A locally adaptive time-stepping algorithm for\ud petroleum reservoir simulations

An algorithm for locally adapting the step-size for large scale finite volume simulations of multi-phase flow in petroleum reservoirs is suggested which allows for an “all-in-one” implicit calculation of behaviour over a very large time scale. Some numerical results for simple two-phase flow in one space dimension illustrate the promise of the algorithm, which has also been applied to very simple 3D cases. A description of the algorithm is presented here along with early results. Further development of the technique is hoped to facilitate useful scaling properties

Update on Multirate Timestepping Methods for Hyperbolic Conservation Laws

This paper constructs multirate time discretizations for hyperbolic conservation laws that allow different timesteps to be used in different parts of the spatial domain. The proposed family of discretizations is second order accurate in time and has conservation and linear and nonlinear stability properties under local CFL conditions. Multirate timestepping avoids the necessity to take small global timesteps (restricted by the largest value of the Courant number on the grid) and therefore results in more efficient algorithms. Numerical results obtained for the advection and Burgers equations confirm the theoretical findings

Space-Time Adaptive Solution of First Order PDES

An explicit time-stepping method is developed for adaptive solution of time-dependent partial differential equations with first order derivatives. The space is partitioned into blocks and the grid is refined and coarsened in these blocks. The equations are integrated in time by a Runge-Kutta-Fehlberg method. The local errors in space and time are estimated and the time and space steps are determined by these estimates. The error equation is integrated to obtain global errors of the solution. The method is shown to be stable if one-sided space discretizations are used. Examples such as the wave equation, Burgers' equation, and the Euler equations in one space dimension with discontinuous solutions illustrate the method