Low dissipative shock-capturing methods using weighted essentially nonoscillatory and central high resolution filters

Author(s):  
V. Weirs
2013 ◽  
Vol 11 (01) ◽  
pp. 1350049
Author(s):  
M. P. RAY ◽  
B. P. PURANIK ◽  
U. V. BHANDARKAR

High-resolution extensions to six Riemann solvers and three flux vector splitting schemes are developed within the framework of a reconstruction-evolution approach. Third-order spatial accuracy is achieved using two different piecewise parabolic reconstructions and a weighted essentially nonoscillatory scheme. A three-stage TVD Runge–Kutta time stepping is employed for temporal integration. The modular development of solvers provides an ease in selecting a reconstruction scheme and/or a Riemann solver/flux vector splitting scheme. The performances of these high-resolution solvers are compared for several one- and two-dimensional test cases. Based on a comprehensive assessment of the solutions obtained with all solvers, it is found that the use of the weighted essentially nonoscillatory reconstruction with the van Leer flux vector splitting scheme provides solutions for a variety of problems with acceptable accuracy.


1994 ◽  
Vol 258 ◽  
pp. 317-333 ◽  
Author(s):  
José Ma Martí ◽  
Ewald Müller

We consider the decay of an initial discontinuity in a polytropic gas in a Minkowski space–time (the special relativistic Riemann problem). In order to get a general analytical solution for this problem, we analyse the properties of the relativistic flow across shock waves and rarefactions. As in classical hydrodynamics, the solution of the Riemann problem is found by solving an implicit algebraic equation which gives the pressure in the intermediate states. The solution presented here contains as a particular case the special relativistic shock-tube problem in which the gas is initially at rest. Finally, we discuss the impact of this result on the development of high-resolution shock-capturing numerical codes to solve the equations of relativistic hydrodynamics.


2014 ◽  
Vol 2014 ◽  
pp. 1-14 ◽  
Author(s):  
Lang Wu ◽  
Dazhi Zhang ◽  
Boying Wu ◽  
Xiong Meng

Fifth-order mapped semi-Lagrangian weighted essentially nonoscillatory (WENO) methods at certain smooth extrema are developed in this study. The schemes contain the mapped semi-Lagrangian finite volume (M-SL-FV) WENO 5 method and the mapped compact semi-Lagrangian finite difference (M-C-SL-FD) WENO 5 method. The weights in the more common scheme lose accuracy at certain smooth extrema. We introduce mapped weighting to handle the problem. In general, a cell average is applied to construct the M-SL-FV WENO 5 reconstruction, and the M-C-SL-FD WENO 5 interpolation scheme is proposed based on an interpolation approach. An accuracy test and numerical examples are used to demonstrate that the two schemes reduce the loss of accuracy and improve the ability to capture discontinuities.


2019 ◽  
Vol 80 (1) ◽  
pp. 555-613 ◽  
Author(s):  
Hassan Yousefi ◽  
Soheil Mohammadi ◽  
Timon Rabczuk

2013 ◽  
Vol 2013 ◽  
pp. 1-11 ◽  
Author(s):  
Chunlei Ruan ◽  
Kunfeng Liang ◽  
Xianjie Chang ◽  
Ling Zhang

Population balance equations (PBEs) are the main governing equations to model the processes of crystallization. Two-dimensional PBEs refer to the crystals that grow anisotropically with the change of two internal coordinates. Since the PBEs are hyperbolic equations, it is necessary to build up high resolution schemes to avoid numerical diffusion and numerical dispersion in order to obtain the accurate crystal size distribution (CSD). In this work, a 5th order weighted essentially nonoscillatory (WENO) method is introduced to compute the two-dimensional PBEs. Several numerical benchmark examples from literatures are carried out; it is found that WENO method has higher resolution than HR method which is well established. Therefore, WENO method is recommended in crystallization simulation when the crystal size distributions are sharp and higher accuracy is needed.


Sign in / Sign up

Export Citation Format

Share Document