A robustness-enhanced method for Riemann solver

Shaolong Guo,Wen-Quan Tao

DOI: https://doi.org/10.1016/j.ijheatmasstransfer.2020.120757

IF: 5.2

2021-01-01

International Journal of Heat and Mass Transfer

Abstract:•An extra term is added to the low diffusion schemes to overcome the shock anomaly.•The revised schemes turn out to be carbuncle-free and shock-stable.•The accuracy of the original schemes in the boundary-layer is not contaminated.

What problem does this paper attempt to address?

An Arbitrary Lagrangian–Eulerian Method Based on the Adaptive Riemann Solvers for General Equations of State

Baolin Tian,Weidong Shen,Song Jiang,Shuanghu Wang,Yan Liu

DOI: https://doi.org/10.1002/fld.1871

IF: 1.847

2008-01-01

International Journal for Numerical Methods in Fluids

Abstract:Approximate or exact Riemann solvers play a key role in Godunov‐type methods. In this paper, three approximate Riemann solvers, the MFCAV, DKWZ and weak wave approximation method schemes, are investigated through numerical experiments, and their numerical features, such as the resolution for shock and contact waves, are analyzed and compared. Based on the analysis, two new adaptive Riemann solvers for general equations of state are proposed, which can resolve both shock and contact waves well. As a result, an ALE method based on the adaptive Riemann solvers is formulated. A number of numerical experiments show good performance of the adaptive solvers in resolving both shock waves and contact discontinuities. Copyright © 2008 John Wiley & Sons, Ltd.
An Improved Meshfree Scheme Based on Radial Basis Functions for Solving Incompressible Navier–Stokes Equations

Yulin Xie,Xizeng Zhao,Matteo Rubinato,Yuemin Yu

DOI: https://doi.org/10.1002/fld.5012

IF: 1.847

2021-01-01

International Journal for Numerical Methods in Fluids

Abstract:In this paper, an improved meshfree scheme based on radial basis functions (RBFs) is provided for solving the incompressible viscous Navier–Stokes equations and two enhancements are proposed to mitigate the typical numerical oscillations. The first one is the combination of the RBFs‐based finite difference (RBF‐FD) method with the semi‐Lagrangian RBFs (SLM‐RBF), with the former being used for the viscous diffusion term and pressure Poisson equation and the latter being used for the advection term. The second enhancement is a regularization term that constructs smooth constraints for RBFs interpolations instead of clipping operations. The capability of the proposed scheme in mitigating numerical fluctuations is demonstrated by validating it against the one‐dimensional (1‐D) advection problem and the advection–diffusion problem with step field functions. The overall performance of the proposed scheme is also validated by the lid‐driven cavity flow and laminar flow around a circular cylinder, showing good agreement with the existing results, indicating that the proposed scheme has good stability in both temporal and spatial domains.
An Optimized Dispersion-Relation-Preserving Combined Compact Difference Scheme to Solve Advection Equations

C. H. Yu,D. Wang,Z. He,T. Paehtz

DOI: https://doi.org/10.1016/j.jcp.2015.07.051

IF: 4.645

2015-01-01

Journal of Computational Physics

Abstract:In this study, we first present an improved version of the classical sixth-order combined compact difference (CCD6) scheme to enhance the convective stability of advection equations through an increased dispersion accuracy. This improved fifth-order dispersion–relation-preserving combined compact difference scheme (DRPCCD5) has been rigorously analyzed through the dispersion, phase speed anisotropy and stability analyses. We then couple the DRPCCD5 scheme with the previous fifth-order compact-reconstruction weighted essentially non-oscillatory (CRWENO5) scheme using a hybrid strategy based on the monotonicity–maintenance criteria. To verify the resulting “optimized” hybrid scheme (ODRPCCD5), several benchmark problems with available exact solution are investigated. The comparison to the previous fifth-order WENO (WENO5) scheme shows that the ODRPCCD5 avoids numerical oscillation around discontinuities, handles large gradients well, and is much faster at the same accuracy because a coarser mesh can be used.
A robust shock-capturing scheme based on rotated Riemann solvers

Yu-Xin Ren

DOI: https://doi.org/10.1016/S0045-7930(02)00114-7

IF: 3.077

2003-01-01

Computers & Fluids

