Skip to content
Licensed Unlicensed Requires Authentication Published by De Gruyter June 17, 2016

A fully-mixed finite element method for the Navier–Stokes/Darcy coupled problem with nonlinear viscosity

  • Sergio Caucao , Gabriel N. Gatica EMAIL logo , Ricardo Oyarzúa and Ivana Šebestová

Abstract

We propose and analyze an augmented mixed finite element method for the coupling of fluid flow with porous media flow. The flows are governed by a class of nonlinear Navier–Stokes and linear Darcy equations, respectively, and the transmission conditions are given by mass conservation, balance of normal forces, and the Beavers–Joseph–Saffman law. We apply dual-mixed formulations in both domains, and the nonlinearity involved in the Navier–Stokes region is handled by setting the strain and vorticity tensors as auxiliary unknowns. In turn, since the transmission conditions become essential, they are imposed weakly, which yields the introduction of the traces of the porous media pressure and the fluid velocity as the associated Lagrange multipliers. Furthermore, since the convective term in the fluid forces the velocity to live in a smaller space than usual, we augment the variational formulation with suitable Galerkin type terms arising from the constitutive and equilibrium equations of the Navier–Stokes equations, and the relation defining the strain and vorticity tensors. The resulting augmented scheme is then written equivalently as a fixed point equation, so that the well-known Schauder and Banach theorems, combined with classical results on bijective monotone operators, are applied to prove the unique solvability of the continuous and discrete systems. In particular, given an integer k ⩾ 0, piecewise polynomials of degree ⩽ k, Raviart–Thomas spaces of order k, continuous piecewise polynomials of degree ⩽ k + 1, and piecewise polynomials of degree ⩽ k are employed in the fluid for approximating the strain tensor, stress, velocity, and vorticity, respectively, whereas Raviart–Thomas spaces of order k and piecewise polynomials of degree ⩽ k for the velocity and pressure, together with continuous piecewise polynomials of degree ⩽ k + 1 for the traces, constitute feasible choices in the porous medium. Finally, several numerical results illustrating the good performance of the augmented mixed finite element method and confirming the theoretical rates of convergence are reported.

Funding

This work was partially supported by CONICYT-Chile through BASAL project CMM, Universidad de Chile, project Anillo ACT1118 (ANANUM), Fondecyt projects 11121347 and 3150047, and the Becas-Chile Programme for Chilean students; by Centro de Investigación en Ingeniería Matemática (CI2MA), Universidad de Concepción; and by Universidad del Bío-Bío through DIUBB project 151408 GI/VC.

References

[1] R. A. Adams and J. J. Fournier, Sobolev Spaces, 140, Academic Press, 2003, Second edition.Search in Google Scholar

[2] D. N. Arnold, F. Brezzi, and J. Douglas, PEERS: a new mixed finite element for plane elasticity, Japan J. Ind. Appl. Math. 1 (1984), 347–367.10.1007/BF03167064Search in Google Scholar

[3] L. Badea, M. Discacciati, and A. Quarteroni, Numerical analysis of the Navier–Stokes/Darcy coupling, Numer. Math. 115 (2010), 195–227.10.1007/s00211-009-0279-6Search in Google Scholar

[4] F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods, 15, Springer Science & Business Media, 2012.10.1007/978-1-4612-3172-1Search in Google Scholar

[5] F. Brezzi, M. Fortin, and L. D. Marini, Mixed finite element methods with continuous stresses, Math. Models Methods Appl. Sci. 3 (1993), 275–287.10.1142/S0218202593000151Search in Google Scholar

[6] M. Cai, M. Mu, and J. Xu, Numerical solution to a mixed Navier–Stokes/Darcy model by the two-grid approach, SIAM J. Numer. Anal. 47 (2009), 3325–3338.10.1137/080721868Search in Google Scholar

[7] Z. Cai, C. Tong, P.S. Vassilevski, and C. Wang, Mixed finite element methods for incompressible flow: stationary Stokes equations, Numer. Meth. Partial Differ. Equ. 26 (2010), 957–978.10.1137/080718413Search in Google Scholar

[8] Z. Cai, C. Wang, and S. Zhang, Mixed finite element methods for incompressible flow: stationary Navier–Stokes equations, SIAM J. Numer. Anal. 48 (2010), 79–94.10.1137/080718413Search in Google Scholar

[9] Z. Cai and S. Zhang, Mixed methods for stationary Navier–Stokes equations based on pseudostress–pressure–velocity formulation, Math. Comp. 81 (2012), 1903–1927.10.1090/S0025-5718-2012-02585-3Search in Google Scholar

[10] J. Camaño, G.N. Gatica, R. Oyarzúa, R. Ruiz-Baier, and P. Venegas, New fully-mixed finite element methods for the Stokes–Darcy coupling, Comput. Methods Appl. Mech. Engrg. 295 (2015), 362–395.10.1016/j.cma.2015.07.007Search in Google Scholar

[11] J. Camaño, G.N. Gatica, R. Oyarzúa, and G. Tierra, An augmented mixed finite element method for the Navier–Stokes equations with variable viscosity, SIAM J. Numer. Anal. 54 (2016), 1069–1092.10.1137/15M1013146Search in Google Scholar

[12] J. Camaño, R. Oyarzúa, and G. Tierra, Analysis of an augmented mixed-FEM for the Navier–Stokes problem, Math. Comp. 86 (2017), 589–615.10.1090/mcom/3124Search in Google Scholar

[13] Y. Cao, Y. Chu, X. He, and M. Wei, Decoupling the stationary Navier–Stokes–Darcy system with the Beavers–Joseph–Saffman interface condition, In: Abstract and Applied Analysis, Hindawi Publishing Corporation, 2013, Artl ID 136483.10.1155/2013/136483Search in Google Scholar

[14] E. Castillo and R. Codina, Stabilized stress–velocity–pressure finite element formulations of the Navier–Stokes problem for fluids with non-linear viscosity, Comput. Methods Appl. Mech. Engrg. 279 (2014), 554–578.10.1016/j.cma.2014.07.003Search in Google Scholar

[15] A. Çeşmelioğlu and B. Rivière, Analysis of time-dependent Navier–Stokes flow coupled with Darcy flow, J. Numer. Math. 16 (2008), 249–280.10.1515/JNUM.2008.012Search in Google Scholar

[16] A. Cesmelioglu, V. Girault, and B. Rivière, Time-dependent coupling of Navier–Stokes and Darcy flows, ESAIM Math. Model. Numer. Anal. 47 (2013), 539–554.10.1051/m2an/2012034Search in Google Scholar

[17] P. Chidyagwai and B. Rivière, On the solution of the coupled Navier–Stokes and Darcy equations, Comput. Methods Appl. Mech. Engrg. 198 (2009), 3806–3820.10.1016/j.cma.2009.08.012Search in Google Scholar

[18] P. G. Ciarlet, Linear and Nonlinear Functional Analysis with Applications, 130, SIAM, 2013.Search in Google Scholar

[19] E. Colmenares, G. N. Gatica, and R. Oyarzúa, Analysis of an augmented mixed-primal formulation for the stationary Boussinesq problem, Numer. Meth. Partial Differ. Equ. 32 (2016), 445–478.10.1002/num.22001Search in Google Scholar

[20] T. A. Davis, Algorithm 832: UMFPACK V4.3 — an unsymmetric-pattern multifrontal method, ACM Trans. Math. Software (TOMS)30 (2004), 196–199.10.1145/992200.992206Search in Google Scholar

[21] M. Discacciati, Domain decomposition methods for the coupling of surface and groundwater flows, Ph.D. thesis, École Polytechnique Fédérale de Lausanne, Switzerland, 2004.Search in Google Scholar

