Magnetic Field Amplification in Galaxy Clusters and its Simulation

We review the present theoretical and numerical understanding of magnetic field amplification in cosmic large-scale structure, on length scales of galaxy clusters and beyond. Structure formation drives compression and turbulence, which amplify tiny magnetic seed fields to the microGauss values that are observed in the intracluster medium. This process is intimately connected to the properties of turbulence and the microphysics of the intra-cluster medium. Additional roles are played by merger induced shocks that sweep through the intra-cluster medium and motions induced by sloshing cool cores. The accurate simulation of magnetic field amplification in clusters still poses a serious challenge for simulations of cosmological structure formation. We review the current literature on cosmological simulations that include magnetic fields and outline theoretical as well as numerical challenges.Comment: 60 pages, 19 Figure

Regularization modeling for large-eddy simulation of homogeneous isotropic decaying turbulence

Inviscid regularization modeling of turbulent flow is investigated. Homogeneous, isotropic, decaying turbulence is simulated at a range of filter widths. A coarse-graining of turbulent flow arises from the direct regularization of the convective nonlinearity in the Navier–Stokes equations. The regularization is translated into its corresponding sub-filter model to close the equations for large-eddy simulation (LES). The accuracy with which primary turbulent flow features are captured by this modeling is investigated for the Leray regularization, the Navier–Stokes-α formulation (NS-α), the simplified Bardina model and a modified Leray approach. On a PDE level, each regularization principle is known to possess a unique, strong solution with known regularity properties. When used as turbulence closure for numerical simulations, significant differences between these models are observed. Through a comparison with direct numerical simulation (DNS) results, a detailed assessment of these regularization principles is made. The regularization models retain much of the small-scale variability in the solution. The smaller resolved scales are dominated by the specific sub-filter model adopted. We find that the Leray model is in general closest to the filtered DNS results, the modified Leray model is found least accurate and the simplified Bardina and NS-α models are in between, as far as accuracy is concerned. This rough ordering is based on the energy decay, the Taylor Reynolds number and the velocity skewness, and on detailed characteristics of the energy dynamics, including spectra of the energy, the energy transfer and the transfer power. At filter widths up to about 10% of the computational domain-size, the Leray and NS-α predictions were found to correlate well with the filtered DNS data. Each of the regularization models underestimates the energy decay rate and overestimates the tail of the energy spectrum. The correspondence with unfiltered DNS spectra was observed often to be closer than with filtered DNS for several of the regularization models

Large-Eddy Simulation and experimental study of cycle-to-cycle variations of stable and unstable operating points in a spark ignition engine

This article presents a comparison between experiments and Large-Eddy Simulation (LES) of a spark ignition engine on two operating points: a stable one characterized by low cycle-to-cycle variations (CCV) and an unstable one with high CCV. In order to match the experimental cycle sample, 75 full cycles (with combustion) are computed by LES. LES results are compared with experiments by means of pressure signals in the intake and exhaust ducts, in-cylinder pressure, chemiluminescence and OH Planar Laser Induced Fluorescence (PLIF). Results show that LES is able to: (1) reproduce the flame behavior in both cases (low and high CCV) in terms of position, shape and timing; (2) distinguish a stable point from an unstable one; (3) predict quantitatively the CCV levels of the two fired operating points. For the unstable case, part of the observed CCV is due to incomplete combustion. The results are then used to analyze the incomplete combustion phenomenon which occurs for some cycles of the unstable point and propose modification of the spark location to control CCV

The Relative Humidity in an Isentropic Advection–Condensation Model: Limited Poleward Influence and Properties of Subtropical Minima

