基于FVCOM和OpenFOAM模型的多尺度水动力耦合方法与流程

本发明涉及潮流与海洋构筑物相互作用多尺度计算技术领域，具体为一种基于fvcom和openfoam模型的多尺度水动力耦合方法。

背景技术：

近年来，海洋灾害频发，沿海地区生命财产损失严重，而数值模拟技术作为一种评估预测工具可有效地对海洋灾害发生过程及结果进行模拟计算。但海洋潮波运动属于大/中尺度或中尺度运动，作用于海洋构筑物时变为小尺度或微尺度运动，例如海啸的运动尺度约1000km，而作用于构筑物后（如桥梁、房屋等）变为0.1~10m的尺度。目前，单一尺度的数值模拟技术，无法详细解析构筑物局部大尺度海洋动力引起的水体运动，并描述构筑物表面水体飞溅、冲击等物理现象。如：单一使用fvcom模型数值模拟技术，可对大/中尺度或中尺度海流运动进行模拟计算，并具有较好的效果，但无法准确模拟海流作用于海洋构筑物时所产生的小尺度或微尺度涡流运动；单一使用openfoam模型数值模拟技术可对海洋构筑物局部涡流运动进行模拟，并具有较高的计算精度，但无法快速模拟大/中尺度或中尺度海流运动。

技术实现要素：

本发明目的在于提供一种基于fvcom和openfoam模型的多尺度水动力耦合方法，以实现大/中尺度与小/微尺度的水域潮流数据的耦合计算，进而评估预测海洋灾害过程。

本发明所述的基于fvcom和openfoam模型的多尺度水动力耦合方法，包括以下步骤：

s1：依据海洋构筑物所处的远端大/中尺度水域作为外域建立fvcom模型，海洋构筑物所处的近区小/微尺度水域作为内域建立openfoam模型；

s2：fvcom模型依托粗分辨率二维非结构网格，模拟三维分层潮流流动，openfoam模型依托小分辨率精细网格，模拟局部三维水体流动；

s3：基于水位初始场和外域潮位边界瞬时值，运用fvcom模型计算获得外域三维分层潮流场，结合三线插值法得到每个网格质心处的潮流值；

s4：利用最邻近点插值方法，首先将潮流数据从非结构网格质心处传递到融合区中间结构化网格节点处，随后进一步将其传递到内域模型openfoam边界节点上；

s5：以到达内域模型边界处的潮流数据作为初始值，通过计算得到潮流与海洋构筑物相互作用结果，解析构筑物近区大尺度海洋动力引起的水体飞溅和冲击现象。

一种基于fvcom和openfoam模型的多尺度水动力耦合方法，通过利用水位初始场和外域潮位边界瞬时值，采用fvcom模型计算获得外域三维分层潮流场，运用三线插值法得到每个网格质心处的潮值。在此基础上，使用最邻近点插值方法，首先将三维分层潮流数据从非结构网格质心处传递到融合区中间结构化网格节点处，随后进一步将其传递到内域模型openfoam边界上作为边界值。以到达内域模型边界处的潮流数据作为初始值，通过计算得到潮流与海洋构筑物相互作用结果，建立一套行之有效的新的耦合方法，可解决大尺度潮流数据与小尺度海洋构筑物相互作用的流固耦合模拟问题，从而捕获复杂的物理现象，如，构筑物表面产生的水体飞溅和迎水面水体爬高。并兼顾模型网格间距和计算效率，以及允许解决方案在远场和近场之间的无缝转换。因此该方法能够同时并直接模拟许多之前无法处理的多尺度流固耦合问题，提高了不同尺度海流运动的高仿真模拟能力。

综上，为了克服原始navier-stokes方程求解过程中的黏性流体计算量大、速度较慢的缺点，本发明结合求解大范围三维分层潮流场的fvcom开源模型和求解局部精细尺度全三维潮流动力的openfoam开源模型建立了多尺度水动力耦合方法。在耦合模型的外域，通过fvcom数值模型计算潮流，获得水位、流速等潮流动力变化。在耦合模型的内域，采用openfoam模拟精细尺度流体动力变化。耦合模型相比于单纯的openfoam模型，可以缩小计算区域，从而减小计算量；相较于单纯的fvcom模型，可以获得海洋构筑物近区全三维精细尺度动力特征，提高模型精度和应用范围。

