Inertial Range Of Turbulence Research Articles

MHD Inertial and Energy-containing Range Turbulence Anisotropy in the Young Solar Wind

We study solar wind turbulence anisotropy in the inertial and energy-containing ranges in the inbound and outbound directions during encounters 1–9 by the Parker Solar Probe (PSP) for distances between ∼21 and 65 R ⊙. Using the Adhikari et al. approach, we derive theoretical equations to calculate the ratio between the 2D and slab fluctuating magnetic energy, fluctuating kinetic energy, and the outward/inward Elsässer energy in the inertial range. For this, in the energy-containing range, we assume a wavenumber k −1 power law. In the inertial range, for the magnetic field fluctuations and the outward/inward Elsässer energy, we consider that (i) both 2D and slab fluctuations follow a power law of k −5/3, and (ii) the 2D and slab fluctuations follow the power laws with k −5/3 and k −3/2, respectively. For the velocity fluctuations, we assume that both the 2D and slab components follow a k −3/2 power law. We compare the theoretical results of the variance anisotropy in the inertial range with the derived observational values measured by PSP, and find that the energy density of 2D fluctuations is larger than that of the slab fluctuations. The theoretical variance anisotropy in the inertial range relating to the k −5/3 and k −3/2 power laws between 2D and slab turbulence exhibits a smaller value in comparison to assuming the same power law k −5/3 between 2D and slab turbulence. Finally, the observed turbulence energy measured by PSP in the energy-containing range is found to be similar to the theoretical result of a nearly incompressible/slab turbulence description.

Turbulence Transitions in Kelvin–Helmholtz Instability “Tube” and “Knot” Dynamics: Vorticity, Helicity, and Twist Waves

We address the sources and dynamics of vorticity and helicity and their relations in transitions to turbulence arising due to Kelvin–Helmholtz instability (KHI) “Tube” and “Knot” (T&K) events. Such events are common in the atmosphere and oceans, and initial numerical simulations reveal that T&K dynamics significantly accelerate turbulence transitions and enhance KHI peak and mean energy dissipation rates. KHI T&K events arise where emerging KH billows exhibit varying wavelengths, phases, amplitudes, and/or discontinuities along their axes. As the KH billows intensify, these regions evolve roughly orthogonal billow cores and induced vortex tubes in close proximity. Their mutual advection as they intensify induces large-amplitude Kelvin vortex waves, or “twist waves”, that arise where locally uniform vortices are distorted by axial or radial advection. The twist waves propagate along, and fragment, the vortex tubes and billow cores, thus accounting for the emergence of helicity and the down-scale energy, enstrophy, and helicity fluxes within the turbulence inertial range. We describe the results of four direct numerical simulations (DNS) addressing KHI T&K dynamics in large and idealized small domains. The large-domain vorticity fields reveal the character and diversity of KHI T&K dynamics, the emergence of twist waves at larger and smaller scales, and their driving of turbulence transitions. Two small-domain DNS exhibit idealized KHI T&K events arising from KH billows that are mis-aligned and that exhibit phase variability along their axes. A third examines the interactions of two vortex tubes in close proximity. These reveal that twist waves drive the character and evolutions of the vorticity and helicity fields.

Stochastic modeling of turbulent mixing based on a hierarchical swapping of fluid parcels