Abstract:This paper presents a robust finite volume shock-capturing scheme based on the rotated approximate Riemann solver. A general framework for constructing the rotated Riemann solver is described and a rotated Roe scheme is discussed in detail. It is found that the robustness of the rotated shock-capturing scheme is closely related to the way in which the direction of upwind differencing is determined. When the upwind direction is determined by the velocity-difference vector, the rotated Roe scheme demonstrates a robust shock-capturing capability and the shock instabilities or carbuncle phenomena can be eliminated completely. The dissipation property associated with the linear field of the rotated flux-difference splitting scheme is analyzed, and several test cases are presented to validate the proposed scheme.
Mechanism-Derived Shock Instability Elimination for Riemann-Solver-Based Shock-Capturing Scheme

Zhiqiang Chen,Xudong Huang,Yu-Xin Ren,Zhifeng Xie,Ming Zhou

DOI: https://doi.org/10.2514/1.j056882

2018-01-01

Abstract:In our study series, the first paper found that the negative pressure dissipation of density flux at the vertical transverse face of a shock leads to shock instability. Here, we present a general eliminating method derived by this mechanism. This method, named transverse pressure dissipation modification, explicitly adds additional flux dissipation at the vertical transverse face of a shock. The stability of the modified scheme was analyzed by a linear stability method, and completely negative values of maximum real part eigenvalues of the linearized system were realized, which indicates the effectiveness of this method. This method was verified and compared in detail with the existing elimination method for typical cases via structured/unstructured mesh and two-dimensional and three-dimensional configurations. Transverse pressure dissipation modification can be connected to any Riemann solver, and the combined scheme is numerically shock stable and contact preserving.
General Procedure for Riemann Solver to Eliminate Carbuncle and Shock Instability

Zhiqiang Chen,Xudong Huang,Yu-Xin Ren,Ming Zhou

DOI: https://doi.org/10.2514/1.j055366

2017-01-01

Abstract:A novel procedure for the Riemann solver flux calculation is proposed in this paper. With this simple normal velocity reconstruction procedure, all of the commonly used flux solvers (such as Godunov, Roe, Harten-Lax-van Leer-Contact, Advection-Upstream-Splitting-Method, etc.) turn out to be carbuncle-free and shock-stable. The normal velocity reconstruction procedure is done by a linear reconstruction of the cell interface normal velocity with the transverse neighbor cells, in consideration of the information transport in the transverse direction, which is neglected in the conventional finite wvolume/difference method. Some typical cases are performed to show that, when the normal velocity reconstruction procedure is used, various schemes (e.g., Roe) become carbuncle-free and shock-stable. In addition, the normal velocity reconstruction procedure has no influence on the contact-preserving property of the original flux solvers, and it adds very little computational cost. The mechanism of the normal velocity reconstruction procedure also is analyzed by a matrix-based stability method, and the results indicate that the normal velocity reconstruction procedure effectively reduces the system's positive eigenvalues that are leading to shock instability.
A shock-stable numerical scheme accurate for contact discontinuities: Applications to 3D compressible flows

Lijun Hu,Xiaohui Wang

DOI: https://doi.org/10.1016/j.cnsns.2023.107602

IF: 4.186

2023-11-04

Communications in Nonlinear Science and Numerical Simulation

Abstract:The low-dissipation Roe scheme is a popular shock capturing scheme but the shock instability, such as the infamous carbuncle phenomenon, seriously affects its application in high Mach number flows. Although considerable research has been undertaken in two dimensions, mechanism analysis of the three-dimensional shock instability is relatively scarce. Numerical simulations of the Sedov blast wave indicate that the shock instability is more prominent in three dimensions, so it has theoretical and practical significance to analyze and heal the numerical instability. The mechanism for three-dimensional shock instability has been thoroughly studied by means of the linearized stability analysis and the dissipation-controlling approach which introduces the required transverse dissipation in the numerical shock layer and is adopted to enhance the Roe scheme's shock stability. Furthermore, the unphysical expansion shock resulting from the violation of entropy condition is eliminated by a simple modification of numerical signal speed. For the improvement of the resolution for contact discontinuities and shear waves, an algebraic method which combines the THINC reconstruction scheme and the BVD algorithm is employed to minimize the density difference in the numerical diffusion term. A series of benchmark numerical experiments fully demonstrate that the proposed scheme is endowed with excellent robustness against the shock instability and high accuracy for contact discontinuities and shear waves.

mathematics, applied, interdisciplinary applications,mechanics,physics, mathematical, fluids & plasmas
Development of accurate and robust genuinely two-dimensional HLL-type Riemann solver for compressible flows

