Premium
An upwind center difference parallel method and numerical analysis for the displacement problem with moving boundary
Author(s) -
Li Changfeng,
Yuan Yirang
Publication year - 2019
Publication title -
numerical methods for partial differential equations
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.901
H-Index - 61
eISSN - 1098-2426
pISSN - 0749-159X
DOI - 10.1002/num.22369
Subject(s) - mathematics , upwind scheme , mathematical analysis , discretization , domain decomposition methods , finite difference method , convection–diffusion equation , finite difference , boundary value problem , nonlinear system , finite element method , physics , quantum mechanics , thermodynamics
A nonlinear coupled mathematical system of two‐phase seepage flow displacement is discussed in this paper including an elliptic equation for the pressure and a convection‐dominated diffusion equation for the saturation. In fact, the boundary of an underground region where the fluid flows through is nonstationary. So a moving boundary should be considered. The saturation equation is convection‐dominated, therefore the method of upwind finite difference is introduced for the accurate computation. The upwind approximation could eliminate numerical oscillation and strong stability is shown. Since the computational work of saturation is larger than the pressure, the authors apply a parallel method, decomposing the whole domain into several nonoverlapping subdomains, to simplify the computation. A domain decomposition method coupled with upwind differences is presented for the saturation. The pressure equation is discretized by a five‐point center finite difference method. By using a transformation and defining new inner products and norms, error estimates in l 2 norm is discussed. Finally, two experimental tests are given to illustrate the efficiency and accuracy of the parallel algorithm.