An Approximate Riemann Solver for Advection-Diffusion Based on the Generalized Riemann Problem

We construct an approximate Riemann solver for scalar advection-diffusion equations with piecewise polynomial initial data.The objective is to handle advection and diffusion simultaneously to reduce the inherent numerical diffusion produced by the usual advection flux calculations.The approximate solution is based on the weak formulation of the Riemann problem and is solved within a space-time discontinuous Galerkin approach with two subregions.The novel generalized Riemann solver produces piecewise polynomial solutions of the Riemann problem.In conjunction with a recovery polynomial,the Riemann solver is then applied to define the numerical flux within a finite volume method.Numerical results for a piecewise linear and a piecewise parabolic approximation are shown.These results indicate a reduction in numerical dissipation compared with the conventional separated flux calculation of advection and diffusion.Also,it is shown that using the proposed solver only in the vicinity of discontinuities gives way to an accurate and efficient finite volume scheme.

Generalized Riemann problem for gas dynamic combustion

The generalized Riemann problem for gas dynamic combustion in a neighborhood of the origin t 0 in the （x, t） plane is considered. Under the modified entropy conditions, the unique solutions are constructed, which are the limits of the selfsimilar Zeldovich-von Neumann-Dring （ZND） combustion model. The results show that, for some cases, there are intrinsical differences between the structures of the perturbed Riemann solutions and the corresponding Riemann solutions. Especially, a strong detonation in the corresponding Riemann solution may be transformed into a weak deflagration coalescing with the pre-compression shock wave after perturbation. Moreover, in some cases, although no combustion wave exists in the corresponding Riemann solution, the combustion wave may occur after perturbation, which shows the instability of the unburnt gases.

ADER Methods for Hyperbolic Equations with a Time-Reconstruction Solver for the Generalized Riemann Problem: the Scalar Case

The ADER approach to solve hyperbolic equations to very high order of accuracy has seen explosive developments in the last few years,including both methodological aspects as well as very ambitious applications.In spite of methodological progress,the issues of efficiency and ease of implementation of the solution of the associated generalized Riemann problem(GRP)remain the centre of attention in the ADER approach.In the original formulation of ADER schemes,the proposed solution procedure for the GRP was based on(i)Taylor series expansion of the solution in time right at the element interface,(ii)subsequent application of the Cauchy-Kowalewskaya procedure to convert time derivatives to functionals of space derivatives,and(iii)solution of classical Riemann problems for high-order spatial derivatives to complete the Taylor series expansion.For realistic problems the Cauchy-Kowalewskaya procedure requires the use of symbolic manipulators and being rather cumbersome its replacement or simplification is highly desirable.In this paper we propose a new class of solvers for the GRP that avoid the Cauchy-Kowalewskaya procedure and result in simpler ADER schemes.This is achieved by exploiting the history of the numerical solution that makes it possible to devise a time-reconstruction procedure at the element interface.Still relying on a time Taylor series expansion of the solution at the interface,the time derivatives are then easily calculated from the time-reconstruction polynomial.The resulting schemes are called ADER-TR.A thorough study of the linear stability properties of the linear version of the schemes is carried out using the von Neumann method,thus deducing linear stability regions.Also,via careful numerical experiments,we deduce stability regions for the corresponding non-linear schemes.Numerical examples using the present simplified schemes of fifth and seventh order of accuracy in space and time show that these compare favourably with conventional ADER methods.This paper is restricted to the one-dimensional scalar case with source term,but preliminary results for the one-dimensional Euler equations indicate that the time-reconstruction approach offers significant advantages not only in terms of ease of implementation but also in terms of efficiency for the high-order range schemes.

A NOTE ON THE GENERALIZED RIEMANN PROBLEM~*

1. Introduction The Riemann problem for the quasilinear hyperbolic system of conservation laws

High-order accurate solutions of generalized Riemann problems of nonlinear hyperbolic balance laws

We provide a systematic study for the generalized Riemann problem(GRP)of the nonlinear hyperbolic balance law,which is critically concerned with the construction of the spatial-temporally coupled high-order Godunov-type scheme.The full analytical GRP solvers up to the third-order accuracy and also a collection of properties of the GRP solution are derived by resolving the elementary waves.The resolution of the rarefaction wave is a crucial point,which relies on the use of the generalized characteristic coordinate(GCC)to analyze the solution at the singularity.From the analysis on the GCC,we derive for the general nonlinear system the evolutionary equations for the derivatives of generalized Riemann invariants.For the nonsonic case,the full set of spatial and temporal derivatives of the GRP solution at the singularity are obtained,whereas for the sonic case the limiting directional derivatives inside the rarefaction wave are derived.In addition,the acoustic approximation of the analytical GRP solver is deduced by estimating the error it introduces.It is shown that the computationally more efficient Toro-Titarev solver can be the approximation of the analytical solver under the suitable condition.Hence this work also provides a theoretical basis of the approximate GRP solver.The theoretical results are illustrated via the examples of the Burgers equation,the shallow water equations and a system for compressible flows under gravity acceleration.Numerical results demonstrate the accuracy of the GRP solvers for both weak and strong discontinuity cases.