AbstractTurbulent mixing is an omnipresent phenomenon that constantly affects our everyday life and plays an important role in a variety of industrial applications. The simulation of turbulent mixing poses great challenges, since the full resolution of all relevant length and time scales is associated with an immense computational effort. This limitation can be overcome by only resolving the large‐scale effects and completely model the sub‐grid scales. The development of an accurate sub‐grid mixing model is therefore a key challenge to capture all interactions in the sub‐grid scales. At this place, the hierarchical parcel‐swapping (HiPS) model formulated by A.R. Kerstein [J. Stat. Phys. 153, 142–161 (2013)] represents a computationally efficient and scale‐resolving turbulent mixing model. HiPS mimics the effects of turbulence on time‐evolving, diffusive scalar fields. In HiPS, the diffusive scalar fields or a state space is interpreted as a binary tree structure, which is an alternative approach compared to the most common mixing models. Every level of the tree represents a specific length and time scale, which is based on turbulence inertial range scaling. The state variables are only located at the base of the tree and are treated as fluid parcels. The effects of turbulent advection are represented by stochastic swaps of sub‐trees at rates determined by turbulent time scales associated with the sub‐trees. The mixing only takes places between adjacent fluid parcels and at rates consistent with the prevailing diffusion time scales. In this work, the HiPS model formulation for the simulation of passive scalar mixing is detailed first. Preliminary results for the mean square displacement, passive scalar probability density function (PDF) and scalar dissipation rate are given and reveal the strengths of the HiPS model considering the reduced order and computational efficiency. These model investigations are an important step of further HiPS advancements. The integrated auxiliary binary tree structure allows HiPS to satisfy a large number of criteria for a good mixing model. From this point of view, HiPS is an attractive candidate for modeling the mixing in transported PDF methods.

Transmission of magnetic island modes across interplanetary shocks: comparison of theory and observations

Interplanetary shock waves are observed frequently in turbulent solar wind. They naturally enhance the temperature/entropy of the plasma through which they propagate. Moreover, many studies have shown that they also act as an amplifier of the fluctuations incident on the shock front. Solar wind turbulent fluctuations can be well described as the superposition of quasi-2D and slab components, the former being energetically dominant. In this paper, we address the interaction of fast forward shocks observed by the Wind spacecraft at 1 AU and quasi-2D turbulent fluctuations in the framework of the Zank et al. (2021) transmission model and we compare model predictions with observations. Our statistical study includes 378 shocks with varying upstream conditions and Mach numbers. We estimate the average ratio of the downstream observed and theoretically predicted power spectra within the inertial range of turbulence. We find that the distributions of this ratio for the whole set and for the subset of shocks that met the assumptions of the model, are remarkably close. We argue that a large statistical spread of the distributions of this ratio is governed by the inherent variation of the upstream conditions. Our findings suggest that the model predicts the downstream fluctuations with a good accuracy and that it may be adopted for a wider class of shocks than it was originally meant for.

Investigation of turbulent mixing using a stochastic hierarchical parcel swapping mixing model

AbstractTurbulent mixing plays an important role in a variety of applications ranging from astrophysics to combustion and even pollutant dispersion. The Direct Numerical Simulation (DNS) that resolves all scales is not feasible for most engineering applications since the flow has a wide range of length and time scales, which yields extremely high resolution requirements. Large Eddy Simulations (LES) overcome this limitation by modeling the sub‐grid scale effects. In transported Probability Density Function (PDF) methods, the key challenge is to develop an accurate mixing model. At this point, the Hierarchical Parcel‐Swapping (HiPS) model, introduced by A.R. Kerstein [J. Stat. Phys. 153, 142‐161 (2013)], is an attractive candidate. It is characterized by a computationally efficient representation of the effects of turbulence on a time‐evolving structure of diffusive scalar fields. The interpretation of the diffusive scalar fields or a state space as a binary tree structure is an alternative approach compared to the most common mixing models. The characteristic feature of HiPS is that every level of the tree corresponds to a specific length and time scale, which is based on turbulence inertial range scaling. The state variables only reside at the base of the tree and are understood as fluid parcels. The effects of turbulent advection are represented by stochastic swaps of sub‐trees at rates determined by prevailing turbulent time scales associated with the sub‐trees. The mixing of adjacent fluid parcels is done either instantaneously or at rates consistent with the corresponding diffusion time scales. In this work, HiPS is detailed for the simulation of passive scalar mixing first. Preliminary results for the scalar power spectra, mean square displacement and scalar dissipation rate are shown and reveal a reasonable agreement with experimental findings. Furthermore, the integrated binary tree structure allows to satisfy a large number of criteria for a good mixing model. Considering the reduced order and computational efficiency, HiPS is an attractive candidate for modeling the mixing in transported PDF methods.

