Physical mechanism of the two-dimensional inverse energy cascade

We study the physical mechanisms of the two-dimensional inverse energy cascade using theory, numerics, and experiment. Kraichnan's prediction of a -5/3 spectrum with constant, negative energy flux is verified in our simulations of 2D Navier-Stokes equations. We observe a similar but shorter range of inverse cascade in laboratory experiments. Our theory predicts, and the data confirm, that inverse cascade results mainly from turbulent stress proportional to small-scale strain rotated by 45 degrees. This "skew-Newtonian" stress is explained by the elongation and thinning of small-scale vortices by large-scale strain which weakens their velocity and transfers their energy upscale.     

Intermittency and non-Gaussian fluctuations of the global energy transfer in fully developed turbulence

We address the experimentally observed non-Gaussian fluctuations for the energy injected into a closed turbulent flow at fixed Reynolds number. We propose that the power fluctuations mirror the internal kinetic energy fluctuations. Using a stochastic cascade model, we construct the excess kinetic energy as the sum over the energy transfers at different levels of the cascade. We find an asymmetric distribution that strongly resembles the experimental data. The asymmetry is an explicit consequence of intermittency and the global measure is dominated by small scale events correlated over the entire system. Our calculation is consistent with the statistical analogy recently made between a confined turbulent flow and a critical system of finite size.     

Inverse energy cascade in three-dimensional isotropic turbulence

We study the statistical properties of homogeneous and isotropic three-dimensional (3D) turbulent flows. By introducing a novel way to make numerical investigations of Navier-Stokes equations, we show that all 3D flows in nature possess a subset of nonlinear evolution leading to a reverse energy transfer: from small to large scales. Up to now, such an inverse cascade was only observed in flows under strong rotation and in quasi-two-dimensional geometries under strong confinement. We show here that energy flux is always reversed when mirror symmetry is broken, leading to a distribution of helicity in the system with a well-defined sign at all wave numbers. Our findings broaden the range of flows where the inverse energy cascade may be detected and rationalize the role played by helicity in the energy transfer process, showing that both 2D and 3D properties naturally coexist in all flows in nature. The unconventional numerical methodology here proposed, based on a Galerkin decimation of helical Fourier modes, paves the road for future studies on the influence of helicity on small-scale intermittency and the nature of the nonlinear interaction in magnetohydrodynamics.     

A theoretical model of turbulent fiber suspension and its application to the channel flow

A theoretical model of turbulent fiber suspension is developed by deriving the equations of Reynolds averaged Navier-Stokes,turbulence kinetic energy and turbulence dissipation rate with the additional term of fibers.In order to close the above equations,the equation of probability distribution function for mean fiber orientation is also derived.The theoretical model is applied to the turbulent channel flow and the corresponding equations are solved numerically.The numerical results are verified by comparisons with the experimental ones.The effects of Reynolds number,fiber concentration and fiber aspect-ratio on the velocity profile,turbulent kinetic energy and turbulent dissipation rate are analyzed.Based on the numerical data,the expression for the velocity profile in the turbulent fiber suspension channel flow,which includes the effect of Reynolds number,fiber concentration and aspect-ratio,is proposed.     

Differential model for 2D turbulence

We present a phenomenological model for 2D turbulence in which the energy spectrum obeys a nonlinear fourth-order differential equation. This equation respects the scaling properties of the original Navier-Stokes equations, and it has both the −5/3 inverse-cascade and the −3 direct-cascade spectra. In addition, our model has Raleigh-Jeans thermodynamic distributions as exact steady state solutions. We use the model to derive a relation between the direct-cascade and the inverse-cascade Kolmogorov constants, which is in good qualitative agreement with the laboratory and numerical experiments. We discuss a steady state solution where both the enstrophy and the energy cascades are present simultaneously, and we discuss it in the context of the Nastrom-Gage spectrum observed in atmospheric turbulence. We also consider the effect of the bottom friction on the cascade solutions and show that it leads to an additional decrease and finite-wavenumber cutoffs of the respective cascade spectra, which agrees with the existing experimental and numerical results. The text was submitted by the authors in English.     

Quasi-3D Navier–Stokes Model for a Rotating Airfoil

A quasi-3D model of the unsteady Navier–Stokes equations in a rotating frame of reference has been developed. The equations governing the flow past a rotating blade are approximated using an order of magnitude analysis on the spanwise derivatives. The model takes into account rotational effects and spanwise outflow at computing expenses in the order of what is typical for similar 2D calculations. Results are presented for both laminar and turbulent flows past blades in pure rotation. In the turbulent case the influence of small-scale turbulence is modelled by the one-equation Baldwin–Barth turbulence model. The computations demonstrate that the main influence of rotation is to increase the maximum lift.     

泡状流运动的湍流模型及液相速度和湍动能的分布预测

