POM模式中文手册及讲解.pdf

POM 模式中文手册 POM(Princeton Ocean Model)是由美国普林斯顿大学于 1977 年共同建立起来的一个三 维斜压原始方程数值海洋模式，后经过多次修改成为今天的样本，是被当今国内外应用较为 广泛的河口、近岸海洋模式。POM 对于中国海域的数值模拟研究有很多文献。POM 在国内 较多人使用，在天津、上海、厦门等多个沿海地区均有人使用POM 模式进行风暴潮的模拟 和预报。 POM 采用蛙跳有限差分格式和分裂算子技术，水平和时间差分格式为显式，垂向差分 格式为隐式，对慢过程(平流项等)和快过程(产生外重力波项)分开，分别用不同的时间步长 积分，快过程的时间步长受严格的 CFL 判据的限制。内外模分离技术比完全三维计算节省 很大计算量。为消除蛙跳格式产生的计算解，POM 在每一时间积分层次上采用了时间滤波。 水平方向采用正交曲线网格，变量空间配置使用“Arakawa C ”网格，可以较好的匹配岸界。 与均匀网格相比，水平曲线正交网格是渐变的，能更好地拟合岸线侧边界，减少“锯齿”效应。 POM 模式在垂向上采用了σ 坐标变换，可体现不规则的海底地形的变化特点，便于引入大 陆架地形 并且引入了干湿网格动边界技术，既可更好地处理三维水动力环境模拟中大量浅 滩的“干出”与“淹没”等难点问题，也可很好地处理复杂地形水域的模拟问题，因此被广泛地 应用于河口近岸海域的潮流数值模拟中。基于POM 模式源程序代码的公开性，便于学者交 流与学习，并可根据实际工作问题的需要进行改进，应用到不同的领域，因而具有很强的生 命力和适用性。 第一章 二维浅水方程数值模式（潮汐、风暴潮数值计算模式） 20 世纪60 年代初期，由于电子计算机计算能力的提高，人们开始了海洋数值预报的 研究和试验工作。1960 年代中期，初步实现了海洋潮汐问题的二维数值模拟，1980 年代中 期二维潮汐数值模拟技术已经成熟。1970 年代中期开始三维海洋温、盐、流数值模拟方法 研究，在1990 年代开始试验性的海洋数值预报（包括温、盐、流和水位）。海浪的数值预 报研究始于 1950 年代后期，1960 年代中期开始进行数值计算实验。 本章讨论二维浅水方程组（也称作深度平均的二维长波运动方程）的数值计算问题，这 种二维数值模式可用于潮汐、风暴潮的模拟和数值预报；加入二维扩散方程后，也可用于浅 海污染物（溶解质）扩散的计算；加入泥沙（悬移质和推移质）运动方程，可用于浅海泥沙 运动的数值模拟。 一、基本方程组 1 在浅海中，由于水平粘性项的重要性不如底摩擦，因此有时忽略此项，或取一个比较任 意的A 值；但是这一项对数值计算的稳定性起到有利的作用。 在边缘海，由于潮汐主要是大洋潮波传入引起的，引潮力的作用相对很小，因此在做边 缘海的数值模拟时一般不考虑引潮力的作用。 二、初始条件与边界条件 2 三、差分方程格式（ADI ） 1、网格设置 空间差分网格采用C 网格，即ζ 、h 位于网格线交叉点上；U 位于x 方向网格线上， 并在相邻的ζ 之间；V 位于y 方向网格线上，并在相邻的ζ 之间。两个同类型的相邻点之 间的距离为ΔS 。 对于具体的海区，则在(x, y)平面上的计算点将具有不同的类型，分为内点和边界点； 边界点进一步可分为闭边界点和开边界点；而开边界点又可进一步分成ζ 点、U 点和V 点。 关于边界点的选取这里还需要特别加以说明。对于闭边界，我们可以相信，海岸线的微小变 动对于整个海区的水体运动的影响将是微小的，只是在离岸很近的区域，才会有较严重的影 响。因此可以在原来的岸线附近假设一条闭边界线，它在U ，V 点通过，而且当它通过U 点时，方向与x 轴方向垂直，当它通过V 点时，方向与y 轴方向垂直。这样，在这些点上 的边界条件就可简单地取作： 一般情况下，只需