Central ADER schemes for hyperbolic conservation laws

In this paper we first briefly review the very high order ADER methods for solving hyperbolic conservation laws. ADER methods use high order polynomial reconstruction of the solution and upwind fluxes as the building block. They use a first order upwind Godunov and the upwind second order weighted average (WAF) fluxes. As well known the upwind methods are more accurate than central schemes. However, the superior accuracy of the ADER upwind schemes comes at a cost, one must solve exactly or approximately the Riemann problems (RP). Conventional Riemann solvers are usually complex and are not available for many hyperbolic problems of practical interest. In this paper we propose to use two central fluxes, instead of upwind fluxes, as the building block in ADER scheme. These are the monotone first order Lax-Friedrich (LXF) and the third order TVD flux. The resulting schemes are called central ADER schemes. Accuracy of the new schemes is established. Numerical implementations of the new schemes are carried out on the scalar conservation laws with a linear flux, nonlinear convex flux and non-convex flux. The results demonstrate that the proposed scheme, with LXF flux, is comparable to those using first and second order upwind fluxes while the scheme, with third order TVD flux, is superior to those using upwind fluxes. When compared with the state of art ADER schemes, our central ADER schemes are faster, more accurate, Riemann solver free, very simple to implement and need less computer memory. A way to extend these schemes to general systems of nonlinear hyperbolic conservation laws in one and two dimensions is presented.     

One-Sided Exact Boundary Null Controllability of Entropy Solutions to a Class of Hyperbolic Systems of Conservation Laws with Characteristics with Constant Multiplicity

This paper proves the local exact one-sided boundary null controllability of entropy solutions to a class of hyperbolic systems of conservation laws with characteristics with constant multiplicity. This generalizes the results in [Li, T. and Yu, L., One-sided exact boundary null controllability of entropy solutions to a class of hyperbolic systems of conservation laws, To appear in Journal de Mathématiques Pures et Appliquées, 2016.] for a class of strictly hyperbolic systems of conservation laws.     

The numerical interface coupling of nonlinear hyperbolic systems of conservation laws: I. The scalar case

Summary We study the theoretical and numerical coupling of two general hyperbolic conservation laws. The coupling preserves in a weak sense the continuity of the solution at the interface without imposing the overall conservativity of the coupled model. In order to analyze the convergence of the coupled numerical scheme, we first revisit the approximation of the boundary value problems. We then prove the convergence and characterize the limit solution of the coupled schemes in a few simple but significative coupling situations. The general coupling problem is analyzed for Riemann initial data and illustrated by numerical simulations. Résumé. Nous nous intéressons à une nouvelle forme de couplage de deux systèmes hyperboliques de lois de conservation. Ce couplage assure de façon faible la continuité de la solution à linterface sans imposer la conservativité du modèle couplé. Pour étudier la convergence du schéma dapproximation numérique, nous commençons par reprendre les résultats concernant lapproximation du problème aux limites. Nous démontrons ensuite la convergence du schéma couplé dans un certain nombre de cas intéressants. Le cas général du couplage est étudié et illustré numériquement pour une donnée initiale de Riemann.Mathematics Subject Classification (2000):65M12, 65M30, 76M12, 35L50, 35L6524 December 2001     

Finite-difference methods for one-dimensional hyperbolic conservation laws

This article contains a survey of some important finite-difference methods for one-dimensional hyperbolic conservation laws. Weak solutions of hyperbolic conservation laws are introduced and the concept of entropy stability is discussed. Furthermore, the Riemann problem for hyperbolic conservation laws is solved. An introduction to finite-difference methods is given for which important concepts such as, e.g., conservativity, stability, and consistency are introduced. Godunov-type methods are elaborated for general systems of hyperbolic conservation laws. Finally, flux limiter methods are developed for the scalar nonlinear conservation law. © 1994 John Wiley & Sons, Inc.     

Relaxation WENO schemes for multidimensional hyperbolic systems of conservation laws

