螺旋桨空化流固耦合
理论与物理
空化FSI概述
为什么螺旋桨的空化需要FSI？
空化（空穴现象）会在螺旋桨叶片表面产生非定常压力波动，引发叶片振动、侵蚀和噪声。叶片的弹性变形会改变空化模式，因此需要流体-结构耦合分析。
控制方程
空化的数学模型是怎样的？
广泛使用均质混合流模型。使用液相和气相的混合密度求解Navier-Stokes方程。空化数定义为：
$p_v$ 是蒸汽压。相变基于Rayleigh-Plesset方程的质量传输模型描述。Schnerr-Sauer模型和Zwart-Gerber-Belamri模型是代表性模型。
$\alpha$ 是气相体积分数，$\dot{m}^+, \dot{m}^-$ 是蒸发和冷凝的质量转移率。
结构侧如何建模？
螺旋桨叶片使用壳单元或实体单元建模。对于复合材料螺旋桨（如CFRP），需要考虑铺层结构的各向异性材料模型。FSI界面上的数据传递在湿表面上进行。
空化数σ——一个无量纲数决定命运
螺旋桨设计者首先计算的就是空化数σ＝(p∞－pv)/(½ρU²)。σ低于1左右时开始产生气泡，跌破0.3则进入螺旋桨完全被蒸汽包裹的“超空化”状态。有趣的是，俄罗斯的“暴风”鱼雷特意利用超空化实现了水下340公里/小时的速度，是普通鱼雷的4倍以上。民用船舶中极力避免的现象，换个思路却能成为武器的秘诀——这种逆向思维也正是空化理论的魅力所在。
各项的物理意义
- 结构-热耦合项：温度变化引起的热膨胀诱发结构变形，变形又影响温度场。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常例子】夏天铁轨伸长导致缝隙变小——温度升高→热膨胀→产生应力的典型例子。电子电路板焊接后翘曲也是不同材料热膨胀系数差异导致的。发动机气缸体因高温区和低温区的温差产生热应力，最坏情况下会导致裂纹。
- 流体-结构耦合（FSI）项：流体压力、剪切力使结构变形，结构变形又改变流体区域，是双向相互作用。【日常例子】强风下悬索桥缆索振动（涡激振动）——风力使结构摇晃，摇晃的结构改变风流，进而放大振动。心脏血流与血管壁的弹性变形、飞机机翼颤振（气动弹性不稳定性）也是典型的FSI问题。有时单向耦合即可满足，但变形较大时双向耦合是必须的。
- 电磁-热耦合项：焦耳发热 $Q = J^2/\sigma$ 引起温度升高，温度变化又改变电阻，形成反馈回路。【日常例子】电炉的镍铬丝通电发热（焦耳热）变红——温度升高电阻改变，电流分布也随之变化。IH电磁炉的涡流发热、输电线路温度升高导致垂度增加也是此类耦合的例子。
- 数据传递项：通过插值解决不同物理场间网格不匹配的问题。【日常例子】天气预报中结合“气温数据”和“风速数据”计算体感温度时，若观测点不同就需要插值——CAE耦合分析中，结构网格和CFD网格通常也不一致，因此界面上的数据传递（插值）精度直接关系到结果的可靠性。
假设条件与适用范围
- 弱耦合假设（单向耦合）：当一方物理场影响另一方而反向影响可忽略时有效
- 需要强耦合的情况：FSI中的大变形、电磁-热耦合中温度依赖性较强的情况
- 时间尺度分离：各物理场特征时间差异较大时，可采用子循环提高效率
- 界面条件一致性：需确认耦合界面上的能量和动量守恒在数值上得到满足
- 不适用的情形：三个及以上物理场同时强耦合时，有时可能需要整体式方法
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 热膨胀系数 $\alpha$ | 1/K | 钢：约12×10⁻⁶，铝：约23×10⁻⁶ |
| 耦合界面力 | N/m²（压力）或N（集中力） | 确认流体侧与结构侧的力平衡 |
| 数据传递误差 | 无量纲（%） | 插值精度取决于网格密度比。建议控制在5%以下 |
数值解法与实现
旋转体的FSI处理
如何处理旋转螺旋桨的FSI？
使用滑移面（sliding mesh / rotating reference frame）设置旋转流体区域。FSI耦合中将螺旋桨叶片表面作为界面，在旋转坐标系下进行结构分析。
Ansys Fluent中的标准做法是先用MRF（Moving Reference Frame）进行稳态计算，再切换到sliding mesh进行非稳态计算。要捕捉空化的非定常特性，必须使用sliding mesh。
空化模型设置
如何设置空化模型的参数？
对于Schnerr-Sauer模型，气泡数密度 $n_b$ 是关键参数。默认值为 $n_b = 10^{13}$ /m³，但有时需要根据实验对比进行调整。
| 参数 | Schnerr-Sauer | Zwart-Gerber-Belamri |
|---|---|---|
| 核密度/核半径 | $n_b = 10^{13}$ /m³ | $R_b = 10^{-6}$ m |
| 蒸发系数 | 1.0 | $F_{vap} = 50$ |
| 冷凝系数 | 1.0 | $F_{cond} = 0.01$ |
| 饱和蒸汽压 | 温度相关 | 温度相关 |
时间步长如何设置？
推荐每旋转一周360~720个时间步（0.5°~1.0°/步）。空化溃灭过程极快，局部区域最好满足CFL < 1。
BEM与CFD——螺旋桨分析的“双刀流”运用
螺旋桨流体分析长期以来并行使用边界元法（BEM）和CFD两种方法。BEM基于势流理论，计算速度极快，适合设计初期阶段成百上千个工况的快速筛选。另一方面，CFD的LES（大涡模拟）能捕捉到空化气泡的生成和溃灭，但每个工况需要数百CPU小时。实际船舶设计采用“BEM粗筛→CFD精算→水池实验验证”的三段式标准流程。最近，将BEM结果作为CFD初始条件的混合方法也在增多，已有案例能将总分析时间减少30~40%。
整体式方法
将所有物理场作为一个联立方程组同时求解。对强耦合问题稳定，但实现复杂，内存消耗大。
分区法（分离迭代法）
各物理场独立求解，在界面交换数据。易于实现，可利用现有求解器。适用于弱耦合。
界面数据传递
最近邻法（最简单但精度低）、投影法（具有守恒性）、RBF插值（对网格不一致鲁棒性强）。需权衡守恒性和精度。
子迭代
在每个耦合步内进行充分迭代，确保界面条件的一致性。残差基准基于各物理场的典型值进行缩放。
Aitken松弛
自动调整耦合迭代的松弛因子。防止过度松弛导致发散，是一种能加速收敛的自适应方法。
稳定性条件
注意附加质量效应（流体-结构耦合中结构密度≈流体密度时）。不稳定时可应用Robin型界面条件或IQN-ILS法。
Aitken松弛的比喻
Aitken松弛类似于“平衡跷跷板”。一方推得太用力，另一方就会弹起，反弹又导致推得更用力——为了抑制这种振荡，自动调整推力大小的就是Aitken松弛。当耦合迭代振荡不收敛时，它会根据前一次的修正量自动调整下一次的修正量，是一种自适应方法。
实践指南
分析步骤
螺旋桨空化FSI分析的实际工作步骤是？
1. 获取螺旋桨叶片形状CAD数据
2. 生成流体网格（旋转区域+静止区域）
3. 创建结构有限元模型（仅叶片。轮毂视为刚体）
4. 在湿表面无空化条件下验证FSI运行
5. 启用空化模型进行正式计算
6. 评估推力系数、扭矩系数、空化模式
网格需要细化到什么程度？
叶片表面 $y^+$ 最好小于1。特别是叶尖区域和前缘附近是空化起始点，需要加密网格。
| 区域 | 单元尺寸参考 |
|---|---|
| 叶片表面边界层 | $y^+ < 1$，15层以上 |
| 叶尖涡区域 | 弦长/100 |
| 尾流区域 | 弦长/50 |
| 远场 | 弦长/5 |
验证有哪些基准问题？
PPTC螺旋桨（Potsdam Propeller Test Case）是ITTC（国际拖曳水池会议）的标准基准。其空化模式、推力系数 $K_T$、扭矩系数 $K_Q$ 的实测数据已公开。
大型油轮的螺旋桨更换——现场“听音”诊断
经验丰富的船舶工程师常说“螺旋桨空化听声音就知道”。正常时是低沉平缓的旋转声，一旦发生空化，船尾就会传来“噼啪”、“嘎吱”的异响。实际上，对于30万吨级的VLCC（超大型油轮），螺旋桨直径超过9米，转速仅每分钟90转左右。即便如此，若发生空化，三个月内叶片表面就可能出现1厘米以上的侵蚀孔。如果能通过仿真预测侵蚀量，优化进坞维修时间，每年可节省数千万日元的成本。
分析流程的比喻
你吹过气球吗？那个瞬间，其实发生了高度复杂的流体-结构耦合。内部空气压力（流体）推展橡胶壁（结构）→扩张的壁改变了内部压力分布→改变的压力进一步使壁变形…在计算中逐步重复这个“抛接球”过程的就是FSI分析。
初学者易犯的错误
“单向耦合足够了吧？”——这个判断失误在耦合分析中最危险。结构变形微小的话，单向耦合确实足够。但对于像心脏瓣膜开闭这样变形会大幅改变流路的情况，单向耦合完全不行。判断标准是“变形量是否超过特征长度的1%”。超过就必须使用双向耦合。如果错误地使用单向耦合，结果会“看似合理实则大错特错”——这是最可怕的情况。
边界条件的思考方式
耦合界面的数据交换就像“边境的出入境管理”。各国（物理场）有自己的法律（控制方程），但如果在边境（界面）不能准确管理人员和物资（力、温度、位移）的往来，两国的经济（能量平衡）就会崩溃。网格不一致时的插值就像“翻译”——误译（插值误差）越小，结果越好。
软件比较
工具比较
可用于螺旋桨空化FSI的软件有哪些？
なった
詳しく
報告