Editorial Type:
Article
Article Type:
Research Article

A Successive Substitution Method for the Evaluation of Trajectories Approximating the Parcel Path by a Linear Function of Space and Time

Valdir Innocentini Instituto Nacional de Pesquisas Espaciais, Sao Jose dos Campos, Brazil

Search for other papers by Valdir Innocentini in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Jul 1999
DOI:
https://doi.org/10.1175/1520-0493(1999)127<1639:ASSMFT>2.0.CO;2
Page(s):
1639–1650
Restricted access

Abstract

Parcel trajectories can be determined by the time integration of the velocity using the Picard successive substitution procedure. However, this scheme can rarely be applied in numerical problems because it implies (in general) an increasing difficulty in solving the analytical solution in each iteration. This strong restriction to its practical use can be overcome by approximating the velocity to a simple function where the analytical integration can be performed as the iteration advances.

The usual procedure adopted in meteorological research is to compute the trajectories using the Petterssen scheme. In this work it is shown that approximating the velocity by a linear time function along the trajectory in the Picard method yields a numerical algorithm formally identical to the Petterssen scheme. An alternative method is proposed following the philosophy of approximating the velocity in the Picard method in order to simplify the integration. In this approach, the parcel path is assumed to be a linear function of space and time, and the velocity along this path at each iteration is computed using an interpolating polynomial. Only one analytical integration is necessary, and the successive integration depends upon the coefficients of the interpolating function adopted for the velocity. For cubic spline, the alternative algorithm can be represented in a suitable vectorial form convenient in designing the numerical code.

The method is tested on idealized situations with high time and space velocity variations using either Lagrangian or Eulerian coordinates. The results for the alternative and Petterssen schemes are similar in the Eulerian coordinate, but the superiority of the alternative scheme is evident in the Lagrangian coordinate.

Corresponding author address: Dr. Valdir Innocentini, Instituto Nacional de Pesquisas Espaciais—DCM, Av. dos Astronautas, 1758, Caixa Postal 515, CEP 12227-010, São Jose dos Campos-SP, Brazil.

Email: valdir@met.inpe.br

Abstract

Parcel trajectories can be determined by the time integration of the velocity using the Picard successive substitution procedure. However, this scheme can rarely be applied in numerical problems because it implies (in general) an increasing difficulty in solving the analytical solution in each iteration. This strong restriction to its practical use can be overcome by approximating the velocity to a simple function where the analytical integration can be performed as the iteration advances.

The usual procedure adopted in meteorological research is to compute the trajectories using the Petterssen scheme. In this work it is shown that approximating the velocity by a linear time function along the trajectory in the Picard method yields a numerical algorithm formally identical to the Petterssen scheme. An alternative method is proposed following the philosophy of approximating the velocity in the Picard method in order to simplify the integration. In this approach, the parcel path is assumed to be a linear function of space and time, and the velocity along this path at each iteration is computed using an interpolating polynomial. Only one analytical integration is necessary, and the successive integration depends upon the coefficients of the interpolating function adopted for the velocity. For cubic spline, the alternative algorithm can be represented in a suitable vectorial form convenient in designing the numerical code.

The method is tested on idealized situations with high time and space velocity variations using either Lagrangian or Eulerian coordinates. The results for the alternative and Petterssen schemes are similar in the Eulerian coordinate, but the superiority of the alternative scheme is evident in the Lagrangian coordinate.

Corresponding author address: Dr. Valdir Innocentini, Instituto Nacional de Pesquisas Espaciais—DCM, Av. dos Astronautas, 1758, Caixa Postal 515, CEP 12227-010, São Jose dos Campos-SP, Brazil.

Email: valdir@met.inpe.br