Lijun Hu,Li Yuan,Kunlei Zhao

DOI: https://doi.org/10.1016/j.compfluid.2020.104719

2020-12-01

Abstract:<p>When using the conventional direction splitting method to calculate multidimensional high speed gasdynamical flows, Riemann solvers capable of resolving contact surface and shear wave accurately will suffer from different forms of shock instability, such as the notorious carbuncle phenomenon. The stability analysis shows that the lack of dissipation in the direction transverse to the shock front leads to the shock instability of low-dissipation HLLEM solver. To overcome this defect, an accurate and carbuncle-free genuinely two-dimensional HLL-type Riemann solver is proposed. Using Zha-Bilgen splitting method, the flux vector of two-dimensional Euler equations is split into the convective flux and pressure flux. An algorithm similar to AUSM+ scheme is adopted to calculate the convective flux and the pressure flux is calculated by the low-dissipation HLLEM scheme. Following Balsara's idea, the genuinely two-dimensional properties of the new solver are achieved by solving the two-dimensional Riemann problem that considers transversal features of the flow at each vertex of the cell interface. Numerical results of benchmark tests demonstrate that the new solver has higher resolution and better robustness than the conventional HLLEM solver implemented in dimension by dimension.</p>

mechanics,computer science, interdisciplinary applications
Mechanism Study of Shock Instability in Riemann-Solver-Based Shock-Capturing Scheme

Zhiqiang Chen,Xudong Huang,Yu-Xin Ren,Zhifeng Xie,Ming Zhou

DOI: https://doi.org/10.2514/1.j056881

IF: 2.624

2018-01-01

AIAA Journal

Abstract:A novel shock instability analysis model named cell-shock-cell system is presented. Shock-related nonunique numerical solutions of the Euler equation are presented in this system, which indicates that this simple system keeps the essential characteristics of numerical shock instability. The stability analysis of this system was conducted, and some important results were obtained. Carbuncle phenomena were shown to originate from multidimensional instability of the upstream shock cell. The numerical shock instability essentially starts from the grid face perpendicular to the shock rather than the grid parallel to it. The analytical eigenvalue analysis of the linearized equations showed that inappropriate pressure dissipation of Riemann solver at the vertical transverse face of a shock was responsible for the carbuncle phenomena and shock instability. Based on the mechanism that was identified, a fundamental cure principle was proposed. These new findings can serve as a basis for understanding the numerical shock instabilities presented in one-dimensional Riemann-solver-based numerical schemes as well as in designing shock-stable schemes.
An efficient, adaptive solver for accurate simulation of multicomponent shock-interface problems for thermally perfect species

Yuqi Wang,Ralf Deiterding,Jianhan Liang

2024-11-20

Abstract:A second-order-accurate finite volume method, hybridized by blending an extended double-flux algorithm and a traditionally conservative scheme, is developed. In this scheme, hybrid convective fluxes as well as hybrid interpolation techniques are designed to ensure stability and accuracy in the presence of both material interfaces and shocks. Two computationally efficient approaches, extended from the original double-flux model, are presented to eliminate the well-known "pressure oscillation" phenomenon at material interfaces observed with the traditional conservative scheme. Numerous verification simulations confirm that the method is capable of handling multi-dimensional shock-interface problems reliably and efficiently, even in the presence of viscous and reactive terms.

Computational Physics
A discussion on numerical shock stability of unstructured finite volume method: Riemann solvers and limiters

Fan Zhang,Zhichao Yuan,Jun Liu

2023-11-29

Abstract:Numerical shock instability is a complexity which may occur in supersonic simulations. Riemann solver is usually the crucial factor that affects both the computation accuracy and numerical shock stability. In this paper, several classical Riemann solvers are discussed, and the intrinsic mechanism of shock instability is especially concerned. It can be found that the momentum perturbation traversing shock wave is a major reason that invokes instability. Furthermore, slope limiters used to depress oscillation across shock wave is also a key factor for computation stability. Several slope limiters can cause significant numerical errors near shock waves, and make the computation fail to converge. Extra dissipation of Riemann solvers and slope limiters can be helpful to eliminate instability, but reduces the computation accuracy. Therefore, to properly introduce numerical dissipation is critical for numerical computations. Here, pressure based shock indicator is used to show the position of shock wave and tunes the numerical dissipation. Overall, the presented methods are showing satisfactory results in both the accuracy and stability.

