Strain, Rotation and Curvature of Non-material Propagating Iso-scalar Surfaces in Homogeneous Turbulence

Abstract

This research aims at gaining some physical insight into the problem of scalar mixing, following the time evolution of propagating iso-surfaces, Y (x, t) = constant, where Y (x, t) stands for any scalar field (e.g., species mass fraction or temperature). First, a rigorous kinematic analysis of non-material line, surface and volume elements, related to propagating iso-scalar surfaces, is presented; this formalism is valid for both constant and variable density flows. Time rates of change of the normal distance and volume between two adjacent iso-surfaces and of area elements, rotation rates of lines and surface elements and an evolution equation for the local mean curvature are obtained. Line and area stretch rates, which encompass additive contributions from the flow and the displacement speed (due to diffusion and reaction), are identified as total strain rates, normal and tangential to the iso-surfaces. Volumetric dilatation rates, addition of line plus area stretch rates, include the mass entrainment rate per unit mass into the non-material volume. Flow and added vorticities, the latter due to gradients of the displacement speed, yield the total vorticity, which provides the real angular velocity of lines and surface elements. A 5123 DNS database for the mixing of inert and reactive scalars in a box of forced statistically stationary and homogeneous turbulence of a constant-density fluid is then examined. A strongly segregated scalar field is prescribed as initial condition. A one-step reaction rate with a characteristic chemical time one order of magnitude greater than the Kolmogorov time micro-scale is used. Data are analyzed at 1.051 large-eddy turnover times after initialization of velocity and scalar fields. Mean negative normal (contractive) and positive tangential (stretching) flow strain rates occur over all mass fractions and scalar-gradient magnitudes. However, means of the total normal strain rate, conditional upon mass fraction, scalar-gradient and mean curvature, are positive everywhere and tend to destroy scalar-gradients for small times. Negative conditioned mean total tangential strain rates (area stretch factor) contract local areas, except for large values of scalar-gradients. Conditional averages of total and added enstrophies are almost identical, which implies a negligible contribution of the flow vorticity to the observed rotation of non-material line and surface elements. The added vorticity is exactly tangential to the iso-surfaces, whereas the flow and total ones are predominantly tangential. Flow sources/sinks of the mean curvature transport equation are much smaller than the added contributions; for this particular DNS database, the local mean curvature development is self-induced by spatial changes of the displacement speed.