The Generalized Riemann Problem for a Chromatography System

The generalized Riemann problem for a model of chromatography system in a conservative form is considered.The unique local solution in the class of piecewise C1functions in a neighborhood of the origin is obtained.The structures of the solutions are similar to the corresponding Riemann problem,which means the Riemann solutions are stable.

Generalized δ-entropy condition to Riemann solution for Chaplygin gas in traffic flow

The Riemann problem for the Aw-Rascle model in the traffic flow with the Delta initial data for the Chaplygin gas is studied. The solutions are constructed globally under the generalized Rankine-Hugoniot relations, the δ-entropy conditions, and the generalized δ-entropy conditions. A new Delta wave, which is called a primary Delta wave, is defined in some solutions. The primary Delta wave satisfies the generalized Rankine- Hugoniot relations and the generalized δ-entropy conditions. It generates initially from the Delta initial data, which either evaluates a Delta wave, whose weight becomes stronger and stronger, or disappears at a finite time.

GLOBAL INSTABILITY OF MULTI-DIMENSIONAL PLANE SHOCKS FOR ISOTHERMAL FLOW

In this paper,we are concerned with the long time behavior of the piecewise smooth solutions to the generalized Riemann problem governed by the compressible isothermal Euler equations in two and three dimensions.A non-existence result is established for the fan-shaped wave structure solution,including two shocks and one contact discontinuity which is a perturbation of plane waves.Therefore,unlike in the one-dimensional case,the multi-dimensional plane shocks are not stable globally.Moreover,a sharp lifespan estimate is established which is the same as the lifespan estimate for the nonlinear wave equations in both two and three space dimensions.

Application of ADER Scheme in MHD Simulation

The Arbitrary accuracy Derivatives Riemann problem method(ADER) scheme is a new high order numerical scheme based on the concept of finite volume integration,and it is very easy to be extended up to any order of space and time accuracy by using a Taylor time expansion at the cell interface position.So far the approach has been applied successfully to flow mechanics problems.Our objective here is to carry out the extension of multidimensional ADER schemes to multidimensional MHD systems of conservation laws by calculating several MHD problems in one and two dimensions: (ⅰ) Brio-Wu shock tube problem,(ⅱ) Dai-Woodward shock tube problem,(ⅲ) Orszag-Tang MHD vortex problem.The numerical results prove that the ADER scheme possesses the ability to solve MHD problem,remains high order accuracy both in space and time,keeps precise in capturing the shock.Meanwhile,the compared tests show that the ADER scheme can restrain the oscillation and obtain the high order non-oscillatory result.

Lie Symmetry Analysis and Conservation Laws of a Generalized Time Fractional Foam Drainage Equation

In this paper, a generalized time fractional nonlinear foam drainage equation is investigated by means of the Lie group analysis method. Based on the Riemann–Liouville derivative, the Lie point symmetries and symmetry reductions of the equation are derived, respectively. Furthermore, conservation laws with two kinds of independent variables of the equation are performed by making use of the nonlinear self-adjointness method.

Local Smooth Solutions to the 3-Dimensional Isentropic Relativistic Euler equations

The authors consider the local smooth solutions to the isentropic relativistic Euler equations in(3+1)-dimensional space-time for both non-vacuum and vacuum cases.The local existence is proved by symmetrizing the system and applying the FriedrichsLax-Kato theory of symmetric hyperbolic systems.For the non-vacuum case,according to Godunov,firstly a strictly convex entropy function is solved out,then a suitable symmetrizer to symmetrize the system is constructed.For the vacuum case,since the coefficient matrix blows-up near the vacuum,the authors use another symmetrization which is based on the generalized Riemann invariants and the normalized velocity.

Selberg's Normal Density Theorem for Automorphic L-Functions for GL_m

Let π be an irreducible unitary cuspidal representation of GLm（AQ） with m ≥ 2, and L（s, Tr） the L-function attached to π. Under the Generalized Riemann Hypothesis for L（s,π）, we estimate the normal density of primes in short intervals for the automorphic L-function L（s, π）. Our result generalizes the corresponding theorem of Selberg for the Riemann zeta-function.

Waves Induced by the Appearance of Single-Point Heat Source in Constant Flow

Heating or cooling one-dimensional inviscid compressible flow can be modeled by the Euler equations with energy sources.A tricky situation is the sudden appearance of a single-point energy source term.This source is discontinuous in both the time and space directions,and results in multiple discontinuous waves in the solution.We establish a mathematical model of the generalized Riemann problem of the Euler equations with source term.Based on the double CRPs coupling method proposed by the authors,we determine the wave patterns of the solution.Theoretically,we prove the existence and uniqueness of solutions to both"heat removal"problem and"heat addition"problem.Our results provide a theoretical explanation for the effect of instantaneous addition or removal of heat on the fluid.