[22] M. Discacciati, E. Miglio, and A. Quarteroni, Mathematical and numerical models for coupling surface and groundwater flows, Appl. Numer. Math. 43 (2002), 57–74.10.1016/S0168-9274(02)00125-3Search in Google Scholar

[23] M. Discacciati and R. Oyarzúa, A conforming mixed finite element method for the Navier–Stokes–Darcy coupled problem, Numer. Math. (2016), 1–36.10.1007/s00211-016-0811-4Search in Google Scholar

[24] M. Discacciati and A. Quarteroni, Navier–Stokes–Darcy coupling: modeling, analysis, and numerical approximation, Rev. Mat. Complut. 22 (2009), 315–426.10.5209/rev_REMA.2009.v22.n2.16263Search in Google Scholar

[25] J. Douglas and J. P. Wang, An absolutely stabilized finite element method for the Stokes problem, Math. Comp. 52 (1989), 495–508.10.1090/S0025-5718-1989-0958871-XSearch in Google Scholar

[26] L. E. Figueroa, G. N. Gatica, and A. Márquez, Augmented mixed finite element methods for the stationary Stokes equations, SIAM J. Sci. Comp. 31 (2008), 1082–1119.10.1137/080713069Search in Google Scholar

[27] L. P. Franca and T. J. Hughes, Two classes of mixed finite element methods, Comput. Methods Appl. Mech. Engrg. 69 (1988), 89–129.10.1016/0045-7825(88)90168-5Search in Google Scholar

[28] L. P. Franca and R. Stenberg, Error analysis of Galerkin least squares methods for the elasticity equations, SIAM J. Numer. Anal. 28 (1991), 1680–1697.10.1137/0728084Search in Google Scholar

[29] G. Gatica, R. Oyarzúa, and F.-J. Sayas, Analysis of fully-mixed finite element methods for the Stokes–Darcy coupled problem, Math. Comp. 80 (2011), 1911–1948.10.1090/S0025-5718-2011-02466-XSearch in Google Scholar

[30] G. N. Gatica, A Simple Introduction to the Mixed Finite Element Method: Theory and Applications, Springer Science & Business Media, 2014.10.1007/978-3-319-03695-3Search in Google Scholar

[31] G. N. Gatica, N. Heuer, and S. Meddahi, On the numerical analysis of nonlinear twofold saddle point problems, IMA J. Numer. Anal. 23 (2003), 301–330.10.1093/imanum/23.2.301Search in Google Scholar

[32] G. N. Gatica, G. C. Hsiao, and S. Meddahi, A coupled mixed finite element method for the interaction problem between an electromagnetic field and an elastic body, SIAM J. Numer. Anal. 48 (2010), 1338–1368.10.1137/090754212Search in Google Scholar

[33] G. N. Gatica, A. Márquez, and S. Meddahi, Analysis of the coupling of primal and dual-mixed finite element methods for a two-dimensional fluid-solid interaction problem, SIAM J. Numer. Anal. 45 (2007), 2072–2097.10.1137/060660370Search in Google Scholar

[34] G. N. Gatica, A. Márquez, R. Oyarzúa, and R. Rebolledo, Analysis of an augmented fully-mixed approach for the coupling of quasi-Newtonian fluids and porous media, Comput. Methods Appl. Mech. Engrg. 270 (2014), 76–112.10.1016/j.cma.2013.11.017Search in Google Scholar

[35] G. N. Gatica, A. Márquez, and W. Rudolph, A priori and a posteriori error analyses of augmented twofold saddle point formulations for nonlinear elasticity problems, Comput. Methods Appl. Mech. Engrg. 264 (2013), 23–48.10.1016/j.cma.2013.05.010Search in Google Scholar

[36] G. N. Gatica, A. Márquez, and M. A. Sánchez, A priori and a posteriori error analyses of a velocity–pseudostress formulation for a class of quasi-Newtonian Stokes flows, Comput. Methods Appl. Mech. Engrg. 200 (2011), 1619–1636.10.1016/j.cma.2011.01.010Search in Google Scholar

