Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAn Efficient, Second Order Accurate, Universal Generalized Riemann Problem Solver Based on the HLLI Riemann Solver
91
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The Riemann problem, and the associated generalized Riemann problem, are increasingly seen as the important building blocks for modern higher order Godunov-type schemes. In the past, building a generalized Riemann problem solver was seen as an intricately mathematical task for complicated physical or engineering problems because the associated Riemann problem is different for each hyperbolic system of interest. This paper changes that situation. The HLLI Riemann solver is a recently-proposed Riemann solver that is universal in that it is applicable to any hyperbolic system, whether in conservation form or with non-conservative products. The HLLI Riemann solver is also complete in the sense that if it is given a complete set of eigenvectors, it represents all waves with minimal dissipation. It is, therefore, very attractive to build a generalized Riemann problem solver version of the HLLI Riemann solver. This is the task that is accomplished in the present paper. We show that at second order, the generalized Riemann problem version of the HLLI Riemann solver is easy to design. Our GRP solver is also complete and universal because it inherits those good properties from original HLLI Riemann solver. We also show how our GRP solver can be adapted to the solution of hyperbolic systems with stiff source terms. Our generalized HLLI Riemann solver is easy to implement and performs robustly and well over a range of test problems. All implementation-related details are presented. Results from several stringent test problems are shown. These test problems are drawn from many different hyperbolic systems, and include hyperbolic systems in conservation form; with non-conservative products; and with stiff source terms. The present generalized Riemann problem solver performs well on all of them.
rate research
Read More
Recently developed concept of dissipative measure-valued solution for compressible flows is a suitable tool to describe oscillations and singularities possibly developed in solutions of multidimensional Euler equations. In this paper we study the convergence of the first-order finite volume method based on the exact Riemann solver for the complete compressible Euler equations. Specifically, we derive entropy inequality and prove the consistency of numerical method. Passing to the limit, we show the weak and strong convergence of numerical solutions and identify their limit. The numerical results presented for the spiral, Kelvin-Helmholtz and the Richtmyer-Meshkov problem are consistent with our theoretical analysis.
A Riemann problem with prescribed initial conditions will produce one of three possible wave patterns corresponding to the propagation of the different discontinuities that will be produced once the system is allowed to relax. In general, when solving the Riemann problem numerically, the determination of the specific wave pattern produced is obtained through some initial guess which can be successively discarded or improved. We here discuss a new procedure, suitable for implementation in an exact Riemann solver in one dimension, which removes the initial ambiguity in the wave pattern. In particular we focus our attention on the relativistic velocity jump between the two initial states and use this to determine, through some analytic conditions, the wave pattern produced by the decay of the initial discontinuity. The exact Riemann problem is then solved by means of calculating the root of a nonlinear equation. Interestingly, in the case of two rarefaction waves, this root can even be found analytically. Our procedure is straightforward to implement numerically and improves the efficiency of numerical codes based on exact Riemann solvers.
We extend our approach for the exact solution of the Riemann problem in relativistic hydrodynamics to the case in which the fluid velocity has components tangential to the initial discontinuity. As in one-dimensional flows, we here show that the wave-pattern produced in a multidimensional relativistic Riemann problem can be predicted entirely by examining the initial conditions. Our method is logically very simple and allows for a numerical implementation of an exact Riemann solver which is both straightforward and computationally efficient. The simplicity of the approach is also important for revealing special relativistic effects responsible for a smooth transition from one wave-pattern to another when the tangential velocities in the initial states are suitably varied. While the content of this paper is focussed on a flat spacetime, the local Lorentz invariance allows its use also in fully general relativistic calculations.
We propose an efficient, accurate and robust implicit solver for the incompressible Navier-Stokes equations, based on a DG spatial discretization and on the TR-BDF2 method for time discretization. The effectiveness of the method is demonstrated in a number of classical benchmarks, which highlight its superior efficiency with respect to other widely used implicit approaches. The parallel implementation of the proposed method in the framework of the deal.II software package allows for accurate and efficient adaptive simulations in complex geometries, which makes the proposed solver attractive for large scale industrial applications.
In [Z. Hu, R. Li, and Z. Qiao. Acceleration for microflow simulations of high-order moment models by using lower-order model correction. J. Comput. Phys., 327:225-244, 2016], it has been successfully demonstrated that using lower-order moment model correction is a promising idea to accelerate the steady-state computation of high-order moment models of the Boltzmann equation. To develop the existing solver, the following aspects are studied in this paper. First, the finite volume method with linear reconstruction is employed for high-resolution spatial discretization so that the degrees of freedom in spatial space could be reduced remarkably without loss of accuracy. Second, by introducing an appropriate parameter $tau$ in the correction step, it is found that the performance of the solver can be improved significantly, i.e., more levels would be involved in the solver, which further accelerates the convergence of the method. Third, Heuns method is employed as the smoother in each level to enhance the robustness of the solver. Numerical experiments in microflows are carried out to demonstrate the efficiency and to investigate the behavior of the new solver. In addition, several order reduction strategies for the choice of the order sequence of the solver are tested, and the strategy $m_{l-1} = lceil m_{l} / 2 rceil$ is found to be most efficient.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
