Universal High Order Subroutine with New Shock Detector for Shock Boundary Layer Interaction

The goal of this work is to develop a new universal high order subroutine for shock boundary layer interaction. First, an effective shock/discontinuity detector has been developed.The detector has two steps.The first step is to check the ratio of the truncation errors on the coarse and fine grids and the second step is to check the local ratio of the left and right slopes. The currently popular shock/discontinuity detectors can detect shock, but mistake high frequency waves and critical points as shock and then damp the physically important high frequency waves.Preliminary results show the new shock/discontinuity detector is very delicate and can detect all shocks including strong, weak and oblique shocks or discontinuity in function and the first, second, and third order derivatives without artificial constants, but never mistake high frequency waves and critical points, expansion waves as shock. This will overcome the bottle neck problem with numerical simulation for the shock-boundary layer interaction, shock-acoustic interaction, image process, porous media flow, multiple phase flow and anywhere the high frequency waves are important, but discontinuity exists and is mixed with high frequency waves. After detecting the shock we can then use one side high order scheme for shocks and high order central compact scheme for the rest if the shock is appropriately located. Then a high order universal subroutine for finite difference method is developed which can be used for any finite difference code for accurate numerical derivatives.

[1]  S. Osher,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes,II , 1989 .

[2]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[3]  Jianzhong Su,et al.  High Order Weighted Compact Scheme for Shock/Boundary Layer Interaction , 2008 .

[4]  M. Pino Martin,et al.  Direct numerical simulation of hypersonic turbulent boundary layers. Part 1. Initialization and comparison with experiments , 2007, Journal of Fluid Mechanics.

[5]  Fermín Navarrina,et al.  Finite volume solvers and movingleast-squares approximations for thecompressible Navier-Stokes equations onunstructured grids , 2007 .

[6]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[7]  S. Lele Compact finite difference schemes with spectral-like resolution , 1992 .

[8]  P. Roe Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[9]  Rick Archibald,et al.  Determining the locations and discontinuities in the derivatives of functions , 2008 .

[10]  Feng He,et al.  A new family of high-order compact upwind difference schemes with good spectral resolution , 2007, J. Comput. Phys..

[11]  Gérard Degrez,et al.  The interaction of an oblique shock wave with a laminar boundary layer revisited. An experimental and numerical study , 1987, Journal of Fluid Mechanics.

[12]  A. Harten High Resolution Schemes for Hyperbolic Conservation Laws , 2017 .

[13]  A. Harten,et al.  The artificial compression method for computation of shocks and contact discontinuities: III. Self , 1978 .

[14]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .