Recent developments in accuracy and stability improvement of nonlinear filter methods for DNS and LES of compressible flows

Abstract

High order methods High order shock-capturing methods Stability improvement of high order methods Accuracy improvement for DNS & LES a b s t r a c t Recent progress in the improvement of numerical stability and accuracy of the Yee and Sjögreen [49] high order nonlinear filter schemes is described. The Yee & Sjögreen adaptive nonlinear filter method consists of a high order non-dissipative spatial base scheme and a nonlinear filter step. The nonlinear filter step consists of a flow sensor and the dissipative portion of a high resolution nonlinear high order shock-capturing method to guide the application of the shock-capturing dissipation where needed. The nonlinear filter idea was first initiated by Yee et al. [54] using an artificial compression method (ACM) of Harten [12] as the flow sensor. The nonlinear filter step was developed to replace high order linear filters so that the same scheme can be used for long time integration of direct numerical simulations (DNS) and large eddy simulations (LES) for both shock-free turbulence and turbulence-shock waves interactions. The improvement includes four major new developments: (a) Smart flow sensors were developed to replace the global ACM flow sensor [21,22,50]. The smart flow sensor provides the locations and the estimated strength of the necessary numerical dissipation needed at these locations and leaves the rest of the flow field free of shock-capturing dissipation. (b) Skew-symmetric splittings were developed for compressible gas dynamics and magnetohydrodynamics (MHD) equations [35,36] to improve numerical stability for long time integration. (c) High order entropy stable numerical fluxes were developed as the spatial base schemes for both the compressible gas dynamics and MHD [37,38]. (d) Several dispersion relation-preserving (DRP) central spatial schemes were included as spatial base schemes in the framework of our nonlinear filter method approach [40]. With these new scheme constructions the nonlinear filter schemes are applicable to a wider class of accurate and stable DNS and LES applications, including forced turbulence simulations where the time evolution of flows might start with low speed shock-free turbulence and develop into supersonic speeds with shocks. Representative test cases for both smooth flows and problems containing discontinuities for compressible flows are included.

On the performance of WENO/TENO schemes to resolve turbulence in DNS/LES of high‐speed compressible flows

“…The blending approach of the abovementioned upwind/central‐WENO hybrid schemes requires a great care to ensure numerical conservation and stability at the interface of different schemes used. Non‐linear filter schemes 28‐31 are developed to avoid this interface deficiency of high‐order hybrid schemes. Specifically, a hybrid scheme is applied to the inviscid fluxes and the nonlinear fluxes are typically evaluated by central discretization of an entropy conservative formulation.…”

Numerical experiments on a hybrid WENO5 filter for shock‐capturing

“…Further works on this approach are the works carried out by Bogey et al [17], Visbal et al [18], Mahmoodi et al [19] and Kang et al [20]. Recently, Kotov et al [21][22][23][24] and also Yee and Sjogreen [25][26][27] developed nonlinear filtering approach for DNS and LES applications. An important advantage of the filtering approach over the other methods such as WENO or hybrid schemes (which require Riemann solvers), is that the filter is applied once a single time step which results in low computational cost when using a multi stage Runge-Kutta scheme for the time discretization.…”