On model equations for particle dispersion in inhomogeneous turbulence

Comparisons are made between the Advection–Diffusion Equation (ADE) approach for particle transport and the two-fluid model approach based on the PDF method. In principle, the ADE approach offers a much simpler way of calculating the inertial deposition of particles in a turbulent boundary layer than that based on the PDF approach. However the ADE equations that have recently been used are only strictly valid for a simple Gaussian process when particle inertia is small. Using a prescribed, but in general non-Gaussian random particle velocity field, it is shown that the net particle mass flux contains a drift term in addition to that from the mean velocity of the particle velocity field, associated with the compressibility of the velocity field. Furthermore the diffusive flux in general depends not only upon the gradient of the mean concentration (true only for a Gaussian random flow field) but also upon higher order derivatives whose relative contribution depends on diffusion coefficients D_ijk… etc. These coefficients depend upon the statistical moments associated with random displacements and compressibility of the particle flow field along particle trajectories which in turn depend upon particle inertia. In contrast the PDF approach offers the advantage of using a simple gradient (Gaussian) approximation in particle phase space which can lead to a non-Gaussian spatial dispersion process when particle inertia is important. Conditions based on the particle mean free path are derived for which a simple ADE is appropriate. Some of the features of particle transport in an inhomogeneous turbulent flow are illustrated by examining particle dispersion in a random flow field composed of pairs of counter rotating vortices which has an rms velocity which increase linearly from a stagnation point.     

Analysis of high-order velocity moments in a strained channel flow

In the current study, model expressions for fifth-order velocity moments obtained from the truncated Gram-Charlier series expansions model for a turbulent flow field probability density function are validated using data from direct numerical simulation (DNS) of a planar turbulent flow in a strained channel. Simplicity of the model expressions, the lack of unknown coefficients, and their applicability to non-Gaussian turbulent flows make this approach attractive to use for closing turbulent models based on the Reynolds-averaged Navier-Stokes equations. The study confirms validity of the model expressions. It also shows that the imposed flow deformation improves an agreement between the model and DNS profiles for the fifth-order moments in the flow buffer zone including when the flow reverses its direction. The study reveals sensitivity of particularly odd velocity moments to the grid resolution. A new length scale is proposed as a criterion for the grid generation near the wall and in the other flow areas dominated by high mean velocity gradients when higher-order statistics have to be collected from DNS.     

Fractal analysis of a circulating flow field with two different velocity laws

The turbulent shear flow around a rotating cylinder in a quiescent flow is a simple case of a rotating turbulent flow field, where centrifugal force works. Two different power-law mean-velocity distributions exist in this flow field. One is U∝1/r and the other is U∝1/r², where r is a distance from the surface of a cylinder. The behavior of chaos and fractal properties for this complex flow field are investigated. The former concerns the dynamic property in fluid flow and the latter is useful to characterize complex flow patterns or the distribution of turbulence quantities. From the instantaneous velocity signal, we defined the iso-velocity set, and its fractal property was investigated both in the U∝1/r and U∝1/r² region. The instantaneous Reynolds stress is found to be a key factor in this fractal property as conceived for the flat plate boundary layer. The intermittency chaos was applied to investigate the turbulent and non-turbulent distribution in the outer region. A simple one-dimensional model could be useful to identify turbulent and non-turbulent distributions even in this complex flow.     

Numerical study of a swirled turbulent flow in the separation zone of a centrifugal air separator

A mathematical model of a swirled turbulent flow in the separation zone of a pneumatic centrifugal device is presented. The model is based on the known k-? model of turbulence proposed by Wilcox. The influence of rotation of the separation-zone walls, input swirl of the gas flow, and other characteristic parameters on the mean velocity field is demonstrated. A comparison of parameters is performed, which reveals good agreement between the numerical and experimental results for a turbulent fluid flow between parallel disks.     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

Turbulent Flow and Dispersion of Inertial Particles in a Confined Jet Issued by a Long Cylindrical Pipe

In this work we examine first the flow field of a confined jet produced by a turbulent flow in a long cylindrical pipe issuing in an abrupt angle diffuser. Second, we examine the dispersion of inertial micro-particles entrained by the turbulent flow. Specifically, we examine how the particle dispersion field evolves in the multiscale flow generated by the interactions between the large-scale structures, which are geometry dependent, with the smaller turbulent scales issued by the pipe which are advected downstream. We use Large-Eddy-Simulation (LES) for the flow field and Lagrangian tracking for particle dispersion. The complex shape of the domain is modelled using the immersed-boundaries method. Fully developed turbulence inlet conditions are derived from an independent LES of a spatially periodic cylindrical pipe flow. The flow field is analyzed in terms of local velocity signals to determine spatial coherence and decay rate of the coherent K–H vortices and to make quantitative comparisons with experimental data on free jets. Particle dispersion is analyzed in terms of statistical quantities and also with reference to the dynamics of the coherent structures. Results show that the particle dynamics is initially dominated by the Kelvin–Helmholtz (K–H) rolls which form at the expansion and only eventually by the advected smaller turbulence scales.     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

Modeling droplet dispersion and interphase turbulent kinetic energy transfer using a new dual-timescale Langevin model

Dispersion of spray droplets and the modulation of turbulence in the ambient gas by the dispersing droplets are two coupled phenomena that are closely linked to the evolution of global spray characteristics, such as the spreading rate of the spray and the spray cone angle. Direct numerical simulations (DNS) of turbulent gas flows laden with sub-Kolmogorov size particles, in the absence of gravity, report that dispersion statistics and turbulent kinetic energy (TKE) evolve on different timescales. Furthermore, each timescale behaves differently with Stokes number, a non-dimensional flow parameter (defined in this context as the ratio of the particle response time to the Kolmogorov timescale of turbulence) that characterizes how quickly a particle responds to turbulent fluctuations in the carrier or gas phase. A new dual-timescale Langevin model (DLM) composed of two coupled Langevin equations for the fluctuating velocities, one for each phase, is proposed. This model possesses a unique feature that the implied TKE and velocity autocorrelation in each phase evolve on different timescales. Consequently, this model has the capability of simultaneously predicting the disparate Stokes number trends in the evolution of dispersion statistics, such as velocity autocorrelations, and TKE in each phase. Predictions of dispersion statistics and TKE from the new model show good agreement with published DNS of non-evaporating and evaporating droplet-laden turbulent flow.     

Temporal Evolution of Vortical Structures in the Wall Region of Turbulent Channel Flow

Hairpin-like vortical structures that form in the wall region of turbulent channel flow are investigated. The analysis is performed by following a procedure in which the Navier-Stokes equations are first integrated by means of a computational code based on a mixed spectral-finite difference technique in the case of the flow in a plane channel. A DNS turbulent-flow database, representing the turbulent statistically steady state of the velocity field through 10 viscous time units, is computed and the vortex-detection method of the imaginary part of the complex eigenvalue pair of the velocity-gradient tensor is applied to the velocity field. As a result, hairpin-like vortical structures are educed. Flow visualizations are provided of the processes of evolution that characterize hairpin vortices in the wall region of turbulent channel flow. The relationship is investigated between vortex dynamics and 2nd- and 4th- quadrant events, showing that ejections and sweeps play a fundamental role in the way the morphological evolution of a hairpin vortex develops with time.     

Subgrid scale formulation of optical flow for the study of turbulent flow

We propose a new formulation of optical flow in the case of passive scalar or solid particles spreading in turbulent flows. The flow equation is defined from the scalar transport equation and a decomposition of our physical quantities (velocity and concentration fields) into two contributions, large and small scales, to account for the lack of spatial resolution in processed images. A subgrid scale model is introduced to model the small scale contribution. Comparisons are made with existing optical flow methods and Particle Image Velocimetry on synthetic and real sequences. The improvement of the estimation of velocity field by the proposed formulation is discussed in the case of a scalar turbulent propagation.     

