A DLM/FD/IB Method for Simulating Compound Cell Interacting with Red Blood Cells in a Microchannel

In this article, a computational model and related methodologies have been tested for simulating the motion of a malaria infected red blood cell (iRBC for short) in Poiseuille flow at low Reynolds numbers. Besides the deformability of the red blood cell membrane, the migration of a neutrally buoyant particle (used to model the malaria parasite inside the membrane) is another factor to determine the iRBC motion. Typically an iRBC oscillates in a Poiseuille flow due to the competition between these two factors. The interaction of an iRBC and several RBCs in a narrow channel shows that, at lower flow speed, the iRBC can be easily pushed toward the wall and stay there to block the channel. But, at higher flow speed, RBCs and iRBC stay in the central region of the channel since their migrations are dominated by the motion of the RBC membrane.     

Transient elastohydrodynamic drag on a particle moving near a deformable wall

We examine the prescribed time-dependent motion of a rigid particle(a sphere or a cylinder) moving in a viscous fluid close toa deformable wall. The fluid motion is described by a nonlinearevolution equation, derived using lubrication theory, whichis solved using numerical and asymptotic methods; a local linearpressure–displacement model describes the wall. When theparticle moves from rest towards the wall, fluid trapping beneaththe particle leads to an overshoot in the normal force on theparticle; a similarity solution is used to describe trappingat early times and a multiregion asymptotic structure describesfluid draining at late times. When the particle is pulled fromrest away from the wall, a peeling process (described by a quasisteadytravelling wave) determines the rate at which fluid can enterthe growing gap between the particle and the wall, leading toa transient adhesive normal force. When a cylinder moves fromrest transversely over the wall, transient peeling motion isagain observed (especially when the wall is initially indented),giving rise to an overshoot in the transverse drag. Simulationsfor a translating sphere show highly nonlinear wall deformationscharacterized by a localized crescent-shaped ridge. Despitegenerating sharp transient deformations, we found no numericalevidence of finite-time choking events.     

An Enriched Biphasic Model for Solute Driven Degradation

Transport of solutes in porous materials plays an important role in many kinds of materials such as biological tissues, porous implants or even soils. In most of the cases the liquid phase in the pores acts as a solvent for one or more solutions. The motion of the solutions is driven by both, the advective and convective transport. The former is related to the fluid phase velocity whereas the letter follows the concentration gradient. The interactions between the solutes and the solid and liquid phase may influence the overall material behavior. Although the solutes often carry electrical charges this paper is focused on neutrally charged solutions. In this contribution the model to describe the solute transport in a fluid saturated porous material is based on the well founded Theory of Porous Media. We will present the basic framework and the governing equations. Finally, we will show a three dimensional numerical example of the solute driven degradation of a skull implant. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Wall effects on the slow steady motion of a particle in a viscous incompressible fluid

In this paper, asymptotic expansions with respect to small Reynolds numbers are proved for the slow steady motion of an arbitrary particle in a viscous, incompressible fluid past a single plane wall. The flow problem is modelled by a certain boundary value problem for the stationary, nonlinear Navier-Stokes equations. The coefficients of these expansions are the solutions of various, linear Stokes problems which can be constructed by single layer potentials and corresponding boundary integral equations on the boundary surface of the particle. Furthermore, some asymptotic estimates at small Reynolds numbers are presented for the slow steady motion of an arbitrary particle in a viscous, incompressible fluid between two parallel, plane walls and in an infinitely long, rectilinear cylinder of arbitrary cross section. In the case of the flow problem with a single plane wall, the paper is based on a thesis, which the author has written under the guidance of Professor Dr. Wolfgang L. Wendland.     

颗粒的有限尺寸对颗粒在均匀各向同性湍流中扩散的影响

