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

湍流运动的两种研究方法湍流运动的两种研究方法1、采用随机过程的统计学方法来反映、采用随机过程的统计学方法来反映大气湍流结构大气湍流结构 平均、方差、标准差、湍强、相关、通平均、方差、标准差、湍强、相关、通 量、应力、湍流动能量、应力、湍流动能 2、解湍流运动控制方程、解湍流运动控制方程（平均运动方（平均运动方程、脉动方程、湍能方程程、脉动方程、湍能方程.）第第三三章章 大气边界层大气边界层支配支配方程方程 3.1 3.1 基本控制方程基本控制方程 3.2 3.2 平均量方程平均量方程 3.3 3.3 湍流脉动量方程湍流脉动量方程 3.4 3.4 湍流方差预报方程湍流方差预报方程 3.5 3.5 湍流通量预报方程湍流通量预报方程 3.6 3.6 闭合理论闭合理论3.1 3.1 基本控制方程基本控制方程为了定量描述和预报边界层状况，需要借助于流为了定量描述和预报边界层状况，需要借助于流体力学来描述大气中气体动力学和热力学方程。体力学来描述大气中气体动力学和热力学方程。描述气体和液体流动的方程组包括：描述气体和液体流动的方程组包括：三个动量守恒方程（三个动量守恒方程（Navier-StokesNavier-Stokes方程）方程）一个质量守恒方程（连续方程）一个质量守恒方程（连续方程）一个热力学能量方程一个热力学能量方程一个状态方程一个状态方程 (水汽及污染物浓度的标量方程类似热量方程水汽及污染物浓度的标量方程类似热量方程)3.1 3.1 基本控制方程基本控制方程3.1.1 控制方程控制方程3.1.2 简化与近似简化与近似Boussinesq近似近似3.1.3 坐标系及其变换坐标系及其变换3.1.1 控制方程控制方程 状态方程状态方程 质量守恒（连续方程）质量守恒（连续方程）动量守恒（牛顿第二定律）动量守恒（牛顿第二定律）热量守恒方程（热力学第一定律）热量守恒方程（热力学第一定律）1）状态方程状态方程vRTP理想气体状态方程：理想气体状态方程：P：气压：气压：湿空气密度：湿空气密度R：干空气气体常数（：干空气气体常数（R=287 J K-1 kg-1）虚温虚温Tv=T(1+0.61q)，q比湿，比湿，T温度温度2）质量守恒（连续方程）质量守恒（连续方程）连续性方程的一般形式：连续性方程的一般形式：0zwyvxut0t不可压缩流体不可压缩流体0zwyvxuzwyvxudtd或或0zwyvxutdtd泰勒假说泰勒假说:运用运用Einstein求和符号惯例，以上的连续方程写成：求和符号惯例，以上的连续方程写成：0zwyvxut0zwyvxuzwyvxudtd0jjxut0iixujjxudtd（j=1、2、3分别代表分别代表x、y、z三个方向）三个方向）3）动量守恒（牛顿第二定律）动量守恒（牛顿第二定律）动量方程的表达式：动量方程的表达式：jji2jij3i3ijijixxuufgxpxuutu1局地时间变化局地时间变化平流项平流项气压梯度力项气压梯度力项重力作用项重力作用项科氏效应科氏效应粘滞应力项粘滞应力项jji2jij3i3ijijixxuufgxpxuutu1等，如任意两个相等，如逆序排列，如顺序排列333322113013221332113122311231ijkijkijkijkjijiij当当01i3为克罗内克符号为克罗内克符号 重力仅在垂直方向起作用重力仅在垂直方向起作用0,3fufvufjij思考：方程右端第三项展开？思考：方程右端第三项展开？4）热量守恒（热力学第一定律）热量守恒（热力学第一定律）热量守恒关系常用位温方程来表示：热量守恒关系常用位温方程来表示：szwyvxuts：引起空气位温变化的热量源引起空气位温变化的热量源(汇汇)项项（分子热传导、与相变有关的潜热释放、净辐射加热）（分子热传导、与相变有关的潜热释放、净辐射加热）水汽及污染物的守恒方程形式与热量守恒形式一致水汽及污染物的守恒方程形式与热量守恒形式一致 关键是要全面准确的了解引起成份变化的关键是要全面准确的了解引起成份变化的源汇项源汇项3.1.2 简化与近似简化与近似Boussinesq近似近似 在一定条件下，控制方程中某些项的量值比其它项小在一定条件下，控制方程中某些项的量值比其它项小得多，以致可以把它略去，使得方程变得较为简单，有助得多，以致可以把它略去，使得方程变得较为简单，有助于方程的求解。于方程的求解。必须考虑地球自转的影响必须考虑地球自转的影响大气密度并非均匀，主要在垂直方向上是大气密度并非均匀，主要在垂直方向上是不均匀的层结流体不均匀的层结流体大气的空间尺度在水平方向远大于铅直方大气的空间尺度在水平方向远大于铅直方向，可视作浅层流体，不可压缩流体向，可视作浅层流体，不可压缩流体大气边界层主要是湍流运动大气边界层主要是湍流运动大气边界层的特点概括有：大气边界层的特点概括有：柯氏力作用柯氏力作用Boussinesq 近似近似Boussinesq 近似的基本假定：近似的基本假定：流体中的动力学粘滞系数流体中的动力学粘滞系数=是常数是常数流体中的分子导热系数流体中的分子导热系数 kT 是常数是常数大气是浅层流体，垂直范围约大气是浅层流体，垂直范围约10km描写流体热力状态的特征量可表示为描写流体热力状态的特征量可表示为),()(0tzyxPzPPd),()(0tzyxzd),()(0tzyxTzTTd000:TP基本状态dddTP:扰动量扰动量远小于基态量：扰动量远小于基态量：|0PPd|0d|0TTd基本状态量假设是静力平衡、绝热的，并且满足基本状态量假设是静力平衡、绝热的，并且满足理想气体状态方程，即：理想气体状态方程，即：静力平衡静力平衡理想气体状态方程理想气体状态方程绝热过程绝热过程gzP00000RTPdzT0Boussinesq 近似下的简化方程：近似下的简化方程：连续方程连续方程 状态方程状态方程 运动方程运动方程 热流量方程热流量方程1）连续方程连续方程 当大气为浅层流体，在该范围内密度的变化当大气为浅层流体，在该范围内密度的变化很小、可以忽略，边界层大气可以近似认为是不很小、可以忽略，边界层大气可以近似认为是不可压缩的，连续方程为：可压缩的，连续方程为：0zwyvxu0iixu或或0zwyvxut2）状态方程状态方程对于干空气，状态方程为对于干空气，状态方程为 P=RT，对该式进行，对该式进行对数微分，可以得到：对数微分，可以得到：TRRTPlnlnlnlnlndTTddPP111如果运动的垂直尺度相对于大气层厚度而言很小，则如果运动的垂直尺度相对于大气层厚度而言很小，则 00PPd取对数：取对数：取微分：取微分：000000TTTPPP00TTdd3）运动方程运动方程Navier-Stokes方程中压力梯度项和重力项可以进方程中压力梯度项和重力项可以进行改写行改写：30030030330300330003111111ididididididdidiidiididiididiigTTxPgxPgxPggxPxPggxPxPgxPjji2jij3i3ijijixxuufgxpxuutu1jjijijididjijixxuufgTTxPxuutu233001将上式代入将上式代入运动方程运动方程得到：得到：扰动温度在重力作用下形成的净浮力项扰动温度在重力作用下形成的净浮力项 0g0dTT4）热流量方程热流量方程边界层支配方程中常用到边界层支配方程中常用到位温位温，Poisson方程：方程：CpRhPPT00对该式进行对数微分，近似有对该式进行对数微分，近似有0hzTd 000000hThzhTTTTTdd此外，有此外，有zzTd因此，有因此，有扰动温度的垂扰动温度的垂直梯度直梯度位温的位温的垂直梯度垂直梯度如果不考虑相变，热源如果不考虑相变，热源 S S 仅考虑分子热传导及辐仅考虑分子热传导及辐射加热，则射加热，则热流量方程热流量方程为：为：jjpjjdTjdjdxRcxxTkxTutT021Rj:j:j方向的方向的辐射热通量辐射热通量分子热传导的加热分子热传导的加热辐射加热辐射加热上式中Td 也可以换成 以上近似处理最早由以上近似处理最早由Boussinesq(1903)Boussinesq(1903)提出提出BoussinesqBoussinesq近似近似 该简化方程假定该简化方程假定流体不可压流体不可压、并限制并限制在一薄层内在一薄层内。适用于研究像积云对流、。适用于研究像积云对流、海陆风环流、边界层急流中的重力波海陆风环流、边界层急流中的重力波活动等发生在浅层内的中尺度运动活动等发生在浅层内的中尺度运动浅水方程浅水方程综上所述，综上所述，Boussinesq近似下的基本方程组为：近似下的基本方程组为：jjpjjdTjdjdxRcxxTkxTutT021zzTdjjijijididjijixxuufgTTxPxuutu2330010zwyvxu0iixu或或00TTddTTPP或或热流量方程热流量方程运动方程运动方程状态方程状态方程连续方程连续方程3.1.2 坐标系及其变换坐标系及其变换 通常我们使用笛卡尔坐标系，在实际应用中，将笛卡通常我们使用笛卡尔坐标系，在实际应用中，将笛卡尔坐标系绕尔坐标系绕z轴旋转，使轴旋转，使x、y轴指向其它方向，能够在处轴指向其它方向，能够在处理问题时更方便。例如，使理问题时更方便。例如，使x轴与平均风向、地转风向、轴与平均风向、地转风向、表面应力方向、垂直于海岸线、山的方向一致，这样就可表面应力方向、垂直于海岸线、山的方向一致，这样就可以简化控制方程中的某些项。以简化控制方程中的某些项。例如，选择例如，选择x轴与平均风方向一致，就可以求得轴与平均风方向一致，就可以求得u=M（平均风速）和（平均风速）和v=0，这时，这时，x轴叫做顺风方向，轴叫做顺风方向，y轴叫做轴叫做侧风方向。侧风方向。3.2 3.2 平均量方程平均量方程以上动力方程组无法求解析解，但可以求得数值解。以上动力方程组无法求解析解，但可以求得数值解。原则上，可以用动力方程组来描述湍流的运动，但是原则上，可以用动力方程组来描述湍流的运动，但是要想囊括所有尺度的湍流运动，计算量太大。要想囊括所有尺度的湍流运动，计算量太大。为了简化，可以截取一定尺度的涡旋，而在这个尺度为了简化，可以截取一定尺度的涡旋，而在这个尺度以下的涡旋用湍流的统计特征来代替。以下的涡旋用湍流的统计特征来代替。在一些中尺度和天气尺度模式中，截取尺度为在一些中尺度和天气尺度模式中，截取尺度为10100公里，而在边界层模式，比如说大涡模拟，截取公里，而在边界层模式，比如说大涡模拟，截取尺度一般为尺度一般为10100米。米。一一 出发方程组出发方程组1.状态方程状态方程2.连续性方程连续性方程3.动量守恒动量守恒4.热量守恒热量守恒5.水汽守恒水汽守恒6.标量守恒标量守恒1 1、状态方程、状态方程vdTRp 干空气气体常数干空气气体常数-1-1kgKJ 04.287 dR2 2、连续方程、连续方程wuvuuuixuii 321 3,2,1 0张量展开：张量展开：0 0z zw wy yv vx xu u 3 3、运动方程、运动方程(动量守恒动量守恒)存储项存储项 平流传输项平流传输项 重力项，仅重力项，仅在垂直方向作用在垂直方向作用 柯氏力项柯氏力项 气压梯度力项气压梯度力项 粘性力项粘性力项 式中：式中：为分子动粘系数，为分子动粘系数，f=2 sin。jjiidjijidjijixxuxPufgTTxuutu2033014 4、热量守恒、热量守恒 存储项存储项 平流传输项平流传输项 热扩散项热扩散项 式中：式中：k为分子热扩散系数，数值为为分子热扩散系数，数值为2.0610-5 m2s-1；Lp为与为与E相相变有关的潜热（变有关的潜热（0C时气液相变取值时气液相变取值2.50106 Jkg-1；液固相变取；液固相变取值值3.34105Jkg-1;气固相变取值气固相变取值2.83106 Jkg-1）；）；cp为湿空气定为湿空气定压比热，与干空气定压比热的关系