Wall modeling via function enrichment within a high‐order DG method for RANS simulations of incompressible flow

We present a novel approach to wall modeling for the Reynolds‐averaged Navier‐Stokes equations within the discontinuous Galerkin method. Wall functions are not used to prescribe boundary conditions as usual, but they are built into the function space of the numerical method as a local enrichment, in addition to the standard polynomial component. The Galerkin method then automatically finds the optimal solution among all shape functions available. This idea is fully consistent and gives the wall model vast flexibility in separated boundary layers or high adverse pressure gradients. The wall model is implemented in a high‐order discontinuous Galerkin solver for incompressible flow complemented by the Spalart‐Allmaras closure model. As benchmark examples, we present turbulent channel flow starting from Re_τ=180 and up to Re_τ=100000 as well as flow past periodic hills at Reynolds numbers based on the hill height of Re_H=10595 and Re_H=19000.     

Prediction of separation flows around a 6：1 prolate spheroid using RANS/LES hybrid approaches

This paper presents hybrid Reynolds-averaged Navier–Stokes (RANS) and large-eddy-simulation (LES) methods for the separated flows at high angles of attack around a 6:1 prolate spheroid. The RANS/LES hybrid methods studied in this work include the detached eddy simulation (DES) based on Spalart–Allmaras (S–A), Menter’s k–ω shear-stress-transport (SST) and k–ω with weakly nonlinear eddy viscosity formulation (Wilcox–Durbin+, WD+) models and the zonal-RANS/LES methods based on the SST and WD+ models. The switch from RANS near the wall to LES in the core flow region is smooth through the implementation of a flow-dependent blending function for the zonal hybrid method. All the hybrid methods are designed to have a RANS mode for the attached flows and have a LES behavior for the separated flows. The main objective of this paper is to apply the hybrid methods for the high Reynolds number separated flows around prolate spheroid at high-incidences. A fourth-order central scheme with fourth-order artificial viscosity is applied for spatial differencing. The fully implicit lower–upper symmetric-Gauss–Seidel with pseudo time sub-iteration is taken as the temporal differentiation. Comparisons with available measurements are carried out for pressure distribution, skin friction, and profiles of velocity, etc. Reasonable agreement with the experiments, accounting for the effect on grids and fundamental turbulence models, is obtained for the separation flows. The project supported by the National Natural Science Foundation of China (10502030 and 90505005).     

A hybrid RANS/LES approach with scale-adaptive capabilities for highly separated flows

A new hybrid RANS/LES approach with scale-adaptive capabilities is developed. The blending function in the SST model is adopted to prevent the invasion of the von Karman length scale to the RANS region, and the compressibility correction proposed by Wilcox is incorporated to produce a realistic shear layer development in compressible flows. The new model is validated for a subcritical flow past a circular cylinder and a supersonic base flow. Time-averaged turbulent statistics predicted by the new model show fairly good agreement with the experimental data, slight improvements over DES simulations, and are much better than SAS results. The main advantage of the new model over the DES method is that the distribution of the blending function reflects local vortex structures instead of grid spacing in the turbulent wake. The sequence of the effect intensity of the compressibility correction from strong to weak is SAS, the new model and DES.     

A high-order nodal discontinuous Galerkin method to solve preconditioned multiphase Euler/Navier-Stokes equations for inviscid/viscous cavitating flows

In this study, a high-order accurate numerical method is applied and examined for the simulation of the inviscid/viscous cavitating flows by solving the preconditioned multiphase Euler/Navier-Stokes equations on triangle elements. The formulation used here is based on the homogeneous equilibrium model considering the continuity and momentum equations together with the transport equation for the vapor phase with applying appropriate mass transfer terms for calculating the evaporation/condensation of the liquid/vapor phase. The spatial derivative terms in the resulting system of equations are discretized by the nodal discontinuous Galerkin method (NDGM) and an implicit dual-time stepping method is used for the time integration. An artificial viscosity approach is implemented and assessed for capturing the steep discontinuities in the interface between the two phases. The accuracy and robustness of the proposed method in solving the preconditioned multiphase Euler/Navier-Stokes equations are examined by the simulation of different two-dimensional and axisymmetric cavitating flows. A sensitivity study is also performed to examine the effects of different numerical parameters on the accuracy and performance of the solution of the NDGM. Indications are that the solution methodology proposed and applied here is based on the NDGM with the implicit dual-time stepping method and the artificial viscosity approach is accurate and robust for the simulation of the inviscid and viscous cavitating flows.