Numerical Analysis,Fluid Dynamics
A Robust High-Order Residual Distribution Type Scheme for Steady Euler Equations on Unstructured Grids

Guanghui Hu,Ruo Li,Tao Tang

DOI: https://doi.org/10.1016/j.jcp.2009.11.002

IF: 4.645

2010-01-01

Journal of Computational Physics

Abstract:A robust high-order algorithm is proposed to solve steady Euler equations on unstructured grids. The main ingredients of the algorithm include a standard Newton method as the outer iterative scheme and a linear multigrid method as the inner iterative scheme with the block lower-upper symmetric Gauss–Seidel (LU-SGS) iteration as its smoother. The Jacobian matrix of the Newton-iteration is regularized by the local residual, instead of using the commonly adopted time-stepping relaxation technique based on the local CFL number. The local Jacobian matrix of the numerical fluxes are computed using the numerical differentiation, which can significantly simplify the implementations by comparing with the manually derived approximate derivatives. The approximate polynomial of solution on each cell is reconstructed by using the values on centroid of the cell, and limited by the WENO hierarchical limiting strategy proposed by Xu et al. [Z.L. Xu, Y.J. Liu, C.W. Shu, Hierarchical reconstruction for discontinuous galerkin methods on unstructured grids with a WENO-type linear reconstruction and partial neighboring cells, Journal of Computational Physics 228 (2009) 2194–2212]. It is found that the proposed algorithm is insensitive to the parameters used. More precisely, in our computations, only one set of the parameters (namely, the proportional constant α for the local residual, the relaxation parameter τ in the Newton-iteration, the weight μ in the WENO scheme and the number of smoothing steps in the multigrid solver) is employed for various geometrical configurations and free-stream configurations. The high-order and robustness of our algorithm are illustrated by considering two-dimensional airfoil problems with different geometrical configurations and different free-stream configurations.
Shock-Stable Roe Scheme Combining Entropy Fix and Rotated Riemann Solver

Xue-Song Li,Xiao-Dong Ren,Chun-Wei Gu,Yu-Hong Li

DOI: https://doi.org/10.2514/1.j058549

IF: 2.624

2019-01-01

AIAA Journal

Abstract:The Roe scheme is an important shock-capturing scheme that is known for its good performance. However, the classical Roe scheme suffers from disastrous problems, such as shock instability for hypersonic flows. The entropy fix and the rotated Riemann solver are common methods used to cure the shock instability. However, these methods may introduce a large and unnecessary numerical dissipation. To cure the shock instability with a minimally increasing numerical dissipation, this study investigates the mechanisms of entropy fix and rotated Roe scheme and combines their complementary advantages. Thus, this paper proposes the Roe-ER scheme, which uses the rotated method to identify where the entropy fix is necessary. Therefore, the possible large dissipation of the rotated Roe scheme is replaced by the relatively small value of entropy fix, and areas that adopt the entropy fix are limited in regions where the shock instability may occur. Numerical cases validate the Roe-ER scheme, which is shock stable, features a low additional numerical dissipation, and can use high-order reconstruction to further reduce areas that activate the entropy fix.
Shock Capturing Scheme Based on Two- Shock Approximation to Riemann Problem

任玉新,徐宏龙

DOI: https://doi.org/10.3321/j.issn:1000-0054.2001.11.005

2001-01-01

Abstract:Riemann problem is a building block of modern high resolution shock capturing schemes. The exact solution of Riemann problem can not be written in close form. Certain iteration procedures are needed to get the explicit solution. Therefore, numerical schemes based on exact Riemann solver are computationally very expensive. In order to improve the computational efficiency, a shock capturing scheme to the original Riemann problem is proposed based on the two shock approximation procedure. In the scheme, all waves in the structure of Riemann problem solution are treated as shock waves. This method keeps non linear feature for Riemann problem, and the entropy does not decrease. Numerical tests indicate that the proposed scheme is efficient and of high resolution to shocks and contact discontinuities.
Numerical experiments using one high-resolution scheme for unsteady inviscid compressible flows

Li XU,Ai-ming YANG,Peng GUAN,Pei-fen WENG

DOI: https://doi.org/10.3969/j.issn.0258-1825.2009.05.018

2009-01-01

