FVCOM模型课题论文.docx

FVCOM模型简介及其应用实例摘要：通过对FVCOM水动力模型基本特点的分析，深入探究了在动力模型中常用的网格选取、方程组离散 求解和选取三维水流模式坐标等方式方法进一步分析了FVCOM在长江、珠江等河口地带对模拟潮汐、盐 度等海洋现象、要素的用途关键词：FVCOM模型有限体积法Q垂向坐标模型应用1 引言FVCOM( Finite Volume Coas tal Ocean Model)是美国 Massachuse tts Dartmouth 州立大 学陈长胜所领导的研究小组于2000年成功建立的海洋环流与生态模型模型包含动量方程、 连续方程、温盐守恒方程以及状态方程，数值模型采用有限体积法(FVM),优点为计算精 确快捷，并且可以较好地拟合海岸线边界和海底地形FVCOM由于其优越性，现在已经成为可以并行计算的模块化的可适用不同需求的模型 本文将着重介绍FVCOM模型的特点及方程，并列举几个应用实例，介绍FVCOM在近海研究中的 主要应用方向2 FVCOM模型特点FVCOM模型采用有限体积法对方程进行离散，综合了现有海洋研究中的有限差分和有限 元模型的优点模型在水平方向上采用无结构化非重叠的三角形网格，在垂向上使用b坐 标或者b -z混合坐标，使用干湿判别法处理潮滩移动边界。

FVCOM模型使用2.5阶湍流闭合 子模型来对控制方程进行封闭，并且分裂外模内模以节省计算时间下面将对这6个特点分 别进行介绍2.1无结构化非重叠的三角形网格 无结构化网格是指网格区域内的内部点不具有相同的毗邻单元即与网格剖分区域内的 不同内点相连的网格数目不同与结构化网格相比，无结构化网格可以方便的拟合复杂的边 界，也可以根据实际需要与进行局部加密，这个优点使其在研究岛屿众多，近岸岸线复杂的 问题时表现尤为突出图 1 FVCOM 三角网格设计[4]2.2离散方法因为计算机不能处理连续问题，所以有必要对海洋原始方程组进行离散求解，目前常用 的离散方法主要为以下三种2.2.1有限差分法有限差分法（FDM）是一种比较成熟和常用的离散方法，简单来说就是将求解域划分为 差分网格，用有限个网格节点代替连续的求解域，在微分方程中用差商代替偏导数，得到相 应的差分方程，通过解差分方程得到微分方程解的近似值2.2.2 有限元法有限元法（FEA）将定解区域看作很多不规则单元，每个单元的顶点被称为节点对每 一单元假定近似解（插值法），推导满足条件，得到包含有限个待定节点参量的简单场函数 通过能量方程然或者加权余量法建立离散方程，最终解出有限元法的数值解。

不考虑整个定 义域的复杂边界2.2.3 有限体积法有限体积法（FVM）将定解区域划分为不重叠的控制体积，每个网格点周围有一个控制 体积，将待解微分方程对每个控制体积积分得到离散方程，方程中的未知数为节点上因变量 的数值FVM得到的离散方程要求因变量在每个控制体积上都有积分守恒，这也就确保了在 整个区域上的积分守恒FDM的好处在于直观，理论成熟，易于编程和并行；FEM适合处理复杂区域但是FDM只 考虑网格点上的值，不考虑值在网格点间的变换；FEA内存和计算量巨大FVM结合了FDM和 FEM的优点，在数值计算中，既可以像FEA 一样使用灵活的几何结构，又具有FDM数值离散上 简单的特点FEA在复杂区域的适应性对FVM毫无优势可言，FVM还具有FEA不具有的物理量守 恒性，物理概念明显的这些优势2.3垂向坐标三维水流模式的垂向坐标一般分为等平面（z）坐标、等密度（p ）坐标和地形拟合（O） 坐标a） z坐标 （b）卩坐标 （0）坐标图2.三维水流模型垂向分层示意[5]2.3.1等平面（z）坐标Z坐标模式下的方程形式简单，便于数值离散海水速度及温度的水平变化 尺度远大于垂直变化尺度，z坐标在一定程度上符合变化规律。

