Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userLarge-eddy simulation of hypersonic flows. Selective procedure to activate the sub-grid model only where small scale turbulence is present
99
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A new method for the localization of the regions where small scale turbulent fluctuations are present in hypersonic flows is applied to the large-eddy simulation (LES) of a compressible turbulent jet with an initial Mach number equal to 5. The localization method used is called selective LES and is based on the exploitation of a scalar probe function $f$ which represents the magnitude of the stretching-tilting term of the vorticity equation normalized with the enstrophy (Tordella et al. 2007). For a fully developed turbulent field of fluctuations, statistical analysis shows that the probability that $f$ is larger than 2 is almost zero, and, for any given threshold, it is larger if the flow is under-resolved. By computing the spatial field of $f$ in each instantaneous realization of the simulation it is possible to locate the regions where the magnitude of the normalized vortical stretching-tilting is anomalously high. The sub-grid model is then introduced into the governing equations in such regions only. The results of the selective LES simulation are compared with those of a standard LES, where the sub-grid terms are used in the whole domain, and with those of a standard Euler simulation with the same resolution. The comparison is carried out by assuming as reference field a higher resolution Euler simulation of the same jet. It is shown that the selective LES modifies the dynamic properties of the flow to a lesser extent with respect to the classical LES. In particular, the prediction of the enstrophy, mean velocity and density distributions and of the energy and density spectra are substantially improved.
rate research
Read More
We present a high-order implicit large-eddy simulation (ILES) approach for simulating transitional turbulent flows. The approach consists of an Interior Embedded Discontinuous Galerkin (IEDG) method for the discretization of the compressible Navier-Stokes equations and a parallel preconditioned Newton-GMRES solver for the resulting nonlinear system of equations. The IEDG method arises from the marriage of the Embedded Discontinuous Galerkin (EDG) method and the Hybridizable Discontinuous Galerkin (HDG) method. As such, the IEDG method inherits the advantages of both the EDG method and the HDG method to make itself well-suited for turbulence simulations. We propose a minimal residual Newton algorithm for solving the nonlinear system arising from the IEDG discretization of the Navier-Stokes equations. The preconditioned GMRES algorithm is based on a restricted additive Schwarz (RAS) preconditioner in conjunction with a block incomplete LU factorization at the subdomain level. The proposed approach is applied to the ILES of transitional turbulent flows over a NACA 65-(18)10 compressor cascade at Reynolds number 250,000 in both design and off-design conditions. The high-order ILES results show good agreement with a subgrid-scale LES model discretized with a second-order finite volume code while using significantly less degrees of freedom. This work shows that high-order accuracy is key for predicting transitional turbulent flows without a SGS model.
MHD turbulence is likely to play an important role in several astrophysical scenarios where the magnetic Reynolds is very large. Numerically, these cases can be studied efficiently by means of Large Eddy Simulations, in which the computational resources are used to evolve the system only up to a finite grid size. The resolution is not fine enough to capture all the relevant small-scale physics at play, which is instead effectively modeled by a set of additional terms in the evolution equations, dubbed as sub-grid-scale model. Here we extend such approach, commonly used in non-relativistic/non-magnetic/incompressible fluid dynamics, applying the so-called gradient model to a general set of balance-law equations, that includes the relevant case in which a non-trivial inversion of conserved to primitive fields is needed. In particular, we focus on the relativistic compressible ideal MHD scenario, providing for the first time (and for any equation of state) all the additional sub-grid-scale terms. As an application, we consider box simulations of the relativistic Kelvin-Helmholtz instability, which is also the first mechanism responsible for the magnetic field amplification in binary neutron star mergers and cannot yet be fully captured by the finest-grid and longest simulations available. The performance of our model is numerically assessed by comparing it to the residuals arising from the filtering of high-resolution simulations. We find that the model can fit very well those residuals from resolutions a few times higher. Although the application shown here explicitly considers the Minkowski metric, it can be directly extended to general relativity, thus settling the basis to implement the gradient sub-grid model in a GRMHD binary merger. Our results suggest that this approach will be potentially able to unveil much better the small-scale dynamics achievable in full GRMHD simulations.
We present a model describing evolution of the small-scale Navier-Stokes turbulence due to its stochastic distortions by much larger turbulent scales. This study is motivated by numerical findings (laval, 2001) that such interactions of separated scales play important role in turbulence intermittency. We introduce description of turbulence in terms of the moments of the k-space quantities using a method previously developed for the kinematic dynamo problem (Nazarenko, 2003). Working with the $k$-space moments allows to introduce new useful measures of intermittency such as the mean polarization and the spectral flatness. Our study of the 2D turbulence shows that the energy cascade is scale invariant and Gaussian whereas the enstrophy cascade is intermittent. In 3D, we show that the statistics of turbulence wavepackets deviates from gaussianity toward dominance of the plane polarizations. Such turbulence is formed by ellipsoids in the $k$-space centered at its origin and having one large, one neutral and one small axes with the velocity field pointing parallel to the smallest axis.
A nonlocal subgrid-scale stress (SGS) model is developed based on the convolution neural network (CNN), a powerful supervised data-driven approach. The CNN is an ideal approach to naturally consider nonlocal spatial information in prediction due to its wide receptive field. The CNN-based models used here only take primitive flow variables as input, then the flow features are automatically extracted without any $priori$ guidance. The nonlocal models trained by direct numerical simulation (DNS) data of a turbulent channel flow at $Re_{tau}=178$ are accessed in both the $priori$ and $posteriori$ test, providing physically reasonable flow statistics (like mean velocity and velocity fluctuations) closing to the DNS results even when extrapolating to a higher Reynolds number $Re_{tau}=600$. In our model, the backscatter is also predicted well and the numerical simulation is stable. The nonlocal models outperform local data-driven models like artificial neural network and some SGS models, e.g. the Smagorinsky model in actual large eddy simulation (LES). The model is also robust since stable solutions can be obtained when examining the grid resolution from one-half to double of the spatial resolution used in training. We also investigate the influence of receptive fields and suggest using the two-point correlation analysis as a quantitative method to guide the design of nonlocal physical models. To facilitate the combination of machine learning (ML) algorithms to computational fluid dynamics (CFD), a novel heterogeneous ML-CFD framework is proposed. The present study provides the effective data-driven nonlocal methods for SGS modelling in the LES of complex anisotropic turbulent flows.
In this essay, we recall the specificities of the transition to turbulence in wall-bounded flows and present recent achievements in the understanding of this problem. The transition is abrupt with laminar-turbulent coexistence over a finite range of Reynolds numbers, the transitional range. The archetypical cases of Poiseuille pipe flow and plane Couette flow are first reviewed at the phenomenological level, together with a few other flow configurations. Theoretical approaches are then examined with particular emphasis on the existence of special nontrivial solutions to the Navier-Stokes equations at finite distance from laminar flow. Dynamical systems theory is most appropriate to analyze their role, in particular with respect to the transient character of turbulence in the lower transitional range. The extensions needed to deal with the prominent spatiotemporal features of the transition are then discussed. Turbulence growth/decay in terms of statistical physics of many-body systems and the relevance of directed percolation as a stochastic process able to account for it are next scrutinized. To conclude, we advocate the recourse to well-designed modeling able to provide us with a conceptually coherent picture of the full transitional range and put forward some open issues.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
