简介面向船舶 CFD 学习者与 Fluent 中高级用户的瞬态仿真案例包使用 Fluent 模拟 Wigley 船体在逆海波浪中的升沉与纵摇响应。案例以 VOF 模型求解明渠流动用明渠波浪边界条件生成浅层波浪并通过动态网格和用户自定义函数UDF将船体运动限制为 4 自由度同时启用数值海滩选项抑制出口处的非物理反射。整套资源共 2 个文件包含 gz 格式的 Fluent 案例文件与 c 格式的 UDF 源码压缩包大小 42.4MB体量轻便适合直接下载运行或按需修改。已有 162 人学习下载。对希望掌握船体耐波性仿真、VOF 波浪建模、动态网格与 UDF 联用的读者来说这份案例包提供了可对照操作的实战参考能帮助快速理清从边界条件设定、动网格更新到数值稳定控制的完整流程是深入理解瞬态 CFD 仿真技巧的实用资料。1. 船体CFD瞬态仿真为什么稳态算不准必须上瞬态做船体水动力仿真的人八成都会遇到同一个困惑明明几何导入没问题、网格质量也过得去稳态算出来的阻力却和试验值差出一大截。原因不复杂——船是在水面上跑的姿态会变自由液面在动船体周围的流场本质上就是非定常的。稳态求解把船“钉死”在一个固定姿态里忽略了浮态变化和兴波的非线性效应结果自然偏。用Fluent做船体CFD瞬态仿真核心价值就是把船放回真实的水动力环境里让船体在重力、浮力和流体动力作用下自由纵摇、垂荡同时捕捉自由液面的瞬态演化。这篇笔记面向的是要用Fluent实际算船的人无论你是做船舶方案论证、螺旋桨干扰分析还是想验证船型优化后的阻力曲线都可以照着这套流程做。我会把物理模型选型、网格策略、瞬态求解参数、避坑经验一次讲透。2. 物理模型选型VOF、湍流模型和六自由度运动2.1 自由液面捕捉VOF多相流模型为什么够用船体CFD仿真首先要决定怎么处理“水-空气分界面”。Fluent里常用的是VOFVolume of Fluid多相流模型它对每个网格单元记录水的体积分数水1、空气0、界面处是0到1的过渡值。选它有三个理由一是对自由液面大变形比如船艏破浪、艉流翻卷的鲁棒性比Level Set好二是能和动网格、重叠网格直接配合三是计算量可控适合瞬态长时程模拟。选VOF时有几个参数必须设对。首先是明渠流动Open Channel选项——如果你关心的是船在无限水域中的航行而不是闭式管道里的流动就需要打开Open Channel Flow并给定自由液面的高度和来流速度。其次是相间作用水和空气的界面张力在船尺度下可以忽略不需要开启Surface Tension。最后是体积分数离散格式建议用Geo-Reconstruct或Compressive界面锐利程度明显优于一阶迎风。; Fluent TUI脚本开启VOF明渠流动 /define/models/multiphase/volume-of-fluid yes /define/models/multiphase/options/open-channel-flow? yes /define/models/multiphase/options/volume-of-fluid-options/interface-modeling-type geo-reconstruct这段TUI脚本是我常用的开局配置。第一行开启VOF模型第二行打开明渠流动第三行把界面重构格式设为Geo-Reconstruct。打开Open Channel Flow后Fluent会自动添加重力项并让你设置水深和波高等参数这是船体绕流和规则波模拟的基础。2.2 湍流模型SST k-omega在船体绕流中的适用边界船体CFD里湍流模型的选择直接影响阻力和边界层厚度预测的准确性。我调试过几个不同方案对比结果是SST k-omega综合表现最好——它在近壁区用k-omega的边界层解析特性在外场用k-epsilon对自由剪切层和尾流的适应性兼顾了船艏高压区和船艉逆压梯度区的计算鲁棒性。但SST k-omega有一个前提要求第一层网格的y值必须控制好。船体表面通常需要y ≈ 1意味着第一层网格厚度只有毫米级甚至更薄具体取决于来流速度和水的运动黏度。雷诺数Re UL/ν以一条12米船、航速3.5米/秒为例Re约4.2×10⁷按平板边界层估算第一层厚度在0.1mm量级。如果你的网格达不到这个分辨率最好退而求其次用Realizable k-epsilon加壁面函数而不是硬上SST算出假边界层。湍流强度入口边界参数我一般设为2%湍流黏度比大约5——这个数值对水槽来流和开阔水域都算合理。如果有实验数据可以直接用实验值没有就按经验值来。2.3 瞬态的核心DFBI六自由度刚体运动船在水面上跑自身姿态是运动的、未知的需要在求解过程中实时算出来。Fluent里对应的模块叫DFBIDynamic Fluid Body Interaction也就是六自由度刚体运动模型。它可以基于流体压力分布积分求出船体受到的力和力矩再结合刚体动力学方程解出船体的垂荡、纵摇角速度然后把姿态更新传给网格。DFBI设置里最容易被忽略的是旋转中心和惯性矩。旋转中心一般取重心位置这个不难但惯性矩如果填错了纵摇和垂荡的耦合频率就全错了横摇阻尼算出来更是离谱。如果找不到准确的惯性矩数据至少要用经验公式估算并记录假设值否则结果没有可比性。/define/models/dynamic-mesh yes /define/models/dynamic-mesh/options/six-dof yes开启DFBI有两种路径一种是靠动态网格的六自由度选项另一种是独立的DFBI界面。上面这两行TUI是走动态网格路径的配置。实践中我推荐用独立的DFBI模型原因是它和重叠网格的配合更顺滑运动幅度可以给得更大而不容易产生负体积。边界条件方面船体表面设为Wall且要在Wall的动量选项中启用运动Moving Wall / Specified Motion让船体表面跟随刚体运动一起移动。底部和侧边界设速度入口或压力出口顶部设压力入口前后方向按来流方向分别设速度入口和压力出口。2.4 为什么DFBI必须和瞬态求解器绑定DFBI解的是刚体运动方程包含速度和位移对时间的积分天然就是时间依赖的。有人问“用稳态求解器加DFBI行不行”——不行稳态求解器没有物理时间步进的概念刚体运动方程根本无法积分。这就是为什么标题里“瞬态”是硬约束不是可选项。另一个容易被忽略的点是DFBI运动需要设置松弛因子和限制器。松弛因子取0.7到0.95之间限制器用来限制单步内允许的最大位移变化防止系统发散。我一般会把最大位移限制在0.1米以内角位移限制在2度以内——宁可多算几步也别让船在第一个时间步就“飞”出去。3. 网格策略重叠网格和船体表面网格的那几条红线3.1 重叠网格船动网格不动的组合方案船体瞬态仿真里船在动而计算域不想跟着动最舒服的做法就是重叠网格Overset Mesh。原理很简单计算域拆成背景网格和组件网格两部分船体放在组件网格里组件网格可以在背景网格里平移旋转两套网格在交界面通过插值交换数据。用重叠网格的收益很明显船动幅度大时组件网格跟着一起平移网格畸变很小不会出现动网格那种越拉越歪、最终负体积极的翻车现场。船体小幅升沉纵摇时甚至可以让组件网格随船刚体运动背景网格完全不动计算稳定性和收敛性都好了不少。上面这个空格不是代码是我画的一个计算域布局思路背景网格沿船长方向取3-5倍船长横向取1.5-2倍船长组件网格包裹船体边界离船面至少0.5个局部船长。重叠区要保证至少有5层网格互相嵌套否则插值误差会直接污染船体表面压力分布。3.2 边界层网格参数怎么定船体表面的边界层网格我一般用18到24层棱柱层增长率1.2第一层高度按y目标反算。目标y1时可以用Fluent的y计算器输入特征长度L和来流速度U就能反推第一层高度。表不同y目标对应的第一层网格高度以12m船、3.5m/s为例y目标第一层高度(mm)壁面处理方式适用湍流模型≤10.08-0.12低雷诺数壁面处理SST k-omega30-803-8壁面函数Realizable k-epsilon这里想提醒一点船体曲面复杂曲率大的地方网格会自动加密第一层高度要按局部流速重新评估。船艏驻点附近流速接近零y自然偏低船艉流动加速区流速偏高同样的第一层高度可能y会暴涨。所以网格质量控制不能只看最小正交质量还要盯船体表面的y分布。3.3 网格无关性验证的最低标准做瞬态仿真最怕网格量不收敛算完后心里没底。我习惯的做法是准备三套网格粗网格约300万单元、基准网格约600万、加密网格约1200万分别算同一个工况的阻力看总阻力系数变化趋势。如果粗网格和基准网格差超过5%说明网格还太粗至少要加密到基准和加密网格之间。要注意的是瞬态仿真的网格无关性验证必须使用相同的物理时间长度和时间步长否则对比的是“不同算例”而不是“不同网格”。这个坑我踩过——第一次做对比时粗网格算了8秒物理时间加密网格只算了4秒阻力值差了一大截还以为网格不收敛后来才发现是时间推进长度不一致。4. 瞬态求解参数时间步长、离散格式和收敛判据4.1 时间步长的确定从CFL条件到物理时间覆盖瞬态仿真的时间步长选择本质上是精度和成本的平衡。CFL条件要求Δt ≤ Δx/U其中Δx是最小网格尺寸U是局部流速。在边界层网格最薄的地方Δx只有0.1mm量级按CFL条件算出来的Δt小得离谱但边界层内流速也低实际的CFL没有那么苛刻。工程上的做法是让船体周围的自由液面网格区域满足CFL≈1到2。以背景网格2 cm、来流3.5m/s为例Δt大约在0.005到0.01秒之间。如果船艏艉兴波区域网格加密到5mmΔt要再缩小。我一般先用相场界面区域的网格尺寸除以流速估计一个初值然后实际迭代看看残差和自由液面是否出现数值振荡。物理时间长度方面至少要覆盖流动“穿越”计算域的时间T L_char / UL_char取3到5倍船长。如果模拟总时间小于这个值船艏兴波还没传到出口边界尾部流场也没到准周期状态阻力系数肯定会一直漂。4.2 离散格式和松弛因子瞬态仿真不能照搬稳态瞬态仿真的空间离散我推荐VOF用Geo-Reconstruct动量方程用二阶迎风压力用PRESTO或Body Force Weighted。这组搭配在船体绕流里表现稳定界面锐利度足够压力梯度精度也够。瞬态时间推进用二阶隐式Second Order Implicit这是瞬态精度的底线。一阶隐式虽然稳但数值耗散大自由液面上小尺度波浪会被抹平而船体兴波里的短波成分恰恰是阻力预测和尾流结构的关键。松弛因子要特别注意——瞬态下压力松弛因子可以比稳态大一点。我常用的设置是压力0.3、动量0.7、k和omega都是0.8。如果你用PISO算法注意把邻近修正Neighbor Correction和偏斜修正Skewness Correction都打开能有效缓解网格略歪时的压力振荡。压力速度耦合我用PISO而不是SIMPLE。原因是在瞬态推进里PISO的分步校正对时间精度更友好SIMPLE在一阶时间精度下容易把时间误差和压力误差混在一起。特别是船体有DFBI运动时压力场每步都得算准否则刚体运动模型不收敛。4.3 自由液面下的收敛判据残差、阻力系数和液面波动一起看瞬态残差设定的标准做法是把所有方程的残差收敛标准设在10⁻³量级。但残差只是必要条件不代表流场已经“准周期稳定”。要看三个东西同时到位一是阻力系数随时间曲线进入周期性或准稳态不再有单调漂移。二是船体纵摇角幅值稳定在某个区间不持续增长或衰减。三是自由液面波高在监控点处波形重复性好周期稳定。这三个信号都到位才能开始取样。表瞬态船体CFD常用监控量及评估标准监控量评估方法达标标志总阻力系数Ct时均值时间长度为3-5个穿越周期时变化1%纵摇角时程曲线出现稳定振荡且幅值不超过1-2度自由液面波高单点波形相邻两周期波形重合度高在做数据取样时只取后半个模拟周期的时均值前半个周期算是“预运行期”让流场从初始条件过渡到真实物理状态。初始条件如果是均匀来流头几个周期的阻力曲线往往会有明显的“启动冲击”这部分绝不能计入统计。4.4 一个可以直接抄的求解设置模板下面是项目里实际用过的设置顺序可以直接照着改参数运行; 求解器设置模板Fluent Scheme脚本片段 ; 时间推进 /define/models/solver/transient yes /define/models/viscous/sst-k-omega yes /solve/set/time-step 0.008 /solve/set/max-iterations-per-time-step 30 ; 离散格式 /solve/set/discretization-scheme/momentum second-order-upwind /solve/set/discretization-scheme/pressure body-force-weighted /solve/set/discretization-scheme/volume-fraction geo-reconstruct ; 收敛残差 /solve/set/eqn-residual-criterion continuity 1e-4 /solve/set/eqn-residual-criterion k 1e-4 /solve/set/eqn-residual-criterion omega 1e-4这几行里最值得说明的是“每步最大迭代30次”这个参数。瞬态求解时每步内部迭代次数不宜设太大——如果某一步老是不收敛说明时间步长可能太大了硬算30步也未必收敛还浪费时间。典型现象是每步都需要28、29次迭代才能打到残差标准那就把Δt缩小一半再看。时间步长0.008秒对应上面的12米船、背景网格2cm工况。如果你网格加密到1cm时间步长顺手减半到0.004保持CFL不变。Changing只改Δt不动网格收敛表现会很差。5. 避坑排查船体瞬态仿真最容易翻车的五个场景5.1 现象计算发散残差飙升到10¹⁰原因压力出口回流量太大自由液面在出口处出现“回流猛击”压力-速度耦合算不下去。最常见是出口边界贴船太近尾流区压力还没有恢复到静压就被出口边界截断。解决把压力出口往船艉下游推远到至少2倍船长如果还不行给出口加一个阻尼区Numerical Beach在阻尼区内逐渐降低水体积分数的物理影响人为吸收反射波。Fluent里可以通过用户自定义源项或在明渠流动选项里开启Damping Zone来实现。5.2 现象船体纵摇角持续增大最终数值崩溃原因DFBI运动的角速度积分有累积误差而限制器没有发挥作用或者是船体重心位置设错导致流体动力产生本来就该有的正反馈力矩。解决先检查重心坐标和旋转中心是否一致再降低DFBI松弛因子到0.5左右试试如果还不行把时间步长减半——纵摇是特征频率很高的自由度时间步长太大时每个步进的角增量过大运动方程数值不稳定。5.3 现象自由液面模糊波面像蒙了一层雾原因网格分辨率不够或者体积分数使用的离散格式耗散太大。一阶迎风下的VOF界面在船艏高压区几乎必然被抹平。解决体积分数离散切换成Geo-Reconstruct或Compressive同时加密自由液面附近的网格。自由液面区域的网格各方向尺寸比控制在1.51以内避免长条形网格让界面被插值拉平。5.4 现象阻力系数波动大时均值无法收敛原因物理模拟时间不够长或者船体运动还处于“过渡期”。瞬态的自由液面流场需要好几轮波在计算域里来回反射才能达到准周期状态尤其是带阻尼区的计算域第一轮反射波特别容易叠加在阻力曲线上。解决延长物理时间至少跑够5个穿越周期或者在数据处理时跳过前20%到30%的过渡段。如果阻力曲线呈现“拍频”现象——两个周期相近的波动叠加——大概率是边界反射波干扰需要加大阻尼区或把侧边界也改成远场条件。5.5 现象网格运动时出现负体积DFBI中途报错原因组件网格和背景网格在重叠区的插值把网格扭曲了或者组件网格边界离背景网格太近运动过程中“顶穿”了背景网格的加密区。解决检查重叠区的最小网格尺寸比组件网格与背景网格尺寸比建议在0.5到2之间把组件网格的尺寸往背景网格靠拢一点或者在DFBI的运动限制器里限制单步最大位移避免一步运动跨过好几个背景网格单元。6. 验证阻力计算结果ITTC公式对比和出口波高检查算出阻力只是第一步验证算得对不对才是真功夫。我最常用两个手段一个是把计算得到的净阻力系数和ITTC 1957公式对比另一个是看出口边界处的波高是否趋近零排除边界反射干扰。ITTC 1957公式给出的是摩擦阻力系数CF 0.075 / (log10(Re) - 2)²其中雷诺数Re UL/ν。以12米船、3.5m/s、水运动黏度1.0×10⁻⁶m²/s为例Re约为4.2×10⁷代入公式算出来CF约为0.00253。船体瞬态仿真算出来的总阻力系数Ct再加上粗糙度补贴和形状因子能跟CF做一个数量级和趋势上的对比——注意量级别指望数值等于CF因为总阻力还包括压差阻力和兴波阻力通常比摩擦阻力大不少。# 阻力结果验证脚本Python示例 Re 4.2e7 cf 0.075 / (math.log10(Re) - 2)**2 ct_sim 0.0038 # Fluent后处理读出的总阻力系数 ratio ct_sim / cf print(fITTC粗糙摩擦系数CF{cf:.5f}, 仿真总阻力系数CT{ct_sim:.5f}, CT/CF{ratio:.2f})这段脚本是拿仿真系数和理论摩擦系数做比值。CT/CF一般在1.2到1.5之间算合理——多出来的部分就是兴波阻力和压差阻力。如果比值小于1.1大概率是自由液面没算好兴波阻力丢失如果大于1.8可能是网格太粗或分离流动没抓住。第二个验证是出口波高。在压力出口处布置一排监测点记录自由液面高度随时间的变化。理想情况下晚期的波高波动幅度应该在来流波高幅值的5%以内。如果波动大说明边界反射波还没衰减干净这时回头加阻尼区或拉长计算域比盲目加密网格更有效。关于出口阻尼区还有一个习惯问题——我每次改动几何或网格后都会强制跑一遍这个波高监测流程哪怕只跑500步物理时间也要确认边界行为没有回弹。从那以后我很少再遇到“换了个网格参数就整体发散”的尴尬。阻力预测验证和边界反射检查建议你第一次做瞬态船体仿真时就养成这个习惯。希望这份拆解对正在用Fluent做船体仿真的人有帮助。照着这套流程把物理模型、网格、时间步长和验证手段都走一遍比你反复试参数要省时得多。本文还有配套的精品资源点击获取