Direct numerical simulation of passive scalar in decaying compressible turbulence

1IntroductionDirectnumericalsimulation(DNS)becomesanimportanttoolinrecentresearchofturbulence[1].DNSofcompressibleturbulenceismoredifficultthanthatoftheincompressibleturbulence.WhentheturbulentMachnumberisgreaterthan0.3theshockletsmayappearinthecompressibleturbulentflowfields.Thereasonandmechanismofshockletsexistencearenotclearyet.TheturbulentMachnumberinDNScannotbeveryhighwiththepresentexistingnumericalmethodsandcomputerresource.Fortheproblemofcompressibleisotropicturbulencewiththeinitia…     

Restricted scaling range models for turbulent velocity and scalar energy transfers in decaying turbulence

The effect of finite Reynolds numbers and/or internal intermittency on the total kinetic energy and scalar energy transfers is examined in detail. For this purpose, two distinct models for velocity and scalar energy transfer are proposed in the specific context of freely decaying isotropic turbulence. The first one extends the already existing dynamical models (hereafter DYM, i.e. based on transport equations originated in Navier–Stokes and advection-diffusion transport equations). The second one relies on the characteristic time of the strain at a specific scale (hereafter SBM). Both models account for the Reynolds number dependence of the scaling exponent of the second-order structure functions, over a range of scales where such exponents may be defined, i.e. a restricted scaling range (RSR). Therefore, the models developed aim at reproducing the energy transfer over the RSR. The predicted energy transfer is very sensible to variations of the scaling exponent, especially at low Reynolds numbers. The approach towards the asymptotic 4/3 law is closely reproduced by the two models. The dynamical model reproduces the experimental results accurately especially in the vicinity of the Taylor microscale, while the SBM agrees almost perfectly with measurements at nearly all scales.     

Batchelor scaling in fast-flowing soap films

The dynamics of a passive scalar such as a dye in the far dissipative range of fluid turbulence is a central problem in nonlinear physics. An important prediction for this problem was made by Batchelor over 40 years ago and is known as Batchelor's scaling law. We here present strong evidence in favor of this law for the thickness fluctuations in the flow of a soap film past a flat plate. The results also capture the dissipative range of the scalar which turns out to have universal features. The probability density function of the scalar increments and their structure functions come out in nice agreement with theoretical predictions.     

Geometry of intensive scalar dissipation events in turbulence

The maxima of the scalar dissipation rate in turbulence appear in the form of sheets and correspond to the potentially most intensive scalar mixing events. Their cross section extension determines a locally varying diffusion scale of the mixing process and extends the classical Batchelor picture of one mean diffusion scale. The distribution of the local diffusion scales is analyzed for different Reynolds and Schmidt numbers with a fast multiscale technique applied to very high-resolution simulation data. The scales always take values across the whole Batchelor range and beyond. Furthermore, their distribution is traced back to the distribution of the contractive short-time Lyapunov exponent of the flow.     

Role of elastic stress in statistical and scaling properties of elastic turbulence

The role of elastic stress in statistical and scaling properties of elastic turbulence in a polymer solution flow between two disks is discussed. The analogy with a small-scale magnetodynamics and a passive scalar turbulent advection in the Batchelor regime is used to explain the experimentally observed statistical properties, the flow structure, and the scaling of elastic turbulence. The emergence of a new length scale, namely, the boundary layer thickness, is observed and studied.     

On continuation of inviscid vortex patches

This investigation concerns solutions of the steady-state Euler equations in two dimensions featuring finite area regions with constant vorticity embedded in a potential flow. Using elementary methods of the functional analysis we derive precise conditions under which such solutions can be uniquely continued with respect to their parameters, valid also in the presence of the Kutta condition concerning a fixed separation point. Our approach is based on the Implicit Function Theorem and perturbation equations derived using shape-differentiation methods. These theoretical results are illustrated with careful numerical computations carried out using the Steklov–Poincaré method which show the existence of a global manifold of solutions connecting the point vortex and the Prandtl–Batchelor solution, each of which satisfies the Kutta condition.     

