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.     

湍流内尺度对光波起伏频谱的影响

在不同湍流内尺度情况下,对大气湍流引起的光波光强闪烁和到达角起伏的时间频谱特征进行了实验研究。由于实际大气湍流的复杂性和不可控性,利用大气湍流模拟箱来生成具有不同内尺度的大气湍流实验环境。利用位置敏感探测器对光波的光强闪烁和到达角起伏进行了同时测量,并反演得到了光传输路径上的湍流内尺度。实验结果显示:湍流内尺度为2.7~5.0 mm,对应同一湍流内尺度,闪烁频谱和到达角起伏频谱在高频段以相同的幂指数关系下降,幂指数的绝对值与湍流内尺度的大小呈线性关系,随着内尺度的增大,频谱指数变化区间逐渐向低频方向移动。     

一种K分布强湍流下的测量设备无关量子密钥分发方案

研究了K分布强湍流下自由空间测量设备无关量子密钥分发协议模型,采用阈值后选择方法来减少大气湍流对密钥生成率的影响,对比分析了使用阈值后选择方法前后协议的密钥率和湍流强度之间的关系.仿真结果表明,使用阈值后选择方法可以有效地提高协议的密钥生成率,尤其是在高损耗和强湍流区域,而且其最佳阈值与湍流强度、信道平均损耗有关,对实际搭建性能较好的自由空间测量设备无关量子密钥分发协议系统具有一定的参考价值.     

Effect of gravity on the development of homogeneous shear turbulence laden with finite-size particles

Fully resolved simulations of homogeneous shear turbulence (HST) laden with sedimenting spherical particles of finite size have been performed to clarify the effects of gravity on the development of particle-laden turbulent shear flows. We consider turbulence in a horizontal flow subjected to vertical or horizontal shear. Numerical results show that the development of HST laden with finite-size particles are significantly altered by gravity. The effects of gravity lead to a slower increase in the Taylor-microscale Reynolds number, whose value is found to be well correlated with the average particle Reynolds number. The gravity also causes a slower increase in the turbulence kinetic energy (TKE) through the enhancement of energy dissipation. The change in the Reynolds shear stress (RSS) due to particles also significantly contributes to the relative change in TKE. In vertically sheared cases, RSS has high values between counter-rotating trailing vortices behind the particles, which causes a transient relative increase in TKE. In horizontally sheared cases, on the other hand, RSS is reduced in the wakes of particles, which contributes to a significant relative reduction in TKE.     

Lens ghost imaging in a non-Kolmogorov slant turbulence atmosphere

Based on the method of computation reference channel and the extended Huygens–Fresnel integral, the model of lens ghost imaging with fully spatially incoherent linear polarization light through a slant turbulent channel has derived. The model shows that the resolution ratio of ghost imaging decreases as the power-law exponent of non-Kolmogorov turbulence increasing or the object location departing from the source. The zenith angle of channel has little influence to the quality of the ghost imaging. The minimum distinguishable centre-separation of double slit decreases as the power-law exponent of non-Kolmogorov turbulence increasing.     

Particle subgrid scale modelling in large-eddy simulations of particle-laden turbulence

This study is concerned with particle subgrid scale (SGS) modelling in large-eddy simulations (LESs) of particle-laden turbulence. Although many particle-laden LES studies have neglected the effect of the SGS on the particles, several particle SGS models have been proposed in the literature. In this research, the approximate deconvolution method (ADM) and the stochastic models of Fukagata et al. (Dynamics of Brownian particles in a turbulent channel flow, Heat Mass Transf. 40 (2004), 715–726) Shotorban and Mashayek (A stochastic model for particle motion in large-eddy simulation, J. Turbul. 7 (2006), 1–13) and Berrouk et al. (Stochastic modelling of inertial particle dispersion by subgrid motion for LES of high Reynolds number pipe flow, J. Turbul. 8 (2007), pp. 1–20) are analysed. The particle SGS models are assessed using both a priori and a posteriori simulations of inertial particles in a periodic box of decaying, homogeneous and isotropic turbulence with an initial Reynolds number of Re_λ = 74. The model results are compared with particle statistics from a direct numerical simulation (DNS). Particles with a large range of Stokes numbers are tested using various filter sizes and stochastic model constant values. Simulations with and without gravity are performed to evaluate the ability of the models to account for the crossing trajectory and continuity effects. The results show that ADM improves results but is only capable of recovering a portion of the SGS turbulent kinetic energy. Conversely, the stochastic models are able to recover sufficient SGS energy, but show a large range of results dependent on the Stokes number and filter size. The stochastic models generally perform best at small Stokes numbers, but are unable to predict preferential concentration.     

