Analysis of an augmented fully-mixed formulation for the coupling of the Stokes and heat equations

Caucao, Sergio; Gatica, Gabriel N.; Oyarzúa, Ricardo

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 52 (2018) no. 5, pp. 1947-1980.

We introduce and analyse an augmented mixed variational formulation for the coupling of the Stokes and heat equations. More precisely, the underlying model consists of the Stokes equation suggested by the Oldroyd model for viscoelastic flow, coupled with the heat equation through a temperature-dependent viscosity of the fluid and a convective term. The original unknowns are the polymeric part of the extra-stress tensor, the velocity, the pressure, and the temperature of the fluid. In turn, for convenience of the analysis, the strain tensor, the vorticity, and an auxiliary symmetric tensor are introduced as further unknowns. This allows to join the polymeric and solvent viscosities in an adimensional viscosity, and to eliminate the polymeric part of the extra-stress tensor and the pressure from the system, which, together with the solvent part of the extra-stress tensor, are easily recovered later on through suitable postprocessing formulae. In this way, a fully mixed approach is applied, in which the heat flux vector is incorporated as an additional unknown as well. Furthermore, since the convective term in the heat equation forces both the velocity and the temperature to live in a smaller space than usual, we augment the variational formulation by using the constitutive and equilibrium equations, the relation defining the strain and vorticity tensors, and the Dirichlet boundary condition on the temperature. The resulting augmented scheme is then written equivalently as a fixed-point equation, so that the well-known Schauder and Banach theorems, combined with the Lax-Milgram theorem and certain regularity assumptions, are applied to prove the unique solvability of the continuous system. As for the associated Galerkin scheme, whose solvability is established similarly to the continuous case by using the Brouwer fixed-point and Lax–Milgram theorems, we employ Raviart–Thomas approximations of order k for the stress tensor and the heat flux vector, continuous piecewise polynomials of order ≤ k + 1 for velocity and temperature, and piecewise polynomials of order ≤ k for the strain tensor and the vorticity. Finally, we derive optimal a priori error estimates and provide several numerical results illustrating the good performance of the scheme and confirming the theoretical rates of convergence.

Classification : 65N30, 65N12, 65N15, 35Q79, 80A20, 76R05, 76D07
Mots clés : Coupling of Stokes and heat equations, stress-velocity formulation, fixed-point theory, augmented fully-mixed formulation, mixed finite element methods, a priori error analysis

[1] R.A. Adams and J.J.F. Fournier, Sobolev Spaces, 2nd edn. Vol 140 of Pure and Applied Mathematics (Amsterdam). Elsevier, Academic Press, Amsterdam (2003) xiv+305. | MR | Zbl

[2] M. Alvarez, G.N. Gatica and R. Ruiz-Baier, An augmented mixed-primal finite element method for a coupled flow-transport problem. ESAIM: M2AN 49 (2015) 1399–1427. | DOI | Numdam | MR | Zbl

[3] M. Amara and J. Baranger, An extra stress-vorticity formulation of Stokes problem for the Oldroyd viscoelastic model. Numer. Math. 94(2003) 603–622. | DOI | MR | Zbl

[4] D.N. Arnold, R.S. Falk and R. Winther, Finite element exterior calculus, homological tecniques, and applications. Acta Numer. 15 (2006) 1–155. | DOI | MR | Zbl

[5] D.N. Arnold and R. Winther, Mixed finite elements for elasticity. Numer. Math. 92 (2002) 401–419. | DOI | MR | Zbl

[6] J. Baranger, C. Guillopé and J.-C. Saut, Mathematical analysis of differential models for viscoelastic fluids. Rheol. Ser. 5 (1996) 199–236. | DOI

[7] J. Baranger and D. Sandri, A formulation of Stokes’s problem and the linear elasticity equations suggested by the Oldroyd model for viscoelastic flow. RAIRO Modél. Math. Anal. Numér. 26 (1992) 331–345. | DOI | Numdam | MR | Zbl

[8] F. Brezzi and M. Fortin, Mixed and Hybrid Finite Element Methods. Vol. 15 of Springer Series in Computational Mathematics. Springer-Verlag, New York (1991) x+350. | MR | Zbl

[9] J. Camaño, G.N. Gatica, R. Oyarzúa and R. Ruiz-Baier, An augmented stress-based mixed finite element method for the steady state Navier–Stokes equations with nonlinear viscosity. Numer. Methods Part. Differ. Equ. 33 (2017) 1692–1725. | DOI | MR | Zbl

[10] A.E. Caola, Y.L. Joo, R.C. Armstrong and R.A. Brown, Highly parallel time integration of viscoelastic flows. J. Non-Newtonian Fluid Mech. 100 (2001) 191–216. | DOI | Zbl

[11] S. Caucao, G.N. Gatica, R. Oyarzúa and I. Šebestová A fully-mixed finite element method for the Navier–Stokes/Darcy coupled problem with nonlinear viscosity. J. Numer. Math. 25 (2017) 55–88. | DOI | MR | Zbl

[12] T. Chinyoka, Modeling of cross-flow heat exchangers with viscoelastic fluids. Nonlinear Anal. Real World Appl. 10 (2009) 3353–3359. | DOI | MR | Zbl

