Motivated by possible generalizations to more complex multiphase multicomponent systems in higher dimensions, we develop an Eulerian-Lagrangian numerical approximation for a system of two conservation laws in one space dimension modeling a simplified two-phase flow problem in a porous medium. The method is based on following tracelines, so it is stable independent of any CFL constraint. The main difficulty is that it is not possible to follow individual tracelines independently. We approximate tracing along the tracelines by using local mass conservation principles and self-consistency. The two-phase flow problem is governed by a system of equations representing mass conservation of each phase, so there are two local mass conservation principles. Our numerical method respects both of these conservation principles over the computational mesh (i.e., locally), and so is a fully conservative traceline method. We present numerical results that demonstrate the ability of the method to handle problems with shocks and rarefactions, and to do so with very coarse spatial grids and time steps larger than the CFL limit. © 2012 Society for Industrial and Applied Mathematics.
|Original language||English (US)|
|Number of pages||1|
|Journal||SIAM Journal on Scientific Computing|
|State||Published - Jan 2012|
Bibliographical noteKAUST Repository Item: Exported on 2020-10-01
Acknowledgements: This author was supported in part by U.S. National Science Foundation grants DMS-0713815 and DMS-0835745, the King Abdullah University of Science and Technology (KAUST) Academic Excellence Alliance program, and the Mathematics Research Promotion Center of Taiwan.This author was supported in part under Taiwan National Science Council grants 96-2115-M-110-002-MY3 and 99-2115-M-110-006-MY3.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.