Mass Correction Applied to Semi-Lagrangian Advection Scheme

Wen-Yih Sun Department of Earth and Atmospheric Sciences, Purdue University, West Lafayette, Indiana

Search for other papers by Wen-Yih Sun in
Current site
Google Scholar
PubMed
Close
and
Ming-Teh Sun University Information Technology Services, Indiana University at Bloomington, Bloomington, Indiana

Search for other papers by Ming-Teh Sun in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

A simple mass correction is proposed for the semi-Lagrangian scheme using forward trajectories. The procedure includes (a) constructing the Lagrangian network induced by the motion of the fluid from the Eulerian network and finding the intersections of the networks by a general interpolation from the irregularly distributed Lagrangian grid to the regularly distributed Eulerian grid, (b) applying the spatial filter to remove the unwanted short waves and the values beyond the constraints, and finally (c) introducing a polynomial or sine function as the correction function that conserves mass while inducing the least modification to the results obtained from (a) and (b). Numerical simulations of pure advection, rotation, and idealized cyclogenesis show that the mass correction function does not degrade the accuracy of the results.

Corresponding author address: Dr. Wen-Yih Sun, Department of Earth and Atmospheric Sciences, Purdue University, 550 Stadium Mall Drive, West Lafayette, IN 47907-2051. Email: wysun@purdue.edu

Abstract

A simple mass correction is proposed for the semi-Lagrangian scheme using forward trajectories. The procedure includes (a) constructing the Lagrangian network induced by the motion of the fluid from the Eulerian network and finding the intersections of the networks by a general interpolation from the irregularly distributed Lagrangian grid to the regularly distributed Eulerian grid, (b) applying the spatial filter to remove the unwanted short waves and the values beyond the constraints, and finally (c) introducing a polynomial or sine function as the correction function that conserves mass while inducing the least modification to the results obtained from (a) and (b). Numerical simulations of pure advection, rotation, and idealized cyclogenesis show that the mass correction function does not degrade the accuracy of the results.

Corresponding author address: Dr. Wen-Yih Sun, Department of Earth and Atmospheric Sciences, Purdue University, 550 Stadium Mall Drive, West Lafayette, IN 47907-2051. Email: wysun@purdue.edu

Save
  • Bates, J. R., 1984: An efficient semi-Lagrangian and alternating- direction implicit method for integrating the shallow-water equations. Mon. Wea. Rev, 112 , 20332047.

    • Search Google Scholar
    • Export Citation
  • Bermejo, R., and A. Staniforth, 1992: The conversion of semi-Lagrangian advection schemes to quasi-monotone schemes. Mon. Wea. Rev, 120 , 26222632.

    • Search Google Scholar
    • Export Citation
  • Bermejo, R., and J. Conde, 2002: A conservative quasi-monotone semi- Lagrangian scheme. Mon. Wea. Rev, 130 , 423430.

  • Chern, J-D., 1994: Numerical simulation of cyclogenesis over the western United States. Ph.D. thesis, Purdue University, 178 pp.

  • Doswell III, C. A., 1984: A kinematic analysis associated with a nondivergent vortex. J. Atmos. Sci, 41 , 12421248.

  • Gravel, S., and A. Staniforth, 1994: A mass conservative semi-Lagrangian scheme for the shallow wave equations. Mon. Wea. Rev, 122 , 243248.

    • Search Google Scholar
    • Export Citation
  • Huang, C-Y., 1994: Semi-Lagrangian advection schemes and Eulerian WKL algorithms. Mon. Wea. Rev, 122 , 16471658.

  • Hólm, E. V., 1995: A fully two-dimensional, nonoscillatory advection scheme for momentum and scalar transport equations. Mon. Wea. Rev, 123 , 536552.

    • Search Google Scholar
    • Export Citation
  • Kuo, H-C., and R. T. Williams, 1990: Semi-Lagrangian solutions to the inviscid Burgers equation. Mon. Wea. Rev, 118 , 12781288.

  • Lin, S-J., and R. B. Rood, 1996: Multimensional flux-form semi- Lagrangian transport schemes. Mon. Wea. Rev, 124 , 20462069.

  • Priestley, A., 1993: A quasi-conservative version of the semi-Lagrangian advection scheme. Mon. Wea. Rev, 121 , 621629.

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

    • Search Google Scholar
    • Export Citation
  • Purser, R. J., and L. M. Leslie, 1994: An efficient semi-Lagrangian scheme using third-order semi-implicit time integration and forward trajectories. Mon. Wea. Rev, 122 , 745756.

    • Search Google Scholar
    • Export Citation
  • Rancic, M., 1992: Semi-Lagrangian piecewise biparabloic scheme for two-dimensional horizontal advection of a passive scalar. Mon. Wea. Rev, 120 , 13941406.

    • Search Google Scholar
    • Export Citation
  • Smolarkiewicz, P. K., and P. J. Rasch, 1991: Monotonic advection on the sphere: An Eulerian versus semi-Lagrangian approach. J. Atmos. Sci, 48 , 793810.

    • Search Google Scholar
    • Export Citation
  • Staniforth, A., and J. Côté, 1991: Semi-Lagrangian integration schemes for atmospheric models—A review. Mon. Wea. Rev, 119 , 22062223.

    • Search Google Scholar
    • Export Citation
  • Sun, W. Y., 2002: Numerical modeling in the atmosphere. J. Korean Meteor. Soc, 38 , 237251.

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

    • Search Google Scholar
    • Export Citation
  • Sun, W. Y., K. S. Yeh, and R-Y. Sun, 1996: A simple semi-Lagrangian scheme for advection equations. Quart. J. Roy. Meteor. Soc, 122 , 12111226.

    • Search Google Scholar
    • Export Citation
  • Wiin-Nielsen, A., 1959: On the application of trajectory methods in numerical forecasting. Tellus, 11 , 180196.

  • Williamson, D. L., and P. J. Rasch, 1989: Two-dimensional semi-Lagrangian transport with shape preserving interpolation. Mon. Wea. Rev, 117 , 102129.

    • Search Google Scholar
    • Export Citation
  • Zalesak, S. T., 1979: Fully multi-dimensional flux-corrected transport. J. Comput. Phys, 31 , 335362.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 86 28 2
PDF Downloads 50 21 3