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

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