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.     

CFD modelling of the flow and reactions in the Olympic Dam flash furnace smelter reaction shaft

The performance of flash furnace burners can be evaluated quickly and efficiently using CFD modelling. Gas flows are modelled using the conventional Eulerian approach, while Lagrangian particle tracking is used to model the flow of solid feed through the burner and into the reaction shaft. A composite particle model has been developed that considers the solid feed to be made up of single particles containing appropriate quantities of concentrate, flux and dust. Solid fuels (such as coal) can also be included in the composite particle. Reactions between the solids and gas are then modelled using standard heat and mass transfer relationships. Results from the modelling process are shown for BHP-Billiton’s Olympic Dam copper flash smelter with the burner that was used from 1998–2003. Flow patterns, temperature and gas composition distributions, particle dispersion and residence time, and overall extent of sulphur removal are predicted and used to evaluate furnace performance. However, results are sensitive to the assumed size of the composite particles, and plant measurements are required to determine the appropriate composite particle size to predict quantitative data.     

A 3-dimensional Eulerian–Eulerian CFD simulation of a hydrocyclone

Hydrocyclones are used in mineral industries for classification and separation of solid particles of different sizes and densities suspended in water medium. In the present study an Eulerian–Eulerian CFD simulation of a solid–liquid hydrocyclone has been carried out taking into account two solid phases and one liquid phase. The average size of the larger particle was 0.6117 and that of the smaller particle was 0.09875 mm. Three separate momentum balance equations for the three phases have been considered unlike that in the mixture model where a single momentum equation is solved for the three phases. Two turbulent models i.e. the Reynolds stress model (RSM) and the standard k–ε model were studied. Comparison of the two turbulence models showed slight variation in prediction of the velocity profile and the separation efficiency. The maximum deviation between the two models was observed near the wall where the stress was maximum for larger size particles.     

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.     

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.     

A 3D non-linear k–ε turbulent model for prediction of flow and mass transport in channel with vegetation

The results from a 3D non-linear k–ε turbulence model with vegetation are presented to investigate the flow structure, the velocity distribution and mass transport process in a straight compound open channel and a curved open channel. The 3D numerical model for calculating flow is set up in non-orthogonal curvilinear coordinates in order to calculate the complex boundary channel. The finite volume method is used to disperse the governing equations and the SIMPLEC algorithm is applied to acquire the coupling of velocity and pressure. The non-linear k–ε turbulent model has good useful value because of taking into account the anisotropy and not increasing the computational time. The water level of this model is determined from 2D Poisson equation derived from 2D depth-averaged momentum equations. For concentration simulation, an expression for dispersion through vegetation is derived in the present work for the mixing due to flow over vegetation. The simulated results are in good agreement with available experimental data, which indicates that the developed 3D model can predict the flow structure and mass transport in the open channel with vegetation.     

Eulerian and Lagrangian predictions of particulate two-phase flows: a numerical study

To predict particulate two-phase flows, two approaches are possible. One treats the fluid phase as a continuum and the particulate second phase as single particles. This approach, which predicts the particle trajectories in the fluid phase as a result of forces acting on particles, is called the Lagrangian approach. Treating the solid as some kind of continuum, and solving the appropriate continuum equations for the fluid and particle phases, is referred to as the Eulerian approach.Both approaches are discussed and their basic equations for the particle and fluid phases as well as their numerical treatment are presented. Particular attention is given to the interactions between both phases and their mathematical formulations. The resulting computer codes are discussed.The following cases are presented in detail: vertical pipe flow with various particle concentrations; and sudden expansion in a vertical pipe flow. The results show good agreement between both types of approach.The Lagrangian approach has some advantages for predicting those particulate flows in which large particle accelerations occur. It can also handle particulate two-phase flows consisting of polydispersed particle size distributions. The Eulerian approach seems to have advantages in all flow cases where high particle concentrations occur and where the high void fraction of the flow becomes a dominating flow controlling parameter.     

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.     

Numerical modeling of 3-D flow on porous broad crested weirs

Gabion weirs with optional design as a broad crested weirs are suitable structures to reduce flash flood with a minimal negative impact on the water environment. In the present study, the 3-D flow was simulated around gabion weirs with respect to free-surface water. The Reynolds-averaged Navier–Stokes equations are solved to predict water surface over the gabion weir. The VOF method with the geometric reconstruction scheme was applied to treat the complex free-surface flow. Simulations were performed using three variants of the k–ε and the RSM models to find the water level and velocity distribution profile and results are compared with several experimental data available in the literature. The structured mesh was used for all domains with high dense mesh near the solid region. A comparison between experimental data and simulations indicates that the k–ε model can be used to predict the complex flow and water level with high accuracy.     