An idealized model of advection and condensation of water vapor is considered as a representation of processes influencing the humidity distribution along isentropic surfaces in the free troposphere. Results are presented for how the mean relative humidity distribution varies in response to changes in the distribution of saturation specific humidity and in the amplitude of a tropical moisture source. Changes in the tropical moisture source are found to have little effect on the relative humidity poleward of the subtropical minima, suggesting a lack of poleward influence despite much greater water vapor concentrations at lower latitudes. The subtropical minima in relative humidity are found to be located just equatorward of the inflection points of the saturation specific humidity profile along the isentropic surface. The degree of mean subsaturation is found to vary with the magnitude of the meridional gradient of saturation specific humidity when other parameters are held fixed. The atmospheric relevance of these results is investigated by comparison with the positions of the relative humidity minima in reanalysis data and by examining poleward influence of relative humidity in simulations with an idealized general circulation model. It is suggested that the limited poleward influence of relative humidity may constrain the propagation of errors in simulated humidity fields

Space-time correlations in turbulent flow: A review

This paper reviews some of the principal uses, over almost seven decades, of correlations, in both Eulerian and Lagrangian frames of reference, of properties of turbulent flows at variable spatial locations and variable time instants. Commonly called space--time correlations, they have been fundamental to theories and models of turbulence as well as for the analyses of experimental and direct numerical simulation turbulence data.Comment: 16 pages, 7 figure

Statistical state dynamics of weak jets in barotropic beta-plane turbulence

Zonal jets in a barotropic setup emerge out of homogeneous turbulence through a flow-forming instability of the homogeneous turbulent state (`zonostrophic instability') which occurs as the turbulence intensity increases. This has been demonstrated using the statistical state dynamics (SSD) framework with a closure at second order. Furthermore, it was shown that for small supercriticality the flow-forming instability follows Ginzburg-Landau (G-L) dynamics. Here, the SSD framework is used to study the equilibration of this flow-forming instability for small supercriticality. First, we compare the predictions of the weakly nonlinear G-L dynamics to the fully nonlinear SSD dynamics closed at second order for a wide ranges of parameters. A new branch of jet equilibria is revealed that is not contiguously connected with the G-L branch. This new branch at weak supercriticalities involves jets with larger amplitude compared to the ones of the G-L branch. Furthermore, this new branch continues even for subcritical values with respect to the linear flow-forming instability. Thus, a new nonlinear flow-forming instability out of homogeneous turbulence is revealed. Second, we investigate how both the linear flow-forming instability and the novel nonlinear flow-forming instability are equilibrated. We identify the physical processes underlying the jet equilibration as well as the types of eddies that contribute in each process. Third, we propose a modification of the diffusion coefficient of the G-L dynamics that is able to capture the asymmetric evolution for weak jets at scales other than the marginal scale (side-band instabilities) for the linear flow-forming instability.Comment: 27 pages, 17 figure

Stirring up trouble: Multi-scale mixing measures for steady scalar sources

The mixing efficiency of a flow advecting a passive scalar sustained by steady sources and sinks is naturally defined in terms of the suppression of bulk scalar variance in the presence of stirring, relative to the variance in the absence of stirring. These variances can be weighted at various spatial scales, leading to a family of multi-scale mixing measures and efficiencies. We derive a priori estimates on these efficiencies from the advection--diffusion partial differential equation, focusing on a broad class of statistically homogeneous and isotropic incompressible flows. The analysis produces bounds on the mixing efficiencies in terms of the Peclet number, a measure the strength of the stirring relative to molecular diffusion. We show by example that the estimates are sharp for particular source, sink and flow combinations. In general the high-Peclet number behavior of the bounds (scaling exponents as well as prefactors) depends on the structure and smoothness properties of, and length scales in, the scalar source and sink distribution. The fundamental model of the stirring of a monochromatic source/sink combination by the random sine flow is investigated in detail via direct numerical simulation and analysis. The large-scale mixing efficiency follows the upper bound scaling (within a logarithm) at high Peclet number but the intermediate and small-scale efficiencies are qualitatively less than optimal. The Peclet number scaling exponents of the efficiencies observed in the simulations are deduced theoretically from the asymptotic solution of an internal layer problem arising in a quasi-static model.Comment: 37 pages, 7 figures. Latex with RevTeX4. Corrigendum to published version added as appendix