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.

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.

Transversal effects of high order numerical schemes for compressible fluid flows

Finite volume schemes for the two-dimensional(2D) wave system are taken to demonstrate the role of the genuine dimensionality of Lax-Wendroff flow solvers for compressible fluid flows. When the finite volume schemes are applied, the transversal variation relative to the computational cell interfaces is neglected, and only the normal numerical flux is used, thanks to the Gauss-Green formula. In order to offset such defects, the Lax-Wendroff flow solvers or the generalized Riemann problem(GRP) solvers are adopted by substituting the time evolution of flows into the spatial variation. The numerical results show that even with the same convergence rate, the error by the GRP2D solver is almost one ninth of that by the multistage Runge-Kutta(RK) method.

计算流体力学的时空观:模型的时空关联性及算法的时空耦合性

流体力学中波的有限传播、粒子的碰撞、各种力之间相互作用,无不体现时空关联效应。本文从计算方法的视角探讨计算流体力学的时空观,即流体力学模型的时空关联性和计算方法的时空耦合性。从流体力学微团法建模出发,明确模型时空关联性的涵义,建立有限体积格式的基本原理,阐述算法时空耦合的必要性,实现流体力学基本控制方程物理建模与有限体积格式数学原理的统一。在实践中,给出时空耦合高精度数值方法设计思路,利用算例比较它与时空解耦方法的差别。期望通过时空观的建立,对未来计算流体力学的算法研究提供帮助。