We present a class of high‐order weighted essentially nonoscillatory (WENO) reconstructions based on relaxation approximation of hyperbolic systems of conservation laws. The main advantage of combining the WENO schemes with relaxation approximation is the fact that the presented schemes avoid solution of the Riemann problems due to the relaxation approach and high‐resolution is obtained by applying the WENO approach. The emphasis is on a fifth‐order scheme and its performance for solving a wide class of systems of conservation laws. To show the effectiveness of these methods, we present numerical results for different test problems on multidimensional hyperbolic systems of conservation laws. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2007     

MUSTA: A multi-stage numerical flux

Numerical methods for conservation laws constructed in the framework of finite volume and discontinuous Galerkin finite elements require, as the building block, a monotone numerical flux. In this paper we present some preliminary results on the MUSTA approach [E.F. Toro, Multi-stage predictor–corrector fluxes for hyperbolic equations, Technical Report NI03037-NPA, Isaac Newton Institute for Mathematical Sciences, University of Cambridge, UK, 17th June, 2003] for constructing upwind numerical fluxes. The scheme may be interpreted as an un-conventional approximate Riemann solver that has simplicity and generality as its main features. When used in its first-order mode we observe that the scheme achieves the accuracy of the Godunov method used in conjunction with the exact Riemann solver, which is the reference first-order method for hyperbolic systems. At least for the scalar model hyperbolic equation, the Godunov scheme is the best of all first-order monote schemes, it has the smallest truncation error. Extensions of the scheme of this paper are realized in the framework of existing approaches. Here we present a second-order TVD (TVD for the scalar case) extension and show numerical results for the two-dimensional Euler equations on non-Cartesian geometries. The schemes find their best justification when solving very complex systems for which the solution of the Riemann problem, in the classical sense, is too complex, too costly or is simply unavailable.     

ON THE CONVERGENCE OF IMPLICIT DIFFERENCE SCHEMES FOR HYPERBOLIC CONSERVATION LAWS

1. IntroductionWe are are interested in the fOllowing Cauchy problem for scalar conservation lawswhere the initial data uo E BV(R) and the flux function f 6 C'(n).It is well known that this problem may not always have a smooth global solution even ifthe i…     

Finite volume schemes for multi-dimensional hyperbolic systems based on the use of bicharacteristics

In this paper we present recent results for the bicharacteristic based finite volume schemes, the so-called finite volume evolution Galerkin (FVEG) schemes. These methods were proposed to solve multi-dimensional hyperbolic conservation laws. They combine the usually conflicting design objectives of using the conservation form and following the characteristics, or bicharacteristics. This is realized by combining the finite volume formulation with approximate evolution operators, which use bicharacteristics of the multi-dimensional hyperbolic system. In this way all of the infinitely many directions of wave propagation are taken into account. The main goal of this paper is to present a self-contained overview on the recent results. We study the L ¹-stability of the finite volume schemes obtained by various approximations of the flux integrals. Several numerical experiments presented in the last section confirm robustness and correct multi-dimensional behaviour of the FVEG methods. This research has been supported under the VW-Stiftung grant I 76 859, by the grant No 201/03 0570 of the Grant Agency of the Czech Republic, by the Deutsche Forschungsgemeinschaft grant GK 431 and partially by the European network HYKE, funded by the EC as contract HPRN-CT-2002-00282.     

一类时空二阶精度高分辨率MmB差分格式的构造及数值试验

1．引言考虑如下二维双曲型守恒律初值问题的数值解．H．M．Wu和S．L．Yang在文山中给出了MmB差分格式的定义如下：给定（．1）M差分格式定义．若则称格式（1．2）为MmB差分格式．这里BmB表示局部MaximumandminimumBounds．由定义可知，若差分格式（1．2）可写为形式且。＼P’三0，＞。：r’一1．则格式（1．4）为MmB差分格式．j＝l文山构造了二维双曲型守恒律的二类二阶精度的MmB差分格式，使构造二维高分辨格式有了新的突破，但他们是从标量线性双曲型守恒律出发，然后把结果推广到非线性情形．本文直接从二维非线性双曲型守恒律…     