水平压强梯度力 很容易求得，并不会引入截断误差但是在海底，等z面会和海底地形相交，在 网格离散时出现许多空洞，在浅水区域，由于采用绝对分层，很难满足必要的垂 直分辨率浅水区潮差比较大时，必须考虑自由面的变化，很容易造成能量和物 质的伪变化z坐标系下也很难反映沿倾斜密度面的物质运动过程2.3.2等密度面（p）坐标由于海洋内部的物质基本沿着等密度面的方向传输，因此等密度模型（即p模型）可以 很好地反映这一动力过程等密度坐标海洋模型对密度流有很高的数值精度,能较好地刻画 海洋中的锋面和层结p模型也很方便地计算水平压强梯度力但是p坐标在混合层和底部 边界层的分辨率较低（冬季浅海），因而在这些区域的表达较差p坐标的使用也会使得非 线性状态方程变得繁琐2 . 3.3地形拟合（o）坐标地形拟合（o）坐标最早出现在大气研究中，后来被引入海洋领域o坐标变换定义如 下：其中E为自静止海面向上起算的海面起伏（即潮位），H为海底到海平面距离，z为垂 向坐标，D为选定水柱的总长度O的变化范围为-1 （底部）到0 （表面）与传统的z坐标系相比，O坐标系统能很好地改变垂向分辨率欠缺的问题，但是并不能 很好的反应表面混合层。

这是因为，在深水区水平相邻网格点之间的垂向网格间距相差较大， 相比较z坐标而言，其分辨率并不高当然，这个缺陷可以通过变O分层（通常称为s坐标） 来改进，这样，即使远离海岸的海域，其表面混合层也可以保证一定的分辨率由于方程经过坐标变换，沿着倾斜密度面的对流项和扩散项变得更加繁琐在地形变化 剧烈的地方，计算水平压强梯度力时，精度往往很难保证因为等 面一般不是水平面，而 水平压强梯度力垂直于重力方向，这样水平压强梯度力在O等值面就会有一个投影在大陆 架海域，o面坡度达到1/100~1/10，在这样的坡度下，o坐标变换下的斜压梯度力会有较大 截断误差这也是o比较明显的缺陷2.4干湿判别法干湿判断法根据计算点的水深以及相邻点的水深和水位值判别该计算点的干湿，用以确 定计算区域由于水位变化产生的动边界单元水深＜0.005m的，为干单元，不参与计算；单 元水深〉O.lm的，为湿单元，参与计算；单元水深在0.005-0.1m之间的，为半湿单元，质 量通量参与计算，动量通量为 02.5 阶湍流闭合子模型在FVCOM的动量方程中，存在一个Km （垂向涡动粘性系数）项，此项在实际探测中很 难测定，就算测定不同地点也会有很大差异。

引入2.5阶湍流闭合子模型就是用来近似垂向 涡动粘性系数，使方程中的未知数个数等于方程个数，从而使方程闭合chi2.6 内外模分离 研究海水运动时，海水表面的重力波称为外重力波，也被称作潮波，对于外重力波的计 算无需考虑垂向坐标海平面以下的重力波称为内重力波，也就是海洋内波对于海洋内波 的计算需要考虑深度因素FVCOM外模式为二维，采用短的时间步长计算水位和垂直平均流速，水位直接提供给内 模态计算使用，垂直积分流速用来校正三维流速，外模态可从内模态获取底应力以及斜压项 和对流项的垂直积分内模式为三维，采用长的时间步长计算三维流速、温度、盐度和湍流 参数这样设置比完全三维计算节省很大计算量可进行2D，3D诊断和3D的斜压模式的计 算3 控制方程控制方程组由动量方程、连续方程、温度、盐度方程和密度方程组成动量方程Po创王u， v， w 是 x， y， z 方向其中在直角坐标系中 x， y， z 分别表示东，北和竖直坐标轴；uii dv的速度分量；T为温度;S为盐度；P为密度；P为压强；f为科氏参量；g为重力加速度;K为垂直旋转粘性系数；K为热量垂直旋转扩散系数；F , F , F ,和F代表水平动量， m h u v T S热量和盐度的扩散项。

