• Arakawa, A., 1972: Design of the UCLA general circulation model: Numerical simulation of weather and climate. Department of Atmospheric Sciences Tech. Rep. 7, University of California, Los Angeles, 116 pp.

  • Baba, Y., K. Takahashi, and T. Sugimura, 2010: Dynamical core of an atmospheric general circulation model on a Yin–Yang grid. Mon. Wea. Rev., 138, 39884005, https://doi.org/10.1175/2010MWR3375.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Beylkin, G., J. M. Keiser, and L. Vozovoi, 1998: A new class of time discretization schemes for the solution of nonlinear PDEs. J. Comput. Phys., 147, 362387, https://doi.org/10.1006/jcph.1998.6093.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bhatt, A., and B. E. Moore, 2017: Structure-preserving exponential Runge–Kutta methods. SIAM J. Sci. Comput., 39, A593A612, https://doi.org/10.1137/16M1071171.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bhatt, A., and B. E. Moore, 2019: Exponential integrators preserving local conservation laws of PDEs with time-dependent damping/driving foresee. J. Comput. Appl. Math., 352, 341351, https://doi.org/10.1016/j.cam.2018.12.003.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Botchev, M. A., 2013: A block Krylov subspace time-exact solution method for linear ordinary differential equation systems. Numer. Linear Algebra Appl., 20, 557574, https://doi.org/10.1002/nla.1865.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Bott, A., 1989: A positive definite advection scheme obtained by nonlinear renormalization of the advective fluxes. Mon. Wea. Rev., 117, 10061016, https://doi.org/10.1175/1520-0493(1989)117<1006:APDASO>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Clancy, C., and J. A. Pudykiewicz, 2013: On the use of exponential time integration methods in atmospheric models. Tellus, 65A, 20898, https://doi.org/10.3402/tellusa.v65i0.20898.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Durran, D. R., 1991: The third-order Adams–Bashforth method: An attractive alternative to leapfrog time difference. Mon. Wea. Rev., 119, 702720, https://doi.org/10.1175/1520-0493(1991)119<0702:TTOABM>2.0.CO;2.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Frommer, A., K. Lund, and D. B. Szyld, 2017: Block Krylov subspace methods for functions of matrices. Electron. Trans. Numer. Anal., 47, 100126.

    • Search Google Scholar
    • Export Citation
  • Gaudreault, S., G. Rainwater, and M. Tokman, 2018: KIOPS: A fast adaptive Krylov subspace solver for exponential integrators. J. Comput. Phys., 372, 236255, https://doi.org/10.1016/j.jcp.2018.06.026.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gong, J., B. Wang, and Z. Ji, 2013a: Three-step difference scheme for solving nonlinear time-evolution partial differential equations. Atmos. Oceanic Sci. Lett., 6, 423427, https://doi.org/10.1080/16742834.2013.11447119.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Gong, J., B. Wang, and Z. Ji, 2013b: Generalized square conservative multistep finite difference scheme incorporating historical observations. Atmos. Oceanic Sci. Lett., 6, 223226, https://doi.org/10.1080/16742834.2013.11447084.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Hochbruck, M., C. Lubich, and H. Selhofer, 1998: Exponential integrators for large systems of differential equations. SIAM J. Sci. Comput., 19, 15521574, https://doi.org/10.1137/S1064827595295337.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Ji, Z., and Q. C. Zeng, 1982: The construction and application of difference scheme of evolution equations. Sci. Atmos. Sin., 6, 8894.

    • Search Google Scholar
    • Export Citation
  • Ji, Z., and B. Wang, 1991: A further discussion on the construction and application of difference scheme of evolution equations (in Chinese). J. Atmos. Sci., 15, 7281.

    • Search Google Scholar
    • Export Citation
  • Koskela, A., 2015: Approximating the matrix exponential of an advection–diffusion operator using the incomplete orthogonalization method. Numerical Mathematics and Advanced Applications—ENUMATH 2013, A. Abdulle et al., Eds., Lecture Notes in Computational Science and Engineering, Vol. 103, Springer, 345–353.

    • Crossref
    • Export Citation
  • Li, D., Y. Cong, and K. Xia, 2016: Flexible exponential integration methods for large systems of differential equations. J. Appl. Math. Comput., 51, 545567, https://doi.org/10.1007/s12190-015-0919-1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Loffeld, J., and M. Tokman, 2013: Comparative performance of exponential, implicit and explicit integrators for stiff systems of ODEs. J. Comput. Appl. Math., 241, 4567, https://doi.org/10.1016/j.cam.2012.09.038.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Moler, C., and C. Van Loan, 1978: Nineteen dubious ways to compute the exponential of a matrix. SIAM Rev., 20, 801836, https://doi.org/10.1137/1020098.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Moler, C., and C. Van Loan, 2003: Nineteen dubious ways to compute the exponential of a matrix, twenty-five years later. SIAM Rev., 45, 349, https://doi.org/10.1137/S00361445024180.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Nair, R. D., and C. Jablonowski, 2008: Moving vortices on the sphere: A test case for horizontal advection problems. Mon. Wea. Rev., 136, 699711, https://doi.org/10.1175/2007MWR2105.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Niesen, J., and W. M. Wright, 2012: Algorithm 919: A Krylov subspace algorithm for evaluating the φ-functions appearing in exponential integrators. ACM Trans. Math. Software, 38, 119, https://doi.org/10.1145/2168773.2168781.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Skamarock, W. C., and M. L. Weisman, 2009: The impact of positive-definite moisture transport on NWP precipitation forecasts. Mon. Wea. Rev., 137, 488494, https://doi.org/10.1175/2008MWR2583.1.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Tokman, M., 2006: Efficient integration of large stiff systems of ODEs with Exponential Propagation Iterative (EPI) methods. J. Comput. Phys., 213, 748776, https://doi.org/10.1016/j.jcp.2005.08.032.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Vo, H. D., and R. B. Sidje, 2017: Approximating the large sparse matrix exponential using incomplete orthogonalization and Krylov subspaces of variable dimension. Numer. Linear Algebra Appl., 24, e2090, https://doi.org/10.1002/nla.2090.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Wang, B., and Z. Z. Ji, 1993: The harmonious dissipative operator and the completely square conservative difference scheme in explicit way. Sci. China, 37B, 462471, https://doi.org/10.1360/yb1994-37-4-462.

    • Search Google Scholar
    • Export Citation
  • Wang, B., Z. Z. Ji, and Q. C. Zeng, 1993: A kind of time-saving explicit numerical computational method. Chin. Sci. Bull., 37, 766768.

    • Search Google Scholar
    • Export Citation
  • Wang, B., Q. Zeng, and Z. Z. Ji, 1995: Square conservation systems and Hamiltonian systems. Sci. China, 38A, 12111219.

  • Wang, B., H. Wan, Z. Ji, X. Zhang, R. Yu, Y. Yu, and H. Liu, 2004: Design of a new dynamical core for global atmospheric models based on some efficient numerical methods. Sci. China, 47A, 421, https://doi.org/10.1360/04za0001.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Williamson, D. L., J. B. Drake, J. J. Hack, R. Jakob, and P. N. Swarztrauber, 1992: A standard test set for numerical approximation to the shallow water equations in spherical geometry. J. Comput. Phys., 102, 211224, https://doi.org/10.1016/S0021-9991(05)80016-6.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Xiao, F., and X. D. Peng, 2004: A convexity preserving scheme for conservative advection transport. J. Comput. Phys., 198, 389402.

  • Zeng, Q. C., and Z. Z. Ji, 1981: On the computational stability of evolution equations. Chin. Comput. Math., 3, 7986.

  • Zerroukat, M., N. Wood, and A. Staniforth, 2004: SLICE-S: A Semi-Lagrangian Inherently conserving and Efficient scheme for transport problems on sphere. Quart. J. Roy. Meteor. Soc., 130, 26492664, https://doi.org/10.1256/qj.03.200.

    • Crossref
    • Search Google Scholar
    • Export Citation
  • Zhang, Y., R. Yu, J. Li, and H. M. Chen, 2013: An implementation of a leaping-point two-step shape-preserving advection scheme in high-resolution spherical latitude-longitude grid (in Chinese). Acta Meteor. Sin., 71, 10891102, https://doi.org/10.11676/qxxb2013.085https://doi.org/10.1007/s10114-013-2274-1.

    • Crossref
    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 42 42 35