入口速度比对射流拟序结构的影响

为了深入了解湍流流动机理以及湍流拟序结构发现过程的影响因素，本文采用大涡模拟方法对不同入口射流伴流速度比的平面湍射流流动进行了数值模拟。采用分步投影法求解动量方程，亚格子项采用标准Smagorinsky亚格子模式模拟，压力泊松方程采用修正的循环消去法快速求解，空间方程采用二阶精度的差分格式，在时间方向上采用二阶精度的显式差分格式。模拟结果给出了平面射流中湍流拟序结构的瞬态发展演变过程，分析了入口速度比对射流拟序结构发展演化过程及宏观流场形态的影响。为进一步研究射流拟序结构及其在湍流流动中的作用提供了基础。     

A fractal analysis of subcooled flow boiling heat transfer

A fractal model for the subcooled flow boiling heat transfer is proposed in this paper. The analytical expressions for the subcooled flow boiling heat transfer are derived based on the fractal distribution of nucleation sites on boiling surfaces. The proposed fractal model for the subcooled flow boiling heat transfer is found to be a function of wall superheat, liquid subcooling, bulk velocity of fluid (or Reynolds number), fractal dimension, the minimum and maximum active cavity size, the contact angle and physical properties of fluid. No additional/new empirical constant is introduced, and the proposed model contains less empirical constants than the conventional models. The proposed model takes into account all the possible mechanisms for subcooled flow boiling heat transfer. The model predictions are compared with the existing experimental data, and fair agreement between the model predictions and experimental data is found for different bulk flow rates.     

LES of Low to High Turbulent Combustion in an Elevated Pressure Environment

A subgrid scale flame surface density combustion model for the Large Eddy Simulation (LES) of premixed combustion is derived and validated. The model is based on fractal characteristics of the flame surface, assuming a self similar wrinkling of the flame between smallest and largest wrinkling length scales. Experimental and direct numerical simulation databases as well as theoretical models are used to derive a model for the fractal parameters, namely the cut-off lengths and the fractal dimension suitable in the LES context. The combustion model is designed with the intent to simulate low as well as high Reynolds number premixed turbulent flame propagation and with a focus on correct scaling with pressure. The combustion model is validated by simulations of turbulent Bunsen flames with methane and propane fuel at pressure levels between 0.1 MPa and 2 MPa and at turbulence levels of $0 < u^{\prime }/s_{L}^{0} < 11$ , conditions typical for spark ignition engines. The predicted turbulent flame speed is in a very good agreement with the experimental data and a smooth transition from resolved flame wrinkling to fully modelled, nearly subgrid-only wrinkling is realized. Evaluating the influence of mesh resolution shows a predicted mean flame surface and turbulent flame speed independent of mesh resolution for cases with 9–86 % resolved flame surface. Additional simulations of a highly turbulent jet flame at 0.1 MPa and 0.5 MPa and the comparison with experimental data in terms of flame shape, velocity field and turbulent fluctuations validates the model also at conditions typical for gas turbines.     

A unified theory of turbulent flow

Summary A critical analysis is made of the assumptions underlying Reynolds' equations for turbulent flow. It is shown that there are regions in a flow field where these assumptions break down, and it is therefore necessary to separate the flow into a turbulent core and a laminar sublayer. The importance of the boundary conditions to be imposed on the mean velocities and Reynolds stresses at the junction is emphasized as this is the way in which the effect of surface roughness enters the theory. A set of equations for calculating turbulent flows is proposed. The distinctive feature is that the turbulent stresses are represented as the difference between viscous terms with a large eddy viscosity and terms satisfying auxiliary differential equations proposed by Broszko. These terms may be associated with the free and wall turbulence respectively. The theory enables the idea of a large eddy viscosity to be applied even where the velocity gradient is large. The results obtained for specific configurations, which will be reported in detail in future papers, are previewed.     