Prandtl number effects in decaying homogeneous isotropic turbulence with a mean scalar gradient

Decaying homogeneous isotropic turbulence with an imposed mean scalar gradient is investigated numerically, thanks to a specific eddy-damped quasi-normal Markovian closure developed recently for passive scalar mixing in homogeneous anisotropic turbulence (BGC). The present modelling is compared successfully with recent direct numerical simulations and other models, for both very large and small Prandtl numbers. First, scalings for the cospectrum and scalar variance spectrum in the inertial range are recovered analytically and numerically. Then, at large Reynolds numbers, the decay and growth laws for the scalar variance and mixed velocity–scalar correlations, respectively, derived in BGC, are shown numerically to remain valid when the Prandtl number strongly departs from unity. Afterwards, the normalised correlation ρ_wθ is found to decrease in magnitude at a fixed Reynolds number when Pr either increases or decreases, in agreement with earlier predictions. Finally, the small scales return to isotropy of the scalar second-order moments is found to depend not only on the Reynolds number, but also on the Prandtl number.     

Soliton in discrete and continuum alpha helical proteins with higher-order excitations

We consider discrete and continuum models of alpha helical proteins with different orders of molecular excitations and interactions. In both the cases integrable models were identified for specific choice of parameters and the elementary excitations expressed in terms of solitons. Also, equivalent integrable Heisenberg spin systems were obtained.     

标量湍流的能量传输特性

应用直接数值模拟数据,从标量湍流传输的三波关系出发,进行湍流及标量湍流传输谱的多尺度分析,研究不同尺度间的能量传输性质,证实标量能量的传输与湍动能传输具有不同性质,大尺度速度脉动对标量传输有较大贡献,尤其是与标量小尺度脉动的相互作用,使标量模拟需要有比速度场更高的网格分辨率;并发现标量湍流的能量传输具有明显的非局部性;另外,定义了能量传输系数,发现在相同的Re数和Pe数条件下,标量湍流的对流惯性较速度脉动的惯性子区宽.     

Decay and growth laws in homogeneous shear turbulence

Homogeneous anisotropic turbulence has been widely studied in the past decades, both numerically and experimentally. Shear flows have received a particular attention because of the numerous physical phenomena they exhibit. In the present paper, both the decay and growth of anisotropy in homogeneous shear flows at high Reynolds numbers are revisited thanks to a recent eddy-damped quasi-normal Markovian closure adapted to homogeneous anisotropic turbulence. The emphasis is put on several aspects: an asymptotic model for the slow part of the pressure–strain tensor is derived for the return to isotropy process when mean velocity gradients are released. Then, a general decay law for purely anisotropic quantities in Batchelor turbulence is proposed. At last, a discussion is proposed to explain the scattering of global quantities obtained in DNS and experiments in sustained shear flows: the emphasis is put on the exponential growth rate of the kinetic energy and on the shear parameter.     

Basic Pattern in Atmospheric Turbulence Model

From the controlling equations of atmosphere motion, Prandtl's mixing length theory is used to derive the atmospheric turbulence models, such as Burgers equation model and Burgers-KdV equation model. And then the projective Riccati equations are applied to solve these atmospheric turbulence models, where much more patterns are obtained, including solitary wave pattern, singular pattern, and so on.     

Shell models of magnetohydrodynamic turbulence

Shell models of hydrodynamic turbulence originated in the seventies. Their main aim was to describe the statistics of homogeneous and isotropic turbulence in spectral space, using a simple set of ordinary differential equations. In the eighties, shell models of magnetohydrodynamic (MHD) turbulence emerged based on the same principles as their hydrodynamic counter-part but also incorporating interactions between magnetic and velocity fields. In recent years, significant improvements have been made such as the inclusion of non-local interactions and appropriate definitions for helicities. Though shell models cannot account for the spatial complexity of MHD turbulence, their dynamics are not over simplified and do reflect those of real MHD turbulence including intermittency or chaotic reversals of large-scale modes. Furthermore, these models use realistic values for dimensionless parameters (high kinetic and magnetic Reynolds numbers, low or high magnetic Prandtl number) allowing extended inertial range and accurate dissipation rate. Using modern computers it is difficult to attain an inertial range of three decades with direct numerical simulations, whereas eight are possible using shell models.     

