Lattice Boltzmann simulations of viscoplastic fluid flows through complex flow channels

We present the results of lattice Boltzmann (LB) simulations for the planar-flow of viscoplastic fluids through complex flow channels. In this study, the Bingham and Casson model fluids are covered as viscoplastic fluid. The Papanastasiou (modified Bingham) model and the modified Casson model are employed in our LB simulations. The Bingham number is an essential physical parameter when considering viscoplastic fluid flows and the modified Bingham number is proposed for modified viscoplastic models. When the value of the modified Bingham number agrees with that of the “normal” Bingham number, viscoplastic fluid flows formulated by modified viscoplastic models strictly reproduce the flow behavior of the ideal viscoplastic fluids. LB simulations are extensively performed for viscoplastic fluid flows through complex flow channels with rectangular and circular obstacles. It is shown that the LB method (LBM) allows us to successfully compute the flow behavior of viscoplastic fluids in various complicated-flow channels with rectangular and circular obstacles. For even low Re and high Bn numbers corresponding to plastic-property dominant condition, it is clearly manifested that the viscosity for both the viscoplastic fluids is largely decreased around solid obstacles. Also, it is shown that the viscosity profile is quite different between both the viscoplastic fluids due to the inherent nature of the models. The viscosity of the Bingham fluid sharply drops down close to the plastic viscosity, whereas the viscosity of the Casson fluid does not rapidly fall. From this study, it is demonstrated that the LBM can be also an effective methodology for computing viscoplastic fluid flows through complex channels including circular obstacles.     

Gas-kinetic numerical method for solving mesoscopic velocity distribution function equation

A gas-kinetic numerical method for directly solving the mesoscopic velocity distribution function equation is presented and applied to the study of three-dimensional complex flows and micro-channel flows covering various flow regimes. The unified velocity distribution function equation describing gas transport phenomena from rarefied transition to continuum flow regimes can be presented on the basis of the kinetic Boltzmann–Shakhov model equation. The gas-kinetic finite-difference schemes for the velocity distribution function are constructed by developing a discrete velocity ordinate method of gas kinetic theory and an unsteady time-splitting technique from computational fluid dynamics. Gas-kinetic boundary conditions and numerical modeling can be established by directly manipulating on the mesoscopic velocity distribution function. A new Gauss-type discrete velocity numerical integration method can be developed and adopted to attack complex flows with different Mach numbers. HPF parallel strategy suitable for the gas-kinetic numerical method is investigated and adopted to solve three-dimensional complex problems. High Mach number flows around three-dimensional bodies are computed preliminarily with massive scale parallel. It is noteworthy and of practical importance that the HPF parallel algorithm for solving three-dimensional complex problems can be effectively developed to cover various flow regimes. On the other hand, the gas-kinetic numerical method is extended and used to study micro-channel gas flows including the classical Couette flow, the Poiseuille- channel flow and pressure-driven gas flows in two-dimensional short micro-channels. The numerical experience shows that the gas-kinetic algorithm may be a powerful tool in the numerical simulation of micro-scale gas flows occuring in the Micro-Electro-Mechanical System (MEMS). The project supported by the National Natural Science Foundation of China (90205009 and 10321002), and the National Parallel Computing Center in Beijing. The English text was polished by Yunming Chen.     

Local vorticity computation approach in double distribution functions based lattice Boltzmann methods for flow and scalar transport

