Attractors for the motion of finite-size particles in a two-sided lid-driven cavity

Particle trajectories in an incompressible cavity flow are investigated experimentally. The flow is driven by two facing walls moving with same speed in opposite directions. When the flow is steady and three-dimensional, finite-size and neutrally-buoyant particles are rapidly attracted to periodic orbits. Three-dimensional particle tracking is used to reconstruct single particle trajectories. Experimental results are presented and compared with numerical simulations. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

螺旋型旋风分离器两相流场的数值模拟

对螺旋型旋风分离器进行了两相流场的三维数值模拟．气体流场通过求解三维N-S方程得到,湍流模型采用了雷诺应力模型．计算结果表明,旋风分离器内部的流场分为两部分:螺旋通道内比较稳定的流场和筒体中心区域的复合涡结构流场．对颗粒运动轨迹的计算表明,颗粒在入口处的初始位置对颗粒分离有比较显著的影响．同时得到了不同入口速度下颗粒的分级效率曲线,并给出了气体流量对旋风分离器性能的影响,结果显示:气体流量的增加会提高分离效率,但同时导致压力损失的急剧增加．     

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.     

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)     

Particle motion in a spatially periodic flow

The motion of small spherical particles in two-dimensional Taylor–Green flow is investigated assuming one-way coupling. The particle dynamics is modeled by a modified Maxey–Riley equation. Trajectories are computed numerically using a Cash–Karp Runge–Kutta scheme which provides an adaptive time-step control. Depending on the Stokes number and the fluid-to-particle density ratio the particle trajectories can either be periodic or chaotic. To clarify the nature of this dynamical system further investigations are carried out by computing the Feigenbaum scenario. It is expected that the particle dynamics in this idealized flow will provide useful informations for an explanation of particle accumulation phenomena in more realistic periodic three-dimensional viscous flows. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

圆孔外Stokes流中球形粒子的运动*

在研究悬浮液入孔问题时,粒子在孔外任意位置所受的力和力矩是所需的最基本数据.本文在严宗毅等(1987)算出的有限个离散数据的基础上.首次给出了球在圆孔外的全部十二个力和力矩系数的近似解析表达式.我们应用这些系数计算零雷诺数下球形粒子入孔时的轨道和旋转角速度,结果与现有实验数据完全一致.分析不同系数的相对重要性表明,在孔口和孔壁附近不能忽略旋转效应.在靠近孔口边缘的局部区域必须计及侧向力效应.以往的理论结果在孔口附近与实验不符,正是由于忽略了这些因素.本文还详细讨论了粒子和孔口尺寸的相对大小以及重力和浮力对于粒子运动轨道、速度分布和旋转的影响,指出用大粒子做中性悬浮实验时,对其密度的要求特别严格.本文所提供的力和力矩系数,考虑因素全面、比较准确,便于计算,为进一步研究各种涉及粒子入口的问题提供了良好的前提.     

Modeling of local particle accumulation in disperse flows with kinematic singularities

The aim of the study is to model the formation of local particle accumulation zones near several typical kinematic singularities. The flows considered are: (i) a steady two–dimensional flow with localized vorticity of the Kelvin cat's eye type (vortex in a mixing layer), (ii) a steady axisymmetric flow formed by a vortex filament normal to a plane in viscous fluid (simple model of tornado), (iii) a neighbourhood of a zero acceleration point in two–dimensional unsteady (harmonic) flow. From parametric numerical calculations, we investigated the inertial mechanisms of forming local particle accumulation zones and found the threshold values of governing parameters separating qualitatively different particle velocity and density patterns. In particular, it is shown that the zero–acceleration point can either “attract” or “scatter” the particles. Zones of concentrated vorticity are typically devoid of particles. In the tornado–like flow, an axisymmetric “cup-shaped” particle accumulation region is formed. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A Meshless Boundary Method for Computation of Stokes Flow with Particles in the Semicircular Canals of the Inner Ear