Aerodynamic focusing of inertial particles in supersonic micronozzles

A particle-laden flow in a supersonic micronozzle is studied using a one-way coupled two-fluid approach. The carrier gas parameters are obtained from the numerical solution of the Navier-Stokes equations, rarefaction effects are taken into account by imposing velocity slip and temperature jump boundary conditions on the nozzle walls. Under conditions considered, the flow around particles is transitional and free-molecular. As a result of numerical solution of the dispersed-phase equations in Lagrangian variables, two types of particle motion in the expanding part of the nozzle are detected: particle spraying and particle accumulation. The particle focusing effect is most pronounced for particles of about 1–2 µm in size. The particle number density fields contain singularities appearing on the envelopes of particle trajectories. However, the model of non-colliding particles remains valid because the mean distance between the particles near the singularities remains much greater than the particle size. The aerodynamic scheme of aerosol particle focusing proposed may be used in various technologies (microthrusters, needle-free drug injection, microfabrication, etc.). (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

液固两相圆柱绕流尾迹内颗粒扩散分布的离散涡数值研究

基于离散涡方法求得的非定常水流场和颗粒的Lagrange运动方程,数值模拟了稀疏液固两相圆柱绕流尾迹内颗粒的扩散分布．获得了流动的涡谱与3种不同St数颗粒（St=0.25,1.0,40）在流场中的分布．通过引入扩散函数来定量表示颗粒在流场中的纵向扩散强度,并计算得到了不同St数颗粒的扩散函数随时间的变化．数值结果揭示出了液固两相圆柱绕流尾迹中的颗粒扩散分布与颗粒的St数和尾涡结构密切相关:1) 中小St数(St=0.25～4.0)颗粒在运动过程中不能进入涡核区,而在旋涡结构的外沿聚集,且颗粒的St数愈大,其越远离涡核区域;2） 在圆柱绕流尾迹区域内,中小St数(St=0.25～4.0)颗粒的纵向扩散强度随其St数的增大而减小．     

Lagrangian coherent structures,transport and chaotic mixing in simple kinematic ocean models

Methods of dynamical system’s theory are used for numerical study of transport and mixing of passive particles (water masses, temperature, salinity, pollutants, etc.) in simple kinematic ocean models composed with the main Eulerian coherent structures in a randomly fluctuating ocean—a jet-like current and an eddy. Advection of passive tracers in a periodically-driven flow consisting of a background stream and an eddy (the model inspired by the phenomenon of topographic eddies over mountains in the ocean and atmosphere) is analyzed as an example of chaotic particle’s scattering and transport. A numerical analysis reveals a non-attracting chaotic invariant set Λ that determines scattering and trapping of particles from the incoming flow. It is shown that both the trapping time for particles in the mixing region and the number of times their trajectories wind around the vortex have hierarchical fractal structure as functions of the initial particle’s coordinates. Scattering functions are singular on a Cantor set of initial conditions, and this property should manifest itself by strong fluctuations of quantities measured in experiments. The Lagrangian structures in our numerical experiments are shown to be similar to those found in a recent laboratory dye experiment at Woods Hole. Transport and mixing of passive particles is studied in the kinematic model inspired by the interaction of a current (like the Gulf Stream or the Kuroshio) with an eddy in a noisy environment. We demonstrate a non-trivial phenomenon of noise-induced clustering of passive particles and propose a method to find such clusters in numerical experiments. These clusters are patches of advected particles which can move together in a random velocity field for comparatively long time. The clusters appear due to existence of regions of stability in the phase space which is the physical space in the advection problem.     

Numerical study of vortex eccentricity in a gas cyclone

