双月刊

ISSN 1006-9585

CN 11-3693/P

+高级检索 English
高阶Runge-Kutta-Li算法对二维线性平流方程的计算检验
作者:
作者单位:

1.中国科学院大气物理研究所大气科学和地球流体力学数值模拟国家重点实验室,北京 100029;2.中国科学院大气物理研究所季风系统研究中心,北京 100190;3.中国海洋大学物理海洋教育部重点实验室,山东青岛266003;4.青岛海洋科学与技术国家实验室区域海洋动力学与数值模拟功能实验室,山东青岛266237

作者简介:

王鹏飞,男,1973年出生,博士,高级工程师,主要从事数值模式、并行计算、非线性可预报性等方面的研究。E-mail: wpf@mail.iap.ac.cn

通讯作者:

基金项目:

国家重点研发计划2018YFA0605904,国家自然科学基金资助项目41530426、41831175、41425019,海洋局国际合作项目GASI-IPOVAI-03


High Order Runge-Kutta-Li Scheme to Solve Two-dimensional Linear Advection Equations
Author:
Affiliation:

1.State Key Laboratory of Numerical Modeling for Atmospheric Sciences and Geophysical Fluid Dynamics, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029;2.Center for Monsoon System Research, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100190;3.Key Laboratory of Physical Oceanography Institute for Advanced Ocean Studies, Ocean University of China, Qingdao, Shandong Province 266003;4.Laboratory for Regional Oceanography and Numerical Modeling, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong Province 266237

Fund Project:

National Key Research and Development Program of China Grant 2018YFA0605904;National Natural Science Foundation of China Grants 41530426 41831175 41425019;SOA International Cooperation Program on Global Change and Air-Sea Interactions Grant GASI-IPOVAI-03National Key Research and Development Program of China (Grant 2018YFA0605904), National Natural Science Foundation of China (Grants 41530426, 41831175, and 41425019), SOA International Cooperation Program on Global Change and Air-Sea Interactions (Grant GASI-IPOVAI-03)

  • 摘要
  • |
  • 图/表
  • |
  • 访问统计
  • |
  • 参考文献
  • |
  • 相似文献
  • |
  • 引证文献
  • |
  • 资源附件
    摘要:

    利用高阶Li空间微分方案(Li, 2005),实现了时间积分为3~6阶Runge-Kutta-Li(RKL)格式的求解算法。二维线性平流方程的试验结果表明:在计算稳定的条件下,各阶算法的计算误差随时间的推移基本上是线性增加的。非转动背景场的平流算例中(高斯型的初值),高阶RKL算法可以取得较好的计算效果。与3、4、5、6阶RK算法配合的Li空间差分方案有效阶数可以达到5、7、9、10阶。RK 算法的阶数为5(6)阶时,总误差控制在10-7(10-8)以内。随RK阶数增加Li微分的有效阶数有增加趋势,且总误差逐渐减小。定常转速的背景场算例中(偏心的高斯型初值),当RK阶数为3时,最优空间差分阶数为10;相应的阶数为4、5、6时对应的空间最优阶为16,22,22,总计算误差可以控制在10-15~10-16。随着精度的提高,误差的绝对值减小很迅速,说明算法是非常有效的。对于圆锥型初值(定常转速的背景场),4、5、6阶RK算法和3阶算法的效果差不多。高阶算法对此类具有导数不连续点的算例,效果不如高斯初始场好,结果不能保持正定,有些地方误差出现下冲和上翘。随着空间差分精度的提高,非正定的解数量和数值减小,误差的绝对值减小,说明了算法在一定程度上是有效的,但并不适合追求极高的算法阶数。这与谱方法中的导数不连续问题有些相似,误差的产生主要源于导数的不连续性,差分类方法仅能获得与导数连续性阶数相当的算法精度。各种算例中,采用恰当的边界条件是必要的,例如旋转背景场算例,比较适合使用无穷远边界条件,否则会出现计算不稳定或无法将计算误差控制到较小的范围内。

    Abstract:

    In this study, aiming to take full advantages of Li’s high-order spatial differential method (Li, 2005), we implement the hybrid Runge-Kutta-Li (RKL) scheme to solve a two-dimensional (2D) linear advection equation. The results indicate that the computation error increased linearly with time. The experiments with a no-rotate background field of the Gaussian initial values by RKL scheme could obtain a precise result. The effective spatial orders were 5, 7, 9, and 10, corresponding to temporal orders of 3, 4, 5, and 6, respectively. The fifth- (sixth-) order Runge-Kutta (RK) integration scheme with the ninth- (tenth-) order Li’s difference scheme in spatial direction controlled the error within 10-7(10-8). The effective order of Li’s scheme (Li, 2005) tended to increase with the increase in the RK order, and the total error gradually decreased.Another rotated background field integrated from an eccentric Gaussian-type initial featured similar results. The effective spatial order was 10 when a third-order RK scheme was applied, and they increased to 16, 22, and 22 when the order of RK scheme changed to 4, 5, and 6, respectively. The computation error could be controlled within 10-15-10-16, and the peak of Gaussian initial values were well maintained. The error decreased sharply while the order of RKL scheme increased, which indicates that the RKL is very effective to deal with such problem. The experiments of RKL scheme to solve a cone initial case (with rotated background filed) indicate that the fourth, fifth, and sixth RK integration obtained almost the same precision result as the third-order RK scheme. The high-order scheme was not as effective as it was for the Gaussian initial condition when it addressed a problem that had discontinuous derivates. The computed solution was not positive in the whole grids, and in some places, the error was downward, while it was upward in some other places. The increase of spatial order could make the error smaller, but the error descent was not very sharp. This result suggests that the high-order spatial difference scheme has some benefits, but we should not expect that an ultra-high-order scheme will lead to an ultra-high precise result. This phenomenon reveals that the high-order scheme is limited by the continuous property of the initial condition, and as a result, the error order is directly proportional to the order of derivates of the initial condition. The proper boundary conditions are important for the above computation cases when the RKL scheme is applied. For instance, in the computation of the rotated background cases, the value outside the grid box tends to 0 at , which is a feasible boundary condition. An improper boundary condition may cause the computation to be unstable, or the error cannot be controlled to an acceptable range.

    参考文献
    相似文献
    引证文献
引用本文

王鹏飞,李建平,黄刚.2019.高阶Runge-Kutta-Li算法对二维线性平流方程的计算检验[J].气候与环境研究,24(4):417-429. WANG Pengfei,LI Jianping,and HUANG Gang.2019.High Order Runge-Kutta-Li Scheme to Solve Two-dimensional Linear Advection Equations[J].Climatic and Environmental Research (in Chinese],24(4):417-429.

复制
分享
文章指标
  • 点击次数:
  • 下载次数:
  • HTML阅读次数:
  • 引用次数:
历史
  • 收稿日期:2018-12-26
  • 最后修改日期:
  • 录用日期:
  • 在线发布日期: 2019-08-08
  • 出版日期: