Mathematical modelling of two‐phase non‐Newtonian flow in a helical pipe

Governing equations for a two‐phase 3D helical pipe flow of a non‐Newtonian fluid with large particles are derived in an orthogonal helical coordinate system. The Lagrangian approach is utilized to model solid particle trajectories. The interaction between solid particles and the fluid that carries them is accounted for by a source term in the momentum equation for the fluid. The force‐coupling method (FCM), developed by M.R. Maxey and his group, is adopted; in this method the momentum source term is no longer a Dirac delta function but is spread on a numerical mesh by using a finite‐sized envelop with a spherical Gaussian distribution. The influence of inter‐particle and particle–wall collisions is also taken into account. Copyright © 2005 John Wiley & Sons, Ltd.     

Direct Numerical Simulation of particulate flow with heat transfer

Numerical simulations of heat transfer in non-isothermal particulate flows are important to better understand the flow pattern. The complexity of numerical algorithms coupling the heat and mass transfer and the considerable computational resources required limit the number of such direct simulations that can be reasonably performed. We suggest a Distributed Lagrange Multiplier/Fictitious Domain (DLM/FD) method to compute the temperature distribution and the heat exchange between the fluid and solid phases. The Boussinesq approximation is considered for the flow/temperature fields coupling. We employ a Finite Element Method (FEM) to solve the fluid flow conservation equations for mass, momentum and energy. The motion of particles is computed by a Discrete Element Method (DEM). On each particle, heat transfer is solved using a FEM. For each class of particles, we generate a single FEM grid and translate/rotate it at each time step to match the physical configuration of each particle. Distributed Lagrange multipliers for both the velocity and temperature fields are introduced to treat the fluid/solid interaction. This work is an extension of the method we proposed in Yu et al. (2006). Two two-dimensional (2D) test cases are proposed to validate the implementation by comparing our computational results with those reported in the literature. Finally, the sedimentation of a single sphere in a semi-infinite channel is presented and the results are discussed.     

No‐slip consistent immersed boundary particle tracking to simulate impaction filtration in porous media

In this paper, we present a new method for simulating the motion of a disperse particle phase in a carrier gas through porous media. We assume a sufficiently dilute particle‐laden flow and compute, independently of the disperse phase, the steady laminar fluid velocity using the immersed boundary method. Given the velocity of the carrier gas, the equations of motion for the particles experiencing the Stokes drag force are solved to determine their trajectories. The ‘no‐slip consistent’ particle tracking algorithm avoids possible numerical filtration of very small particles due to the nonzero velocity field at the solid–fluid interface introduced by the immersed boundary method. This physically consistent tracking allows a reliable estimation of the filtration efficiency of porous filters due to inertial impaction. We illustrate and test our new approach for model porous media consisting of a structured array of aligned rectangular fibers, arranged in line and staggered. In the staggered geometry, the effect of the residual velocity at the solid–fluid interface is significant for particles with low inertia. Without adopting the developed no‐slip consistent numerical method, an artificial numerical filtration is observed, which becomes dominant for small enough particles. For both the in line and the staggered geometries, the filtration rate depends quite strongly and non monotonically on the particle inertia. This is expressed most clearly in the staggered arrangement in which a very strong increase in the filtration efficiency is observed at a well‐defined critical droplet size, corresponding to a qualitative change in the dominant particle paths in the porous medium. Copyright © 2013 John Wiley & Sons, Ltd.     

THE DERIVATION OF REYNOLDS TRANSPORT EQUATION FOR ARBITRARY MOTION CONTROL VOLUME BASED ON THE DYNAMIC BOUNDARY CALCULUS

严格而言,流体力学中所有守恒定律均是针对物质体系的(或称流体系统),如质量、动量、动量矩和能量等守恒定律。如果跟随物质体系描述和表征流体质点系的运动行为,即为Lagrange描述方法;如果把物质体系的运动和守恒定律转换到空间坐标系中,即为人们常说的Euler描述方法。因此,对于具体考察(跟随的)的流体物质系统而言,各守恒定律存在由物质体系表征到空间体系表征的转换,这个转换关系就是著名的Reynolds输运方程。本文从动边界微积分关系式出发,系统推导了在不同运动速度控制体上的雷诺输运方程,并通过讨论进一步阐明各种不同形式输运方程的物理意义。     

