螺旋桨CFD分析
螺旋桨CFD理论基础
螺旋桨性能的基础
老师，如何用CFD预测船舶螺旋桨的性能，要看哪些参数呢？
螺旋桨性能用推力系数 $K_T$、转矩系数 $K_Q$、效率 $\eta_0$ 评估。
其中 $T$ 是推力，$Q$ 是转矩，$n$ 是转速，$D$ 是螺旋桨直径，$J = V_A/(nD)$ 是前进系数。开水面特性（$K_T$-$J$ 曲线和 $K_Q$-$J$ 曲线）代表螺旋桨的基本性能。
CFD计算如何与实验数据比较？
与ITTC（国际船舶试验委员会）或SMPF（船舶性能推进委员会）的标准螺旋桨（如DTMB 4119、KCS）的水池试验数据比较。如果CFD与实验的 $K_T$ 和 $K_Q$ 差异在2%以内则很好，5%以内为实用水平。
数值建模的方法
螺旋桨CFD有哪些建模方法？
主要有三种。(1) MRF（移动参照系）：在旋转坐标系中围绕螺旋桨进行定常计算。成本最低。(2) 滑动网格（Rigid Body Motion）：实际旋转螺旋桨的非定常计算。船体-螺旋桨干涉更精确。(3) 重叠网格（Chimera法）：用背景网格和重叠处理运动。广泛用于STAR-CCM+。
应该选择哪种方法？
开水面特性计算用MRF足够。船尾伴流中的螺旋桨（自航状态）需要滑动网格或重叠网格。空化分析采用非定常滑动网格作为标准。
船舶螺旋桨理论的历史——Rankine-Froude运动量理论（1865年）
螺旋桨推力的最初理论描述是由William Rankine（1865）和W.J.M. Froude（1878）提出的"运动量理论（Actuator Disk Theory）"。螺旋桨被理想化为假想的无限薄圆盘，对流体赋予动量，推导出简单公式 T = ρA(V+va)×2va（va为诱导速度）。这个理论虽然简化了，但抓住了产生推力物理本质，计算简便，至今仍用于设计初期的功率和效率估算。螺旋桨效率理论上限（Betz极限 η=1/(1+va/V)）也由此理论导出，与风力涡轮机的Betz系数（η_max≈59.3%）有相同的数学结构。130年后的CFD能详细建模叶片，进行完整翼型分析，但Rankine理论的清晰思路至今在教学中被传承。
螺旋桨CFD数值计算方法
网格设计
螺旋桨网格生成需要注意什么？
叶片表面附近的边界层网格最关键。用SST k-ω模型时，配置 $y^+ \approx 1$ 的prism层来直接求解壁面。叶前缘和后缘需要局部细化，至少要在前缘曲率半径范围内布置10个单元。为捕捉叶尖涡，叶尖附近网格也需细化。
网格数量的推荐值是多少？
每个叶片50万到200万单元为标准。利用周期边界条件，可以只计算1个叶片，将成本降低到1/$Z$（$Z$为叶数）。但在自航计算中，由于上升流不均匀，必须计算全叶模型。STAR-CCM+的Trimmed Mesh和Fluent的Polyhedral Mesh在品质和自动化上较优。
空化分析
空化的CFD模型是怎样的？
Schnerr-Sauer模型和Zwart-Gerber-Belamri模型为主流。用VOF法求解气液两相流，局部压力低于蒸气压的区域发生气泡生成（蒸发），压力恢复的区域发生消灭（凝聚）。
在Fluent中通过Multiphase > VOF > Schnerr and Sauer激活。STAR-CCM+中在Physics设置Cavitation Model。OpenFOAM用内置Schnerr-Sauer模型的interPhaseChangeFoam。
空化分析的注意事项？
(1) 时间步长极小（$\Delta t \sim 10^{-5}$～$10^{-4}$s）。(2) 网格要能分辨空化层厚度（数毫米）。(3) 在SST k-ω模型上加Reboud修正（密度修正）可改进叶片和云空化的行为。
为什么Schnerr-Sauer模型在空化中受欢迎
在螺旋桨CFD空化分析中，OpenFOAM用户最常用的是Schnerr-Sauer模型。虽然有Zwart-Gerber-Belamri等其他选项，但Schnerr-Sauer的优势在于蒸发和凝聚系数无需调参，只要设定泡核密度就能运行。这种易用性得到广泛支持。但真正再现空化坍塌（内爆）的细节需要LES或超高密度网格，计算规模往往超过10亿单元。现场中"全叶面均匀6层prism层"不崩溃的网格质管理往往是第一关难点。
螺旋桨CFD实务应用
自航分析（Self-Propulsion）
什么是自航分析？
将船体和螺旋桨放在同一计算域中，模拟螺旋桨推进船舶的状态。求出推力与阻力平衡的自航点（Self-Propulsion Point），可预测实船所需功率。
计算步骤如何进行？
(1) 先计算仅船体的阻力（光船阻力）。(2) 再单独计算螺旋桨的开水面特性。(3) 在船体+螺旋桨的滑动网格计算中，调整转速找到推力=阻力的点。STAR-CCM+的Propeller Performance功能可自动搜索自航点。
如何评估船尾伴流的影响？
在螺旋桨面处计算伴流分布（wake fraction）$w = 1 - V_A / V_S$。需评估名义伴流（无螺旋桨）和有效伴流（有螺旋桨）。伴流的非均匀性引起螺旋桨荷载变化、振动和空化，需要详细分析周向分布。
验证用基准算例
有没有可用的螺旋桨CFD验证数据？
有代表性的基准数据如下。
| 螺旋桨 | 叶数 | 特性 | 数据来源 |
|---|---|---|---|
| DTMB 4119 | 3叶 | 非空化验证 | ITTC |
| PPTC VP1304 | 5叶 | 空化/压力脉动 | SVA Potsdam |
| KCS | 5叶 (KP505) | 自航分析。与KCS船体配合 | NMRI/SIMMAN |
这些数据通过SIMMAN研讨会和ITTC基准案例对外发布。
LNG船的省能源螺旋桨——压舱和满载条件兼容设计
液化天然气（LNG）运输船在压舱（空载）和满载状态下吃水差异大，推进效率需求完全不同。过去设计通常优化满载条件，但国际海事组织（IMO）的CII（碳强度指标）法规趋严，两个条件的燃料效率都变得重要。CFD（RANS+自由表面VOF法+螺旋桨转动MRF模型）的多工况优化中，叶片展向桨距分布和弦翼面在两条件间寻找折中。某造船厂的CFD形状优化结果显示，通过这种方法年燃油消耗量降3.8%，同时CO₂排放和CII等级改善，成果在造船学会发表。
螺旋桨CFD软件比较
STAR-CCM+中的螺旋桨分析
听说STAR-CCM+在螺旋桨分析中很强。
STAR-CCM+在海洋工程领域应用广泛。Rigid Body Motion（相当于滑动网格）和Overset Mesh稳定，与DFBI（动态流体结构相互作用）整合优秀。螺旋桨分析推荐配置为Segregated Flow + SST k-ω + VOF（空化时）+ Rigid Body Motion。Polyhedral Mesh + Prism Layer的组合在质量上被推荐。
Ansys Fluent中的螺旋桨分析
Fluent怎么样？
Fluent用滑动网格（Mesh Motion）处理螺旋桨旋转。先用MRF（Frame Motion）计算开水面特性，再用滑动网格进行自航分析的两阶段方法效率高。压力耦合求解器配SST k-ω为标准。空化用Multiphase VOF + Schnerr-Sauer。Fluent Meshing中的polyhedral meshing能生成高质量的螺旋桨叶面网格。
OpenFOAM中的螺旋桨分析
OpenFOAM呢？
MRF用simpleFoam或pimpleFoam加MRFZone配置。滑动网格用pimpleFoam + dynamicMeshDict的solidBodyMotionFvMesh指定旋转。网格用snappyHexMesh生成，但螺旋桨这样复杂形状用商用网格生成器（Pointwise、ANSA等）导入能更稳定保证质量。
OpenFOAM也能做空化分析吗？
interPhaseChangeFoam内置Schnerr-Sauer模型。但与商用求解器相比，空化分析稳定性略弱，对时间步长和网格质量敏感度高。OpenFOAM v2212及以后版本在这方面有改进。
船舶螺旋桨CFD工具比较——MPUF-3A vs StarCCM+ vs OpenFOAM的KT和KQ精度
在螺旋桨推力系数KT和转矩系数KQ预测精度上比较CFD工具，专业海洋码（MPUF-3A等涡升力面法）计算快但对叶尖涡和空化细节再现不足。StarCCM+的MRF和滑动网格组合在船舶领域很成熟，ITTC对比试验中报告KT误差±3～5%，KQ误差±5～8%。OpenFOAM的pimpleFoam + AMI（任意网格接口）经开放验证，在螺旋桨开水面特性基准计算中与商用工具精度相当，公开论文有报道。但在设置方便度和自动化上商用工具依然占优。
螺旋桨CFD前沿研究
水下辐射噪声（URN）
螺旋桨的噪声预测也用CFD做吗？
由于IMO（国际海事组织）水下噪声低减指南的影响，螺旋桨水下辐射噪声（Underwater Radiated Noise）预测变得重要。噪声主要来自空化，空化泡坍塌时产生高压脉冲。
CFD中用Ffowcs Williams-Hawkings（FW-H）方程预测远场声学。Fluent和STAR-CCM+都内置FW-H求解器，可从非定常螺旋桨计算结果得出噪声频谱。但准确预测空化噪声需要网格分辨率和时间分辨率都极高。
能源节减装置（ESD）
为了应对EEDI（能源效率设计指标），螺旋桨周围的装置受关注。
螺旋桨前的Pre-Swirl Device（定子、导管）和后的Post-Swirl Device（PBCF：螺旋桨中心盖翅、Hubvane）是代表。这些装置通过优化螺旋桨周围流场，推进效率可提高2～6%。
ESD的CFD评估方法？
分别进行装置装着和未装着的自航计算，量化所需功率的差异。对ESD形状优化用参数化研究或伴随法。STAR-CCM+的Design Manager和Fluent的Adjoint Solver可探索ESD翼角和弦长的最优值。
电动推进和轮毂驱动
电动推进的螺旋桨设计也用CFD吗？
用的。方向推进器（ABB的Azipod、Rolls-Royce的Mermaid等）和轮毂驱动（rim-driven thruster）的CFD中，电动机部分的间隙流和电磁力耦合成为追加课题。吊舱（nacelle）和支撑杆的形状优化也属CFD范畴。
螺旋桨空化侵蚀预测——CFD与侵蚀模型耦合的最前沿
船舶螺旋桨的空化侵蚀是气泡坍塌时产生的局部冲击压力（数GPa）反复敲击金属表面造成。CFD（RANS平均+空化模型）求得气泡生成和坍塌分布，再转换为侵蚀风险指标（Erosion Risk Indicator: ERI）的耦合方法是前沿。特别是"微喷流气泡坍塌（Microjet）"的压力方向性（垂直壁面还是斜向）决定侵蚀强度，单纯空化体积分数无法准确预测侵蚀部位。DNV GL公开的基准试验中，CFD空化模型的侵蚀预测与实验侵蚀迹线空间一致需要历经5代改良。
螺旋桨CFD故障排除
$K_T$ 与实验值偏差超过10%
开水面计算中 $K_T$ 与实验值不符。
逐项检查。(1) 叶片表面网格是否充分分辨前缘曲率（前缘处至少10个单元）。(2) 是否达到 $y^+ < 1$。特别SST k-ω模型中 $y^+$ 为5～30的中间值会精度恶化。(3) 叶尖后流是否充分求解。进行网格敏感性分析确认收敛。
MRF和滑动网格的结果不一样？
在一致流入的开水面条件下，MRF和滑动网格的定常平均值应基本相同。差异超5%则可能是MRF接口网格不匹配。Fluent中检查Mesh InterfaceInterpolation设为Conservative，确认没有orphan cells。
滑动网格压力振荡
滑动网格计算中压力大幅振荡，接近发散。
时间步长过大。螺旋桨每步转过的角度应不超1度为目标。举例：$n = 10$ rps时，$\Delta t < 1/(360 \times 10) \approx 2.8 \times 10^{-4}$s。还有滑动接口处旋转域和固定域的网格尺寸差异过大会导致插值误差。接口面的单元尺寸应一致。
空化分析不稳定
空化分析数转后发散。
加空化模型在数值上易不稳。(1) 先用非空化条件充分收敛流场作初值。(2) 进一步减小时间步长（$\Delta t \sim 10^{-5}$s）。(3) 降低VOF的under-relaxation至0.3左右。(4) 科学数控制在0.5以下。STAR-CCM+的自动CFL Ramping功能能自动克服初期不稳定。
螺旋桨CFD推力与实验差15%——叶轮模型设置和旋转区域边界处理
螺旋桨CFD推力与实验偏差超15%的典型原因有两个。第一是"MRF（移动参照系）区域设置错误"——旋转区域边界离叶尖太近会在旋转到静止坐标系的速度变换中放大误差。MRF边界推荐离螺旋桨半径R至少1.5倍以上。第二是"滑流速度处理不当"——螺旋桨后流强旋转（涡旋成分）传播到出口边界时，若出口距离不足3～5D会产生反流，推力计算变不稳定。有自由表面时，水面换气（Air Ingestion）可能导致螺旋桨转速大幅波动，需与VOF耦合。