Double diffusive convection

File:Salt fingers.png

Numerical simulations results by Singh and Srinivasan show concentration fields at different Rayleigh numbers for fixed value of R_ρ = 6. The parameters are: (a) Ra_T = 7×10⁸ , t=1.12×10⁻², (b) Ra_T =3.5×10⁸, t=1.12×10⁻², (c) Ra_T =7×10⁶, t=1.31×10⁻², (d) Ra_T=7×10⁵, t=3.69×10⁻². It is seen from the figure that finger characteristics such as width, evolution pattern are a function of Rayleigh numbers.

Double diffusive convection is a fluid dynamics phenomenon that describes a form of convection driven by two different density gradients, which have different rates of diffusion.

Convection in fluids is driven by density variations within them. These density variations may be caused by gradients in the composition of the fluid, or by differences in temperature (through thermal expansion). Thermal and compositional gradients can often diffuse with time, reducing their ability to drive the convection, and requiring that gradients in other regions of the flow exist in order for convection to continue. A common example of double diffusive convection is in oceanography, where heat and salt concentrations exist with different gradients and diffuse at differing rates. An effect that affects both of these variables is the input of cold freshwater from an iceberg.

Double diffusive convection is important in understanding the evolution of a number of systems that have multiple causes for density variations. These include convection in the Earth's oceans (as mentioned above), in magma chambers, and in the sun (where heat and helium diffuse at differing rates).

Two quite different types of fluid motion exist—and therefore are classified accordingly—depending on whether the stable stratification is provided by the density-affecting component with the lowest or the highest molecular diffusivity. If the stratification is provided by the component with the lower molecular diffusivity (for example in case of a stable salt-stratified ocean perturbed by a thermal gradient due to an iceberg), the stratification is called to be of "diffusive" type (see external link below), otherwise it is of "finger" type, occurring frequently in oceanographic studies as salt-fingers. These long fingers of rising and sinking water occur when hot saline water lies over cold fresh water of a higher density. A perturbation to the surface of hot salty water results in an element of hot salty water surrounded by cold fresh water. This element loses its heat more rapidly than its salinity because the diffusion of heat is faster than of salt; this is analogous to the way in which just unstirred coffee goes cold before the sugar has diffused to the top. Because the water becomes cooler but remains salty, it becomes denser than the fluid layer beneath it. This makes the perturbation grow and causes the downward extension of a salt finger. As this finger grows, additional thermal diffusion accelerates this effect.

Role of salt fingers in oceans

Double diffusion convection plays a significant role in upwelling of nutrients and vertical transport of heat and salt in oceans. Salt fingering contributes to vertical mixing in the oceans. Such mixing helps regulate the gradual overturning circulation of the ocean, which control the climate of the earth. Apart from playing an important role in controlling the climate, fingers are responsible for upwelling of nutrients which supports flora and fauna. The most significant aspect of finger convection is that they transport the fluxes of heat and salt vertically, which has been studied extensively during the last five decades.

Governing Equations

The conservation equations for vertical momentum, heat and salinity equations (under Boussinesq’s approximation) have following form for double diffusive salt fingers:

${\nabla}\cdot U =0$

$\frac{\partial U}{\partial t} + U\cdot\nabla U = \nu {\nabla}^2 U - g(\beta \Delta S-\alpha \Delta T)\mathbf{k}$

$\frac{\partial T}{\partial t} + U\cdot\Delta T = k_T {\nabla}^2 T$

$\frac{\partial S}{\partial t} + U\cdot\Delta S = k_s {\nabla}^2 s$

Where, U and W are velocity components in horizontal (x axis) and vertical (z axis) direction; k is the unit vector in the Z-direction, k_T is molecular diffusivity of heat, k_s is molecular diffusivity of salt, α is coefficient of thermal expansion at constant pressure and salinity and β is the coefficient of saline expansion at constant pressure and temperature. The above set of conservation equations governing the two-dimensional finger-convection system is non-dimensionalised using the following scaling: the depth of the total layer height H is chosen as the characteristic length, velocity (U, W), salinity (S), temperature (T) and time (t) are non-dimensionalised as

$x=\frac{X}{H} , z=\frac{Z}{H}, u=\frac{U}{k_T /H}, w= \frac{W}{k_T /H}, S^*= \frac{S-S_B}{S_T-S_B} , T^*= \frac{T-T_B}{T_T-T_B}, t^* = \frac{t}{H^2 /k_T }.$

Where, (T_T, S_T) and (T_B, S_B) are the temperature and concentration of the top and bottom layers respectively. On introducing the above non-dimensional variables, the above governing equations reduce to the following form:

${\nabla}\cdot u =0$

$\frac{\partial u}{\partial t^*} + u\cdot\nabla u = Pr {\nabla}^2 u - \left[Pr Ra_T (\frac{S^*}{R_\rho}-T^*)\right] \mathbf{k}$

$\frac{\partial T^*}{\partial t^*} + u\cdot\Delta T^* = {\nabla}^2 T^*$

$\frac{\partial S^*}{\partial t^*} + u\cdot\Delta S^* = \frac{Pr}{Sc} {\nabla}^2 S^*$

Where, R_ρ is the density stability ratio, Ra_T is the thermal Rayleigh number, Pr is the Prandtl number, Sc is the Schmidt number which are defined as

$R_\rho = \frac{\alpha \Delta T }{\beta \Delta S}, Ra_T = \frac{g\beta \Delta S H^3}{\nu k_T}, Pr=\frac{\nu}{\alpha}, Sc=\frac{\nu}{\beta} .$

Figure 1(a-d) shows the evolution of salt fingers in heat-salt system for different Rayleigh numbers at a fixed R_ρ. It can be noticed that thin and thick fingers form at different Ra_T. Fingers flux ratio, growth rate, kinetic energy, evolution pattern, finger width etc. are found to be the function of Rayleigh numbers and R_ρ.Where, flux ratio is another important non-dimensional parameter. It is the ratio of heat and salinity fluxes, defined as,

$R_f=\frac{\alpha T'}{\beta S'}.$

Applications

Double diffusion convection holds importance in natural processes and engineering applications. The effect of double diffusive convection is not limited to oceanography but also in various fields like geology, astrophysics, metallurgy.

See also

• Rayleigh number
• Prandtl number

External links

• Oceanic Double-diffusion: Introduction
• Double Diffusion in Oceanography
• Diffusive-mode Double Diffusive Convection, Stability and Density-Driven Flows
• Stockman, H.W; Li, C.; Cooper, C.; 1997. Practical application of lattice-gas and lattice Boltzmann methods to dispersion problems. InterJournal of Complex Systems, manuscript no. 90.