Numerical study of concentrated fluid–particle suspension flow in a wavy channel

The particle migration effects and fluid–particle interactions occurring in the flow of highly concentrated fluid–particle suspension in a spatially modulated channel have been investigated numerically using a finite volume method. The mathematical model is based on the momentum and continuity equations for the suspension flow and a constitutive equation accounting for the effects of shear‐induced particle migration in concentrated suspensions. The model couples a Newtonian stress/shear rate relationship with a shear‐induced migration model of the suspended particles in which the local effective viscosity is dependent on the local volume fraction of solids. The numerical procedure employs finite volume method and the formulation is based on diffuse‐flux model. Semi‐implicit method for pressure linked equations has been used to solve the resulting governing equations along with appropriate boundary conditions. The numerical results are validated with the analytical expressions for concentrated suspension flow in a plane channel. The results demonstrate strong particle migration towards the centre of the channel and an increasing blunting of velocity profiles with increase in initial particle concentration. In the case of a stenosed channel, the particle concentration is lowest at the site of maximum constriction, whereas a strong accumulation of particles is observed in the recirculation zone downstream of the stenosis. The numerical procedure applied to investigate the effects of concentrated suspension flow in a wavy passage shows that the solid particles migrate from regions of high shear rate to low shear rate with low velocities and this phenomenon is strongly influenced by Reynolds numbers and initial particle concentration. Copyright © 2008 John Wiley & Sons, Ltd.     

A parallel computational method for simulating two‐phase gel dynamics on a staggered grid

We develop a parallel computational algorithm for simulating models of gel dynamics where the gel is described by two phases, a networked polymer and a fluid solvent. The models consist of transport equations for the two phases, two coupled momentum equations, and a volume‐averaged incompressibility constraint. Multigrid with Vanka‐type box‐relaxation scheme is used as preconditioner for the Krylov subspace solver (GMRES) to solve the momentum and incompressibility equations. Through numerical experiments of a model problem, the efficiency, robustness and scalability of the algorithm are illustrated. Copyright © 2008 John Wiley & Sons, Ltd.     

An SPH model for free surface flows with moving rigid objects

This paper presents a computational model for free surface flows interacting with moving rigid bodies. The model is based on the SPH method, which is a popular meshfree, Lagrangian particle method and can naturally treat large flow deformation and moving features without any interface/surface capture or tracking algorithm. Fluid particles are used to model the free surface flows which are governed by Navier–Stokes equations, and solid particles are used to model the dynamic movement (translation and rotation) of moving rigid objects. The interaction of the neighboring fluid and solid particles renders the fluid–solid interaction and the non‐slip solid boundary conditions. The SPH method is improved with corrections on the SPH kernel and kernel gradients, enhancement of solid boundary condition, and implementation of Reynolds‐averaged Navier–Stokes turbulence model. Three numerical examples including the water exit of a cylinder, the sinking of a submerged cylinder and the complicated motion of an elliptical cylinder near free surface are provided. The obtained numerical results show good agreement with results from other sources and clearly demonstrate the effectiveness of the presented meshfree particle model in modeling free surface flows with moving objects. Copyright © 2013 John Wiley & Sons, Ltd.     

An immersed‐boundary finite‐volume method for direct simulation of flows with suspended paramagnetic particles