Multi‐Scale Kelvin‐Helmholtz Instability Dynamics Observed by PMC Turbo on 12 July 2018: 2. DNS Modeling of KHI Dynamics and PMC Responses

AbstractKjellstrand et al. (2022), https://10.1029/2021JD036232 describes the evolution and dynamics of a strong, large‐scale Kelvin‐Helmholtz instability (KHI) event observed in polar mesospheric clouds (PMCs) on 12 July 2018 by high‐resolution imagers aboard the PMC Turbulence (PMC Turbo) stratospheric long‐duration balloon experiment. The imaging provides evidence of KH billow interactions and instabilities that are strongly influenced by gravity waves at larger scales. Specific features include initially separated regions of KHI, secondary convective and KH instabilities of individual billows, and “tubes” and “knots” that arise where billow cores are mis‐aligned or discontinuous along their axes. This study describes a direct numerical simulation of KH billow interactions in a periodic domain seeded with random initial noise that enables excitation of multiple KH billows exhibiting variable phase structures that capture multiple features of the observed KHI dynamics. Variable KH billow phases along their axes yield initial vortex tubes having diagonal alignments that link adjacent, but mis‐aligned, billow cores. Weak initial vortex tubes and billow cores having nearly orthogonal alignments amplify, interact strongly, and drive intense vortex knots at these sites. These vortex tube and knot (T&K) dynamics excite “twist waves” that unravel the initial vortex tubes, and drive increasingly strong vortex interactions and a cascade of energy and enstrophy to successively smaller scales in the turbulence inertial range. The implications of T&K dynamics are much more rapid and intense breakdown and decay of the KH billows, and significantly enhanced energy dissipation rates, where these interactions occur.

High-latitude Observations of Inertial-range Turbulence by the Ulysses Spacecraft During the Solar Minimum of 1993–96

Abstract We have examined Ulysses magnetic field measurements for the years 1993 through 1996 as the spacecraft moved sunward from 5 au at high southern latitudes, passing through perihelion during the first fast-latitude scan to achieve high northern latitudes, and finally returning to 5 au. These years represent near-solar-minimum activity, providing a clear measure of high-latitude solar-wind turbulence. We apply a series of tests to the data, examining both the magnetic variance anisotropy and the underlying wavevector anisotropy, finding them to be consistent with past 1 au observations. The variance anisotropy depends upon both the thermal proton temperature parameter and the amplitude of the magnetic power spectrum, while the underlying wavevector anisotropy is dominated by the component perpendicular to the mean magnetic field. We also examine the amplitude of the magnetic power spectrum as well as the associated turbulent transport of energy to small scales that results in the heating of the thermal plasma. The measured turbulence is found to be stronger than that seen at low latitudes by the Voyager spacecraft as it traverses the distance from 1 to 5 au during the years approaching solar maximum. If the high- and low-latitude sources are comparable, this would indicate that while the heating processes are active in both regions, the turbulence has had less decay time in the transport of energy to small scales. Alternatively, it may also be that the high-latitude source is stronger.

Scaling Laws for Partially Developed Turbulence

We formulate multifractal models for velocity differences and gradients which describe the full range of length scales in turbulent flow, namely: laminar, dissipation, inertial, and stirring ranges. The models subsume existing models of inertial range turbulence. In the localized ranges of length scales in which the turbulence is only partially developed, we propose multifractal scaling laws with scaling exponents modified from their inertial range values. In local regions, even within a fully developed turbulent flow, the turbulence is not isotropic nor scale invariant due to the influence of larger turbulent structures (or their absence). For this reason, turbulence that is not fully developed is an important issue which inertial range study can not address. In the ranges of partially developed turbulence, the flow can be far from universal, so that standard inertial range turbulence scaling models become inapplicable. The model proposed here serves as a replacement. Details of the fitting of the parameters for the τp and ζp models in the dissipation range are discussed. Some of the behavior of ζp for larger p is unexplained. The theories are verified by comparing to high resolution simulation data.

