Numerical Study of the Onset of Double-Diffusive Cellular Convection Due to a Uniform Lateral Heat Flux

1982 ◽  
Vol 104 (4) ◽  
pp. 649-655 ◽  
Author(s):  
S. Takao ◽  
M. Tsuchiya ◽  
U. Narusawa

When a fluid with a vertical solute gradient of (−dS/dy)0 is heated laterally, roll cells start to form at the boundary, developing into a series of convective layers. Numerical experiments were performed to investigate the onset of the abovementioned double-diffusive convection under the application of a uniform lateral heat flux. The paper reports the results and discussion of the following aspects of the stability of double-diffusive convection; (i) the relationship between the critical value, (Ra/Rs)c, above which convection cells form along the vertical wall and the nondimensional slot width, (d/L), (ii) the effect of the Lewis number on (Ra/Rs)c. It was also confirmed that values of (Ra/Rs)c as well as H/L (the nondimensional vertical size of incipient cells) obtained in this numerical experiment for wide slot widths (d/L>∼30), agreed well with those obtained previously by physical experiments.

2017 ◽  
Vol 374 ◽  
pp. 1-17 ◽  
Author(s):  
Mohamed Amine Medebber ◽  
Nourddine Retiel

This article reports a numerical study of double-diffusive convection within the annular region of two concentric vertical cylinders. The outer vertical wall is maintained at lower uniform temperature and concentration, while the inner vertical wall is maintained at higher uniform temperature and concentration. The top and bottom horizontal walls are adiabatic and impermeable to mass transfer. The resulting governing equations are solved using a finite volume method. The coupling between the continuity and momentum equations is solved using the SIMPLER algorithm. The compilations have been obtained for Prandtl numbers (Pr) equal to 7.0, and Lewis number (Le) equal to 100. The thermal Rayleigh number (RaT) and height ratio (X) are, respectively, varied in the range 103≤RaT≤106 and 0.0≤X≤1.0. The influence of physical and geometrical parameters on the streamlines, isotherms, isoconcentrations, average Nusselt and Sherwood numbers has been numerically investigated in detail.


2016 ◽  
Vol 802 ◽  
pp. 667-689 ◽  
Author(s):  
Yantao Yang ◽  
Roberto Verzicco ◽  
Detlef Lohse

Direct numerical simulations are conducted for double diffusive convection (DDC) bounded by two parallel plates. The Prandtl numbers, i.e. the ratios between the viscosity and the molecular diffusivities of scalars, are similar to the values of seawater. The DDC flow is driven by an unstable salinity difference (here across the two plates) and stabilized at the same time by a temperature difference. For these conditions the flow can be in the finger regime. We develop scaling laws for three key response parameters of the system: the non-dimensional salinity flux $\mathit{Nu}_{S}$ mainly depends on the salinity Rayleigh number $\mathit{Ra}_{S}$, which measures the strength of the salinity difference and exhibits a very weak dependence on the density ratio $\unicode[STIX]{x1D6EC}$, which is the ratio of the buoyancy forces induced by two scalar differences. The non-dimensional flow velocity $Re$ and the non-dimensional heat flux $\mathit{Nu}_{T}$ are dependent on both $\mathit{Ra}_{S}$ and $\unicode[STIX]{x1D6EC}$. However, the rescaled Reynolds number $Re\unicode[STIX]{x1D6EC}^{\unicode[STIX]{x1D6FC}_{u}^{eff}}$ and the rescaled convective heat flux $(\mathit{Nu}_{T}-1)\unicode[STIX]{x1D6EC}^{\unicode[STIX]{x1D6FC}_{T}^{eff}}$ depend only on $\mathit{Ra}_{S}$. The two exponents are dependent on the fluid properties and are determined from the numerical results as $\unicode[STIX]{x1D6FC}_{u}^{eff}=0.25\pm 0.02$ and $\unicode[STIX]{x1D6FC}_{T}^{eff}=0.75\pm 0.03$. Moreover, the behaviours of $\mathit{Nu}_{S}$ and $Re\unicode[STIX]{x1D6EC}^{\unicode[STIX]{x1D6FC}_{u}^{eff}}$ agree with the predictions of the Grossmann–Lohse theory which was originally developed for the Rayleigh–Bénard flow. The non-dimensional salt-finger width and the thickness of the velocity boundary layers, after being rescaled by $\unicode[STIX]{x1D6EC}^{\unicode[STIX]{x1D6FC}_{u}^{eff}/2}$, collapse and obey a similar power-law scaling relation with $\mathit{Ra}_{S}$. When $\mathit{Ra}_{S}$ is large enough, salt fingers do not extend from one plate to the other and horizontal zonal flows emerge in the bulk region. We then show that the current scaling strategy can be successfully applied to the experimental results of a heat–copper–ion system (Hage & Tilgner, Phys. Fluids, vol. 22, 2010, 076603). The fluid has different properties and the exponent $\unicode[STIX]{x1D6FC}_{u}^{eff}$ takes a different value $0.54\pm 0.10$.


Sign in / Sign up

Export Citation Format

Share Document