Eulerian–Lagrangian multiscale methods for solving scalar equations – Application to incompressible two-phase flows

The present article proposes a new hybrid Eulerian–Lagrangian numerical method, based on a volume particle meshing of the Eulerian grid, for solving transport equations. The approach, called Volume Of Fluid Sub-Mesh method (VOF-SM), has the advantage of being able to deal with interface tracking as well as advection–diffusion transport equations of scalar quantities. The Eulerian evolutions of a scalar field could be obtained on any orthogonal curvilinear grid thanks to the Lagrangian advection and a redistribution of particles on the Eulerian grid. The Eulerian concentrations result from the projection of the volume and scalar informations handled by the particles. The particle velocities are interpolated from the Eulerian velocity field. The VOF-SM method is validated on several scalar interface tracking and transport problems and is compared to existing schemes within the literature. It is finally coupled to a Navier–Stokes solver and applied to the simulation of two free-surface flows, i.e. the two-dimensional buckling of a viscous jet during the filling of a square mold and the three-dimensional dam-break flow in a tank.     

Intermittency of velocity time increments in turbulence

We analyze the statistics of turbulent velocity fluctuations in the time domain. Three cases are computed numerically and compared: (i) the time traces of Lagrangian fluid particles in a (3D) turbulent flow (referred to as the dynamic case); (ii) the time evolution of tracers advected by a frozen turbulent field (the static case); (iii) the evolution in time of the velocity recorded at a fixed location in an evolving Eulerian velocity field, as it would be measured by a local probe (referred to as the virtual probe case). We observe that the static case and the virtual probe cases share many properties with Eulerian velocity statistics. The dynamic (Lagrangian) case is clearly different; it bears the signature of the global dynamics of the flow.     

Lagrangian averages,averaged Lagrangians,and the mean effects of fluctuations in fluid dynamics

We begin by placing the generalized Lagrangian mean (GLM) equations for a compressible adiabatic fluid into the Euler-Poincare (EP) variational framework of fluid dynamics, for an averaged Lagrangian. This is the Lagrangian averaged Euler-Poincare (LAEP) theorem. Next, we derive a set of approximate small amplitude GLM equations (glm 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 glm equations uses the linearized relations between Eulerian and Lagrangian fluctuations, in the tradition of Lagrangian stability analysis for fluids. The glm derivation also uses the method of averaged Lagrangians, in the tradition of wave, mean flow interaction. Next, the new glm EP motion equations for incompressible ideal fluids are compared with the Euler-alpha turbulence closure equations. An alpha model is a GLM (or glm) fluid theory with a Taylor hypothesis closure. Such closures are based on the linearized fluctuation relations that determine the dynamics of the Lagrangian statistical quantities in the Euler-alpha equations. Thus, by using the LAEP theorem, we bridge between the GLM equations and the Euler-alpha closure equations, through the small-amplitude glm approximation in the EP variational framework. We conclude by highlighting a new application of the GLM, glm, and alpha-model results for Lagrangian averaged ideal magnetohydrodynamics. (c) 2002 American Institute of Physics.     

气管支气管树内气流和颗粒运动的大涡模拟

为更精确地描述真实人体呼吸道内的空气流动,明晰颗粒的运动沉积规律,本文从直接医学CT扫描得到的原始数据出发,利用图像辨识技术,重构了一个男性真人气管支气管树前三级的三维几何模型.采用大涡模拟的方法计算了非规则几何曲面结构内的气体流动现象,并在拉格朗日框架下跟踪颗粒的运动规律.数值计算得到了气流场的三维分布,以及颗粒的运动轨迹情况,结果表明现有基于Weibel的对称模型与真实人体的几何结构有较大的差异,而几何结构对流动影响较大;受非对称复杂结构影响,在不同截面的二次气流速度的分布规律不同;分叉后颗粒进入左右支气管的数量有明显的不同.     

Remapping-free ALE-type kinetic method for flow computations

Based on the integral form of the fluid dynamic equations, a finite volume kinetic scheme with arbitrary control volume and mesh velocity is developed. Different from the earlier unified moving mesh gas-kinetic method [C.Q. Jin, K. Xu, An unified moving grid gas-kinetic method in Eulerian space for viscous flow computation, J. Comput. Phys. 222 (2007) 155–175], the coupling of the fluid equations and geometrical conservation laws has been removed in order to make the scheme applicable for any quadrilateral or unstructured mesh rather than parallelogram in 2D case. Since a purely Lagrangian method is always associated with mesh entangling, in order to avoid computational collapsing in multidimensional flow simulation, the mesh velocity is constructed by considering both fluid velocity (Lagrangian methodology) and diffusive velocity (Regenerating Eulerian mesh function). Therefore, we obtain a generalized Arbitrary-Lagrangian–Eulerian (ALE) method by properly designing a mesh velocity instead of re-generating a new mesh after distortion. As a result, the remapping step to interpolate flow variables from old mesh to new mesh is avoided. The current method provides a general framework, which can be considered as a remapping-free ALE-type method. Since there is great freedom in choosing mesh velocity, in order to improve the accuracy and robustness of the method, the adaptive moving mesh method [H.Z. Tang, T. Tang, Adaptive mesh methods for one-and two-dimensional hyperbolic conservation laws, SIAM J. Numer. Anal. 41 (2003) 487–515] can be also used to construct a mesh velocity to concentrate mesh to regions with high flow gradients.     

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

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

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

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

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.     

Turbulent transport of material particles: an experimental study of finite size effects

We present experimental Lagrangian statistics of finite sized, neutrally bouyant, particles transported in an isotropic turbulent flow. The particle's diameter is varied over turbulent inertial scales. Finite size effects are shown not to be trivially related to velocity intermittency. The global shape of the particle's acceleration probability density functions is not found to depend significantly on its size while the particle's acceleration variance decreases as it becomes larger in quantitative agreement with the classical k(-7/3) scaling for the spectrum of Eulerian pressure fluctuations in the carrier flow.     

An immersed boundary energy-based method for incompressible viscoelasticity

Incompressible viscoelastic materials are prevalent in biological applications. In this paper we present a method for incompressible viscoelasticity in which the elasticity of the material is described in Lagrangian form (i.e. in material coordinates), and Eulerian (spatial) coordinates are used for the equations of motion and to enforce the incompressibility condition. The elastic forces are computed directly from an energy functional without the use of stress tensors, and the immersed boundary method is used to communicate between Lagrangian and Eulerian variables. The method is first applied to a warm-up problem, in which a viscoelastic incompressible material fills a two-dimensional periodic domain. For this problem, we study convergence of the velocity field, the deformation map, and the Eulerian force density. The numerical results indicate that the velocity field and deformation map converge strongly at second order and the Eulerian force density converges weakly at second order. Incompressibility is well maintained, as indicated by area conservation in this 2D problem. Finally, the method is applied to a three-dimensional fluid–structure interaction problem with two different materials: an isotropic neo-Hookean model and an anisotropic fiber-reinforced model.     

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.