Direct numerical simulation of elastic turbulence and its mixing-enhancement effect in a straight channel flow

In this paper,we present a direct numerical simulation(DNS) of elastic turbulence of viscoelastic fluid at vanishingly low Reynolds number(Re = 1) in a three-dimensional straight channel flow for the first time,using the Giesekus constitutive model for the fluid.In order to generate and maintain the turbulent fluid motion in the straight channel,a sinusoidal force term is added to the momentum equation,and then the elastic turbulence is numerically realized with an initialized chaotic velocity field and a stretched conformation field.Statistical and structural characteristics of the elastic turbulence therein are analyzed based on the detailed information obtained from the DNS.The fluid mixing enhancement effect of elastic turbulence is also demonstrated for the potential applications of this phenomenon.     

大气湍流对激光通信系统的影响

从分析激光在大气湍流场中的传输方程出发,忽略系统中的其它噪声,仅考虑由大气喘流引起的系统误码率,讨论了激光信号在传输过程中的振幅起伏以及强度起伏;推导出由大气湍流引起的光强起伏和系统误码率的关系,结果表明:在弱起伏条件下,对于系统误码率为10－9以下的要求,光强起伏应小于0.67;随着湍流强度C2n的增大,误码率增加很快;采用长波长的激光进行传输可以有效地降低系统误码率.     

Capacity of wander and spread beams in log-normal distribution non-Kolmogorov turbulence optical links

We model the average channel capacity of optical wireless communication systems in weak turbulence horizontal channels, using the log-normal distribution models. The effects of beam wander and spread, pointing errors, turbulence inner scale, turbulence outer scale and the spectral index of non-Kolmogorov turbulence on system's performance are included. The model can evaluate the influence of the atmospheric turbulence conditions in the performance of a ground-to-train optical wireless communication system.     

On the role of pressure in elasto-inertial turbulence

The dynamics of elasto-inertial turbulence is investigated numerically from the perspective of the coupling between polymer dynamics and flow structures. In particular, direct numerical simulations of channel flow with Reynolds numbers ranging from 1000 to 6000 are used to study the formation and dynamics of elastic instabilities and their effects on the flow. Based on the splitting of the pressure into inertial and polymeric contributions, it is shown that the polymeric pressure is a non-negligible component of the total pressure fluctuations, although the rapid inertial part dominates. Unlike Newtonian flows, the slow inertial part is almost negligible in elasto-inertial turbulence. Statistics on the different terms of the Reynolds stress transport equation also illustrate the energy transfers between polymers and turbulence and the redistributive role of pressure. Finally, the trains of cylindrical structures around sheets of high polymer extension that are characteristics of elasto-inertial turbulence are shown to be correlated with the polymeric pressure fluctuations.     

Influence of oceanic turbulence on propagation characteristics of Gaussian array beams

Based on the power spectrum of oceanic turbulence proposed by Nikishov, the influence of oceanic optical turbulence on propagation of Gaussian array beams is studied in detail by using five propagation characteristic parameters (i.e., mean-squared beam width, Rayleigh range, turbulence distance, power in the bucket and Strehl ratio). It is shown that the influence of oceanic optical turbulence on propagation of laser beams becomes stronger as τ (ratio of temperature to salinity contribution to the refractive index spectrum) and χ_T (rate of dissipation of mean-squared temperature) increase and ɛ (rate of dissipation of kinetic energy per unit mass of fluid) decreases, which is in agreement with the further analysis of oceanic turbulence shown in this paper. Furthermore, it is concluded that propagation of laser beams is more affected by oceanic optical turbulence in abyssal region than that in active region or oceanic surface.     

Investigation of anomalous very fast decay regimes in homogeneous isotropic turbulence

