Results 1 
5 of
5
Numerical Solution Of Reservoir Flow Models Based On Large Time Step Operator Splitting Algorithms
 FILTRATION IN POROUS MEDIA AND INDUSTRIAL APPLICATIONS, LECTURE NOTES IN MATHEMATICS
, 1999
"... During recent years the authors and collaborators have been involved in an activity related to the construction and analysis of large time step operator splitting algorithms for the numerical simulation of multiphase flow in heterogeneous porous media. The purpose of these lecture notes is to revie ..."
Abstract

Cited by 25 (14 self)
 Add to MetaCart
During recent years the authors and collaborators have been involved in an activity related to the construction and analysis of large time step operator splitting algorithms for the numerical simulation of multiphase flow in heterogeneous porous media. The purpose of these lecture notes is to review some of this activity. We illustrate the main ideas behind these novel operator splitting algorithms for a basic twophase flow model. Special focus is posed on the numerical solution algorithms for the saturation equation, which is a convection dominated, degenerate convectiondiffusion equation. Both theory and applications are discussed. The general background for the reservoir flow model is reviewed, and the main features of the numerical algorithms are presented. The basic mathematical results supporting the numerical algorithms are also given. In addition, we present some results from the BV solution theory for quasilinear degenerate parabolic equations, which provides the correct ...
An Unconditionally Stable Method For The Euler Equations
 J. COMPUT. PHYS
, 1999
"... We discuss how to combine a front tracking method with dimensional splitting to solve numerically systems of conservation laws in two space dimensions. In addition we present an adaptive grid refinement strategy. The method is unconditionally stable and allows for moderately high cfl numbers (typ ..."
Abstract

Cited by 7 (4 self)
 Add to MetaCart
We discuss how to combine a front tracking method with dimensional splitting to solve numerically systems of conservation laws in two space dimensions. In addition we present an adaptive grid refinement strategy. The method is unconditionally stable and allows for moderately high cfl numbers (typically 14), and thus it is highly efficient. The method is
Unconditionally Stable Methods For HamiltonJacobi Equations
, 2000
"... We present new numerical methods for constructing approximate solutions to the Cauchy problem for HamiltonJacobi equations of the form u t + H(Dxu) = 0. The methods are based on dimensional splitting and front tracking for solving the associated (nonstrictly hyperbolic) system of conservation laws ..."
Abstract

Cited by 5 (5 self)
 Add to MetaCart
We present new numerical methods for constructing approximate solutions to the Cauchy problem for HamiltonJacobi equations of the form u t + H(Dxu) = 0. The methods are based on dimensional splitting and front tracking for solving the associated (nonstrictly hyperbolic) system of conservation laws p t + DxH(p) = 0, where p = Dxu. In particular, our methods depends heavily on a front tracking method for onedimensional scalar conservation laws with discontinuous coecients. The proposed methods are unconditionally stable in the sense that the time step is not limited by the space discretization and they can be viewed as "large time step" Godunov type (or front tracking) methods. We present several numerical examples illustrating the main features of the proposed methods. We also compare our methods with several methods from the literature.
Numerical Solution Of The Polymer System By Front Tracking
, 1999
"... . The paper describes the application of front tracking to the polymer system, an example of a nonstrictly hyperbolic system. Front tracking computes piecewise constant approximations based on approximate Riemann solutions and exact tracking of waves. It is wellknown that the front tracking method ..."
Abstract

Cited by 5 (4 self)
 Add to MetaCart
. The paper describes the application of front tracking to the polymer system, an example of a nonstrictly hyperbolic system. Front tracking computes piecewise constant approximations based on approximate Riemann solutions and exact tracking of waves. It is wellknown that the front tracking method may introduce a blowup of the initial total variation for initial data along the curve where the two eigenvalues of the hyperbolic system are identical. It is demonstrated by numerical examples that the method converges to the correct solution after a finite time that decreases with the discretization parameter. For multidimensional problems, front tracking is combined with dimensional splitting and numerical experiments indicate that large splitting steps can be used without loss of accuracy. Typical CFL numbers are in the range 10 to 20, and comparisons with Riemann free, highresolution methods confirm that the high efficiency of front tracking. The polymer system, coupled with an ellip...
Front tracking for scalar balance equations
 J. Hyperbolic Differ. Equ
"... 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 f ..."
Abstract

Cited by 3 (2 self)
 Add to MetaCart
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 nonstiff 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.