基于势流理论和OpenFOAM的耦合模型对多向不规则-水科学进展.PDF

第２７卷 第 １期 水 科 学 进 展 Ｖｏｌ．２７，Ｎｏ．１ ２０１６年１月　 　 ＡＤＶＡＮＣＥＳ ＩＮ ＷＡＴＥＲ ＳＣＩＥＮＣＥ Ｊａｎ．，２０１６ ＤＯＩ： １０１４０４２／ ｊｃｎｋｉ３２１３０９２０１６０１０１０ 基于势流理论和ＯｐｅｎＦＯＡＭ 的 耦合模型对多向不规则波浪的模拟 １，２ １ ２，３ ２ 季新然 ，柳淑学 ，Ｂｏ Ｔ ＰＡＵＬＳＥＮ ，Ｈａｒｒｙ Ｂ ＢＩＮＧＨＡＭ （１ 大连理工大学海岸和近海工程国家重点实验室， 辽宁 大连　 １１６０２４； ２． 丹麦技术大学机械学院， 灵比　 ２８００； ３ 代尔夫特三角洲研究中心， 代尔夫特　 ２６２９ＨＤ） 摘要： 基于势流理论的数值水池可以快速计算波浪的传播及其对建筑物的作用， 但是势流理论是基于波浪的无黏性 假设的， 而在工程中， 通常需要在固体边界及波浪破碎的区域考虑黏性效应。 针对基于求解Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程的黏性 水池计算量较大、 速度较慢的缺点， 采用耦合的方法模拟多向不规则波浪的传播， 即在外域通过基于势流理论的数值 水池产生多向不规则波浪， 内域采用求解Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程和流体体积法（ＶＯＦ方法）对自由表面进行追踪， 通过外 域所提供的边界波浪， 内域计算可以在较小区域进行计算，从而达到减少计算量、 提高计算效率的目的。 关键词： 势流理论； ＯｐｅｎＦＯＡＭ； 耦合； 多向不规则波浪 中图分类号： Ｏ３５３２　 　 　 文献标志码： Ａ　 　 　 文章编号： １００１⁃６７９１（２０１６）０１⁃００８８⁃１２ 随着高性能计算机的发展和数值模拟技术日趋成熟， 数值水池成为海岸及海洋工程中波浪问题的重要研 究手段。 根据解决问题的不同， 一般将水池分为两种： 一种是基于无旋、 无黏假设的完全非线性势流理论 （ｆｕｌｌｙ ｎｏｎｌｉｎｅａｒ ｐｏｔｅｎｔｉａｌ ｆｌｏｗ ｍｏｄｅｌ）； 另一种是基于考虑黏性的Ｎａｖｉｅｒ⁃Ｓｔｏｋｅｓ方程。 势流理论建立的数值水池是基于无旋、 无黏的势流假设， 一般采用边界元法离散拉普拉斯方程和非线性 的自由面边界条件， 采用混合欧拉⁃拉格朗日方法来更新自由面的运动。 这类水池一般计算量比较小， 理论 ［１］ 和数值技术比较成熟。 Ｔａｎｉｚａｗａ等 在数值水池中模拟了波浪、 船和抗摇舱内水体间的耦合运动， 将实验、 ［２］ 理论分析和数值模拟３种手段相结合， 以优化抗摇舱的设计。 Ｇｒｉｌｌｉ等 建立了高精度的三维数值水池， 并 模拟了孤立波的爬升翻卷， 还利用势流理论水池和黏性流体水池各自的优点， 结合边界元法（ＢＥＭ）模型和 ［３］ ［４］ 流体体积法（ＶＯＦ方法）研究了波浪的破碎 。 Ｂａｉ和Ｔａｙｌｏｒ 建立了高阶边界元模型， 采用混合的欧拉⁃拉格 朗日法对数值波浪水槽内直立圆柱与非线性规则波浪和聚焦波作用时波浪的绕射及传播进行了研究。 房克照 ［５］ 和邹志利 建立了基于高阶Ｂｏｕｓｓｉｎｅｓｑ 水波方程的一维波浪破碎数值模型， 对波浪在均匀海岸上的浅化以及 ［６⁃７］ ［８］ 破碎过程进行了模拟。 Ｍａ等 、 Ｗａｎｇ和Ｗｕ 采用有限元的方法建立了完全非线性的数值水池， 模拟了 ［９］ ［１０］ 非线性波浪与圆柱的作用。 Ｅｎｇｓｉｇ⁃Ｋａｒｕｐ等 将Ｂｉｎｇｈａｍ和Ｚｈａｎｇ 提出的二维变阶有限差分方法的完全非 线性