基于三次插值函数算法的时间积分方案与二阶时空 - 中国气象局广州 .pdf

27 卷 5 期 热 带 气 象 学 报 Vol.27 ，No.5 20 11 年 10 月 JOURNAL OF TROPICAL METEOROLOGY Oct. ，20 11 辜旭赞，张兵，王明欢. 基于三次插值函数算法的时间积分方案与二阶时空余差数值模式 以原始大气运动方程与理想全球模拟个例为例[J]. 热 带气象学报，2011 ，27(5): 669-678. 文章编号 1004-4965(2011)05-0669-10 基于三次插值函数算法的时间积分方案与二阶时空 余差数值模式 以原始大气运动方程与理想全球模拟个例为例 辜旭赞，张兵，王明欢 (中国气象局武汉暴雨研究所，湖北 武汉 430074) 摘 要: 从大气运动原始方程和欧拉算符出发，用泰勒级数展开，给出二阶时空微商余项预报方程。进 而讨论用三次插值函数—— 三次曲面拟合求上游点的准拉格朗日时间积分方案与相应的二阶时空余差数 值模式—— “ 三次模式”。则 三次模式是通过实现各个大气物理量场的二阶可导，从而可对预报方程 做空间非线性(“三次”)时间离散积分，成为“ 三次曲面拟合——时间步积分—— 三次曲面拟合——……” 一种新算法数值模式。讨论 三次数值模式的数学基础：三次插值函数及其数值分析极性定律用于数值模 式。指出： 三次模式和谱模式都具有数学 “收敛性”；而 Coons 三次曲面具有对变量场拟合二阶可导 “最优性”；和Hermite 三次曲面片具有对网格变量场二阶可导运算 “等价性”。又指出：有限差分模式 的中央差近似斜率和曲率，分别是三次样条斜率和曲率作“三点平滑”。 三次模式适合采用原始大气运 动方程，适合采用准拉格朗日时间积分方案，并给出一个理想全球模拟个例。因大气运动本质上是非线性 的，理论上可按变量场 三次曲面曲率判断，以采用符合物理诠释的局域或单点平滑，以保持模式时间积 分稳定性。且未来容易实现全球多重/时变套网格 三次数值模式。 关 键 词 数值预报与模拟；二阶时空余差数值模式；三次插值函数；Navier-Stokes 大气运动方程 组；准拉格朗日时间积分；欧拉向前差分时间积分；全球模拟个例 中图分类号 P435 文献标识码 A Doi ：10.3969/j.issn.1004-4965.2011.05.007 完整描述与实现。 1 言 本文重点讨论采用 三次曲面插值的准拉格 朗日时间积分方案[7 －14] 、时空微商余项分别为 数值模式的动力框架核心是时间积分方案， 2 2 2 O( Δt)和O( Δx , Δy , Δz) 的新算法数值模式， 它决定了模式大气动力系统的物理机制、数学精 可称为“ 三次数值模式”，简称“ 三次模式”， 度和计算误差，同时也决定模式描述大气运动的 或“三次模式”。它在计算精度上高于所谓“双线性 物理保守性(如“质量守恒”、“ 能量守恒”等) 、数 插值”[13－14]准拉格朗