Abstract:<p>The immersed boundary method (IBM) with direct forcing is very popular in the simulation of rigid particulate flows. In the IBM, an interaction force is introduced at the interface between fluid and particle in order to approximate the no-slip boundary condition. The interaction force is calculated through dividing the velocity difference (or error) between fluid and particle at the interface by the time step. Here, a dynamic equation for the velocity difference is derived. Analyses on the dynamic equation provide a few new findings: (i) The interaction force is the solution of a least-squares error problem, with the direct implication that the Lagrangian marker distribution has no effect on the large scale flow structure once the distribution of Lagrangian markers become saturated along the interface (i.e., each marker remains properly correlated with all its neighbors); (ii) The Lagrangian volume-weight is a relaxation factor to control how fast the velocity error decays to the ideal value of zero; (iii) The optimal choice of the Lagrangian volume-weight is the largest value permissible by a stability condition. A comprehensive convergence analysis with regard to the spatial and temporal resolution is presented for the velocity error and also for the shear-stress and surface pressure. In three simple canonical problems, it is analytically and numerically shown that the IBM results converge to the theoretical solutions obtained with precise imposition of no-slip and no-penetration boundary conditions. It is observed that it is not necessary to match the Lagrangian marker volume-weight to that of the local Eulerian cell volume and in fact this matching leads to lower than optimal computational efficiency. However, it is found that extremely high Eulerian grid resolution and small time step have to be used to obtain high precision simulation results. Especially, the time step should be inversely proportional to the particle Reynolds number for low Reynolds number flows. For high frequency oscillation problems, the grid size needs to be reduced by a factor of the square root of the frequency, and the time step to be reduced by a factor of the frequency. The theoretical findings here can be used to alleviate the technical difficulties in simulating non-spherical particles by not requiring the Lagrangian marker distribution to match the Eulerian grids and also in the implementation of IBM on non-uniform Eulerian grids. The present work also provides simple practical guidance on the choice of temporal and spatial resolution so as to control the simulation error a priori.</p>

Is Lagrangian weight crucial in the direct forcing immersed boundary method?

An analysis of the spatio-temporal resolution of the immersed boundary method with direct forcing

An Improved Direct-Forcing Immersed Boundary Method with Inward Retraction of Lagrangian Points for Simulation of Particle-Laden Flows

A LB - DF/FD Method for Particle Suspensions

Inflow and Outflow Numerical Simulation Using Least-Square Moving Particle Semi-Implicit Method on GPU

An improved direct-forcing immersed boundary method for simulations of flow and heat transfer in particle-laden flows

Direct Numerical Simulation of Multiple Particles Sedimentation at an Intermediate Reynolds Number

A boundary thickening-based direct forcing immersed boundary method for fully resolved simulation of particle-laden flows

A Fictitious Domain Method for Particulate Flows of Arbitrary Density Ratio

A Parallel Fictitious Domain Method for the Interface-Resolved Simulation of Particle-Laden Flows and Its Application to the Turbulent Channel Flow

A Modified Immersed Boundary Method for Simulations of Fluid–particle Interactions

Full-Scale Solutions to Particle-Laden Flows: Multidirect Forcing and Immersed Boundary Method

Simulation of Flows Around Two Tandem Rotating Cylinders in a Channel Using LB-DF/FD Method

Improvement of Unresolved CFD-DEM by Velocity Field Reconstruction on Unstructured Grids

Simulation of flow around two rotating cylinders in a side-by-side arrangement at long channel using LB-DF/FD method

A Three-Dimensional Fictitious Domain Method for the Simulation of Fluid-Structure Interactions

Interface-resolved Direct Numerical Simulations of the Interactions Between Neutrally Buoyant Spheroidal Particles and Turbulent Channel Flows

Numerical Simulation of Hydrodynamic Focusing of Particles in Straight Channel Flows with the Immersed Boundary-Lattice Boltzmann Method

Interface-resolved Numerical Simulations of Particle-Laden Turbulent Flows in a Vertical Channel Filled with Bingham Fluids

A New Modification of the Immersed‐boundary Method for Simulating Flows with Complex Moving Boundaries

Flow Modulation by Finite-Size Neutrally Buoyant Particles in A Turbulent Channel Flow