Abstract
The incompressible miscible displacement of three-dimensional Darcy–Forchheimer flow is discussed in this paper, and the mathematical model is formulated by two partial differential equations, a Darcy–Forchheimer flow equation for the pressure and a convection-diffusion equation for the concentration. The pressure plays an important role and determines the Darcy–Forchheimer velocity. A conservative mixed finite element is used to approximate the pressure and the velocity, and the computational accuracy of Darcy–Forchheimer velocity is improved by one order. A second-order upwind fractional step difference scheme is adopted to obtain the concentration, and numerical oscillation and dispersion are eliminated. Computational work is reduced greatly by decomposing the whole computation into successive one-dimensional subcomputations, where the speedup solvers are used. Applying energy-norm method, operator decomposition, speedup algorithm and induction hypotheses, we derive an optimal second-order estimates in \(L^2\) norm. Numerical experiments are given to show the theoretical accuracy and the applicability for solving actual problems.
Similar content being viewed by others
References
Aziz, K., Settari, A.: Petroleum Reservoir Simulation. Applied Science Publishers, London (1979)
Ewing, R.E., Lazarov, R.D., Lyons, S.L., Papavassiliou, D.V., Pasciak, J., Qin, G.: Numerical well model for non-darcy flow through isotropic porous media. Comput. Geosci. 3(3–4), 185–204 (1999)
Ruth, D., Ma, H.: On the derivation of the Forchheimer equation by means of the averaging theorem. Transp. Porous Media 7(3), 255–264 (1992)
Douglas Jr., J.: Finite difference methods for two-phase incompressible flow in porous media. SIAM J. Numer. Anal. 20(4), 681–696 (1983)
Douglas Jr., J.: Simulation of miscible displacement in porous media by a modified method of characteristic procedure. In: Numerical Analysis, Dundee, Lecture Notes in Mathematics, vol 912, pp. 1982. Springer-Verlag, Berlin (1981)
Ewing, R.E., Russell, T.F., Wheeler, M.F.: Convergence analysis of an approximation of miscible displacement in porous media by mixed finite elements and a modified method of characteristics. Comput. Methods Appl. Mech. Eng. 47(1–2), 73–92 (1984)
Douglas Jr., J., Yuan, Y. R.: Numerical simulation of immiscible flow in porous media based on combining the method of characteristics with mixed finite element procedure. In: Numerical Simulation in Oil Recovery, pp. 119–132. Springer, New York (1986)
Fabrie, P.: Regularity of the solution of Darcy-Forchheimer’s equation. Nonlinear Anal. Theory. Methods Appl. 13(9), 1025–1049 (1989)
Girault, V., Wheeler, M.F.: Numerical discretization of a Darcy Forchheimer model. Numer. Math. 110(2), 161–198 (2008)
Lopez, H., Molina, B., Salas, J.J.: Comparison between different numerical discretization for a Darcy-Forchheimer model. Electron. Trans. Numer. Anal. 34, 187–203 (2009)
Pan, H., Rui, H.: Mixed element method for two-dimensional Darcy-Forchheimer model. J. Sci. Comput. 52, 563–587 (2012)
Pan, H., Rui, H.: A mixed element method for Darcy-Forchheimer incompressible miscible displacement problem. Comput. Methods Appl. Mech. Eng. 264, 1–11 (2013)
H. Rui, H. Pan. A block-centered finite difference method for compressible Darcy Forchheimer model. Submitted for publication
Rui, H., Pan, H.: A block-centered finite difference method for the Darcy-Forchheimer model. SIAM. J Numer. Anal. 50(5), 2612–2631 (2012)
J. Douglas Jr., P. J. Paes-Leme, T. Giorgi. Generalized Forchheimer flow in porous media. In: Boundary value problems for partial differential equations and applications. In: RMA Research Notes Applied Mathematics, vol. 29, pp. 99–111. Masson, Paris (1993)
Park, J.E.: Mixed finite element method for generalized Forchheimer flow in porous media. Numer. Methods Part. Differ. Equ. 21(2), 213–228 (2005)
Axelsson, O., Gustafasson, I.: A modified upwind scheme for convective transport equations and the use of a conjugate gradient method for the solution of non-symmetric systems of equations. J. Inst. Math. Appl. 23, 321–337 (1979)
Ewing, R.E., Lazarov, R.D., Vassilev, A.T.: Finite difference scheme for parabolic problems on a composite grids with refinement in time and space. SIAM J. Numer. Anal. 31(6), 1605–1622 (1994)
Lazarov, R.D., Mischev, I.D., Vassilevski, P.S.: Finite volume methods for convection-diffusion problems. SIAM J. Numer. Anal. 33(1), 31–55 (1996)
Ewing, R.E.: The Mathematics of Reservior Simulation. SIAM, Philadelphia (1983)
Peaceman, D.W.: Fundamental of Numerical Reservoir Simulation. Elsevier, Amsterdam (1980)
Yuan, Y.R.: Theory and Application of Reservoir Numerical Simulation, p. 6. Science Press, Beijing (2013)
Douglas Jr., J., Gann, J.E.: Two order correct difference analogues for the equation of multidimensional heat flow. Math. Comp. 17(81), 71–80 (1963)
Douglas Jr., J., Gunn, J.E.: A general formulation of alternation methods, part I. Parabolic and hyperbolic problems. Numer. Math. 9(5), 428–453 (1964)
Yuan, Y.R.: Time stepping along characteristics for the finite element approximation of compressible miscible displacement in porous media. Math. Numer. Sinica 14(4), 385–400 (1992)
Yuan, Y.R.: Finite difference methods for a compressible miscible displacement problem in porous media. Math. Numer. Sinica 15(1), 16–28 (1993)
Acknowledgements
The authors would like to thank the reviewers for their helpful suggestions.
Author information
Authors and Affiliations
Corresponding author
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Project supported by the National Natural Science Foundation of China (Grant No. 11271231), Natural Science Foundation of Shandong Province (Grant No. ZR2016AM08)
Rights and permissions
About this article
Cite this article
Yuan, Y., Li, C. & Yang, Q. Mixed Finite Element-Second Order Upwind Fractional Step Difference Scheme of Darcy–Forchheimer Miscible Displacement and Its Numerical Analysis. J Sci Comput 86, 24 (2021). https://doi.org/10.1007/s10915-020-01393-9
Received:
Revised:
Accepted:
Published:
DOI: https://doi.org/10.1007/s10915-020-01393-9
Keywords
- Darcy–Forchheimer flow
- Miscible displacement
- Mixed finite element-upwind fractional step difference
- Second order error estimates in \(L^2\)-norm
- Numerical computation