Modulation of homogeneous shear turbulence laden with finite-size particles

The dynamics of homogeneous shear turbulence laden with spherical finite-size particles is investigated using fully resolved numerical simulations to understand how the presence of particles modulates turbulent shear flows. We focus on a dilute flow laden with non-sedimenting particles whose diameter is slightly smaller than or comparable with those of vortex cores in turbulence. An immersed boundary method is adopted to represent a spherical finite-size particle. Numerical results show that the presence of particles augments the viscous dissipation of turbulence kinetic energy, which leads to a slower increase in the turbulence energy. Although the augmentation of energy dissipation occurs predominantly inside viscous layers surrounding particles in an initial period, the contribution from their outside becomes more significant due to the modification of turbulence structures as turbulence develops. It is found that the particles exhibit weak tendency to accumulate in vortex layers. The particles approaching and colliding with vortex layers induce large velocity fluctuations, which leads to the generation and shedding of thin vortex tubes. Newly generated vortex tubes interact with developed vortex tubes and layers, and modify the entire structure of the vorticity field.     

Fractal dimension of superfluid turbulence

Superfluid turbulence consists of a disordered tangle of quantized vortex filaments which interact with each other and with the normal fluid. We develop a kinematic model of normal-fluid turbulence to study superfluid vortex tangles at finite temperatures and show by numerical simulation that the system of filaments has a fractal dimension larger than one. We find that the fractal dimension is directly related to the vortex-line density and is independent of temperature over a wide range.     

Theoretical study of the decaying convective turbulence in a shear-buoyancy PBL

The derivation of a theoretical model for the decaying convective turbulence in a shear-buoyancy planetary boundary layer is considered. The model is based on the dynamical equation for the energy density spectrum in which the buoyancy, mechanical and inertial transfer terms are retained. The parameterization for the buoyancy and mechanical terms is provided by the flux Richardson number. Regarding the inertial term an approach employing Heisenberg’s spectral transfer theory is used to describe the turbulence friction, caused by small eddies, responsible for the energy dissipation of the large eddies. Therefore, a novelty in this study is to utilize the Adomian decomposition method to solve directly without linearization the energy density spectrum equation, with this the nonlinear nature of the problem is preserved. Therefore, the errors found are only due to the parameterization used. Comparison of the theoretical model is performed against large-eddy simulation data for a decaying convective turbulence in a shear-buoyancy planetary boundary layer. The results show that the existence of a mechanical turbulent driving mechanism reduces in an accentuated way the energy density spectrum and turbulent kinetic energy decay generated by the decaying convective production in a shear-buoyancy planetary boundary layer.     

Equilibrium statistics of a vortex filament with applications

The thermodynamic functions and scaling exponents (including the Kolmogorov and Flory exponents) of a vortex filament in thermal equilibrium are calculated, giving a quantitative content to earlier qualitative analyses. The numerical results uncover a percolation property of vortex filaments near the maximum entropy state. The implications of the results for the onset of turbulence, for the structure of its inertial range, and for superfluid vortices are discussed. In particular, it is shown that vortex stretching pushes a vortex system to a polymeric state and a Kolmogorov spectrum.This work was supported in part by the Applied Mathematical Sciences subprogram of the Office of Energy Research, U.S. Department of Energy, under contract DE-AC03-76SF-00098, and in part by the National Science Foundation under grant number DMS89-19074     

各向同性湍流中颗粒引起湍流变动的直接模拟

本文通过直接数值模拟对均匀各向同性湍流中颗粒对湍流的变动作用进行了研究.颗粒相的体积分数很小而质量载荷足够大,以至于颗粒之间的相互作用可以忽略不计,而重点考虑颗粒与湍流间能量的交换。颗粒对湍流的反向作用使得湍动能的耗散率增强,以至于湍动能的衰减速率增大.湍动能的衰减速率随颗粒惯性的增大而增大。三维湍动能谱显示,颗粒对湍动能的影响在不同的尺度上是不均匀的。在低波数段,流体带动颗粒,而高波数段则相反.     

