Numerical simulation of the performance of a snow fence with airfoil snow plates by FVM

The aim of this work is to analyze the efficiency of a snow fence with airfoil snow plates to avoid the snowdrift formation, to improve visibility and to prevent blowing snow disasters on highways and railways. In order to attain this objective, it is necessary to solve particle transport equations along with the turbulent fluid flow equations since there are two phases: solid phase (snow particles) and fluid phase (air). In the first place, the turbulent flow is modelled by solving the Reynolds-averaged Navier-Stokes (RANS) equations for incompressible viscous flows through the finite volume method (FVM) and then, once the flow velocity field has been determined, representative particles are tracked using the Lagrangian approach. Within the particle transport models, we have used a particle transport model termed as Lagrangian particle tracking model, where particulates are tracked through the flow in a Lagrangian way. The full particulate phase is modelled by just a sample of about 15,000 individual particles. The tracking is carried out by forming a set of ordinary differential equations in time for each particle, consisting of equations for position and velocity. These equations are then integrated using a simple integration method to calculate the behaviour of the particles as they traverse the flow domain. Finally, the conclusions of this work are exposed.     

Steady state numerical simulation of the particle collection efficiency of a new urban sustainable gravity settler using design of experiments by FVM

The aim of this work is to analyze the efficiency of a new sustainable urban gravity settler to avoid the solid particle transport, to improve the water waste quality and to prevent pollution problems due to rain water harvesting in areas with no drainage pavement. In order to get this objective, it is necessary to solve particle transport equations along with the turbulent fluid flow equations since there are two phases: solid phase (sand particles) and fluid phase (water). In the first place, the turbulent flow is modelled by solving the Reynolds-averaged Navier-Stokes (RANS) equations for incompressible viscous flows through the finite volume method (FVM) and then, once the flow velocity field has been determined, representative particles are tracked using the Lagrangian approach. Within the particle transport models, a particle transport model termed as Lagrangian particle tracking model is used, where particulates are tracked through the flow in a Lagrangian way. The full particulate phase is modelled by just a sample of about 2,000 individual particles. The tracking is carried out by forming a set of ordinary differential equations in time for each particle, consisting of equations for position and velocity. These equations are then integrated using a simple integration method to calculate the behaviour of the particles as they traverse the flow domain. The entire FVM model is built and the design of experiments (DOE) method was used to limit the number of simulations required, saving on the computational time significantly needed to arrive at the optimum configuration of the settler. Finally, conclusions of this work are exposed.     

A numerical study on statistical temporal scales in inertia particle dispersion

The statistical temporal scales involved in inertia particle dispersion are analyzed numerically. The numerical method of large eddy simulation, solving a filtered Navier-Stokes equation, is utilized to calculate fully developed turbulent channel flows with Reynolds numbers of 180 and 640, and the particle Lagrangian trajectory method is employed to track inertia particles released into the flow fields. The Lagrangian and Eulerian temporal scales are obtained statistically for fluid tracer particles and three different inertia particles with Stokes numbers of 1, 10 and 100. The Eulerian temporal scales, decreasing with the velocity of advection from the wall to the channel central plane, are smaller than the Lagrangian ones. The Lagrangian temporal scales of inertia particles increase with the particle Stokes number. The Lagrangian temporal scales of the fluid phase ‘seen’ by inertia particles are separate from those of the fluid phase, where inertia particles travel in turbulent vortices, due to the particle inertia and particle trajectory crossing effects. The effects of the Reynolds number on the integral temporal scales are also discussed. The results are worthy of use in examining and developing engineering prediction models of particle dispersion.     

Simulation of free turbulent particle-laden jet using Reynolds-stress gas turbulence model

Free two-phase flows occur in many practical applications, such as sprays or particle drying and combustion. This paper deals with mathematical modelling of a free turbulent two-phase jet. A steady, axisymmetric, dilute, monodisperse, particle-laden, turbulent jet injected into a still environment, has been considered. The model treats the gas-phase from an Eulerian standpoint and the motion of particles from a Lagrangian one. Closure of the system of time averaged transport equations has been accomplished by using a Reynolds-stress turbulence model. The particles–fluid interaction has been considered by the PSI-Cell concept. Both the effect of interphase slip and the effect of particle dispersion have been taken into account.     

