A new approach of a limiting process for multi-dimensional flows

An enhanced Multi-dimensional Limiting Process (e-MLP) is developed for the accurate and efficient computation of multi-dimensional flows based on the Multi-dimensional Limiting Process (MLP). The new limiting process includes a distinguishing step and an enhanced multi-dimensional limiting process. First, the distinguishing step, which is independent of high order interpolation and flux evaluation, is newly introduced. It performs a multi-dimensional search of a discontinuity. The entire computational domain is then divided into continuous, linear discontinuous and nonlinear discontinuous regions. Second, limiting functions are appropriately switched according to the type of each region; in a continuous region, there are no limiting processes and only higher order accurate interpolation is performed. In linear discontinuous and nonlinear discontinuous regions, TVD criterion and MLP limiter are respectively used to remove oscillation. Hence, e-MLP has a number of advantages, as it incorporates useful features of MLP limiter such as multi-dimensional monotonicity and straightforward extensionality to higher order interpolation. It is applicable to local extrema and prevents excessive damping in a linear discontinuous region through application of appropriate limiting criteria. It is efficient because a limiting function is applied only to a discontinuous region. In addition, it is robust against shock instability due to the strict distinction of the computational domain and the use of regional information in a flux scheme as well as a high order interpolation scheme. This new limiting process was applied to numerous test cases including one-dimensional shock/sine wave interaction problem, oblique stationary contact discontinuity, isentropic vortex flow, high speed flow in a blunt body, planar shock/density bubble interaction, shock wave/vortex interaction and, particularly, magnetohydrodynamic (MHD) cloud-shock interaction problems. Through these tests, it was verified that e-MLP substantially enhances the accuracy and efficiency with both continuous and discontinuous multi-dimensional flows.

[1]  P. Janhunen,et al.  A Positive Conservative Method for Magnetohydrodynamics Based on HLL and Roe Methods , 2000 .

[2]  G. Tóth The ∇·B=0 Constraint in Shock-Capturing Magnetohydrodynamics Codes , 2000 .

[3]  Chongam Kim,et al.  Multi-dimensional limiting process for three-dimensional flow physics analyses , 2008, J. Comput. Phys..

[4]  M. Brio,et al.  An upwind differencing scheme for the equations of ideal magnetohydrodynamics , 1988 .

[5]  P. Sweby High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws , 1984 .

[6]  Sanghoon Han,et al.  Accurate and Robust Pressure Weight Advection Upstream Splitting Method for Magnetohydrodynamics Equations , 2009 .

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

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

[9]  C. Hirsch,et al.  Numerical Computation of Internal and External Flows. By C. HIRSCH. Wiley. Vol. 1, Fundamentals of Numerical Discretization. 1988. 515 pp. £60. Vol. 2, Computational Methods for Inviscid and Viscous Flows. 1990, 691 pp. £65. , 1991, Journal of Fluid Mechanics.

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

[11]  Kyu Hong Kim,et al.  Accurate, efficient and monotonic numerical methods for multi-dimensional compressible flows Part II: Multi-dimensional limiting process , 2005 .

[12]  Chongam Kim,et al.  Cures for the shock instability: development of a shock-stable Roe scheme , 2003 .

[13]  Elaine S. Oran,et al.  The interaction of a shock with a vortex: Shock distortion and the production of acoustic waves , 1995 .

[14]  Neil D. Sandham,et al.  Low-Dissipative High-Order Shock-Capturing Methods Using Characteristic-Based Filters , 1999 .

[15]  Chongam Kim,et al.  Accurate, efficient and monotonic numerical methods for multi-dimensional compressible flows Part I Spatial discretization , 2005 .

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

[17]  Oh-Hyun Rho,et al.  Methods for the accurate computations of hypersonic flows: I. AUSMPW + scheme , 2001 .

[18]  Sungtae Kim,et al.  Wavenumber-extended high-order oscillation control finite volume schemes for multi-dimensional aeroacoustic computations , 2008, J. Comput. Phys..

[19]  J. P. Boris,et al.  Vorticity generation by shock propagation through bubbles in a gas , 1988, Journal of Fluid Mechanics.

[20]  Y. Hattori,et al.  Sound generation by shock–vortex interactions , 1999, Journal of Fluid Mechanics.

[21]  Chi-Wang Shu TVB uniformly high-order schemes for conservation laws , 1987 .

[22]  Charles Hirsch,et al.  Numerical computation of internal & external flows: fundamentals of numerical discretization , 1988 .

[23]  K Ravi,et al.  Some aspects of high-speed blunt body flow computations with Roe scheme , 1995 .