In the present study, we have proposed an immersed‐boundary finite‐volume method for the direct numerical simulation of flows with inertialess paramagnetic particles suspended in a nonmagnetic fluid under an external magnetic field without the need for any model such as the dipole–dipole interaction. In the proposed method, the magnetic field (or force) is described by the numerical solution of the Maxwell equation without current, where the smoothed representation technique is employed to tackle the discontinuity of magnetic permeability across the particle–fluid interface. The flow field, on the other hand, is described by the solution of the continuity and momentum equations, where the discrete‐forcing‐based immersed‐boundary method is employed to satisfy the no‐slip condition at the interface. To validate the method, we performed numerical simulations on the two‐dimensional motion of two and three paramagnetic particles in a nonmagnetic fluid subjected to an external uniform magnetic field and then compared the results with the existing finite‐element and semi‐analytical solutions. Comparison shows that the proposed method is robust in the direct simulation of such magnetic particulate flows. This method can be extended to more general flows without difficulty: three‐dimensional particulate flows, flows with a great number of particles, or flows under an arbitrary external magnetic field. Copyright © 2010 John Wiley & Sons, Ltd.     

Prediction of turbulent gas–solids flow in curved ducts using the Eulerian–Lagrangian method

The flow of particulate two‐phase flow mixtures occur in several components of solid fuel combustion systems, such as the pressurised fluidised bed combustors (PFBC) and suspension‐fired coal boilers. A detailed understanding of the mixture characteristics in the conveying component can aid in refining and optimising its design. In this study, the flow of an isothermal, dilute two‐phase particulate mixture has been examined in a high curvature duct, which can be representative of that transporting the gas–solid mixture from the hot clean‐up section to the gas turbine combustor in a PFBC plant. The numerical study has been approached by utilising the Eulerian–Lagrangian methodology for describing the characteristics of the fluid and particulate phases. By assuming that the mixture is dilute and the particles are spherical, the governing particle momentum equations have been solved with appropriately prescribed boundary conditions. Turbulence effects on the particle dispersion were represented by a statistical model that accounts for both the turbulent eddy lifetime and the particle transit time scales. For the turbulent flow condition examined it was observed that mixtures with small particle diameters had low interphase slip velocities and low impaction probability with the pipe walls. Increasing the particle diameters (>50 μm) resulted in higher interphase slip velocities and, as expected, their impaction probability with the pipe walls was significantly increased. The particle dispersion is significant for the smaller sizes, whereas the larger particles are relatively insensitive to the gas turbulence. The main particle impaction region, and locations most prone to erosion damage, is estimated to be within an outer duct length of two to six times the duct diameter, when the duct radius of curvature to the duct diameter ratio is equal to unity. Copyright © 1999 John Wiley & Sons, Ltd.     

A high-resolution finite-difference method for simulating two-fluid,viscoelastic gel dynamics

An important class of gels are those composed of a polymer network and fluid solvent. The mechanical and rheological properties of these two-fluid gels can change dramatically in response to temperature, stress, and chemical stimulus. Because of their adaptivity, these gels are important in many biological systems, e.g. gels make up the cytoplasm of cells and the mucus in the respiratory and digestive systems, and they are involved in the formation of blood clots. In this study we consider a mathematical model for gels that treats the network phase as a viscoelastic fluid with spatially and temporally varying material parameters and treats the solvent phase as a viscous Newtonian fluid. The dynamics are governed by a coupled system of time-dependent partial differential equations which consist of transport equations for the two phases, constitutive equations for the viscoelastic stresses, two coupled momentum equations for the velocity fields of the two fluids, and a volume-averaged incompressibility constraint. We present a numerical method based on a staggered grid, second order finite-difference discretization of the momentum equations and a high-resolution unsplit Godunov method for the transport equations. The momentum and incompressibility equations are solved in a coupled manner with the Generalized Minimum Residual (GMRES) method using a multigrid preconditioner based on box-relaxation. We present results on the accuracy and robustness of the method together with an illustration of the interesting behavior of this gel model for the four-roll mill problem.     