Save
Cover Monthly Weather Review
Monthly Weather Review

  • Bermejo, R., 1990: On the equivalence of semi-Lagrangian schemes and particle-in-cell finite element methods. Mon. Wea. Rev.,118, 979–987.

  • Boyce, W. E., and R. C. DiPrima, 1977: Elementary Differential Equations. John Wiley and Sons, 451 pp.

  • Coddington, E. A., 1961: An Introduction to Ordinary Differential Equations. Prentice-Hall, 292 pp.

  • Draxler, R. R., 1991: The accuracy of trajectories during ANATEX calculated using dynamic model analyses versus rawinsonde observations. J. Appl. Meteor.,30, 1446–1467.

  • Gear, C. W., 1971: Numerical Initial Value Problems in Ordinary Differential Equations. Prentice-Hall, 253 pp.

  • Innocentini, V., and E. S. Caetano Neto, 1992: A numerical study of the role of humidity in the updraft driven by moist slantwise convection. J. Atmos. Sci.,49, 1092–1106.

  • Jakobs, H. J., H. Fieldmann, H. Hass, and M. Memmesheimer, 1995:The use of nested models for air pollution studies: An application of the FURAD model to a SANA episode. J. Appl. Meteor.,34, 1301–1319.

  • Kreider, D. L., R. G. Kuller, and D. R. Ostberg, 1968: Elementary Differential Equations. Addison-Wesley, 485 pp.

  • Kuo, Y.-H., M. Skumanich, P. L. Haagenson, and J. S. Chang, 1985:The accuracy of trajectory models as revealed by observing system simulation experiments. Mon. Wea. Rev.,113, 1852–1867.

  • ——, R. J. Reed, and S. Low-Nam, 1992: Thermal structure and airflow in a model simulation of an occluded marine cyclone. Mon. Wea. Rev.,120, 2280–2297.

  • Leslie, L. M., and R. J. Purser, 1995: Three-dimensional mass-conserving semi-Lagrangian scheme employing forward trajectories. Mon. Wea. Rev.,123, 2551–2566.

  • Makar, P. A., and S. R. Karpik, 1996: Basis-spline interpolation on the sphere: Applications to semi-Lagrangian advection. Mon. Wea. Rev.,124, 182–199.

  • Maryon, R. H., and C. C. Heasman, 1988: The accuracy of plume trajectories forecast using the U. K. Meteorological Office operational forecasting models and their sensitivity to calculation schemes. Atmos. Environ.,22, 259–272.

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

  • Petterssen, S., 1940: Motion and Motion Systems. Vol. I, Weather Analysis and Forecasting, McGraw-Hill, 505 pp.

  • Press, H. P., S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery, 1992: Numerical Recipes in Fortran. 2d. ed. Cambridge University Press, 963 pp.

  • Pudykiewicz, J., and A. Staniforth, 1984: Some properties and comparative performance of the semi-Lagrangian method of Robert in the solution of the advection-diffusion equation. Atmos.–Ocean,22, 283–308.

  • ——, R. Benoit, and A. Staniforth, 1985: Preliminary results from a partial LRTAP model based on an existing meteorological forecast model. Atmos.–Ocean,23, 257–303.

  • Purnell, D. K., 1976: Solution of the advective equation by upstream interpolation with cubic splines. Mon. Wea. Rev.,104, 42–48.

  • Purser, R. J., and L. M. Leslie, 1991: An efficient interpolation procedure for high-order three-dimensional semi-Lagrangian models. Mon. Wea. Rev.,119, 2492–2498.

  • Saylor, R. D., and G. D. Ford, 1995: On the comparison of numerical methods for the integration of kinetic equations in atmospheric chemistry and transport models. Atmos. Environ.,29, 2585–2593.

  • Seibert, P., 1993: Convergence and accuracy of numerical methods for trajectory calculations. J. Appl. Meteor.,32, 558–566.

  • Smolarkiewicz, P. K., and J. A. Pudykiewicz, 1992: A class of semi-Lagrangian approximation for fluids. J. Atmos. Sci.,49, 2082–2096.

  • Staniforth, A., and J. Côté, 1991: Semi-Lagrangian integration schemes for atmospheric models—A review. Mon. Wea. Rev.,119, 2206–2223.

  • Sun, W.-Y., and K.-S. Yeh, 1997: A general semi-Lagrangian advection scheme employing forward trajectories. Quart. J. Roy. Meteor. Soc.,123, 2463–2476.

  • Zalesch, S. T., 1979: Fully multidimensional flux-corrected transport. J. Comput. Phys.,31, 335–352.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 878 165 24
PDF Downloads 549 99 14