Operator-split Runge-Kutta-Rosenbrock methods for nonhydrostatic atmospheric model

Monthly Weather Review

Volumn 140, Issue 4, 2012, Pages 1257-1284

  Ullrich, Paul ^a   Jablonowski, Christiane ^a  

^a UNIVERSITY OF MICHIGAN   (United States)

Author keywords

Channel flows; Nonhydrostatic models; Numerical analysis modeling; Numerical weather prediction forecasting; Vertical coordinates

Indexed keywords

A-STABLE; ATMOSPHERIC FLOWS; ATMOSPHERIC MODEL; CARTESIAN GEOMETRY; FINITE-VOLUME; GOVERNING EQUATIONS; GRID-SPACINGS; HIGH-ORDER; LARGE ASPECT RATIO; METHOD OF LINES; MODEL FORMULATION; NONHYDROSTATIC MODEL; NONLINEAR ORDINARY DIFFERENTIAL EQUATION; NUMERICAL WEATHER PREDICTION; ROSENBROCK METHODS; RUNGE-KUTTA FORMULA; SEMI-IMPLICIT METHODS; SPATIAL DISCRETIZATIONS; THIRD-ORDER; TIME STEP; TIME STEPPING METHOD; TIME-STEPPING; VERTICAL COMPONENT; VERTICAL COORDINATES; VERTICAL GRIDS; WAVE SPEED;

ASPECT RATIO; ATMOSPHERIC MOVEMENTS; CHANNEL FLOW; EULER EQUATIONS; NONLINEAR EQUATIONS; ORDINARY DIFFERENTIAL EQUATIONS;

RUNGE KUTTA METHODS;

ATMOSPHERIC MODELING; CHANNEL FLOW; CLIMATE PREDICTION; EULERIAN ANALYSIS; NUMERICAL MODEL; WEATHER FORECASTING;

EID: 84859773201     PISSN: 00270644     EISSN: None     Source Type: Journal    
DOI: 10.1175/MWR-D-10-05073.1     Document Type: Article

References (46)