Study of a new macro-particle model for the low-temperature pyrolysis of dried wood chips

 This paper presents a transient one-dimensional mathematical model which simulates the pyrolysis of a single dried wood particle. The porous wood particle is considered as a two-phase system: the solid phase consisting of wood and char and the gas phase consisting of volatiles and tar. Conservation equations for mass, momentum and energy are developed for each phase. Chemical processes are described through an existing one-stage three-reactions scheme, leading separately to char, tar and volatiles evolution dynamics. The variation of transport and physical properties with temperature and with composition is presented by algebraic equations. The model presented in this paper is more advanced than the current models found in literature, since it contains physical effects not included in past models, such as cross diffusion, differing solid and gas phase temperatures and a transient gas phase momentum equation incorporating the wall friction experienced by a fluid flowing through a porous medium. Furthermore, an adequate reference system for enthalpy, based on temperature dependent reaction heats, is used. The mathematical equations with initial and boundary conditions are solved numerically by means of a commercial CFD code (PHOENICS). The validity of the pyrolysis kinetics scheme is examined through comparison with experimental data. Furthermore, the macro-particle model is validated by (1) examining the limitations and importance of the newly-modelled effects (different solid phase and gas phase temperature, cross diffusion and wall friction) and (2) making a comparison between predicted and experimental results for large particles. Received on 18 December 2000     

Simulating 3D deformable particle suspensions using lattice Boltzmann method with discrete external boundary force

A method for direct numerical analysis of three‐dimensional deformable particles suspended in fluid is presented. The flow is computed on a fixed regular ‘lattice’ using the lattice Boltzmann method (LBM), where each solid particle is mapped onto a Lagrangian frame moving continuously through the domain. Instead of the bounce‐back method, an external boundary force (EBF) is used to impose the no‐slip boundary condition at the fluid–solid interface for stationary or moving boundaries. The EBF is added directly to the lattice Boltzmann equation. The motion and orientation of the particles are obtained from Newtonian dynamics equations. The advantage of this approach is outlined in comparison with the standard and higher‐order interpolated bounce‐back methods as well as the LBM immersed‐boundary and the volume‐of‐fluid methods. Although the EBF method is general, in this application, it is used in conjunction with the lattice–spring model for deformable particles. The methodology is validated by comparing with experimental and theoretical results. Copyright © 2009 John Wiley & Sons, Ltd.     

A multiphase DNS approach for handling solid particles motion with heat transfer

In the current work we propose a multiphase DNS method capable of resolving the motion of solid particles coupled with heat transfer effects. The method is based on solving a shared set of momentum and energy balance equations for the carrier phase and the particulate phase. Individual particles are tracked using a number of volume fraction advection equations. The proposed method is in very good agreement with the available data in the literature for the following cases: isothermal particle motion (in the presence of walls and other particles), natural convection around a stationary particle and solid particles motion accompanied with heat transfer effects. In addition, we show that the method is inherently capable of handling deformable particles (i.e. droplets and bubbles) co-existing with solid particles. The method is thus well suited to deal with challenging multiphase systems, such as diesel spray combustion with soot formation, spray drying with particle nucleation, and biological treatment of waste water.     

A discrete forcing method dedicated to moving bodies in two‐phase flow

The numerical simulation of interaction between structures and two‐phase flows is a major concern for many industrial applications. To address this challenge, the motion of structures has to be tracked accurately. In this work, a discrete forcing method based on a porous medium approach is proposed to follow a nondeformable rigid body with an imposed velocity by using a finite‐volume Navier‐Stokes solver code dedicated to multiphase flows and based on a two‐fluid approach. To deal with the action reaction principle at the solid wall interfaces in a conservative way, a porosity is introduced allowing to locate the solid and insuring no diffusion of the fluid‐structure interface. The volumetric fraction equilibrium is adapted to this novelty. Mass and momentum balance equations are formulated on a fixed Cartesian grid. Interface tracking is addressed in detail going from the definition of the porosity to the changes in the discretization of the momentum balance equation. This so‐called time‐ and space‐dependent porosity method is then validated by using analytical and elementary test cases.     

Numerical simulation of 3D free surface flows,with multiple incompressible immiscible phases. Applications to impulse waves

