Validation of a novel very large eddy simulation method for simulation of turbulent separated flow

The paper describes the validation of a newly developed very LES (VLES) method for the simulation of turbulent separated flow. The new VLES method is a unified simulation approach that can change seamlessly from Reynolds‐averaged Navier–Stokes to DNS depending on the numerical resolution. Four complex test cases are selected to validate the performance of the new method, that is, the flow past a square cylinder at Re = 3000 confined in a channel (with a blockage ratio of 20%), the turbulent flow over a circular cylinder at Re = 3900 as well as Re = 140,000, and a turbulent backward‐facing step flow with a thick incoming boundary layer at Re = 40,000. The simulation results are compared with available experimental, LES, and detached eddy simulation‐type results. The new VLES model performs well overall, and the predictions are satisfactory compared with previous experimental and numerical results. It is observed that the new VLES method is quite efficient for the turbulent flow simulations; that is, good predictions can be obtained using a quite coarse mesh compared with the previous LES method. Discussions of the implementation of the present VLES modeling are also conducted on the basis of the simulations of turbulent channel flow up to high Reynolds number of Re_τ = 4000. The efficiency of the present VLES modeling is also observed in the channel flow simulation. From a practical point of view, this new method has considerable potential for more complex turbulent flow simulations at relative high Reynolds numbers. Copyright © 2013 John Wiley & Sons, Ltd.     

Closure modeling in bridging regions of variable-resolution (VR) turbulence computations

Optimal utilization of computational resources mandates spatio-temporal variation in resolution for computing complex engineering flows. Closure modeling in regions bridging between different resolutions is rendered difficult due to changing interactions between resolved and unresolved fields. We develop a closure model for the bridging region based on energy conservation principles. Then we proceed to provide a proof of concept in decaying isotropic turbulence with temporally varying resolution. The simplicity of the flow permits a thorough examination of various aspects of the proposed closure not feasible in more complex flows. The results demonstrate the potential promise of the approach, but more validation studies need to be performed.While the present development is in the context of partially averaged Navier–Stokes (PANS) method, the closure principle should apply for other variable-resolution (VR) approaches.     

Energy flow along the medium-induced parton cascade

We discuss the dynamics of parton cascades that develop in dense QCD matter, and contrast their properties with those of similar cascades of gluon radiation in vacuum. We argue that such cascades belong to two distinct classes that are characterized respectively by an increasing or a constant (or decreasing) branching rate along the cascade. In the former class, of which the BDMPS, medium-induced, cascade constitutes a typical example, it takes a finite time to transport a finite amount of energy to very soft quanta, while this time is essentially infinite in the latter case, to which the DGLAP cascade belongs. The medium induced cascade is accompanied by a constant flow of energy towards arbitrary soft modes, leading eventually to the accumulation of the initial energy of the leading particle at zero energy. It also exhibits scaling properties akin to wave turbulence. These properties do not show up in the cascade that develops in vacuum. There, the energy accumulates in the spectrum at smaller and smaller energy as the cascade develops, but the energy never flows all the way down to zero energy. Our analysis suggests that the way the energy is shared among the offsprings of a splitting gluon has little impact on the qualitative properties of the cascades, provided the kernel that governs the splittings is not too singular.     

湍流大气传输中跟踪抖动对远场影响的仿真研究

分析了跟踪抖动对湍流大气传输远场光斑的影响。基于麦克斯韦电磁场理论,采用大气相干长度对大气湍流进行描述,推导了发射光束因跟踪抖动导致光轴偏离的远场表达式。在此基础上,利用相位屏法模拟抖动引起的倾斜相位和大气折射率起伏引起的相位调制,并采用低频补偿的功率谱反演法对传输过程进行了数值仿真。分析了不同跟踪抖动、湍流强度条件下远场光斑质心脱靶量的变化,以及不同尺寸模拟目标的回波概率。分析结果表明,在传输距离为10 km时,强湍流造成的远场光斑脱靶量可达几十μrad;当跟踪抖动较大时,湍流强弱对脱靶量影响差别很小。最后,对一定尺寸的模拟目标,从探测回波概率的角度给出了发射系统跟踪抖动量的控制范围。     

A mathematical and physical Study of multiscale deconvolution models of turbulence

This work presents a rigorous analysis of mathematical and physical properties for solutions of multiscale deconvolution turbulence models. We show that solutions of these models exactly conserve model quantities for the integral invariants of fundamental physical importance: kinetic energy, helicity, and (in two dimensions) enstrophy. The kinetic energy conservation is the key that allows us to next apply the phenomenology of homogeneous, isotropic turbulence to establish the existence of a model energy cascade and, in particular, that the cascade exhibits enhanced energy dissipation in a secondary accelerated cascade, which ends at the model's microscale (which we establish is larger than the Kolmogorov microscale). We also prove that the model dissipates energy at the same rate as true turbulent flow, ～ O(U³ ∕ L), independent of Reynolds number. Lastly, we prove the existence of global attractors for the model solutions; the proof of which also shows that solutions are actually one degree of regularity higher than previously known. Copyright © 2012 John Wiley & Sons, Ltd.     