１前言在核能、航天等技术领域和能源、动力、石油化工等工业过程都存在泡状流现象,其中绝大多数的泡状流流动状态为湍流,因此对湍流泡状流进行深入研究极为必要。近年来Ｌｅｅ［１］、Ｌｏｐｅｚ　ｄｅ　Ｂｅｒｔｏｄａｎｏ［２,３］等提出了一些多维的湍流输运的双流...     

Local Collocation Approach for Solving Turbulent Combined Forced and Natural Convection Problems

An application of the meshless Local Radial Basis Function Collocation Method (LRBFCM) [22, 30–33] in solution of incompressible turbulent combined forced and natural convection is for the first time explored in the present paper. The turbulent flow equations are described by the low-Re number $k − ε$ model with Launder and Sharma [23] and Abe et al. [1] closure coefficients. The involved temperature, velocity, pressure, turbulent kinetic energy and dissipation fields are represented on overlapping 5-noded sub-domains through the collocation by using multiquadrics Radial Basis Functions (RBF). The involved first and second order partial derivatives of the fields are calculated from the respective derivatives of the RBF's. The involved equations are solved through the explicit time stepping. The pressure-velocity coupling is based on Chorin's fractional step method [11]. The adaptive upwinding technique, proposed by Lin and Atluri [27], is used because of the convection dominated situation. The solution procedure is represented for a 2D upward channel flow with differentially heated walls. The results have been assessed by achieving a reasonable agreement with the direct numerical simulation of Kasagi and Nishimura [20] for Reynolds number 4494, based on the channel width, and Grashof number 9.6 × 10⁵. The advantages of the represented mesh-free approach are its simplicity, accuracy, similar coding in 2D and 3D, and straightforward applicability in non-uniform node arrangements.     

Energy spectrum of quasigeostrophic turbulence

We consider the energy spectrum of a quasigeostrophic model of forced, rotating turbulent flow. We provide a rigorous a priori bound E(k)相似文献   

Comparison of sound power radiation from isolated airfoils and cascades in a turbulent flow

An analytical model of the sound power radiated from a flat plate airfoil of infinite span in a 2D turbulent flow is presented. The effects of stagger angle on the radiated sound power are included so that the sound power radiated upstream and downstream relative to the fan axis can be predicted. Closed-form asymptotic expressions, valid at low and high frequencies, are provided for the upstream, downstream, and total sound power. A study of the effects of chord length on the total sound power at all reduced frequencies is presented. Excellent agreement for frequencies above a critical frequency is shown between the fast analytical isolated airfoil model presented in this paper and an existing, computationally demanding, cascade model, in which the unsteady loading of the cascade is computed numerically. Reasonable agreement is also observed at low frequencies for low solidity cascade configurations.     

Turbulence: Does Energy Cascade Exist?