Influence of Mass Loading on Particle-Laden Turbulent Pipe Flow

A CFD code in the framework of OpenFOAM was validated for simulations of particle-laden pipe and channel flows at low to intermediate mass loadings. The code is based on an Eulerian two-fluid approach with Reynolds-averaged conservation equations, including turbulence modeling and four-way coupling. Pipe flow simulations of particles in air against gravity were conducted at Reynolds numbers up to 50000. The particle mass loading was varied and its effect on the mean velocities and turbulent fluctuations of the two phases was studied. Special attention was paid to the influence of mass loading on the centerline velocity and the wall shear velocity of the fluid phase for various flow parameters and particle properties. Empirical correlations were established between these two quantities and the flow Reynolds number, particle Reynolds number, Stokes number and particle to fluid density ratio for a range of particle mass loadings. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An Eulerian–Lagrangian mixed discrete least squares meshfree method for incompressible multiphase flow problems

Mixed discrete least squares meshfree (MDLSM) method has been developed as a truly meshfree method and successfully used to solve single-phase flow problems. In the MDLSM, a residual functional is minimized in terms of the nodal unknown parameters leading to a set of positive-definite system of algebraic equations. The functional is defined using a least square summation of the residual of the governing partial differential equations and its boundary conditions at all nodal points discretizing the computational domain. Unlike the discrete least squares meshfree (DLSM) which uses an irreducible form of the governing equations, the MDLSM uses a mixed form of the original governing equations allowing for direct calculation of the gradients leading to more accurate computational results. In this study, an Eulerian–Lagrangian MDLSM method is proposed to solve incompressible multiphase flow problems. In the Eulerian step, the MDLSM method is used to solve the governing phase averaged Navier–Stokes equations discretized at fixed nodal points to get the velocity and pressure fields. A Lagrangian based approach is then used to track different flow phases indexed by a set of marker points. The velocities of marker points are calculated by interpolating the velocity of fixed nodal points using a kernel approximation, which are then used to move the marker points as Lagrangian particles to track phases. To avoid unphysical clustering and dispersing of the marker points, as a common drawback of Lagrangian point tracking methods, a new approach is proposed to smooth the distribution of marker points. The hybrid Eulerian and Lagrangian characteristics of the approach used here provides clear advantages for the proposed method. Since the nodal points are static on the Eulerian step, the time-consuming moving least squares (MLS) approximation is implemented only once making the proposed method more efficient than corresponding fully Lagrangian methods. Furthermore, phases can be simply tracked using the Lagrangian phase tracking procedure. Efficiency of the proposed MDLSM multiphase method is evaluated using several benchmark problems and the results are presented and discussed. The results verify the efficiency and accuracy of the proposed method for solving multiphase flow problems.     

密相两相流基本方程组的封闭性和分散相在壁面处的剪应力

本文对将连续相和分散相分别作连续介质考虑的两相流基本方程组,通过分析密相下两相流动的物理机制,提出了分散相的应力满足宾汉流体的本构方程,并用分散相在壁面处剪应力作为分散相的一个边界条件,从而得到了一个封闭的处理密相系统的两相流数学模型.最后,还讨论了分散相在壁面处剪切力的表达式.     

宾汉流体稠密两相湍流流动中的二阶矩模型

建立的Bingham流体稠密两相流动的二阶矩-颗粒动力论湍流模型（USM-theta模型）既体现了两相的作用,又体现了屈服应力所引起的附加项,并提出了USM-theta模型下考虑浓度修正值影响的两相湍流流动的算法．利用该模型对圆管内Bingham流体的单相湍流流动、稠密液固两相的湍流流动进行了计算,并和五方程湍流模型进行了比较,结果表明该模型的预测效果更好．利用USM-theta模型对含颗粒的Bingham流体的两相湍流流动进行了模拟,随着屈服应力的增加,Bingham流体相与颗粒相在管道中心附近的主流速度减小．液固两相湍流和Bingham流体两相湍流的计算结果表明屈服应力引起的附加项对流动有很重要的影响．     

