Results 1  10
of
10
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 46 (22 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 ...
Wellbalanced finite volume schemes of arbitrary order of accuracy for shallow water flows
, 2006
"... ..."
Unconditionally stable methods for HamiltonJacobi equations
 J. Comput. Phys
, 2002
"... We present new numerical methods for constructing approximate solutions to the Cauchy problem for HamiltonJacobi equations of the form u t + H (D x u) = 0. The methods are based on dimensional splitting and front tracking for solving the associated (nonstrictly hyperbolic) system of conservation ..."
Abstract

Cited by 12 (8 self)
 Add to MetaCart
(Show Context)
We present new numerical methods for constructing approximate solutions to the Cauchy problem for HamiltonJacobi equations of the form u t + H (D x u) = 0. The methods are based on dimensional splitting and front tracking for solving the associated (nonstrictly hyperbolic) system of conservation laws p t + D x H ( p) = 0, where p = D x u. In particular, our methods depend heavily on a front tracking method for onedimensional scalar conservation laws with discontinuous coefficients. 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 "largetimestep" Godunovtype (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. c 2002 Elsevier Science (USA)
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 10 (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
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 8 (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 5 (2 self)
 Add to MetaCart
(Show Context)
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.
A Second Order Front Tracking Solution of the Euler Equations
"... A second order front tracking method is developed for solving the Euler equations of inviscid fluid dynamics numerically. Front tracking methods are usually limited to first order accuracy, since they are based on a piecewise constant approximation of the solution. Here the second order convergence ..."
Abstract
 Add to MetaCart
A second order front tracking method is developed for solving the Euler equations of inviscid fluid dynamics numerically. Front tracking methods are usually limited to first order accuracy, since they are based on a piecewise constant approximation of the solution. Here the second order convergence is achieved by building a piecewise linear reconstruction of the piecewise constant front tracking solution in a postprocessing step. The linearization is performed by decomposing the piecewise constant solution of the hyperbolic system into its wave components and by linearizing the wave solutions separately. In order to achieve a physically correct linearization, the front types of the previously developed improved front tracking method are employed. It is illustrated numerically for the onedimensional unsteady interacting blast waves problem and a twodimensional supersonic airfoil flow validation study that the proposed front tracking method can achieve second order convergence also in the presents of strong discontinuities. I.
Hydropolis Dam Failure: Impact Assessments
"... By approving this article, the publisher recognizes that the U.S. Government retains nonexclusive, royaltyfree license to publish or reproduce the published form of this contribution, or to allow others to do so, for U.S. Government purposes. Los Alamos National Laboratory requests that the publish ..."
Abstract
 Add to MetaCart
(Show Context)
By approving this article, the publisher recognizes that the U.S. Government retains nonexclusive, royaltyfree license to publish or reproduce the published form of this contribution, or to allow others to do so, for U.S. Government purposes. Los Alamos National Laboratory requests that the publisher identify this article as work performed under the auspices of the U.S. Departmentof Energy. Los Alamos National Laboratory strongly supports academic freedom and a researcher's right to publish; as an institution, however, the Laboratory does not endorse the viewpoint of a publication or guarantee its technical correctness.