Full Text Views 11 11 7
PDF Downloads 17 17 11

A New Algorithm with Square Conservative Exponential Integral for Transport Equation

View More View Less
  • 1 Zhejiang Meteorological Observatory, Hangzhou, Zhejiang, China
  • 2 Heavy Rain and Drought-Flood Disasters in Plateau and Basin Key Laboratory of Sichuan Province, Chengdu, SiChuan, China
  • 3 Key Laboratory of Cloud-Precipitation Physics and Severe Storms (LACS), Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing, China
  • 4 Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science and Technology, Nanjing, China
  • 5 University of Chinese Academy of Sciences, Beijing, China
  • 6 NOAA/NCEP/Environmental Modeling Center, I. M. System Group, College Park, Maryland
© Get Permissions
Restricted access

Abstract

The square conservative exponential integral method (SCEIM) is proposed for transport problems on the sphere. The method is a combination of the square conservation algorithm and the exponential integral method. The main emphasis in the development of SCEIM is on conservation, positive-definite, and reversibility as well as achieving comparable accuracy to other published schemes. The most significant advantage of SCEIM is to change the forward model to the backward model by setting a negative time step, and the backward model can be used to solve the inverse problem. Moreover, the polar problem is significantly improved by using a simple effective central skip-point difference scheme without major penalty on the overall effectiveness of SCEIM. To demonstrate the effectiveness and generality of the SCEIM, this method is evaluated by standard cosine bell tests and deformational flow tests. The numerical results show that SCEIM is a time-convergence method as well as a grid-convergence method, and has a strong shape-preserving ability. In the tests of the inverse problem, the sharp fronts are successfully regressed back into their initial weak fronts and the cosine bells move against the wind direction and return to the initial position with high accuracy. The numerical results of forward simulations are compared with those of published schemes, the total mass conservation, and error norms are competitive in term of accuracy.