两相流中柱状固粒对流体湍动特性影响的研究

对含柱状固粒的两相流场,建立了包含柱状固粒对流场影响的流体脉动速度方程,在求解脉动速度方程的基础上,经平均得到流体的湍流强度和雷诺应力.将该方法用于槽流湍流场的求解,并与单相流实验结果进行了比较.计算中变化柱状固粒的参数,给出了固粒的体积分数、长径比、松驰时间对流场湍动特性的影响,说明粒子对流场的湍动特性起着抑制作用,其抑制的程度与粒子的体积分数、长径比成正比,与粒子的松弛时间成反比.     

Modeling of two-phase flows with surface tension by finite pointset method (FPM)

A meshfree method for two-phase immiscible incompressible flows including surface tension is presented. The continuum surface force (CSF) model is used to include the surface tension force. The incompressible Navier–Stokes equation is considered as the mathematical model. Application of implicit projection method results in linear second-order partial differential equations for velocities and pressure. These equations are then solved by the finite pointset method (FPM), which is a meshfree and Lagrangian method. The fluid is represented as finite number of particles and the immiscible fluids are distinguished by the color of each particle. The interface is tracked automatically by advecting the color functions for each particle. Two test cases, Laplace's law and the Rayleigh–Taylor instability in 2D have been presented. The results are found to be consistent with the theoretical results.     

固液两相流中的K-ε双方程湍流模式*

建立了固液两相流中更一般的K-ε双方程湍流模式。模化了固相和液相的连续方程、动量方程及K方程和ε方程。该湍流模型考虑了固液两相间速度的滑移,颗粒间的作用及相间作用。使用本文所建立的湍流模型,数值预测了一管湍流中的沙水混合流动,其预测结果与实验结果比较一致。     

Computation of dilute two-phase flow in a pump

This paper is a report on a joint project between academia and industry which is concerned with computation of dilute two-phase flow through a pump in turbulent condition. The flow field for the continuous phase is computed using the Reynolds averaged Navier–Stokes equations together with mixing length turbulence modeling. The dispersed phase is treated using the Lagrangian approach by tracking it's trajectory along which the information is passed. It is found that the bubbles and small solid particles flow out of the chamber (between the rotating impeller and the casing wall) with the conveying fluid. The solid particles of relatively bigger sizes accumulate at the low pressure zones near the cashing wall or the rotating shaft.     

Dispersion of solid particles in turbulant flow through pipe expansions

Studies of the effects upon gas-liquid two-phase flows of pipefittings such as expansions, contractions, bends, and valveshave usually concentrated upon pressuredrop correlations andhave not attempted to determine changes in the distributionsof the gas and liquid phases caused by the fitting. However,it is known that such information is important if, for example,flow separators, which divide the gas and liquid phases in avariety of industrial processes, are to function efficiently.It is therefore important to gain an understanding of the influenceupon phase distributions of the common pipe fittings mentionedabove, which will be found in almost any industrial pipeworksystem. As a first step, the dispersion of solid particles carriedby turbulent gas flows through a pipe expansion has been modellednumerically. The commercial fluid-flow code CFDS-FLOW{smalltilde}hDas been used to model the gas flow, together with aneddy interaction model for determination of the motion of thesolid particles. Mean particle velocities and root-mean-squarevalues of the particle velocity fluctuations, as well as particleconcentrations, are evaluated and compared with recent experimentalresults. The influence of different eddy-length and eddy-lifetimespecifications upon the dispersion of particles of various sizesis investigated. It is found that the different eddy characteristicshave little effect on predicted mean particle velocities, whereasfluctuations in particle velocities and particle concentrationare sensitive to the changes made. By comparing the resultswith experimental data, it is possible to draw conclusions aboutthe relative merits of the different eddy specifications.     

Eulerian Simulation of Gas‐Solid Particle Flow by BDIM

The numerical scheme based on the boundary domain integral method (BDIM) for the numerical simulation of twophase two‐component flows is presented. A program is being developed to model the hydrodynamics of fluidized bed systems by using the Eulerian approach in terms of velocity‐vorticity variables formulation. With the vorticity vector both phases motion computation scheme is partitioned into its kinematic and kinetic aspect. Influence of the drag coefficient on the two‐phase two‐component flow field is studied on the two‐phase gas‐solid particles vertical channel flow.     

