数值预报复习要点.docx

数值预报复习要点.docx

《数值预报复习要点.docx》由会员分享，可在线阅读，更多相关《数值预报复习要点.docx（14页珍藏版）》请在冰豆网上搜索。

数值预报复习要点

第一章

大气运动方程

大气运动方程应遵循牛顿第二定律（运动方程）、质量守恒定律（连续方程）、能量守恒定律（热力学或能量方程）、气体试验定律（状态方程）、水汽守恒定律（水汽方程）

运动方程

连续方程

状态方程

热力学方程

水汽方程

球坐标下得大气运动方程组

球坐标经度λ、纬度φ,地心到空间点得距离r

球坐标下δx就是沿纬圈得微小位移、δy就是沿经圈得微小位移、δz就是垂直方向得微小位移

δx=r*cosφ*δλ

δy=r*δφ

δz=δr

u就是纬向速度u=r*cosφ*dλ/dt

v就是经向速度v=r*dφ/dt

w就是垂直速度w=dr/dt

球坐标下加速度得展开

局地直角坐标系下坐标与球坐标得关系

局地直角坐标系就是球坐标系得简化形式,保留了球坐标得框架,忽略了球面曲率得影响。

区分球坐标方程与局地直角指标方程

P坐标系下大气运动方程组

P坐标系通常应用于天气尺度得大气运动,具有准静力平衡得特点,满足方程

P坐标系下得垂直速度ω=dp/dt

区分P坐标方程与Z坐标方程

σ坐标

σ坐标系就是与气压相联系得坐标系,具有下边界简单,便于引进地形得动力作用等特点。

σ坐标得定义

pT就是模式层顶得气压,ps就是地面气压

σ坐标得边界条件

区分p坐标方程与σ坐标方程

状态方程

垂直运动方程

数值模式得分类

过滤模式只能模拟准地转演变过程,而原始方程模式既能模拟准地转演变过程又能模拟地转适应过程。

原始方程模式分为正压原始方程模式（垂直方向一层）与斜压原始方程模式（垂直方向有多层）

地图投影得概念

地图投影就是按照一定得数学条件,把球形得地球表面展绘于平面地图上。

或者说把地球表面投影到一个简单得曲面上。

能够识别出光源、目标物、投影面

投影得误差

距离误差、面积误差、角度（形状）误差

投影类型

等角投影、等面积投影、任意投影

等角投影中,经过投影后地球表面得任意两条交线得夹角保持不变,且在投影面任意一点得各个方向上长度放大或缩小得倍数相等,投影之后不产生角度或者形状得误差。

按地图投影面得性质分类型

方位投影或平面投影、圆锥投影、圆柱投影

地图投影得基本概念及几个重要因子

映像面:

投影得投射面、投影面

映像平面:

映像面沿某一条经线切开所展成得平面

地图:

映像平面按地图比例尺缩小后得图

切投影:

映像面与地球表面相切于某一点得投影

割投影:

映像面与地球表面相割得投影

标准纬度:

映像面与地球表面相交得纬度（标准纬度上,映像面得距离等于地球表面上相应得距离）

映像比例尺m:

映像平面上得距离除以地球表面上相应得距离。

又称地图放大因子。

标准纬度上m=1

缩小比例尺:

地图上任意纬度上得距离除以映像平面上相应得距离

实际比例尺:

地图上任意纬度上得距离除以地球表面上相应得距离

正形投影

正形投影得光源位于球心,映像面为圆锥面,映像面圆锥角为α,标准纬度为φ0

地图放大系数得计算

其中k为单位经度所张得圆锥角,表示了圆锥得几何特征,成为圆锥常数;θ0就是标准纬度得余角。

极射赤面投影就是一种正形割投影,其光源位于南极,映像面为一个与地球相割于北纬60度得平面,标准纬度为60°N。

根据网格坐标计算放大系数

柯氏参数得计算

兰伯托投影就是一种正形投影,其光源位于地球球心,映像面为一个与地球表面相割与30°N与60°N得圆锥面,圆锥角为90°。

麦卡托投影光源位于球心,映像面就是与地球表面相割于南北纬22。

5°得圆柱面,标准纬度为22、5°N与22、5°S

投影后,经线为等距平行得直线,纬线为与经线垂直得直线,正形圆锥投影得极限情形。

