Search

Abstract: A numerical scheme for fully resolved simulation of particle–fluid systems with freely moving rigid particles is developed. The approach is based on a fictitious domain method wherein the entire particle–fluid domain is assumed to be an incompressible fluid but with variable density. The flow inside the particle domain is constrained to be a rigid body motion using an additional rigidity constraint in a fractional step scheme. The rigidity constraint force is obtained based on the fast computation technique proposed by Sharma and Patankar (2005) [1]. The particle is assumed to be made up of material points moving on a fixed background mesh where the fluid flow equations are solved. The basic finite-volume solver is based on a co-located grid incompressible but variable density flow. The incompressibility constraint is imposed by solving a variable-coefficient pressure equation. Use of density-weighted reconstruction of the pressure gradients was found to give a stable scheme for high density ratio particle–fluid systems. Various verification and validation test cases on fixed and freely moving particles are performed to show that the numerical approach is accurate and stable for a wide range ( – ) of particle–fluid density ratios. [Copyright &y& Elsevier]

Abstract: A fictitious-domain based formulation for fully resolved simulations of arbitrary shaped, freely moving rigid particles in unsteady flows is presented. The entire fluid–particle domain is assumed to be an incompressible, but variable density, fluid. The numerical method is based on a finite-volume approach on a co-located, Cartesian grid together with a fractional step method for variable density, low-Mach number flows. The flow inside the fluid region is constrained to be divergence-free for an incompressible fluid, whereas the flow inside the particle domain is constrained to undergo rigid body motion. In this approach, the rigid body motion constraint is imposed by avoiding the explicit calculation of distributed Lagrange multipliers and is based upon the formulation developed by Patankar [N. Patankar, A formulation for fast computations of rigid particulate flows, Center for Turbulence Research Annual Research Briefs 2001 (2001) 185–196]. The rigidity constraint is imposed and the rigid body motion (translation and rotational velocity fields) is obtained directly in the context of a two-stage fractional step scheme. The numerical approach is applied to both imposed particle motion and fluid–particle interaction problems involving freely moving particles. Grid and time-step convergence studies are performed to evaluate the accuracy of the approach. Finally, simulation of rigid particles in a decaying isotropic turbulent flow is performed to study the feasibility of simulations of particle-laden turbulent flows. [Copyright &y& Elsevier]

• A general framework for the undisturbed flow in point-particle approach. • Two approaches to estimate the particle disturbance fields. • Applicable to arbitrary grid shapes, unbounded and wall-bounded flows. • Any particle size-to-grid ratio, particle density, and Reynolds numbers at low volume loadings. A general, two-way coupled, point-particle formulation that accounts for the disturbance created by the dispersed particles in obtaining the undisturbed fluid flow field needed for accurate computation of the force closure models is presented. Specifically, equations for the disturbance field created by the presence of particles are first derived based on the inter-phase momentum coupling force in a finite-volume formulation. Solution to the disturbance field is obtained using two approaches: (i) direct computation of the disturbance velocity and pressure using the reaction force due to particles at computational control volumes, and (ii) a linearized, approximate computation of the disturbance velocity field, specifically applicable for low Reynolds number flows. In both approaches, the computed disturbance field is used to obtain the undisturbed fluid velocity necessary to model the aerodynamic forces on the particle. The two approaches are thoroughly evaluated for a single particle in an unbounded and wall-bounded flow on uniform, anisotropic, as well as unstructured grids to show accurate computation of the particle motion and inter-phase coupling. The approach is straightforward and can be applied to any Euler-Lagrange formulations of particle-laden flows. [ABSTRACT FROM AUTHOR]

The accuracy of Euler-Lagrange point-particle models employed in particle-laden fluid flow simulations depends on accurate estimation of the particle force through closure models. Typical force closure models require computation of the slip velocity at the particle location, which in turn requires accurate estimation of the undisturbed fluid velocity. Such an undisturbed velocity is not readily available when the fluid and particle phases are two-way coupled, due to the disturbance created by the particle's force in the nearby fluid velocity field. A common practice is to use the disturbed velocity to compute the particle force which can result in errors as much as 100% in predicting the particle dynamics. In this work, a correction scheme is developed that facilitates accurate estimation of the undisturbed fluid velocity in particle-laden fluid flows with and without no-slip walls. The model is generic and can handle particles of different size and density, arbitrary interpolation and projection functions, anisotropic grids with large aspect ratios, and wall-bounded flows. The present correction scheme is motivated by the recent work of Esmaily & Horwitz (JCP, 2018) on unbounded particle-laden flows. Modifications necessary for wall-bounded flows are developed such that the undisturbed fluid velocity at any wall distance is accurately recovered, asymptotically approaching the result of unbounded schemes for particles far away from walls. A detailed series of verification tests was conducted on settling velocity of a particle in parallel and perpendicular motions to a no-slip wall. A range of flow parameters and grid configurations; involving anisotropic rectilinear grids with aspect ratios typically encountered in particle-laden turbulent channel flows was considered in detail. When the wall effects are accounted for, the present correction scheme reduces the errors in predicting the near-wall particle motion by one order of magnitude smaller values compared to the unbounded correction schemes. [ABSTRACT FROM AUTHOR]