第3章大气边界层支配方程之1课件

单击以编辑母版标题样式,单击以编辑母版文本样式,第二级,第三级,第四级,第五级,*,第,三,章 大气边界层,支配,方程,3.1 基本控制方程,3.2 平均量方程,3.3 湍流脉动量方程,3.4 湍流方差预报方程,3.5 湍流通量预报方程,3.6 闭合理论,第三章 大气边界层支配方程 3.1 基本控制方程,1,3.1 基本控制方程,为了定量描述和预报边界层状况，需要借助于流体力学来描述大气中气体动力学和热力学方程。,描述气体和液体流动的方程组包括：,三个动量守恒方程（Navier-Stokes方程）,一个质量守恒方程（连续方程）,一个热力学能量方程,一个状态方程,(水汽及污染物浓度的标量方程类似热量方程),3.1 基本控制方程为了定量描述和预报边界层状况，需要借助于,2,3.1 基本控制方程,3.1.1 控制方程,3.1.2 简化与近似Boussinesq近似,3.1.3 坐标系及其变换,3.1 基本控制方程3.1.1 控制方程,3,3.1.1 控制方程,状态方程,质量守恒（连续方程）,动量守恒（牛顿第二定律）,热量守恒方程（热力学第一定律）,3.1.1 控制方程状态方程,4,1）状态方程,理想气体状态方程：,P：气压,：湿空气密度,R：干空气气体常数（R=287 J K,-1,kg,-1,）,虚温,T,v,=T(1+0.61q)，q比湿，T温度,1）状态方程理想气体状态方程：P：气压,5,2）质量守恒（连续方程）,连续性方程的一般形式：,不可压缩流体,或,泰勒假说:,2）质量守恒（连续方程）连续性方程的一般形式：不可压缩流,6,运用Einstein求和符号惯例，以上的连续方程写成：,（j=1、2、3分别代表x、y、z三个方向）,运用Einstein求和符号惯例，以上的连续方程写成：（j=,7,3）动量守恒（牛顿第二定律）,动量方程的表达式：,局地时间变化,平流项,气压梯度力项,重力作用项,科氏效应,粘滞应力项,3）动量守恒（牛顿第二定律）动量方程的表达式：局地时间变,8,i3,为克罗内克符号,重力仅在垂直方向起作用,思考：方程右端第三项展开？,i3为克罗内克符号 重力仅在垂直方向起作用思考：方程右端第,9,4）热量守恒（热力学第一定律）,热量守恒关系常用位温方程来表示：,s,：,引起空气位温变化的热量源(汇)项,（分子热传导、与相变有关的潜热释放、净辐射加热）,水汽及污染物的守恒方程形式与热量守恒形式一致,关键是要全面准确的了解引起成份变化的,源汇项,4）热量守恒（热力学第一定律）热量守恒关系常用位温方程来,10,3.1.2 简化与近似Boussinesq近似,在一定条件下，控制方程中某些项的量值比其它项小得多，以致可以把它略去，使得方程变得较为简单，有助于方程的求解。,必须考虑地球自转的影响,大气密度并非均匀，主要在垂直方向上是不均匀的层结流体,大气的空间尺度在水平方向远大于铅直方向，可视作浅层流体，不可压缩流体,大气边界层主要是湍流运动,大气边界层的特点概括有：,柯氏力作用,Boussinesq 近似,3.1.2 简化与近似Boussinesq近似,11,Boussinesq 近似的基本假定：,流体中的动力学粘滞系数,=,是常数,流体中的分子导热系数,k,T,是常数,大气是浅层流体，垂直范围约10km,描写流体热力状态的特征量可表示为,Boussinesq 近似的基本假定：流体中的动力学粘滞系数,12,扰动量远小于基态量：,基本状态量假设是静力平衡、绝热的，并且满足理想气体状态方程，即：,静力平衡,理想气体状态方程,绝热过程,扰动量远小于基态量：基本状态量假设是静力平衡、绝热的，并且满,13,Boussinesq 近似下的简化方程：,连续方程,状态方程,运动方程,热流量方程,Boussinesq 近似下的简化方程：连续方程,14,1）连续方程,当大气为浅层流体，在该范围内密度的变化很小、可以忽略，边界层大气可以近似认为是不可压缩的，连续方程为：,或,1）连续方程 当大气为浅层流体，在该范围内,15,2）状态方程,对于干空气，状态方程为,P=,RT,，对该式进行对数微分，可以得到：,如果运动的垂直尺度相对于大气层厚度而言很小，则,取对数：,取微分：,2）状态方程对于干空气，状态方程为 P=RT，对该式,16,3）运动方程,Navier-Stokes方程中压力梯度项和重力项可以进行改写,：,将上式代入,运动方程,得到：,扰动温度在重力作用下形成的净浮力项,3）运动方程Navier-Stokes方程中压力梯度项和,17,4）热流量方程,边界层支配方程中常用到,位温,，Poisson方程：,对该式进行对数微分，近似有,此外，有,因此，有,扰动温度的垂直梯度,位温的,垂直梯度,4）热流量方程边界层支配方程中常用到位温，Poisso,18,如果不考虑相变，热源,S,仅考虑分子热传导及辐射加热，则,热流量方程,为：,R,j,:j方向的辐射热通量,分子热传导的加热,辐射加热,上式中T,d,也可以换成,如果不考虑相变，热源 S 仅考虑分子热传导及辐射加热，则热,19,以上近似处理最早由Boussinesq(1903)提出,Boussinesq近似,该简化方程假定,流体不可压,、,并限制在一薄层内,。适用于研究像积云对流、海陆风环流、边界层急流中的重力波活动等发生在浅层内的中尺度运动,浅水方程,以上近似处理最早由Boussinesq(1903)提出Bo,20,综上所述，Boussinesq近似下的基本方程组为：,或,或,热流量方程,运动方程,状态方程,连续方程,综上所述，Boussinesq近似下的基本方程组为：或或热流,21,3.1.2 坐标系及其变换,通常我们使用笛卡尔坐标系，在实际应用中，将笛卡尔坐标系绕,z,轴旋转，使,x、y,轴指向其它方向，能够在处理问题时更方便。例如，使x轴与平均风向、地转风向、表面应力方向、垂直于海岸线、山的方向一致，这样就可以简化控制方程中的某些项。,例如，选择,x,轴与平均风方向一致，就可以求得,u,=M（平均风速）和,v,=0，这时，,x,轴叫做顺风方向，,y,轴叫做侧风方向。,3.1.2 坐标系及其变换 通常我们,22,3.2 平均量方程,以上动力方程组无法求解析解，但可以求得数值解。,原则上，可以用动力方程组来描述湍流的运动，但是要想囊括所有尺度的湍流运动，计算量太大。,为了简化，可以截取一定尺度的涡旋，而在这个尺度以下的涡旋用湍流的统计特征来代替。,在一些中尺度和天气尺度模式中，截取尺度为10100公里，而在边界层模式，比如说大涡模拟，截取尺度一般为10100米。,3.2 平均量方程以上动力方程组无法求解析解，但可以求得数值,23,一 出发方程组,1.状态方程,2.连续性方程,3.动量守恒,4.热量守恒,5.水汽守恒,6.标量守恒,一 出发方程组1.状态方程,24,1、状态方程,干空气气体常数,1、状态方程干空气气体常数,25,2、连续方程,张量展开：,2、连续方程张量展开：,26,3、运动方程(动量守恒),存储项,平流传输项,重力项，仅在垂直方向作用,柯氏力项,气压梯度力项,粘性力项,式中：,为分子动粘系数，,f,=2,sin,。,3、运动方程(动量守恒)存储项 平流传输项 重力项，,27,4、热量守恒,存储项,平流传输项,热扩散项,式中：,k,为分子热扩散系数，数值为,2.0610,-5,m,2,s,-1,；,L,p,为与,E,相变有关的潜热（0C时气液相变取值2.5010,6,Jkg,-1,；液固相变取值3.3410,5,Jkg,-1,;气固相变取值2.8310,6,Jkg,-1,）；,c,p,为湿空气定压比热，与干空气定压比热的关系为,c,p,=,c,pd,(1+0.86,q,)；,c,pd,取值1004.07 Jkg,-1,K,-1,；,E,为蒸发量,辐射加热项,与相变有关的潜热释放,4、热量守恒 存储项 平流传输项 热扩散项 式中：,28,5、水汽守恒,存储项,平流传输项,水汽扩散项,式中：v,q,为空气中水汽分子扩散率，,S,q,是方程中不含有其余过程时的净水体源项（源-汇），单位是：单位时间单位体积的总水体质量。,其它过程的净水体源项,5、水汽守恒 存储项 平流传输项 水汽扩散项 式中,29,6、标量守恒,存储项,平流传输项,C的扩散项,其它过程的体源项,式中：v,c,为空气中C的分子扩散率，,S,c,是不存在于方程中的其余过程的体源项，例如化学反应等。,6、标量守恒 存储项 平流传输项 C的扩散项,30,在湍流运动的大气边界层中，上述方程组还不能完整地描述边界层中的全部过程，应将上述的主要变量转换成平均量和脉动量相加。即：,平均场方程描述长时间过程，,脉动场方程描述短时间过程。,二 湍流中平均变量方程,在湍流运动的大气边界层中，上述方程组还不能,31,推导思路：,出发方程：Boussinesq 近似方程组,采用雷诺平均的方法，将任意一个物理量表示成,平均量和脉动量之和,，代入方程组，然后再,取平均,。,大气边界层平均量控制方程,推导思路：出发方程：Boussinesq 近似方程组,32,雷诺平均规则,Stull.书 P41-44,雷诺平均规则Stull.书 P41-44,33,1,状态方程,进行雷诺平均后：,最后一项很小、略去不计,平均量的状态方程,1 状态方程进行雷诺平均后：最后一项很小、略去不计平均量的,34,2 连续方程,湍流脉动连续方程,湍流平均量连续方程,2 连续方程湍流脉动连续方程湍流平均量连续方程,35,3 动量方程,3 动量方程,表示为雷诺应力对平均运动的影响,湍流应力或雷诺应力,重要！,再进行雷诺平均，得到：,3 动量方程3 动量方程表示为雷诺应力对平均运动的影响,36,假设，定常状态，即,假设，略去下沉，即,假设，水平均匀性，即,P.S.：定常、水平均匀性、下沉,展开任一平均变量,的全导数：,象位温,、湍流动能e等,平均变量，垂直变化大,、水平变化很小；,风速,相反，,u、v量级,m/s，而w量级mm/s；因此方程中,项,多数情况下量级几乎,相等,。,水平平流,垂直,假设，定常状态，即 假设，略去下沉，即假设，水平均匀性，即,37,4、热量方程,湍流热通量的输送对温度变化的影响,：湍流热通量,再进行雷诺平均，得到：,分子热传导引起的热通量，通常忽略,4、热量方程湍流热通量的输送对温度变化的影响：湍流热通量再进,38,5、水汽守恒,6、标量守恒,5、水汽守恒6、标量守恒,39,上面平均方程组均出现了湍流通量散度项，表现出湍流通量对平均场动量、热量和水汽含量增减的贡献。,上面平均方程组均出现了湍流通量散度项，表现出,40,湍流中平均变量方程概要,（略去分子扩散和粘性）,状态:,连续:,动量:,热量:,水汽:,标量:,湍流中平均变量方程概要（略去分子扩散和粘性）状态:连续:,41,【例 子 1】,设湍流热通量按 随高度线性递减，其中 a=0.3(K m s,-1,)和 b=310,-4,(K s,-1,)。如果初始位温廓线是任意形状（即选择某一形状），那么1小时后廓线的最终形状是什么样子？略去下沉、辐射、潜热加热，并假设水平均匀。,-【解法】-,略去下沉、辐射、潜热加热，位温平均量方程为：,假设水平均匀，略去,x,和,y,导数，得到：,1 小时的增温是,【讨论】,ML，所有高度空气以相同速率增温，廓线形状不变,【例 子 1】设湍流热通量按,42,【例 子 2】,如果10m/s的风速把干空气平流到某一区域，该区水汽水平梯度为(5g水/kg气)/100km，那么要保持定常状态的比湿，边界层湍流水汽通量的垂直梯度多大？假设所有水都是气态，且不存在水汽体源。,问题中没提到下沉或水平通量梯度，假设为零，得到：,-【解法】-,定常即 ，选择,x,坐标与平均风向一致，平均量方程:,【讨论】,梯度大小相等于在垂直距离1km上 减少0.5(g/kg)(m/s),【例 子 2】如果10m/s的风速把干空气平流到,43,【作 业】,令 C 是空气中hockipuculis