Benign paroxismal positional vertigo (BPPV) is modeled by introducing free-floating particles (canaliths) which settle inside the semicircular canals (SCC). The Stokes flow induced by a canalith is evaluated by coupling the force coupling method (FCM) to the method of fundamental solutions (MFS). The proposed methodology results in a straightforward meshless boundary method for the simulation of bounded Stokes flow with finite-size particles. (© 2010 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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.     

Numerical study of gas–solid flow in a cyclone separator

This paper presents a numerical study of the gas–powder flow in a typical Lapple cyclone. The turbulence of gas flow is obtained by the use of the Reynolds stress model. The resulting pressure and flow fields are verified by comparing with those measured and then used in the determination of powder flow that is simulated by the use of a stochastic Lagrangian model. The separation efficiency and trajectory of particles from simulation are shown to be comparable to those observed experimentally. The effects of particle size and gas velocity on separation efficiency are quantified and the results agree well with experiments. Some factors which affect the performance of cyclone were identified. It is shown that the collision between gas streams after running about a circle and that just entering occurred around the junction of the inlet duct and the cylinder of the cyclone, resulting in a short-circuiting flow. The combination of flow source and sink was distributed near the axis of cyclone, forming a flow dipole at axial section. Particles entering at different positions gave different separation efficiency. A particle with size exceeding a critical diameter, which was condition-dependant, would stagnate on the wall of cyclone cone. This was regarded as one of the main reasons for the deposition on the inner conical surface in such cyclones used in the cement industry.     

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

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

柱状固粒两相边界层流场稳定性研究

从运动方程和本构方程出发,推导得到了含柱状粒子两相流场的修正Orr-Sommerfeld方程,然后在边界层流场中,采用数值计算方法,得到了含柱状粒子流场的稳定性中性曲线,给出了流场失稳的临界雷诺数.结果表明在所述情况下,柱状粒子对流场起着抑制失稳的作用,而且抑制的程度随着柱状粒子体积分数和长径比的增加而提高.     

Hybrid Simulation of Space Plasmas: Models with Massless Fluid Representation of Electrons. IV. Kelvin–Helmholtz Instability

This is a survey of the literature on hybrid simulation of the Kelvin–Helmholtz instability. We start with a brief review of the theory: the simplest model of the instability—a transition layer in the form of a tangential discontinuity; compressibility of the medium; finite size of the velocity shear region; pressure anisotropy. We then describe the electromagnetic hybrid model (ions as particles and electrons as a massless fluid) and the main numerical schemes. We review the studies on two-dimensional and three-dimensional hybrid simulation of the process of particle mixing across the magnetopause shear layer driven by the onset of a Kelvin–Helmholtz instability. The article concludes with a survey of literature on hybrid simulation of the Kelvin–Helmholtz instability in finite-size objects: jets moving across the magnetic field in the middle of the field reversal layer; interaction between a magnetized plasma flow and a cylindrical plasma source with zero own magnetic field.     

A Lagrangian stochastic model for nonpassive particle diffusion in turbulent flows

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

Existence of an infinite particle limit of stochastic ranking process

We study a stochastic particle system which models the time evolution of the ranking of books by online bookstores (e.g., Amazon.co.jp). In this system, particles are lined in a queue. Each particle jumps at random jump times to the top of the queue, and otherwise stays in the queue, being pushed toward the tail every time another particle jumps to the top. In an infinite particle limit, the random motion of each particle between its jumps converges to a deterministic trajectory. (This trajectory is actually observed in the ranking data on web sites.) We prove that the (random) empirical distribution of this particle system converges to a deterministic space–time-dependent distribution. A core of the proof is the law of large numbers for dependent random variables.     

Numerical Investigation of Particulate Flow over a Backward‐Facing Step

Flow separation and recirculation caused by a sudden expansion in the channel geometry in the form of a backwardfacing step (BFS) appear in numerous practical applications. Additionally, BFS flow has been used as a generic test case to study fundamental flow properties, such as separation or re‐attachment. In the present work, BFS flow laden with dispersed particles is investigated by numerical simulations using a spectral element method [1]. The motion of the dispersed particles is computed by Lagrangian particle tracking. In a first step, only the influence of the flow on the particles is accounted for, while possible effects of the particle motion on the flow are neglected. Spatial distribution of the particles is investigated, and effects of different wall‐particle interaction models on the computational results are examined.     

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 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.