Computation of vorticity, or the skew-symmetric velocity gradient tensor, in conjunction with the strain rate tensor, plays an important role in the flow classification, in vortical structure identification and in the modeling of various complex fluids and flows. For the simulation of flows accompanied by the advection-diffusion transport of a scalar field (e.g., temperature), double distribution functions (DDF) based lattice Boltzmann (LB) methods, involving a pair of LB schemes are commonly used. We present a new local vorticity computation approach by introducing an intensional anisotropy of the scalar flux in the third order, off-diagonal moment equilibria of the LB scheme for the scalar field, and then combining the second order non-equilibrium components of both the LB methods. As such, any pair of lattice sets in the DDF formulation that can independently support the third order off-diagonal moments would enable local determination of the complete flow kinematics, with the LB methods for the fluid motion and the transport of the passive scalar respectively providing the necessary moment relationships to determine the symmetric and skew-symmetric components of the velocity gradient tensor. Since the resulting formulation is completely local and do not rely on any finite difference approximations for velocity derivatives, it is by design naturally suitable for parallel computation. As an illustration of our approach, we formulate a DDF-LB scheme for local vorticity computation using a pair of multiple relaxation times (MRT) based collision approaches on two-dimensional, nine velocity (D2Q9) lattices, where the necessary moment relationships to determine the velocity gradient tensor and the vorticity are established via a Chapman-Enskog analysis. Simulations of various benchmark flows demonstrate good accuracy of the predicted vorticity fields using our approach against available solutions, including numerical results, with a second order convergence. Furthermore, extensions of our formulation for a variety of collision models, including those based on cascaded and non-cascaded central moments, to enable local vorticity computation are presented.     

Application of the lattice‐Boltzmann method to study flow and dispersion in channels with and without expansion and contraction geometry

The lattice‐Boltzmann (LB) method, derived from lattice gas automata, is a relatively new technique for studying transport problems. The LB method is investigated for its accuracy to study fluid dynamics and dispersion problems. Two problems of relevance to flow and dispersion in porous media are addressed: (i) Poiseuille flow between parallel plates (which is analogous to flow in pore throats in two‐dimensional porous networks), and (ii) flow through an expansion–contraction geometry (which is analogous to flow in pore bodies in two‐dimensional porous networks). The results obtained from the LB simulations are compared with analytical solutions when available, and with solutions obtained from a finite element code (FIDAP) when analytical results are not available. Excellent agreement is found between the LB results and the analytical/FIDAP solutions in most cases, indicating the utility of the lattice‐Boltzmann method for solving fluid dynamics and dispersion problems. Copyright © 1999 John Wiley & Sons, Ltd.     

Numerical study of the properties of the central moment lattice Boltzmann method

Central moment lattice Boltzmann method (LBM) is one of the more recent developments among the lattice kinetic schemes for computational fluid dynamics. A key element in this approach is the use of central moments to specify the collision process and forcing, and thereby naturally maintaining Galilean invariance, an important characteristic of fluid flows. When the different central moments are relaxed at different rates like in a standard multiple relaxation time (MRT) formulation based on raw moments, it is endowed with a number of desirable physical and numerical features. Because the collision operator exhibits a cascaded structure, this approach is also known as the cascaded LBM. While the cascaded LBM has been developed sometime ago, a systematic study of its numerical properties, such as the accuracy, grid convergence, and stability for well‐defined canonical problems is lacking, and the present work is intended to fulfill this need. We perform a quantitative study of the performance of the cascaded LBM for a set of benchmark problems of differing complexity, viz., Poiseuille flow, decaying Taylor–Green vortex flow, and lid‐driven cavity flow. We first establish its grid convergence and demonstrate second‐order accuracy under diffusive scaling for both the velocity field and its derivatives, that is, the components of the strain rate tensor, as well. The method is shown to quantitatively reproduce steady/unsteady analytical solutions or other numerical results with excellent accuracy. The cascaded MRT LBM based on the central moments is found to be of similar accuracy when compared with the standard MRT LBM based on the raw moments, when a detailed comparison of the flow fields are made, with both reproducing even the small scale vortical features well. Numerical experiments further demonstrate that the central moment MRT LBM results in significant stability improvements when compared with certain existing collision models at moderate additional computational cost. Copyright © 2015 John Wiley & Sons, Ltd.     

Unsteady three-dimensional MHD flow of the micropolar fluid over an oscillatory disk with Cattaneo-Christov double diffusion

