Boundary Conditions for Semi-Lagrangian Schemes: Testing Some Alternatives in One-Dimensional Models

A. McDonald Met Éireann, Dublin, Ireland

Search for other papers by A. McDonald in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

The one-dimensional advection equation and the one-dimensional advection adjustment equation with rotation are used as test beds for the design of well-posed boundary conditions for the initial-boundary-value problem using semi-Lagrangian discretization. Three options are found to be stable in experimental tests: trajectory truncation, time interpolation, and a well-posed buffer zone. Stability is proved for all three for the one-dimensional advection equation when linear interpolation is used for the interpolation associated with the semi-Lagrangian discretization.

Corresponding author address: Aidan McDonald, Met Éireann, Glasnevin Hill, Dublin 9, Ireland.

Email: Aidan.McDonald@met.ie

Abstract

The one-dimensional advection equation and the one-dimensional advection adjustment equation with rotation are used as test beds for the design of well-posed boundary conditions for the initial-boundary-value problem using semi-Lagrangian discretization. Three options are found to be stable in experimental tests: trajectory truncation, time interpolation, and a well-posed buffer zone. Stability is proved for all three for the one-dimensional advection equation when linear interpolation is used for the interpolation associated with the semi-Lagrangian discretization.

Corresponding author address: Aidan McDonald, Met Éireann, Glasnevin Hill, Dublin 9, Ireland.

Email: Aidan.McDonald@met.ie

Save
  • Bates, J. R., and A. McDonald, 1982: Multiply-upstream, semi-Lagrangian advective schemes: Analysis and application to a multi-level primitive equation model. Mon. Wea. Rev.,110, 1831–1842.

  • Cats, G. J., and O. Ã…kesson, 1983: An investigation into a marked difference between two successive ECMWF forecasts of September, 1982. Beitr. Phys. Atmos.,56, 440–451.

  • Davies, H. C., 1976: A lateral boundary formulation for multi-level prediction models. Quart. J. Roy. Meteor. Soc.,102, 405–418.

  • Durran, D. R., 1999: Numerical Methods for Wave Equations in Geophysical Fluid Dynamics. Springer-Verlag, 465 pp.

  • Elvius, T., and A. Sundström, 1973: Computationally efficient schemes and boundary conditions for a fine mesh barotropic model based on the shallow water equations. Tellus,25, 132–156.

  • Engquist, B., and A. Majda, 1977: Absorbing boundary conditions for the numerical simulation of waves. Math. Comput.,31, 629–651.

  • Falcone, M., and R. Ferretti, 1998: Convergence analysis for a class of high-order semi-Lagrangian advection schemes. SIAM J. Numer. Anal.,35, 909–940.

  • Gill, A. E., 1982: Atmosphere–Ocean Dynamics. Academic Press, 662 pp.

  • Gustafsson, B., H.-O. Kreiss, and A. Sundström, 1972: Stability theory of difference approximations for mixed initial boundary value problems. Math. Comput.,26, 649–686.

  • Gustafsson, N., E. Källén, and S. Thorsteinsson, 1998: Sensitivity of forecast errors to initial and lateral boundary conditions. Tellus,50A, 167–185.

  • Le Dimet, F. X., and O. Talagrand, 1986: Variational algorithms for analysis and assimilation of meteorological observations: Theoretical aspects. Tellus,38, 97–110.

  • McDonald, A., 1984: Accuracy of multiply-upstream, semi-Lagrangian advective schemes. Mon. Wea. Rev.,112, 1267–1275.

  • ——, 1998: Conservation of mass in HIRLAM. HIRLAM Newsletter, Vol. 32, 39–44. [Available from A. McDonald, Irish Meteorological Service, Glasnevin Hill, Dublin 9, Ireland.].

  • ——, 1999a: A review of lateral boundary conditions for limited area forecast models. Proc. Ind. Nat. Sci. Acad.,65, 91–105.

  • ——, 1999b: Well-posed boundary conditions for semi-Lagrangian schemes: The one-dimensional case. HIRLAM Tech. Rep. 43, HIRLAM Group, 26 pp. [Available from A. McDonald, Irish Meteorological Service, Glasnevin Hill, Dublin 9, Ireland.].

  • ——, 2000: Well-posed boundary conditions for semi-Lagrangian schemes: The one-dimensional case; part II. HIRLAM Tech. Rep. 44, HIRLAM Group, 19 pp. [Available from A. McDonald, Irish Meteorological Service, Glasnevin Hill, Dublin 9, Ireland.].

  • Oliger, J., and A. Sundström, 1978: Theoretical and practical aspects of some initial boundary value problems in fluid dynamics. SIAM J. Appl. Math.,35, 419–446.

  • Robert, A., and E. Yakimiw, 1986: Identification and elimination of an inflow boundary computational solution in limited area model integrations. Atmos.–Ocean,24, 369–385.

  • SkÃ¥lin, R., and I. Lie, 1995: On the stability of semi-Lagrangian advection with interpolation in time. DNMI, 24 pp. [Available from Ivar Lie, DNMI, Box 43, Blindern, N-0313 Oslo, Norway.].

  • Sundström, A., and T. Elvius, 1979: Computational problems related to limited area modelling. Numerical Methods Used in Atmospheric Models, Vol. 2, GARP Publ. Series, No. 17, 379–416.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 249 98 6
PDF Downloads 195 74 3