To answer the question whether a cascade of energy exists or not in turbulence, we propose a set of correlation functions able to test if there is an irreversible transfert of energy, step by step, from large to small structures. These tests are applied to real Eulerian data of a turbulent velocity flow, taken in the wind grid tunnel of Modane, and also to a prototype model equation for wave turbulence. First we demonstrate the irreversible character of the flow by using multi-time correlation function at a given point of space. Moreover the unexpected behavior of the test function leads us to connect irreversibility and finite time singularities (intermittency). Secondly we show that turbulent cascade exists, and is a dynamical process, by using a test function depending on time and frequency. The cascade shows up only in the inertial domain where the kinetic energy is transferred more rapidly (on average) from the wavenumber \(k_{1}\) to \(k_{2}\) than from \(k_{1}\) to \(k'_{2}\) larger than \(k_{2}\).     

Turbulence decay rate as a measure of flow dimensionality

The dimensionality of turbulence in fluid layers determines their properties. We study electromagnetically driven flows in finite-depth fluid layers and show that eddy viscosity, which appears as a result of three-dimensional motions, leads to increased bottom damping. The anomaly coefficient, which characterizes the deviation of damping from the one derived using a quasi-two-dimensional model, can be used as a measure of the flow dimensionality. Experiments in turbulent layers show that when the anomaly coefficient becomes high, the turbulent inverse energy cascade is suppressed. In the opposite limit turbulence can self-organize into a coherent flow.     

Vortex burst as a source of turbulence

An important issue in turbulence theory is to understand what kinds of elementary flow structures are responsible for the part of the turbulent energy spectrum described by Kolmogorov's celebrated k(-5/3) law. A model for such structure has been proposed by Lundgren [Phys. Fluids 25, 2193-2203 (1982)]] in the form of a vortex with spiral structure subjected to an axially straining field. We report experimental results of a vortex burst in a laminar-flow environment showing that this structure is responsible for a k(-5/3) part in the energy spectrum. If there are many experimental evidences of the existence of vortices with spiral structures in turbulent flows, it is the first time that such an elementary structure is experimentally shown to be responsible for the turbulent energy cascade.     

有液滴卷吸时的环状流换热系数的预测模型

本文对垂直上升光管中环状流流动沸腾的理论模型进行了分析,以液膜紊流的动量方程和能量方程为基础,推导了环状流的数学模型,通过求解动量方程和能量方程,获得了流动沸腾换热系数的预测模型,并对该预测模型进行了数值求解,将预测的换热系数同实验值作了比较。     

Grid and subgrid-scale interactions in viscoelastic turbulent flow and implications for modelling

Using direct numerical simulations of turbulent plane channel flow of homogeneous polymer solutions, described by the Finitely Extensible Nonlinear Elastic-Peterlin (FENE-P) rheological constitutive model, a-priori analyses of the filtered momentum and FENE-P constitutive equations are performed. The influence of the polymer additives on the subgrid-scale (SGS) energy is evaluated by comparing the Newtonian and the viscoelastic flows, and a severe suppression of SGS stresses and energy is observed in the viscoelastic flow. All the terms of the transport equation of the SGS kinetic energy for FENE-P fluids are analysed, and an approximated version of this equation for use in future large eddy simulation closures is suggested. The terms responsible for kinetic energy transfer between grid-scale (GS) and SGS energy (split into forward/backward energy transfer) are evaluated in the presence of polymers. It is observed that the probability and intensity of forward scatter events tend to decrease in the presence of polymers.     

Numerical simulation of pulsation processes in hydraulic turbine based on 3D model of cavitating flow

A new approach was proposed for simulation of unsteady cavitating flow in the flow passage of a hydraulic power plant. 1D hydro-acoustics equations are solved in the penstock domain. 3D equations of turbulent flow of isothermal compressible liquid-vapor mixture are solved in the turbine domain. Cavitation is described by a transfer equation for liquid phase with a source term which is responsible for evaporation and condensation. The developed method was applied for simulation of pulsations in pressure, discharge, and total energy propagating along the flow conduit of the hydraulic power plant. Simulation results are in qualitative and quantitative agreement with experiment. The influence of key physical and numerical parameters like discharge, cavitation number, penstock length, time step, and vapor density on simulation results was studied.     

Turbulence of incompressible liquid flow near local equilibrium and the principle of secondary maximum of entropy

The variational principle of maximum entropy is used to describe the dynamics of weakly nonequilibrium turbulence using the theory of Reynolds stresses for viscous incompressible liquid flow. From this principle, equations closing the theory of Reynolds stresses and also equations describing mean flow-turbulence interaction for 3D turbulent flows are derived. The theory is reduced to 2D flows and weak turbulence. Thermodynamic analogues and an example of Couette flow are considered.     

不可压湍流的控制方程

使用侧偏统计平均方法理性推导了不可压湍流的封闭方程组．侧偏统计平均提供了湍流脉动量的一阶统计信息，在引入加权漂移速度对称性及正交各向异性后，使用动量传输链概念模化封闭了整个方程组．方程组不含任何经验系数，不使用壁面函数，保留了ＮＳ方程的均化的非线性特性．其级数形式的能量方程与非线性现象多尺度层次现象相对应，具备了描述湍流统计平均流动及拟序结构流动的双重功能．用平面射流，圆射流，层流－湍流转捩及后向台阶流等算例初步验证了方程对真实湍流的适定性．     

A causal continuous-time stochastic model for the turbulent energy cascade in a helium jet flow

We discuss continuous cascade models and their potential for modelling the energy dissipation in a turbulent flow. Continuous cascade processes, expressed in terms of stochastic integrals with respect to Lévy bases, are examples of ambit processes. These models are known to reproduce experimentally observed properties of turbulence: the scaling and self-scaling of the correlators of the energy dissipation and of the moments of the coarse-grained energy dissipation. We compare three models: a normal model, a normal inverse Gaussian model, and a stable model. We show that the normal inverse Gaussian model is superior to both, the normal and the stable models, in terms of reproducing the distribution of the energy dissipation; and that the normal inverse Gaussian model is superior to the normal model and competitive with the stable model in terms of reproducing the self-scaling exponents. Furthermore, we show that the presented analysis is parsimonious in the sense that the self-scaling exponents are predicted from the one-point distribution of the energy dissipation, and that the shape of these distributions is independent of the Reynolds number.     

Near Continuum Velocity and Temperature Coupled Compressible Boundary Layer Flow over a Flat Plate

The problem of a compressible gas flows over a flat plate with the velocity-slip and temperature-jump boundary conditions are being studied. The standard single- shooting method is applied to obtain the exact solutions for velocity and temperature profiles when the momentum and energy equations are weakly coupled. A double-shooting method is applied if these two equations are closely coupled. If the temperature affects the velocity directly, more significant velocity slip happens at locations closer to the plate’s leading edge, and inflections on the velocity profiles appear, indicating flows may become unstable. As a consequence, the temperature-jump and velocity-slip boundary conditions may trigger earlier flow transitions from a laminar to a turbulent flow state.