本文指出固体颗粒对流体湍流运动的响应有不同的机理,颗粒受大涡的粘性拖动,但受小涡的随机碰撞.基于这种原理,本文计算了有限尺寸的固体颗粒在均匀各向同性湍流中的扩散.结果显示存在着二种相互抵消的效应:颗粒的惯性使颗粒长期扩散系数上升,而颗粒尺寸使颗粒的长期扩散系数下降.     

Brownian dynamics of rigid body by fluctuating hydrodynamic equations

In this work, Brownian dynamics of rigid body in an incompressible fluid with fluctuating hydrodynamic equations is presented. To demonstrate the Brownian motion of rigid body, fluctuating hydrodynamic equations have been coupled with equations of motion of rigid body. Thermal fluctuation is included in the fluid equations via random stress terms unlike the random terms in the conventional Brownian dynamics type approach. Calculation of random stress terms in the fluid is easier in comparison to the random terms in the particle motion. Direct numerical simulation for the Brownian motion of rigid body with a meshfree framework is analysed. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Rotation of an axisymmetric particle in a plane flow

The rotation of an inertialess ellipsoidal particle in a shear flow of a Newtonian fluid has been firstly analyzed by Jeffery [1]. He found that the particle rotates such that the end of its symmetry axis describes a closed periodic orbit. Based on the balance equation of the angular momentum we derived the equation of rotational motion of a cylindrical particle, that is suspended in a plane shear flow field of a viscous fluid, and solved numerically. The rotary inertia is taken into account. The solution is compared with the rotation of a slender particle. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Direct particle motion and interaction modeling method applied to simulate propellant burn

A direct particle motion and particle interaction modeling method was developed to provide an alternative means of capturing the fundamental phenomena occurring during the burning of propellant grains. Individual propellant grains and other moving components are directly incorporated into the computational domain, removing the need for correlations for particle drag and interaction effects. The motion of the individual particles is calculated from the locally acting fluid induced and collision effect forces and moments. Particle/object interactions are handled through a soft particle collision algorithm. Localized mass and energy sources, accompanied by a shrinking particle size, simulate the effects of the combustion process.     

Numerical simulation of particles movement in cellular flows under the influence of magnetic forces