Cattaneo-Christov heat and mass flux models are considered rather than Fourier and Fick laws due to the presence of thermal and concentration transport hyperbolic phenomena. The generalized form of the Navier-Stokes model is considered in hydromagnetic flow. Three-dimensional(3 D) unsteady fluid motion is generated by the periodic oscillations of a rotating disk. Similarity transformations are used to obtain the normalized fluid flow model. The successive over relaxation(SOR) method with finite difference schemes are accomplished for the numerical solution of the obtained partial differential non-linear system. The flow features of the velocity, microrotation, temperature,and concentration fields are discussed in pictorial forms for various physical flow parameters. The couple stresses and heat and mass transfer rates for different physical quantities are explained via tabular forms. For better insight of the physical fluid model, 3 D fluid phenomena and two-dimensional(2 D) contours are also plotted. The results show that the micropolar fluids contain microstructure having non-symmetric stress tensor and are useful in lubrication theory. Moreover, the thermal and concentration waves in CattaneoChristov models have a significance role in the laser heating and enhancement in thermal conductivity.     

An adaptive finite element method for viscoplastic flows in a square pipe with stick–slip at the wall

This paper presents the numerical solution of non-linear yield stress phenomena by using a new mixed anisotropic auto-adaptive finite element method. The Poiseuille flow of a Bingham fluid with slip yield boundary condition at the wall is considered. Despite its practical interest, for instance for pipeline flows of yield-stress fluids such as concrete and cements, this problem has not been addressed yet to our knowledge. The case of a pipe with a square section has been investigated in detail. The computations cover the full range of the two main dimensionless numbers and exhibit complex flow patterns: all the different flow regimes are completely identified.     

Lattice Boltzmann methods for shallow water flow applications

We apply the lattice Boltzmann (LB) method for solving the shallow water equations with source terms such as the bed slope and bed friction. Our aim is to use a simple and accurate representation of the source terms in order to simulate practical shallow water flows without relying on upwind discretization or Riemann problem solvers. We validate the algorithm on problems where analytical solutions are available. The numerical results are in good agreement with analytical solutions. Furthermore, we test the method on a practical problem by simulating mean flow in the Strait of Gibraltar. The main focus is to examine the performance of the LB method for complex geometries with irregular bathymetry. The results demonstrate its ability to capture the main flow features. Copyright © 2007 John Wiley & Sons, Ltd.     

Simple,coupled algorithms for solving creeping flows and their application to electro‐osmotic flows

In this study we developed simple, coupled algorithms for solving low‐Reynolds‐number flows applicable to micro‐scale flows such as electro‐osmotic flows. The most popular scheme, i.e. the projection method, is not suitable for such flows because of its undesirable slip effect on boundaries at low‐Reynolds‐numbers. In our method, the velocity and pressure are strongly coupled, and the momentum and pressure equations are solved iteratively by using the successive over relaxation (SOR) method while exchanging the unknown variables as soon as they have been updated. The developed methods are applied to a model flow for evaluating their performance. It was found that the coupled schemes are indeed superior to a projection method, i.e. the fractional‐step method, in both numerical accuracy and CPU time. The code is then applied to a dc electro‐osmotic flow within a cavity driven by electrical force acting on the ions spread in the fluid. In this application, the system of equations for the fluid flow and that for the ion transport are solved in a decoupled way, but each system is solved by using fully implicit schemes. From the simulations and by introducing the concept of vorticity source, we can identify two roles of the body force, one contributing to build‐up of the osmotic pressure and the other to the fluid flow. The interesting reverse flow occurring after the external potentials applied on the electrodes have been shut off is also investigated in terms of the vorticity source. Copyright © 2010 John Wiley & Sons, Ltd.     

Lie symmetry solutions for two-phase mass flows