A numerical method for the solution to the density‐dependent incompressible Navier–Stokes equations modeling the flow of N immiscible incompressible liquid phases with a free surface is proposed. It allows to model the flow of an arbitrary number of liquid phases together with an additional vacuum phase separated with a free surface. It is based on a volume‐of‐fluid approach involving N indicator functions (one per phase, identified by its density) that guarantees mass conservation within each phase. An additional indicator function for the whole liquid domain allows to treat boundary conditions at the interface between the liquid domain and a vacuum. The system of partial differential equations is solved by implicit operator splitting at each time step: first, transport equations are solved by a forward characteristics method on a fine Cartesian grid to predict the new location of each liquid phase; second, a generalized Stokes problem with a density‐dependent viscosity is solved with a FEM on a coarser mesh of the liquid domain. A novel algorithm ensuring the maximum principle and limiting the numerical diffusion for the transport of the N phases is validated on benchmark flows. Then, we focus on a novel application and compare the numerical and physical simulations of impulse waves, that is, waves generated at the free surface of a water basin initially at rest after the impact of a denser phase. A particularly useful application in hydraulic engineering is to predict the effects of a landslide‐generated impulse wave in a reservoir. Copyright © 2014 John Wiley & Sons, Ltd.     

Dynamic simulation of sedimentation of solid particles in an Oldroyd-B fluid

In this paper we present a two-dimensional numerical study of the viscoelastic effects on the sedimentation of particles in the presence of solid walls or another particle. The Navier-Stokes equations coupled with an Oldroyd-B model are solved using a finite-element method with the EVSS formalism, and the particles are moved according to their equations of motion. In a vertical channel filled with a viscoelastic fluid, a particle settling very close to one side wall experiences a repulsion from the wall; a particle farther away from the wall is attracted toward it. Thus a settling particle will approach an eccentric equilibrium position, which depends on the Reynolds and Deborah numbers. Two particles settling one on top of the other attract and form a doublet if their initial separation is not too large. Two particles settling side by side approach each other and the doublet also rotates till the line of centers is aligned with the direction of sedimentation. The particle-particle interactions are in qualitative agreement with experimental observations, while the wall repulsion has not been documented in experiments. The driving force for lateral migrations is shown to correlate with the pressure distribution on the particle's surface. As a rule, viscoelasticity affects the motion of particles by modifying the pressure distribution on their surface. The direct contribution of viscoelastic normal stresses to the force and torque is not important.     

Study of water waves with submerged obstacles using a vortex method with Helmholtz decomposition

The present study develops a 2‐D numerical scheme that combines the vortex method and the boundary integral method by a Helmholtz decomposition to investigate the interaction of water waves with submerged obstacles. Viscous effects and generation of vorticity on the free surface are neglected. The second kind of Fredholm integral equations that govern the strengths of vortex sheets along boundaries are solved iteratively. Vorticity is convected and diffused in the fluid via a Lagrangian vortex (blob) method with varying cores, using the particle strength exchange method for diffusion, with particle redistribution. A grid‐convergence study of the numerical method is reported. The inviscid part of the method and the simulation of the free‐surface motion are tested using two calculations: solitary wave propagation in a uniform channel and a moving line vortex in the fluid. Finally, the full model is verified by simulating periodic waves travelling over a submerged rectangular obstacle using nonuniform vortex blobs with a mapping of the redistribution lattice. Overall, the numerical model predicts the vortices' evolution and the free‐surface motion reasonably well. Copyright © 2008 John Wiley & Sons, Ltd.     

Constrained moving least‐squares immersed boundary method for fluid‐structure interaction analysis

A numerical method is presented for the analysis of interactions of inviscid and compressible flows with arbitrarily shaped stationary or moving rigid solids. The fluid equations are solved on a fixed rectangular Cartesian grid by using a higher‐order finite difference method based on the fifth‐order WENO scheme. A constrained moving least‐squares sharp interface method is proposed to enforce the Neumann‐type boundary conditions on the fluid‐solid interface by using a penalty term, while the Dirichlet boundary conditions are directly enforced. The solution of the fluid flow and the solid motion equations is advanced in time by staggerly using, respectively, the third‐order Runge‐Kutta and the implicit Newmark integration schemes. The stability and the robustness of the proposed method have been demonstrated by analyzing 5 challenging problems. For these problems, the numerical results have been found to agree well with their analytical and numerical solutions available in the literature. Effects of the support domain size and values assigned to the penalty parameter on the stability and the accuracy of the present method are also discussed.