Numerical simulation of double‐diffusive finger convection

Water Resources Research
By: , and 



A hybrid finite element, integrated finite difference numerical model is developed for the simulation of double‐diffusive and multicomponent flow in two and three dimensions. The model is based on a multidimensional, density‐dependent, saturated‐unsaturated transport model (SUTRA), which uses one governing equation for fluid flow and another for solute transport. The solute‐transport equation is applied sequentially to each simulated species. Density coupling of the flow and solute‐transport equations is accounted for and handled using a sequential implicit Picard iterative scheme. High‐resolution data from a double‐diffusive Hele‐Shaw experiment, initially in a density‐stable configuration, is used to verify the numerical model. The temporal and spatial evolution of simulated double‐diffusive convection is in good agreement with experimental results. Numerical results are very sensitive to discretization and correspond closest to experimental results when element sizes adequately define the spatial resolution of observed fingering. Numerical results also indicate that differences in the molecular diffusivity of sodium chloride and the dye used to visualize experimental sodium chloride concentrations are significant and cause inaccurate mapping of sodium chloride concentrations by the dye, especially at late times. As a result of reduced diffusion, simulated dye fingers are better defined than simulated sodium chloride fingers and exhibit more vertical mass transfer.

Additional publication details

Publication type Article
Publication Subtype Journal Article
Title Numerical simulation of double‐diffusive finger convection
Series title Water Resources Research
DOI 10.1029/2003WR002777
Volume 41
Issue 1
Year Published 2005
Language English
Publisher American Geophysical Union
Contributing office(s) Toxic Substances Hydrology Program
Description W01019; 16 p.
Google Analytic Metrics Metrics page