根据实际情况，有时也可以用位温3代替温度T水平动量，热量和 盐度的扩散项为摩擦修正项还可以使用坐标变化得到O坐标下的控制方程，在此不再列举下面简要介绍水平动量方程及连续方程，温盐方程的构成形式与动量方程类似3.1 水平动量方程叟项为一局地微商，也称欧拉加速度，指固定空间点上的速度时间变化率,otu Ou + vOu + wOu项为平流变化的负值，表示由于沿流点运动方向速度分布不均导致的速 ox oy ozdu度变化，这两项的和为-，即拉格朗日加速度，表示流点的速度变化率fv项表示科氏1 OP力在x轴上的分量，- 表示x方向上的水平压强梯度力，后两项表示垂向粘性和水P O x0平摩擦对运动的影响dP3.2 连续方程dp (Ou Ou Ou) 一dt连续方程来源于公式击+P忖+石+疋J，以水团微元为研究对象时’默认在微Ou Ou Ou兀中密度保持不变，有+ = + = 0连续方程表示流体的散度为零，也就是流体是Ox Oy Oz不可压缩的4 FVCOM 模型应用实例FVCOM 水动力模型能够对近岸的一些基本海洋现象，如潮汐、环流、风暴潮、海浪以及 由此衍生的泥沙输运和海底地形变化进行模拟，并将其应用在海洋工程、环境保护和海洋资 源利用上，具有重要的作用。

4.1 FVCOM 在中国近海区域的应用4.1.1 FVCOM 模型在长江口杭州湾海域的一次台风过程计算曾经将此模型应用到了实际海区长江口—杭州湾附近的海域，考察了一次台风过程中浪、 流、泥沙之间的相互作用，并且讨论了耦合模型对于风暴潮数值模拟的作用此例选取了 1997年11号9711号台风“Winnie”，它于8月19日从浙江登陆为了充分考虑边界的影响并且节省计算时间，模型采用了嵌套的形式，包括一个大区网格和一个小区网格大网格 在长江口附近约为l-2km,小网格在长江口区域为200-300m大网格是为小网格海洋和海浪 模式提供开边界,没有用耦合型，具体来说分别用FVC0M和SWAVE计算了台风在小区开边界的 水位以及波浪在耦合模型中水平方向采用了球坐标垂向为坐标，均匀的分了 10层然 后分别对浓度场、浪场、泥沙浓度进行了模拟和分析浪、流、泥沙的相互影响,相互作用, 当其中的每一个要素被改变时，同时会影响到其他要素的计算那么在使用FVCOM模型时， 能够较为合理的模拟近岸的浪、流、潮汐、泥沙输运等现象同时存在的动力过程,是真正的 各要素同步计算并且从动力上耦合的模型2ft130.5122 122.5 I 牛 15.15 124 125l.unpiludc S"f:)图4 长江口网格【4】图5 东海网格【4】4.1.2 无结构网格耦合模型在渤海中应用与验证渤海是中国北部的一个接近封闭的海湾,只有东部通过渤海海峡与黄海相通,平均水深只 有18m。

模型水平方向采用球坐标,垂向为O坐标，均匀的分6层，计算时间共15天通过 FVC0M水动力模型、FVC0M波浪模型和耦合模型计算了风海流、潮流和波浪通过计算波浪 场和水位场的差值，浪流耦合对波浪和水位模拟影响最大的地方在渤海湾内，特别是靠近岸 边的地方，因此正确模拟浪流相互作用需要模型能够更详尽的分辨出近岸复杂的地形，这说 明在渤海的浪流耦合数值模拟中使用无结构网格是非常重要的40七m翳:;冷曲图6 渤海网格【4】120 (214.1.3 FVCOM 模型在象山港三维潮流盐度计算中的应用象山港是浙江省三大养殖基地之一,由于象山港的地形和水动力特点,其水体运动受到潮 汐、径流和地形的共同影响,湾内的流动具有很强的三维结构,任何一种二维模式（水平或垂 向）都不可能正确反映这个结构浙江省水利河口研究院的研究院利用了 FVCOM模型，建立 了象山港海域的三维有限体积潮流盐度模型水平方向上选用三角形非结构网格和垂直方向 上选用O。