Effects of multi-Schell beams on orbital angular momentum of entangled signal photons in low-order turbulence aberration channels

We model the detection and crosstalk probability of orbital angular momentum (OAM) states of the entangled signal photon in the Kolmogorov channels of the low-order turbulence aberrations and by the Rytov approximation. The results show that lower OAM mode number of signal photons and larger sub-beam number of multi-Gaussian Schell-model pump beam, the less susceptible of the detection probability of the signal photon to spatial coherence of source and turbulence aberrations is achieved. The maximum crosstalk probability is decrease as the decreasing of the sub-beam number of multi-Gaussian Schell-model. Enlarging OAM difference value or decreasing sub-beam number of multi-Gaussian Schell-model pump beam results in a lower crosstalk probability of the OAM of entangled signal photons.     

Direct and adjoint global stability analysis of turbulent transonic flows over a NACA0012 profile

In this work, various turbulent solutions of the two‐dimensional (2D) and three‐dimensional compressible Reynolds averaged Navier–Stokes equations are analyzed using global stability theory. This analysis is motivated by the onset of flow unsteadiness (Hopf bifurcation) for transonic buffet conditions where moderately high Reynolds numbers and compressible effects must be considered. The buffet phenomenon involves a complex interaction between the separated flow and a shock wave. The efficient numerical methodology presented in this paper predicts the critical parameters, namely, the angle of attack and Mach and Reynolds numbers beyond which the onset of flow unsteadiness appears. The geometry, a NACA0012 profile, and flow parameters selected reproduce situations of practical interest for aeronautical applications. The numerical computation is performed in three steps. First, a steady baseflow solution is obtained; second, the Jacobian matrix for the RANS equations based on a finite volume discretization is computed; and finally, the generalized eigenvalue problem is derived when the baseflow is linearly perturbed. The methodology is validated predicting the 2D Hopf bifurcation for a circular cylinder under laminar flow condition. This benchmark shows good agreement with the previous published computations and experimental data. In the transonic buffet case, the baseflow is computed using the Spalart–Allmaras turbulence model and represents a mean flow where the high frequency content and length scales of the order of the shear‐layer thickness have been averaged. The lower frequency content is assumed to be decoupled from the high frequencies, thus allowing a stability analysis to be performed on the low frequency range. In addition, results of the corresponding adjoint problem and the sensitivity map are provided for the first time for the buffet problem. Finally, an extruded three‐dimensional geometry of the NACA0012 airfoil, where all velocity components are considered, was also analyzed as a Triglobal stability case, and the outcoming results were compared to the previous 2D limited model, confirming that the buffet onset is well detected. Copyright © 2014 John Wiley & Sons, Ltd.     

Eulerian spatial and temporal autocorrelations: assessment of Taylor's hypothesis and a model

We present measurements of Eulerian longitudinal velocity autocorrelations in homogeneous, isotropic, high-intensity (～9%) free-stream turbulence behind an active grid. Spatial correlations are measured using particle image velocimetry as well as with two-point hot-wire anemometry (HWA), while temporal correlations are measured using HWA. The temporal correlations are transformed into spatial correlations by using Taylor's ‘frozen’ hypothesis with both the mean as well as instantaneous velocities. A model relating Eulerian spatial and temporal autocorrelations is also used for this purpose. The differences from the measured spatial correlation resulting from the use of Taylor's hypothesis on the temporal correlation is quantified; even at this moderately high level of turbulent intensity, the result from the use of the instantaneous velocity as convection velocity is practically indistinguishable from that obtained using the mean velocity. Use of the model produces a good agreement between the estimates of the spatial correlation function. A relation between Eulerian spatial and temporal integral scales is also derived.     

Analysis of marine atmospheric turbulence effects on infrared imaging system by angle of arrival fluctuations

For the long-range infrared imaging system, the marine atmospheric turbulence degrades seriously the probability of object recognition and tracking. In this study, the angle of arrival fluctuations of an optical wave, which describes the distortion effects of marine atmospheric turbulence on an infrared optical imaging system, is investigated in detail both analytically and numerically. Analytic expressions of the angle of arrival fluctuations are derived for optical plane and spherical waves propagating through weak marine atmospheric turbulence with horizontal path, and they consider simultaneously finite turbulence inner scale, turbulence outer scale, wavelength, and aperture diameters. Numerical calculations are conducted to analyze the influence of marine weak turbulence on the infrared imaging. The results are useful for understanding the potential impact of deviations from the terrestrial turbulence.