一种“准拉格朗日法”和“欧拉法”统一算法时间积分方案

从欧拉算符出发,用泰勒级数展开,给出二阶时空微商余项预报方程,进而讨论一种数值分析新算法———双三次曲面拟合(插值)的准拉格朗日时间积分方案与数值模式。它是将大气运动描述成为非线性的三次运动,即是通过对原始大气运动方程中包括标量、矢量的压、温、湿、风、以及旋转地球上广义牛顿力加速度场和散度场等,做双三次曲面拟合,实现对各个大气运动变量场的二阶可导,即限定气块上游点在各个不同双三次曲面(片)、具有斜率、曲率和挠率的非线性三次变量场上活动,从而可对各个大气运动方程做时间离散积分,即为双三次曲面拟合—时间步积分—双三次曲面拟合……,实现成为一种新动力框架数值模式。由于双三次曲面具有数学定律收敛性和二阶可导最优性,故选用双三次曲面插值求算二阶余差上游点,具有充分必要的数学理由:它包含了大气运动变量场之斜率、曲率和挠率。因此,埃尔米特双三次曲面片具有对网格变量场二阶可导运算等价性、及其数学收敛性与最佳曲率最优性,并且将准拉格朗日法与欧拉法,以及柯朗-弗里德里希斯-列维判据统一起来。容易实现全球网格变量场的双三次曲面拟合,和可按双三次曲面变量场的斜率、曲率或挠率判断,作变量场局域或单点平滑,以此保持三次模式的时间积分稳定...

A new quasi-Lagrangian time integration scheme with the interpolation of fitting bicubic surface

Starting from the Eulerian operator,the forecasting equations with 2-order time and space differential remainder are derived in light of the Taylor series expansion.The quasi-Lagrangian integration scheme with an algorithm of numerical analysis the-bicubic surface interpolation-for a new meteorological numerical model is then suggested.It describes atmospheric motions just as non-linear cubic waves,with fitting different bicubic surfaces to physical scalar or vector variable fields in the model atmosphere,i...