Direct methods fixed step size

In the IBM, the presence of the solid boundary (fixed or moving) in the fluid can be represented by a virtual body force field -rp( ) applied on the computational grid at the vicinity of solid-flow interface. Considering the stability and efficiency in a 3-D simulation, the direct forcing scheme is adopted in this model. Details of this scheme are introduced in Section II.B. In this study, a new velocity interpolation method is developed based on the particle level-set function (p), which is shown in Fig. 20. At each time step of the simulation, the fluid-particle boundary condition (no-slip or free-slip) is imposed on the computational cells located in a small band across the particle surface. The thickness of this band can be chosen to be equal to 3A, where A is the mesh size (assuming a uniform mesh is used). If a grid point (like p and q in Fig. 20), where the velocity components of the control volume are defined, falls into this band, that is... 

---