Averaged Lagrangians and the mean effects of fluctuations in ideal fluid dynamics

We begin by placing the generalized Lagrangian mean (GLM) equations for a compressible adiabatic fluid into the Euler–Poincaré (EP) variational framework of fluid dynamics, for an averaged Lagrangian. We then state the EP Averaging Result—that GLM equations arise from GLM Hamilton’s principles in the EP framework. Next, we derive a new set of approximate small-amplitude GLM equations (gℓm equations) at second order in the fluctuating displacement of a Lagrangian trajectory from its mean position. These equations express the linear and nonlinear back-reaction effects on the Eulerian mean fluid quantities by the fluctuating displacements of the Lagrangian trajectories in terms of their Eulerian second moments. The derivation of the gℓm equations uses the linearized relations between Eulerian and Lagrangian fluctuations, in the tradition of Lagrangian stability analysis for fluids. The gℓm derivation also uses the method of averaged Lagrangians, in the tradition of wave, mean flow interaction (WMFI). The gℓm EP motion equations for compressible and incompressible ideal fluids are compared with the Euler-alpha turbulence closure equations. An alpha model is a GLM (or gℓm) fluid theory with a Taylor hypothesis closure (THC). Such closures are based on the linearized fluctuation relations that determine the dynamics of the Lagrangian statistical quantities in the Euler-alpha closure equations. We use the EP Averaging Result to bridge between the GLM equations and the Euler-alpha closure equations. Hence, combining the small-amplitude approximation with THC yields in new turbulence closure equations for compressible fluids in the EP variational framework.     

Averaged Lagrangians and MacCallum-Taub's limit in macroscopic gravity