We apply Lie symmetry method to a set of non-linear partial differential equations, which describes a two-phase rapid gravity mass flow as a mixture of solid particles and viscous fluid down a slope (Pudasaini, J. Geophys. Res. 117 (2012) F03010, 28 pp [1]). In order to systematically explore the mathematical structure and underlying physics of the two-phase mixture flow, we generate several similarity forms in general form and construct self-similar solutions. Our analysis generalizes the results, obtained by applying the Lie symmetry method to relatively simple single-phase pressure-driven gravity mass flows, to the two-phase mass flows that include several dominant driving forces and strong phase-interactions. Analytical and numerical solutions are presented for the symmetry-reduced homogeneous and non-homogeneous systems of equations. Analytical and numerical results show that the new models presented here can adequately describe the dynamics of two-phase debris flows, and produce observable phenomena that are consistent with the physics of the flow. The solutions are strongly dependent on the choice of the symmetry-reduced model, as characterized by the group parameters, and the physical parameters of the flows. These solutions reveal strong non-linear and distinct dynamic evolutions, and phase-interactions between the solid and fluid phases, namely the phase-heights and phase-velocities.     

基于Boltzmann模型方程的气体运动论统一算法研究

模型方程出发,研究确立含流态控制参数可描述不同流域气体流动特征的气体分子速度分布函数方程; 研究发展气体运动论离散速度坐标法, 借助非定常时间分裂数值计算方法和NND差分格式, 结合DSMC方法关于分子运动与碰撞去耦技术, 发展直接求解速度分布函数的气体运动论耦合迭代数值格式; 研制可用于物理空间各点宏观流动取矩的离散速度数值积分方法, 由此提出一套能有效模拟稀薄流到连续流不同流域气体流动问题统一算法. 通过对不同Knudsen数下一维激波内流动、二维圆柱、三维球体绕流数值计算表明, 计算结果与有关实验数据及其它途径研究结果(如DSMC模拟值、N-S数值解)吻合较好, 证实气体运动论统一算法求解各流域气体流动问题的可行性. 尝试将统一算法进行HPF并行化程序设计, 基于对球体绕流及类``神舟'返回舱外形绕流问题进行HPF初步并行试算, 显示出统一算法具有很好的并行可扩展性, 可望建立起新型的能有效模拟各流域飞行器绕流HPF并行算法研究方向. 通过将气体运动论统一算法推广应用于微槽道流动计算研究, 已初步发展起可靠模拟二维短微槽道流动数值算法; 通过对Couette流、Poiseuille流、压力驱动的二维短槽道流数值模拟, 证实该算法对微槽道气体流动问题具有较强的模拟能力, 可望发展起基于Boltzmann模型方程能可靠模拟MEMS微流动问题气体运动论数值计算方法研究途径.      

Laminar flow of micropolar fluid in rectangular microchannels

Compared with the classic flow on macroscale, flows in microchannels have some new phenomena such as the friction increase and the flow rate reduction. Papautsky and co-workers explained these phenomena by using a micropolar fluid model where the effects of micro-rotation of fluid molecules were taken into account. But both the curl of velocity vector and the curl of micro-rotation gyration vector were given incorrectly in the Cartesian coordinates and then the micro-rotation gyration vector had only one component in the z-direction. Besides, the gradient term of the divergence of micro-rotation gyration vector was missed improperly in the angular moment equation. In this paper, the governing equations for laminar flows of micropolar fluid in rectangular microchannels are reconstructed. The numerical results of velocity profiles and micro-rotation gyrations are obtained by a procedure based on the Chebyshev collocation method. The micropolar effects on velocity and micro-rotation gyration are discussed in detail.The project was supported by the National Natural Science Foundation of China (10472054). The English text was polished by Boyi Wang.     

A stabilized Nitsche‐type extended embedding mesh approach for 3D low‐ and high‐Reynolds‐number flows

