简介：一份基于Ansys Fluent的船体CFD瞬态仿真案例，聚焦船体在逆海波浪中移动时的升沉与纵摇响应。内容采用VOF模型求解明渠流动，通过明渠波浪边界条件生成浅层波浪，并借助动态网格模拟Wigley船体运动，同时利用UDF将运动限制为4个自由度，且启用数值海滩选项抑制出口波浪反射。适用于船舶CFD工程师、相关专业学生，以及需要开展船舶耐波性瞬态仿真分析的研发人员。整个文件包共2个文件，包括Fluent案例文件（.c）和配套的网格/数据文件（.gz），压缩包约42.4MB，整体结构清晰，便于直接导入运行。已有162人浏览学习。下载后可对照案例文件完整理解波浪生成、动态网格运动、自由度约束、无反射边界等关键设置的实现方式，快速上手同类船舶波浪耦合仿真问题，也可在此基础上修改参数开展进一步研究。
1. 船体CFD瞬态仿真：为什么定常结果到了波浪里总翻车
做船体CFD仿真的人大多有过这个经历：用定常求解算出来的阻力、纵倾角，放到真实海况里一对比，对不上。某模拟项目X里，我拿定常结果的兴波阻力系数去估一条排水艇的顶浪失速，结果预报值比实际偏了15%以上。问题不在网格也不在湍流模型，而是定常求解根本模拟不了“船在波浪中不断变化的水线面、不断变化的浮态”这回事。
瞬态仿真就是为这个场景准备的。它把时间项、自由液面、船体六自由度运动全部耦合在一起，船在波峰上被抬起、在波谷里下坠、纵摇俯仰、垂荡沉浮，阻力和姿态随时间真实变化。不过，瞬态仿真把网格、时间步长、动网格策略的难度成倍放大，新手最容易在这三块里卡死。
这篇笔记面向已经会跑Fluent稳态算例、但想转到波浪中船体运动仿真的工程师。目标是让你在上手之前，先搞清楚船体瞬态仿真的模型选择逻辑、网格策略、时间步长控制，再把最常见的翻车点提前避开。从模型选型讲到后处理验证，每一节都能直接落到你的项目里。
2. 模型选型的思路：VOF、DFBI与湍流模型的取舍
船体瞬态仿真涉及三组核心模型选择：自由液面用什么捕捉、船体运动怎么解耦、湍流用什么闭合。这三组选型互相耦合，后面任何一个参数出问题，都可能是最开始模型选错埋下的雷。
2.1 自由液面模型：VOF与界面压缩的配合
船体CFD里的自由液面，通常用VOF（Volume of Fluid）模型来捕捉。它的基本思路很简单：每个网格单元记录一个水相体积分数，0代表全空气、1代表全水、0到1之间代表两者交界面。Fluent里设置路径是Models → Multiphase → Volume of Fluid，需要开启Implicit Body Force，否则重力与压力梯度的平衡靠数值硬扛，会在自由液面附近产生虚假速度。
VOF模型的关键参数不在这两行菜单里，而在离散格式。压力插值要选Body Force Weighted（PRESTO!也行，但Body Force Weighted在自由液面重力主导的问题里更稳）；动量方程里的体积分数插值，建议开Compressive界面压缩格式。对比一下，默认的Geo-Reconstruct（几何重构）精度高，但每步迭代里要重建界面，计算成本大；Compressive的数值扩散稍大，却能在同样的网格密度下维持更稳定的自由液面形状。对于船体这种动边界不断变化的问题，我一般选Compressive，稳定优先。
需要认识的一个深层问题是：VOF对网格的依赖比对离散格式的依赖更大。自由液面处的网格必须接近水平分辨率才能保住“波浪不被数值抹平”。一条垂向网格尺寸是水平尺寸两倍的网格，经过一个静水衰减算例就能明显看到波浪高度被压掉。做瞬态船体仿真之前，先用一个单独的静水网格跑恒定波高衰减验证，这是最便宜的自由液面质量检查。
设置时还要注意两相各自的物理属性。空气和水都要设置参考密度，默认参考密度在1 kg/m³以下时压力基求解器会出现可压缩性伪影。Fluent的Operating Conditions里把Operating Density设置成空气密度，而不是默认0，能大幅改善压力收敛和自由液面处速度。
2.2 船体运动自由度：DFBI模型与约束策略
瞬态船体仿真要让船动起来，Fluent里通常用DFBI（Dynamic Fluid Body Interaction）模型。它把船当成刚体，在每一个时间步先求解流场，再通过表面压力/剪切力积分得到力和力矩，解刚体运动方程更新位置和姿态。路径是Dynamic Mesh → Dynamic Mesh Zones → Six DOF（或新版菜单里的DFBI）。船体自由度的完整释放包括纵荡、横荡、垂荡、横摇、纵摇、艏摇，但在大多数规则波顶浪算例里，真正起作用的是垂荡和纵摇两个自由度，其他自由度受对称性和航行稳定性约束，一般锁死。
这里有一个容易忽略的点：DFBI的旋转轴位置。很多人直接用船体重心坐标做旋转中心，但注意你要输入的是船体质量、惯性矩，而且惯性矩是相对质心的。如果船体几何原点不在重心，惯性矩张量的坐标也需要平移到几何原点上，不换算的话纵摇频率会偏差很大，误差直接体现在纵摇周期的偏移上。
自由度的约束不是一次性设置完的。收敛性差的算例里，可以先用四个自由度（垂荡、纵摇、纵荡、艏摇）试试，然后逐步加载全部自由度。DFBI还有阻尼选项，在Motion Attributes里能附加转动阻尼，这个选项不是用来模拟物理阻尼的，而是用来让数值稳定性差的算例“拖一拖”，让它不至于在初始几个时间步直接发散。阻尼系数不能给太大，否则船体运动响应被扭曲。常见做法是先给一个很小的阻尼系数，比如临界阻尼的约5%，等瞬态稳定后再逐段递减到0。
2.3 湍流模型：从k-ε到SST k-ω的现实考量
湍流模型选择不是越大越好。很多做船体仿真的人一上来就用LES或DES，结果算力根本扛不住，而且非定常的LES要统计足够长的时域才有意义。对瞬态船体仿真，SST k-omega是工业上的首选。
SST k-omega的优势在于边界层内用k-omega的鲁棒性，边界层外和自由剪切层用k-epsilon的优势，并在两者之间通过混合函数平滑过渡。船体表面附近的逆压梯度、流动分离模拟都需要这种过渡能力，用标准k-epsilon做船尾伴流场时分离点偏晚、阻力偏小，这在阻力和耐波性联合预报里很致命。
在Fluent里选择SST k-omega之后，需要特别关注Near-Wall Treatment。Fluent里默认是Automatic Wall Treatment，它试图在不同y+区域自动适配壁面函数。但船体仿真里适合的是Low-Re Correction（低雷诺数修正），强制启用后壁面网格y+在1左右，直接用粘性底层解析。代价是壁面网格必须足够密，近壁第一层网格高度要到微米量级（跟船体尺度、来流速度相关，后面给出估算方法）。
另外一个容易忽略的选项是Curvature Correction（曲率修正）。SST模型默认不感知流线弯曲，水面船艏部曲率变化剧烈，开启后对艏部流动分离模拟有明显改善，计算量增加很小，建议默认开启。有人在这里用Transition SST模型，它在小尺度滑行艇上效果不错，实尺度排水船反而经常因为转捩位置不对而预测偏航。所以明确一下：常规排水船瞬态仿真选SST k-omega就够，不要盲目上转捩模型。
3. 从几何到网格：瞬态仿真需要什么样的计算配置
网格是瞬态仿真的地基，地基歪了，后面求解器再稳也是空中楼阁。船体几何处理、网格拓扑、近壁面网格估算、动网格策略这四件事必须在前处理阶段一次做对。
3.1 船体几何清理与计算域尺寸设定
船体模型导入Fluent前需要把几何拓扑修整干净。来自CAD的船体模型常见问题包括：曲面间微小的间隙/重叠面、船艉的压载水舱或附体小面、螺旋桨轴和舵系的非水密结构。这些在Fluent Meshing里可以用Share Topology或Intersect/Subtract处理。船体表面的网格决定水线捕捉质量，不能有穿透或间隙。
计算域的设计有讲究。三维船体在波浪中的兴波会向各个方向传播，域外边界反射波会污染算例。顶浪规则波算例，入口距船艏至少1.5倍船长，出口距船艉至少3倍船长，两侧各2倍船长，顶部离静水面至少1倍船长，底部水深至少2倍船长。倍数是针对最大波长和船体尺度取较大的一个。底部用Speed Inlet或压力出口都会产生边界反射，原则上自由液面区域的六个边界都要用能吸收波的边界条件。
实际项目中我常用数值消波功能。Fluent的Wave BC插件（或用户自定义函数）可以在出口和侧边界的自由液面区域添加动量源项，用于吸收外行波。如果不用这个，出口处会看到驻波反射上游，纵摇幅值会出现周期性脉动叠加，数据后期也很难完全滤掉。
3.2 近壁面网格：y+估算与第一层网格高度计算
第一层网格高度是由边界层解析需求决定的。船体CFD中，y+目标在1左右，需要估算摩擦速度uτ，第一步估算摩擦系数Cf，再用Cf计算壁面剪应力，进而得uτ。公式路径为：Cf ≈ 0.075 / (log10(Re) - 2)^2，τ_w = 0.5 * ρ * V^2 * Cf，uτ = sqrt(τ_w / ρ)，y = y+ * μ / (ρ * uτ)。这个公式处理的是平板边界层近似，船艏和船艉的局部流动偏差较大，但对于第一层网格高度估算已经足够。
下面是一个估算脚本，可以用Python快速算不同航速下的第一层网格高度：
import math rho = 1025.0 # 海水密度 kg/m^3 mu = 0.0012 # 动力黏度 Pa·s V = 10.0 # 航速 m/s L = 50.0 # 船长 m y_plus = 1.0 # 目标 y+ Re = rho * V * L / mu Cf = 0.075 / (math.log10(Re) - 2.0) ** 2 tau_w = 0.5 * rho * V * V * Cf u_tau = math.sqrt(tau_w / rho) y = y_plus * mu / (rho * u_tau) print(f"Re = {Re:.2e}, Cf = {Cf:.4f}") print(f"u_tau = {u_tau:.4f} m/s, 第一层网格高度 = {y * 1000:.4f} mm")这段脚本的逻辑是：从雷诺数出发估算摩擦系数，再算出壁面剪应力和摩擦速度，最后用y+定义反推第一层网格物理高度。结果出来通常是0.5到2毫米的量级（取决于航速和船长）。如果输出的第一层高度小于0.1毫米，说明边界层网格数会非常大，要考虑分区划分或局部加密，而不是整体铺满。
Fluent Meshing里设置边界层的方式可以选Thin/Prism层。船体表面需要5到10层边界层，增长率1.1到1.2。还要注意边界层网格在自由液面附近不要与液面网格垂直交错，否则界面捕捉和边界层同时抢网格，体积分数收敛会互相干扰。
3.3 动网格策略：重叠网格与动态层更新的取舍
船体大幅垂荡和纵摇时，网格必须跟着动。两个主流选择是重叠网格（Overset）和动态层更新（Dynamic Layering）。
重叠网格的逻辑是把计算域分成背景网格和船体随动网格，两块网格各自独立划分，再通过插值耦合。船体运动幅值大、姿态变化剧烈时，重叠网格的优势明显，因为船体附近的边界层网格可以刚性跟随船体运动，不会因为大幅变形而扭曲。代价是数据交换层的插值误差。
动态层更新则是在船体边界的相邻区域动态生成/销毁网格层。船体下沉时上部网格层压缩并合并，船体抬升时下部网格层分裂增长。对垂荡幅值小的算例它足够，但船体纵摇角度超过几度，网格单元很容易被拉伸到极扁，单元质量急剧下降，发散风险很大。
从实际角度，常规规则波中船体运动幅度在0.1倍船长以内时，动态层更新足够。大波高、强非线性碎波场景，直接选重叠网格。要留意的是重叠网格的插值区域不能在自由液面大梯度区太薄，否则速度插值会把界面抹平。给一个量级判断：插值边界与船体表面的距离，至少要大于边界层总厚度，同时避开自由液面剧烈波动的区域。
4. 瞬态求解设置：从时间步长到收敛判据
模型选好、网格做好，接下来是瞬态求解器的核心调度：压力速度耦合、离散格式、时间步长、收敛判据。这一节有大量“差一点就全线发散”的细节，值得多花时间看。
4.1 压力速度耦合与瞬态离散格式选择
船体CFD是低速不可压流动，压力基求解器（Pressure-Based）是默认选择。压力速度耦合在Fluent里通常用Coupled Scheme，它的优势在于同时求解动量和压力修正方程，在时间步长较大的情况下比SIMPLE系更稳定。代价是每步迭代内存消耗增大，但对瞬态船体问题来说稳定优先。
时间离散格式，用户手册里写的是First Order和Second Order Implicit。瞬态船体问题第一选择是Second Order Implicit。它具备无条件稳定特性，但实际使用中在积分初始阶段要小心——开始几个时间步内压力冲击大，二阶格式对初始波动响应敏感，可能震荡。我一般会在前20个时间步把时间离散降为一阶，等压力场稳定后再切回二阶，这样能避免“一起步就发散”的常见问题。
对流项的离散格式选择也要同步调整。动量方程用Second Order Upwind，湍流方程也用Second Order Upwind，VOF的体积分数方程用Compressive。压力插值用Body Force Weighted。不要把默认值直接沿用，默认的中心差分格式会在对流占主导的船体绕流里引入数值震荡。
4.2 时间步长与CFL数：经验值先行
时间步长是瞬态仿真最重要的旋钮，它直接决定库朗数（CFL）：
CFL = V * dt / dx
V是局部流速，dt是时间步长，dx是网格尺寸。船体CFD里自由液面处的CFL尤其关键。经验法则是自由液面区域的CFL小于1，最大不超过2。网格最小尺寸是水面处垂直方向加密后的尺寸，如果网格在自由液面处是0.1米，船速5米/秒，那么dt约为0.02秒。这会导致每秒钟要跑50步，一步假设10到20次迭代，一分钟就是数千次迭代——瞬态仿真的算力消耗就是这么被抬起来的。
另一种基于波浪的时间步策略：规则波周期T，一个波周期建议至少采样200到400个时间步来捕捉波形演化。船的垂荡纵摇响应频率约等于遭遇频率，为了捕捉响应幅值和相位，每周期至少100步。两个条件取更严格的那个。
Fluent的Adaptive Time Stepping可以自动调整dt，但我不建议一上来就用。理由很简单：船体运动刚开始瞬态响应时，速度变化快，自适应会频繁进退，产生额外的插值误差。先用固定步长确定流动特征的稳定趋势，再开启自适应来优化中段效率。一般来说，自适应之后网格尺寸较小区域的CFL保守值设在0.5，会有比较平稳的推进效果。
4.3 收敛判定：看残差不如看积分量
瞬态仿真的收敛判定和定常非常不同。每一时间步内会迭代若干次，看的是每一步内的残差是否降到目标。常见做法是设残差下降3个量级或达到1e-4。但船体CFD里残差降到多少并不是终极标准，监测量的稳定才是。
在每个时间步内，要同时监测船体阻力的变化曲线、船体垂荡位移、纵摇角的变化。如果这三个量在该时间步内已经不再变化，说明流场和刚体运动已经耦合收敛。有一个经验值：单步内船体运动监测量的相对变化小于0.1%后再进入下一步。Fluent里有几类瞬态求解推进方式，其中可以选择时间步内自动迭代，比如每步最多迭代20次，但最好是设置一个收敛判据，满足即提前进入下一步。
开启自动保存和自动记录：Data File Quantities里要选上Wall Shear Stress、Pressure Coefficient、Volume Fraction和DFBI的刚体运动量。自动保存在瞬态里通常设置为每10步或每0.1秒保存一次，以便回查故障点和重启动。磁盘开销建议提前规划，单算例几十GB很正常，要预留空间。
5. 常见问题与排查：船体瞬态仿真的典型翻车点
这里整理的是几个反复出现的典型问题。以下每条都按现象、原因、解决的逻辑展开，都是我实际调试时硬碰出来的经验。
5.1 初始几步直接发散：压力和速度像是爆炸一样增长
现象：算例在刚开始的十几步内出现NaN，或者压力场出现非物理的极大值，流场直接崩溃。
原因：初始流场和船体运动不一致，压力速度耦合在初始状态下产生强冲击。另一个常见原因是船体在初始位置有局部几何不良——船艏底部网格有负体积，或者自由液面初始体积分数分布没有与船体水线面对齐。
解决：把初始化的船体位置放在静水水线面上，确保体积分数初始化的水面与船体水线面吻合。用Fluent的Adapt → Region来强制船体内部区域设定为水相。前20到50个时间步，时间离散用一阶、时间步长缩到目标步长的四分之一到二分之一慢慢放量。如果仍然发散，检查船体附近网格单元质量，特别是最小正交质量低于0.1的区域要重新处理。
5.2 波浪在传播中快速衰减：波幅越传越小甚至消失
现象：入口给定规则波，但波往计算域深处传播时波高明显下降，到船体时已不足给定的70%。
原因：网格分辨率不足以承载波浪。自由液面垂向网格尺寸过大，或水平方向波浪每波长网格数不足。一般规则波需要每波长至少40到60个网格单元，波高方向至少10到15个网格单元。如果网格够但波仍衰减，检查入口的波边界条件参数，比如波高和周期是否输入正确、入口边界是否启用了数值消波。
解决：在船体周围的自由液面区域做局部加密。如果加密后网格数无法接受（总单元数到千万级），考虑缩短计算域，把消波区压缩，让有效传播区域内的网格尽可能保留在波高分辨率上。检查数值格式：VOF用Compressive，时间离散用二阶，一阶时间离散会造成明显数值耗散。
5.3 船体运动发散：垂荡幅度越来越大直到撞出计算域
现象：船体在波里像蹦床一样越跳越高，垂荡位移单调增大，最终冲破动网格更新范围。
原因：初始状态与波作用力不平衡。船体浮态设置错误，重心位置不对，或者垂荡刚度（水线面面积决定的恢复力）与实际不符。另一个常见原因是时间步长过大，船体在每个时间步内位移过大，动网格更新产生的网格变形过大，导致力计算失真，进一步放大运动。
解决：回查三件事——初始浮态是不是静浮平衡；重心位置和转动惯量是否与排水量匹配；时间步长的CFL是否高于2。然后把阻尼系数临时加上去，比如在纵摇通道加5%临界阻尼跑200个时间步，让船体运动先稳定下来，再逐步解除。垂荡幅值收敛到稳定极限环之后，再加参数减阻尼。模拟中途船体位移过大导致网格更新失败，要回到前一个保存点，缩短时间步再试。
5.4 质量不守恒：进出口流量不平衡，阻力曲线跳变
现象：后处理发现进出口流量不平衡程度超过1%，阻力曲线在没有物理原因的情况下出现台阶式跳变。
原因：自由液面区域的进出口边界条件设置不匹配，或者VOF的界面在进出口处发生了数值回流。出口压力边界给定静水压力时，波面运动导致压力动态变化，出口边界反馈错误流量。
解决：进出口处采用Wave BC，且出口距船体尾部足够远。检查体积分数是否在所有边界上正确初始化，特别是入口和水面边界。另一个技巧是把质量守恒残差设成监测量，而不是只测全局残差。每步的网体积分数不平衡量应该低于0.1%。如果跳变发生在某特定时刻，回查那个时间步的保存文件，定位是哪一区域出现回流或界面破碎。
5.5 网格更新失败信息反复出现
现象：运行过程中Fluent报错Dynamic Mesh Update Failed，程序中止。
原因：动网格更新时，船体局部网格变形过大。尤其是艏部局部加密区、船体与自由液面交界处的边界层网格，这些区域的网格在变形后发生负体积翻转。另一个原因是动态层更新设置了过大的层高比，压缩/分裂不协调。
解决：如果是动态层更新，缩小Split Factor和Collapse Factor的触发值；如果是重叠网格，检查插值边界与船体表面的间距。更有效的方式是给船体运动区域单独包一层刚性运动网格区域，外部再通过插值或变形网格过渡。这个分层策略能让运动区域的网格完全刚性跟随船体，避免贴合船体的网格产生变形。
6. 进阶：做完网格无关性验证再谈结果可信
上一节解决了“算得动”的问题，这一节解决“算得准”的问题。在做大量工况参数扫描之前，必须先做网格无关性验证。
网格无关性验证的操作方法：保持几何、物理模型、湍流模型、时间步长不变，只改变网格密度。三套网格分别是Coarse（基准网格的0.5倍数量）、Medium（基准网格数）和Fine（基准网格数的2倍数量）。比较量选择三个：总阻力系数、垂荡幅值、纵摇幅值。如果Medium和Fine的偏差在2%以内，可以认为网格无关性满足。
具体步骤是：先跑Medium网格算一个规则波周期，记录垂荡和纵摇的时历；再跑Fine网格重复同样算例。比较两者的垂荡幅值相对偏差。如果偏差超过5%，说明网格太粗，需要加密自由液面网格而不是全计算域加密——重点加密波高方向的网格密度。注意时间步长也要在网格加密后重新检查CFL，因为dx变小了，CFL会相应增大。
后处理提取结果的逻辑也有讲究。瞬态仿真输出的是随时间变化的阻力和姿态，不能直接拿单一时步的结果说话。规则波里要取稳定后的多个周期做统计分析：用FFT或最小二乘拟合提取幅值和相位，还要看是否稳定（没有持续增长或衰减）。阻力系数是瞬时值与平均值并存的关系：平均值用于总阻力评估，瞬时波动用于分析推进系统的载荷谱。
稳定的瞬态结果可用后，导出数据。导出的格式要包含时间、阻力、垂荡位移、纵摇角。Python的pandas可以直接读Fluent导出的CSV做快速统计：
import pandas as pd df = pd.read_csv("ship_motion.csv") # 假设数据从第200秒后才是稳定段 steady = df[df["time"] > 200.0] heave_amp = 0.5 * (steady["heave"].max() - steady["heave"].min()) pitch_amp = 0.5 * (steady["pitch"].max() - steady["pitch"].min()) print(f"稳定段垂荡幅值: {heave_amp:.4f} m") print(f"稳定段纵摇幅值: {pitch_amp:.4f} deg")这段脚本先取稳定段数据，再用极差的一半估算幅值。要注意的是，极差法对噪声敏感，数据中有少量尖峰就会高估幅值。更稳妥的做法是先对稳定段做带通滤波再取幅值，但工程快评场景下极差法够用，只要确认时历曲线没有异常尖峰。
从那以后，每接手一个船体瞬态仿真项目，我不管前面的算例多着急，都强制走一遍静水衰减验证、网格无关性验证、固定步长稳定化这三道前置门禁。看起来磨蹭，实际上是在给后面所有工况试错兜底。先把基准算例跑稳了再谈参数扫描，比反复踩同一个坑要快得多。希望这些经验能帮到正在被瞬态仿真折磨的你，少走几步弯路。
本文还有配套的精品资源，点击获取