Erratum: “Dependence of 3D Self-correlation Level Contours on the Scales in the Inertial Range of Solar Wind Turbulence” (2019, ApJL, 883, L9)

Erratum: “Dependence of 3D Self-correlation Level Contours on the Scales in the Inertial Range of Solar Wind Turbulence” (2019, ApJL, 883, L9)

Evolution of localized magnetic field perturbations and the nature of turbulent dynamo

Kinematic dynamo in incompressible isotropic turbulent flows with high magnetic Prandtl number is considered. The approach interpreting an arbitrary magnetic field distribution as a superposition of localized perturbations (blobs) is developed. We derive a general relation between stochastic properties of an isolated blob and a stochastically homogenous distribution of magnetic field advected by the same stochastic flow. This relation allows us to investigate the evolution of a localized blob at a late stage when its size exceeds the viscous scale. It is shown that in three-dimensional flows, the average magnetic field of the blob increases exponentially in the inertial range of turbulence, as opposed to the late-batchelor stage when it decreases. Our approach reveals the mechanism of dynamo generation in the inertial range both for blobs and homogenous contributions. It explains the absence of dynamo in the two-dimensional case and its efficiency in three dimensions. We propose a way to observe the mechanism in numerical simulations.

Advective structure functions in anisotropic two-dimensional turbulence

Abstract

Laboratory model for plastic fragmentation in the turbulent ocean

We study the fragmentation of deformable and brittle fibers in the inertial range of turbulence using laboratory experiments and numerical simulations. The fragmentation process is shown to be limited at small scales by a physical cut-off length due to fluid-structure interactions of the object with turbulence, and thus independent of the fiber brittleness. This scenario, comprehensively modeled by an evolution equation, leads to the accumulation of fragments slightly longer than the cut-off scale, as smaller fragments are too short to be deformed and broken by the turbulence. This result may improve our understanding of microplastic formation in the ocean.

Nanoflare Theory Revisited

Abstract At any scale l in the turbulent inertial range, the magnetic field can be divided up into a large-scale component and a small-scale, high spatial frequency component which undergoes magnetic reversals. Such local reconnections, i.e., on any inertial scale l, seem to be an inseparable part of magnetohydrodynamic (MHD) turbulence, whose collective outcome can lead to global reconnection with a rate independent of the small-scale physics dominant at dissipative scales. We show that this picture, known as stochastic reconnection, is intimately related to nanoflare theory, proposed long ago to explain solar coronal heating. We argue that, due to stochastic flux freezing, a generalized version of magnetic flux freezing in turbulence, the field follows the flow in a statistical sense. Turbulence bends and stretches the field, increasing its spatial complexity. Strong magnetic shears associated with such a highly tangled field can trigger local reversals and field annihilations on a wide range of inertial scales which convert magnetic energy into kinetic and thermal energy. The former may efficiently enhance turbulence and the latter heat generation. We support this theoretical picture using scaling laws of MHD turbulence and also recent analytical and numerical studies which suggest a statistical correlation between magnetic spatial complexity and energy dissipation. Finally, using an MHD numerical simulation, we show that the time evolution of the magnetic complexity is statistically correlated with the rate of kinetic energy injection and/or magnetic-to-thermal energy conversion, in agreement with our proposed theoretical picture.

Power and complexity in stochastic reconnection

Previous work has invoked kinetic and magnetic spatial complexities, associated with velocity and magnetic fields u(x,t) and B(x,t), respectively, in order to study magnetic reconnection and diffusion in turbulent and magnetized fluids. In this paper, using the coarse-grained momentum equation, we argue that the fluid jets associated with magnetic reconnection events at an arbitrary scale l in the turbulence inertial range are predominantly driven by the Lorentz force Nl=(j×B)l−jl×Bl. This force is induced by the subscale currents and is analogous to the turbulent electromotive force El=(u×B)l−ul×Bl in dynamo theories. Typically, high (low) magnetic complexities during reconnection imply large (small) spatial gradients for the magnetic field, i.e., strong (weak) Lorentz forces Nl. Reconnection launches jets of fluid, hence the rate of change of kinetic complexity is expected to strongly correlate with the power injected by the Lorentz force Nl. We test this prediction using an incompressible, homogeneous magnetohydrodynamic (MHD) simulation and associate it with previous results. It follows that the stronger (weaker) the turbulence, the more (less) complex the magnetic field and the stronger (weaker) the driving Lorentz forces and thus the ensuing reconnection.

