[论文解读] A Numerical Formulation for Simulating Free-Surface Hydrodynamics

Thomas T. O’Shea, Kyle A. Brucker+2 more|arXiv (Cornell University)|Oct 7, 2014

Fluid Dynamics Simulations and Interactions参考文献 7被引用 23

一句话总结

本文提出一种基于结构化笛卡尔网格的有限体积法，结合浸入边界法与体积分数法（VOF），用于模拟船舶周围自由表面流体动力学，特别关注船尾平板形尾流。该方法能准确捕捉波浪破碎、飞溅及空气卷入现象，数值预测结果在多个拖曳速度下与实验室测量的阻力和波形剖面数据高度一致。

ABSTRACT

Cartesian-grid methods in combination with immersed-body and volume-of-fluid methods are ideally suited for simulating breaking waves around ships. A surface panelization of the ship hull is used as input to impose body-boundary conditions on a three-dimensional cartesian grid. The volume-of-fluid portion of the numerical algorithm is used to capture the free-surface interface, including the breaking of waves. The numerical scheme is implemented on a parallel computer. Recent improvements to the numerical scheme are discussed, including implementation of a new multigrid procedure and conversion to MPI communication. Numerical predictions are compared to laboratory measurements of a towed transom-stern model.

研究动机与目标

• 开发一种鲁棒、即插即用的CFD框架，用于模拟船舶周围复杂的自由表面流，包括波浪破碎与空气卷入现象。
• 通过使用带有浸入边界条件的笛卡尔网格，消除对边界拟合网格的需求。
• 通过新型多重网格求解器与基于MPI的通信机制，提升数值效率与可扩展性，以支持大规模模拟。
• 将数值预测结果与不同拖曳速度下平板形尾模型的详细实验室测量数据进行验证。
• 实现对大规模海军流体动力学问题的高分辨率模拟，包括十亿网格点量级的应用。

提出的方法

• 采用结构化交错笛卡尔网格求解不可压Navier-Stokes方程，时间积分采用分数步长格式。
• 利用体积分数函数在单元内追踪自由表面界面，通过平滑阶跃函数表示清晰的液气界面。
• 通过体力量方法施加浸入边界条件，利用表面面板化输入数据，在船体表面强制实现无滑移与无穿透边界条件。
• 实现一种新型多重网格求解器，以加速分数步长方法中压力泊松方程的收敛。
• 集成基于MPI的通信例程，提升高性能计算系统上的并行可扩展性。
• 采用亚格子尺度模型处理湍流与黏性应力，其中SGS应力张量采用动态过程建模。

实验结果

研究问题

• RQ1 基于VOF与浸入边界技术的笛卡尔网格方法能否准确模拟船舶模型周围的波浪破碎与船尾平板形尾流？
• RQ2 在不同拖曳速度下，数值预测的阻力与波形剖面与实验室测量结果的吻合程度如何？
• RQ3 新型多重网格求解器在大规模自由表面流CFD模拟中对收敛速率的提升程度如何？
• RQ4 数值预测的自由表面高程与水流粒子速度与 plunging breaking wave 的实验数据相比如何？
• RQ5 在解析高频表面谱特征方面存在哪些局限性？这些局限性与混叠或模型简化有何关联？

主要发现

• 平板形尾模型的阻力数值预测结果在四个拖曳速度下与DTMB实验室测量值高度一致。
• NFA代码预测的波形剖面在所有四个速度及横向位置均与实验数据极为吻合，尤其在船尾区域保持了稳定的精度。
• 对 plunging breaking wave 的模拟所得自由表面高程在形状与量级上均与实验室测量结果高度一致。
• 该模型成功捕捉了波浪破碎、飞溅形成与空气卷入等复杂流动特征，尤其在中等速度（7–8节）时表现突出。
• R/V Athena船后方的鸡尾浪区域，平均自由表面高程预测与实验结果吻合良好，但高频成分存在偏差，可能源于混叠或数值模型中缺少附体结构。
• 新型多重网格求解器与MPI通信机制的引入显著提升了求解器收敛速度与并行可扩展性，使模拟分辨率接近十亿网格单元。