A Novel Framework for Integration of Random-Walk Particle-Tracking Simulation in Subsurface Multi-phase Flow Modeling

Coarse-scale models are generally preferred in the numerical simulation of multi-phase flow due to computational constraints. However, capturing the effects of fine-scale heterogeneity on flow and isolating the impacts of numerical (artificial) dispersion, which increases with scale, are not trivial. In this paper, a particle-tracking method is devised and integrated in a scale-up workflow to estimate the conditional probability distributions of multi-phase flow functions, which can be considered as inputs in coarse-scale simulations with existing commercial packages. First, a novel particle-tracking method is developed to solve the saturation transport equation. The transport calculation is coupled with a velocity update, following the implicit pressure, explicit saturation framework, to solve the governing equations of two-phase immiscible flow. Each phase particle is advanced in a deterministic convection step according to the phase velocity, as well as in a stochastic dispersion step based on the random Brownian motion. A kernel-based formulation is proposed for computation of fluid saturation in accordance with the phase particle distribution. A novel aspect is that this method employs the kernel approach to construct saturation from phase particle distribution, which is an important improvement to the conventional box method that necessitates a large number of particles per grid cell for consistent saturation interpolation. The model is validated against various analytical solutions. Finally, the validated model is integrated in a statistical scale-up procedure to calibrate effective, or “pseudo,” multi-phase flow functions (e.g., relative permeability functions). The proposed scale-up framework does not impose any length scale requirement regarding the distribution of sub-grid heterogeneities.     

再悬浮底泥中非吸附性污染物释放的数值模拟

环境水动力学中, 湖库底泥中污染物释放是人们研究的主要问题之一. 在水动力学条件作用下, 污染底泥再悬浮使大量污染物被重新释放出来, 造成水体的二次污染. 本文基于水槽实验研究提供的大量实测数据, 建立上覆水体-底泥-污染物的耦合力学模型. 在上覆水体不同流速条件下, 数值模拟底泥起动再悬浮过程以及污染物释放过程. 分析流场特性和污染物浓度之间的关系, 得到速度、颗粒体积分数、污染物浓度、湍动能以及时间等参数之间的定量关系. 研究表明, 底泥再悬浮污染物释放过程, 是由上覆水体-底泥-污染物构成的耦合过程. 迅速进入上覆水体的底泥颗粒, 影响了上覆水体流动特性, 进而影响到污染物的释放. 对于非吸附性污染物, 底泥起动后复杂的流场特性是底泥再悬浮污染物释放的主要影响因素. 当流场特性(如雷诺数)改变时, 对流和湍流扩散作用在污染物输运过程贡献不同. 建立上覆水体-底泥-污染物的耦合模型, 研究水动力学条件与底泥污染物释放规律的定量化关系, 可为构建湖库区域水污染模型提供支撑.      

Numerical analysis of turbulent flow in a rotating U-bend with roughened walls

A numerical analysis has been performed for a developing turbulent flow in a rotating U-bend of strong curvature with rib-roughened walls using an anisotropic turbulent model. In this calculation, an algebraic Reynolds stress model is used to precisely predict Reynolds stresses, and a boundary-fitted coordinate system is introduced as a method of coordinate transformation to set the exact boundary conditions along the complicated shape of U-bend with rib-roughened walls. Calculated results for mean velocity and Reynolds stresses are compared to the experimental data in order to validate the proposed numerical method and the algebraic Reynolds stress model. Although agreement is certainly not perfect in all details, the present method can predict characteristic velocity profiles and reproduce the separated flow generated near the outer wall, which is located just downstream of the curved duct. The Reynolds stresses predicted by the proposed turbulent model agree well with the experimental data, except in regions of flow separation.     

Turbulent particle dispersion in arbitrary wall-bounded geometries: A coupled CFD-Langevin-equation based approach