Switchbacks in the Near-Sun Magnetic Field: Long Memory and Impact on the Turbulence Cascade

Abstract One of the most striking observations made by Parker Solar Probe during its first solar encounter is the omnipresence of rapid polarity reversals in a magnetic field that is otherwise mostly radial. These so-called switchbacks strongly affect the dynamics of the magnetic field. We concentrate here on their macroscopic properties. First, we find that these structures are self-similar, and have neither a characteristic magnitude, nor a characteristic duration. Their waiting time statistics show evidence of aggregation. The associated long memory resides in their occurrence rate, and is not inherent to the background fluctuations. Interestingly, the spectral properties of inertial range turbulence differ inside and outside of switchback structures; in the latter the 1/f range extends to higher frequencies. These results suggest that outside of these structures we are in the presence of lower-amplitude fluctuations with a shorter turbulent inertial range. We conjecture that these correspond to a pristine solar wind.

Hierarchical Parcel‐Swapping: An efficient mixing model for turbulent reactive flows

AbstractHierarchical Parcel‐Swapping (HiPS) developed by A.R. Kerstein [J. Stat. Phys. 153, 142‐161 (2013)] is a computationally efficient and novel model for the effects of turbulence on time‐evolving, diffusive scalar fields. The characteristic feature of HiPS is the interpretation of the one‐dimensional flow domain or a state space as a binary tree structure. Every tree level corresponds to a specific length and time scale, which is based on a turbulence inertial range scaling. The state variables reside at the base of the tree and are interpreted as fluid parcels. The effects of turbulent advection are represented by stochastic swaps of sub‐trees at rates determined by turbulent time scales associated with the sub‐trees. The mixing of adjacent fluid parcels is done at rates consistent with the prevailing diffusion time scales. In this work, we investigate the influence of turbulent time scale variations on an isothermal series‐parallel reaction scheme. The production of a desired chemical species is evaluated by means of a defined selectivity and is strongly affected by the underlying mixing time scales.

Scaling laws in Hall inertial-range turbulence

Abstract. There is an increasing amount of observational evidence in space plasmas for the breakdown of inertial-range spectra of magnetohydrodynamic (MHD) turbulence on spatial scales smaller than the ion-inertial length. Magnetic energy spectra often exhibit a steepening, which is reminiscent of dissipation of turbulence energy, for example in wave–particle interactions. Electric energy spectra, on the other hand, tend to be flatter than those of MHD turbulence, which is indicative of a dispersive process converting magnetic into electric energy in electromagnetic wave excitation. Here we develop a model of the scaling laws and the power spectra for the Hall inertial range in plasma turbulence. In the present paper we consider a two-dimensional geometry with no wave vector component parallel to the magnetic field as is appropriate in Hall MHD. A phenomenological approach is taken. The Hall electric field attains an electrostatic component when the wave vectors are perpendicular to the mean magnetic field. The power spectra of Hall turbulence are steep for the magnetic field with a slope of -7/3 for compressible magnetic turbulence; they are flatter for the Hall electric field with a slope of -1/3. Our model for the Hall turbulence gives a possible explanation for the steepening of the magnetic energy spectra in the solar wind as an indication of neither the dissipation range nor the dispersive range but as the Hall inertial range. Our model also reproduces the shape of energy spectra in Kelvin–Helmholtz turbulence observed at the Earth's magnetopause.

Dependence of 3D Self-correlation Level Contours on the Scales in the Inertial Range of Solar Wind Turbulence