k=0所以不能采用普遍得正形投影中得关系式来对之进行讨论。

而就是从地图放大系数入手求有关表达式

Je为网格点相对于赤道得坐标。

放大系数就是关于赤道成纬向轴对称。

普遍正交曲线坐标系中得方程组

qj就是正交曲线得坐标,dqj就是相应得坐标变元,dlj就是空间点沿坐标线所移动得距离,称为坐标线元

dlj=Hj*dqj

其中Hj成为拉密系数。

坐标线元不等于坐标变元而就是等于坐标变元与拉密系数得乘积

正交曲线坐标下得常用关系式

气压梯度力

涡度

散度

风速矢量平流

绝对温度平流

普遍地图投影坐标系中得方程组

设X与Y轴地图投影放大系数为m与n,Z方向得地图投影放大系数为1

拉密系数

要求可以利用给出得关系式得到普遍地图投影坐标系中得大气方程组表达式

例如根据连续方程表达式与散度在正交曲线坐标系下得表达式,得到地图投影坐标下得连续方程表达式,要求将求与符号展开成各项相加得形式

第一步将H得表达式代入散度表达式,写出

第二步写出

第三步写出

差分方法与差分格式

离散化得概念

u（x,t）就是连续函数,u（i∆x,n∆t）就是u（x,t）经离散化后得形式。

所谓离散化,即把连续得x以i∆x代替,连续得t以n∆t代替,其中i与n为整数。

均就是以一维线性平流方程为例

差分格式基本都就是通过泰勒展开式来构造得。

前差格式:

后差格式:

中央差格式:

二阶微分得差分格式

拉普拉斯得差分格式

拉普拉斯得差分格式涉及到得格点

截断误差

上面差分格式中得R被成为截断误差。

意思就是用差商来近似代替偏微商时,将会因舍去R所代表得项而造成得误差。

一般用R中最大得项来表示截断误差得大小。

如果R中最大得项就是∆x,则R=O（∆x）。

如果就是∆x2,则R=O（∆x2）。

（注:

所谓最大得项指得就是偏导阶数最小得一项,一阶偏导项大于二阶偏导项）

R反映了差分方程代替微分方程时得截断误差,它在一定程度上代表了差分格式得精度,R得阶次越高,则差分格式得精度越高,误差越小。

这个得精度就就是R=O（∆x2）

相容性（一致性）

当空间步长∆x与时间步长∆t很小时,差分方程就是否逼近微分方程,这就就是差分格式得相容性（一致性）问题。

收敛性

在一定得定解条件下,差分方程得解就是否逼近微分方程得解得问题,称之为差分格式得收敛性问题。

稳定性

在时间积分过程中,由于舍入误差得影响,差分解得误差就是否随时间增长得问题,即差分格式得计算稳定性问题。

拉克斯（Lax）等价定理:

如果差分方程逼近微分方程,即差分格式与微分方程就是相容得,或者差分格式满足相容性条件,差分格式得稳定性,保证了其收敛性（计算稳定性就是收敛性得充分必要条件）。

用Von-Neumann稳定性判别方法来证明差分格式得计算稳定性时得主要步骤为:

1设解得波动形式,代入差分方程。

2得出其对应得增幅因子G。

3讨论时得情况。

4判断格式稳定性及满足格式稳定性得条件。

CFL判据

增幅因子G:

其中An+1与An分别就是n+1时刻与n时刻得振幅。

微分方程波动形式解:

差分方程形式解:

各个差分格式稳定性

时间前差,空间后差

条件稳定

时间前差,空间前差

条件稳定

时间前差,空间中央差

任何情况下不成立,所以为绝对不稳定

时间积分格式分类:

1、二时间层得积分格式（非迭代格式）

a、欧拉格式:

绝对不稳定

b、后差格式（隐式格式）,绝对稳定格式

c、梯形格式（隐式格式）,中性格式,振幅不变

2、二时间层得积分格式（迭代格式）

特点:

先通过前差格式算一个n+1时刻得粗略预报值,然后在通过这个预报值用后差格式计算n+1时刻得准确预报值

a、欧拉—后差格式（显式格式）:

条件稳定格式

b、赫恩格式:

绝对不稳定格式

3、三个时间层得积分格式

中央差格式,又称跳背格式

条件稳定

