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.

The evolution of the invariants (R and Q) of the velocity gradient tensor in homogeneous isotropic turbulence is investigated using data from direct numerical simulation (DNS). The concepts of conditional average time rate of change of the invariants and conditional mean trajectories (CMT) in invariant phase space are introduced to study the dynamics of this flow. The resulting dynamical system in the (R,Q) phase space is a clockwise spiral with a stable focus at the origin, illustrating that in the mean, the cyclic sequence of topological evolution following a fluid particle is unstable-node/saddle/saddle (UN/S/S)→stable-node/saddle/saddle (SN/S/S)→stable-focus/stretching (SF/S)→unstable-focus/contracting (UF/C). The mean rates of change of R and Q, i.e., Ṙ, Q̇, are found to be negligible near the right branch of the null discriminant (D=0) curve, indicating that this curve is an attractor in the (R,Q) space. The effects of both the diffusion term and the anisotropic part of the pressure Hessian term on the dynamics of the invariants have also been analyzed using the conditional averages. Both contributions are found to be important in the dynamics of the velocity gradient invariants. Based on these results the extent of the validity of the model equations governing the evolution of R and Q proposed by Cantwell [Phys. Fluids A 4, 782 (1992)] and Dopazo et al. [“Velocity gradients in turbulent flows. Stochastic models,” Ninth Symposium on “Turbulent Shear Flows,” Kyoto, Japan, 1993, pp. 26-2-1–26-2-5] are discussed.

No abstract

A three-dimensional direct numerical simulation of a propagating turbulent premixed flame is performed using one-step Arrhenius model chemistry. The interaction of the flame thermochemical processes with the local geometries of the scalar field and flow topologies is studied. Four regions (“fresh reactants,” “preheating,” “burning,” and “hot products”), characterized by their reaction rate and mass fraction values, are examined. Thermochemical processes in the “preheating” and “burning” regions smooth out highly contorted iso-scalar surfaces, present in the “fresh reactants,” and annihilate large curvatures. Positive volumetric dilatation rates, −P = ∇ · u, display maxima for elliptic concave and minima for convex scalar micro-structures. Constant average tangential strain rates, aT, with large fluctuations, occur throughout the flow domain, whereas normal strain rates, aN, follow the trends of volumetric dilatation rates. Focal topologies, present in the “fresh reactants,” tend to disappear in favor of nodal structures as moving towards the “hot products.” The vorticity vector is predominantly tangential to the iso-scalar surfaces. The Unstable Node/Saddle/Saddle and Stable Focus/Stretching topologies, present in the “fresh reactants,” correlate with large values of aN and aT providing hints on the flow topologies fostering scalar mixing.