海洋数值模型的理论和应用.ppt

第一节　海洋数值模型的发展概况 观测，试验 理论分析 数值模型 海洋的研究手段 海洋数值模型发展的客观要求 海水运动方程推导过程中的动力学假设，比如海水静力学假设，都是建立在大部分空间尺度大于10km，时间尺度大于惯性周期（2π/f）的情况下的 海水运动方程不存在精确解：（1）偏微分方程本身是非线性的，（2）将一个实际问题具体化时所需要的环境场如海底地形、海岸线几何形状、表面驱动力等都无法给出解析函数表达式，只能在离散化的空间和时间区域取得 因此，求解海水运动方程只能通过近似求解的办法来确定，而这些数值解又可能引入明显的近似误差 海洋数值模型发展的历史 海洋模型到按其水平网格的离散方式以及所使用的垂向坐标系的不同大致经历了如下几个发展阶段 最早出现并且还在使用的海洋模型是Bryan等人开发的基于原始方程的低阶精度的有限差分模型，它在水深方向采用z坐标系 目前常见的HOPS（ Harvard Ocean Prediction System ），MOM（ GFDL Modular Ocean Model），POP (Parallel Ocean Program)，NCOM（NCAR Community Ocean Model）模型在某种程度上都可以认为是该模型改进版 海洋数值模型发展的历史 上个世纪70年代，sigma坐标系开始应用于海洋模型在水深方向，比如目前被广泛使用的POM（Princeton Ocean Model）、ECOM （Estuarine Coastal and Ocean Model） 、ROMS（Regional Ocean Modeling System）模型都属于这种类型的模型 与传统的 z 坐标系相比，sigma坐标可以更好地贴合海底地形变化 海洋数值模型发展的历史 自上世纪90年代以来，非结构化网格技术开始在海洋模型中得到应用，相应的也出现了基于非结构化网格离散方式的有限元模型，如SEOM模型（ Spectral Finite Element Ocean Model），和有限体积模型，如FVCOM 虽然与传统的结构化网格相比,非结构化网格可以更好地拟合陆地边界，但是代码实现上的困难以及计算稳定性的问题使其迄今还没有得到非常广泛的应用 新一代的海洋数值模型 新一代的海洋模型广泛采用随地坐标系（ terrain-following coordinates），进而促进了有关时间步长，对流项和压力梯度项等数值算法的改进 进入新世纪以来，下一代的海洋数值动力模型正在紧锣密鼓的研制中，代表性的是TOMS (Terrain-following Ocean Modeling System)，它融合了目前最先进的物理知识、数值方法和数据同化技术 /WWWPUBLIC/htdocs.pom/TOMS.htm 海洋模型的商业软件 就海流的仿真建模研究而言，现阶段国外已形成了不少具有代表性的、具有强大前后处理功能及核心仿真建模技术的软件或系统（商业软件） 如加拿大Dalhousie University研究的DALCOAST河口海岸预报系统；丹麦的DHI Water Environment 机构开发的的MIKE系列软件系统；荷兰的WL Delft Hydraulics开发的的Delft－3D软件 第二节　基于POM模型的一般介绍 水动力模型基本特点 垂直方向的 𝝈 坐标变换 水平网格采用正交曲线坐标和ArakawaC差分格式 水平时间差分采用显格式，而垂直差分为隐格式 自由表面可以模拟水位变化 垂直和水平方向的混合扩散分别采用2.5阶的Mellor-Yamada湍流闭合模式和Smagorinski模式 内外模态分别处理速度较慢的内重力波和速度较快的外重力波以提高整个模式计算效率 包含了海水的热动力过程 Sigma坐标变换： 模型采用（a）垂向sigma坐标系和（b）水平Arakawa C网格 (1) 控制方程 连续方程 海水运动方程 (2) (3) (4) 温度方程 湍流动能方程 盐度方程 湍流混合长度方程 (5) (6) (7) (8) 其中，ω是基于sigma坐标系的垂向流速矢量，它与笛卡儿坐标系下的垂向流速W之间的关系可表示为 另外，方程中 分别表示海水的实际密度，Boussinesq 近似密度以及密度扰动 其中， 是实际水深 (9) M—Y湍流闭合模型（profq） 在湍流动能及混合长度方程中 分别表示湍流混合系数，热扩散系数和湍流动能的垂直扩散系数； 墙近似函数 ，其中 ，κ=0.4 是von Karman 常数 湍流动能和混合长度方程组由下列等式关系闭合 其中， 是稳定函