Browsing by Author "Vlaykov, Dimitar G."

Compressible magnetohydrodynamic (MHD) turbulence is ubiquitous in astrophysical phenomena ranging from the intergalactic to the stellar scales. In studying them, numerical simulations are nearly inescapable, due to the large degree of nonlinearity involved. However, the dynamical ranges of these phenomena are much larger than what is computationally accessible. In large eddy simulations (LESs), the resulting limited resolution effects are addressed explicitly by introducing to the equations of motion additional terms associated with the unresolved, subgrid-scale dynamics. This renders the system unclosed. We derive a set of nonlinear structural closures for the ideal MHD LES equations with particular emphasis on the effects of compressibility. The closures are based on a gradient expansion of the finite-resolution operator [W. K. Yeo (CUP, 1993)] and require no assumptions about the nature of the flow or magnetic field. Thus, the scope of their applicability ranges from the sub-to the hyper-sonic and -Alfvenic regimes. The closures support spectral energy cascades both up and down-scale, as well as direct transfer between kinetic and magnetic resolved and unresolved energy budgets. They implicitly take into account the local geometry, and in particular, the anisotropy of the flow. Their properties are a priori validated in Paper II [P. Grete et al., Phys. Plasmas 23, 062317 (2016)] against alternative closures available in the literature with respect to a wide range of simulation data of homogeneous and isotropic turbulence. Published by AIP Publishing.

Even though compressible plasma turbulence is encountered in many astrophysical phenomena, its effect is often not well understood. Furthermore, direct numerical simulations are typically not able to reach the extreme parameters of these processes. For this reason, large-eddy simulations (LES), which only simulate large and intermediate scales directly, are employed. The smallest, unresolved scales and the interactions between small and large scales are introduced by means of a subgridscale (SGS) model. We propose and verify a new set of nonlinear SGS closures for future application as an SGS model in LES of compressible magnetohydrodynamics. We use 15 simulations (without explicit SGS model) of forced, isotropic, homogeneous turbulence with varying sonic Mach number M-s = 0.2-20 as reference data for the most extensive a priori tests performed so far in literature. In these tests, we explicitly filter the reference data and compare the performance of the new closures against the most widely tested closures. These include eddy-viscosity and scale-similarity type closures with different normalizations. Performance indicators are correlations with the turbulent energy and cross-helicity flux, the average SGS dissipation, the topological structure and the ability to reproduce the correct magnitude and the direction of the SGS vectors. We find that only the new nonlinear closures exhibit consistently high correlations (median value >0.8) with the data over the entire parameter space and outperform the other closures in all tests. Moreover, we show that these results are independent of resolution and chosen filter scale. Additionally, the new closures are effectively coefficient-free with a deviation of less than 20%. Published by AIP Publishing.

Large-eddy simulations (LES) are a powerful tool in understanding processes that are inaccessible by direct simulations due to their complexity, for example, in the highly turbulent regime. However, their accuracy and success depends on a proper subgrid-scale (SGS) model that accounts for the unresolved scales in the simulation. We evaluate the applicability of two traditional SGS models, namely the eddy-viscosity (EV) and the scale-similarity (SS) models, and one recently proposed nonlinear (NL) SGS model in the realm of compressible magnetohydrodynamic (MHD) turbulence. Using 209 simulations of decaying, supersonic (initial sonic Mach number M-s approximate to 3) MHD turbulence with a shock-capturing scheme and varying resolution, SGS model, and filter, we analyze the ensemble statistics of kinetic and magnetic energy spectra and structure functions. Furthermore, we compare the temporal evolution of lower-and higher-order statistical moments of the spatial distributions of kinetic and magnetic energy, vorticity, current density, and dilatation magnitudes. We find no statistical influence on the evolution of the flow by any model if grid-scale quantities are used to calculate SGS contributions. In addition, the SS models, which employ an explicit filter, have no impact in general. On the contrary, both the EV and NL models change the statistics if an explicit filter is used. For example, they slightly increase the dissipation on the smallest scales. We demonstrate that the nonlinear model improves higher-order statistics already with a small explicit filter, i.e., a three-point stencil. The results of, e.g., the structure functions or the skewness and kurtosis of the current density distribution are closer to the ones obtained from simulations at higher resolution. In addition, no additional regularization to stabilize the model is required. We conclude that the nonlinear model with a small explicit filter is suitable for application in more complex scenarios when higher-order statistics are important.

Turbulence in compressible plasma plays a key role in many areas of astrophysics and engineering. The extreme plasma parameters in these environments, e.g. high Reynolds numbers, supersonic and super-Alfvenic flows, however, make direct numerical simulations computationally intractable even for the simplest treatment-magnetohydrodynamics (MHD). To overcome this problem one can use subgrid-scale (SGS) closures-models for the influence of unresolved, subgrid-scales on the resolved ones. In this work we propose and validate a set of constant coefficient closures for the resolved, compressible, ideal MHD equations. The SGS energies are modeled by Smagorinsky-like equilibrium closures. The turbulent stresses and the electromotive force (EMF) are described by expressions that are nonlinear in terms of large scale velocity and magnetic field gradients. To verify the closures we conduct a priori tests over 137 simulation snapshots from two different codes with varying ratios of thermal to magnetic pressure (beta(p) = 0.25, 1, 2.5, 5, 25) and sonic Mach numbers (M-s = 2, 2.5, 4). Furthermore, we make a comparison to traditional, phenomenological eddy-viscosity and alpha- beta - gamma closures. We find only mediocre performance of the kinetic eddy-viscosity and alpha - beta - gamma closures, and that the magnetic eddy-viscosity closure is poorly correlated with the simulation data. Moreover, three of five coefficients of the traditional closures exhibit a significant spread in values. In contrast, our new closures demonstrate consistently high correlations and constant coefficient values over time and over the wide range of parameters tested. Important aspects in compressible MHD turbulence such as the bi-directional energy cascade, turbulent magnetic pressure and proper alignment of the EMF are well described by our new closures.

We derive oscillatory signals in correlation functions in two-field open inflation by means of the in-in formalism; such signatures are caused by resonances between oscillations in the tunnelling field and fluctuations in the inflaton during the curvature dominated, intermediate and subsequent inflationary regime. While amplitudes are model-dependent, we find distinct oscillations in the power and bi-spectrum that can act as a direct probe of the curvature dominated phase and thus, indirectly, strengthen the claim of the string landscape if they were observed. We comment on the prospects of detecting these tell-tale signs in current experiments, which is challenging, but not impossible. At the technical level, we pay special attention to the applicability conditions for truncating fluctuations to the light (inflaton) field and derive upper limits on the oscillation amplitude of the heavy field. A violation of these bounds requires a multi-field analysis at the perturbed level.