MarineFlow 使用路线：从船体模型到流场和水动力结果 / A MarineFlow Workflow from Hull Geometry to Hydrodynamic Results
📅 创建时间：2026-07-20 🏷️ 标签：#MarineFlow #CFD工作流 #船体阻力 #螺旋桨 📚 前置知识：[[04-cfd-equations]]
1. MarineFlow 能力边界
根据安装目录，可以确认 MarineFlow 集成了几何、网格、OpenFOAM、NaViiX、MPI、CGNS/HDF5 和 VTK 类工具。本篇描述的是可迁移的标准使用方法，具体菜单名称以实际版本为准。
2. 创建项目之前
先准备算例说明：
text
目标：计算什么？
模型：船体、螺旋桨还是组合体？
工况：速度、转速、吃水、流体属性？
物理：稳态/瞬态、自由液面、空化、湍流？
输出：阻力、推力、扭矩、压力还是波面？
验证：和试验、经验公式或公开算例怎样比较？3. 几何导入与检查
检查内容
- 船体尺度和单位
- 坐标原点和航向
- 水线位置
- 表面法向
- 是否存在破面和缝隙
- 是否有重复面
- 螺旋桨轴线和旋转方向
流体几何的特殊要求
结构有限元可以接受开放曲面，但生成封闭流体域通常需要水密表面。一个很小的缝隙也可能让体网格生成失败。
4. 建立计算域
以船体阻力为例：
text
顶部空气边界
┌─────────────────────┐
入口 → │ 自由液面 │ → 出口
│ 船体 │
│ │
└─────────────────────┘
底部和侧边界边界太近可能：
- 限制流线
- 改变压力
- 反射波浪
- 夸大阻塞效应
可以通过扩大计算域进行敏感性检查。
5. 网格生成
表面网格
要准确描述：
- 船首曲率
- 船尾尖锐区域
- 螺旋桨前缘和尾缘
- 水线附近几何
体网格
需要在关键区域加密：
- 近壁面
- 自由液面
- 船尾尾流
- 螺旋桨旋转域
- 预计空化区域
边界层网格
壁面函数计算通常通过 y+ 判断首层网格是否合适。不同湍流处理要求不同 y+ 范围，不能只追求越小越好。
6. 定义边界
常见边界：
- inlet：来流速度或质量流量
- outlet：压力或零梯度
- hull：无滑移壁面
- symmetry：对称面
- atmosphere：大气压力
- seabed/farfield：底部和远场
- rotating interface：旋转区域接口
边界命名应在网格和求解配置之间保持稳定。
7. 选择物理模型
阻力快速评估
可以从稳态、单相或简化自由液面开始，先获得初始流场和基本趋势。
兴波阻力
需要空气-水两相和重力，通常使用 VOF。
螺旋桨平均性能
可以从旋转参考系方法开始，成本低于真实滑移网格。
叶片周期载荷和压力脉动
需要瞬态滑移网格或等价非定常方法。
空化
需要相变模型、蒸汽压、足够时间分辨率和局部网格。
8. 设置求解控制
需要理解而不是照抄的参数：
- 时间步或伪时间步
- 最大 Courant 数
- 离散格式
- 压力求解器
- 速度求解器
- 预条件器
- 欠松弛因子
- 终止条件
- 输出时间间隔
输出过于频繁会显著增加磁盘占用和并行 I/O 时间。
9. 初始化
好的初始场可以显著减少迭代：
- 先使用势流或简化稳态结果
- 先计算低速工况，再逐步提高速度
- 先关闭空化，获得稳定流场后再开启
- 先使用一阶格式，稳定后切换二阶
这些是延续法思想：先求容易问题，再逐步接近目标问题。
10. 并行计算
典型过程：
text
检查完整网格
→ 按区域分解
→ 启动多个 MPI 进程
→ 每个进程求解自己的子域
→ 交换交界面数据
→ 计算结束后重构结果进程数选择要考虑：
- 每个分区是否有足够单元
- 分区边界是否过多
- 机器内存
- CPU 核心数
- 磁盘和网络能力
小算例使用过多进程可能更慢。
11. 计算监控
数值指标
- 各方程残差
- 连续性误差
- Courant 数
- 线性求解器迭代次数
工程指标
- 阻力
- 推力
- 扭矩
- 自由液面波高
- 关键点压力
- 空化体积
残差下降但工程监控量仍在漂移，通常不能认为结果已经稳定。
12. 后处理
船体阻力
- 检查力的坐标方向
- 区分压力和黏性贡献
- 转换为无量纲系数
- 做时间平均
螺旋桨
- 推力、扭矩和效率
- 叶面压力
- 叶尖涡
- 尾流速度分布
- 单叶周期载荷
自由液面
α=0.5等值面- 波高切线
- 波浪传播和反射
空化
- 汽相体积分数等值面
- 空泡随时间的演化
- 叶面低压区域
13. 一个结果是否可用
至少完成：
text
□ 网格质量检查
□ 质量守恒检查
□ 工程监控量稳定
□ 网格敏感性分析
□ 时间步敏感性分析（瞬态）
□ 计算域敏感性分析
□ 与试验或公开数据对比下一篇：[[06-structural-physics]]