@article{farthing_kees_coffey_kelley_miller_2003, title={Efficient steady-state solution techniques for variably saturated groundwater flow}, volume={26}, ISSN={["0309-1708"]}, DOI={10.1016/S0309-1708(03)00076-9}, abstractNote={We consider the simulation of steady-state variably saturated groundwater flow using Richards’ equation (RE). The difficulties associated with solving RE numerically are well known. Most discretization approaches for RE lead to nonlinear systems that are large and difficult to solve. The solution of nonlinear systems for steady-state problems can be particularly challenging, since a good initial guess for the steady-state solution is often hard to obtain, and the resulting linear systems may be poorly scaled. Common approaches like Picard iteration or variations of Newton’s method have their advantages but perform poorly with standard globalization techniques under certain conditions. Pseudo-transient continuation has been used in computational fluid dynamics for some time to obtain steady-state solutions for problems in which Newton’s method with standard line-search strategies fails. Here, we examine the use of pseudo-transient continuation as well as Newton’s method combined with standard globalization techniques for steady-state problems in heterogeneous domains. We investigate the methods’ performance with direct and preconditioned Krylov iterative linear solvers. We then make recommendations for robust and efficient approaches to obtain steady-state solutions for RE under a range of conditions.}, number={8}, journal={ADVANCES IN WATER RESOURCES}, author={Farthing, MW and Kees, CE and Coffey, TS and Kelley, CT and Miller, CT}, year={2003}, month={Aug}, pages={833–849} }
 @article{coffey_mcmullan_kelley_mcrae_2003, title={Globally convergent algorithms for nonsmooth nonlinear equations in computational fluid dynamics}, volume={152}, ISSN={["0377-0427"]}, DOI={10.1016/S0377-0427(02)00697-0}, abstractNote={In this paper we report on a computational study in which a nonsmooth discretization of the Euler equations for flow in a nozzle is solved with splitting method which is in turn globalized with the method of pseudo-transient continuation.}, number={1-2}, journal={JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS}, author={Coffey, T and McMullan, RJ and Kelley, CT and McRae, DS}, year={2003}, month={Mar}, pages={69–81} }