1. Adcroft, A., C. Hill, and J. Marshall, 1997: Representation of topography by shaved cells in a height coordinate ocean model. Mon. Wea. Rev., 125, 2293-2315.
2. Ahmad, N., and J. Lindeman, 2007: Euler solutions using fluxbased wave decomposition. Int. J. Numer. Methods Fluids, 54, 47-72, doi:10.1002/fld.1392.
3. Arakawa, A., and C. S. Konor, 2009: Unification of the anelastic and quasi-hydrostatic systems of equations. Mon. Wea. Rev., 137, 710-726.
4. Ascher, U. M., S. J. Ruuth, and R. J. Spiteri, 1997: Implicit-explicit Runge-Kutta methods for time-dependent partial differential equations. Appl. Numer. Math., 25 (2-3), 151-167, doi:10.1016/ S0168-9274(97)00056-1.
5. Barad, M., and P. Colella, 2005: A fourth-order accurate local refinement method for Poisson's equation. J. Comput. Phys., 209, 1-18, doi:10.1016/j.jcp.2005.02.027.
6. Bonaventura, L., 2000: A semi-implicit semi-Lagrangian scheme using the height coordinate for a nonhydrostatic and fully elastic model of atmospheric flows. J. Comput. Phys., 158, 186-213, doi:10.1006/jcph.1999.6414.
7. Crouzeix, M., 1980: Une méthode multipas implicite-explicite pour l'approximation des équations d'évolution paraboliques. Numerische Math., 35, 257-276, doi:10.1007/BF01396412.
8. Davies, T., A. Staniforth, N. Wood, and J. Thuburn, 2003: Validity of anelastic and other equation sets as inferred from normalmode analysis. Quart. J. Roy. Meteor. Soc., 129, 2761-2775, doi:10.1256/qj.02.1951.
9. Dudhia, J., 1993: A nonhydrostatic version of the Penn State NCAR MesoscaleModel:Validation tests and simulation of anAtlantic cyclone and cold front. Mon. Wea. Rev., 121, 1493-1513.
10. Durran, D. R., 1989: Improving the anelastic approximation. J. Atmos. Sci., 46, 1453-1461.
11. Durran, D. R., 1999: Numerical Methods for Wave Equations in Geophysical Fluid Dynamics. Springer-Verlag, 489 pp.
12. Gadd, A. J., 1978: A split explicit integration scheme for numerical weather prediction. Quart. J. Roy. Meteor. Soc., 104, 569-582, doi:10.1002/qj.49710444103.
13. Gal-Chen, T., and R. C. J. Somerville, 1975: On the use of a coordinate transformation for the solution of the Navier-Stokes equations. J. Comput. Phys., 17, 209-228, doi:10.1016/0021-9991(75)90037-6.
14. Giraldo, F. X., 2005: Semi-implicit time-integrators for a scalable spectral element atmospheric model. Quart. J. Roy. Meteor. Soc., 131, 2431-2454, doi:10.1256/qj.03.218.
15. Giraldo, F. X., and M. Restelli, 2008: A study of spectral element and discontinuous Galerkin methods for the Navier-Stokes equations in nonhydrostatic mesoscale atmo spheric modeling: Equation sets and test cases. J. Comput. Phys., 227, 3849-3877, doi:10.1016/j.jcp.2007.12.009.
16. Gottlieb, S., C.-W. Shu, and E. Tadmor, 2001: Strong stabilitypreserving high-order time discretization methods. SIAM Rev., 43, 89-112, doi:10.1137/S003614450036757X.
17. Jablonowski, C., and D. L. Williamson, 2006:Abaroclinic instabilitiy test case for atmospheric model dynamical cores. Quart. J. Roy. Meteor. Soc., 132, 2943-2975.
18. Jameson, A., W. Schmidt, and E. Turkel, 1981: Numerical solution of the Euler equations by finite volume methods using Runge Kutta time stepping schemes. Proc. 14th Fluid and Plasma Dynamics Conf., Palo Alto, CA, AIAA, 15 pp.
19. Jebens, S., O. Knoth, and R. Weiner, 2011: Partially implicit peer methods for the compressible Euler equations. J. Comput. Phys., 230, 4955-4974, doi:10.1016/j.jcp.2011.03.015.
20. Klein,R.,U.Achatz,D.Bresch,O.M.Knio, andP. K.Smolarkiewicz, 2010: Regime of validity of soundproof atmospheric flowmodels. J. Atmos. Sci., 67, 3226-3237.
21. Klemp, J. B., W. C. Skamarock, and O. Fuhrer, 2003: Numerical consistency of metric terms in terrain-following coordinates. Mon. Wea. Rev., 131, 1229-1239.
22. Kwizak,M., and A. J. Robert, 1971: A semi-implicit scheme for grid point atmospheric models of the primitive equations. Mon. Wea. Rev., 99, 32-36.
23. Lanser, D., J. G. Blom, and J. G. Verwer, 2001: Time integration of the shallow water equations in spherical geometry. J. Comput. Phys., 171, 373-393, doi:10.1006/jcph.2001.6802.
24. Lin, S., 2004: A "vertically Lagrangian" finite-volume dynamical core for global models. Mon. Wea. Rev., 132, 2293-2307.
25. Liou, M.-S., 2006: A sequel to AUSM. Part II: AUSM1-up for all speeds. J. Comput. Phys., 214, 137-170, doi:10.1016/j.jcp.2005. 09.020.
26. McCorquodale, P., and P. Colella, 2011: A high-order finite-volume method for conservation laws on locally refined grids. Commun. Appl. Math. Comput. Sci., 6, 1-25.
27. Norman, M. R., R. D. Nair, and F. H. M. Semazzi, 2011: A low communication and large time step explicit finite-volume solver for non-hydrostatic atmospheric dynamics. J. Comput. Phys., 230, 1567-1584, doi:10.1016/j.jcp.2010.11.022.
28. Nørsett, S. P., and A. Wolfbrandt, 1979: Order conditions for Rosenbrock type methods. Numerische Math., 32, 1-15, doi:10.1007/BF01397646.
29. Ogura, Y., and N. A. Phillips, 1962: Scale analysis of deep and shallow convection in the atmosphere. J. Atmos. Sci., 19, 173-179.
30. Restelli, M., and F. X. Giraldo, 2009:A conservative discontinuous Galerkin semi-implicit formulation for the Navier-Stokes equations in nonhydrostatic mesoscale modeling. SIAM J. Sci. Comput., 31, 2231-2257, doi:10.1137/070708470.
31. Robert, A., 1993: Bubble convection experiments with a semiimplicit formulation of the Euler equations. J. Atmos. Sci., 50, 1865-1873.
32. Roe, P. L., 1981: Approximate Riemann solvers, parameter vectors, and difference schemes. J. Comput. Phys., 43, 357-372, doi:10.1016/0021-9991(81)90128-5.
33. Rosenbrock, H., 1963: Some general implicit processes for the numerical solution of differential equations. Comput. J., 5, 329-330.
34. Rusanov, V., 1961: Calculation of intersection of non-steady shock waves with obstacles. J. Comput. Math. Phys. USSR, 1, 267-279.
35. Sandu, A., J. G. Verwer, J. G. Blom, E. J. Spee, G. R. Carmichael, and F. A. Potra, 1997: Benchmarking stiff ode solvers for atmospheric chemistry problems II: Rosenbrock solvers. Atmos. Environ., 31, 3459-3472, doi:10.1016/S1352-2310(97)83212-8.
36. Schär, C., D. Leuenberger, O. Fuhrer, D. Lüthi, and C. Girard, 2002: A new terrain-following vertical coordinate formulation for atmospheric predictionmodels. Mon.Wea. Rev., 130, 2459-2480.
37. Schiesser, W. E., 1991: The Numerical Method of Lines: Integration of Partial Differential Equations. Academic Press, 326 pp.
38. Schiesser, W. E., and G. W. Griffiths, 2009: A Compendium of Partial Differential Equation Models: Method of Lines Analysis with Matlab. Cambridge University Press, 490 pp.
39. Skamarock, W. C., and J. B. Klemp, 2008: A time-split nonhydrostatic atmospheric model for weather research and forecasting applications. J. Comput. Phys., 227, 3465-3485, doi:10.1016/j.jcp.2007.01.037.
40. St-Cyr, A., and D. Neckels, 2009: A fully implicit Jacobian-free high-order discontinuousGalerkinmesoscale flow solver. Proc. Ninth Int. Conf. on Computational Science, Baton Rouge, LA, Intel, 243-252.
41. Straka, J. M., R. B. Wilhelmson, L. J. Wicker, J. R. Anderson, and K. K. Droegemeier, 1993: Numerical solutions of a nonlinear density current: A benchmark solution and comparisons. Int. J. Numer. Methods Fluids, 17, 1-22, doi:10.1002/ fld.1650170103.
42. Taylor, M. A., J. Edwards, and A. St Cyr, 2008: Petascale atmospheric models for the Community Climate System Model: New developments and evaluation of scalable dynamical cores. J. Phys. Conf. Ser., 125, 012023, doi:10.1088/1742-6596/ 125/1/012023.
43. Ullrich, P. A., C. J. Jablonowski, and B. L. van Leer, 2010: Highorder finite-volume models for the shallow-water equations on the sphere. J. Comput. Phys., 229, 6104-6134.
44. Varah, J. M., 1980: Stability restrictions on second order, three level finite difference schemes for parabolic equations. SIAM J. Numer. Anal., 17, 300-309, doi:10.1137/0717025.
45. Verwer, J. G., E. J. Spee, J. G. Blom, and W. Hundsdorfer, 1999: A second-order Rosenbrock method applied to photochemical dispersion problems. SIAM J. Sci. Comput., 20, 1456-1480, doi:10.1137/S1064827597326651.
46. Williamson, D. L., 2008: Equivalent finite volume and Eulerian spectral transform horizontal resolutions established from aqua-planet simulations. Tellus, 60A, 839-847, doi:10.1111/ j.1600-0870.2008.00340.x.