×

Numerical analysis of one dimensional motion of magma without mass forces. (English) Zbl 1490.86001

Summary: In this paper, we present numerical treatment for a compaction-driven Darcian flow viscoelastic rock magma model. This problem is a strongly coupled system of one quasi-linear parabolic equation and one integro-differential equation for the density and the porosity of the flow. The numerical discretization uses cell-centered finite difference method, combined with semi-implicit and implicit time stepping. Implicit-explicit schemes, as well as implicit-explicit iterative algorithms have been developed to solve the corresponding discrete problems. Some properties (positivity, boundness, conservation) of the numerical solutions are investigated. Convergence study of the iteration processes is also presented. The efficiency and the accuracy of the proposed methods are illustrated numerically by test examples with near-real data.

MSC:

86-08 Computational methods for problems pertaining to geophysics
65M06 Finite difference methods for initial value and initial-boundary value problems involving PDEs
76S05 Flows in porous media; filtration; seepage
86A60 Geological problems
PDFBibTeX XMLCite
Full Text: DOI

References:

[1] Bear, J., Dynamics of Fluids in Porous Media (1972), Elsevier: Elsevier New York · Zbl 1191.76001
[2] Papin, A. A.; Tokareva, M., On local solvability of the system of the equations of one dimensional motion of magma, J. Sib. Fed. Univ. Math. Phys., 10, 1, 385-395 (2017) · Zbl 07325365
[3] Simpson, M.; Spiegelman, M.; Weinstein, M. I., Degenerate dispersive equations arising in the study steady of magma dynamics, Nonlinearity, 20, 21-40 (2007) · Zbl 1177.86013
[4] Chen, Z.; Huan, G.; Ma, Y., (Computational Methods for Multiphase Flows in Porous Media. Computational Methods for Multiphase Flows in Porous Media, Computational Science and Engineering Series, vol. 2 (2006), SIAM: SIAM Philadelphia, PA) · Zbl 1092.76001
[5] Akhmerova, I. G.; Papin, A. A.; Tokareva, M. A., Mathematical Models of Heterogeneous Media, Part 1 (2012), Altai Gos. Univ. Barnaul, (in Russian)
[6] J. A.D. Connolly Y. Y.Podladchikov, I. G., Compaction-driven fluid flow in viscoelastic rock, Geodin. Acta, 11, 55-84 (1998)
[7] Arbogast, T.; Obeyeskere, M.; Wheeler, M. F., Numerical methods for the simulation of flow in root-soil systems, SIAM J. Numer. Anal., 30, 6, 1677-1702 (1993) · Zbl 0792.76103
[8] Arbogast, T.; Hesse, M. A.; Taicher, A. L., Mixed methods for two-phase Darcy-Stokes mixtures of partially melted materials with regions of zero porosity, SIAM J. Sci. Comput., 39, 2, B375-B402 (2017) · Zbl 1383.76330
[9] Radu, F. A.; Kumar, K.; Nordbotten, J. M.; Pop, I. S.; robust, A., Mass conservative scheme for two-phase flow in porous media including Hoelder continuous nonlinearities, IMA J. Numer. Anal., 38, 2, 884-920 (2018) · Zbl 1477.76055
[10] Arbogast, T.; Taicher, A. L., A linear degenerate elliptic equation arising from two-phase mixtures, SIAM J. Numer. Anal., 54, 5, 3105-3122 (2016) · Zbl 1457.65177
[11] Arbogast, T.; Taicher, A. L., A cell-centered finite difference method for a degenerate elliptic equation arising from two-phase mixtures, Comput. Geosci., 21, 4, 701-712 (2017) · Zbl 1396.76065
[12] Brenner, K.; Cancès, C.; Hilhorst, D., Finite volume approximation for an immiscible two-phase flow in porous media with discontinuous capillary pressure, Comput. Geosci., 17, 573-597 (2013) · Zbl 1392.76035
[13] Cao, X.; Nemadjieu, S. F.; Pop, S., Convergence of a MPFA finite volume scheme for a two phase porous media flow model with dynamic capillarity, IMA J. Numer. Anal. (2018), (published online 27 January 2018), http://dx.doi.org/10.1093/imanum/drx078 · Zbl 1483.76041
[14] Radu, F. A.; Nordbotten, J. M.; Pop, I. S.; Kumar, K., A robust linearization scheme for finite volume based discretizations for simulation of two-phase flow in porous media, J. Comput. Appl. Math., 289, 134-141 (2015) · Zbl 1320.76084
[15] Morency, C.; Huismans, R. S.; Beaumont, C.; Fullsack, P., A numerical model for coupled fluid flow and matrix deformation with applications to disequilibrium compaction and delta stability, J. Geophys. Res., 112, B10407 (2007)
[16] C. Etchegaray, N. Meunier, Numerical solutions of a 2D fluid problem coupled to a nonlinear non-local reaction-advection-diffusion problem for cell crawling migration in a discoidal domain, in: R. Anguelov, M. Lachowicz (Eds.), Mathematical Methods and Models in Biosciences, pp. 122-139, http://dx.doi.org/10.11145/texts.2018.03.113; C. Etchegaray, N. Meunier, Numerical solutions of a 2D fluid problem coupled to a nonlinear non-local reaction-advection-diffusion problem for cell crawling migration in a discoidal domain, in: R. Anguelov, M. Lachowicz (Eds.), Mathematical Methods and Models in Biosciences, pp. 122-139, http://dx.doi.org/10.11145/texts.2018.03.113 · Zbl 1404.92034
[17] Papin, A. A.; Akhmerova, I. G., Sovability of the system of equations of one-dimensional motion of a heat-conducting two-phase mixture, Math. Notes, 87, 2, 230-243 (2010) · Zbl 1201.35003
[18] Samarskii, A. A., The Theory of Difference Schemes (2001), Marcel Dekker Inc · Zbl 0971.65076
[19] Peña, J. M., M-matrices whose inverses are totally positive, Linear Algebra Appl., 221, 189-193 (1995) · Zbl 0827.15022
[20] Faragó, I.; Horváth, R., Discrete maximum principle and adequate discretizations of linear parabolic problems, SIAM J. Sci. Comput., 28, 2313-2336 (2006) · Zbl 1130.65086
[21] Nikolaevskiy, V. N., Geomechanics and Fluidodynamics with Applications to Reservoir Engineering (1996), Springer Netherlands, 352 p., ISBN 978-90-481-4638-3 (Nedra, Moskow, 1996, 447p., ISBN 5-247-03675-1, original Russian edition)
[22] Turcotte, D. L.; Schubert, G., Geodynamics: Application of Continuum Physics to Geological Problems (1982), Wiley, 450p
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.