© 2021 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Dr. ShifengHao, shifenghao@aliyun.com

Abstract

The square conservative exponential integral method (SCEIM) is proposed for transport problems on the sphere. The method is a combination of the square conservation algorithm and the exponential integral method. The main emphasis in the development of SCEIM is on conservation, positive-definite, and reversibility as well as achieving comparable accuracy to other published schemes. The most significant advantage of SCEIM is to change the forward model to the backward model by setting a negative time step, and the backward model can be used to solve the inverse problem. Moreover, the polar problem is significantly improved by using a simple effective central skip-point difference scheme without major penalty on the overall effectiveness of SCEIM. To demonstrate the effectiveness and generality of the SCEIM, this method is evaluated by standard cosine bell tests and deformational flow tests. The numerical results show that SCEIM is a time-convergence method as well as a grid-convergence method, and has a strong shape-preserving ability. In the tests of the inverse problem, the sharp fronts are successfully regressed back into their initial weak fronts and the cosine bells move against the wind direction and return to the initial position with high accuracy. The numerical results of forward simulations are compared with those of published schemes, the total mass conservation, and error norms are competitive in term of accuracy.

© 2021 American Meteorological Society. For information regarding reuse of this content and general copyright information, consult the AMS Copyright Policy (www.ametsoc.org/PUBSReuseLicenses).

Corresponding author: Dr. ShifengHao, shifenghao@aliyun.com
Save