A Lagrangian stochastic model for tetrad dispersion

We present results for tetrad (four-particle) dispersion in homogeneous isotropic turbulence by means of a simple Lagrangian stochastic model with a focus on the inertial subrange. We show that for appropriate values of C ₀, the constant of proportionality in the second-order Lagrangian velocity structure function, the shape statistics agree well with equivalent results from a direct numerical simulation (DNS) of turbulence. Moreover, we show that the shape statistics are independent of C ₀ for a wide range of C ₀-values. We also show that the parameters which characterise the shape of the tetrads can be approximately related to appropriate ratios of the growth rates of the mean square separation, the mean square area and the mean square volume of the tetrads. By means of exit times, we are able to estimate the equivalent values for the DNS data. We also consider the statistics of four-point velocity differences (via a diffusion tensor) which agree well with the DNS. We show that the nature of the velocity field experienced by the tetrad varies significantly with C ₀.     

Quantized Superfluid Vortex Filaments Induced by the Axial Flow Effect

We report the quantized superfluid vortex filaments induced by the axial flow effect, which exhibit intriguing loop structures on helical vortexes. Such new vortex filaments correspond to a series of soliton excitations including the multipeak soliton, W-shaped soliton, and anti-dark soliton, which have no analogue when the axial flow effect is absent. In particular, we show that the vortex filaments induced by the multipeak soliton and W-shaped soliton arise from the dual action of bending and twisting of the vortex, while the vortex filament induced by the anti-dark soliton is caused only by the bending action, which is consistent with the case of the standard bright soliton. These results will deepen our understanding of breather-induced vortex filaments and will be helpful for controllable ring-like excitations on vortices.     

On the Global Evolution of Vortex Filaments,Blobs, and Small Loops in 3D Ideal Flows

We consider a wide class of approximate models of evolution of singular distributions of vorticity in three dimensional incompressible fluids and we show that they have global smooth solutions. The proof exploits the existence of suitable Hamiltonian functions. The approximate models we analyze (essentially discrete and continuous vortex filaments and vortex loops) are related to some problem of classical physics concerning turbulence and also to the numerical approximation of flows with very high Reynolds number. Finally, we apply our strategy to discrete models for filaments used in numerical methods.     

惯性颗粒所见被动标量统计特性的DNS与模型研究

本文采用DNS方法,对惯性颗粒所见各向同性湍流中具有平均标量梯度的被动标量场统计特性进行了研究。结果表明：惯性对颗粒温度脉动强度,两相温度关联,自相关特性以及颗粒热流与两相交叉热流的统计特性具有明显的影响。在PDF方法的框架下,系统地推导了非等温气固两相流的PDF方程,且基于朗之万随机体系对方程进行了封闭,并利用前面的...     

Statistics and structures of pressure and density in compressible isotropic turbulence

We study statistics and structures of pressure and density in the presence of large-scale shock waves in a forced compressible isotropic turbulence using high-resolution numerical simulation. The spectra for pressure and density exhibit a ?2 scaling over an operational definition of the inertial range. Both the numerical simulation and a heuristic PDF model reveal that the PDFs of pressure increment exhibit a ?2 power law region for the separation in the operational definition of inertial range, quantitatively similar to the PDF of pressure gradient, which also displays a ?2 power law region. Moreover, the statistical relation between density increment and pressure increment has been investigated through a shock-relation model. There is a positive correlation between the vorticity magnitude and pressure, which is different from the case of incompressible turbulence. We argue that this difference is due to large-scale shock waves, another type of intermittent structures in addition to vortex structures in incompressible turbulence.     

A subgrid model for clustering of high-inertia particles in large-eddy simulations of turbulence

