omega调报告.docVIP

omega调研报告 1.前言 Omega是可调整网格多尺度环境模式的简称，是一个新的大气模拟系统，用运了全面的非流体静力学方程组这一先进的计算流体力学方法。它基于非结构的三棱柱网格，该网格的水平网格尺度为100km到1km，垂直方向上尺度从边界层中的几米到自由大气中的1km。因此，它能够模拟全尺度的危险物释放。它代表着天气预报中一个非常先进的网格处理方法，它能够迅速地适应包含有地面和地形特征网格；它代表了一种很先进的大气释放的紧急响应；它代表了一种很先进的浮质输送和扩散的模式。它较多的运用于大气释放、台风预报以及点的天气预报等方面。 2.正文 2.1 OMEGA网格的产生 OMEGA基于一个非结构的三棱柱网格。这种网格在水平方向上是非结构的但在垂直方向上时结构化的。使用的是旋转笛卡尔坐标系（图1）。这种网格产生的过程如下： 第一步：初始化网格时，使用一个二十面体(图2)，该二十面体有十二个节点三十条边，让每一个节点在球面上。将十二个节点按照从北向南的顺序标号，第一个点在北极，第二个点在x轴上，依次标号，则这十二个点的坐标（表1）可给出。 图1 OMEGA的笛卡尔坐标系及垂直格点系列 图2 左面为二十面体的正面 右面为二十面体的反面 表1 节点的坐标（设地球为正球体，半径为1，为正二十面体） Node i X（i） Y（i） Z（i） 1 0 0 0.5*sqrt(5.0) 2,3,4,5,6 cos{(i-2)*2*π/5} sin{(i-2)*2*π/5} 0.5 7,8,9,10,11 cos{π/5 +(i-7)*2*π/5} sin{π/5 +(i-7)*2*π/5} -0.5 12 0 0 -0.5*sqrt(5.0) 第二步：以每个三角形边的中点为节点，将每个面分为四个三角形（图3），再将这三个节点投影到球面上。经过六次叠代，可达81920个网格（20*4n）。网格面积约为15.5*108m2，最大周长为131km，最小周长为110km，边长在30km到40km。 图3 叠代产生网格 第三步：调整 调整是在第二部之后根据一定条件让格点合理化的过程，包括四个步骤： 1） 网格三等分 在三角形中心增加一个节点，则会有两个额外的三角形产生。三等分会产生钝角三角形，所以要经过修正（边缘调整）。 2） 边对分 主要用在不协调的边界上，节点的插入不严格，可根据长宽率调整。 3） 节点松弛 通过周围节点将某个节点移动到多边形中心，是个简单的平滑过程。 4）重接（边缘交换） 在删除顶点后，有一个紧接着的重接过程。在三等分后有一个紧接着的重接过程。用于消除钝角三角形，三角形角度可以计算出来。重接步骤包括计算每个单元看是否有低的长宽比和移动边及在连接顶点。 OMEGA网格加密和放大条件： 加密条件：①水陆边界 ②坡度陡峭处 ③相邻单元太小 放大条件：①长宽比不合理 ②平坦地区 ③陆地网格周围是陆地或水面网格周围是水面（性质相同） OMEGA在垂直方向上的划分： 在水平的非结构网格产生后，垂直方向的分层就会从地球中心开始向外分层。垂直层的高度可以是常数或可变高度。在边界层，它比较合理的接近地面。垂直向的高度通过指定的第一层的高度和伸长比率决定，以确保在近地面比较合理。垂直网格在地面时沿着地形的，但在计算区域顶层上是固定的（图4）。非结构网格在垂直方向上增加了当地的确定性，所以能够更好地捕获地形或者重要的大气环流和云动力学的物理量特征。三棱柱网格点能很精确的代表地形信息，与MM5传统网格相比较有三大优点。第一个优点为OMEGA没有要求像传统中尺度模式MM5样的2-4的地形平滑；第二个优点为OMEGA静力适应能够很好的解决在复杂地形及水陆混合边界（complex land/water boundaries）的物理过程；第二个优点为动力适应增加了当需要时（例如区域中的天气系统或超过区域的烟羽的路径）的空间物理过程的处理能力，它能自动调整网格。 图4 OMEGA在垂直方向上的划分 图5 一个半球面的omega网格 2.2 OMEGA使用的坐标系及插值方法 OMEGA模式用的是旋转笛卡尔坐标系。在旋转笛卡尔坐标系中，xoy平面在赤道平面上，让x轴与本初子午线垂直，z轴垂直于xoy平面穿过北极点，y轴和x轴和z轴垂直（图1）。每组格点元素为三棱柱网格中心点的值（图6），这些元素在垂直方向上按统一方式堆叠着。在所用单元产生后，所有物理量用P-L平均权重方案插值。 图6 OMEGA格点元素 一个顶点上P-L平均权重方案插值方法公式如下： （1） （2） 需是个稳定的权重，要尽量一致 由（1）式和（2）式可得： （3） （4） 定义为功能函数，用最小功能函数的最优化问题给出约束方程（3）和（4） 用拉