Multi‐resolution simulation of double‐diffusive convection in porous media

Abstract

PurposeThe purpose of this paper is to consider double‐diffusive convection in a heated porous medium saturated with a fluid. Of particular interest is the case where the fluid has a stabilizing concentration gradient and small diffusivity.Design/methodology/approachA fully‐coupled stabilized finite element scheme and adaptive mesh refinement (AMR) methodology are introduced to solve the resulting coupled multiphysics application and resolve fine scale solution features. The code is written on top of the open source finite element library LibMesh, and is suitable for parallel, high‐performance simulations of large‐scale problems.FindingsThe stabilized adaptive finite element scheme is used to compute steady and unsteady onset of convection in a generalized Horton‐Rogers‐Lapwood problem in both two and three‐dimensional domains. A detailed study confirming the applicability of AMR in obtaining the predicted dependence of solutal Nusselt number on Lewis number is given. A semi‐permeable barrier version of the generalized HRL problem is also studied and is believed to present an interesting benchmark for AMR codes owing to the different boundary and internal layers present in the problem. Finally, some representative adaptive results in a complex 3D heated‐pipe geometry are presented.Originality/valueThis work demonstrates the feasibility of stabilized, adaptive finite element schemes for computing simple double‐diffusive flow models, and it represents an easily‐generalizable starting point for more complex calculations since it is based on a highly‐general finite element library. The complementary nature of h‐adaptivity and stabilized finite element techniques for this class of problem is demonstrated using particularly simple error indicators and stabilization parameters. Finally, an interesting double‐diffusive convection benchmark problem having a semi‐permeable barrier is suggested.