Clustering (or preferential concentration) of inertial particles suspended in a homogeneous, isotropic turbulent flow is strongly influenced by the smallest scales of the turbulence. In particle-laden large-eddy simulations (LES) of turbulence, these small scales are not captured by the grid and hence their effect on particle motion needs to be modelled. In this paper, we use a subgrid model based on kinematic simulations of turbulence (Kinematic Simulation based SubGrid Model or KSSGM), for the first time in the context of predicting the clustering and the relative velocity statistics of inertial particles. This initial study focuses on the special case of inertial particles in the absence of gravitational settling. We show that the KSSGM gives excellent predictions for clustering in a priori tests for inertial particles with St ≥ 2.0, where St is the Stokes number, defined as the ratio of the particle response time to the Kolmogorov time-scale. To the best of our knowledge, the KSSGM represents the first model that has been shown to capture the effect of the subgrid scales on inertial particle clustering for St ≥ 2.0. We also show that the mean inward radial relative velocity between inertial particles (?w_r?^(?), which enters into the formula for the collision kernel) is accurately predicted by the KSSGM for all St. We explain why the model captures clustering at higher St?but not for lower St?, and provide new insights into the key statistical parameters of turbulence that a subgrid model would have to describe, in order to accurately predict clustering of low-St?particles in an LES.     

Physics of multiple level hairpin vortex structures in turbulence

Previous experimental and numerical studies have revealed that the hairpin vortex is a basic flow element of transitional boundary layer. The hairpin vortex is believed to have legs, necks and a ring head. Based on our DNS study, the legs and the ring head are generated separately by different mechanisms. The legs function like an engine to generate low speed zones by rotation, create shear layers with surrounding high speed neighbor fluids, and further cause vortex ring formation through shear layer instability. In addition, the ring head is ?-shaped and separated from quasi-streamwise legs from the beginning. Contrary to the classical concept of "vortex breakdown", we believe transition from laminar flow to turbulence is a "buildup" process of multiple level vortical structures. The vortex rings of first level hairpins are mostly responsible for positive spikes, which cause new vorticity rollup, second level vortex leg formation and finally smaller second level vortex ring generation. The third and lower level vortices are generated following the same mechanism. In this paper, the physical process from ?-vortex to multi-level hairpin vortices is described in detail.     

Lagrangian measurements of inertial particle accelerations in grid generated wind tunnel turbulence

We describe Lagrangian measurements of water droplets in grid generated wind tunnel turbulence at a Taylor Reynolds number of R(lambda)=250 and an average Stokes number (St) of approximately 0.1. The inertial particles are tracked by a high speed camera moving along the side of the tunnel at the mean flow speed. The standardized acceleration probability density functions of the particles have spread exponential tails that are narrower than those of a fluid particles (St approximately 0) and there is a decrease in the acceleration variance with increasing Stokes number. A simple vortex model shows that the inertial particles selectively sample the fluid field and are less likely to experience regions of the fluid undergoing the largest accelerations. Recent direct numerical simulations compare favorably with these first measurements of Lagrangian statistics of inertial particles in highly turbulent flows.     

混合层流场中涡结构对流速度的特性

基于大涡模拟和光线追踪方法, 对光线穿越流场后的光程分布与混合层流场中涡结构之间的关系进行了分析, 提出了一种基于涡核位置提取的涡结构瞬时对流速度定量计算方法, 并使用直接几何测量数据进行了验证. 通过对不同尺寸的涡结构、涡-涡配对及融合过程中的涡结构和强压缩性流场中涡结构瞬时对流速度的定量数值计算, 揭示了混合层流场中涡结构对流速度的特性: 对单个涡结构而言, 其瞬时对流速度具有脉动特性, 且脉动幅度随涡结构尺寸和流场压缩性而变化; 在涡-涡配对及融合过程中, 涡对中各个涡结构的瞬时对流速度都表现出类似正弦波动的特点. 针对混合层流场中涡结构对流速度的特性, 给出了其背后的物理原因.     

Emission of phonons by vortex filaments and the decay of a turbulent state in the Bose condensate

