Premium
Insights from a pseudospectral approach to the Elder problem
Author(s) -
van Reeuwijk Maarten,
Mathias Simon A.,
Simmons Craig T.,
Ward James D.
Publication year - 2009
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/2008wr007421
Subject(s) - discretization , steady state (chemistry) , benchmark (surveying) , mathematics , bifurcation , pseudo spectral method , mathematical optimization , computer science , mathematical analysis , physics , geology , nonlinear system , fourier transform , fourier analysis , chemistry , geodesy , quantum mechanics
The aim of this paper is to clarify and circumvent the issue of multiple steady state solutions in the Elder problem. A pseudospectral method is used to avoid numerical error associated with spatial discretization. The pseudospectral method is verified by comparison to an analytical solution at Rayleigh number, Ra = 0, and by reproducing the three stable steady state solutions that are known to exist at Ra = 400. A bifurcation diagram for 0 < Ra < 400, which is free of discretization error, confirms that multiple steady states are indeed an intrinsic characteristic of the Elder problem. The existence of multiple steady states makes the Ra = 400 Elder problem less suitable for benchmarking numerical models. To avoid the multiple steady states, we propose a benchmark at Ra = 60. The results for this Low Rayleigh Number Elder Problem are presented and compared to simulations with the commercial groundwater modeling package FEFLOW. Correspondence between the pseudospectral model and FEFLOW is excellent.