《数值预报复习要点》由会员分享，可在线阅读，更多相关《数值预报复习要点（16页珍藏版）》请在金锄头文库上搜索。
1、第一章大气运动方程大气运动方程应遵循牛顿第二定律运动方程、质量守恒定律连续方程、能量守恒定律热力学或能量方程、气体试验定律状态方程、水汽守恒定律水汽方程运动方程连续方程状态方程热力学方程水汽方程球坐标下的大气运动方程组球坐标 经度、纬度，地心到空间点的距离r球坐标下 x是沿纬圈的微小位移、y是沿经圈的微小位移、z是垂直方向的微小位移x=r*cos*y=r*z=ru是纬向速度u=r* cos*d/dtv是经向速度v= r*d/dtw是垂直速度w=dr/dt球坐标下加速度的展开局地直角坐标系下坐标及球坐标的关系局地直角坐标系是球坐标系的简化形式，保存了球坐标的框架，忽略了球面曲率的影响。区分球坐标
2、方程与局地直角指标方程P坐标系下大气运动方程组P坐标系通常应用于天气尺度的大气运动，具有准静力平衡的特点，满足方程P坐标系下的垂直速度=dp/dt区分P坐标方程及Z坐标方程坐标坐标系是及气压相联系的坐标系，具有下边界简单，便于引进地形的动力作用等特点。坐标的定义pT是模式层顶的气压，ps是地面气压坐标的边界条件区分p坐标方程与坐标方程 状态方程 垂直运动方程数值模式的分类过滤模式只能模拟准地转演变过程，而原始方程模式既能模拟准地转演变过程又能模拟地转适应过程。原始方程模式分为正压原始方程模式垂直方向一层与斜压原始方程模式垂直方向有多层地图投影的概念地图投影是按照一定的数学条件，把球形的地球外表
3、展绘于平面地图上。或者说把地球外表投影到一个简单的曲面上。能够识别出光源、目标物、投影面投影的误差距离误差、面积误差、角度形状误差投影类型等角投影、等面积投影、任意投影等角投影中，经过投影后地球外表的任意两条交线的夹角保持不变，且在投影面任意一点的各个方向上长度放大或缩小的倍数相等，投影之后不产生角度或者形状的误差。按地图投影面的性质分类型方位投影或平面投影、圆锥投影、圆柱投影地图投影的根本概念及几个重要因子映像面：投影的投射面、投影面映像平面：映像面沿某一条经线切开所展成的平面地图：映像平面按地图比例尺缩小后的图切投影：映像面及地球外表相切于某一点的投影割投影：映像面及地球外表相割的投影标准
4、纬度：映像面及地球外表相交的纬度标准纬度上，映像面的距离等于地球外表上相应的距离映像比例尺m：映像平面上的距离除以地球外表上相应的距离。又称地图放大因子。标准纬度上m=1缩小比例尺：地图上任意纬度上的距离除以映像平面上相应的距离实际比例尺：地图上任意纬度上的距离除以地球外表上相应的距离正形投影正形投影的光源位于球心，映像面为圆锥面，映像面圆锥角为，标准纬度为0地图放大系数的计算其中k为单位经度所张的圆锥角，表示了圆锥的几何特征，成为圆锥常数；0是标准纬度的余角。极射赤面投影是一种正形割投影，其光源位于南极，映像面为一个及地球相割于北纬60度的平面，标准纬度为60N。根据网格坐标计算放大系数柯氏
5、参数的计算兰伯托投影是一种正形投影，其光源位于地球球心，映像面为一个及地球外表相割及30N与60N的圆锥面，圆锥角为90。麦卡托投影光源位于球心，映像面是及地球外表相割于南北纬22。5的圆柱面，标准纬度为22.5N与22.5S投影后，经线为等距平行的直线，纬线为及经线垂直的直线，正形圆锥投影的极限情形。k=0所以不能采用普遍的正形投影中的关系式来对之进展讨论。而是从地图放大系数入手求有关表达式Je为网格点相对于赤道的坐标。放大系数是关于赤道成纬向轴对称。普遍正交曲线坐标系中的方程组qj是正交曲线的坐标，dqj是相应的坐标变元，dlj是空间点沿坐标线所移动的距离，称为坐标线元dlj=Hj*dqj
6、其中Hj成为拉密系数。坐标线元不等于坐标变元而是等于坐标变元及拉密系数的乘积正交曲线坐标下的常用关系式气压梯度力涡度散度风速矢量平流绝对温度平流普遍地图投影坐标系中的方程组设X与Y轴地图投影放大系数为m与n，Z方向的地图投影放大系数为1拉密系数要求可以利用给出的关系式得到普遍地图投影坐标系中的大气方程组表达式例如根据连续方程表达式与散度在正交曲线坐标系下的表达式，得到地图投影坐标下的连续方程表达式，要求将求与符号展开成各项相加的形式第一步将H的表达式代入散度表达式，写出第二步写出第三步写出差分方法与差分格式离散化的概念u(x,t)是连续函数，u(ix,nt)是u(x,t)经离散化后的形式。所谓
7、离散化，即把连续的x以ix代替，连续的t以nt代替，其中i与n为整数。均是以一维线性平流方程为例 差分格式根本都是通过泰勒展开式来构造的。前差格式：后差格式：中央差格式：二阶微分的差分格式拉普拉斯的差分格式拉普拉斯的差分格式涉及到的格点截断误差上面差分格式中的R被成为截断误差。意思是用差商来近似代替偏微商时，将会因舍去R所代表的项而造成的误差。一般用R中最大的项来表示截断误差的大小。如果R中最大的项是x，那么R=Ox。如果是x2，那么R=Ox2。注：所谓最大的项指的是偏导阶数最小的一项，一阶偏导项大于二阶偏导项R反映了差分方程代替微分方程时的截断误差，它在一定程度上代表了差分格式的精度，R的阶
8、次越高，那么差分格式的精度越高，误差越小。这个的精度就是R=Ox2相容性一致性当空间步长 Dx 与时间步长Dt 很小时，差分方程是否逼近微分方程，这就是差分格式的相容性一致性问题。收敛性在一定的定解条件下，差分方程的解是否逼近微分方程的解的问题，称之为差分格式的收敛性问题。稳定性在时间积分过程中，由于舍入误差的影响，差分解的误差是否随时间增长的问题，即差分格式的计算稳定性问题。拉克斯Lax等价定理：如果差分方程逼近微分方程，即差分格式及微分方程是相容的，或者差分格式满足相容性条件，差分格式的稳定性，保证了其收敛性计算稳定性是收敛性的充分必要条件。 用Von-Neumann稳定性判别方法来证明差
9、分格式的计算稳定性时的主要步骤为：1设解的波动形式，代入差分方程。2得出其对应的增幅因子G。3讨论时的情况。4判断格式稳定性及满足格式稳定性的条件。CFL判据增幅因子G：其中A n+1与An分别是n+1时刻与n时刻的振幅。微分方程波动形式解：差分方程形式解：各个差分格式稳定性 时间前差，空间后差条件稳定时间前差，空间前差条件稳定时间前差，空间中央差任何情况下不成立，所以为绝对不稳定时间积分格式分类：1、二时间层的积分格式非迭代格式 a、欧拉格式：，绝对不稳定b、后差格式隐式格式，绝对稳定格式c、梯形格式隐式格式，中性格式，振幅不变2、二时间层的积分格式迭代格式特点：先通过前差格式算一个n+1时
10、刻的粗略预报值，然后在通过这个预报值用后差格式计算n+1时刻的准确预报值a、欧拉后差格式显式格式：条件稳定格式b、赫恩格式：绝对不稳定格式3、三个时间层的积分格式中央差格式，又称跳背格式条件稳定这种格式，可以看出不仅需要一个具有物理意义的初值u0，同样还需要一个出于计算要求的初值u1，前者称为物理初值，后者称为计算初值。差分格式误差相关内容：实际工作中，不可能任意缩小步长由于计算量过大等原因，实际计算是在有限的网格下进展的，一定程度的误差是不可防止的。中央差的计算解问题：使用中央差会产生两个波解，其中一个有物理意义，另一个不具有物理意义，是计算过程中产生的虚假波形，称之为计算解。中央差分格式的
11、数值解为两个波动的叠加。实际工作中通常采用提高网格分辩能力来抑制与减小计算解带来的影响。时间的截断误差频率误差当使用中央差格式显式时，原频率=f由替代当使用梯形格式时隐式，原频率=f由替代当t很小时，两种差分格式及=f差异很小当t增加时，差分格式解偏离原频率幅度增大ft1时显式格式会出现不稳定，隐式格式的频率解随ft增加而减小显式格式及隐式格式的频率解及真值的比值空间的截断误差波数误差用中央差格式展开时，波数k的数值解为：差分近似精度随k或x的减小而增大，这也就是说对于波长较短的波，其产生的波数误差较大；而波长较长的波，那么差分方程可以比拟准确地表示其空间微商，波数误差很小，精度较高。相速度与
12、群速度误差空间差分格式的波数误差与时间积分格式的频率误差会造成相速度与群速度的误差，从而引起计算频散。差分格式在波的移动与能量传播方面均可造成误差。而且： 1由于相速度误差，减慢了平流过程；2造成虚假的计算频散，且对短波尤为明显。1波长越长，误差越小；波长减小，其误差也就更为严重；2提高网格分辩率，使x取得足够小，可以提高相速度的准确率。差分格式误差特征总结：1、三层时间积分格式存在计算解问题： 计算解对差分解的影响依赖于网格分辨率与波长。2、时间积分格式引起频率误差：显式格式使其频率明显增加，振动加快；隐式格式使其频率明显减小，振动减慢。3、空间差分格式引起波数误差： 高阶差分格式所引起的波
13、数误差要比低阶格式小； 波长较短的波，误差尤为严重。4、空间差分格式会引起计算频散： 尤其对于短波，相速度与群速度均会产生很大的误差。 通常可采用提高网格分辩率的方法减小各种误差。非线性不稳定对于非线性偏微分方程，线性偏微分方程的稳定性条件，只能给出其计算稳定性的必要条件，即使满足这一条件，也可能会因为差分方程的边界条件与非线性项的不正确表示而产生计算的不稳定现象，我们把这种由于非线性作用而产生的不稳定，称为非线性不稳定。混淆误差差分方法是用有限的自由度系统来代替原来的连续介质系统的。而有限格点上的函数值只能分解有限的波数，其最短波长为2x，对于非线性作用产生的波长小于2x的波动，网格系统不能
14、正确地将它表示出来，而把它错误地表示成为某种波长大于2x的波，从而产生了误差，我们把这种波的误差称之为混淆误差。对于网格数为I的网格，假设ui中包含两个波，波数分别为k1与k2，当这两个波发生非线性作用而产生k1+k2的波时，如果k1+k2I/2，那么网格将把这个波的波数识别为S=I-k1+k2自激反应假设k1与k2相互作用，能量反应到上k1，那么有：k1=I-k1+k2，或2k1=I-k2。由于k2=I/4。因此，发生能量反应的波其波长必在2x与4x之间，也就是说只有对波长很短的波才发生能量反应；非线性不稳定的产生主要是由于短波能量的虚假增长所造成的。抑制与克制非线性计算不稳定的做法通常有： 空间与时间平滑，滤去短波分量滤除波长小于4倍格距的波动； 在方程中参加扩散项； 构造具有隐式平滑与某种选择性衰减作用的差分格式 构造守衡的差分格式，使差分方程尽可能保持原来的 物理规律与能量关系； 采用谱变换方法 可防止非线性不稳定问题。正压原始方程模式以下在3个假设的根底上导出正压原始方程模式的预报方程组第一、假设大气均匀不可压，密度为一常数。第二、假设模式满足静力平衡。第三 、假设模式大气是正压的，初始时刻水平风速不随气压变化正压原始方程模式的预报方