This paper presents a stabilized extended finite element method (XFEM) based fluid formulation to embed arbitrary fluid patches into a fixed background fluid mesh. The new approach is highly beneficial when it comes to computational grid generation for complex domains, as it allows locally increased resolutions independent from size and structure of the background mesh. Motivating applications for such a domain decomposition technique are complex fluid‐structure interaction problems, where an additional boundary layer mesh is used to accurately capture the flow around the structure. The objective of this work is to provide an accurate and robust XFEM‐based coupling for low‐ as well as high‐Reynolds‐number flows. Our formulation is built from the following essential ingredients: Coupling conditions on the embedded interface are imposed weakly using Nitsche's method supported by extra terms to guarantee mass conservation and to control the convective mass transport across the interface for transient viscous‐dominated and convection‐dominated flows. Residual‐based fluid stabilizations in the interior of the fluid subdomains and accompanying face‐oriented fluid and ghost‐penalty stabilizations in the interface zone stabilize the formulation in the entire fluid domain. A detailed numerical study of our stabilized embedded fluid formulation, including an investigation of variants of Nitsche's method for viscous flows, shows optimal error convergence for viscous‐dominated and convection‐dominated flow problems independent of the interface position. Challenging two‐dimensional and three‐dimensional numerical examples highlight the robustness of our approach in all flow regimes: benchmark computations for laminar flow around a cylinder, a turbulent driven cavity flow at Re = 10000 and the flow interacting with a three‐dimensional flexible wall. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical modelling of 3D stratified free surface flows: a case study of sediment dumping

A three‐dimensional numerical model has been developed to simulate stratified flows with free surfaces. The model is based on the Reynolds‐averaged Navier–Stokes (RANS) equations with variable fluid density. The equations are solved in a transformed σ‐coordinate system with the use of operator‐splitting method (Int. J. Numer. Meth. Fluids 2002; 38 :1045–1068). The numerical model is validated against the one‐dimensional diffusion problem and the two‐dimensional density‐gradient flow. Excellent agreements are obtained between numerical results and analytical solutions. The model is then used to study transport phenomena of dumped sediments into a water body, which has been modelled as a strongly stratified flow. For the two‐dimensional problem, the numerical results compare well with experimental data in terms of mean particle falling velocity and spreading rate of the sediment cloud for both coarse and medium‐size sediments. The model is also employed to study the dumping of sediments in a three‐dimensional environment with the presence of free surface. It is found that during the descending process an annulus‐like cloud is formed for fine sediments whereas a plate‐like cloud for medium‐size sediments. The model is proven to be a good tool to simulate strongly stratified free surface flows. Copyright © 2005 John Wiley & Sons, Ltd.     

A coupled lattice Boltzmann model for the shallow water‐contamination system

A coupled numerical method for the direct simulation of shallow water dynamics and pollutant transport is formulated and implemented. The conservation equations of shallow water dynamics equations and the convection–diffusion equations are solved using the lattice Boltzmann (LB) method. The local equilibrium distribution of the pollutant has no terms of second order in flow velocity. And the relaxation time of the pollutant deviates from a constant for the flows with variable free surface water depth. The numerical tests show that this scheme strictly obeys the conservation law of mass and momentum. Excellent agreement is obtained between numerical predictions and analytical solutions in the pure diffusion problem and convection–diffusion problem. Furthermore, the influences on the accuracy of the lattice size and the diffusivity are also studied. The results indicate that the variation in the free surface water depth cannot affect the conservation of the model, and the model has the ability to simulate the complex topography problem. The comparison shows that the LB scheme has the capacity to solve the complex convection–diffusion problem in shallow water. Copyright © 2008 John Wiley & Sons, Ltd.     

Transport and deposition of neutral particles in magnetohydrodynamic turbulent channel flows at low magnetic Reynolds numbers

The effect of Lorentz force on particle transport and deposition is studied by using direct numerical simulation of turbulent channel flow of electrically conducting fluids combined with discrete particle simulation of the trajectories of uncharged, spherical particles. The magnetohydrodynamic equations for fluid flows at low magnetic Reynolds numbers are adopted. The particle motion is determined by the drag, added mass, and pressure gradient forces. Results are obtained for flows with particle ensembles of various densities and diameters in the presence of streamwise, wall-normal or spanwise magnetic fields. It is found that the particle dispersion in the wall-normal and spanwise directions is decreased due to the changes of the underlying fluid turbulence by the Lorentz force, while it is increased in the streamwise direction. The particle accumulation in the near-wall region is diminished in the magnetohydrodynamic flows. In addition, the tendency of small inertia particles to concentrate preferentially in the low-speed streaks near the walls is strengthened with increasing Hartmann number. The particle transport by turbophoretic drift and turbulent diffusion is damped by the magnetic field and, consequently, particle deposition is reduced.     