A turbulence decay mechanism in the Bose condensate (superfluid liquid) that is associated with the emission of acoustic waves upon reconnection of vortex filaments is studied. An approach based on the master equation for the length distribution of vortex loops is developed. Evolution equations for the density of the vortex filaments in the quasistationary and non-quasistationary regimes are derived. The decay law is compared to the numerical simulation results.     

确定分布的展向Lorentz力调制下的槽道湍流涡结构

采用直接数值模拟方法,对槽道湍流中确定分布的Lorentz力的流动控制与减阻问题进行研究.讨论了Lorentz力作用于槽道湍流后,流场的特性和涡结构的特性,并对此类Lorentz力对槽道湍流的控制与减阻机理进行了讨论.研究发现:1)Lorentz力诱导的层流流场壁面附近存在梯度极大的展向速度剪切层,该剪切层容易形成流向涡结构;2)在给定合适参数的确定分布的Lorentz力作用下,湍流流场仅剩周期分布的准流向涡;3)与未控制流场相比,控制后的流场中,准流向涡的抬升高度大大降低,从而减小猝发强度,使壁面阻力下降.     

Large scale dissipation and filament instability in two-dimensional turbulence

Coherent vortices in two-dimensional turbulence induce far-field effects that stabilize vorticity filaments and inhibit the generation of new vortices. We show that the large-scale energy sink often included in numerical simulations of statistically stationary two-dimensional turbulence reduces the stabilizing role of the vortices, leading to filament instability and to continuous formation of new coherent vortices. This counterintuitive effect sheds new light on the mechanisms responsible for vortex formation in forced-dissipated two-dimensional turbulence, and it has significant impact on the temporal evolution of the vortex population in freely decaying turbulence. The time dependence of vortex statistics in the presence of a large-scale energy sink can be approximately described by a modified version of the scaling theory developed for small-scale dissipation.     

Modeling of classical turbulence by quantized vortices

Quantum turbulence in superfluids appears as a stochastic tangle of quantized vortex lines. Interest to this system extends beyond the field of superfluid helium to include a large variety of topics both fundamental and engineering problem. In the article we present a discussion of the hot topic, which is undoubtedly mainstream in this field, and which deals with the quasi-classical properties of quantum turbulence. The idea that classical turbulence can be modeled by a set of slim vortex tubes (or vortex sheets) has been discussed for quite a long time. In classical fluids, the concept of thin vortex tubes is a rather fruitful mathematical model. Quantum fluids, where the vortex filaments are real objects, give an excellent opportunity for the study of the question, whether the dynamics of a set of vortex lines is able to reproduce (at least partially) the properties of real hydrodynamic turbulence. The main goal of this article is to discuss the current state of this activity. We cover such important topics as theoretical justification of this model, and experimental and numerical evidence for quasi-classical turbulence.     

Lagrangian coherent structures analysis of gas-liquid flow in a bubble column

The objective of this paper is to apply a new identifying method to investigating the gas-liquid two-phase flow behaviors in a bubble column with air injected into water.In the numerical simulations,the standard k-?turbulence model is employed to describe the turbulence phenomenon occurring in the continuous fluid.The Finite-Time Lyapunov Exponent(FTLE)and Lagrangian Coherent Structures(LCS)are applied to analyze the vortex structures in multiphase flow.Reasonable agreements are obtained between the numerical and experimental data.The numerical results show that the evolution of gas-liquid in the column includes initial and periodical developing stages.During the initial stage,the bubble hose is forming and extending along the vertical direction with the vortex structures formed symmetrically.During the periodical developing stage,the bubble hose starts to oscillate periodically,and the vortexes move along the bubble hose to the bottom of column alternately.Compared to the Euler-system-based identification criterion of a vortex,the FTLE field presents the boundary of a vortex without any threshold defined and the LCS represents the divergence extent of infinite neighboring particles.During the initial stage,the interfaces between the forward and backward flows are highlighted by the LCS.As for the periodical developing stage,the LCS curls near the vortex centers,providing a method of analyzing a flow field from a dynamical system perspective.