Hostname: page-component-cd9895bd7-hc48f Total loading time: 0 Render date: 2024-12-28T13:37:55.205Z Has data issue: false hasContentIssue false
  • English
  • Français

The Mass-Preserving S-DDM Scheme for Two-Dimensional Parabolic Equations

Part of: Flows in porous media; filtration; seepage Partial differential equations, initial value and time-dependent initial-boundary value problems

Published online by Cambridge University Press:  01 February 2016

Zhongguo Zhou  and
Dong Liang
Zhongguo Zhou
Affiliation:
School of Mathematics, Shandong University, Jinan, Shandong, 250100, China
Dong Liang*
Affiliation:
School of Mathematics, Shandong University, Jinan, Shandong, 250100, China Department of Mathematics and Statistics, York University, Toronto, Ontario, M3J1P3, Canada
*
*Corresponding author. Email addresses:zhouzhongguo@mail.sdu.edu.cn (Z.Zhou), dliang@mathstat.yorku.ca (D. Liang)
Get access

Abstract

In the paper, we develop and analyze a new mass-preserving splitting domain decomposition method over multiple sub-domains for solving parabolic equations. The domain is divided into non-overlapping multi-bock sub-domains. On the interfaces of sub-domains, the interface fluxes are computed by the semi-implicit (explicit) flux scheme. The solutions and fluxes in the interiors of sub-domains are computed by the splitting one-dimensional implicit solution-flux coupled scheme. The important feature is that the proposed scheme is mass conservative over multiple non-overlapping sub-domains. Analyzing the mass-preserving S-DDM scheme is difficult over non-overlapping multi-block sub-domains due to the combination of the splitting technique and the domain decomposition at each time step. We prove theoretically that our scheme satisfies conservation of mass over multi-block non-overlapping sub-domains and it is unconditionally stable. We further prove the convergence and obtain the error estimate in L2-norm. Numerical experiments confirm theoretical results.

Keywords

Parabolic equationsS-DDMmass-preservingunconditional stabilityconvergencemulti-blocksnon-overlapping

MSC classification

Secondary: 65M06: Finite difference methods 65M12: Stability and convergence of numerical methods 65M55: Multigrid methods; domain decomposition 76S05: Flows in porous media; filtration; seepage
Type
Research Article
Information
Communications in Computational Physics , Volume 19 , Issue 2 , February 2016 , pp. 411 - 441
Copyright
Copyright © Global-Science Press 2016 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)

References

[1]Aziz, K. and Settari, A., Petroleum Reservoir Simulation, Applied Science Publisher, Ltd., London, 1979.Google Scholar
[2]Bear, J., Hydraulics of Groundwater, McGraw-Hill, New York, 1978.Google Scholar
[3]Chen, Z. X., Huan, G. R. and Ma, Y. L., Computational Methods for Multiphase Flows in Porous Media, Computational Science and Engineering Series, SIAM, Philadelphia, 2 (2006).Google Scholar
[4]Dawson, C. N., Du, Q. and Dupont, T. F., A finite difference domain decomposition algorithm for numerical solution of heat equations, Math. Comp., 57 (1991), 63–71.Google Scholar
[5]Du, C. and Liang, D., An Efficient S-DDM Iterative Approach for Compressible Contamination Fluid Flows in Porous Media, J. Comput. Phys., 229 (2010), 4501–4521.Google Scholar
[6]Du, Q., Mu, M. and Wu, Z. N., Efficient parallel algorithms for parabolic problems, SIAM J. Numer. Anal., 39 (2001), 1469–1487.CrossRefGoogle Scholar
[7]Dryja, M. and Tu, X., A domain decomposition discretization of parabolic problems, Numer. Math., 107 (2007), 625–640.CrossRefGoogle Scholar
[8]Gaiffe, S., Glowinski, R. and Lemonnier, P., Domain decomposition and splitting methods for parabolic equations via a mixed formula, in the 12th international conference on Domain Decomposition, Chiba, Japan, 1999.Google Scholar
[9]Kuznetsov, Y. A., New algorithms for approximate realization of implicit difference scheme, Soviet J. Numer. Anal. Math. Modelling, 3 (1988), 99–114.CrossRefGoogle Scholar
[10]Liang, D. and Du, C., The efficient S-DDM scheme and its analysis for solving parabolic equations, J. Comput. Phys., 272 (2014), 46–69.Google Scholar
[11]Lions, P., On the Schwarz alternating method I., In Domain Decomposition Methods for Partial Differential Equations, Glowinski, R., Golub, G. H., Meurant, G. A. and Periaux, J., eds., SIAM, Philadelphia, 1988, 1–42.Google Scholar
[12]Lions, P., On the Schwarz alternating method II., In Domain Decomposition Methods for Partial Differential Equations, Chan, T., Glowinski, R., Meurant, G. A., Periaux, J. and Widlund, O., eds., SIAM, Philadelphia, 1989, 47–70.Google Scholar
[13]Rivera, W., Zhu, J. and Huddleston, D., An efficient parallel algorithm with application to computational fluid dynamics, Comput. Math. Appl., 45 (2003), 165–188.Google Scholar
[14]Shi, H. S. and Liao, H. L., Unconditional stability of corrected explicit/implicit domain decomposition algorithms for parallel approximation of heat equations, SIAMJ. Numer. Anal., 44 (2006), 1584–1611.Google Scholar
[15]Sheng, Z. Q., Yuan, G. W. and Hang, X. D., Unconditional stability of parallel difference schemes with second order accuracy for parabolic equation, Appl. Math. Comput., 184 (2007), 1015–1031.Google Scholar
[16]Smith, B. F., Bjørstad, P. and Gropp, W., Domain Decomposition: Parallel Multilevel Methods for Elliptic Partial Differential Equations, Cambridge University Press, 1996.Google Scholar
[17]Vabishchevich, P. N., Domain decomposition methods with overlapping subdomains for the time-dependent problems of mathematical physics, Comput. Methods Appl. Math., 8 (2008), 393–405.Google Scholar
[18]Zhu, S. H., Conservative domain decomposition procedure with unconditional stability and second-order accuracy, Appl. Math. Comput., 216 (2010), 3275–3282.Google Scholar
[19]Zhuang, Y. and Sun, X. H., Stabilitized explicit-implicit domain decomposition methods for the numerical solution of parabolic equations, SIAMJ. Sci. Comput., 24 (2002), 335–358.Google Scholar