Abstract:One high-resolution scheme is developed to compute unsteady, inviscid, compressible flows. In this paper,we are concerned about AUSMPW+Riemann solver with the weighting function f based on pressure by considering monotonicity and compare it with Roe Riemann solver. For better accuracy, fifth-order WENO schemes are adopted to interpolate higher order left and right states across a cell interface. Third-order TVD Runge-Kutta method is used to advance the solution in time. The performance of the schemes has been assessed in four typical unsteady inviscid compressible flow problems, moving contact discontinuity problem, Sod's problem, regular shock reflection, and double Math reflection of a strong shock. AUSMPW+scheme with fifth-order WENO schemes interpolation proves to have the great capabilities to capture shock and contact discontinuities with high resolution and be more efficient than Roe scheme.
A Robust Riemann Solver for Multiple Hydro-Elastoplastic Solid Mediums

Ruo Li,Yanli Wang,Chengbao Yao

DOI: https://doi.org/10.48550/arXiv.1902.05145

2019-02-03

Abstract:We propose a robust approximate solver for the hydro-elastoplastic solid material, a general constitutive law extensively applied in explosion and high speed impact dynamics, and provide a natural transformation between the fluid and solid in the case of phase transitions. The hydrostatic components of the solid is described by a family of general Mie-Grüneisen equation of state (EOS), while the deviatoric component includes the elastic phase, linearly hardened plastic phase and fluid phase. The approximate solver provides the interface stress and normal velocity by an iterative method. The well-posedness and convergence of our solver are proved with mild assumptions on the equations of state. The proposed solver is applied in computing the numerical flux at the phase interface for our compressible multi-medium flow simulation on Eulerian girds. Several numerical examples, including Riemann problems, shock-bubble interactions, implosions and high speed impact applications, are presented to validate the approximate solver.

Numerical Analysis
A shock-capturing scheme with a novel limiter for compressible flows solved by smoothed particle hydrodynamics

Zi-Fei Meng,A-Man Zhang,Ping-Ping Wang,Fu-Ren Ming

DOI: https://doi.org/10.1016/j.cma.2021.114082

IF: 6.588

2021-12-01

Computer Methods in Applied Mechanics and Engineering

Abstract:There are many compressible flows involving strong discontinuities in nature and engineering applications. Their behaviors are complicated by the existence of shock waves, rarefaction waves and contact discontinuities. To investigate these flows, a shock-capturing scheme based on smoothed particle hydrodynamics (SPH) is proposed. In this scheme, Roe’s approximate Riemann solver cooperated with a novel limiter is embedded into the SPH governing equations to capture shocks. This limiter is simple and effective to control numerical dissipations, and the use of this limiter eliminates the tunable artificial viscosity as required by the conventional SPH. Additionally, to restore the accuracy limited by heterogeneous particle distribution that is usually encountered in the simulation of strongly-compressible flows, the gradient operator in the continuity equation is corrected by a renormalization procedure. Last but not least, for initial particle distribution, unlike the equal mass particle distribution commonly used in compressible SPH simulations, the equal spacing particle distribution is adopted in this scheme, which makes it much easier to model multidimensional problems. The present scheme has been verified by a set of benchmark tests involving contact discontinuities, extreme shock waves and strong rarefaction waves, some of which are firstly simulated by the SPH method.

mechanics,engineering, multidisciplinary,mathematics, interdisciplinary applications
Improved SIP scheme for numerical solutions of 3-D Euler and N-S equations and its application

Qiusheng Liu,Baoguo Wang,Mengyu Shen

1996-01-01

Abstract:The scalar scheme of Stone's strongly implicit procedure (SIP) was extended to general systems of multidimensional conservation laws. The systems contain 3-D Euler equations or 3-D Navier-Stokes equations. This study focuses on the simulation of fluid flows with strong shocks in three dimensions, and presents an efficient and high-resolution implicit hybrid scheme combining improved SIP with the NND scheme, which is proved suitable for shock-capturing calculation in gas dynamics. Numerical experiments suggest that the improved SIP scheme not only be able to have a fairly rapid convergence rate, but also can generate a highly resolved approximation to the steady-state solution. The scheme leads to an effective nonoscillatory shock capturing method for steady flow.
Robust and accurate Roe-type Riemann solver with compact stencil: Rotated-RoeM scheme

Seongyu Choi,Donguk Kim,Jaehyong Park,Jin Seok Park

DOI: https://doi.org/10.1016/j.jcp.2024.112913

IF: 4.645

2024-03-10

Journal of Computational Physics