The derivation of MacCallum-Taub's averaged high-frequency Lagrangian [1] is analysed with special attention paid to the assumptions made along the derivation. It is shown that averaged high-frequency Lagrangians of the same form as MacCallum-Taub's Lagrangian can be derived by applying the Brill-Hartle and macroscopic gravity averaging schemes. A procedure for the derivation of a Lagrangian of macroscopic gravity (an averaged Hilbert action) is proposed and its high-frequency limit (namely, its high-frequency perturbation expansion up to the second order terms in perturbations, which is referred to as MacCallum-Taub's limit) is calculated. There is disagreement [2] in the expressions for MacCallum-Taub's averaged high-frequency Lagrangian and the high-frequency limit of the macroscopic gravity Lagrangian. Possible reasons for such disagreement are analysed. The origin of the difference is shown to consist in using the propagation equation for perturbations, i.e. the linearized Ricci tensor vanishes, during the derivation (averaging) carried out in [1]. A new derivation of an averaged high-frequency Lagrangian without assuming the propagation equation to hold and by taking into account the proper correlation functions is given. The newly derived expression is shown to coincide with MacCallum-Taub's limit of the macroscopic gravity Lagrangian, which resolves the disagreement.     

Two-fluid equations for phonons without momentum

In the customary microscopic derivation of the two-fluid equations of liquid He II explicit use is made of the assumption that an elementary excitation of wave vector k carries a momentum ?k. In this paper it is shown that phonons in a liquid can be defined as carrying momentum ?k (Eulerian phonons) or zero momentum (Lagrangian phonons). A careful analysis—in particular of the concept “velocity of the bearer fluid”—shows that the two-fluid equations turn out to be the same in either case.     

ALE计算中基于声速分布的网格优化技术

基于声速分布,提出一种拉氏流体力学计算中大变形网格优化的数值技术.该方法不但可以优化网格的几何形状且可以提高拉氏流体计算的时间步长.介绍基于声速分布的网格松弛泛函、修正梯度流方程的推导、离散和求解方法,启动/终止网格优化过程的条件,及基于这种网格优化方法的ALE算法.给出Rayleigh-Taylor不稳定性问题等数值算例,用以证明该方法的有效性.     

Geometric gradient-flow dynamics with singular solutions

The gradient-flow dynamics of an arbitrary geometric quantity is derived using a generalization of Darcy’s Law. We consider flows in both Lagrangian and Eulerian formulations. The Lagrangian formulation includes a dissipative modification of fluid mechanics. Eulerian equations for self-organization of scalars, 1-forms and 2-forms are shown to reduce to nonlocal characteristic equations. We identify singular solutions of these equations corresponding to collapsed (clumped) states and discuss their evolution.     

湍流两相流动有燃烧颗粒相概率密度函数输运方程理论

由有燃烧的湍流气粒两相流动的瞬态方程和统计力学概率密度函数概念出发,推导了有燃烧颗粒相的质量－动量－能量联合概率密度函数（ＰＤＦ）输运方程,并对方程中条件期望项用梯度模拟概念进行了模拟封闭。封闭后的ＰＤＦ方程可作为建立颗粒拟流体模型方程和封闭二阶矩模型的基础,也可以通过Ｍｏｎｔｅ－Ｃａｒｌｏ 法求解用以直接计算颗粒雷诺应力和湍流动能,以便和二阶 矩模型的结果相对照,改善二阶矩模型。     

Numerical simulations of instabilities in the implosion process of inertial confined fusion in 2D cylindrical coordinates

In this paper, we introduce a multi-material arbitrary Lagrangian and Eulerian method for the hydrodynamic radiative multi-group diffusion model in 2D cylindrical coordinates. The basic idea in the construction of the method is the following: In the Lagrangian step, a closure model of radiation-hydrodynamics is used to give the states of equations for materials in mixed cells. In the mesh rezoning step, we couple the rezoning principle with the Lagrangian interface tracking method and an Eulerian interface capturing scheme to compute interfaces sharply according to their deformation and to keep cells in good geometric quality. In the interface reconstruction step, a dual-material Moment-of-Fluid method is introduced to obtain the unique interface in mixed cells. In the remapping step, a conservative remapping algorithm of conserved quantities is presented. A number of numerical tests are carried out and the numerical results show that the new method can simulate instabilities in complex fluid field under large deformation,and are accurate and robust.     

The shallow water equations in Lagrangian coordinates

Recent advances in the collection of Lagrangian data from the ocean and results about the well-posedness of the primitive equations have led to a renewed interest in solving flow equations in Lagrangian coordinates. We do not take the view that solving in Lagrangian coordinates equates to solving on a moving grid that can become twisted or distorted. Rather, the grid in Lagrangian coordinates represents the initial position of particles, and it does not change with time. We apply numerical methods traditionally used to solve differential equations in Eulerian coordinates, to solve the shallow water equations in Lagrangian coordinates. The difficulty with solving in Lagrangian coordinates is that the transformation from Eulerian coordinates results in solving a highly nonlinear partial differential equation. The non-linearity is mainly due to the Jacobian of the coordinate transformation, which is a precise record of how the particles are rotated and stretched. The inverse Jacobian must be calculated, thus Lagrangian coordinates cannot be used in instances where the Jacobian vanishes. For linear (spatial) flows we give an explicit formula for the Jacobian and describe the two situations where the Lagrangian shallow water equations cannot be used because either the Jacobian vanishes or the shallow water assumption is violated. We also prove that linear (in space) steady state solutions of the Lagrangian shallow water equations have Jacobian equal to one. In the situations where the shallow water equations can be solved in Lagrangian coordinates, accurate numerical solutions are found with finite differences, the Chebyshev pseudospectral method, and the fourth order Runge–Kutta method. The numerical results shown here emphasize the need for high order temporal approximations for long time integrations.     

Variational symmetries and Lagrangian multiforms

Letters in Mathematical Physics - By considering the closure property of a Lagrangian multiform as a conservation law, we use Noether’s theorem to show that every variational symmetry of a...     

An Improved Method to Determine Particle Dispersion Width for Efficient Modeling of Turbulent Two‐Phase Flows

An improved approach is presented for the hybrid Eulerian‐Lagrangian modeling of turbulent two‐phase flows. The hybrid model consists of a nonlinear k–ε model for the fluid flow and an efficient Lagrangian trajectory model for the particulate flow. The improved approach avoids an empirical correlation required to determine the dispersion width for the existing Stochastic‐Probabilistic Efficiency Enhanced Dispersion (SPEED) model. The improved SPEED model is validated using experimental data for a poly‐dispersed water spray interacting with a turbulent annular air jet behind a bluff‐body. Numerical results for the number‐mean and Sauter‐mean droplet diameters, as well as mean and fluctuating droplet velocities are compared with the experimental data and with the predictions of other dispersion models. It is demonstrated that higher computational efficiency and smoother profiles of Sauter‐mean diameter can be obtained with the improved stochastic‐probabilistic model than with the eddy‐interaction model.     

非旋转性前进波Euler与Lagrange解间的转换

对等深水中非旋转性的前进重力波动场,以求得的Euler与Lagrange两种形式至第三阶的解,按照同一流体质点在相同时间与位置处其流速唯一与质量守恒性及在自由表面水位处Euler形式解与Lagrange形式解为同一值的特性,来推导二者可相互转换.由连续的Taylor级数展开,考虑波动场中各流体质点的运动轨迹与运动周期,将已知的Euler形式解转换成完全未知的Lagrange形式解,解决了以往成果中出现含时间的不合理的共振项,以及无法得到与Euler系统不同的Lagrange形式的流体质点运动频率与平均运动 关键词： 非旋转性前进波 Euler-Lagrange转换 质点运动轨迹 质点运动频率     

Euler's fluid equations: Optimal control vs optimization

An optimization method used in image-processing (metamorphosis) is found to imply Euler's equations for incompressible flow of an inviscid fluid, without requiring that the Lagrangian particle labels exactly follow the flow lines of the Eulerian velocity vector field. Thus, an optimal control problem and an optimization problem for incompressible ideal fluid flow both yield the same Euler fluid equations, although their Lagrangian parcel dynamics are different. This is a result of the gauge freedom in the definition of the fluid pressure for an incompressible flow, in combination with the symmetry of fluid dynamics under relabeling of their Lagrangian coordinates. Similar ideas are also illustrated for SO(N) rigid body motion.     

钝体后湍流预混燃烧的PDF模拟

本文采用PDF方法对矩形燃烧室内钝体后的湍流预混火焰进行了数值模拟。脉动速度-频率-标量联合的PDF输运方程用Monte Carlo方法求解,质量、动量和能量的平均值由基于无结构网格的有限体积法求解,压力通过状态方程获得。PDF方程中所需的平均密度、平均速度和压力由有限体积法提供,并将用Monte Carlo方法求出的雷诺应力、化学反应源项和比热比传递给有限体积法。本文对丙烷和空气燃烧的不同简化化学反应机理进行了研究,并与实验结果进行比较,获得满意的结果。     

Closure Approximations for Passive Scalar Turbulence: A Comparative Study on an Exactly Solvable Model with Complex Features

Some standard closure approximations used in turbulence theory are analyzed by examining systematically the predictions these approximations produce for a passive scalar advection model consisting of a shear flow with a fluctuating cross sweep. This model has a general geometric structure of a jet flow with transverse disturbances, which occur in a number of contexts, and it encompasses a wide variety of possible spatio-temporal statistical structures for the velocity field, including strong long-range correlations. Even though the Eulerian and Lagrangian velocity statistics are not equal and the passive scalar statistics exhibit broader-than-Gaussian intermittency, this model is nevertheless simple enough so that many passive scalar statistics can be computed exactly and compared systematically with the predictions of the closure approximations. Our comparative study illustrates the strength and weaknesses of the closure approximations and points out the physical phenomena that these approximations are able or not able to describe properly. In particular it is shown that the direct interaction approximation (DIA), one of the most sophisticated closure approximations available, fails to reproduce adequately the statistical features of the scalar and may even lead to absdurd predictions, even though the equations it produces are rather complicated and difficult to analyze. Two alternative closure approximations, the Modified DIA (MDIA) and the Renormalized Lagrangian Approximation (RLA), with different levels of sophistication, both are simpler to use than the DIA and perform better. In particular, it is shown that both closure approximations always reproduce exactly the second order statistics for the scalar and that the MDIA is even able to capture intermittency effects.     

Hydrodynamics, Particle Relabelling and Relativity

Using the wave equation as an example, it is shown how to extend the hydrodynamic Lagrangian-picture method of constructing field evolution using a continuum of trajectories to second-order theories. The wave equation is represented through Eulerian-picture models that are distinguished by their Lorentz transformation properties. Introducing the idea of the relativity of the particle label, it is demonstrated how the corresponding trajectory models are compatible with the relativity principle. It is also shown how the Eulerian variational formulation may be obtained by canonical transformation from the Lagrangian picture, and how symmetries in the Lagrangian picture may be used to generate Eulerian conserved charges.     

非旋转性自由表面重力驻波的Euler与Lagrange方法解间的转换

对于等深水中的非旋转性重力驻波流场,本文用Euler与Lagrange两种方法求得其至三阶的解,根据同一粒流体质点在相同时间与位置处其流速值为唯一与质量守恒及在自由表面水位的Euler形式解与Lagrange形式解相同等特性,来推导其间互可转换.由一系列连续的Taylor级数展开,在考虑波动场中各流体质点的运动轨迹与运动周期条件下,将已知的Euler解转换成完全未知的Lagrange形式解.接着再将所得的Lagrange解转换成对应的Euler形式,均可得到完全相同的结果.由此可得知,在考虑波动场各流体质 关键词： 重力驻波 Euler与Lagrange解间的转换 质点运动轨迹