A Lagrangian continuous random walk (CRW) model is developed to predict turbulent particle dispersion in arbitrary wall-bounded flows with prevailing anisotropic, inhomogeneous turbulence. The particle tracking model uses 3D mean flow data obtained from the Fluent CFD code, as well as Eulerian statistics of instantaneous quantities computed from DNS databases. The turbulent fluid velocities at the current time step are related to those of the previous time step through a Markov chain based on the normalized Langevin equation which takes into account turbulence inhomogeneities. The model includes a drift velocity correction that considerably reduces unphysical features common in random walk models. It is shown that the model satisfies the well-mixed criterion such that tracer particles retain approximately uniform concentrations when introduced uniformly in the domain, while their deposition velocity is vanishingly small, as it should be. To handle arbitrary geometries, it is assumed that the velocity rms values in the boundary layer can locally be approximated by the DNS data of fully developed channel flows. Benchmarks of the model are performed against particle deposition data in turbulent pipe flows, 90° bends, as well as more complex 3D flows inside a mouth-throat geometry. Good agreement with the data is obtained across the range of particle inertia.     

雷诺应力各向异性涡黏模型的层析TRPIV测量

利用层析TRPIV测量水洞中平板湍流边界层3D-3C速度场的高分辨率时间序列数据库. 提出了空间局部平均多尺度速度结构函数的新概念, 描述湍流多尺度涡结构的空间拉伸、压缩、剪切变形和旋转. 用空间局部平均多尺度速度结构函数对湍流脉动速度进行了空间多尺度分解. 用空间流向局部平均多尺度速度结构函数, 根据湍流多尺度涡结构在流向的拉伸和压缩物理特征, 提出了新的湍流相干结构条件采样方法, 检测并提取了层析TRPIV数据中相干结构“喷射”和“扫掠”事件中的脉动速度、平均速度变形率、雷诺应力等物理量的空间拓扑形态. 通过研究平均速度变形率各分量与雷诺应力各分量之间的空间相位差异,肯定了壁湍流相干结构雷诺应力各向异性复涡黏模型的合理性.      

Scalar transport from a point source in flows over wavy walls

Simultaneous measurements of the velocity and concentration field in fully developed turbulent flows over a wavy wall are described. The concentration field originates from a low-momentum plume of a passive tracer. PLIF and digital particle image velocimetry are used to make spatially resolved measurements of the structure of the scalar distribution and the velocity. The measurements are performed at three different Reynolds numbers of Re _b = 5,600, Re _b = 11,200 and Re _b = 22,400, respectively, based on the bulk velocity u _b and the total channel height 2h. The velocity field and the scalar field are investigated in a water channel with an aspect ratio of 12:1, where the bottom wall of the test section consists of a train of sinusoidal waves. The wavy wall is characterized by the amplitude to wavelength ratio α = 0.05 and the ratio β between the wave amplitude and the half channel height where β = 0.1. The scalar is released from a point source at the wave crest. For the concentration measurements, Rhodamine B is used as tracer dye. At low to moderate Reynolds number, the flow field is characterized through a recirculation zone which develops after the wave crest. The recirculation zone induces high intensities of the fluctuations of the streamwise velocity and wall-normal velocity. Furthermore, large-scale structures are apparent in the flow field. In previous investigations it has been shown that these large-scale structures meander laterally in flows over wavy bottom walls. The investigations show a strong effect of the wavy bottom wall on the scalar mixing. In the vicinity of the source, the scalar is transported by packets of fluid with a high scalar concentration. As they move downstream, these packets disintegrate into filament-like structures which are subject to strong gradients between the filaments and the surrounding fluid. The lateral scale of the turbulent plume is smaller than the lateral scale of the large-scale structures in the flow field and the plume dispersion is dominated by the structures in the flow field. Due to the lateral meandering of the large-scale structures of the flow field, also the scalar plume meanders laterally. Compared to turbulent plumes in plane channel flows, the wavy bottom wall enhances the mixing effect of the turbulent flow and the spreading rate of the scalar plume is increased.