Citation: | Haochen Liu, Pengzhan Huang. A TWO-GRID DECOUPLED FINITE ELEMENT METHOD FOR THE STATIONARY CLOSED-LOOP GEOTHERMAL SYSTEM[J]. Journal of Applied Analysis & Computation, 2023, 13(4): 1837-1851. doi: 10.11948/20220258 |
A two-grid decoupled finite element method is proposed and analyzed for the stationary closed-loop geothermal model, which is coupled by the Navier-Stokes/Darcy equations and the heat equations with some interface conditions. The main idea of the proposed method is to solve the nonlinear problem on a coarse grid to obtain an initial approximation, then solve the decoupled, linear problem on a fine grid. Hence, the original problem is solved by two subsystems using the two-grid technique, which will save computational time. Moreover, the stability of the proposed method is proved, and numerical examples are presented.
[1] | R. A. Adams and J. J. F. Fournier, Sobolev spaces, Academic Press, New York, 2003. |
[2] | E. Burman and M. Fernández, Stabilization of explicit coupling in fluid-structure interaction involving fluid incompressibility, Comput. Methods Appl. Mech. Engrg., 2009, 198(5–8), 766–784. doi: 10.1016/j.cma.2008.10.012 |
[3] | G. Du, Q. Li and Y. Zhang, A two-grid method with backtracking for the mixed Navier–Stokes/Darcy model, Numer. Meth. Part. Differ. Equs., 2020, 36(6), 1601–1610. doi: 10.1002/num.22493 |
[4] | B. Du, H. Su and X. Feng, Two-level variational multiscale method based on the decoupling approach for the natural convection problem, Int. Commun. Heat Mass Transf., 2015, 61, 128–139. doi: 10.1016/j.icheatmasstransfer.2014.12.004 |
[5] | J. Fang, P. Huang and Y. Qin, A two-level finite element method for the steady-state Navier-Stokes/Darcy model, J. Korean Math. Society, 2020, 57, 915–933. |
[6] | M. Fernández, J. Gerbeau and S. Smaldone, Explicit coupling schemes for a fluid-fluid interaction problem arising in hemodynamics, SIAM J. Sci. Comput., 2014, 36(6), A2557–A2583. doi: 10.1137/130948653 |
[7] | Y. He, Y. Zhang, Y. Shang and H. Xu, Two-level newton iterative method for the 2D/3D steady Navier-Stokes equations, Numer. Meth. Part. Differ. Equs., 2012, 28(5), 1620–1642. doi: 10.1002/num.20695 |
[8] | X. Hu, P. Huang and X. Feng, Two-grid method for Burgers' equation by new mixed finite element schemes, Math. Model. Anal., 2014, 19, 1–17. doi: 10.3846/13926292.2014.892902 |
[9] | P. Huang, A two-level stabilized Oseen iterative method for stationary conduction-convection equations, Math. Rep., 2014, 16, 285–293. |
[10] | P. Huang, An efficient two-level finite element algorithm for the natural convection equations, Appl. Numer. Math., 2017, 118, 75–86. doi: 10.1016/j.apnum.2017.02.012 |
[11] | P. Huang and X. Feng, Error estimates for two-level penalty finite volume method for the stationary Navier-Stokes equations, Math. Meth. Appl. Sci., 2013, 36, 1918–1928. doi: 10.1002/mma.2736 |
[12] | P. Huang, X. Feng and Y. He, Two-level defect-correction Oseen iterative stabilized finite element methods for the stationary Navier-Stokes equations, Appl. Math. Model., 2013, 37, 728–741. doi: 10.1016/j.apm.2012.02.051 |
[13] | P. Huang, X. Feng and D. Liu, Two-level stabilized method based on three corrections for the stationary Navier-Stokes equations, Appl. Numer. Math., 2012, 62, 988–1001. doi: 10.1016/j.apnum.2012.03.006 |
[14] | P. Huang, X. Feng and D. Liu, Two-level stabilized method based on Newton iteration for the steady Smagorinsky model, Nonlinear Anal. Real World Appl., 2013, 14, 1795–1805. doi: 10.1016/j.nonrwa.2012.11.011 |
[15] | P. Huang, X. Feng and H. Su, Two-level defect-correction locally stabilized finite element method for the steady Navier-Stokes equations, Nonlinear Anal. Real World Appl., 2013, 14, 1171–1181. doi: 10.1016/j.nonrwa.2012.09.008 |
[16] | P. Huang, Y. He and X. Feng, Two-level stabilized finite element method for Stokes eigenvalue problem, Appl. Math. Mech. (English Edition), 2012, 33(5), 621–630. doi: 10.1007/s10483-012-1575-7 |
[17] | P. Huang, Y. He and X. Feng, Convergence and stability of two-level penalty mixed finite element method for the stationary Navier-Stokes equations, Frontiers Math. China, 2013, 8, 837–854. doi: 10.1007/s11464-013-0257-2 |
[18] | W. Layton and W. Lenferink, Two-level Picard and modified Picard methods for the Navier-Stokes equations, Appl. Math. Comput., 1995, 69(2–3), 263–274. |
[19] | W. Layton, F. Schieweck and I. Yotov, Coupling fluid flow with porous media flow, SIAM J. Numer. Anal., 2002, 40(6), 2195–2218. doi: 10.1137/S0036142901392766 |
[20] | H. Liu, P. Huang and Y. He, Well-posedness and finite element approximation for the steady-state closed-loop geothermal system, 2022, Submitted. |
[21] | M. Mahbub, X. He, N. Nasu, C. Qiu, Y. Wang and H. Zheng, A coupled multiphysics model and a decoupled stabilized finite element method for the closed-loop geothermal system, SIAM J. Sci. Comput., 2020, 42(4), B951–B982. doi: 10.1137/19M1293533 |
[22] | M. Mu and J. Xu, A two-grid method of a mixed Stokes-Darcy model for coupling fluid flow with porous media flow, SIAM J. Numer. Anal., 2007, 45(5), 1801–1813. doi: 10.1137/050637820 |
[23] | L. Wang, J. Li and P. Huang, An efficient two-level algorithm for the 2D/3D stationary incompressible magnetohydrodynamics based on the finite element method, Int. Commun. Heat Mass Transf., 2018, 98, 183–190. doi: 10.1016/j.icheatmasstransfer.2018.02.019 |
[24] | J. Xu, A novel two-grid method for semilinear elliptic equations, SIAM J. Sci. Comput., 1994, 15(1), 231–237. doi: 10.1137/0915016 |
[25] | J. Xu, Two-grid discretization techniques for linear and nonlinear PDEs, SIAM J. Numer. Anal., 1996, 33(5), 1759–1777. doi: 10.1137/S0036142992232949 |
[26] | Y. Zhang, L. Shan and Y. Hou, Well-posedness and finite element approximation for the convection model in superposed fluid and porous layers, SIAM J. Numer. Anal., 2020, 58(1), 541–564. doi: 10.1137/19M1241532 |
[27] | T. Zhang, X. Zhao and P. Huang, Decoupled two level finite element methods for the steady natural convection problem, Numer. Algor., 2015, 68(4), 837–866. doi: 10.1007/s11075-014-9874-4 |
Temperature distribution with different Darcy numbers. Left:
Temperature distribution with different horizontal pipelines. Left: length =4; Middle: length =2; Right: length =1.
Temperature distribution with different injection temperatures. Left: injection temperature