scholarly journals A closure for Lagrangian velocity gradient evolution in turbulence using recent-deformation mapping of initially Gaussian fields

2016 ◽  
Vol 804 ◽  
pp. 387-419 ◽  
Author(s):  
Perry L. Johnson ◽  
Charles Meneveau

The statistics of the velocity gradient tensor in turbulent flows is of both theoretical and practical importance. The Lagrangian view provides a privileged perspective for studying the dynamics of turbulence in general, and of the velocity gradient tensor in particular. Stochastic models for the Lagrangian evolution of velocity gradients in isotropic turbulence, with closure models for the pressure Hessian and viscous Laplacian, have been shown to reproduce important features such as non-Gaussian probability distributions, skewness and vorticity strain-rate alignments. The recent fluid deformation (RFD) closure introduced the idea of mapping an isotropic Lagrangian pressure Hessian as the upstream initial condition using the fluid deformation tensor. Recent work on a Gaussian fields closure, however, has shown that even Gaussian isotropic velocity fields contain significant anisotropy for the conditional pressure Hessian tensor due to the inherent velocity–pressure couplings, and that assuming an isotropic pressure Hessian as the upstream condition may not be realistic. In this paper, Gaussian isotropic field statistics is used to generate more physical upstream conditions for the recent fluid deformation mapping. In this new framework, known isotropy relations can be satisfied by tuning the free model parameters and the original Gaussian field coefficients can be directly used without direct numerical simulation (DNS)-based re-adjustment. A detailed comparison of results from the new model, referred to as the recent deformation of Gaussian fields (RDGF) closure, with existing models and DNS shows the improvements gained, especially in various single-time statistics of the velocity gradient tensor at moderate Reynolds numbers. Application to arbitrarily high Reynolds numbers remains an open challenge for this type of model, however.

2020 ◽  
Vol 198 ◽  
pp. 104389 ◽  
Author(s):  
Ryosuke Nagata ◽  
Tomoaki Watanabe ◽  
Koji Nagata ◽  
Carlos B. da Silva

2018 ◽  
Vol 848 ◽  
pp. 876-905 ◽  
Author(s):  
Christopher J. Keylock

The velocity gradient tensor for turbulent flow contains crucial information on the topology of turbulence, vortex stretching and the dissipation of energy. A Schur decomposition of the velocity gradient tensor (VGT) is introduced to supplement the standard decomposition into rotation and strain tensors. Thus, the normal parts of the tensor (represented by the eigenvalues) are separated explicitly from non-normality. Using a direct numerical simulation of homogeneous isotropic turbulence, it is shown that the norm of the non-normal part of the tensor is of a similar magnitude to the normal part. It is common to examine the second and third invariants of the characteristic equation of the tensor simultaneously (the$\unicode[STIX]{x1D64C}{-}\unicode[STIX]{x1D64D}$diagram). With the Schur approach, the discriminant function separating real and complex eigenvalues of the VGT has an explicit form in terms of strain and enstrophy: where eigenvalues are all real, enstrophy arises from the non-normal term only. Re-deriving the evolution equations for enstrophy and total strain highlights the production of non-normality and interaction production (normal straining of non-normality). These cancel when considering the evolution of the VGT in terms of its eigenvalues but are important for the full dynamics. Their properties as a function of location in$\unicode[STIX]{x1D64C}{-}\unicode[STIX]{x1D64D}$space are characterized. The Schur framework is then used to explain two properties of the VGT: the preference to form disc-like rather than rod-like flow structures, and the vorticity vector and strain alignments. In both cases, non-normality is critical for explaining behaviour in vortical regions.


2010 ◽  
Vol 2010 ◽  
pp. 1-8 ◽  
Author(s):  
Luc Florack ◽  
Hans van Assen

Myocardial deformation and strain can be investigated using suitably encoded cine MRI that admits disambiguation of material motion. Practical limitations currently restrict the analysis to in-plane motion in cross-sections of the heart (2D + time), but the proposed method readily generalizes to 3D + time. We propose a new, promising methodology, which departs from a multiscale algorithm that exploits local scale selection so as to obtain a robust estimate for the velocity gradient tensor field. Time evolution of the deformation tensor is governed by a first-order ordinary differential equation, which is completely determined by this velocity gradient tensor field. We solve this matrix-ODE analytically and present results obtained from healthy volunteers as well as from patient data. The proposed method requires only off-the-shelf algorithms and is readily applicable to planar or volumetric tagging MRI sampled on arbitrary coordinate grids.


1999 ◽  
Vol 381 ◽  
pp. 141-174 ◽  
Author(s):  
ANDREW OOI ◽  
JESUS MARTIN ◽  
JULIO SORIA ◽  
M. S. CHONG

Since the availability of data from direct numerical simulation (DNS) of turbulence, researchers have utilized the joint PDFs of invariants of the velocity gradient tensor to study the geometry of small-scale motions of turbulence. However, the joint PDFs only give an instantaneous static representation of the properties of fluid particles and dynamical Lagrangian information cannot be extracted. In this paper, the Lagrangian evolution of the invariants of the velocity gradient tensor is studied using conditional mean trajectories (CMT). These CMT are derived using the concept of the conditional mean time rate of change of invariants calculated from a numerical simulation of isotropic turbulence. The study of the CMT in the invariant space (RA, QA) of the velocity-gradient tensor, invariant space (RS, QS) of the rate-of-strain tensor, and invariant space (RW, QW) of the rate-of-rotation tensor show that the mean evolution in the (Σ, QW) phase plane, where Σ is the vortex stretching, is cyclic with a characteristic period similar to that found by Martin et al. (1998) in the cyclic mean evolution of the CMT in the (RA, QA) phase plane. Conditional mean trajectories in the (Σ, QW) phase plane suggest that the initial reduction of QW in regions of high QW is due to viscous diffusion and that vorticity contraction only plays a secondary role subsequent to this initial decay. It is also found that in regions of the flow with small values of QW, the local values of QW do not begin to increase, even in the presence of self-stretching, until a certain self-stretching rate threshold is reached, i.e. when Σ≈0.25 〈QW〉1/2. This study also shows that in regions where the kinematic vorticity number (as defined by Truesdell 1954) is low, the local value of dissipation tends to increase in the mean as observed from a Lagrangian frame of reference. However, in regions where the kinematic vorticity number is high, the local value of enstrophy tends to decrease. From the CMT in the (−QS, RS phase plane, it is also deduced that for large values of dissipation, there is a tendency for fluid particles to evolve towards having a positive local value of the intermediate principal rate of strain.


Sign in / Sign up

Export Citation Format

Share Document