Lipo Wang Institut für Technische Verbrennung RWTH-Aachen, Germany TMB-2009, Trieste 2009.07.28

1. Local and nonlocal conditional strain rates along gradient trajectories from various scalar fields in turbulence Lipo Wang Institut für Technische Verbrennung RWTH-Aachen, Germany TMB-2009, Trieste 2009.07.28

2. Content • Content • Introduction: background and general analysis • Local statistics • strain rate along scalar gradients • interaction of different gradient trajectories • Non-local statistics • a new two-point velocity difference structure function • Conclusions

3. Background Mathematically, dependent variables to describe turbulence can be: (a) in scalar form: the passive scalar , kinetic energy … (b) in vector form: the velocity vector , the vorticity vector … (c) In tensor form: rate-of-strain, the stress tensor … Only for those in (a), their presentations are independent of the choice of coordinate systems. In the Cartesian coordinate system, vector and tensor variables consist of projections (along constant directions) in scalar form. Scalar field variables: variables in scalar form + spatial projections of vector variables.

4. Background The scalars of interest and importance: the passive scalar, projections of the velocity and vorticity vector, kinetic energy and energy dissipation. Question: are there any generic properties or relations among these different scalars?

5. Background • Hints from previous studies: • Tennekes & Lumley (1972): turbulence transports passive contaminants in much the same way as momentum. • Tsinober (2001): under the action of strain tensor, local statistics of the passive vector behaves qualitatively the same as the active ones. • Wang & Peters (2008): dissipation elements from different scalar fields evolve in a similar fashion.

6. General analysis Formally the governing equation is: For the momentum equations ( ) : at large enough Reynolds numbers, the viscous term and pressure gradient term are relative small; the convection term and time-dependent term approximately balance.

7. General analysis Observation: in if instantaneously the magnitude of is relatively small (at high enough Re), different scalars may assume some generalities. • for , , therefore instantaneously small at high Re; • for , at high Re, both the pressure gradient and  are small; • large  tends to be perpendicular to the directions of largest (positive and negative) strain rates. The product of  and ‘could’ be small?

8. Local statistics Object parameter(s): for different scalars, their convection terms contain the same (turbulent) velocity; thus the velocity-related properties may statistically be similar: Laminar diffusion analysis by Batchelor (1959): locally under the action of a constant velocity field (x,y,z), the governing equation of the passive scalar becomes Thus the gradient vector evolves from any arbitrary orientation to Ashurst et al. (1987): the passive vector locally tends to align with the most compressive strain direction.

9. Local statistics Numerical verification from DNS (homogenous shear flow) A finer resolution of x/<1 is important for diagnosing fine structures.

10. Local statistics PDF of (the passive scalar) Case 2: PDF of

11. Local statistics DNS result: the mean strain rate from different scalar gradients is an index showing the alignment between the scalar gradient vector and the principle directions of the strain rate. Question: how will the orientation of different scalar gradient vectors behave?

12. Local statistics Orientation conditioned on the magnitude of gradients P1: from small gradient points P2: from large gradient points

13. Local statistics The sixth moment of the orientation PDF. case 1 case 2

14. Local statistics Conclusion: in regions of large scalar gradients, gradient vectors tend to align with each other, while if gradients are of small magnitudes, they orientate irrelevantly. Physically, large scalar gradients are due to compressive strain.

15. Local statistics Illustration of the interaction of various gradient trajectories.

16. Nonlocal statistics Local statistics, although informative, are not enough to describe turbulent fields. There is a strong need to develop knowledge about the nonlocal properties. Existing theories: Karman-Howarth equation; Kolmogorov’s structure function… Two-point structure function conditioned along gradient trajectories?

17. Nonlocal statistics: theory The passive scalar equation: The passive gradient equation: Along a same gradient trajectory, the two-point correlation equation of the scalar gradient:

18. Nonlocal statistics: theory The scalar gradient correlation in the Cartesian frame: Normalized from:

19. Nonlocal statistics: theory The scalar gradient correlation along gradient trajectories: Normalized from: In the inertial range:

20. Nonlocal statistics: theory From the two-point scalar gradient correlation equation we obtain (a) in the inertial range: (b) In the viscous diffusive range:

21. a s b Nonlocal statistics: theory Turbulent signals: primary variables correlate differently from derivatives b a b a structure of a large eddy primary variables (ui, ) derivatives Because derivatives are not correlated at scales >>, then (*)

22. Nonlocal statistics: theory Consequently, in the inertial range: (Reference: Lipo Wang, PRE 79, 046325 (2009)) In the viscous range: Negative velocity difference is from

23. Nonlocal statistics: numerical results Numerical verification: the Reynolds number effect in the viscous range case 1 case 2

24. Nonlocal statistics: theory Differently from the Kolmogorov 1/3 scaling in the Cartesian, the two-point velocity difference yields a linear scaling when conditioned along gradient trajectories. Physically: extensive strain elongates gradient trajectories, therefore at large separation arclengths gradient trajectories selectively proceed through extensive strain rate regions; while compressive regions are likely to be occupied by shorter gradient trajectories. In the Cartesian, information can be mixed and equally partitioned into each axis to lead to a smaller 1/3 scaling.

25. Nonlocal statistics: application Application: at small scales dissipation elements are compressed, while stretched at large scales. (JFM 608, 113-138 (2008))

26. Nonlocal statistics: extension to other scalars Numerical results of other scalars: the linear scaling holds as well case 1 case 2

27. Nonlocal statistics: extension to other scalars The gradient correlation equation of other scalar : At scales >>:

28. Nonlocal statistics: extension to other scalars DNS results: the effect from Re case 1 case 2

29. Conclusions • Conclusions: • Because of the relative small magnitude of the source terms, there are some generic properties with respect to the strain rate conditioned on gradients of different scalars. Similar to the passive scalar, the mean conditional strain rate of other scalars are negative as well. • The alignment relation of scalar gradients behaves qualitatively different in different regimes; based on which it can be expected that gradient trajectories tend to be parallel along compressive axes, while in extensive planes under the stretching action, these trajectories becomes uncorrelated.

30. Conclusions • Conclusions (continued): • For the passive scalar, along gradient trajectories the two-point velocity difference is proportional to the separation arclength of the gradient trajectories. This property holds for other scalar in a similar way if the Reynolds number is high enough. • This linear scaling is qualitatively different from the classic Kolmogorov’s 1/3 scaling. Physically it can be explained by the selectivity of large gradient trajectories under the stretching action of turbulent velocity.

31. Velocity difference structure function The effect of the weighting factor on statistics Differently from the Cartesian system, in which reference points are equally weighted, the weighting factor of sample points varies along gradient trajectories. sample points (a) Cartesian points (b) trajectory points