Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA systematic comparison of two-equation RANS turbulence models applied to shock-cloud interactions
72
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Turbulence models attempt to account for unresolved dynamics and diffusion in hydrodynamical simulations. We develop a common framework for two-equation Reynolds-Averaged Navier-Stokes (RANS) turbulence models, and we implement six models in the Athena code. We verify each implementation with the standard subsonic mixing layer, although the level of agreement depends on the definition of the mixing layer width. We then test the validity of each model into the supersonic regime, showing that compressibility corrections can improve agreement with experiment. For models with buoyancy effects, we also verify our implementation via the growth of the Rayleigh-Taylor instability in a stratified medium. The models are then applied to the ubiquitous astrophysical shock-cloud interaction in three dimensions. We focus on the mixing of shock and cloud material, comparing results from turbulence models to high-resolution simulations (up to 200 cells per cloud radius) and ensemble-averaged simulations. We find that the turbulence models lead to increased spreading and mixing of the cloud, although no two models predict the same result. Increased mixing is also observed in inviscid simulations at resolutions greater than 100 cells per radius, which suggests that the turbulent mixing begins to be resolved.
rate research
Read More
We study density fluctuations in supersonic turbulence using both theoretical methods and numerical simulations. A theoretical formulation is developed for the probability distribution function (PDF) of the density at steady state, connecting it to the conditional statistics of the velocity divergence. Two sets of numerical simulations are carried out, using either a Riemann solver to evolve the Euler equations or a finite-difference method to evolve the Navier-Stokes (N-S) equations. After confirming the validity of our theoretical formulation with the N-S simulations, we examine the effects of dynamical processes on the PDF, showing that the nonlinear term in the divergence equation amplifies the right tail of the PDF and reduces the left one, the pressure term reduces both the right and left tails, and the viscosity term, counter-intuitively, broadens the right tail of the PDF. Despite the inaccuracy of the velocity divergence from the Riemann runs, as found in our previous work, we show that the density PDF from the Riemann runs is consistent with that from the N-S runs. Taking advantage of their much higher effective resolution, we then use the Riemann runs to study the dependence of the PDF on the Mach number, $mathcal{M}$, up to $mathcal{M}sim30$. The PDF width, $sigma_{s}$, follows the relation $sigma_{s}^2 = ln (1+b^2 {mathcal M}^2)$, with $bapprox0.38$. However, the PDF exhibits a negative skewness that increases with increasing $mathcal{M}$, so much of the growth of $sigma_{s}$ is accounted for by the growth of the left PDF tail, while the growth of the right tail tends to saturate. Thus, the usual prescription that combines a lognormal shape with the standard variance-Mach number relation greatly overestimates the right PDF tail at large $mathcal{M}$, which may have a significant impact on theoretical models of star formation.
The interaction of a shock with a cloud has been extensively studied in the literature, where the effects of magnetic fields, radiative cooling and thermal conduction have been considered. However, the formation of fully developed turbulence has often been prevented by the artificial viscosity inherent in hydrodynamical simulations, and a uniform post-shock flow has been assumed in all previous single-cloud studies. In reality, the flow behind the shock is also likely to be turbulent, with non-uniform density, pressure and velocity structure created as the shock sweeps over inhomogenities upstream of the cloud. To address these twin issues we use a sub-grid compressible k-epsilon turbulence model to estimate the properties of the turbulence generated in shock-cloud interactions and the resulting increase in the transport coefficients that the turbulence brings. A detailed comparison with the output from an inviscid hydrodynamical code puts these new results into context. We find that cloud destruction in inviscid and k-epsilon models occurs at roughly the same speed when the post-shock flow is smooth and when the density contrast between the cloud and inter-cloud medium is less than 100. However, there are increasing and significant differences as this contrast increases. Clouds subjected to strong ``buffeting by a highly turbulent post-shock environment are destroyed significantly quicker. Additional calculations with an inviscid code where the post-shock flow is given random, grid-scale, motions confirms the more rapid destruction of the cloud. Our results clearly show that turbulence plays an important role in shock-cloud interactions, and that environmental turbulence adds a new dimension to the parameter space which has hitherto been studied (abridged).
The rich structure that we observe in molecular clouds is due to the interplay between strong magnetic fields and supersonic (turbulent) velocity fluctuations. The velocity fluctuations interact with the magnetic field, causing it too to fluctuate. Using numerical simulations, we explore the nature of such magnetic field fluctuations, $vec{delta B}$, over a wide range of turbulent Mach numbers, $mathcal{M} = 2 - 20$ (i.e., from weak to strong compressibility), and Alfven Mach numbers, $mathcal{M}_{text{A}0} = 0.1 - 100$ (i.e., from strong to weak magnetic mean fields, $B_0$). We derive a compressible quasi-static fluctuation model from the magnetohydrodynamical (MHD) equations and show that velocity gradients parallel to the mean magnetic field give rise to compressible modes in sub-Alfvenic flows, which prevents the flow from becoming two-dimensional, as is the case in incompressible MHD turbulence. We then generalise an analytical model for the magnitude of the magnetic fluctuations to include $mathcal{M}$, and find $|vec{delta B}| = delta B = c_ssqrt{pirho_0}mathcal{M}mathcal{M}_{text{A}0}$, where $c_s$ is the sound speed and $rho_0$ is the mean density of gas. This new relation fits well in the strong $B$-field regime. We go on to study the anisotropy between the perpendicular ($ B_{perp}$) and parallel ($ B_{parallel}$) fluctuations and the mean-normalised fluctuations, which we find follow universal scaling relations, invariant of $mathcal{M}$. We provide a detailed analysis of the morphology for the $delta B_{perp}$ and $delta B_{parallel}$ probability density functions and find that eddies aligned with $B_0$ cause parallel fluctuations that reduce $B_{parallel}$ in the most anisotropic simulations. We discuss broadly the implications of our fluctuation models for magnetised gases in the interstellar medium.
Using the adaptive mesh refinement code MG, we perform hydrodynamic simulations of the interaction of a shock with a molecular cloud evolving due to thermal instability and gravity. To explore the relative importance of these processes, three case studies are presented. The first follows the formation of a molecular cloud out of an initially quiescent atomic medium due to the effects of thermal instability and gravity. The second case introduces a shock whilst the cloud is still in the warm atomic phase, and the third scenario introduces a shock once the molecular cloud has formed. The shocks accelerate the global collapse of the clouds with both experiencing local gravitational collapse prior to this. When the cloud is still atomic, the evolution is shock dominated and structures form due to dynamical instabilities within a radiatively cooled shell. While the transmitted shock can potentially trigger the thermal instability, this is prevented as material is shocked multiple times on the order of a cloud crushing time-scale. When the cloud is molecular, the post-shock flow is directed via the pre-existing structure through low-density regions in the inter-clump medium. The clumps are accelerated and deformed as the flow induces clump-clump collisions and mergers that collapse under gravity. For a limited period, both shocked cases show a mixture of Kolmogorov and Burgers turbulence-like velocity and logarithmic density power spectra, and strongly varying density spectra. The clouds presented in this work provide realistic conditions that will be used in future feedback studies.
We revisit the problem on the inner structure of shock waves in simple gases modelized by the Boltzmann kinetic equation. In cite{pomeau1987shock}, a self-similarity approach was proposed for infinite total cross section resulting from a power law interaction, but this self-similar form does not have finite energy. Motivated by the work of Pomeau, Bobylev and Cercignani started the rigorous study of the solutions of the spatial homogeneous Boltzmann equation, focusing on those which do not have finite energy cite{bobylev2002self,bobylev2003eternal}. In the present work, we provide a correction to the self-similar form, so that the solutions are more physically sound in the sense that the energy is no longer infinite and that the perturbation brought by the shock does not grow at large distances of it on the cold side in the soft potential case.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
