A Fully-Mixed Finite Element Method for the n-Dimensional Boussinesq Problem with Temperature-Dependent Parameters

In this paper, we introduce and analyze a high-order, fully-mixed finite element method for the free convection of -dimensional fluids, , with temperature-dependent viscosity and thermal conductivity. The mathematical model is given by the coupling of the equations of continuity, momentum (Navier–St...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Journal of computational methods in applied mathematics Jg. 20; H. 2; S. 187 - 213
Hauptverfasser: Almonacid, Javier A., Gatica, Gabriel N.
Format: Journal Article
Sprache:Englisch
Veröffentlicht: Minsk De Gruyter 01.04.2020
Walter de Gruyter GmbH
Schlagworte:
ISSN:1609-4840, 1609-9389
Online-Zugang:Volltext
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Beschreibung
Zusammenfassung:In this paper, we introduce and analyze a high-order, fully-mixed finite element method for the free convection of -dimensional fluids, , with temperature-dependent viscosity and thermal conductivity. The mathematical model is given by the coupling of the equations of continuity, momentum (Navier–Stokes) and energy by means of the Boussinesq approximation, as well as mixed thermal boundary conditions and a Dirichlet condition on the velocity. Because of the dependence on the temperature of the fluid properties, several additional variables are defined, thus resulting in an augmented formulation that seeks the rate of strain, pseudostress and vorticity tensors, velocity, temperature gradient and pseudoheat vectors, and temperature of the fluid. Using a fixed-point approach, smallness-of-data assumptions and a slight higher-regularity assumption for the exact solution provide the necessary well-posedness results at both continuous and discrete levels. In addition, and as a result of the augmentation, no discrete inf-sup conditions are needed for the well-posedness of the Galerkin scheme, which provides freedom of choice with respect to the finite element spaces. In particular, we suggest a combination based on Raviart–Thomas, Lagrange and discontinuous elements for which we derive optimal a priori error estimates. Finally, several numerical examples illustrating the performance of the method and confirming the theoretical rates of convergence are reported.
Bibliographie:ObjectType-Article-1
SourceType-Scholarly Journals-1
ObjectType-Feature-2
content type line 14
ISSN:1609-4840
1609-9389
DOI:10.1515/cmam-2018-0187