On Wavewise Entropy Inequalities for High-Resolution Schemes. I: The Semidiscrete Case

We develop a new approach, the method of wavewise entropy inequalities for the numerical analysis of hyperbolic conservation laws. The method is based on a new extremum tracking theory and Volpert's theory of BV solutions. The method yields a sharp convergence criterion which is used to prove the convergence of generalized MUSCL schemes and a class of schemes using flux limiters previously discussed in 1984 by Sweby.

     

双曲型守恒律方程的熵稳定格式的一些讨论

本文讨论双曲型守恒律方程的熵稳定格式.对于给定的熵对,格式所满足的熵条件中的数值熵通量是不唯一的.Tadmor的充分条件可以唯一地确定标量方程的熵守恒通量,但不能唯一确定方程组的熵守恒通量,却可以给出方程组的空间一阶精度的熵守恒格式.也讨论了在熵守恒通量上添加数值粘性得到的显式熵稳定格式需要满足的条件及常见的时间离散对熵守恒和熵稳定的影响.     

Multirate Runge–Kutta schemes for advection equations

Explicit time integration methods can be employed to simulate a broad spectrum of physical phenomena. The wide range of scales encountered lead to the problem that the fastest cell of the simulation dictates the global time step. Multirate time integration methods can be employed to alter the time step locally so that slower components take longer and fewer time steps, resulting in a moderate to substantial reduction of the computational cost, depending on the scenario to simulate [S. Osher, R. Sanders, Numerical approximations to nonlinear conservation laws with locally varying time and space grids, Math. Comput. 41 (1983) 321–336; H. Tang, G. Warnecke, A class of high resolution schemes for hyperbolic conservation laws and convection-diffusion equations with varying time and pace grids, SIAM J. Sci. Comput. 26 (4) (2005) 1415–1431; E. Constantinescu, A. Sandu, Multirate timestepping methods for hyperbolic conservation laws, SIAM J. Sci. Comput. 33 (3) (2007) 239–278]. In air pollution modeling the advection part is usually integrated explicitly in time, where the time step is constrained by a locally varying Courant–Friedrichs–Lewy (CFL) number. Multirate schemes are a useful tool to decouple different physical regions so that this constraint becomes a local instead of a global restriction. Therefore it is of major interest to apply multirate schemes to the advection equation. We introduce a generic recursive multirate Runge–Kutta scheme that can be easily adapted to an arbitrary number of refinement levels. It preserves the linear invariants of the system and is of third order accuracy when applied to certain explicit Runge–Kutta methods as base method.     

A new nonlinear functional for general scalar conservation laws

The approach based on the construction of some nonlinear functionals was proved to be robust in the study of the well-posedness theories of hyperbolic conservation laws, especially in one space dimensional case. In particular, a generalized entropy functional was constructed in [T.-P. Liu, T. Yang, A new entropy functional for scalar conservation laws, Comm. Pure Appl. Math. 52 (1999) 1427-1442] for the L¹ stability of weak solutions. However, this generalized functional is so far only defined for scalar equations with convex flux function. In this paper, we introduce a new nonlinear functional which is motivated by the new Glimm functional introduced in [J.-L. Hua, Z.-H. Jiang, T. Yang, A new Glimm functional and convergence rate of Glimm scheme for general systems of hyperbolic conservation laws, preprint] for general scalar conservation laws. This functional improves the one given in [H.-X. Liu, T. Yang, A nonlinear functional for general scalar hyperbolic conservation laws, J. Differential Equations 235 (2) (2007) 658-667] and it can be viewed as a better attempt for the generalized entropy functional for general equations.     

On formulation of stress-strain relations for soils at failure