[37] G. N. Gatica, S. Meddahi, and R. Oyarzúa, A conforming mixed finite-element method for the coupling of fluid flow with porous media flow, IMA J. Numer. Anal. 29 (2009), 86–108.10.1093/imanum/drm049Search in Google Scholar

[38] G. N. Gatica, R. Oyarzúa, and F.-J. Sayas, A twofold saddle point approach for the coupling of fluid flow with nonlinear porous media flow, IMA J. Numer. Anal. 32 (2012), 845–887.10.1093/imanum/drr020Search in Google Scholar

[39] G. N. Gatica and W. L. Wendland, Coupling of mixed finite elements and boundary elements for linear and nonlinear elliptic problems, Applicable Analysis63 (1996), 39–75.10.1080/00036819608840495Search in Google Scholar

[40] V. Girault and P.-A. Raviart, Finite Element Methods for Navier–Stokes Equations: Theory and Algorithms, 5, Springer Science & Business Media, 2012.Search in Google Scholar

[41] V. Girault and B. Rivière, DG approximation of coupled Navier–Stokes and Darcy equations by Beaver–Joseph–Saffman interface condition, SIAM J. Numer. Anal. 47 (2009), 2052–2089.10.1137/070686081Search in Google Scholar

[42] P. Hansbo and M. Juntunen, Weakly imposed Dirichlet boundary conditions for the Brinkman model of porous media flow, Appl. Numer. Math. 59 (2009), 1274–1289.10.1016/j.apnum.2008.07.003Search in Google Scholar

[43] F. Hecht, New development in FreeFem++, J. Numer. Math. 20 (2012), 251–266.10.1515/jnum-2012-0013Search in Google Scholar

[44] R. Hiptmair, Finite elements in computational electromagnetism, Acta Numerica11 (2002), 237–339.10.1017/CBO9780511550140.004Search in Google Scholar

[45] C. O. Horgan, Korn’s inequalities and their applications in continuum mechanics, SIAM Rev. 37 (1995), 491–511.10.1137/1037123Search in Google Scholar

[46] J. S. Howell, Dual-mixed finite element approximation of Stokes and nonlinear Stokes problems using trace-free velocity gradients, J. Comp. Appl. Math. 231 (2009), 780–792.10.1016/j.cam.2009.05.002Search in Google Scholar

[47] J. S. Howell and N. J. Walkington, Dual-mixed finite element methods for the Navier–Stokes equations, ESAIM Math. Model. Numer. Anal. 47 (2013), 789–805.10.1051/m2an/2012050Search in Google Scholar

[48] A. Márquez, S. Meddahi, and F.-J. Sayas, Strong coupling of finite element methods for the Stokes–Darcy problem, IMA J. Numer. Anal. 35 (2015), 969–988.10.1093/imanum/dru023Search in Google Scholar

[49] W. C. H. McLean, Strongly Elliptic Systems and Boundary Integral Equations, Cambridge University Press, 2000.Search in Google Scholar

[50] A. Quarteroni and A. Valli, Numerical Approximation of Partial Differential Equations, 23, Springer Science & Business Media, 2008.Search in Google Scholar

[51] L. Zuo and Y. Hou, Numerical analysis for the mixed Navier–Stokes and Darcy Problem with the Beavers–Joseph interface condition, Numer. Meth. Partial Differ. Equ. 31 (2015), 1009–1030. 10.1002/num.21933Search in Google Scholar

Received: 2015-11-13
Revised: 2016-6-8
Accepted: 2016-6-18
Published Online: 2016-6-17
Published in Print: 2017-7-26

© 2017 by Walter de Gruyter Berlin/Boston

Downloaded on 1.5.2024 from https://www.degruyter.com/document/doi/10.1515/jnma-2015-0121/html
Scroll to top button