Generalised scale-by-scale energy-budget equations and large-eddy simulations of anisotropic scalar turbulence at various Schmidt numbers

A necessary condition for the accurate prediction of turbulent flows using large-eddy simulation (LES) is the correct representation of energy transfer between the different scales of turbulence in the LES. For scalar turbulence, transfer of energy between turbulent length scales is described by a transport equation for the second moment of the scalar increment. For homogeneous isotropic turbulence, the underlying equation is the well-known Yaglom equation. In the present work, we study the turbulent mixing of a passive scalar with an imposed mean gradient by homogeneous isotropic turbulence. Both direct numerical simulations (DNS) and LES are performed for this configuration at various Schmidt numbers, ranging from 0.11 to 5.56. As the assumptions made in the derivation of the Yaglom equation are violated for the case considered here, a generalised Yaglom equation accounting for anisotropic effects, induced by the mean gradient, is derived in this work. This equation can be interpreted as a scale-by-scale energy-budget equation, as it relates at a certain scale r terms representing the production, turbulent transport, diffusive transport and dissipation of scalar energy. The equation is evaluated for the conducted DNS, followed by a discussion of physical effects present at different scales for various Schmidt numbers. For an analysis of the energy transfer in LES, a generalised Yaglom equation for the second moment of the filtered scalar increment is derived. In this equation, new terms appear due to the interaction between resolved and unresolved scales. In an a-priori test, this filtered energy-budget equation is evaluated by means of explicitly filtered DNS data. In addition, LES calculations of the same configuration are performed, and the energy budget as well as the different terms are thereby analysed in an a-posteriori test. It is shown that LES using an eddy viscosity model is able to fulfil the generalised filtered Yaglom equation for the present configuration. Further, the dependence of the terms appearing in the filtered energy-budget equation on varying Schmidt numbers is discussed.     

Irradiance scintillation index for Gaussian beam based on the generalized anisotropic von Karman spectrum

Experiments and theoretical investigations have shown that the atmosphere turbulence exhibits both anisotropic and non-Kolmogorov properties. In this paper, based on the anisotropic generalized von Karman spectrum and the Rytov approximation theory, new expression for the irradiance scintillation index of optical waves is derived for Gaussian beam propagating through weak anisotropic non-Kolmogorov turbulence. Compared with previously published results, it considers simultaneously the asymmetry property of turbulence cells or eddies in the orthogonal xy-plane, the general spectral power law in the range 3–4 instead of constant value of 11/3 for the Kolmogorov turbulence, and the finite turbulence inner and outer scales. Two anisotropic factors are introduced to parameterize the anisotropy of turbulence cells or eddies in horizontal and vertical directions. In the special cases of these two anisotropic factors equaling one and the finite turbulence inner and outer scales equaling separately zero and infinite, the derived expression can reduce correctly to the previously published results. Calculations are performed to analyze the derived results.     

Flavour dynamics with general scalar fields

We consider a spontaneously broken gauge theory based on the standard model (SM) group with scalar fields that carry arbitrary representations of G, and we investigate some general properties of the charged and neutral current involving these fields. In particular we derive the conditions for having real or complex couplings of the Z boson to two different neutral or charged scalar fields, and for the existence of CP-violating Z-scalar-scalar couplings. Moreover, we study models with the same fermion content as in the SM, with one SU(2) Higgs singlet, and an arbitrary number of Higgs doublets. We show that the structure of the Z-Higgs boson and of the Yukawa couplings in these models can be such that CP-violating form factors which conserve chirality are induced at the one-loop level. Received: 18 December 1998 / Published online: 22 March 1999     

Speckle statistical properties of Gaussian beam from a semi-rough target in the atmospheric turbulence