This paper presents a numerical study of the vortex eccentricity in a gas cyclone and its effect on the performance of the cyclone. The gas flow in the cyclone was modeled as an unsteady flow by the Navier–Stokes equations with the Reynolds Stress Model (RSM) as the turbulence model. The particles were modelled by the Lagrangian particle tracking (LPT) approach in an unsteady gas flow. Gas cyclones with the same dimensions and total flow rates but different numbers of inlets were simulated with the inlet velocity varying from 12 to 20 m/s. The vortex eccentricities in different cases were analyzed in terms of radial deviation and angular deviation. In addition, the frequency of the precessing vortex core (PVC) was analyzed by the fast Fourier transform (FFT). The results show that the vortex center in the single inlet cyclone has a great eccentricity and its precession center is also different from the geometric center, which reduces the particle collection efficiency. The increase in the symmetry of the inlet causes only a very small increase in the pressure drop in the simulated cases, but it can significantly reduce the vortex eccentricity, particularly by eliminating the eccentricity of the PVC center. The improvement of the vortex eccentricity can generally increase the collection efficiency for particles greater than 2.0 µm. The increase of the collection efficiency is mainly because the symmetrical vortex can restrain the short-circuiting flow of particles. The results can improve the understanding of the vortex flow in gas cyclones which may guide the optimization of gas cyclones.     

CFD of turbulent transport of particles behind a backward-facing step using a new model—k-ε-Sp

The k-ε-S_p model, describing two-dimensional gas–solid two-phase turbulent flow, has been developed. In this model, the diffusion flux and slip velocity of solid particles are introduced to represent the particle motion in two-phase flow. Based on this model, the gas–solid two-phase turbulent flow behind a vertical backward-facing step is simulated numerically and the turbulent transport velocities of solid particles with high density behind the step are predicted. The numerical simulation is validated by comparing the results of the numerical calculation with two other two-phase turbulent flow models (k-ε-A_p, k-ε-k_p) by Laslandes and the experimental measurements. This model, not only has the same virtues of predicting the longitudinal transport of the solid particles as the present practical two-phase flow models, but also can predict the lateral transport of the solid particles correctly.     

A Lagrangian simulation of particle deposition in a turbulent boundary layer in the presence of thermophoresis

Knowledge of particle deposition in turbulent flows is often required in engineering situations. Examples include fouling of turbine blades, plate-out in nuclear reactors and soot deposition. Thus it is important for numerical simulations to be able to predict particle deposition. Particle deposition is often principally determined by the forces acting on the particles in the boundary layer. The particle tracking facility in the CFD code uses the eddy lifetime model to simulate turbulent particle dispersion, no specific boundary layer being modelled. The particle tracking code has been modified to include a boundary layer. The non-dimensional yplus, y⁺, distance of the particle from the wall is determined and then values for the fluid velocity, fluctuating fluid velocity and eddy lifetime appropriate for a turbulent boundary layer used. Predictions including the boundary layer have been compared against experimental data for particle deposition in turbulent pipe flow. The results giving much better agreement. Many engineering problems also involve heat transfer and hence temperature gradients. Thermophoresis is a phenomena by which small particles experience a force in the opposite direction to the temperature gradient. Thus particles will tend to deposit on cold walls and be repulsed by hot walls. The effect of thermophoresis on the deposition of particles can be significant. The modifications of the particle tracking facility have been extended to include the effect of thermophoresis. A preliminary test case involving the deposition of particles in a heated pipe has been simulated. Comparison with experimental data from an extensive experimental programme undertaken at ISPRA, known as STORM (Simplified Tests on Resuspension Mechanisms), has been made.     

Numerical investigation of blood flow. Part I: In microvessel bifurcations

In some diseases there is a focal pattern of velocity in regions of bifurcation, and thus the dynamics of bifurcation has been investigated in this work. A computational model of blood flow through branching geometries has been used to investigate the influence of bifurcation on blood flow distribution. The flow analysis applies the time-dependent, three-dimensional, incompressible Navier–Stokes equations for Newtonian fluids. The governing equations of mass and momentum conservation were solved to calculate the pressure and velocity fields. Movement of blood flow from an arteriole to a venule via a capillary has been simulated using the volume of fluid (VOF) method. The proposed simulation method would be a useful tool in understanding the hydrodynamics of blood flow where the interaction between the RBC deformation and blood flow movement is important. Discrete particle simulation has been used to simulate the blood flow in a bifurcation with solid and fluid particles. The fluid particle method allows for modeling the plasma as a particle ensemble, where each particle represents a collective unit of fluid, which is defined by its mass, moment of inertia, and translational and angular momenta. These kinds of simulations open a new way for modeling the dynamics of complex, viscoelastic fluids at the micro-scale, where both liquid and solid phases are treated with discrete particles.     

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.     

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.