A Characteristic Finite Difference Domain Decomposition Substructuring Method for Immiscible Displacement Problems in Three-Dimensional Porous Media

Abstract:

Two-phase immiscible displacement in porous media is described by a coupled nonlinear system of an elliptic equation (for the pressure) and a parabolic equation (for the saturation). For the saturation changes much rapidly than the pressure, a more accurate solution (in both time and space) should be illustrated in practical numerica simulaiton for the former unknown. In this paper we present a seven-point central finite difference scheme to simulate the pressure and a characteristic finite difference combinng with domain decomposition method for the saturation equation. This method consists of reduced two-dimensional computation on the subdomain interface boundaries and fully implicit computation parallelly in subdomains. Aparallel algorithm is outlined and an error estimate in discrete norm is derived by introducing new inner products and norms. At the end of this paper, numerical experiments are presented in order to demonstrate theoretical results and the efficiency.

• Periodical: Advanced Materials Research (Volume 629)
• Pages: 915-919
• DOI: https://doi.org/10.4028/www.scientific.net/AMR.629.915
• Online since: December 2012
• Authors: Chang Feng Li
• Keywords: Characteristic Difference, Domain Decomposition, Error Estimates, Parallel Computation, Two-Phase Displacement

References

[1] Douglas Jr, J. Finite difference methods for two-phase incompressible flow in porous media. SIAM J. Numer. Anal., 1983; 20 (4): 681-696.

DOI: 10.1137/0720046

[2] Russell TF. Time stepping along characteristics with incomplete iteration for a Galerkin approximation of miscible displacement in porous media. SIAM J. Numer. Anal., 1985; 22 (5): 970-1013.

DOI: 10.1137/0722059

[3] Ewing RE, Russell TF, Wheeler MF. Convergence analysis of approximation of miscible displacement in porous media by mixed finite elements and a modified method of characteristics. Computational Methods in Applied Mechanics and Engineering, 1984; 47: 73-92.

DOI: 10.1016/0045-7825(84)90048-3

[4] Douglas Jr., J, Russell TF. Numerical Methods for Convection-dominated Diffusion Problems Based on Combining the Method of Characteristics with Finite Element or Finite Difference Procedures. SIAM J. Numer. Anal., 1982; 19 (5): 871-885.

DOI: 10.1137/0719063

[5] Dawson CN, Du Q, Dupont TF. A finite difference domain decomposition algorithm for numerical solution of the heat equation. Math. Comp., 1991; 57 (195): 63-71.

DOI: 10.1090/s0025-5718-1991-1079011-4

[6] Dawson CN, Dupont TF. Explicit/Implicit conservative Galerkin domain decomposition procedures for parabolic problems. Math. Comp., 1992; 58 (197): 21-34.

DOI: 10.1090/s0025-5718-1992-1106964-9

[7] DU Q, MU M, WU ZN. Efficient parallel algorithms for parabolic problems. SIAM J. Numer. Anal., 2001; 39 (5): 1469-1487.

DOI: 10.1137/s0036142900381710

[8] Yuan YR. Finite difference method and analysis for three-dimensional semiconductor device of heat conduction. Science in China (Series A), 1996; 26 (11): 973-983. (in Chinese).

[9] Li CF, Yuan YR. Domain decomposition with characteristic finite difference method for two-phase displacement problems. Applied Numerical Mathematics, 2008; 58 (9): 1262-1273.

DOI: 10.1016/j.apnum.2007.07.006

[10] Li CF, Yuan YR. Nonoverlapping Domain Decomposition Characteristic Finite Differences for Three-Dimensional Convection-Diffusion Equations. Numerical Methods for Partial Differential Equations, 2012; 28 (1): 17-37.

DOI: 10.1002/num.20605