A numerical model of particle motion in fluid flow under the influence of hydrodynamic and magnetic forces is presented. As computational tool, a flow solver based on the Boundary Element Method is used. The Euler-Lagrange formulation of multiphase flow is considered. In the case of a particle with a magnetic moment in a nonuniform external magnetic field, the Kelvin body force acts on a single particle. The derived Lagrangian particle tracking algorithm is used for simulation of dilute suspensions of particles in viscous flows taking into account gravity, buoyancy, drag, pressure gradient, added mass and magnetophoretic force. As a benchmark test case the magnetite particle motion in cellular flow field of water is computed with and without the action of the magnetic force. The effect of the Kelvin force on particle motion and separation from the main flow is studied for a predefined magnetic field and different values of magnetic flux density. (© 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of scour around a submarine pipeline using computational fluid dynamics and discrete element method

Scour under a submarine pipeline can lead to structural failure; hence, a good understanding of the scour mechanism is paramount. Various numerical methods have been proposed to simulate scour, such as potential flow theory and single-phase and two-phase turbulent models. However, these numerical methods have limitations such as their reliance on calibrated empirical parameters and inability to provide detailed information. This paper investigates the use of a coupled computational fluid dynamics-discrete element method (CFD-DEM) model to simulate scour around a pipeline. The novelty of this work is to use CFD-DEM to extract detailed information, leading to new findings that enhance the current understanding of the underlying mechanisms of the scour process. The simulated scour evolution and bed profile are found to be in good agreement with published experimental results. Detailed results include the contours of the fluid velocity and fluid pressure, particle motion and velocity, fluid forces on the particles, and inter-particle forces. The sediment transport rate is calculated using the velocity of each single particle. The quantitative analysis of the bed load layer is also presented. The numerical results reveal three scour stages: onset of scour, tunnel erosion, and lee-wake erosion. Particle velocity and force distributions show that during the tunnel erosion stage, the particle motion and particle–particle interactive forces are particularly intense, suggesting that single-phase models, which are unable to account for inter-particle interactions, may be inadequate. The fluid pressure contours show a distinct pressure gradient. The pressure gradient force is calculated and found to be comparable with the drag force for the onset of scour and the tunnel erosion. However, for the lee-wake erosion, the drag force is shown to be the dominant mechanism for particle movements.     

A Lagrangian stochastic model for nonpassive particle diffusion in turbulent flows

In this paper, we present a Lagrangian stochastic model for heavy particle dispersion in turbulence. The model includes the equation of motion for a heavy particle and a stochastic approach to predicting the velocity of fluid elements along the heavy particle trajectory. The trajectory crossing effect of heavy particles is described by using an Ito type stochastic differential equation combined with a fractional Langevin equation. The comparison of the predicted dispersion of four heavy particles with the observations shows that the model is potentially useful but requires further development.     

Numerical modeling of the movement of a rigid particle in viscous fluid

Modeling the movement of a rigid particle in viscous fluid is a problem physicists and smathematicians have tried to solve since the beginning of this century. A general model for an ellipsoidal particle was first published by Jeffery in the twenties. We exploit the fact that Jeffery was concerned with formulae which can be used to compute numerically the velocity field in the neighborhood of the particle during his derivation of equations of motion of the particle. This is our principal contribution to the subject. After a thorough check of Jeffery's formulae, we coded software for modeling the flow around a rigid particle based on these equations. Examples of its applications are given in conclusion. A practical example is concerned with the simulation of sigmoidal inclusion trails in porphyroblast.     

Coloured noise for dispersion of contaminants in shallow waters

In this article, we explore the application of a set of stochastic differential equations called particle model in simulating the advection and diffusion of pollutants in shallow waters. The Fokker–Planck equation associated with this set of stochastic differential equations is interpreted as an advection–diffusion equation. This enables us to derive an underlying particle model that is exactly consistent with the advection–diffusion equation. Still, neither the advection–diffusion equation nor the related traditional particle model accurately takes into account the short-term spreading behaviour of particles. To improve the behaviour of the model shortly after the deployment of contaminants, a particle model forced by a coloured noise process is developed in this article. The use of coloured noise as a driving force unlike Brownian motion, enables to us to take into account the short-term correlated turbulent fluid flow velocity of the particles. Furthermore, it is shown that for long-term simulations of the dispersion of particles, both the particle due to Brownian motion and the particle model due to coloured noise are consistent with the advection–diffusion equation.     

A numerical model for magnetic chromatography

In this paper, we present details of a mathematical model for magnetic chromatography (MC) systems where strong distorted magnetic fields are used to separate particles from a colloidal mixture. The model simulates the effect of magnetic field gradients on particle motion, and includes calculation of the fluid flow, magnetic field, and particle concentration field. It is based on the finite-volume method (FVM) and uses an expanding-grid technique to handle domains with large aspect ratios. The model has been validated against the results from an analytical model. The numerical model has been used to simulate the performance of a real MC system under various operating conditions.     

Turbulent particle dispersion in an electrostatic precipitator

The behaviour of charged particles in turbulent gas flow in electrostatic precipitators (ESPs) is crucial information to optimise precipitator efficiency. This paper describes a strongly coupled calculation procedure for the rigorous computation of particle dynamics during ESP taking into account the statistical particle size distribution. The turbulent gas flow and the particle motion under electrostatic forces are calculated by using the commercial computational fluid dynamics (CFD) package FLUENT linked to a finite volume solver for the electric field and ion charge. Particle charge is determined from both local electrical conditions and the cell residence time which the particle has experienced through its path. Particle charge density and the particle velocity are averaged in a control volume to use Lagrangian information of the particle motion in calculating the gas and electric fields. The turbulent particulate transport and the effects of particulate space charge on the electrical current flow are investigated. The calculated results for poly-dispersed particles are compared with those for mono-dispersed particles, and significant differences are demonstrated.     

The Maxey–Riley equation: Existence,uniqueness and regularity of solutions

The Maxey–Riley equation describes the motion of an inertial (i.e., finite-size) spherical particle in an ambient fluid flow. The equation is a second-order, implicit integro-differential equation with a singular kernel, and with a forcing term that blows up at the initial time. Despite the widespread use of the equation in applications, the basic properties of its solutions have remained unexplored. Here we fill this gap by proving local existence and uniqueness of mild solutions. For certain initial velocities between the particle and the fluid, the results extend to strong solutions. We also prove continuous differentiability of the mild and strong solutions with respect to their initial conditions. This justifies the search for coherent structures in inertial flows using the Cauchy–Green strain tensor.     

Thermocapillary motion of deformable fluid particles in tubes

The boundary integral technique is used to study the effect of deformation on the steady, creeping, thermocapillary migration of a fluid particle under conditions of axisymmetry, negligible thermal convection and an insulated tube wall. The spherical radius of the fluid particle (i.e. the radius as if the particle were a sphere, a ′= (3V _p ^′ /4π)^1/3, V _p ^′ is the particle volume) and that of the tube are denoted, respectively, by a′and b′. For small capillary numberCa = 0.05, only for a large fluid particle (a′/b′ = 0.8) is deformation significant. Fora′/b′= 0.8, hydrodynamic stresses squeeze the particle, reduce the interaction of the particle with the wall and thereby increase the terminal velocity. For small particles a′/b′< 0.8 and Ca = 0.05 the fluid particles translate as spheres, due to the fact that the fluid particle is too far away from the wall to be subject to distending hydrodynamic stresses. The deformable particle moves faster than a spherical one in the thermocapillary migration. The increase in velocity with capillary number is larger for thermocapillary motion than for buoyancy.     

湍流边界层中固体小颗粒湍流运动的Lagrangian模型

给出了固体小颗粒在边界层中的Lagrangian运动方程,方程中包括受壁面影响的粘性阻力,Saffman升力及Magus升力等.使用频谱法,得到了颗粒响应流体的Lagrangian能谱的表达式,使用这些结果研究了各种响应特性.本文的结果清楚地表明了固体个颗粒在湍流扩散过程中,其湍流扩散是可能大于流体的.     

Effects of Poiseuille flow on peristaltic transport of a particulate suspension

The problem of peristaltic transport induced by sinusoidal waves of a particle-fluid mixture in the presence of a Poiseuille flow, is analysed. The governing equations of motion resulting from the Navier-Stokes equations for both the fluid and particle phases are solved and closed form solutions are obtained for limiting values of Reynolds number, wave number and the Poiseuille flow parameter while the method of Frobenius series solution is used for the general case. It is found that the mean flow is strongly dependent on the Poiseuille flow parameter. The effects of particle concentration in the fluid is well discharged throughout the analysis and the results are compared with the other studies in the literature.     

Settling and asymptotic motion of aerosol particles in a cellular flow field

Summary This paper presents a proof that given a dilute concentration of aerosol particles in an infinite, periodic, cellular flow field, arbitrarily small inertial effects are sufficient to induce almost all particles to settle. It is shown that when inertia is taken as a small parameter, the equations of particle motion admit a slow manifold that is globally attracting. The proof proceeds by analyzing the motion on this slow manifold, wherein the flow is a small perturbation of the equation governing the motion of fluid particles. The perturbation is supplied by the inertia, which here occurs as a regular parameter. Further, it is shown that settling particles approach a finite number of attracting periodic paths. The structure of the set of attracting paths, including the nature of possible bifurcations of these paths and the resulting stability changes, is examined via a symmetric one-dimensional map derived from the flow.