Based on extended Huygens–Fresnel principle and Goodman model for target surface, analytic expression is developed for the mutual coherent function (MCF) of a reflected Gaussian-beam from a semi-rough target in single-pass atmospheric turbulence. Then according to the MCF we derive expression about the mean intensity and average speckle size at the receiver. The analysis indicates that the mean intensity is closely related to the ratio of root mean square (rms) height to the lateral correlation length; in addition, the speckle size is associated with turbulence strength and roughness of target. Our results agree with well-known the limiting cases of perfectly smooth and Lambertian target.     

A priori evaluation of large eddy simulation subgrid-scale scalar flux models in isotropic passive-scalar and anisotropic buoyancy-driven homogeneous turbulence

Direct numerical simulation (DNS) of passive (non-buoyant) and active (buoyant) scalar homogeneous turbulence is carried out using a standard pseudo-spectral numerical method. The flow settings simulated include stationary forced and decaying passive-scalar turbulence, as well as decaying anisotropic active-scalar turbulence. The Schmidt number is unity in all cases. The results are compared with, and are found to be in very good agreement with, previous similar DNS studies. The well-validated DNS data are divided into 19 sets, and are employed to study different large eddy simulation (LES) subgrid-scale (SGS) models for the SGS scalar flux. The models examined include three eddy-viscosity-type models (Smagorinsky, Vreman and Sigma with a constant SGS Schmidt number), a Dynamic Structure model and two versions of the Gradient (Gradient and Modulated Gradient) model. The models are investigated with respect to their ability to predict the orientation, and the magnitude, of the SGS scalar flux. Eddy-viscosity models are found to predict the magnitude of the SGS scalar flux accurately, but are poor at predicting the orientation of the SGS scalar flux. The Dynamic Structure and Gradient models are better than eddy-viscosity models at predicting both the magnitude and direction. However, neither of them can be realised in an actual LES, without carrying additional transport equations. Based on these observations, four new models are proposed – combining directions from Dynamic Structure and Gradient models, and magnitudes from Smagorinsky and Vreman eddy-viscosity models. These models are expected to be better than eddy-viscosity and Modulated Gradient models, and this is confirmed by preliminary a posteriori tests.     

Propagation factor of electromagnetic concentric rings Schell-model beams in non-Kolmogorov turbulence

We derive an analytical expression for the propagation factor(known as M2-factor) of electromagnetic concentric rings Schell-model(EM CRSM) beams in non-Kolmogorov turbulence by utilizing the extended Huygens–Fresnel diffraction integral formula and the second-order moments of the Wigner distribution function(WDF). Our results show that the EM CRSM beam has advantage over the scalar CRSM beam for reducing the turbulence-induced degradation under suitable conditions. The EM CRSM beam with multi-rings far-fields in free space is less affected by the turbulence than the one with dark-hollow far-fields or the electromagnetic Gaussian Schell-model(EGSM) beam. The dependence of the M2-factor on the beam parameters and the turbulence are investigated in detail.     

A Closure for Isotropic Turbulence Based on Extended Scale Similarity Theory in Physical Space

The closure of a turbulence field is a longstanding fundamental problem, while most closure models are introduced in spectral space. Inspired by Chou's quasi-normal closure method in spectral space, we propose an analytical closure model for isotropic turbulence based on the extended scale similarity theory of the velocity structure function in physical space. The assumptions and certain approximations are justified with direct numerical simulation. The asymptotic scaling properties are reproduced by this new closure method, in comparison to the classical Batchelor model.     

On the dynamics of a continuum spin system

For a one-dimensional system of classical spins with nearest neighbour Heisenberg interaction we derive the equation of motion for each three-dimensional spin vector. In the continuum limit where the spins lie dense on a line this set of equations reduces to a nonlinear partial differential equation. In addition to spin-wave solutions we obtain some other special solutions of this equation. In particular we find solitary waves having total energy localised in a finite region, with velocity of propagation inversely proportional to the width of this region. Solutions of still another type are shown to have a diffusive character. The stability of such solutions and the possibility of interaction of two or more solitary waves have not yet been studied.