Mathematical modelling of bottom evolution by particle deposition and resuspension

A two-dimensional numerical model for the evolution of a bottom due to particle deposition and resuspension by a fluid flow is here presented. A computational fluid dynamic approach is used to calculate the flow field and a Lagrangian particle tracking technique is applied to solve the dispersed phase. The evolution of the lower boundary is simulated taking into account the mass conservation of the solid phase and the geotechnical properties of the granular material. The model is characterized by two important features. First, fluid dynamics are coupled with the bottom evolution due to particle deposition and resuspension. This permits to use the model to simulate complex flow fields as well as complex time-evolving geometries. Second, the dispersed phase is calculated by a Lagrangian approach, which retains the discrete information of the individual particles of the granular bottom which may be of interest for some industrial processes (coating) and environmental flows (sediment stratification). First consistency checks have been performed for some deposition and resuspension test cases with fluid at rest. The model has also been tested by comparison with a physical experiment of deposition inside a cavity. Finally, as an example of possible applications of industrial and environmental interest, the model has been applied to investigate particle deposition in rectangular cavities and the evolution of a sand heap by a fluid flow.     

垂直管道中密相两相流动的解析解

本文依据文献[1]的密相两相流动的数学模型,对垂直圆管中密相两相流动进行了解析求解,分别得到了连续相和分散相的速度解析表达式.在相间阻力与相间速度差成比例时,除了在离管壁面很近的薄区之外,管道流动规律与达西渗流定律完全一致.本文验证了文献[1]的密相两相流动数学模型的假定在本文情形下是合理的.     

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.     

二相流计算的一种差分算法

考虑两相流的力学行为,忽略相间的耗散作用,建立了Euler型的基本控制方程．状态方程采用刚性状态方程．基于Abgrall提出的准则,在流动区域内,对可压两相流提出了一个精度较高的Euler型数值方法,数值格式是Godunov型格式,对守恒型和非守恒型方程采用HLLC型和Lax-Friedrichs型近似Riemann解算器,引入了速度驰豫和压强驰豫过程来代替两相间的相互作用．在一维情形下给出数值算例,并且和Saurel的算例进行了比较,结果表明该算法不但精确而且稳定,且在间断处没有数值振荡．     

Analysis of the one-dimensional Euler–Lagrange equation of continuum mechanics with a Lagrangian of a special form

The equations describing the flow of a one-dimensional continuum in Lagrangian coordinates are studied in this paper by the group analysis method. They are reduced to a single Euler–Lagrange equation which contains two undetermined functions (arbitrary elements). Particular choices of these arbitrary elements correspond to different forms of the shallow water equations, including those with both, a varying bottom and advective impulse transfer effect, and also some other motions of a continuum. A complete group classification of the equations with respect to the arbitrary elements is performed.One advantage of the Lagrangian coordinates consists of the presence of a Lagrangian, so that the equations studied become Euler–Lagrange equations. This allows us to apply Noether’s theorem for constructing conservation laws in Lagrangian coordinates. Not every conservation law in Lagrangian coordinates has a counterpart in Eulerian coordinates, whereas the converse is true. Using Noether’s theorem, conservation laws which can be obtained by the point symmetries are presented, and their analogs in Eulerian coordinates are given, where they exist.     

Development of boundary transfer method in simulation of gas–solid turbulent flow of a riser

Computational fluid dynamics (CFD) is used to simulate the behavior of two phase gas solid in a fluid catalytic cracking (FCC) riser. Gas and particle phases are considered as separate fully interpenetrating continuous media within each control volume. Each phase described in terms of its own separate mass and momentum conservation equations. Simple k–epsilon (k_g–?_g) turbulence model is used for the gas phase and the solid phase is handled with the kinetic theory of granular flows. Source terms are used to account for the influence of hydrodynamic drag on the production, dissipation and exchange of turbulent kinetic energy between the phases. For the particles partial slip condition is considered at the wall.