Abstract The self-correlation level contours at the 1010 cm scale reveal a 3D isotropic feature in the slow solar wind and a quasi-anisotropic feature in the fast solar wind. However, the 1010 cm scale is approximately near the low-frequency break (outer scale of turbulence cascade), especially in the fast wind. How the self-correlation level contours behave with dependence on the scales in the inertial range of solar wind turbulence remains unknown. Here we present the 3D self-correlation function level contours and their dependence on the scales in the inertial range for the first time. We use data at 1 au from instruments on the Wind spacecraft in the period 2005–2018. We show the 3D isotropic self-correlation level contours of the magnetic field in the inertial range of both slow and fast solar wind turbulence. We also find that the self-correlation level contours of the velocity in the inertial range present 2D anisotropy with an elongation in the perpendicular direction and 2D isotropy in the plane perpendicular to the mean magnetic field. These results indicate differences between the magnetic field and the velocity, providing new clues to interpret the solar wind turbulence on the inertial scale.

A population balance model for large eddy simulation of polydisperse droplet evolution

In the context of many applications of turbulent multi-phase flows, knowledge of the dispersed phase size distribution and its evolution is critical to predicting important macroscopic features. We develop a large eddy simulation (LES) model that can predict the turbulent transport and evolution of size distributions, for a specific subset of applications in which the dispersed phase can be assumed to consist of spherical droplets, and occurring at low volume fraction. We use a population dynamics model for polydisperse droplet distributions specifically adapted to a LES framework including a model for droplet breakup due to turbulence, neglecting coalescence consistent with the assumed small dispersed phase volume fractions. We model the number density fields using an Eulerian approach for each bin of the discretized droplet size distribution. Following earlier methods used in the Reynolds-averaged Navier–Stokes framework, the droplet breakup due to turbulent fluctuations is modelled by treating droplet–eddy collisions as in kinetic theory of gases. Existing models assume the scale of droplet–eddy collision to be in the inertial range of turbulence. In order to also model smaller droplets comparable to or smaller than the Kolmogorov scale we extend the breakup kernels using a structure function model that smoothly transitions from the inertial to the viscous range. The model includes a dimensionless coefficient that is fitted by comparing predictions in a one-dimensional version of the model with a laboratory experiment of oil droplet breakup below breaking waves. After initial comparisons of the one-dimensional model to measurements of oil droplets in an axisymmetric jet, it is then applied in a three-dimensional LES of a jet in cross-flow with large oil droplets of a single size being released at the source of the jet. We model the concentration fields using $N_{d}=15$ bins of discrete droplet sizes and solve scalar transport equations for each bin. The resulting droplet size distributions are compared with published experimental data, and good agreement for the relative size distribution is obtained. The LES results also enable us to quantify size distribution variability. We find that the probability distribution functions of key quantities such as the total surface area and the Sauter mean diameter of oil droplets are highly variable, some displaying strong non-Gaussian intermittent behaviour.

Dust Polarization Maps from TIGRESS: E/B Power Asymmetry and TE Correlation

Abstract We present the first large set of all-sky synthetic dust polarization maps derived directly from a self-consistent magnetohydrodynamics simulation using the TIGRESS framework. Turbulence in this simulation is predominantly driven by supernova explosions, with rates that are self-consistently regulated by feedback loops. The simulation covers both the outer scale and inertial range of turbulence with uniformly high resolution. The shearing-box utilized in the simulation, in concert with resolved supernova-driven turbulence, enables the capturing of generation, growth, and saturation of both turbulent and mean magnetic fields. We construct polarization maps at 353 GHz, as seen by observers inside a model of the multiphase, turbulent, magnetized interstellar medium (ISM). To fully sample the simulated ISM state, we use 350 snapshots spanning over (more than six feedback loops) and nine representative observers. The synthetic skies show a prevalent E/B power asymmetry ( ) and positive TE correlation in broad agreement with observations by the Planck satellite. However, the ranges of and are generally lower than those measured by Planck. We find large fluctuations of E/B asymmetry and TE correlation depending on the observer’s position and temporal fluctuations of ISM properties due to bursts of star formation. The synthetic maps are made publicly available to provide novel models of the microwave sky.