第3章大气边界层支配方程之1.ppt

《第3章大气边界层支配方程之1.ppt》由会员分享，可在线阅读，更多相关《第3章大气边界层支配方程之1.ppt（47页珍藏版）》请在课桌文档上搜索。

1、湍流运动的两种研究方法,1、采用随机过程的统计学方法来反映大气湍流结构 平均、方差、标准差、湍强、相关、通 量、应力、湍流动能 2、解湍流运动控制方程（平均运动方程、脉动方程、湍能方程.）,第三章 大气边界层支配方程,3.1 基本控制方程 3.2 平均量方程 3.3 湍流脉动量方程 3.4 湍流方差预报方程 3.5 湍流通量预报方程 3.6 闭合理论,3.1 基本控制方程,为了定量描述和预报边界层状况，需要借助于流体力学来描述大气中气体动力学和热力学方程。描述气体和液体流动的方程组包括：,(水汽及污染物浓度的标量方程类似热量方程),3.1 基本控制方程,3.1.1 控制方程3.1.2 简化与近

2、似Boussinesq近似3.1.3 坐标系及其变换,3.1.1 控制方程,状态方程质量守恒（连续方程）动量守恒（牛顿第二定律）热量守恒方程（热力学第一定律）,1）状态方程,理想气体状态方程：,P：气压：湿空气密度R：干空气气体常数（R=287 J K-1 kg-1）虚温Tv=T(1+0.61q)，q比湿，T温度,2）质量守恒（连续方程）,连续性方程的一般形式：,不可压缩流体,或,泰勒假说:,运用Einstein求和符号惯例，以上的连续方程写成：,（j=1、2、3分别代表x、y、z三个方向）,3）动量守恒（牛顿第二定律）,动量方程的表达式：,局地时间变化,平流项,气压梯度力项,重力作用项,科氏

3、效应,粘滞应力项,i3为克罗内克符号,重力仅在垂直方向起作用,思考：方程右端第三项展开？,4）热量守恒（热力学第一定律）,热量守恒关系常用位温方程来表示：,s：引起空气位温变化的热量源(汇)项（分子热传导、与相变有关的潜热释放、净辐射加热）,水汽及污染物的守恒方程形式与热量守恒形式一致 关键是要全面准确的了解引起成份变化的源汇项,3.1.2 简化与近似Boussinesq近似,在一定条件下，控制方程中某些项的量值比其它项小得多，以致可以把它略去，使得方程变得较为简单，有助于方程的求解。,必须考虑地球自转的影响,大气密度并非均匀，主要在垂直方向上是不均匀的层结流体,大气的空间尺度在水平方向远大于

4、铅直方向，可视作浅层流体，不可压缩流体,大气边界层主要是湍流运动,大气边界层的特点概括有：,Boussinesq 近似,Boussinesq 近似的基本假定：,流体中的动力学粘滞系数=是常数,流体中的分子导热系数 kT 是常数,大气是浅层流体，垂直范围约10km,描写流体热力状态的特征量可表示为,扰动量远小于基态量：,基本状态量假设是静力平衡、绝热的，并且满足理想气体状态方程，即：,Boussinesq 近似下的简化方程：,连续方程状态方程运动方程热流量方程,1）连续方程,当大气为浅层流体，在该范围内密度的变化很小、可以忽略，边界层大气可以近似认为是不可压缩的，连续方程为：,或,2）状态方程,

5、对于干空气，状态方程为 P=RT，对该式进行对数微分，可以得到：,如果运动的垂直尺度相对于大气层厚度而言很小，则,取对数：,取微分：,3）运动方程,Navier-Stokes方程中压力梯度项和重力项可以进行改写：,将上式代入运动方程得到：,扰动温度在重力作用下形成的净浮力项,4）热流量方程,边界层支配方程中常用到位温，Poisson方程：,对该式进行对数微分，近似有,此外，有,因此，有,扰动温度的垂直梯度,位温的垂直梯度,如果不考虑相变，热源 S 仅考虑分子热传导及辐射加热，则热流量方程为：,Rj:j方向的辐射热通量,分子热传导的加热,辐射加热,上式中Td 也可以换成,以上近似处理最早由Bou

6、ssinesq(1903)提出Boussinesq近似该简化方程假定流体不可压、并限制在一薄层内。适用于研究像积云对流、海陆风环流、边界层急流中的重力波活动等发生在浅层内的中尺度运动浅水方程,综上所述，Boussinesq近似下的基本方程组为：,或,或,热流量方程,运动方程,状态方程,连续方程,3.1.2 坐标系及其变换,通常我们使用笛卡尔坐标系，在实际应用中，将笛卡尔坐标系绕z轴旋转，使x、y轴指向其它方向，能够在处理问题时更方便。例如，使x轴与平均风向、地转风向、表面应力方向、垂直于海岸线、山的方向一致，这样就可以简化控制方程中的某些项。例如，选择x轴与平均风方向一致，就可以求得u=M（平

