第二章：基本方程.pptVIP

未知量：二阶矩量 三阶矩量 耗散 初看增加方程组中最后两个方程对我们妨害而不是有利，因为又增加了６个未知量。那我们为什么要这么做呢？理由是：认识ＴＫＥ方程和温度方差方程可以给我们提供湍流强度和湍流效果的数据，从而使我们在改进涡度扩散率参数化中能应用这些信息，参数化方案如下： * a．? 二重相关矩近似为平均值梯度的函数 此项的意义在于即使没有平均梯度也能输送热量。 b．? 三重相关矩近似为二重相关矩项梯度的函数 温度方差的垂直通量沿温度方差局部垂直梯度由高值向低值输送。 * 其中?是经验长度尺度，是在使模拟流场与实际实验或实际流场相匹配的尝试中根据实验数据和误差选择的。 Ｋ的表达式比较复杂，但我们可以用下式近似: 上述方程组太复杂，需进行数值求解。 这种闭合方法中，?的取值多少有些随意，增加了预报精度的难度。为此引进另一种闭合?e??闭合，除了ＴＫＥ方程外，还加进了高程度参数化的耗散率预报方程，因此避免了?误差。 ? * 改进的一阶半闭合（TKE）的边界层模式： * * (a)风速u分量廓线计算值(左边)和观测值(右边)的比较 u分量的模拟与实测符合较好,只有16:00时低层200m、300m范围内模式值偏小 模拟结果 * (b)风速v分量廓线计算值(左边)和观测值(右边)的比较 v分量较u分量差一些,14:00、18:00模拟值偏大约1~1.5m/s * (c)位温廓线计算值(左边)和观测值(右边)的比较 (c)中实测廓线是每隔3小时记录一次,模式值每2小时输出一次，12:00 的模拟结果与实测很一致； 14:00和18:00时计算值略低,但不超过1OC；18:00的逆温层的高度略低一些； 实测廓线中PBL顶处很薄的层次中位温发生很大变化,呈跃变状“而模拟结果中，PBL顶的跃变层不如实测结果明显”，原因是垂直格距较大的缘故。 * * E-?闭合：对TKE模式的发展，引入?预报方程 * 高阶闭合优点： 模拟出均匀混合层在白天随时间增厚的情况 夜间，稳定边界层发展，夜间急流形成的证据 白天湍流强度可增大到很大的数值，而夜间较小 低阶闭合无法给出湍流强度、脉动温度方差及平均梯度为0时的充分发展的混合层。 * 例题 （1）设湍流热通量随高度线性递减 其中a=0.3(Kms-1)和b=3×10-4(Ks-1)。如果初始位温廓线是任意形状，那么1小时后廓线的最终形状是什么样子？略去下沉、辐射和潜热加热并假设水平均匀。 这个结果不是z的函数，所以探测时每个高度上的空气以相同的速率增温，对时间t积分得1小时后的增温为1.08K。 * （2）如果10m/s水平风速把干空气平流到某一区域，该区水汽水平梯度为(5g水/kg气)/100km，那么要保持定常状态的比湿，边界层湍流水汽通量的垂直梯度是多大？假设所有的水都是汽态，而且不存在水汽体源。务必说明你的补充假设。 忽略下沉和水平通量梯度的数据有： * 第三节 湍流通量和方差预报方程 上面介绍的各个量如温度、速度等平均值预报方程中，均含有 等协方差项。这一节我们推导方差和协方差方程，方差方程向我们提供湍流能量和湍流强度的信息，而协方差方程则是描述运动湍流通量的。 一、方差方程 1．? 速度方差方程（运动方程） 以2ui?乘以速度脉动方程： * 利用乘积的微分规则做类似 的变换，并对整个方程求雷诺平均： 方程右端最后一项的雷诺平均为零； 方程左端最后一项：将湍流连续方程乘以 与这一项相加的通量形式 耗散项： * 则 I II III 第2项为速度方差的分子扩散，它含有方差的曲率，从混合层10-6s-2到近地层10-2s-2，乘以?后，量级为10-11?10-7m2s-3；第3项，如穿过一个直径为1cm的湍涡，其湍涡速度变化为0.1m/s，则穿过这个湍涡的瞬时切变为10s-1，对于更小的湍涡来讲，这个切变就更大了。该值平方乘以?后，其值为10-6?10-2m2s-3之间，混合层典型值10-4?10-3m2s-3。 定义粘滞耗散： * 气压扰动项：