森林下垫面热动力过程及其物理机制的大涡模拟研讨.pdfVIP

森林下垫面热动力过程及其物理机制的大涡模拟研究* 苗世光1,2 蒋维楣2 1中国气象局北京城市气象研究所 北京 100089 2南京大学大气科学系 南京 210093 摘 要 在采用各向异性湍流动能闭合方案和3阶Runge-Kutta时间积分方案的大涡模式中 引入由森林冠 层粗糙元造成的动量拖曳项 热量输入项和TKE耗散项 以模拟森林冠层和森林边界层的气象场 通过中 性和不稳定层结条件下不同叶面积指数算例的模拟分析及与已有观测结果的比较表明 本文所建大涡模式 对森林冠层和森林边界层有较好的模拟效果 进一步研究表明 不稳定层结条件下较稠密的森林冠层中特 有的KinkingPairing湍涡结构与森林边界层中湍流的大涡运动相互作用 形成了森林冠层附近的温度斜 坡型结构 关键词 森林冠层 森林边界层 大涡模拟 温度斜坡型结构 KinkingPairing湍涡结构 1 引言 大气与陆面之间的动量 热量和物质交换 水汽 二氧化碳等 对诸多的气象 环境问 题 如微气象 植被层通量 尘埃的飞扬等都是至关重要的 在大气边界层的研究中 有时 将植被简化成单层下垫面 对于研究农林牧有关的微气象问题 显然是过于粗糙的 并且这 样做往往不能正确预测陆气界面的通量 从而影响模式的可靠性 因此 对植被层湍流的研 究是十分必要的 在数值模拟的三种方法中 直接模拟(DNS)方法直接求解适合所研究问题初始条件和边 界条件的Navier -Stokes NS 方程组 不对其作任何重要简化和假设 计算量较大 由于大 气运动是高Reynolds 数运动 用直接模拟方法研究时 计算量巨大 难以实现 系综平均 是应用最早的一种模拟方法 在 研究上得到广泛的应用 但是 寻求系综平均 (RANS) PBL 的一种普遍适用的闭合方案是很困难的 更重要的是计算结果丢失了瞬时的信息 不能给出 流场的瞬时结构 所以无从解释因大涡结构而导致的下扫 喷发和温度斜坡型特征的机理 因而也无法改进通量计算的精度 大涡模拟(LES)从一定程度上可以避开这些困难 LES技 术是直接模拟的原理和计算量的限制在某种程度上的一种折中方法 已成为研究湍流问题的 主要手段之一 [1] [2] [3] Baldocchi 等 和Gao 等 对森林冠层的湍流结构进行了系统的观测研究 Shaw 等 首 [4] 次用LES 方法研究了森林冠层的湍流流场 Su 等 分析了中性层结稀疏森林冠层大涡模拟 和外场观测的湍流通计量 然而 对森林全边界层的研究尚未见诸于国内外文献 本文作了 这方面的初步尝试 以期为非均匀下垫面边界层湍流结构和流场特征研究及污染物扩散研究 提供研究思路和高效能的技术手段 文中采用各向异性湍流动能闭合方案和3 阶Runge-Kutta 时间积分方案的大涡模式[5] 引入由森林冠层粗糙元造成的动量拖曳项 热量输入项和TKE 耗散项 以模拟森林冠层和 森林边界层的气象场 通过中性和不稳定层结条件下不同叶面积指数算例的模拟 分析森林 冠层和森林边界层的湍流统计特征 以及湍流流场的空间结构 并与已有的外场观测 同类 数值模拟结果和理论研究作比较 2 大涡模拟中森林冠层的描述 * 本项研究得到国家自然科学基金项共同资助 sgmiao@ium.cn; jwm@nju.edu.cn 1 [5] 苗世光等 所建大涡模式的方程组由连续方程 动量方程 热流量方程组成 采用先进 的各向异性TKE 闭合方案和精度较高的3 阶Runge-Kutta 时间积分方案 模式在水平方向 和垂直方向都采用等距网格 水平方向采用伪谱方法 垂直方向采用交错网格和二阶有限差 分方法 模式的四个侧边界取周期性的边界条件 上边界采用辐射边界条件,下边界采用M-O