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.
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
@article{M2AN_2018__52_5_1947_0,
author = {Caucao, Sergio and Gatica, Gabriel N. and Oyarz\'ua, Ricardo},
title = {Analysis of an augmented fully-mixed formulation for the coupling of the {Stokes} and heat equations},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis },
pages = {1947--1980},
publisher = {EDP-Sciences},
volume = {52},
number = {5},
year = {2018},
doi = {10.1051/m2an/2018027},
zbl = {1426.65172},
mrnumber = {3885702},
language = {en},
url = {http://www.numdam.org/articles/10.1051/m2an/2018027/}
}
TY - JOUR AU - Caucao, Sergio AU - Gatica, Gabriel N. AU - Oyarzúa, Ricardo TI - Analysis of an augmented fully-mixed formulation for the coupling of the Stokes and heat equations JO - ESAIM: Mathematical Modelling and Numerical Analysis PY - 2018 SP - 1947 EP - 1980 VL - 52 IS - 5 PB - EDP-Sciences UR - http://www.numdam.org/articles/10.1051/m2an/2018027/ DO - 10.1051/m2an/2018027 LA - en ID - M2AN_2018__52_5_1947_0 ER -
%0 Journal Article %A Caucao, Sergio %A Gatica, Gabriel N. %A Oyarzúa, Ricardo %T Analysis of an augmented fully-mixed formulation for the coupling of the Stokes and heat equations %J ESAIM: Mathematical Modelling and Numerical Analysis %D 2018 %P 1947-1980 %V 52 %N 5 %I EDP-Sciences %U http://www.numdam.org/articles/10.1051/m2an/2018027/ %R 10.1051/m2an/2018027 %G en %F M2AN_2018__52_5_1947_0
Caucao, Sergio; Gatica, Gabriel N.; Oyarzúa, Ricardo. Analysis of an augmented fully-mixed formulation for the coupling of the Stokes and heat equations. ESAIM: Mathematical Modelling and Numerical Analysis , Tome 52 (2018) no. 5, pp. 1947-1980. doi : 10.1051/m2an/2018027. http://www.numdam.org/articles/10.1051/m2an/2018027/
[1] and , Sobolev Spaces, 2nd edn. Vol 140 of Pure and Applied Mathematics (Amsterdam). Elsevier, Academic Press, Amsterdam (2003) xiv+305. | MR | Zbl
[2] , and , An augmented mixed-primal finite element method for a coupled flow-transport problem. ESAIM: M2AN 49 (2015) 1399–1427. | DOI | Numdam | MR | Zbl
[3] and , An extra stress-vorticity formulation of Stokes problem for the Oldroyd viscoelastic model. Numer. Math. 94(2003) 603–622. | DOI | MR | Zbl
[4] , and , Finite element exterior calculus, homological tecniques, and applications. Acta Numer. 15 (2006) 1–155. | DOI | MR | Zbl
[5] and , Mixed finite elements for elasticity. Numer. Math. 92 (2002) 401–419. | DOI | MR | Zbl
[6] , and , Mathematical analysis of differential models for viscoelastic fluids. Rheol. Ser. 5 (1996) 199–236. | DOI
[7] and , 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] and , Mixed and Hybrid Finite Element Methods. Vol. 15 of Springer Series in Computational Mathematics. Springer-Verlag, New York (1991) x+350. | MR | Zbl
[9] , , and , 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] , , and , Highly parallel time integration of viscoelastic flows. J. Non-Newtonian Fluid Mech. 100 (2001) 191–216. | DOI | Zbl
[11] , , and 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] , Modeling of cross-flow heat exchangers with viscoelastic fluids. Nonlinear Anal. Real World Appl. 10 (2009) 3353–3359. | DOI | MR | Zbl
[13] , Linear and Nonlinear Functional Analysis with Applications. Society for Industrial and Applied Mathematics, Philadelphia, PA (2013) xiv+832. | MR | Zbl
[14] , and , 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] , and , 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] , and , An augmented fully-mixed finite element method for the stationary Boussinesq problem. Calcolo 54 (2017) 167–205. | DOI | MR | Zbl
[17] , and , Finite element approximation of the non-isothermal Stokes–Oldroyd equations. Int. J. Numer. Anal. Model. 4 (2007) 425–440. | MR | Zbl
[18] and , Non-isothermal flows of viscoelastic incompressible fluids. Nonlinear Anal. 7 (2001) 919–942. | DOI | MR | Zbl
[19] , Algorithm 832: UMFPACK V4.3 – an unsymmetric-pattern multifrontal method. ACM Trans. Math. Softw. 30 (2004) 196–199. | DOI | MR | Zbl
[20] and , A dual mixed formulation for non-isothermal Oldroyd–Stokes problem. MMNP 6 (2011) 130–156. | MR | Zbl
[21] , A Simple Introduction to the Mixed Finite Element Method: Theory and Applications. SpringerBriefs in Mathematics. Springer, Cham (2014). | MR | Zbl
[22] , and , 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] , and , 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] and , 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] , 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] , New development in FreeFem++. J. Numer. Math. 20 (2012) 251–265. | DOI | MR | Zbl
[27] , , , , and , 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] , and , Analytical investigation for cooling turbine disks with a non-Newtonian viscoelastic fluid. Comput. Math. Appl. 61 (2011) 1728–1738. | DOI | MR | Zbl
[29] , and 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] and , Optimal control of non-isothermal viscoelastic fluid flow. J. Non-Newtonian Fluid Mech. 88 (2000) 261–301. | DOI | Zbl
[31] , Strongly Elliptic Systems and Boundary Integral Equations. Cambridge University Press, Cambridge (2000) xiv+357. | MR | Zbl
[32] , , , and , 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] and , Modelling of non-isothermal viscoelastic flows. J. Non-Newtonian Fluid Mech. 68 (1997) 205–224. | DOI
[34] and , Numerical Approximation of Partial Differential Equations. Vol. 23 of Springer Series in Computational Mathematics. Springer-Verlag, Berlin (1994) xvi+543. | MR | Zbl
[35] and , 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] and , Mixed and hybrid methods, Handbook of Numerical Analysis. Vol. II of Handb. Numer. Anal. II. North-Holland, Amsterdam (1991) 523–639. | MR | Zbl
Cité par Sources :