附图说明

图1为本发明fvcom与openfoam模型多尺度水动力耦合方法流程图；

图2为本发明模拟fvcom和openfoam模型多尺度潮流传播耦合示意图；

图3本发明fvcom与openfoam模型网格分布示意图；

图4本发明fvcom与openfoam模型耦合数据传递示意图。

具体实施方式

如图1所示，本发明所述的基于fvcom和openfoam模型的多尺度水动力耦合方法，包括以下步骤：

依据海洋构筑物所处的远端大/中尺度水域作为外域建立fvcom模型，海洋构筑物所处的近区小/微尺度水域作为内域建立openfoam模型；fvcom模型依托粗分辨率二维非结构网格，模拟三维分层潮流流动，openfoam模型依托小分辨率精细网格，模拟局部三维水体流动；基于水位初始场和外域潮位边界瞬时值，运用fvcom模型计算获得外域三维分层潮流场，结合三线插值法得到每个网格质心处的潮流值；利用最邻近点插值方法，首先将三维分层潮流数据从非结构网格质心处传递到融合区中间结构化网格节点处，随后进一步将其传递到内域模型openfoam边界节点上；以到达内域模型边界处的潮流数据作为初始值，通过计算得到潮流与海洋构筑物相互作用结果。

在模型外域设置粗分辨率二维非结构网格，网格间距介于2m~200m之间，采用中尺度fvcom模型模拟大范围三维分层潮流运动，如图3所示。fvcom的原始方程主要包含动量方程、质量连续方程以及温度、盐度和密度方程等。在物理和数学上用mellor-yamada2.5阶垂向湍流闭合模型及smagorinsky水平湍流闭合模型对方程组进行闭合。在垂向上使用σ坐标系或通用垂向坐标系对不规则底部地形进行拟和，在水平上利用非结构三角形网格对水平计算区域进行空间离散。模型计算完成后，对三维分层非结构化网格节点数据进行三线插值，得到每个网格质心处的潮流数据，包括流速u、水位h、湍流动能k等。

融合区设置中间网格,将fvcom模型计算得到的三维分层潮流数据通过最邻近点插值方法传递到中间网格上，中间网格为结构性网格，即耦合作用区，如图4所示，从点c(u,h,k)到点c(u',h',k')过程。

对openfoam模型控制的小尺度内域进行精细化处理，三维网格的每个网格间距离尺度为0.1m~10m之间。为了更好的捕捉构筑物附近流体动力变化，一般可根据海洋构筑物的特征，将三维网格划分为结构性网格或非结构性网格，如图3所示，因此三维网格可以是四面体、六面体，或者它们之间的组合。

采用最邻近点插值方法将潮流数据从中间网格传递到openfoam模型控制的内域边界上，即图4中点c(u',h',k')到点c(u'',h'',k'')过程。利用swak4foam工具将达到三维网格边界的潮流数据进行初始化赋值，并采用interfoam作为内域模型求解器。

在fvcom模型与openfoam模型潮流数据传递过程中，中间网格的引入，提高了fvcom模型与openfoam模型间解决方案的传递效率，避免了如果fvcom模型为非结构性网格，openfoam模型同样为非结构性网格时，两个模型间直接数据传递引起的过程复杂、时间过长问题。

内域模型openfoam以有限体积法fvm离散模型空间网格，以eular方法离散时间项，将计算域离散成单元体后，物理量分布在单元网格的质心上。将单元体上的积分方程转化为用单元体中心物理量表示的代数方程，最后得到方程组数值解。

同时，使用流体体积函数（vof）方法捕捉复杂的水面变化，并采用gaussmuscl格式求解体积分数方程中的对流项，此外为更好的重构水气界面，模型中引入了gaussinterfacecompression格式的界面压缩项。潮流传播过程中的速度压强的耦合计算采用pimple算法，离散后的线性方程组采用预条件双共轭梯度法求解。

通过内域模型openfoam计算，得到潮流与海洋构筑物的相互作用结果，解析大范围海洋动力作用下，局部小尺度构筑物附近的水体飞溅和冲击现象。

以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。