Preconditioned Conjugate-Residual Solvers for Helmholtz Equations in Nonhydrostatic Models

Numerical integration of the compressible nonhydrostatic equations using semi-implicit techniques is complicated by the need to solve a Helmholtz equation at each time step. The authors present an accurate and efficient technique for solving the Helmholtz equation using a conjugate-residual (CR) algorithm that is accelerated by ADI preconditioners. These preconditioned CR solvers possess four distinct advantages over most other solvers that have been used with the Helmholtz equations that arise in compressible nonhydrostatic semi-implicit atmospheric models: the preconditioned CR methods 1) can solve Helmholtz equations containing variable coefficients, alleviating the need to prescribe a reference state in order to simplify the elliptic problem; 2) transparently include the cross-derivative terms arising from terrain transformations; 3) are efficient and accurate for nonhydrostatic models used across a broad range of scales, from cloud scales to synoptic-global scales; and 4) are easy to formulate and program. These features of the CR solver allow semi-implicit formulations that are unconstrained by the form of the Helmholtz equations, and the authors propose a formulation that is more consistent than those most often used in that it includes implicit treatment of all terms associated with the pressure gradients and divergence. This formulation is stable for nonhydrostatic-scale simulations involving steep terrain, whereas the more common semi-implicit formulation is not. The ADI preconditioners are presented for use in simulations of both hydrostatic and nonhydrostatic scale flows. These simulations demonstrate the efficiency and accuracy of the preconditioned CR method and the overall stability of the model formulation. The simulations also suggest a general convergence criteria for the iterative algorithm in terms of the solution divergence.

[1]  W. Skamarock,et al.  Adaptive Grid Refinement for Two-Dimensional and Three-Dimensional Nonhydrostatic Atmospheric Flow , 1993 .

[2]  Jeffrey S. Scroggs,et al.  A global nonhydrostatic semi-Lagrangian atmospheric model without orography , 1995 .

[3]  H. Kapitza,et al.  The non-hydrostatic mesoscale model GESIMA. I : Dynamical equations and tests , 1992 .

[4]  Len G. Margolin,et al.  On Forward-in-Time Differencing for Fluids: Stopping Criteria for Iterative Solutions of Anelastic Pressure Equations , 1997 .

[5]  Josef Stoer,et al.  Solution of Large Linear Systems of Equations by Conjugate Gradient Type Methods , 1982, ISMP.

[6]  M. Cullen.  A test of a semi‐implicit integration technique for a fully compressible non‐hydrostatic model , 1990 .

[7]  P. Roache.  QUANTIFICATION OF UNCERTAINTY IN COMPUTATIONAL FLUID DYNAMICS , 1997 .

[8]  J. J. Douglas.  Alternating direction methods for three space variables , 1962 .

[9]  John Marshall,et al.  Application of a parallel Navier-Stokes model to ocean circulation , 1996 .

[10]  Chris Snyder,et al.  An Analysis of Frontogenesis in Numerical Simulations of Baroclinic Waves , 1994 .

[11]  E. Stiefel,et al.  Relaxationsmethoden bester Strategie zur Lösung linearer Gleichungssysteme , 1955 .

[12]  D. Durran,et al.  A Compressible Model for the Simulation of Moist Mountain Waves , 1983 .

[13]  Monique Tanguay,et al.  A Semi-implicit Send-Lagrangian Fully Compressible Regional Forecast Model , 1990 .

[14]  B. Golding.  An efficient non-hydrostatic forecast model , 1992 .

[15]  William C. Skamarock,et al.  Efficiency and Accuracy of the Klemp-Wilhelmson Time-Splitting Technique , 1994 .

[16]  O. Axelsson.  Iterative solution methods , 1995 .

[17]  Richard C. J. Somerville,et al.  On the use of a coordinate transformation for the solution of the Navier-Stokes equations , 1975 .

[18]  P. Bernardet,et al.  The Pressure Term in the Anelastic Model: A Symmetric Elliptic Solver for an Arakawa C Grid in Generalized Coordinates , 1995 .

[19]  J. Klemp,et al.  The Simulation of Three-Dimensional Convective Storm Dynamics , 1978 .

[20]  S. Eisenstat,et al.  Variational Iterative Methods for Nonsymmetric Systems of Linear Equations , 1983 .

[21]  Jean Côté,et al.  A Generalized Family of Schemes that Eliminate the Spurious Resonant Response of Semi-Lagrangian Schemes to Orographic Forcing , 1995 .

[22]  M. Hestenes,et al.  Methods of conjugate gradients for solving linear systems , 1952 .

[23]  Motohki Ikawa,et al.  Comparison of Some Schemes for Nonhydrostatic Models with Orography , 1988 .

[24]  W. Skamarock,et al.  The stability of time-split numerical methods for the hydrostatic and the nonhydrostatic elastic equations , 1992 .