这种格式,可以瞧出不仅需要一个具有物理意义得初值u0,同样还需要一个出于计算要求得初值u1,前者称为物理初值,后者称为计算初值。

差分格式误差相关内容:

实际工作中,不可能任意缩小步长（由于计算量过大等原因）,实际计算就是在有限得网格下进行得,一定程度得误差就是不可避免得。

中央差得计算解问题:

使用中央差会产生两个波解,其中一个有物理意义,另一个不具有物理意义,就是计算过程中产生得虚假波形,称之为计算解。

中央差分格式得数值解为两个波动得叠加。

实际工作中通常采用提高网格分辩能力来抑制与减小计算解带来得影响。

时间得截断误差（频率误差）

当使用中央差格式（显式）时,原频率ω=f由替代

当使用梯形格式时（隐式）,原频率ω=f由替代

当∆t很小时,两种差分格式与ω=f差异很小

当∆t增加时,差分格式解偏离原频率幅度增大

f∆t>1时显式格式会出现不稳定,隐式格式得频率解随f∆t增加而减小

显式格式与隐式格式得频率解与真值得比值

空间得截断误差（波数误差）

用中央差格式展开时,波数k得数值解为:

差分近似精度随k或∆x得减小而增大,这也就就是说对于波长较短得波,其产生得波数误差较大;而波长较长得波,则差分方程可以比较精确地表示其空间微商,波数误差很小,精度较高。

相速度与群速度误差

空间差分格式得波数误差与时间积分格式得频率误差会造成相速度与群速度得误差,从而引起计算频散。

差分格式在波得移动与能量传播方面均可造成误差。

而且:

（1）由于相速度误差,减慢了平流过程;

（2）造成虚假得计算频散,且对短波尤为明显。

（1）波长越长,误差越小;波长减小,其误差也就更为严重;

（2）提高网格分辩率,使∆x取得足够小,可以提高相速度得准确率。

差分格式误差特征总结:

1、三层时间积分格式存在计算解问题:

计算解对差分解得影响依赖于网格分辨率与波长。

2、时间积分格式引起频率误差:

显式格式使其频率明显增加,振动加快;

隐式格式使其频率明显减小,振动减慢。

3、空间差分格式引起波数误差:

高阶差分格式所引起得波数误差要比低阶格式小;

波长较短得波,误差尤为严重。

4、空间差分格式会引起计算频散:

尤其对于短波,相速度与群速度均会产生很大得误差。

通常可采用提高网格分辩率得方法减小各种误差。

非线性不稳定

对于非线性偏微分方程,线性偏微分方程得稳定性条件,只能给出其计算稳定性得必要条件,即使满足这一条件,也可能会因为差分方程得边界条件与非线性项得不正确表示而产生计算得不稳定现象,我们把这种由于非线性作用而产生得不稳定,称为非线性不稳定。

混淆误差

差分方法就是用有限得自由度系统来代替原来得连续介质系统得。

而有限格点上得函数值只能分解有限得波数,其最短波长为2∆x,对于非线性作用产生得波长小于2∆x得波动,网格系统不能正确地将它表示出来,而把它错误地表示成为某种波长大于2∆x得波,从而产生了误差,我们把这种波得误差称之为混淆误差。

对于网格数为I得网格,假设ui中包含两个波,波数分别为k1与k2,当这两个波发生非线性作用而产生k1+k2得波时,如果k1+k2>I/2,则网格将把这个波得波数识别为S=I-（k1+k2）

自激反馈

假设k1与k2相互作用,能量反馈到上k1,则有:

k1=I-（k1+k2）,或2k1=I-k2。

由于k2<=I/2,所以k1>=I/4。

因此,发生能量反馈得波其波长必在2∆x与4∆x之间,也就就是说只有对波长很短得波才发生能量反馈;非线性不稳定得产生主要就是由于短波能量得虚假增长所造成得。

抑制与克服非线性计算不稳定得做法通常有:

▪空间与时间平滑,滤去短波分量（滤除波长小于4倍格距得波动）;

▪在方程中加入扩散项;

▪构造具有隐式平滑与某种选择性衰减作用得差分格式

▪构造守衡得差分格式,使差分方程尽可能保持原来得

物理规律与能量关系;

▪采用谱变换方法（可避免非线性不稳定问题）。

正压原始方程模式

以下在3个假设得基础上导出正压原