螺旋桨空蚀FSI
螺旋桨空蚀FSI理论基础
空蚀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倍多。这种在民用船舶中必须避免的现象，在武器应用中却成了制胜的秘诀——这种逆向思维正是空蚀理论的精妙之处。
螺旋桨空蚀FSI数值计算方法
旋转体FSI处理
旋转螺旋桨的FSI如何处理？
使用滑移面（滑动网格/旋转参考系）来设置旋转流体区域。FSI耦合中，螺旋桨叶面是界面，在旋转坐标系中执行结构分析。
在Ansys Fluent中，通常先用MRF（移动参考系）进行定常计算，然后切换到滑动网格进行非定常计算。要捕捉空蚀的非定常特性，必须使用滑动网格。
空蚀模型设置
空蚀模型参数如何设置？
对于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%。
螺旋桨空蚀FSI实际应用
分析步骤
螺旋桨空蚀FSI分析的实际步骤是什么？
1. 获取螺旋桨叶片形状的CAD数据
2. 生成流体网格（旋转区域+静止区域）
3. 创建结构有限元模型（仅叶片；毂盘为刚体）
4. 在无空蚀湿润条件下验证FSI功能
5. 启用空蚀模型进行主计算
6. 评估推力系数、扭矩系数和空蚀模式
网格需要多细？
叶面的$y^+$应小于1。特别是叶尖和前缘周围易发生空蚀，应加密网格。
| 区域 | 单元尺寸指南 |
|---|---|
| 叶面边界层 | $y^+ < 1$，15层及以上 |
| 翼尖涡流域 | 弦长/100 |
| 尾流域 | 弦长/50 |
| 远场 | 弦长/5 |
有哪些基准问题用于验证？
PPTC螺旋桨（波茨坦螺旋桨测试用例）是ITTC（国际水池会议）的标准基准。空蚀模式、推力系数$K_T$和扭矩系数$K_Q$的实验值都已公开。
大型油轮的螺旋桨更换——现场"声音"诊断
资深的船舶工程师常说"通过声音可以识别螺旋桨空蚀"。正常情况下转速声低沉温和，但出现空蚀时，船尾会发出"吧嗒吧嗒"或"卡嚓卡嚓"的异音。实际上，排水量30万吨级的超大油轮（VLCC）螺旋桨直径超过9米，转速仅每分钟90转左右。即使这样，空蚀发生时3个月内就能在叶片表面形成1厘米以上的侵蚀孔。若能通过仿真预测侵蚀量，优化干坞期，每年可节省数千万元成本。
螺旋桨空蚀FSI软件比较
工具比较
螺旋桨空蚀FSI有哪些可用软件？
STAR-CCM+在海事领域人气为什么那么高？
西门子能源（前CD-adapco）在海事CFD领域历史悠久。STAR-CCM+的重叠网格和自动网格生成在螺旋桨分析中易用。同时支持与Abaqus的强耦合FSI。
复合材料螺旋桨怎么处理？
需要能建模异向铺层结构的结构求解器。Abaqus的铺层壳元（S8R等）和Nastran的PCOMP/PCOMPG定义都很适合。Ansys Mechanical ACP（复合材料前后处理）也便于铺层定义。
造船大国的分析工具——韩国、中国、日本的应用差异
争夺全球造船份额的韩国现代重工和中国中船集团，在螺旋桨空蚀分析上用什么工具？韩国趋向于STAR-CCM+和ANSYS Fluent并行，积累了与水池试验相关的海量数据。中国对OpenFOAM进行改进，开发了独有的分析栈。日本三菱、川崎、今治则倾向于采用FINE/Marine等船舶专用求解器的高可信度。有趣的是，即使采用相同空蚀模型（Schnerr-Sauer法），因网格策略和时间步不同，结果可能相差10～20%，各公司都守着这些诀窍。
螺旋桨空蚀FSI前沿研究
复合材料螺旋桨自适应变形
为什么复合材料螺旋桨备受关注？
CFRP（碳纤维增强塑料）的各向异性特性可实现荷载相关的叶片桨距角变化——弯扭耦合（BTC），在宽运行范围内维持高效率的自适应螺旋桨成为可能。
铺层方向与BTC系数的关系由式表示
当$D_{16}$非零时（如[+45/-45]非对称铺层），会产生BTC效应。
空蚀侵蚀预测
能通过仿真预测空蚀对叶面的侵蚀吗？
可以估算气泡崩溃时的冲击压力，与材料耐侵蚀性比较。Patela（2016）的侵蚀强度指标方法正逐步推广。通过FSI，叶片变形改变空蚀模式，进而影响侵蚀分布——这种间接耦合也很重要。
噪声预测的集成
空蚀噪声预测与FSI有关吗？
空蚀是船舶水下辐射噪声的主要声源。应用FW-H声学类比到CFD结果以预测远场噪声。叶片变形改变空蚀体积波动，影响噪声频谱。Ansys Fluent和STAR-CCM+都内置FW-H求解器。
量子计算与空蚀——下一代仿真的前沿
空蚀-FSI耦合的计算成本异常高昂。忠实追踪单个气泡崩溃的DNS计算，即使用现代超算也需数周才能计算螺旋桨一转。因此量子计算加速格子Boltzmann法备受关注。IBM与DLR的联合研究中，用量子线路表示流体的相空间，理论上可实现相比经典计算100倍的加速（2023年）。虽然实用化还远，但"用量子比特求解Navier-Stokes方程"这一理念本身正在根本改变FSI分析的视野。
螺旋桨空蚀FSI故障处理
空蚀未出现
计算中未能再现实验观测到的空蚀。
整理核查清单。
| 确认项 | 详情 |
|---|---|
| 蒸汽压设置 | 是否为对应温度的正确$p_v$ |
| 参考压力 | operating pressure设置是否正确 |
| 乱流模型 | 推荐SST k-omega。k-epsilon过大估计粘度，抑制空蚀 |
| 网格密度 | 叶面低压区是否充分解析 |
| 时间步长 | 定常计算中空蚀无法稳定 |
推力和扭矩不匹配
$K_T$或$K_Q$与实验值偏差很大。
FSI耦合不收敛如何处理？
旋转体FSI中，离心力引起的初始变形较大，应先对结构单独进行离心力分析，以变形后的形状作为FSI计算的初始形状。同时将耦合弛豫系数设小（0.2～0.5），稳定后逐步增大。
"分析无问题"却船体振动——常见陷阱
实务中有个著名失败案例："空蚀-FSI分析显示船体无振动问题"，但投入使用后船尾甲板激烈震动，船员不适。原因查明后发现：分析假设了均匀来流，但实际上螺旋桨前方有船体，产生不均匀伴流，流速周向变化达30%。这种不均匀性导致空蚀发生时刻错开，周期压力脉动被放大。"均匀来流假设"简化了计算，但在振动问题上尤需留意。
价值
详细
错误