Testing density-dependent groundwater models: Two-dimensional steady state unstable convection in infinite, finite and inclined porous layers

Weatherill, D.; Simmons, C. T.; Voss, C. I.; Robinson, N. I.

Abstract:

This study proposes the use of several problems of unstable steady state convection with variable fluid density in a porous layer of infinite horizontal extent as two-dimensional (2-D) test cases for density-dependent groundwater flow and solute transport simulators. Unlike existing density-dependent model benchmarks, these problems have well-defined stability criteria that are determined analytically. These analytical stability indicators can be compared with numerical model results to test the ability of a code to accurately simulate buoyancy driven flow and diffusion. The basic analytical solution is for a horizontally infinite fluid-filled porous layer in which fluid density decreases with depth. The proposed test problems include unstable convection in an infinite horizontal box, in a finite horizontal box, and in an infinite inclined box. A dimensionless Rayleigh number incorporating properties of the fluid and the porous media determines the stability of the layer in each case. Testing the ability of numerical codes to match both the critical Rayleigh number at which convection occurs and the wavelength of convection cells is an addition to the benchmark problems currently in use. The proposed test problems are modelled in 2-D using the SUTRA [SUTRA-A model for saturated-unsaturated variable-density ground-water flow with solute or energy transport. US Geological Survey Water-Resources Investigations Report, 02-4231, 2002. 250 p] density-dependent groundwater flow and solute transport code. For the case of an infinite horizontal box, SUTRA results show a distinct change from stable to unstable behaviour around the theoretical critical Rayleigh number of 4??2 and the simulated wavelength of unstable convection agrees with that predicted by the analytical solution. The effects of finite layer aspect ratio and inclination on stability indicators are also tested and numerical results are in excellent agreement with theoretical stability criteria and with numerical results previously reported in traditional fluid mechanics literature. ?? 2004 Elsevier Ltd. All rights reserved.

Genre:

Article

ProdID:

70027189

Citation Author:

Weatherill, D.; Simmons, C. T.; Voss, C. I.; Robinson, N. I.

Citation Contributing Office:

Citation Datum:

Citation Day:

Citation Edition:

Citation Editor:

Citation End Page:

562

Citation Issue:

Citation Keywords:

Citation Language:

English

Citation Larger Work Title:

Advances in Water Resources

Citation LatN:

Citation LatS:

Citation LonE:

Citation LonW:

Citation Month:

Citation No Pagination:

Citation Number Of Pages:

Citation Online Only Flag:

Citation Phsyical Description:

Citation Projection:

Citation Public Comments:

Citation Publisher:

Citation Series:

Citation Series Code:

Citation Series Number:

Citation Search Results Text:

Testing density-dependent groundwater models: Two-dimensional steady state unstable convection in infinite, finite and inclined porous layers; 2004; Article; Journal; Advances in Water Resources; Weatherill, D.; Simmons, C. T.; Voss, C. I.; Robinson, N. I.