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.

Simulation of sedimentation of two spheres with different densities in a square tube

Flow Patterns in the Sedimentation of an Elliptical Particle

Freely Rising or Falling of a Sphere in a Square Tube at Intermediate Reynolds Numbers

Dynamic Simulation of Sphere Motion in a Vertical Tube

Discontinuity in the Sedimentation System with Two Particles Having Different Densities in a Vertical Channel.

Settling Behavior of Two Particles with Different Densities in a Vertical Channel

Numerical Simulation Of Two Elliptical Particles' Sedimentation With Inertial Using Lattice Boltzmann Method

Study on the Interaction of Sedimenting Cylindrical Particles in Still Fluid

Numerical simulation of the sedimentation of cylindrical pollutant particles in fluid.

Sedimentation of a Single Particle Between Two Parallel Walls

Sedimentation of an Ellipsoidal Particle in Narrow Tubes

Settling of two rigidly connected spheres

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

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

Sedimentation of rectangular particle in Oldroyd-B fluid

SEDIMENTATION OF CIRCULAR PARTICLES IN OLDROYDB FLUID

Interaction Between Two Unequal Particles at Intermediate Reynolds Numbers: A Pattern of Horizontal Oscillatory Motion.

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

Sedimentation of a spherical squirmer in a square tube under gravity

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

Numerical Simulation for Flow Over two Circular Cylinders in Side-by-side Arrangement with Lattice Boltzmann Method