Abstract:In this study, we present a fictitious domain method for the dynamic simulation of particle motion in a Bingham viscoplastic fluid at moderate Bingham numbers. Our method is built on the framework established by Glowinski and his coworkers, in the sense that we use their formulation and their operator-splitting idea to simplify the computation, but differs from their method in both spatial and temporal discretizations of the governing equations. Concerning the space scheme, we use a finite-difference method to discretize equations and a collocation-element method to enforce the rigid-body motion constraint inside the particle boundaries. Concerning the time scheme, the combined system is decoupled into three sub-systems: a Navier–Stokes problem, a plasticity problem and a rigid-body-motion problem, and we solve the Navier–Stokes problem with the classic projection method. The plasticity and rigid-body-motion multipliers at the previous time-level are retained in the Navier–Stokes problem to reduce the splitting error. In addition, the present study shows that retaining a viscous diffusion term in the plasticity problem is not necessary for the convergence of the iteration. The method is verified by comparing our results on the lid-driven cavity flow, the drag coefficient for a sphere settling in a tube and the hydrodynamic interactions between two spheres translating along the tube axis to the data available in the literature. Our results confirm that the drag coefficient for a sphere settling in a Bingham fluid at non-zero Reynolds numbers can be well correlated with an effective Reynolds number. For two approaching spheres, there exists a critical separation distance above which a drag-reduction is observed and below which a drag-enhancement takes place, compared to the drag for a single sphere. The drag-reduction however does not happen to a sphere falling towards a solid wall at the moderate Bingham numbers we studied. These observations are explained by a consideration of the competition between a shear-thinning plastic force and a repulsive lubrication force on the sphere occurring in the squeezing flow.

Research on Aerodynamic Characteristics of the Sphere in Transitional Region by Bridging-Function and DSMC Method

Simulation of Transitional Flow in a 2D Groove by DSMC Method

The Overall Numerical Simulation of Both High-Altitude Plume Flowfield and Inflow Region with DSMC Method

Engineering Calculation of Heat Flux at the Stagnation Point of Cylinder in Transitional Region

Modeling and numerical simulation of gas-solid jet in the transitional-rarefied regime

Numerical study of uniform flow over sphere at intermediate Reynolds numbers

Study on Flow Past Two Spheres in Tandem Arrangement Using a Local Mesh Refinement Virtual Boundary Method

The study of hybrid DSMC/EPSM method of supersonic flat plate

Numerical investigations of wake and force for flow past a sphere

Study on transitions for flow past two spheres in tandem arrangement

Scale-adaptive simulation of the separated flow past a 90°-inclined prolate spheroid

Interaction between a stationary sphere and turbulent flow in a boundary layer

Dynamic Simulation of Sphere Motion in a Vertical Tube

Fluidization of the Spherocylindrical Particles: Comparison of Multi-Sphere and Bond-Sphere Models

Comparison of Transitional Turbulence Models for Numerical Simulation of Flow Past a Sphere in a Pipe

Investigation of Supersonic Turbulent Flows over a Sphere by Fully Resolved Direct Numerical Simulation

Wake Structures of Two Spheres in Tandem Arrangement at Various Gaps for Re = 300

New Implement of Pressure Boundary Conditions for DSMC

The Influence of Advection Schemes and Turbulence Closure Models on Drag Coefficient Calculation Around a Circular Cylinder at High Reynolds Number

Computational Analysis of Sphere Wakes in a Linearly Stratified Fluid

A Fictitious Domain Method for Dynamic Simulation of Particle Sedimentation in Bingham Fluids