No Arabic abstract
Previous work (S. Davidovits and N. J. Fisch, Sudden viscous dissipation of compressing turbulence, Phys. Rev. Lett., 116(105004), 2016) demonstrated that the compression of a turbulent field can lead to a sudden viscous dissipation of turbulent kinetic energy (TKE), and suggested this mechanism could potentially be used to design new fast-ignition schemes for inertial confinement fusion. We expand on previous work by accounting for finite Mach numbers, rather than relying on a zero-Mach-limit assumption as previously done. The finite-Mach-number formulation is necessary to capture a self-consistent feedback mechanism in which dissipated TKE increases the temperature of the system, which in turn modifies the viscosity and thus the TKE dissipation itself. Direct numerical simulations with a tenth-order accurate Pade scheme were carried out to analyze this self-consistent feedback loop for compressing turbulence. Results show that, for finite Mach numbers, the sudden viscous dissipation of TKE still occurs, both for the solenoidal and dilatational turbulent fields. As the domain is compressed, oscillations in dilatational TKE are encountered due to the highly-oscillatory nature of the pressure dilatation. An analysis of the source terms for the internal energy shows that the mechanical work term dominates the viscous turbulent dissipation. As a result, the effect of the suddenly dissipated TKE on temperature is minimal for the Mach numbers tested. Moreover, an analytical expression is derived that confirms the dissipated TKE does not significantly alter the temperature evolution for low Mach numbers, regardless of compression speed. The self-consistent feedback mechanism is thus quite weak for subsonic turbulence, which could limit its applicability for inertial confinement fusion.
We present a simple model for the turbulent kinetic energy behavior of subsonic plasma turbulence undergoing isotropic three-dimensional compression, such as may exist in various inertial confinement fusion experiments or astrophysical settings. The plasma viscosity depends on both the temperature and the ionization state, for which many possible scalings with compression are possible. For example, in an adiabatic compression the temperature scales as $1/L^2$, with $L$ the linear compression ratio, but if thermal energy loss mechanisms are accounted for, the temperature scaling may be weaker. As such, the viscosity has a wide range of net dependencies on the compression. The model presented here, with no parameter changes, agrees well with numerical simulations for a range of these dependencies. This model permits the prediction of the partition of injected energy between thermal and turbulent energy in a compressing plasma.
Intense fluctuations of energy dissipation rate in turbulent flows result from the self-amplification of strain rate via a quadratic nonlinearity, with contributions from vorticity (via the vortex stretching mechanism) and the pressure Hessian tensor, which we analyze here using direct numerical simulations of isotropic turbulence in periodic domains of up to $12288^3$ grid points, and Taylor-scale Reynolds numbers in the range $140-1300$. We extract the statistics of various terms involved in amplification of strain and additionally condition them on the magnitude of strain. We find that strain is overall self-amplified by the quadratic nonlinearity, and depleted via vortex stretching; whereas pressure Hessian acts to redistribute strain fluctuations towards the mean-field and thus depleting intense strain. Analyzing the intense fluctuations of strain in terms of its eigenvalues reveals that the net amplification is solely produced by the third eigenvalue, resulting in strong compressive action. In contrast, the self-amplification terms acts to deplete the other two eigenvalues, whereas vortex stretching acts to amplify them, both effects canceling each other almost perfectly. The effect of the pressure Hessian for each eigenvalue is qualitatively similar to that of vortex stretching, but significantly weaker in magnitude. Our results conform with the familiar notion that intense strain is organized in sheet-like structures, which are in the vicinity of, but never overlap with regions of intense vorticity due to fundamental differences in their amplifying mechanisms.
Recent experiments and simulations have shown that unsteady turbulent flows, before reaching a dynamic equilibrium state, display a universal behaviour. We show that the observed universal non-equilibrium scaling can be explained using a non-equilibrium correction of Kolmogorovs energy spectrum. Given the universality of the experimental and numerical observations, the ideas presented here lay the foundation for the modeling of a wide class of unsteady turbulent flows.
We present a new mesoscale model for ionic liquids based on a low Mach number fluctuating hydrodynamics formulation for multicomponent charged species. The low Mach number approach eliminates sound waves from the fully compressible equations leading to a computationally efficient incompressible formulation. The model uses a Gibbs free energy functional that includes enthalpy of mixing, interfacial energy, and electrostatic contributions. These lead to a new fourth-order term in the mass equations and a reversible stress in the momentum equations. We calibrate our model using parameters for [DMPI+][F6P-], an extensively-studied room temperature ionic liquid (RTIL), and numerically demonstrate the formation of mesoscopic structuring at equilibrium in two and three dimensions. In simulations with electrode boundaries the measured double layer capacitance decreases with voltage, in agreement with theoretical predictions and experimental measurements for RTILs. Finally, we present a shear electroosmosis example to demonstrate that the methodology can be used to model electrokinetic flows.
A systematic study of the influence of the viscous effect on both the spectra and the nonlinear fluxes of conserved as well as non conserved quantities in Navier-Stokes turbulence is proposed. This analysis is used to estimate the helicity dissipation scale which is shown to coincide with the energy dissipation scale. However, it is shown using the decomposition of helicity into eigen modes of the curl operator, that viscous effects have to be taken into account for wave vector smaller than the Kolomogorov wave number in the evolution of these eigen components of the helicity.