Interaction of temperature- and salinity-driven natural convection in homogeneous porous media

<p>Onset of thermal and haline convection was studied separately by Lapwood (1948) and Wooding (1956) in theoretical models using analytical methods. They established that the buoyancy force caused by difference in temperature (Δ<em>T</em>) or concentration (Δ<em>c</em>) can induce natural convection depending on the model properties (e.g. geometry, permeability, etc.). In the course of further numerical simulations, the thermal (<em>Ra<sub>T</sub></em>) and the haline Rayleigh number (<em>Ra<sub>H</sub></em>) proved itself useful to characterise the type, the intensity and the form of the natural convection (e.g. Diersch and Kolditz, 2002). The main purpose of our study was to examine numerically the combined effect of temperature- and salinity-driven natural convection in a two-dimensional homogeneous porous medium.</p><p>Two-dimensional finite element base model was set up in agreement with the Elder problem (Wooding, 1956) in order to verify the numerical calculation. First, it was established that (1) the critical Rayleigh numbers are mathematically equivalent in the two separated cases (<em>Ra</em><sub><em>Tcr</em></sub>=<em>Ra<sub>Hcr</sub></em>=4π<sup>2</sup>), and (2) time-dependent thermal or haline natural convections evolve, when the Rayleigh number lies within the range of 300–600. Numerical simulations were accomplished to investigate the interaction of the temperature- and salinity-driven natural convection. Non-dimensional thermal expansion and haline concentration were increased from αΔ<em>T</em>=0.01 to 1 and from βΔ<em>c</em>=10<sup>-5</sup> to 10<sup>-3</sup>, respectively, while the variation of the Darcy flux, the temperature, the concentration, the Nusselt and the Sherwood numbers was computed. The main points of this study were that (1) how the onset of the thermohaline convection is facilitated by the positive interaction of the thermal and haline effects (<em>Ra<sub>THcr</sub></em>); (2) under what conditions time-dependent flow evolves in the theoretical models; (3) whether a new non-dimensional number can be defined instead of the two separated Rayleigh numbers in order to characterise the behaviour of the thermonaline convection. These simulations draw attention to the importance of understanding the physical background of thermohaline convection, for instance, at the margin of confined and unconfined carbonate systems (e.g. Buda Thermal Karst), or in the case of groundwater flow induced by water pumping/injection of deep geothermal power plants.</p><p>The project was supported by the ÚNKP-19-3 and ÚNKP-19-4 New National Excellence Program of the Ministry for Innovation and Technology, the Hungarian Research Fund (K 129279) and the János Bolyai Scholarship of the Hungarian Academy of Science. This research is a part of a project that has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement No 810980.</p><p>References:</p><p>Diersch, H.-J.G., Kolditz, O. (2002), Variable-density flow and transport in porous media: approaches and challenges. Advances in Water Resources, 25, 899-944.</p><p>Lapwood, E.R. (1948), Convection of a fluid in a porous medium. Mathematical Proceedings of the Cambridge Philosophical Society, 44, 508-521.</p><p>Wooding, R.A. (1956), Steady state free thermal convection of liquid in a saturated permeable medium. Journal of Fluid Mechanics, 2, 273-285.</p>

Causal Heat Conduction Contravening the Fading Memory Paradigm

We propose a causal heat conduction model based on a heat kernel violating the fading memory paradigm. The resulting transport equation produces an equation for the temperature. The model is applied to the discussion of two important issues such as the thermohaline convection and the nuclear burning (in)stability. In both cases, the behaviour of the system appears to be strongly dependent on the transport equation assumed, bringing out the effects of our specific kernel on the final description of these problems. A possible relativistic version of the obtained transport equation is presented.