A Fourier Series Solution for Transient Three‐Dimensional Thermohaline Convection in Porous Enclosures. Issue 11 (6th November 2020)

• Record Type: Journal Article
• Authors: Tabrizinejadas, Sara; Fahs, Marwan; Ataie‐Ashtiani, Behzad; Simmons, Craig T.; di Chiara Roupert, Raphaël; Younes, Anis
• Is Part Of: Water resources research. Volume 56:Issue 11(2020)
• Page Start: n/a
• Page End: n/a
• Publication Date: 2020-11-06
• Subjects: thermohaline convection -- 3‐D semianalytical solution -- crossed thermal and salinity gradients -- Fourier series -- benchmarking -- COMSOL multiphysics; Hydrology -- Periodicals; 333.91
• DOI: 10.1029/2020WR028111
• Languages: English
• ISSNs: 0043-1397

Abstract

Abstract: Thermohaline convection (THC) in porous media is frequently investigated using the problem of porous enclosure. Most of the existing modeling‐based studies are limited to 2‐D simulations, because 2‐D assumption is widely used to deal with computational requirement of 3‐D numerical solutions. Analytical solutions serve as an alternative to deal with computational requirement of numerical solutions. Existing analytical solutions of THC are mostly limited to 2‐D and also under steady‐state regime. In this work, we develop a meshless 3‐D semianalytical solution for the problem of THC in a porous box under crossed thermal and solute gradients, for both steady state and transient regimes. The semianalytical solution is developed using the Fourier series (FS) method applied to the vector potential form of the governing equations. The extension to transient solutions represents an important technical feature of this work, as the applications of the FS method to density‐driven problems have been limited to steady‐state conditions. The FS solution is validated against a finite element solution obtained using COMSOL Multiphysics. Numerical experiments show the worthiness of the developed FS solution as a benchmark because it clearly allows making distinction between different numerical techniques. The effects of governing parameters on three‐dimensional THC have not been investigated previously. We perform a detailed parameter sensitivity analysis to address this gap. A vortex convective flow is observed and the orientation and intensity of the flow is sensitive to the gravity number. The increase in the temperature gradient reduces the salinity flux. Key Points: Existing numerical and analytical studies of thermohaline convection (THC) are limited to 2‐D and steady‐state conditions A semianalytical solution is developed for THC in 3‐D box with crossed heat‐salinity gradients at both transient and steady‐state regimes Observations: a single vortex flow with orientation dependent on the Gravity number and different transient and steady‐state behaviors