The emergence of anomalous fast decay regimes in homogeneous isotropic turbulence (HIT) decay is investigated via both theoretical analysis and eddy-damped quasi-normal Markovian simulations. The work provides new insight about a fundamental issue playing a role in HIT decay, namely the influence of non-standard shapes of the energy spectrum, in particular in the large energetic scale region. A detailed analysis of the kinetic energy spectrum E(k) and the non-linear energy transfer T(k) shows that anomalous decay regimes are associated with the relaxation of initial energy spectra which exhibit a bump at energetic scales. This feature induces an increase in the energy cascade rate, toward solutions with a smooth shape at the spectrum peak. Present results match observations reported in wind-tunnel experiments dealing with turbulence decay in the wake of grids and bluff bodies, including scaling laws for the dissipation parameter C_?. They also indicate that the ratio between the initial eddy turnover time and the advection time determines of how fast anomalous regimes relax toward classical turbulence free-decay. This parameter should be used for consistent data comparison and it opens perspectives for the control of multiscale effects in industrial applications.     

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.     

基于拉盖尔-高斯光束的通信系统在非Kolmogorov湍流中传输的系统容量

涡旋波束在大气湍流中的传输有非常重要的理论研究和实际应用意义. 本文基于利托夫近似和广义惠更斯-菲涅耳原理, 推导出拉盖尔-高斯(LG)光束在非Kolmogorov湍流中斜程传输时的螺旋谱, 并进一步推导出系统的容量. 对基于LG光束的通信系统容量进行了数值计算, 并对指数参数、光束波长、天顶角、湍流内尺度、外尺度、结构常数对系统容量的影响进行了分析比较. 本文的结论能够为LG光束在非Kolmogorov湍流中的通信提供一定的参考价值. 关键词： 拉盖尔-高斯光束 非Kolmogorov湍流 平均容量     

Degree of polarization for quantum light field propagating through non-Kolmogorov turbulence

Nonclassical polarization properties of a quantum field propagating through non-Kolmogorov turbulence are studied in a turbulent atmosphere paraxial channel. The analytic equation for the quantum degree of polarization of linearly polarized light is obtained. It is shown by numerical simulation that the polarization fluctuations of the quantum field are a function of the turbulent strength, the photon number, the propagation distance, the fractal constant α and the coherence length ρ₀.     

汉克-贝塞尔光束在海洋湍流信道中的螺旋相位谱分析

携带轨道角动量的汉克-贝塞尔(Hankel-Bessel,HB)光束具有无衍射和自聚焦特性,用来作为信息传输的载体有望增大信息传输容量.基于Rytov近似理论,推导得到了HB涡旋光束经过海洋水平弱湍流信道后的螺旋相位谱的解析表达式,并利用数值仿真方法研究了海洋湍流参数对轨道角动量模式探测概率的影响.结果表明,海洋湍流导致发射轨道角动量模式的探测概率下降,出现模式串扰和螺旋相位谱扩展.海洋湍流对HB涡旋光束的负面影响随着轨道角动量模式数、传输距离、温度方差耗散率的增加而增强,随湍流动能耗散率的增加而减弱.HB涡旋光束受以盐度波动驱动的海洋湍流的负面影响更大.另外,在弱湍流及几十米传输距离条件下,HB涡旋光束的传输性能要差于最佳束腰大小设置的拉盖尔-高斯涡旋光束.这些结果有望为海洋环境水下光通信链路的实现提供一定的参考价值.     

间歇湍流的分数阶动力学

间歇湍流意味着湍流涡旋并不充满空间,其维数介于2和3之间.湍流扩散为超扩散,且概率密度分布具有长尾特征.本文将流体力学的Navier-Stokes(NS)方程中的黏性项用分数阶的拉普拉斯算子表达.分析表明,分数阶拉普拉斯的阶数α和间歇湍流的维数D相联系.对于均匀各向同性的Kolmogorov湍流α=2,即用整数阶NS方程描述.而对于间歇性湍流,一定用分数阶的NS方程来描述.对于Kolmogorov湍流,扩散方差正比于t3,即Richardson扩散.而对于间歇性湍流,扩散方差要比Richardson扩散更强.     

Wave turbulent statistics in non-weak wave turbulence

In wave turbulence, which is made by nonlinear interactions among waves, it has been believed that statistical properties are well described by the weak turbulence theory, where separation of linear and nonlinear time scales derived from weak nonlinearity is assumed. However, the separation of the time scales is often violated. To get rid of this inconsistency, closed equations are derived in wave turbulence without assuming the weak nonlinearity according to Direct-Interaction Approximation (DIA), which has been successful in Navier-Stokes turbulence. The DIA equations is a natural extension of the conventional kinetic equation to not-necessarily-weak wave turbulence.     

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.