Abstract:In this study, we aim to develop a robust and accurate Riemann solver to resolve strong shock waves with compact stencil information by applying the rotated hybrid concept to the RoeM scheme. The original RoeM scheme offers good robustness from the perspective of avoiding shock instability, also called the carbuncle phenomenon, in hypersonic flows by introducing Mach-number-based functions. These functions rely on multi-dimensional shock discontinuity sensing terms to capture grid-aligned shock waves, thus necessitating sharing of values of cells adjacent to each vertex. Consequently, constructing a data structure that can handle such a wide stencil is necessary, particularly for implementations on unstructured grids. Moreover, communication of this stencil information on parallel computation can increase computational costs. To overcome these complexities, the rotated hybrid concept is applied to the RoeM scheme, to replace the multi-dimensional terms while still resolving the grid-aligned shock instability. In addition, a Mach-number-based weighting function is introduced to distinguish between the flow around the boundary layer and the onset of shock instability. The damping characteristic of the proposed scheme is compared to that of the original RoeM scheme by conducting a linear perturbation analysis. The findings reveal that the proposed scheme effectively controls spurious oscillations across strong shock waves. Numerous test cases are presented to demonstrate the robustness and accuracy of the proposed approach. Notably, the proposed scheme successfully resolves thermal boundary layers and accurately predicts the heat transfer rate around stagnation points in hypersonic flows.

physics, mathematical,computer science, interdisciplinary applications
An Efficient Modified WENO Scheme Based on the Identification of Inflection Points

Delin Chai,Guang Xi,Zhongguo Sun,Zhiheng Wang,Zhu Huang

DOI: https://doi.org/10.1016/j.compfluid.2018.04.036

IF: 3.077

2018-01-01

Computers & Fluids

Abstract:In this paper, a new modified fifth-order weighted essentially non-oscillatory (WENO) scheme is proposed by introducing a switching principle based on the identification of inflection points. The switching principle adjusts the scheme between the linear upwind scheme and the WENO-M scheme automatically. If there is one or more inflection points located in the large reconstruction stencil, the WENO-M procedure is applied to reconstruct the numerical flux, instead of the classical WENO scheme; otherwise, the linear upwind scheme is used. The new scheme is more efficient due to the identification of the inflection points, and has higher numerical accuracy and resolution due to the usage of WENO-M procedure. A set of benchmark problems are conducted and analyzed to verify the performance enhancement. The extreme shock tube problem of Leblanc is computed to demonstrate the robustness of the new scheme.

A robustness-enhanced method for Riemann solver

An Arbitrary Lagrangian–Eulerian Method Based on the Adaptive Riemann Solvers for General Equations of State

An Improved Meshfree Scheme Based on Radial Basis Functions for Solving Incompressible Navier–Stokes Equations

An Optimized Dispersion-Relation-Preserving Combined Compact Difference Scheme to Solve Advection Equations

A robust shock-capturing scheme based on rotated Riemann solvers

Mechanism-Derived Shock Instability Elimination for Riemann-Solver-Based Shock-Capturing Scheme

General Procedure for Riemann Solver to Eliminate Carbuncle and Shock Instability

A shock-stable numerical scheme accurate for contact discontinuities: Applications to 3D compressible flows

Development of accurate and robust genuinely two-dimensional HLL-type Riemann solver for compressible flows

Mechanism Study of Shock Instability in Riemann-Solver-Based Shock-Capturing Scheme

An efficient, adaptive solver for accurate simulation of multicomponent shock-interface problems for thermally perfect species

A discussion on numerical shock stability of unstructured finite volume method: Riemann solvers and limiters

A Robust High-Order Residual Distribution Type Scheme for Steady Euler Equations on Unstructured Grids

Shock-Stable Roe Scheme Combining Entropy Fix and Rotated Riemann Solver

Shock Capturing Scheme Based on Two- Shock Approximation to Riemann Problem

Numerical experiments using one high-resolution scheme for unsteady inviscid compressible flows

A Robust Riemann Solver for Multiple Hydro-Elastoplastic Solid Mediums

A shock-capturing scheme with a novel limiter for compressible flows solved by smoothed particle hydrodynamics

Improved SIP scheme for numerical solutions of 3-D Euler and N-S equations and its application

Robust and accurate Roe-type Riemann solver with compact stencil: Rotated-RoeM scheme

An Efficient Modified WENO Scheme Based on the Identification of Inflection Points