Implicit large-eddy simulation of compressible flows using the Interior Embedded Discontinuous Galerkin method

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.Comment: 54th AIAA Aerospace Sciences Meeting, AIAA SciTech, 201

Flow separation control over a rounded ramp with spanwise alternating wall actuation

An implicit large-eddy simulation is carried out to study turbulent boundary-layer separation from a backward-facing rounded ramp with active wall actuation control. This method, called spanwise alternating distributed strips control, is imposed onto the flat plate surface upstream of a rounded ramp by alternatively applying out-of-phase control and in-phase control to the wall-normal velocity component in the spanwise direction. As a result, the local turbulence intensity is alternatively suppressed and enhanced, leading to the creation of vertical shear-layers, which is responsible for the presence of large-scale streamwise vortices. These vortices exert a predominant influence on the suppression of the flow separation. The interaction between the large-scale vortices and the downstream recirculation zone and free shear-layer is studied by examining flow statistics. It is found that in comparison with the non-controlled case the flow separation is delayed, the reattachment point is shifted upstream, and the length of the mean recirculation zone is reduced up to 8.49%. The optimal control case is achieved with narrow in-phase control strips. An in-depth analysis shows that the delay of the flow separation is attributed to the activation of the near-wall turbulence by the in-phase control strips and the improvement of the reattachment location is mainly due to the large-scale streamwise vortices, which enhance the momentum transport between the main flow and separated region

Comparison of Subgrid-scale Viscosity Models and Selective Filtering Strategy for Large-eddy Simulations

Explicitly filtered large-eddy simulations (LES), combining high-accuracy schemes with the use of a selective filtering without adding an explicit subgrid-scales (SGS) model, are carried out for the Taylor-Green-vortex and the supersonic-boundary-layer cases. First, the present approach is validated against direct numerical simulation (DNS) results. Subsequently, several SGS models are implemented in order to investigate if they can improve the initial filter-based methodology. It is shown that the most accurate results are obtained when the filtering is used alone as an implicit model, and for a minimal cost. Moreover, the tests for the Taylor-Green vortex indicate that the discretization error from the numerical methods, notably the dissipation error from the high-order filtering, can have a greater influence than the SGS models

Scale-adaptive simulation of unsteady cavitation around a naca66 hydrofoil

Distances between consecutive aftershocks are analysed by means of mono- and multifractal theory with the aim of quantifying the complexity of the physical mechanism governing them, as well as their predictability and predictive instability. Hausdorff, Ha, and Hurst, H, exponents are determined by semivariograms and rescaled analysis, respectively. The exponent ß of the power law describing power spectral contents is also quantified. These three parameters permit a generation of fractional Gaussian noise, fGn, simulating distances. The complexity and predictive instability of physical mechanism generating the series of distances is quantified by means of the correlation dimension, µ*, the Kolmogorov entropy, ¿, and the Lyapunov exponents, ¿i, which are based on the reconstruction theorem formulation. Additionally, the multifractal detrended fluctuation analysis, MF-DFA, contributes with a different point of view to quantify the complexity of the series, in terms of fractal spectral width, W, spectral asymmetry, B, and the critical Hölder exponent, a0. By one hand, the MF-DFA is applied to the complete set of distances characterising the whole aftershock process. By the other hand, the MF-DFA is applied to segments of the series of distances with the aim of determining the evolution of the complexity since the mainshock up to the end of the stress relaxation process. Finally, an ARIMA multilinear regression process is applied to obtain some improvements, in comparison with fGn simulations, on the prediction of distances. The database for this analysis is obtained from the Southern California Seismic Network (SCSN) catalogue. Three series of aftershocks equalling to or exceeding magnitudes of 2.0, assuring seismic catalogue completeness, and associated with Landers (06/28/1992), Northridge (01/17/1994) and Hector Mine (10/16/1999) mainshocks are obtained. It is worth mentioning that common mono-multifractal behaviour for the three aftershocks series is not detected, whatever aftershock periods or segments of them are considered.Postprint (published version

Wall-resolved large eddy simulation over NACA0012 airfoil

The work presented here forms part of a project on Large-Eddy Simulation (LES) of aeroengine aeroacoustic interactions. In this paper we concentrate on LES of near-field flow over an isolated NACA0012 airfoil at zero angle of attack with Rec=2e5. The predicted unsteady pressure/velocity field is used in an analytically-based scheme for far-field trailing edge noise prediction. A wall resolved implicit LES or so-callednumerical Large Eddy Simulation (NLES) approach is employed to resolve streak-like structure in the near-wall flow regions. The mean and RMS velocity and pressure profile on airfoil surface and in wake are validated against experimental data and computational results from other researchers. The results of the wall-resolved NLES method are very encouraging. The effects of grid-refinement and higher-order numerical scheme on the wall-resolved NLES approach are also discussed

Large-eddy simulation and wall modelling of turbulent channel flow

We report large-eddy simulation (LES) of turbulent channel flow. This LES neither resolves nor partially resolves the near-wall region. Instead, we develop a special near-wall subgrid-scale (SGS) model based on wall-parallel filtering and wall-normal averaging of the streamwise momentum equation, with an assumption of local inner scaling used to reduce the unsteady term. This gives an ordinary differential equation (ODE) for the wall shear stress at every wall location that is coupled with the LES. An extended form of the stretched-vortex SGS model, which incorporates the production of near-wall Reynolds shear stress due to the winding of streamwise momentum by near-wall attached SGS vortices, then provides a log relation for the streamwise velocity at the top boundary of the near-wall averaged domain. This allows calculation of an instantaneous slip velocity that is then used as a ‘virtual-wall’ boundary condition for the LES. A Kármán-like constant is calculated dynamically as part of the LES. With this closure we perform LES of turbulent channel flow for Reynolds numbers Re_τ based on the friction velocity u_τ and the channel half-width δ in the range 2 × 10^3 to 2 × 10^7. Results, including SGS-extended longitudinal spectra, compare favourably with the direct numerical simulation (DNS) data of Hoyas & Jiménez (2006) at Re_τ = 2003 and maintain an O(1) grid dependence on Re_τ