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.^[2]

Convection in fluids is driven by density variations within them under the influence of gravity. 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. Another example of double diffusion is the formation of false bottoms at the interface of sea ice and under-ice meltwater layers.^[3] A good discussion of many of these processes is in Stewart Turner's monograph "Buoyancy effects in fluids".^[4]

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,^[5] and in the sun (where heat and helium diffuse at differing rates). Sediment can also be thought as having a slow Brownian diffusion rate compared to salt or heat, so double diffusive convection is thought to be important below sediment laden rivers in lakes and the ocean.^[6][7]

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—a density ratio between 0 and 1), 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.^[8] 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.^[9]

Governing equations

The conservation equations for vertical momentum, heat and salinity equations (under Boussinesq's approximation) have the following form for double diffusive salt fingers:^[10]

$${\displaystyle {\nabla }\cdot U=0}$$

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

$${\displaystyle {\frac {\partial T}{\partial t))+U\cdot \nabla T=k_{T}{\nabla }^{2}T}$$

$${\displaystyle {\frac {\partial S}{\partial t))+U\cdot \nabla 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 haline contraction coefficient 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^[11]

$${\displaystyle 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))}.}$$

$${\displaystyle {\nabla }\cdot u=0}$$

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

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

$${\displaystyle {\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

$${\displaystyle R_{\rho }={\frac {\alpha \Delta T}{\beta \Delta S)),Ra_{T}={\frac {g\alpha \Delta TH^{3)){\nu k_{T))},Pr={\frac {\nu }{k_{T))},Sc={\frac {\nu }{k_{S))}.}$$

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,

$${\displaystyle R_{f}={\frac {\alpha T'}{\beta S')).}$$

Applications

Double diffusive convection holds importance in natural processes and engineering applications.^[12][13] The effect of double diffusive convection is not limited to oceanography, also occurring in geology,^[14] astrophysics,^[15] and metallurgy.^[16]