Résumé Une forme générale des relations contrainte-vitesse de déformation est étudiée pour un mouvement non visqueux. Le tenseur contrainte est représenté par une fonction tensorielle isotrope du tenseur vitesse de déformation: le nombre et la forme des fonctions indépendantes décrivant le matériau sont déterminés.Le principe de la vitesse de dissipation maximale est employé pour obtenir des restrictions supplémentaires pour la relation constitutive. De plus, ce principe permet d'expliquer les contradictions provenant de la forme standard de la loi d'écoulement du potentiel plastique appliquée aux sols. Des formes particulières des relations tensorielles linéaires, dans le cas d'incompressibilité et de dissipation constante sont discutés.     

ON THE CENTRAL RELAXING SCHEME Ⅱ: SYSTEMS OF HYPERBOLIC CONSERVATION LAWS

1. IntroductionWe are interested in construction of the central reltalng sChemes for system of noIilinearhyperbolic conservation lawswith initial data U(0, x) = Uo(x), x = (x1 ? ...! xd), based on the local relaJxation approkimationof Eq.(1.1) [2, 3, 6, 8, 9, 12].To i11ustrate the basic idea of the relaalng schemes, for the sake of simplicity in the presentation, we restrict our attention to onedimensional scalar conservaioll lawsFirst, introduce a linear hyperbollc system with a stiff sourc…     

On the frame-dependence of stress and heat flux in polar fluids

Summary A special theory of polar fluids is developed with a view to providing some insight into the questions of the frame-dependence of stress and heat flux, arising out of a recent paper by Müller (1).

---

Résumé Le but de cet article est de développer une théorie spéciale des fluides à molécules polaires afin d'éclaircir deux questions soulevées dans un récent article de Muller (1), à savoir: le flux de chaleur et une possible relation de dépendance entre la tension et le système de coordonnées choisi.

     

Nonstrictly hyperbolic systems with stiff relaxation terms

In this paper, a zero factor idea is introduced to extend the convergence framework in [G.-Q. Chen, C.D. Levermore, T.-P. Liu, Hyperbolic conservation laws with stiff relaxation terms and entropy, Comm. Pure Appl. Math. 47 (1994) 787-830] for the singular limits of stiff relaxation from general 2×2 hyperbolic conservation systems to nonstrictly hyperbolic systems and an application of this framework on the so-called system of extended traffic flow is obtained.     

Adjoint IMEX-based schemes for control problems governed by hyperbolic conservation laws

Starting from relaxation schemes for hyperbolic conservation laws we derive continuous and discrete schemes for optimization problems subject to nonlinear, scalar hyperbolic conservation laws. We discuss properties of first- and second-order discrete schemes and show their relations to existing results. In particular, we introduce first and second-order relaxation and relaxed schemes for both adjoint and forward equations. We give numerical results including tracking type problems with non-smooth desired states.     

Construction Des Applications Harmoniques Non Rigides D'un Tore Dans La Sphère

Zusammenfassung Le but de Particle present est l'étude des déformations harmoniques géodésiques. f d'un tore dans la sphere avec la densité d'énergie 1/2 . Nous montrons que tout les deformations harmoniques sont isométriques, c'est-à-dire que l'application f est rigide, si et seulement si f est totalement géodésique. De plus nous calculons les espaces de deformations harmoniques pour la 3-sphère.     

一类交错网格的Gauss型格式

本文在交错网格的情况下 ,利用 Gauss型求积公式构造了一类不需解 Riemann问题的求解一维单个双曲守恒律的二阶显式 Gauss型差分格式 ,证明了该格式在CFL条件限制下为 TVD格式 ,并证明了这类格式的收敛性 ,然后将格式推广到方程组的情形 .由于在交错网格的情况下构造的这类差分格式 ,不需要求解 Riemann问题 ,因此这类格式与诸如 Harten等的 TVD格式相比具有如下优点 :由于不需要完整的特征向量系 ,因此可用于求解弱双曲方程组 ,计算更快、编程更加简便等 .