POM模式中文手册.docVIP

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