Adaptive Finite Element Simulation of Double-Diffusive Convection

Jorge Milhazes and Pedro J. Coelho ()
Additional contact information
Jorge Milhazes: Romax Technology Ltd., Ergo House, Mere Way, Ruddington Fields Business Park, Nottingham NG11 6JS, UK
Pedro J. Coelho: IDMEC, Mechanical Engineering Department, Instituto Superior Técnico, Universidade de Lisboa, Av. Rovisco Pais, 1049-001 Lisboa, Portugal

Energies, 2023, vol. 16, issue 4, 1-19

Abstract: Double-diffusive convection plays an important role in many physical phenomena of practical importance. However, the numerical simulation of these phenomena is challenging since fine meshes are often required to capture the flow physics. Hence, several different numerical methods have been employed in the past. This work reports the development and application of an adaptive finite element method for the simulation of these phenomena, thereby avoiding the need for the use of very fine meshes over the whole domain. The weak formulation of the conservation equations for mass, momentum, energy and species concentration is used. The Boussinesq approximation relates the density of the fluid to the temperature and/or the species concentration. A second-order backward difference method is used for time discretization and the Galerkin method is employed for spatial discretization. Both adaptive time step and grid refinement techniques are employed, and the code is parallelized using MPI. Three different stabilization methods of the convective-diffusion equations are compared; namely, the streamline upwind Petrov–Galerkin (SUPG) method, and two modified methods aimed at diminishing spurious oscillations that include an artificial diffusion term. This diffusion term may be either isotropic or orthogonal to the streamlines. The addition of artificial isotropic diffusion to the SUPG method provides enhanced stability. The method is applied to double-diffusive finger convection in a sucrose-salt aqueous mixture and a stratified salt solution heated from below. The method accurately reproduces the experimentally observed temporal evolution of the salt fingers in the former case and the location of the interfaces between convective and non-convective zones in the latter.

Keywords: finite elements; double-diffusion; finger convection; stabilization methods; adaptive methods (search for similar items in EconPapers)
JEL-codes: Q Q0 Q4 Q40 Q41 Q42 Q43 Q47 Q48 Q49 (search for similar items in EconPapers)
Date: 2023
References: View complete reference list from CitEc
Citations:

Downloads: (external link)
https://www.mdpi.com/1996-1073/16/4/2010/pdf (application/pdf)
https://www.mdpi.com/1996-1073/16/4/2010/ (text/html)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:gam:jeners:v:16:y:2023:i:4:p:2010-:d:1072277

Access Statistics for this article

Energies is currently edited by Ms. Agatha Cao

More articles in Energies from MDPI
Bibliographic data for series maintained by MDPI Indexing Manager ().