. Conservation laws with source terms often have steady states in which the flux gradients are nonzero but exactly balanced by source terms. Many numerical methods (e.g., fractional step methods) have difficulty preserving such steady states and cannot accurately calculate small perturbations of such states. Here a variant of the wave-propagation algorithm is developed which addresses this problem by introducing a Riemann problem in the center of each grid cell whose flux difference exactly cancels the source term. This leads to modified Riemann problems at the cell edges in which the jump now corresponds to perturbations from the steady state. Computing waves and limiters based on the solution to these Riemann problems gives high-resolution results. The 1D and 2D shallow water equations for flow over arbitrary bottom topography are use as an example, though the ideas apply to many other systems. The method is easily implemented in the software package clawpack. Keywords: Godunov meth...

Abstract. This paper investigates the behavior of numerical schemes for nonlinear conservation laws with source terms. We concentrate on two significant examples: relaxation approximations and genuinely nonhomogeneous scalar laws. The main tool in our analysis is the extensive use of weak limits and nonconservative products which allow us to describe accurately the operations achieved in practice when using Riemann-based numerical schemes. Some illustrative and relevant computational results are provided. 1.

This paper deals with the diusion limit of a kinetic equation where the collisions are modeled by a Lorentz type operator. The main aim is to construct a discrete scheme to approximate this equation which gives for any value of the Knudsen number, and in particular at the diusive limit, the right discrete diusion equation with the same value of the diusion coecient as in the continuous case. We are also naturally interested with a discretization which can be used with few velocity discretization points, in order to reduce the cost of computation.

Abstract not found

We are concerned with efficient numerical simulation of the radiative transfer equations...

Abstract. This paper is concerned with the development of high order methods for the numerical approximation of one-dimensional nonconservative hyperbolic systems. In particular, we are interested in high order extensions of the generalized Roe methods introduced by I. Toumi in 1992, based on WENO reconstruction of states. We also investigate the well-balanced properties of the resulting schemes. Finally, we will focus on applications to shallow-water systems. 1.

incline, when the Froude number is above two. The goal of this paper is to analyze the behavior of numerical approximations to a model roll wave equation u(x; 0) = u 0 (x); which arises as a weakly nonlinear approximation of the shallow water equations. The main difficulty associated with the numerical approximation of this problem is its linear instability. Numerical round-off error can easily overtake the numerical solution and yields false roll wave solution at the steady state.

We derive and study Well-Balanced schemes for quasimonotone discrete kinetic models. By means of a rigorous localization procedure, we reformulate the collision terms as nonconservative products and solve the resulting Riemann problem whose solution is self-similar. The construction of an Asymptotic Preserving (AP) Godunov scheme is straightforward and various compactness properties are established within different scalings. At last, some computational results are supplied to show that this approach is realizable and ecient on concrete 2 × 2 models.

We describe the extension of high resolution finite volume methods and adaptive refinement for the shallow water equations in the context of tsunami modeling. Godunov-type methods have been used extensively for modeling the shallow water equations in many contexts, however, tsunami modeling presents some unique challenges that must be overcome. We describe some of the specific difficulties associated with tsunami modeling, and summarize some numerical approaches that we have used to overcome these challenges. For instance, we have developed a well-balanced Riemann solver that is appropriate in the deep ocean regime as well as robust in near-shore and dry regions. Additionally, we have extended adaptive refinement algorithms to this application. We briefly describe some of the modifications necessary for using these adaptive methods for tsunami modeling.

Abstract. We propose and prove convergence of a front tracking method for scalar conservation laws with source term. The method is based on writing the single conservation law as a 2 × 2 quasilinear system without a source term, and employ the solution of the Riemann problem for this system in the front tracking procedure. In this way the source term is processed in the Riemann solver, and one avoids using operator splitting. Since we want to treat the resonant regime, classical arguments for bounding the total variation of numerical solutions do not apply here. Instead compactness of a sequence of front tracking solutions is achieved using a variant of the singular mapping technique invented by Temple [69]. The front tracking method has no CFL–condition associated with it, and it does not discriminate between stiff and non-stiff source terms. This makes it an attractive approach for stiff problems, as is demonstrated in numerical examples. In addition, the numerical examples show that the front tracking method is able to preserve steady–state solutions (or achieving them in the long time limit) with good accuracy. 1.