No Arabic abstract
The advection and mixing of a scalar quantity by fluid flow is an important problem in engineering and natural sciences. If the fluid is turbulent, the statistics of the passive scalar exhibit complex behavior. This paper is concerned with two Lagrangian scalar turbulence models based on the recent fluid deformation model that can be shown to reproduce the statistics of passive scalar turbulence for a range of Reynolds numbers. For these models, we demonstrate how events of extreme passive scalar gradients can be recovered by computing the instanton, i.e., the saddle-point configuration of the associated stochastic field theory. It allows us to both reproduce the heavy-tailed statistics associated with passive scalar turbulence, and recover the most likely mechanism leading to such extreme events. We further demonstrate that events of large negative strain in these models undergo spontaneous symmetry breaking.
Turbulent fluid flows are ubiquitous in nature and technology, and are mathematically described by the incompressible Navier-Stokes equations (INSE). A hallmark of turbulence is spontaneous generation of intense whirls, resulting from amplification of the fluid rotation-rate (vorticity) by its deformation-rate (strain). This interaction, encoded in the non-linearity of INSE, is non-local, i.e., depends on the entire state of the flow, constituting a serious hindrance in turbulence theory and in establishing regularity of INSE. Here, we unveil a novel aspect of this interaction, by separating strain into local and non-local contributions utilizing the Biot-Savart integral of vorticity in a sphere of radius R. Analyzing highly-resolved numerical turbulent solutions to INSE, we find that when vorticity becomes very large, the local strain over small R surprisingly counteracts further amplification. This uncovered self-attenuation mechanism is further shown to be connected to local Beltramization of the flow, and could provide a direction in establishing the regularity of INSE.
We study transport of a weakly diffusive pollutant (a passive scalar) by thermoconvective flow in a fluid-saturated horizontal porous layer heated from below under frozen parametric disorder. In the presence of disorder (random frozen inhomogeneities of the heating or of macroscopic properties of the porous matrix), spatially localized flow patterns appear below the convective instability threshold of the system without disorder. Thermoconvective flows crucially effect the transport of a pollutant along the layer, especially when its molecular diffusion is weak. The effective (or eddy) diffusivity also allows to observe the transition from a set of localized currents to an almost everywhere intense global flow. We present results of numerical calculation of the effective diffusivity and discuss them in the context of localization of fluid currents and the transition to a global flow. Our numerical findings are in a good agreement with the analytical theory we develop for the limit of a small molecular diffusivity and sparse domains of localized currents. Though the results are obtained for a specific physical system, they are relevant for a broad variety of fluid dynamical systems.
The reduction of dimensionality of physical systems, specially in fluid dynamics, leads in many situations to nonlinear ordinary differential equations which have global invariant manifolds with algebraic expressions containing relevant physical information of the original system. We present a method to identify such manifolds, and we apply it to a reduced model for the Lagrangian evolution of field gradients in homogeneous and isotropic turbulence with a passive scalar.
In scalar turbulence it is sometimes the case that the scalar diffusivity is smaller than the viscous diffusivity. The thermally-driven turbulent convection in water is a typical example. In such a case the smallest scale in the problem is the Batchelor scale $l_b$, rather than the Kolmogorov scale $l_k$, as $l_b = l_k/Sc^{1/2}$, where Sc is the Schmidt number (or Prandtl number in the case of temperature). In the numerical studies of such scalar turbulence, the conventional approach is to use a single grid for both the velocity and scalar fields. Such single-resolution scheme often over-resolves the velocity field because the resolution requirement for scalar is higher than that for the velocity field, since $l_b<l_k$ for $Sc>1$. In this paper we put forward an algorithm that implements the so-called multiple-resolution method with a finite-volume code. In this scheme, the velocity and pressure fields are solved in a coarse grid, while the scalar field is solved in a dense grid. The central idea is to implement the interpolation scheme on the framework of finite-volume to reconstruct the divergence-free velocity from the coarse to the dense grid. We demonstrate our method using a canonical model system of fluid turbulence, the Rayleigh-Benard convection. We show that, with the tailored mesh design, considerable speed-up for simulating scalar turbulence can be achieved, especially for large Schmidt (Prandtl) numbers. In the same time, sufficient accuracy of the scalar and velocity fields can be achived by this multiple-resolution scheme. Although our algorithm is demonstrated with a case of an active scalar, it can be readily applied to passive scalar turbulent flows.
Lagrangian properties obtained from a Particle Tracking Velocimetry experiment in a turbulent flow at intermediate Reynolds number are presented. Accurate sampling of particle trajectories is essential in order to obtain the Lagrangian structure functions and to measure intermittency at small temporal scales. The finiteness of the measurement volume can bias the results significantly. We present a robust way to overcome this obstacle. Despite no fully developed inertial range we observe strong intermittency at the scale of dissipation. The multifractal model is only partially able to reproduce the results.