《第1章 大气运动的基本方程组PPT讲稿.ppt》由会员分享，可在线阅读，更多相关《第1章 大气运动的基本方程组PPT讲稿.ppt（60页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第1章 大气运动的基本方程组第1页，共60页，编辑于2022年，星期日 1.1运动学基础运动学基础一、标量场的空间变化一、标量场的空间变化1、位置矢量位置矢量：空间上的任一点空间上的任一点M(x，y，z)的位置则可用一个的位置则可用一个位置矢量表示：位置矢量表示：空间上位置的变化可用空间上位置的变化可用位置矢量位置矢量的改变量（位移矢量）的改变量（位移矢量）表示：表示：2第2页，共60页，编辑于2022年，星期日2、标量场的梯度、标量场的梯度 任一标量场任一标量场(以气压场以气压场p为例为例)可表为空间点和时间的函数：可表为空间点和时间的函数：考虑某一指定时刻考虑某一指定时刻(t=t0)气压气
2、压p在某一点的邻域的空间变化，在某一点的邻域的空间变化，则则p可视为只是空间变量的函数，其空间微分可表为：可视为只是空间变量的函数，其空间微分可表为：定义：定义：为气压梯度。为气压梯度。梯度算子（符）：梯度算子（符）：3第3页，共60页，编辑于2022年，星期日 气压气压p沿沿 （或（或）方向）方向的方向导数可表为：的方向导数可表为：此式清楚地表明了气压空间此式清楚地表明了气压空间变化与气压梯度的关系：变化与气压梯度的关系：当当 0时，时，即即p的方向导数取得的方向导数取得最大正值；当最大正值；当 /2时，时，4第4页，共60页，编辑于2022年，星期日三、三、场变量的时间变化场变量的时间变化
3、1、局地时间变化率、局地时间变化率 当考察在空间某个固定点（位置矢为当考察在空间某个固定点（位置矢为 ）上一个场变量随时间）上一个场变量随时间t的变化时，所测得的变化时，所测得(或观测到或观测到)的变化称为该场变量在该地点上的局的变化称为该场变量在该地点上的局地（时间）变化。地（时间）变化。场变量场变量F（，t）在点）在点 上的局地变化率上的局地变化率(单位时间内的变化量单位时间内的变化量)可定量地表为：可定量地表为：2、个别时间变化率、个别时间变化率 个别变化率是指跟随某个个别变化率是指跟随某个“动点动点”（如移动的飞机、车、（如移动的飞机、车、5第5页，共60页，编辑于2022年，星期日
4、船、空气质点或天气系统中的特性点等）在运动过程中所历经（或船、空气质点或天气系统中的特性点等）在运动过程中所历经（或测得）某物理量测得）某物理量F 随时间的变化率。随时间的变化率。其数学表达式可写为：其数学表达式可写为：与局地变化率不同，它是物理量在不同地点、不同时刻的变化与局地变化率不同，它是物理量在不同地点、不同时刻的变化率。率。3、平流变化率、平流变化率 改写个别变化率的表达式：改写个别变化率的表达式：6局地时间变化率局地时间变化率第6页，共60页，编辑于2022年，星期日右边第二项的分子可表为：右边第二项的分子可表为：取取 (同时有同时有 )的极限：的极限：平流变化率平流变化率 其中：
5、其中：为动点的为动点的位置矢的时间变化率，位置矢的时间变化率，也称为平流速度。也称为平流速度。F的个别变化率等于其局地变化率与平流变化率之的个别变化率等于其局地变化率与平流变化率之和。和。当动点就是空气质点时，气象上通常用当动点就是空气质点时，气象上通常用 表示空气运动的速表示空气运动的速度：度：7第7页，共60页，编辑于2022年，星期日其中：其中：这样有：这样有：上式可以看成是上式可以看成是“个别微分算子个别微分算子”：作用于场变量作用于场变量F F 的结果。若令的结果。若令F=TF=T,T T为气温为气温,则由上式有则由上式有：这是局地温度的预报方程。左边代表局地的温度变化率，右边的项这
6、是局地温度的预报方程。左边代表局地的温度变化率，右边的项可视为影响局地温度变化的强迫因子。可视为影响局地温度变化的强迫因子。8第8页，共60页，编辑于2022年，星期日上式右边第二项（上式右边第二项（）称为温度平流。）称为温度平流。当当 0 时时,称为暖平流称为暖平流，可造成升温：，可造成升温：T-T+T-T+冷平流冷平流 暖平流暖平流例题例题20第9页，共60页，编辑于2022年，星期日三、速度场的散度和涡度三、速度场的散度和涡度1、速度散度和连续方程、速度散度和连续方程1 1）速度散度）速度散度考虑表面积为考虑表面积为S、体积为、体积为的空气块的空气块（如图），由于其表面上各点的速（如图）
7、，由于其表面上各点的速度分布不均匀而引起的体积变化率度分布不均匀而引起的体积变化率 高斯公式高斯公式其中：其中：速度散度速度散度第10页，共60页，编辑于2022年，星期日考虑气块体积趋于零有考虑气块体积趋于零有：速度散度的物理意义速度散度的物理意义：空气微团体积的相对变化率。空气微团体积的相对变化率。当垂直速度为零时，当垂直速度为零时，空气运动为水平运动，空气微团的体空气运动为水平运动，空气微团的体积变化率退化为水平面积积变化率退化为水平面积(A)的变化率：的变化率：辐散辐散辐合辐合第11页，共60页，编辑于2022年，星期日2）连续方程）连续方程3）气压倾向方程（）气压倾向方程（P10）第
8、12页，共60页，编辑于2022年，星期日2、速度场的涡度、速度场的涡度 1）涡度：是用来描述空气微团的旋转特性：涡度：是用来描述空气微团的旋转特性：x分量分量y分量分量z分量分量第13页，共60页，编辑于2022年，星期日第14页，共60页，编辑于2022年，星期日 对于大尺度运动，垂直方向的涡度分量是主要的，天气学上常对于大尺度运动，垂直方向的涡度分量是主要的，天气学上常常主要考虑垂直涡度分量常主要考虑垂直涡度分量z，并且约定：并且约定：在北半球：在北半球：0，称之为气旋式涡度，称之为气旋式涡度，0(0(C 0)时，气象上称之为气旋式（反气旋时，气象上称之为气旋式（反气旋式）环流式）环流。
9、3）速度环流与涡度的关系）速度环流与涡度的关系 计算沿围线计算沿围线ABCD的速度环流：的速度环流：AB:BC:CD:DA:第16页，共60页，编辑于2022年，星期日沿围线沿围线ABCD的总速度环流则为的总速度环流则为:由此有由此有:为矩形中心点处的铅直涡度分量，为矩形为矩形中心点处的铅直涡度分量，为矩形ABCD的面积。当矩形面积趋于零，取极限则有的面积。当矩形面积趋于零，取极限则有 可见，铅直涡度分量可见，铅直涡度分量可解释为水平围线上的速度环流在面积趋于可解释为水平围线上的速度环流在面积趋于零时的极限，或者说是单位面积上的速度环流。零时的极限，或者说是单位面积上的速度环流。第17页，共6
10、0页，编辑于2022年，星期日1.2 旋转坐标系中的大气运动方程旋转坐标系中的大气运动方程1惯性坐标系与非惯性坐标系惯性坐标系与非惯性坐标系 牛顿第二定律只适用于某种特定的坐标系（或参照系）。牛顿第二定律只适用于某种特定的坐标系（或参照系）。按牛顿第二定律是否成立，可将坐标系分成：按牛顿第二定律是否成立，可将坐标系分成：绝对绝对 （静止）坐标系（静止）坐标系：能使牛顿第二定律成立的坐标系。在这种坐：能使牛顿第二定律成立的坐标系。在这种坐标系中，牛顿惯性定律亦成立，故又称之为标系中，牛顿惯性定律亦成立，故又称之为惯性坐标系惯性坐标系。相对于相对于惯性坐标系作匀速直线运动的坐标系仍是惯性坐标系。惯
11、性坐标系作匀速直线运动的坐标系仍是惯性坐标系。相对坐标系相对坐标系：相对于惯性坐标系作加速运动的坐标系，也称：相对于惯性坐标系作加速运动的坐标系，也称非非惯性坐标系惯性坐标系。第18页，共60页，编辑于2022年，星期日在气象学上：在气象学上：1）通常将相对于恒星静止、不随地球自转的坐标系称为绝对坐标）通常将相对于恒星静止、不随地球自转的坐标系称为绝对坐标系（惯性坐标系或系（惯性坐标系或“静止静止”坐标系），在绝对坐标系中观测到坐标系），在绝对坐标系中观测到的大气运动称为的大气运动称为绝对运动绝对运动；并且，通常略去地球绕太阳公转；并且，通常略去地球绕太阳公转引起的加速度引起的加速度(610-
12、3m/s2)；2）将固定于地球上、跟随地球自转一起转动的坐标系称为相对坐标）将固定于地球上、跟随地球自转一起转动的坐标系称为相对坐标系（旋转坐标系），它是一种非惯性坐标系系（旋转坐标系），它是一种非惯性坐标系,在此坐标系中观测在此坐标系中观测到的大气运动称为到的大气运动称为相对运动相对运动。2 惯性坐标系中的运动方程惯性坐标系中的运动方程 在惯性坐标系中，按牛顿第二运动定律，单位质量空气徽团在惯性坐标系中，按牛顿第二运动定律，单位质量空气徽团的运动方程可表的运动方程可表 第19页，共60页，编辑于2022年，星期日 上式是在绝对坐标系中，单位质量空气微团所遵从的运动方程，上式是在绝对坐标系中，
13、单位质量空气微团所遵从的运动方程，有时称为有时称为绝对运动方程绝对运动方程。但是，由于在地球坐标系中（例如地。但是，由于在地球坐标系中（例如地球上的测站）无法直接观测到球上的测站）无法直接观测到绝对速度绝对速度和和绝对加速度绝对加速度，只能观，只能观测到测到相对速度和相对加速度相对速度和相对加速度。因此，上式并不能直接用于研究。因此，上式并不能直接用于研究地球大气运动。地球大气运动。找出绝对速度与相对速度以及绝对加速度与相对加速度的找出绝对速度与相对速度以及绝对加速度与相对加速度的关系？关系？3 两种坐标系中的速度和加速度的关系两种坐标系中的速度和加速度的关系 第20页，共60页，编辑于202
14、2年，星期日绝对位移绝对位移相对位移相对位移第21页，共60页，编辑于2022年，星期日 若将由于地球自转引起若将由于地球自转引起p点的移动速度（称为牵连速度）记为点的移动速度（称为牵连速度）记为 ，则，则p点的牵连位移为点的牵连位移为 空气微团的绝对位移等于其相对位移与牵连位移之向量和空气微团的绝对位移等于其相对位移与牵连位移之向量和 若用若用 除上式两端，并取除上式两端，并取 趋于零的极限，则有趋于零的极限，则有 即即 绝对速度与相对绝对速度与相对速度的关系速度的关系第22页，共60页，编辑于2022年，星期日位于纬度位于纬度 处的空气处的空气质点的牵连速度就是质点的牵连速度就是该质点随地
15、球自转时该质点随地球自转时在纬圈平面上以角速在纬圈平面上以角速度度 作匀速圆周运动作匀速圆周运动的线速度的线速度：（推导）（推导）于是有：于是有：和和 从个别变化率的定义出发，可直接证明，对于任一标量从个别变化率的定义出发，可直接证明，对于任一标量F,绝对坐绝对坐标系中的个别变化率等于相对坐标系中的个别变化率：标系中的个别变化率等于相对坐标系中的个别变化率：空气质点的牵连速度空气质点的牵连速度作业：习题作业：习题21，29第23页，共60页，编辑于2022年，星期日证明：对任意矢量证明：对任意矢量 ，成立。成立。第24页，共60页，编辑于2022年，星期日令：令：有：有：（推导）（推导）即表述
16、即表述绝对加速度绝对加速度与与相对加速度相对加速度关系的定量关系式。关系的定量关系式。绝对加速绝对加速度度等于等于相对加速度相对加速度加上两个由于坐标系旋转而引起的加上两个由于坐标系旋转而引起的附加加速度：附加加速度：1）科里奥利（科里奥利（Coriolis）加速度：）加速度：2）向心加速度：）向心加速度：通过已经确定了旋转坐标系中的通过已经确定了旋转坐标系中的相对速度与绝对速度相对速度与绝对速度以及以及相对相对加速度与绝对加速度加速度与绝对加速度的定量关系，的定量关系，那么我们就完全可以通过那么我们就完全可以通过地球上探测到的风速（相对速度）来定量地表述绝对速度和绝地球上探测到的风速（相对速
17、度）来定量地表述绝对速度和绝对加速度。因而，我们可以进一步导出便于直接用于研究地球对加速度。因而，我们可以进一步导出便于直接用于研究地球大气运动规律的运动方程大气运动规律的运动方程旋坐标系中的运动方程（相对运旋坐标系中的运动方程（相对运动方程）。动方程）。第25页，共60页，编辑于2022年，星期日4、相对运动方程相对运动方程 在等号的不同边，在等号的不同边，或称或称“惯性力惯性力”，或称或称 “加速度加速度”第26页，共60页，编辑于2022年，星期日5、作用于空气微团上的作用力、作用于空气微团上的作用力1）气压梯度力）气压梯度力 考虑右图空气考虑右图空气块所受压力：块所受压力：于是，对单位
18、质量的空气小体而言，在于是，对单位质量的空气小体而言，在x方向上所受压力的合力为方向上所受压力的合力为：第27页，共60页，编辑于2022年，星期日Y方向的合力：方向的合力：Z方向的合力：方向的合力：综合起来，单位质量空气小体所受总的压力合力即气压梯度综合起来，单位质量空气小体所受总的压力合力即气压梯度力为力为气压梯度力的性质：气压梯度力的性质：1 1）气压梯度力的方向与气压梯度的方向相反，即与等压面）气压梯度力的方向与气压梯度的方向相反，即与等压面（线）垂直、指向气压降低的方向；（线）垂直、指向气压降低的方向；2 2）气压梯度力的大小与气压梯度的大小成正比，与空气密度）气压梯度力的大小与气压
19、梯度的大小成正比，与空气密度成反比成反比。第28页，共60页，编辑于2022年，星期日2)科里奥利（科里奥利（Coriolis）力）力:它的存在条件：它的存在条件：0（旋转）；（旋转）；0，即有相对于地球的运动。，即有相对于地球的运动。科氏力（地转偏向力）的性质：科氏力（地转偏向力）的性质：i)即科氏力一定在纬圈平面上；即科氏力一定在纬圈平面上；ii)科氏力只会改变运动的方向，不改变其大科氏力只会改变运动的方向，不改变其大小，即对空气微团不做功。小，即对空气微团不做功。iii)对于北半球的水平运动，科氏力总是指向运动前进方对于北半球的水平运动，科氏力总是指向运动前进方向的右方（观测者面向运动前
20、进方向），南半球的情形则相向的右方（观测者面向运动前进方向），南半球的情形则相反，指向运动前进方向的左方。反，指向运动前进方向的左方。第29页，共60页，编辑于2022年，星期日例如，例如，在北半球，向南（北）流的水流会受到指向西（东）的科氏在北半球，向南（北）流的水流会受到指向西（东）的科氏力的作用，引起水流向南（北）流的河床的西（东）岸受到更为力的作用，引起水流向南（北）流的河床的西（东）岸受到更为严重的冲刷严重的冲刷。又如，。又如，在北半球运行的远程火箭，当它铅直上升在北半球运行的远程火箭，当它铅直上升（下降）时，其轨道要向西（东）偏移；当它在水平面方向（下降）时，其轨道要向西（东）偏移
21、；当它在水平面方向向东（西）飞行时，其轨道要向南并向上（向北并向下）偏向东（西）飞行时，其轨道要向南并向上（向北并向下）偏移。移。3 3 重力重力：地心引力与惯性离心力的合力地心引力与惯性离心力的合力单位质量空气微团所受地心引力可表为单位质量空气微团所受地心引力可表为 第30页，共60页，编辑于2022年，星期日 地心引力是指向地心的有势力，设地心引力是指向地心的有势力，设 为地心引力势，则地心为地心引力势，则地心引力可表为：引力可表为：势力的性质势力的性质：沿闭合路线积分为零，无旋。（证明）：沿闭合路线积分为零，无旋。（证明）若假定极地海平面（若假定极地海平面（）上的地心引力势为零）上的地心
22、引力势为零 （），则可以求得地心引力势为：），则可以求得地心引力势为：物理图像：物理图像：地心引力的等势面为以地心为球心的同心球面族地心引力的等势面为以地心为球心的同心球面族。第31页，共60页，编辑于2022年，星期日 惯性离心力惯性离心力：垂直于地轴、从地轴指向外，也是一种有垂直于地轴、从地轴指向外，也是一种有势力。设其势函数为势力。设其势函数为 ，则可将惯性离心力表为：，则可将惯性离心力表为：（证明）（证明）若假设地轴上（若假设地轴上（）的惯性离心力势为零（）的惯性离心力势为零（），则可求），则可求得得 物理图像：物理图像：惯性离心力的等势面为以地轴为轴的同轴圆柱面族。惯性离心力的等势面
23、为以地轴为轴的同轴圆柱面族。这样，重力加速度可表为这样，重力加速度可表为：重力位势重力位势：第32页，共60页，编辑于2022年，星期日 由地心引力和惯性离心力的性质可知，由地心引力和惯性离心力的性质可知，重力的大小重力的大小随所随所处的处的纬度纬度和和高度高度不同而不同。不同而不同。在同一海拔高度上，重力加速度在同一海拔高度上，重力加速度随纬度增大而增大，赤道上最小，极地最大；在同一纬度，重随纬度增大而增大，赤道上最小，极地最大；在同一纬度，重力加速度随海拔高度增大而减小力加速度随海拔高度增大而减小。计算重力的近似公式（。计算重力的近似公式（1930年国际重力公式）可表为：年国际重力公式）可
24、表为：这里这里 为纬度，为纬度，z为海拔高度（以米为单位）。为海拔高度（以米为单位）。由于重力加速度随纬度和高度的变化很小，它与地心引力的由于重力加速度随纬度和高度的变化很小，它与地心引力的 夹角也很小，故气象上一般将重力加速度视为常数，取夹角也很小，故气象上一般将重力加速度视为常数，取g9.81 米米/秒秒2；并近似地视地球为半径；并近似地视地球为半径a6371公里公里的球形，视等重力位的球形，视等重力位势面为同心球面族。势面为同心球面族。单位？单位？第33页，共60页，编辑于2022年，星期日重力位势则可近似估算为（假定重力位势则可近似估算为（假定 ）：）：（推导）（推导）4 分子粘性力（
25、内摩擦力）分子粘性力（内摩擦力）其中，其中，为运动学粘性系数，为运动学粘性系数，为动力粘性系数。为动力粘性系数。大气是一种低粘流体，除了贴近地面几厘米厚度的薄层，因为空大气是一种低粘流体，除了贴近地面几厘米厚度的薄层，因为空气运动速度垂直梯度很大，必须考虑分子粘性作用的影响外，一般气运动速度垂直梯度很大，必须考虑分子粘性作用的影响外，一般都可忽略分子粘性力的作用。都可忽略分子粘性力的作用。第34页，共60页，编辑于2022年，星期日最后，矢量形式的相对运动方程可改写为最后，矢量形式的相对运动方程可改写为 第35页，共60页，编辑于2022年，星期日1.3 运动方程的分量形式运动方程的分量形式
26、1 球坐标系：球坐标系：原点原点O却在地心却在地心 经度；度；纬度度 r 距球心距离距球心距离 为沿纬圈指向东为沿纬圈指向东 沿经圈指向北沿经圈指向北 沿径向指向天顶的沿径向指向天顶的 单位向量。单位向量。第36页，共60页，编辑于2022年，星期日沿各方向的坐标线元分别为：沿各方向的坐标线元分别为：在球坐标中，若将空气微团的运动速度矢表为在球坐标中，若将空气微团的运动速度矢表为则有：则有：，第37页，共60页，编辑于2022年，星期日2 球坐标系的个别微分算子球坐标系的个别微分算子 对任意一个场变量对任意一个场变量 ，有：，有：两边除以两边除以dt，取极限，并利用球坐标各速度分量的表达式有：
27、，取极限，并利用球坐标各速度分量的表达式有：所以：球坐标个别微分算子：所以：球坐标个别微分算子：第38页，共60页，编辑于2022年，星期日梯度算子：梯度算子：要求得运动方程在球坐标系中的分量方程，就得先将运动方程要求得运动方程在球坐标系中的分量方程，就得先将运动方程中向量形式的各项分别分解到球坐标系的各个方向上，即求出各中向量形式的各项分别分解到球坐标系的各个方向上，即求出各项的分量式。项的分量式。3 球坐标系加速度的分量形式球坐标系加速度的分量形式 (1)相对加速度可写为：相对加速度可写为：问题归结为求单位向量问题归结为求单位向量 、和和 的个别变化率的分量表达的个别变化率的分量表达式。式
28、。第39页，共60页，编辑于2022年，星期日 由于由于 所以：所以：考虑单位矢量考虑单位矢量 随经度随经度的变化：的变化：第40页，共60页，编辑于2022年，星期日 大小：大小：于是：于是：的方向（即的方向（即 的方向）与的方向）与 方向相反，指向地轴。方向相反，指向地轴。第41页，共60页，编辑于2022年，星期日于是，有于是，有 （还有简便的推导还有简便的推导）第42页，共60页，编辑于2022年，星期日其中，包含因子其中，包含因子 1/r 的项是由于地球的球面性引起的曲率项，的项是由于地球的球面性引起的曲率项，称为称为“曲率加速度曲率加速度”。4 各作用力在球坐标系中的分量各作用力在
29、球坐标系中的分量(1)气压梯度力的分量形式气压梯度力的分量形式 利用球坐标系中的梯度算子，可得气压梯度力在球坐标系中分解式利用球坐标系中的梯度算子，可得气压梯度力在球坐标系中分解式为为:(2)(2)科氏力科氏力 其中，其中，地转角速度可表为：地转角速度可表为：第43页，共60页，编辑于2022年，星期日各分量可分别表为：各分量可分别表为：于是，在球坐标系中，科氏力可表为于是，在球坐标系中，科氏力可表为其中其中称为科氏参数。称为科氏参数。第44页，共60页，编辑于2022年，星期日（3）重力和分子粘性力）重力和分子粘性力重力只在垂直方向有分量重力只在垂直方向有分量：分子粘性力可形式上表为分子粘性
30、力可形式上表为：5 球坐标系的运动方程的分量形式球坐标系的运动方程的分量形式，第45页，共60页，编辑于2022年，星期日 等式右边的含因子等式右边的含因子1/r的项源于球坐标系中的的项源于球坐标系中的“曲率加速度曲率加速度”，可称为可称为“曲率惯性力曲率惯性力”，可以证明该力与速度垂直，即其对可以证明该力与速度垂直，即其对空气微团不做功空气微团不做功。由于。由于90%以上的大气都集中在离地面以上的大气都集中在离地面20公里以公里以下的薄层内，有时采用郭晓岚称谓的下的薄层内，有时采用郭晓岚称谓的“薄层近似薄层近似”，即在上述方程，即在上述方程中，当中，当r 作为系数出现时，近似地取作为系数出现
31、时，近似地取 于是薄层近似下的运动方程可表为于是薄层近似下的运动方程可表为:第46页，共60页，编辑于2022年，星期日6 球坐标系的速度散度和涡度球坐标系的速度散度和涡度球坐标系中的单位向量的基本微分公式可表为：球坐标系中的单位向量的基本微分公式可表为：及及 第47页，共60页，编辑于2022年，星期日利用上述微分关系，可得到球坐标系中速度场的散度：利用上述微分关系，可得到球坐标系中速度场的散度：类似地，速度场的涡度可表为类似地，速度场的涡度可表为球坐标系中的拉普拉斯算子可表为球坐标系中的拉普拉斯算子可表为 第48页，共60页，编辑于2022年，星期日6 球坐标系得连续方程球坐标系得连续方程
32、 将球坐标系中的散度表达式代入上连续方程，可将球坐标系中将球坐标系中的散度表达式代入上连续方程，可将球坐标系中的质量连续方程写为：的质量连续方程写为：或：或：另一种推导方法：另一种推导方法：小体积得空气质量增加小体积得空气质量增加净的质量流入净的质量流入 第49页，共60页，编辑于2022年，星期日 1.4 局地直角坐标系中的运动方程与连续方程局地直角坐标系中的运动方程与连续方程 1 局地直角坐标系局地直角坐标系 局地直角坐标系的坐标局地直角坐标系的坐标原点原点位于海平面上指定地点，位于海平面上指定地点，x x轴轴沿纬线指向东，沿纬线指向东，y y轴沿经线指轴沿经线指向北，向北，z z轴指向天
33、顶轴指向天顶。三个坐标轴向的单位向量分三个坐标轴向的单位向量分别为和，它们与球坐标系中沿别为和，它们与球坐标系中沿和方向的单位矢的指向完全相和方向的单位矢的指向完全相同，所不同的是，现在，在坐标原点附近的同，所不同的是，现在，在坐标原点附近的“局部地区局部地区”（研究问（研究问题的区域），单位矢量将视为题的区域），单位矢量将视为不随地点的改变而改变的单位矢量不随地点的改变而改变的单位矢量。（当范围不大时，即局地看成一个平面，而不是球面）（当范围不大时，即局地看成一个平面，而不是球面）第50页，共60页，编辑于2022年，星期日2 局地直角坐标系中的运动方程与连续方程局地直角坐标系中的运动方程与
34、连续方程 在球坐标系的基础上，假定：在球坐标系的基础上，假定：1）略去球面性所产生的曲率项。）略去球面性所产生的曲率项。2）形式上，令）形式上，令空气运动的速度分量可表为：空气运动的速度分量可表为：，第51页，共60页，编辑于2022年，星期日局地直角坐标系中的运动方程可表为：局地直角坐标系中的运动方程可表为：局地直角坐标系中的连续方程为局地直角坐标系中的连续方程为：或或 第52页，共60页，编辑于2022年，星期日3-平面近似平面近似 采采用用局局地地直直角角坐坐标标系系，忽忽略略地地球球的的球球面面性性，可可以以使使方方程程形形式式变变得得简简单单，但但是是，以以后后将将会会看看到到，如如
35、果果完完全全略略去去地地球球的的球球面面性性，取取科科氏氏参参数数为为常常数数（即即所所谓谓的的f-f-平平面面近近似似），除除了了精精度度降降低低外外，还还会会带带来来一一个个严严重重的的动动力力学学缺缺陷陷，即即消消除除了了气气象象学学中中非非常常重重要要的的一一种种波波动动罗罗斯斯贝贝（Rossby）波波的的生生存存条条件件。为为了了弥弥补补这这种种缺缺陷陷，可采用所谓的可采用所谓的-平面近似平面近似，部分地保留地球球面性的影响。，部分地保留地球球面性的影响。设设局局地地直直角角坐坐标标系系原原点点所所在在纬纬度度为为 ，科科氏氏参参数数 f 可可在在的的邻域展为如下泰勒级数：邻域展为如
36、下泰勒级数：其中，其中，第53页，共60页，编辑于2022年，星期日若只保留泰勒级数右边第一项，略去其它项，则得如下最低阶若只保留泰勒级数右边第一项，略去其它项，则得如下最低阶近似：近似：常数常数 此近似称为此近似称为“-平面近似平面近似”。当取泰勒级数右边第一和第二项，略去其它项时，有：当取泰勒级数右边第一和第二项，略去其它项时，有：其中，其中，。以上近似即所谓的。以上近似即所谓的“-平面近似平面近似”。球坐标系和局地直角坐标系的比较球坐标系和局地直角坐标系的比较第54页，共60页，编辑于2022年，星期日1.5 大气运动的闭合方程组与初、边条件大气运动的闭合方程组与初、边条件 1 闭合方程
37、组闭合方程组 闭合方程组闭合方程组：独立方程的个数与未知函数的个数相同的方程组。：独立方程的个数与未知函数的个数相同的方程组。描述大气运动的基本方程有：描述大气运动的基本方程有：运动方程（含三个分量方程）、连续运动方程（含三个分量方程）、连续方程、热力学方程、状态方程和水份守恒方程等共七个方程方程、热力学方程、状态方程和水份守恒方程等共七个方程：第55页，共60页，编辑于2022年，星期日 并非是一个精确的已知量；并非是一个精确的已知量；为对单位质量空气的为对单位质量空气的非绝热加热率，包括辐射加热、水的相变加热和分子的或湍非绝热加热率，包括辐射加热、水的相变加热和分子的或湍流的热传导等。每一
38、种非绝热过程都是非常复杂的过程，至流的热传导等。每一种非绝热过程都是非常复杂的过程，至少目前人们还不能用严格而精确的解析式来表达它们，即这少目前人们还不能用严格而精确的解析式来表达它们，即这一项也并不是严格意义上的已知函数。一项也并不是严格意义上的已知函数。所以，要闭合上述方程组的主要困难是缺乏精确的表达式来描述所以，要闭合上述方程组的主要困难是缺乏精确的表达式来描述大气中的重要源、汇项大气中的重要源、汇项 和和 等。在本课程范围内，我们将主等。在本课程范围内，我们将主要考虑干空气（要考虑干空气（=0）的绝热（）的绝热（=0）过程，这时，上述方）过程，这时，上述方程组变为包含六个方程和六个未知
39、函数，程组变为包含六个方程和六个未知函数，u,v,w，p，T,和和 的闭的闭合方程组。合方程组。第56页，共60页，编辑于2022年，星期日2 初始条件与边界条件初始条件与边界条件 初始条件即未知函数初始条件即未知函数 ，P，T和和等在初始时刻等在初始时刻()的值的值:边界条件则比较复杂，与边界的几何与物理性质有关。对边界条件则比较复杂，与边界的几何与物理性质有关。对于全球大气，主要是考虑铅直方向上的于全球大气，主要是考虑铅直方向上的上、下边界条件上、下边界条件。第57页，共60页，编辑于2022年，星期日1)下边界条件下边界条件 a)无地形：无地形：（粘性大气）（粘性大气）（运动学边界）（运
40、动学边界）（法向速度为零）（法向速度为零）（动力学边界）：（动力学边界）：b)有地形：则运动学和动力学下边界条件可分别表为有地形：则运动学和动力学下边界条件可分别表为 第58页，共60页，编辑于2022年，星期日2）上边界条件）上边界条件 在大气上边界处（在大气上边界处（，任一物理量（气象要素）都应，任一物理量（气象要素）都应有界。常见的条件有：假定在大气上界处，空气密度（有界。常见的条件有：假定在大气上界处，空气密度（）、）、风速（风速（V）、大气压力（）、大气压力（P）等为零及大气与外界没有动量和动）等为零及大气与外界没有动量和动能的交换。能的交换。3）自由边界条件）自由边界条件 如果假定大气是均质（密度为常数）大气，则模式大气顶是一个自如果假定大气是均质（密度为常数）大气，则模式大气顶是一个自由面：由面：与外界无质量交换与外界无质量交换与外界无动能交换与外界无动能交换自由面高度第59页，共60页，编辑于2022年，星期日假定在自由面上的压力为常数（假定在自由面上的压力为常数（），则自由面条件可表为：），则自由面条件可表为：和和 如果研究的不是全球大气，而是某个有限区域，则还须适当给如果研究的不是全球大气，而是某个有限区域，则还须适当给定区域侧边界上的条件，称为侧边界条件。定区域侧边界上的条件，称为侧边界条件。第60页，共60页，编辑于2022年，星期日