7、均风速）和v=0，这时，x轴叫做顺风方向，y轴叫做侧风方向。,3.2 平均量方程,以上动力方程组无法求解析解，但可以求得数值解。原则上，可以用动力方程组来描述湍流的运动，但是要想囊括所有尺度的湍流运动，计算量太大。为了简化，可以截取一定尺度的涡旋，而在这个尺度以下的涡旋用湍流的统计特征来代替。在一些中尺度和天气尺度模式中，截取尺度为10100公里，而在边界层模式，比如说大涡模拟，截取尺度一般为10100米。,一 出发方程组,1.状态方程2.连续性方程3.动量守恒4.热量守恒5.水汽守恒6.标量守恒,1、状态方程,干空气气体常数,2、连续方程,张量展开：,3、运动方程(动量守恒),存储项,平流传

8、输项,重力项，仅在垂直方向作用,柯氏力项,气压梯度力项,粘性力项,式中：为分子动粘系数，f=2sin。,4、热量守恒,存储项,平流传输项,热扩散项,式中：k为分子热扩散系数，数值为2.0610-5 m2s-1；Lp为与E相变有关的潜热（0C时气液相变取值2.50106 Jkg-1；液固相变取值3.34105Jkg-1;气固相变取值2.83106 Jkg-1）；cp为湿空气定压比热，与干空气定压比热的关系为cp=cpd(1+0.86q)；cpd取值1004.07 Jkg-1K-1；E为蒸发量,辐射加热项,与相变有关的潜热释放,5、水汽守恒,存储项,平流传输项,水汽扩散项,式中：vq为空气中水汽分

9、子扩散率，Sq是方程中不含有其余过程时的净水体源项（源-汇），单位是：单位时间单位体积的总水体质量。,其它过程的净水体源项,6、标量守恒,存储项,平流传输项,C的扩散项,其它过程的体源项,式中：vc为空气中C的分子扩散率，Sc是不存在于方程中的其余过程的体源项，例如化学反应等。,在湍流运动的大气边界层中，上述方程组还不能完整地描述边界层中的全部过程，应将上述的主要变量转换成平均量和脉动量相加。即：,平均场方程描述长时间过程，脉动场方程描述短时间过程。,二 湍流中平均变量方程,推导思路：,出发方程：Boussinesq 近似方程组 采用雷诺平均的方法，将任意一个物理量表示成平均量和脉动量之和，代

10、入方程组，然后再取平均。大气边界层平均量控制方程,雷诺平均规则,Stull.书 P41-44,1 状态方程,进行雷诺平均后：,2 连续方程,湍流脉动连续方程,湍流平均量连续方程,3 动量方程,3 动量方程,湍流应力或雷诺应力,重要！,P.S.：定常、水平均匀性、下沉,展开任一平均变量的全导数：,象位温、湍流动能e等平均变量，垂直变化大、水平变化很小；风速相反，u、v量级m/s，而w量级mm/s；因此方程中项多数情况下量级几乎相等。,水平平流 垂直,4、热量方程,：湍流热通量,5、水汽守恒,6、标量守恒,上面平均方程组均出现了湍流通量散度项，表现出湍流通量对平均场动量、热量和水汽含量增减的贡献。

11、,湍流中平均变量方程概要（略去分子扩散和粘性）,状态:,连续:,动量:,热量:,水汽:,标量:,【例 子 1】,设湍流热通量按 随高度线性递减，其中 a=0.3(K m s-1)和 b=310-4(K s-1)。如果初始位温廓线是任意形状（即选择某一形状），那么1小时后廓线的最终形状是什么样子？略去下沉、辐射、潜热加热，并假设水平均匀。,-【解法】-略去下沉、辐射、潜热加热，位温平均量方程为：,假设水平均匀，略去 x 和 y 导数，得到：,1 小时的增温是,【讨论】ML，所有高度空气以相同速率增温，廓线形状不变,【例 子 2】,如果10m/s的风速把干空气平流到某一区域，该区水汽水平梯度为(5

12、g水/kg气)/100km，那么要保持定常状态的比湿，边界层湍流水汽通量的垂直梯度多大？假设所有水都是气态，且不存在水汽体源。,问题中没提到下沉或水平通量梯度，假设为零，得到：,-【解法】-定常即，选择x坐标与平均风向一致，平均量方程:,【讨论】梯度大小相等于在垂直距离1km上 减少0.5(g/kg)(m/s),【作 业】,令 C 是空气中hockipuculis细菌浓度，大家知道这种传染病菌，每年冬天都要掠过美国北部一些州，当湖上结冰时它便增加。Sieve研究所发现了空气中hockipuculis的下列守恒方程:,式中 a 是常数，假设水平均匀、无下沉，求出湍流大气中的 守恒方程。,【个 例 分 析】,