Interaction of a deformable solid with two-phase flows: An Eulerian-based numerical model for fluid-structure interaction using the level contour reconstruction method

We describe the formulation of a method for fluid-structure interaction involving the coupling of moving and/or flexible solid structures with multiphase flows in the framework of the Level Contour Reconstruction Method. We present an Eulerian-based numerical procedure for tracking the motion and interaction of a liquid-gas interface with a fluid-solid interface in the Lagrangian frame together with the evaluation of the fluid transport equations coupled to those for the solid transport, namely the left Cauchy-Green strain tensor field, in the Eulerian frame. To prevent excessive dissipation due to the convective nature of the solid transport equation, a simple incompressibility constraint for the strain field is enforced. A single grid structure is used for both the fluid and solid phases which allows for a simple and natural coupling of the fluid and solid dynamics. Several benchmark tests are performed to show the accuracy of the numerical method and which demonstrate accurate results compared to several of those in the existing literature. In particular we show that surface tension effects including contact line dynamics on the deforming solid phase can be properly simulated. The three-phase interaction of a droplet impacting on a flexible cantilever is investigated in detail. The simulations follow the detailed motion of the droplet impact (and subsequent deformation, breakup, and fall trajectory) along with the motion of the deformable solid cantilever due to its own weight as well as due to the force of the droplet impact.     

Numerical modeling and simulation of PEM fuel cells: Progress and perspective

This paper provides a comprehensive review on the research and development in multi-scale numerical modeling and simulation of PEM fuel cells. An overview of recent progress in PEM fuel cell modeling has been provided. Fundamental transport phenomena in PEM fuel cells and the corresponding mathematical formulation of macroscale models are analyzed. Various important issues in PEM fuel cell modeling and simulation are examined in detail, including fluid flow and species transport, electron and proton transport, heat transfer and thermal management, liquid water transport and water management, transient response behaviors, and cold-start processes. Key areas for further improvements have also been discussed.     

Numerical simulation of cavitating flow in 2D and 3D inducer geometries

A computational method is proposed to simulate 3D unsteady cavitating flows in spatial turbopump inducers. It is based on the code FineTurbo, adapted to take into account two‐phase flow phenomena. The initial model is a time‐marching algorithm devoted to compressible flow, associated with a low‐speed preconditioner to treat low Mach number flows. The presented work covers the 3D implementation of a physical model developed in LEGI for several years to simulate 2D unsteady cavitating flows. It is based on a barotropic state law that relates the fluid density to the pressure variations. A modification of the preconditioner is proposed to treat efficiently as well highly compressible two‐phase flow areas as weakly compressible single‐phase flow conditions. The numerical model is applied to time‐accurate simulations of cavitating flow in spatial turbopump inducers. The first geometry is a 2D Venturi type section designed to simulate an inducer blade suction side. Results obtained with this simple test case, including the study of its general cavitating behaviour, numerical tests, and precise comparisons with previous experimental measurements inside the cavity, lead to a satisfactory validation of the model. A complete three‐dimensional rotating inducer geometry is then considered, and its quasi‐static behaviour in cavitating conditions is investigated. Numerical results are compared to experimental measurements and visualizations, and a promising agreement is obtained. Copyright © 2004 John Wiley & Sons, Ltd.     

饱和-非饱和土壤中污染物运移过程的数值模拟

本文提出了一个模拟饱和 非饱和土壤中溶和污染物运移过程的数值模型．模拟的控制污染物运移的物理 化学现象包括：对流，机械逸散，分子弥散，吸附，蜕变，不动水效应．发展了一个修正的特征线Ｇａｌｅｒｋｉｎ方法以离散污染物运移过程的控制方程并导出了一个用于有限元方程求解的显式算法．数值例题结果表明所提出模型和算法的功能