Numerical analysis of conservative unstructured discretisations for low Mach flows

This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving. https://authorservices.wiley.com/author-resources/Journal-Authors/licensing-and-open-access/open-access/self-archiving.htmlUnstructured meshes allow easily representing complex geometries and to refine in regions of interest without adding control volumes in unnecessary regions. However, numerical schemes used on unstructured grids have to be properly defined in order to minimise numerical errors. An assessment of a low-Mach algorithm for laminar and turbulent flows on unstructured meshes using collocated and staggered formulations is presented. For staggered formulations using cell centred velocity reconstructions the standard first-order method is shown to be inaccurate in low Mach flows on unstructured grids. A recently proposed least squares procedure for incompressible flows is extended to the low Mach regime and shown to significantly improve the behaviour of the algorithm. Regarding collocated discretisations, the odd-even pressure decoupling is handled through a kinetic energy conserving flux interpolation scheme. This approach is shown to efficiently handle variable-density flows. Besides, different face interpolations schemes for unstructured meshes are analysed. A kinetic energy preserving scheme is applied to the momentum equations, namely the Symmetry-Preserving (SP) scheme. Furthermore, a new approach to define the far-neighbouring nodes of the QUICK scheme is presented and analysed. The method is suitable for both structured and unstructured grids, either uniform or not. The proposed algorithm and the spatial schemes are assessed against a function reconstruction, a differentially heated cavity and a turbulent self-igniting diffusion flame. It is shown that the proposed algorithm accurately represents unsteady variable-density flows. Furthermore, the QUICK schemes shows close to second order behaviour on unstructured meshes and the SP is reliably used in all computations.Peer ReviewedPostprint (author's final draft

Modelling the rapid spherical compression of isotropic turbulence

Rapid distortion analysis is used to modify the form of the closure model for the dissipation rate of the turbulent kinetic energy. The modification is such that the evolution of the dissipation rate during a rapid compression is predicted exactly; good agreement between the model prediction and direct simulation data is obtained. Previous closure proposals fail to properly predict the rapid compression case. The reason for the difference between the present and previous models is traced to the fact that previous workers neglected variations of kinematic viscosit

Adaptively refined large eddy simulations of clusters

We present a numerical scheme for modelling unresolved turbulence in cosmological adaptive mesh refinement codes. As a first application, we study the evolution of turbulence in the intra-cluster medium and in the core of a galaxy cluster. Simulations with and without subgrid scale model are compared in detail. Since the flow in the ICM is subsonic, the global turbulent energy contribution at the unresolved length scales is smaller than 1% of the internal energy. We find that the production of turbulence is closely correlated with merger events occurring in the cluster environment, and its dissipation locally affects the cluster energy budget. Because of this additional source of dissipation, the core temperature is larger and the density is smaller in the presence of subgrid scale turbulence than in the standard adiabatic run, resulting in a higher entropy core value.Comment: Submitted to ApJ, 14 pages, 14 figures, 3 table

A Fluid-Dynamical Subgrid Scale Model for Highly Compressible Astrophysical Turbulence

We formulate and implement the Euler equations with SGS dynamics and provide numerical tests of an SGS turbulence energy model that predicts the turbulent pressure of unresolved velocity fluctuations and the rate of dissipation for highly compressible turbulence. We test closures for the turbulence energy cascade by filtering data from high-resolution simulations of forced isothermal and adiabatic turbulence. Optimal properties and an excellent correlation are found for a linear combination of the eddy-viscosity closure that is employed in LES of weakly compressible turbulence and a term that is non-linear in the Jacobian matrix of the velocity. Using this mixed closure, the SGS turbulence energy model is validated in LES of turbulence with stochastic forcing. It is found that the SGS model satisfies several important requirements: 1. The mean SGS turbulence energy follows a power law for varying grid scale. 2. The root mean square (RMS) Mach number of the unresolved velocity fluctuations is proportional to the RMS Mach number of the resolved turbulence, independent of the forcing. 3. The rate of dissipation and the turbulence energy flux are constant. Moreover, we discuss difficulties with direct estimates of the turbulent pressure and the dissipation rate on the basis of resolved flow quantities that have recently been proposed. In combination with the energy injection by stellar feedback and other unresolved processes, the proposed SGS model is applicable to a variety of problems in computational astrophysics. Computing the SGS turbulence energy, the treatment of star formation and stellar feedback in galaxy simulations can be improved. Further, we expect that the turbulent pressure on the grid scale affects the stability of gas against gravitational collapse.Comment: 19 pages, 16 figures, submitted to A&

Toward the large-eddy simulation of compressible turbulent flows

New subgrid-scale models for the large-eddy simulation of compressible turbulent flows are developed and tested based on the Favre-filtered equations of motion for an ideal gas. A compressible generalization of the linear combination of the Smagorinsky model and scale-similarity model, in terms of Favre-filtered fields, is obtained for the subgrid-scale stress tensor. An analogous thermal linear combination model is also developed for the subgrid-scale heat flux vector. The two dimensionless constants associated with these subgrid-scale models are obtained by correlating with the results of direct numerical simulations of compressible isotropic turbulence performed on a 96(exp 3) grid using Fourier collocation methods. Extensive comparisons between the direct and modeled subgrid-scale fields are provided in order to validate the models. A large-eddy simulation of the decay of compressible isotropic turbulence (conducted on a coarse 32(exp 3) grid) is shown to yield results that are in excellent agreement with the fine grid direct simulation. Future applications of these compressible subgrid-scale models to the large-eddy simulation of more complex supersonic flows are discussed briefly

The analysis and simulation of compressible turbulence

Compressible turbulent flows at low turbulent Mach numbers are considered. Contrary to the general belief that such flows are almost incompressible, (i.e., the divergence of the velocity field remains small for all times), it is shown that even if the divergence of the initial velocity field is negligibly small, it can grow rapidly on a non-dimensional time scale which is the inverse of the fluctuating Mach number. An asymptotic theory which enables one to obtain a description of the flow in terms of its divergence-free and vorticity-free components has been developed to solve the initial-value problem. As a result, the various types of low Mach number turbulent regimes have been classified with respect to the initial conditions. Formulae are derived that accurately predict the level of compressibility after the initial transients have disappeared. These results are verified by extensive direct numerical simulations of isotropic turbulence