[13] P.G. Ciarlet, Linear and Nonlinear Functional Analysis with Applications. Society for Industrial and Applied Mathematics, Philadelphia, PA (2013) xiv+832. | MR | Zbl

[14] E. Colmenares, G.N. Gatica and R. Oyarzúa, Fixed point strategies for mixed variational formulations of the stationary Boussinesq problem. C. R. Math. Acad. Sci. Paris 354 (2016) 57–62. | DOI | MR | Zbl

[15] E. Colmenares, G.N. Gatica and R. Oyarzúa, Analysis of an augmented mixed-primal formulation for the stationary Boussinesq problem. Numer. Methods Part. Differ. Equ. 32 (2016) 445–478. | DOI | MR | Zbl

[16] E. Colmenares, G.N. Gatica and R. Oyarzúa, An augmented fully-mixed finite element method for the stationary Boussinesq problem. Calcolo 54 (2017) 167–205. | DOI | MR | Zbl

[17] C. Cox, H. Lee and D. Szurley, Finite element approximation of the non-isothermal Stokes–Oldroyd equations. Int. J. Numer. Anal. Model. 4 (2007) 425–440. | MR | Zbl

[18] S. Damak and C. Guillopé, Non-isothermal flows of viscoelastic incompressible fluids. Nonlinear Anal. 7 (2001) 919–942. | DOI | MR | Zbl

[19] T.A. Davis, Algorithm 832: UMFPACK V4.3 – an unsymmetric-pattern multifrontal method. ACM Trans. Math. Softw. 30 (2004) 196–199. | DOI | MR | Zbl

[20] M. Farhloul and A. Zine, A dual mixed formulation for non-isothermal Oldroyd–Stokes problem. MMNP 6 (2011) 130–156. | MR | Zbl

[21] G.N. Gatica, A Simple Introduction to the Mixed Finite Element Method: Theory and Applications. SpringerBriefs in Mathematics. Springer, Cham (2014). | MR | Zbl

[22] 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. Eng. 264 (2013) 23–48. | DOI | MR | Zbl

[23] 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. Eng. 200 (2011) 1619–1636. | DOI | MR | Zbl

[24] V. Girault and P.-A. Raviart, Finite element methods for Navier–Stokes equations: Theory and algorithms. Vol. 5 of Springer Series in Computational Mathematics. Springer-Verlag, Berlin (1986) x+374. | MR | Zbl

[25] B. Hjertager, Multi-fluid CFD Analysis of Chemical Reactors, Multiphase Reacting Flows: Modelling and Simulation. CISM Courses and Lectures, SpringerWienNewYork, Vienna (2007) Vol. 492, 125–179. | MR | Zbl

[26] F. Hecht, New development in FreeFem++. J. Numer. Math. 20 (2012) 251–265. | DOI | MR | Zbl

[27] Y.L. Joo, J. Sun, M.D. Smith, R.C. Armstrong, R.A. Brown and R.A. Ross, Two-dimensional numerical analysis of non-isothermal melt spinning with and without phase transition. J. Non-Newtonian Fluid Mech. 102 (2002) 37–70. | DOI | Zbl

[28] F. Khani, M.T. Darvishi and R. Gorla, Analytical investigation for cooling turbine disks with a non-Newtonian viscoelastic fluid. Comput. Math. Appl. 61 (2011) 1728–1738. | DOI | MR | Zbl

[29] A. Kufner, O. Jhon and S. Fučík Function Spaces. Monographs and Textbooks on Mechanics of Solids and Fluids; Mechanics: Analysis. Noordhoff International Publishing, Academia, Leyden, Prague (1977) xv+454. | MR | Zbl

[30] K. Kunisch and X. Marduel, Optimal control of non-isothermal viscoelastic fluid flow. J. Non-Newtonian Fluid Mech. 88 (2000) 261–301. | DOI | Zbl

[31] W. Mclean, Strongly Elliptic Systems and Boundary Integral Equations. Cambridge University Press, Cambridge (2000) xiv+357. | MR | Zbl

[32] Y. Mu, G. Zhao, X. Wu, L. Hang and H. Chu, Continuous modeling and simulation of flow-swell-crystallization behaviors of viscoelastic polymer melts in the hollow profile extrusion process. Appl. Math. Model. 39 (2015) 1352–1368. | DOI | MR | Zbl

[33] G.W.M. Peters and F.P.T. Baaijens, Modelling of non-isothermal viscoelastic flows. J. Non-Newtonian Fluid Mech. 68 (1997) 205–224. | DOI

[34] A. Quarteroni and A. Valli, Numerical Approximation of Partial Differential Equations. Vol. 23 of Springer Series in Computational Mathematics. Springer-Verlag, Berlin (1994) xvi+543. | MR | Zbl

[35] P.-A. Raviart and J.-M. Thomas, Introduction à l’analyse numérique des équations aux dérivées partielles, Collection Mathématiques Appliquées pour la Maîtrise. Masson, Paris (1983) 224. | MR | Zbl

[36] J.E. Roberts and J.-M. Thomas, Mixed and hybrid methods, Handbook of Numerical Analysis. Vol. II of Handb. Numer. Anal. II. North-Holland, Amsterdam (1991) 523–639. | MR | Zbl