简介这份专题资料面向水动力学、船舶工程与流体力学方向的研究人员和学生围绕NACA0015二维水翼的空化现象及其控制机理展开整理可帮助读者理解空化数值计算的建模思路与参数选取依据。内容涵盖计算域设置弦长c0.1m攻角6°~10°进口距前缘5倍弦长、出口距尾缘10倍弦长、上下边界相距12倍弦长、空化数定义与取值空化数0.9、进口速度12m/s、汽化压强0.044bar、时间步长10⁻⁵s并延伸至Clark-Y翼型弦长70mm、三维展宽68mm及700mm×190mm矩形二维计算域的替代方案还涉及NACA66(MOD)的f/c0.020、a0.8、t/c0.9等翼型参数与汽化压强、空化数关系。资源为单个docx文档压缩包约332KB结构紧凑便于检索查阅已有112人学习。适合需要参考翼型空化建模、工况设定与气泡动力学分析的读者作为专题笔记使用。1. 从一份 NACA0015 空化计算笔记说起算例能跑通靠的是边界条件不是网格数很多人第一次做水翼空化仿真卡住的地方不是网格画不细而是算出来的空化数对不上、云空化脱落周期跟文献差一个量级。这份 2021-2022 年的专题笔记里NACA0015 二维水翼算例其实已经把关键量给全了弦长 c0.1m攻角 6°~10°弧形进口边距前缘 5c出口边距尾缘 10c上下边界相距 12c翼型居中空化数取 0.9进口速度 12m/s汽化压强 0.044bar时间步长 1e-5 s。这些数字看着琐碎但每一个都直接决定空化数是否算得准、回流水汽是否混得进去。这份材料适合两类人一类是做船舶推进、水下航行器、水泵诱导轮空化分析的研究生和工程师另一类是已经会点 CFD 但总在空化数守恒和云空化非定常捕捉上翻车的人。Clark-Y 翼型弦长 70mm、三维展向 68mm、NACA66(MOD) 的 f/c0.020、a0.8、t/c0.9 这些参数也在笔记里出现说明作者是在做多翼型对比。2. NACA0015 二维计算域搭建与空化数标定2.1 计算域几何为什么按 5c/10c/12c 布笔记里的域布置不是随手写的。进口放到前缘上游 5 倍弦长是为了让来流在到达翼型前有足够距离发展成均匀流避免把入口的数值扰动直接怼到前缘出口放到尾缘下游 10 倍弦长是给尾迹和脱落的空泡云留出耗散空间否则出口边界会把未衰减的涡结构反射回计算域云空化脱落频率就假了。上下边界相距 12 倍弦长、翼型居中得到 6c 的半无限感阻塞比对 NACA0015 在 6°~10° 攻角下阻塞效应可以忽略。边界位置类型备注进口前缘上游 5c速度入口v12 m/s出口尾缘下游 10c压力出口由空化数反算上/下翼型上下各 6c滑移/对称视求解器而定翼型弦长 c0.1m无滑移壁面y 需匹配湍流模型弧形进口边是笔记里特意提的细节常见做法是把入口段做成以翼型前缘为圆心的圆弧使来流方向沿径向指向前缘攻角变化时不用重画域只改入口速度矢量方向即可做 6°、9° 多攻角扫描时省事。2.2 空化数定义与出口压力反算空化数是这一整套算例的锚点笔记给的表达式是$$\sigma \frac{p_\infty - p_v}{\tfrac{1}{2}\rho v_\infty^2}$$其中 v∞ 是进口速度ρ 是水密度pv 是汽化压强笔记取 0.044barp∞ 是参考压强工程上常用出口压力代替。给定 σ0.9、v∞12 m/s就能反算出口压力。水密度取 998.2 kg/m³pv0.044bar4400 Parho 998.2 # 水密度 kg/m^3 v_inf 12.0 # 进口速度 m/s p_v 4400.0 # 汽化压强 Pa (0.044 bar) sigma 0.9 # 目标空化数 q 0.5 * rho * v_inf**2 # 动压 p_out p_v sigma * q # 由 sigma 定义反算出口压力 print(f动压 q {q:.1f} Pa) print(f出口压力 p_out {p_out:.1f} Pa {p_out/1e5:.4f} bar)逻辑说明空化数定义里的 p∞ 用的是远场参考压力二维外流里最方便的就是把压力出口设成这个值进口只给速度。参数说明σ 越小越容易空化0.9 属于中等空化强度能稳定出现片状空化并周期性脱落把 σ 从 0.9 调到 1.5 左右通常片空化会缩短甚至消失调到 0.5 以下容易发展成超空化。这一步算错了后面看到的所有云空化图像都不可信。提示pv 强烈依赖水温。笔记给的 0.044bar 对应约 30℃ 的水20℃ 时 pv 约 2340 Pa。如果仿真里水温设错空化数会系统性偏移脱落频率跟着一起偏。2.3 时间步长与库朗数约束笔记给的时间步长是 1e-5 s。这个值不是拍脑袋定的它要满足库朗数 Co v·Δt/Δx ≤ 1最好控制在 0.5 附近。进口 12 m/s若前缘附近最小网格 0.5mm则 Co 12×1e-5/5e-4 0.24属于比较稳的范围。若把网格加密到 0.1mm 而时间步不变Co 会跳到 1.2非定常空化脱落就会数值发散或频率失真。# 以 OpenFOAM interCondensatingEvaporatingFoam 为例检查时间步与库朗数 # controlDict 关键项 deltaT 1e-5; # 笔记给定时间步长 adjustTimeStep no; # 空化非定常算例一般不自动调步 maxCo 0.5; # 若开启 adjustTimeStep 时的目标库朗数 writeInterval 0.01; # 脱落周期约几十 ms取足够密参数说明adjustTimeStep 关掉是为了让脱落的频谱分析有均匀采样自动变步会让 FFT 出现频率涂抹。writeInterval 建议让一个脱落周期内至少存 50 帧否则 PSD 峰位置会糊。3. 求解器选型与空化模型参数落地3.1 均质平衡流模型与输运方程模型怎么选水翼云空化是典型的两相可压缩问题主流两条路一是均质平衡流正压状态方程二是输运方程类空化模型Kunz、Schnerr-Sauer、Zwart-Gerber-Belamri。均质平衡流把混合物密度写成压力的函数算得快、对片空化稳定但蒸发冷凝的滞后效应表达不了云空化脱落的细节偏弱。输运方程模型显式求解汽相体积分数输运能捕捉脱落和回射流代价是要调蒸发/冷凝系数。笔记的场景是研究控制机理需要看脱落和回射流所以我一般选 Schnerr-Sauer 或 Zwart。Schnerr-Sauer 里的成核点密度参数很关键# Schnerr-Sauer 关键参数constant/phaseProperties 或 cavitationProperties Cv 1.0; // 冷凝系数 Cc 1.0; // 蒸发系数 n 1.0e13; // 单位体积成核点数量 1/m^3 dNuc 1.0e-6; // 成核点初始直径 m参数说明n 和 dNuc 控制初始汽核含量n 调大空化更容易触发、片空化更靠前Cv、Cc 分别控制冷凝和蒸发速率抖得厉害时把 Cc 降到 0.5 往往能压住。Zwart 模型对应的 Cvap、Ccond 通常在 50 和 0.01 量级起步再调。模型蒸发/冷凝系数适用云空化脱落Schnerr-SauerCc/Cv ~1.0片空化、中高 σ可捕捉ZwartCvap50, Ccond0.01通用可捕捉Kunz率常数调稳定计算稍弱均质平衡流无快速扫参偏弱3.2 湍流模型与壁面处理空化区域是强逆压梯度加相变SST k-ω 是这类算例的常规选择因为它对逆压梯度和壁面附近处理比标准 k-ε 好。第一层网格高度按 y 定SST 用壁面函数时 y 落在 30~100 即可若要做低雷诺数解析则 y1网格量会涨一个量级。二维 NACA0015 弦长 100mm我把第一层高度定在 0.05mm 左右对应 y≈40既能抓住前缘吸力面的片空化又不至于让网格在三万以内下不来。# snappyHexMesh 局部加密片段前缘与吸力面加密 refinementRegions { frontRegion { mode distance; levels ((0.02 4) (0.005 6)); // 前缘 2cm 内 4 级, 5mm 内 6 级 } suctionSide { mode distance; levels ((0.01 5)); // 吸力面周围 1cm 内 5 级 } }逻辑说明空化主要在吸力面前半段和尾缘附近发生全局加密性价比低用距离场在翼型周围做局部加密能在网格量可控的前提下把空化界面捕捉住。参数说明levels 里前一个数是距离阈值单位 m后一个数是加密级数每加一级网格尺度减半。前缘最外 2cm 给 4 级、5mm 内给 6 级是二维算例比较省的配法。注意空化区汽液界面处的网格分辨率直接决定脱落形态。若只在前缘加密而不在尾缘加密脱落的空泡云会过早耗散看起来像不脱落。4. Clark-Y、NACA66(MOD) 多翼型算例与工况对比4.1 Clark-Y 弦长 70mm、三维展向 68mm 的域怎么摆笔记里 Clark-Y 的设定是弦长 70mm、三维展宽 68mm弦展比接近 1属于低展弦比矩形翼。二维域可以照 NACA0015 的方法按 5c/10c 缩放也可以用笔记给的矩形域长 700mm、高 190mm翼型前端距进口 130mm、高度居中、攻角 0顺流平放。这个域的阻塞比是 70/190≈0.37比 NACA0015 的域大不少做零攻角基准工况够用但做 9° 攻角时上边界离吸力面太近阻塞会抬高空化数阈值。参数NACA0015Clark-Y弦长 c0.1 m70 mm域长×高—700×190 mm前缘距进口5c130 mm攻角6°~10°0°平放展向二维68 mm三维三维计算要注意展向端面。若不做周期性或对称处理端面回流会从翼尖附近倒灌进来污染空化区。常见做法是端面设成对称面或自由滑移再把展向取到足够长让端面影响远离主流空化区。4.2 NACA66(MOD) 参数含义与建模注意笔记给的 NACA66(MOD) 参数是 f/c0.020、a0.8、t/c0.9。这里要小心三个量的常规含义分别是最大弯度比 f/c、中弧线相对位置 a、厚度比 t/c。但 t/c0.9 对水翼来说是极厚的不太可能一般水翼 t/c 在 0.05~0.15怀疑笔记里 t/c 记的是别的量或笔误。实际建模前应该回到原始翼型坐标库核对NACA66 mod 属于 NACA 6 系列修正型常用中弧线 a0.8 指中弧线在 80% 弦长附近达到最大弯度。# NACA 6 系列中弧线近似a0.8, f/c0.020示例 import numpy as np def naca6_meanline(xc, a0.8, fc0.020): # xc: 弦向坐标 0~1 yc np.zeros_like(xc) # a0.8 时最大弯度位置约在 xca for i, x in enumerate(xc): if x a: yc[i] fc / a**2 * (2*a*x - x**2) else: yc[i] fc / (1-a)**2 * ((1-2*a) 2*a*x - x**2) return yc x np.linspace(0, 1, 101)逻辑说明这段是给建模脚本用的中弧线生成把 f/c 和 a 代进去就能得到弯度分布再叠加厚度分布得到完整翼型。参数说明a 越大最大弯度越靠后前缘载荷越低、空化起始越靠后跟 NACA0015 对比正好能看出弯度位置对空化起始的影响。所有坐标生成完建议先画出来跟翼型库叠一下确认没跑偏再进网格。4.3 多工况扫描的批处理思路做 0°、6°、9° 攻角加多个空化数的矩阵扫描一算就是十几个算例。手动改文件不现实常规做法是脚本化生成 casefor aoa in 0 3 6 9; do for sig in 0.5 0.9 1.5; do case${aoa}deg_sigma${sig} cp -r template_case $case # 改入口速度方向为攻角 sed -i s/ALPHA/$aoa/ $case/0/U # 改出口压力为对应空化数反算值 python set_sigma.py --case $case --sigma $sig done done参数说明aoa 控制入口速度矢量方向注意改的是矢量方向而不是网格网格在攻角扫描时保持不变最省事sig 通过 set_sigma.py 按 2.2 的公式反算出口压力写进 p 文件。每个 case 跑完统一后处理把吸力面压力系数分布和汽相体积分数时程抽出来对比才能看出攻角和空化数哪个对脱落频率影响更大。提示多工况扫描时保持时间步、网格、湍流模型完全一致只变攻角和空化数否则结论会互相污染没法归因。5. 云空化非定常结果的验证与两个易踩的坑云空化算完不能只看几张漂亮的汽相云图得做定量的非定常验证。第一步是提取监测点的压力脉动或某截面汽相体积分数时程做 FFT 拿脱落频率再折算斯特劳哈尔数 St f·c/v∞。NACA0015 在 σ≈0.9、攻角 8° 附近文献 St 通常在 0.1~0.2 量级如果算出来差三五倍八成不是物理问题是数值问题。import numpy as np # signal: 某监测点汽相体积分数时程, dt: 采样间隔 freq np.fft.rfftfreq(len(signal), ddt) amp np.abs(np.fft.rfft(signal - signal.mean())) peak freq[np.argmax(amp[1:])1] St peak * c / v_inf print(f脱落频率 {peak:.2f} Hz, St {St:.3f})逻辑说明先减去均值去掉直流再取实FFT看主峰。参数说明dt 是写盘间隔不是时间步要保证一个脱落周期内至少 50 个采样点c 和 v∞ 用各自算例的值别套错翼型。第一个坑是回射流没捕捉到就下结论说云空化机理。回射流是云空化脱落的关键需要足够的流向分辨率和合适的输运系数Cv 给太小回射流会被抹平看云图会觉得空泡直接溃灭。把吸力面后半段也加密、Cv 提到 1.5 左右再跑一轮对比一般能明显看到反向流动。第二个坑是空化数用了进口压力而不是出口压力做参考。进口给速度、出口给压力才闭合若两头都给压力进口速度会自己找平衡实际 v∞ 就不是 12 m/s 了空化数跟着漂。验证方法很简单跑稳定后从后处理里读进口截面的平均速度看是不是 12 m/s不是就回到边界条件重设。最后一招是把翼型吸力面最小压力系数 Cp_min 读出来跟空化数比较Cp_min 的绝对值小于 σ 时理论上不该有明显空化用这个做算例成立性自检比看云图靠谱。本文还有配套的精品资源点击获取 NACA0015水翼云空化CFD仿真：边界条件与空化模型选型
简介这份专题资料面向水动力学、船舶工程与流体力学方向的研究人员和学生围绕NACA0015二维水翼的空化现象及其控制机理展开整理可帮助读者理解空化数值计算的建模思路与参数选取依据。内容涵盖计算域设置弦长c0.1m攻角6°~10°进口距前缘5倍弦长、出口距尾缘10倍弦长、上下边界相距12倍弦长、空化数定义与取值空化数0.9、进口速度12m/s、汽化压强0.044bar、时间步长10⁻⁵s并延伸至Clark-Y翼型弦长70mm、三维展宽68mm及700mm×190mm矩形二维计算域的替代方案还涉及NACA66(MOD)的f/c0.020、a0.8、t/c0.9等翼型参数与汽化压强、空化数关系。资源为单个docx文档压缩包约332KB结构紧凑便于检索查阅已有112人学习。适合需要参考翼型空化建模、工况设定与气泡动力学分析的读者作为专题笔记使用。1. 从一份 NACA0015 空化计算笔记说起算例能跑通靠的是边界条件不是网格数很多人第一次做水翼空化仿真卡住的地方不是网格画不细而是算出来的空化数对不上、云空化脱落周期跟文献差一个量级。这份 2021-2022 年的专题笔记里NACA0015 二维水翼算例其实已经把关键量给全了弦长 c0.1m攻角 6°~10°弧形进口边距前缘 5c出口边距尾缘 10c上下边界相距 12c翼型居中空化数取 0.9进口速度 12m/s汽化压强 0.044bar时间步长 1e-5 s。这些数字看着琐碎但每一个都直接决定空化数是否算得准、回流水汽是否混得进去。这份材料适合两类人一类是做船舶推进、水下航行器、水泵诱导轮空化分析的研究生和工程师另一类是已经会点 CFD 但总在空化数守恒和云空化非定常捕捉上翻车的人。Clark-Y 翼型弦长 70mm、三维展向 68mm、NACA66(MOD) 的 f/c0.020、a0.8、t/c0.9 这些参数也在笔记里出现说明作者是在做多翼型对比。2. NACA0015 二维计算域搭建与空化数标定2.1 计算域几何为什么按 5c/10c/12c 布笔记里的域布置不是随手写的。进口放到前缘上游 5 倍弦长是为了让来流在到达翼型前有足够距离发展成均匀流避免把入口的数值扰动直接怼到前缘出口放到尾缘下游 10 倍弦长是给尾迹和脱落的空泡云留出耗散空间否则出口边界会把未衰减的涡结构反射回计算域云空化脱落频率就假了。上下边界相距 12 倍弦长、翼型居中得到 6c 的半无限感阻塞比对 NACA0015 在 6°~10° 攻角下阻塞效应可以忽略。边界位置类型备注进口前缘上游 5c速度入口v12 m/s出口尾缘下游 10c压力出口由空化数反算上/下翼型上下各 6c滑移/对称视求解器而定翼型弦长 c0.1m无滑移壁面y 需匹配湍流模型弧形进口边是笔记里特意提的细节常见做法是把入口段做成以翼型前缘为圆心的圆弧使来流方向沿径向指向前缘攻角变化时不用重画域只改入口速度矢量方向即可做 6°、9° 多攻角扫描时省事。2.2 空化数定义与出口压力反算空化数是这一整套算例的锚点笔记给的表达式是$$\sigma \frac{p_\infty - p_v}{\tfrac{1}{2}\rho v_\infty^2}$$其中 v∞ 是进口速度ρ 是水密度pv 是汽化压强笔记取 0.044barp∞ 是参考压强工程上常用出口压力代替。给定 σ0.9、v∞12 m/s就能反算出口压力。水密度取 998.2 kg/m³pv0.044bar4400 Parho 998.2 # 水密度 kg/m^3 v_inf 12.0 # 进口速度 m/s p_v 4400.0 # 汽化压强 Pa (0.044 bar) sigma 0.9 # 目标空化数 q 0.5 * rho * v_inf**2 # 动压 p_out p_v sigma * q # 由 sigma 定义反算出口压力 print(f动压 q {q:.1f} Pa) print(f出口压力 p_out {p_out:.1f} Pa {p_out/1e5:.4f} bar)逻辑说明空化数定义里的 p∞ 用的是远场参考压力二维外流里最方便的就是把压力出口设成这个值进口只给速度。参数说明σ 越小越容易空化0.9 属于中等空化强度能稳定出现片状空化并周期性脱落把 σ 从 0.9 调到 1.5 左右通常片空化会缩短甚至消失调到 0.5 以下容易发展成超空化。这一步算错了后面看到的所有云空化图像都不可信。提示pv 强烈依赖水温。笔记给的 0.044bar 对应约 30℃ 的水20℃ 时 pv 约 2340 Pa。如果仿真里水温设错空化数会系统性偏移脱落频率跟着一起偏。2.3 时间步长与库朗数约束笔记给的时间步长是 1e-5 s。这个值不是拍脑袋定的它要满足库朗数 Co v·Δt/Δx ≤ 1最好控制在 0.5 附近。进口 12 m/s若前缘附近最小网格 0.5mm则 Co 12×1e-5/5e-4 0.24属于比较稳的范围。若把网格加密到 0.1mm 而时间步不变Co 会跳到 1.2非定常空化脱落就会数值发散或频率失真。# 以 OpenFOAM interCondensatingEvaporatingFoam 为例检查时间步与库朗数 # controlDict 关键项 deltaT 1e-5; # 笔记给定时间步长 adjustTimeStep no; # 空化非定常算例一般不自动调步 maxCo 0.5; # 若开启 adjustTimeStep 时的目标库朗数 writeInterval 0.01; # 脱落周期约几十 ms取足够密参数说明adjustTimeStep 关掉是为了让脱落的频谱分析有均匀采样自动变步会让 FFT 出现频率涂抹。writeInterval 建议让一个脱落周期内至少存 50 帧否则 PSD 峰位置会糊。3. 求解器选型与空化模型参数落地3.1 均质平衡流模型与输运方程模型怎么选水翼云空化是典型的两相可压缩问题主流两条路一是均质平衡流正压状态方程二是输运方程类空化模型Kunz、Schnerr-Sauer、Zwart-Gerber-Belamri。均质平衡流把混合物密度写成压力的函数算得快、对片空化稳定但蒸发冷凝的滞后效应表达不了云空化脱落的细节偏弱。输运方程模型显式求解汽相体积分数输运能捕捉脱落和回射流代价是要调蒸发/冷凝系数。笔记的场景是研究控制机理需要看脱落和回射流所以我一般选 Schnerr-Sauer 或 Zwart。Schnerr-Sauer 里的成核点密度参数很关键# Schnerr-Sauer 关键参数constant/phaseProperties 或 cavitationProperties Cv 1.0; // 冷凝系数 Cc 1.0; // 蒸发系数 n 1.0e13; // 单位体积成核点数量 1/m^3 dNuc 1.0e-6; // 成核点初始直径 m参数说明n 和 dNuc 控制初始汽核含量n 调大空化更容易触发、片空化更靠前Cv、Cc 分别控制冷凝和蒸发速率抖得厉害时把 Cc 降到 0.5 往往能压住。Zwart 模型对应的 Cvap、Ccond 通常在 50 和 0.01 量级起步再调。模型蒸发/冷凝系数适用云空化脱落Schnerr-SauerCc/Cv ~1.0片空化、中高 σ可捕捉ZwartCvap50, Ccond0.01通用可捕捉Kunz率常数调稳定计算稍弱均质平衡流无快速扫参偏弱3.2 湍流模型与壁面处理空化区域是强逆压梯度加相变SST k-ω 是这类算例的常规选择因为它对逆压梯度和壁面附近处理比标准 k-ε 好。第一层网格高度按 y 定SST 用壁面函数时 y 落在 30~100 即可若要做低雷诺数解析则 y1网格量会涨一个量级。二维 NACA0015 弦长 100mm我把第一层高度定在 0.05mm 左右对应 y≈40既能抓住前缘吸力面的片空化又不至于让网格在三万以内下不来。# snappyHexMesh 局部加密片段前缘与吸力面加密 refinementRegions { frontRegion { mode distance; levels ((0.02 4) (0.005 6)); // 前缘 2cm 内 4 级, 5mm 内 6 级 } suctionSide { mode distance; levels ((0.01 5)); // 吸力面周围 1cm 内 5 级 } }逻辑说明空化主要在吸力面前半段和尾缘附近发生全局加密性价比低用距离场在翼型周围做局部加密能在网格量可控的前提下把空化界面捕捉住。参数说明levels 里前一个数是距离阈值单位 m后一个数是加密级数每加一级网格尺度减半。前缘最外 2cm 给 4 级、5mm 内给 6 级是二维算例比较省的配法。注意空化区汽液界面处的网格分辨率直接决定脱落形态。若只在前缘加密而不在尾缘加密脱落的空泡云会过早耗散看起来像不脱落。4. Clark-Y、NACA66(MOD) 多翼型算例与工况对比4.1 Clark-Y 弦长 70mm、三维展向 68mm 的域怎么摆笔记里 Clark-Y 的设定是弦长 70mm、三维展宽 68mm弦展比接近 1属于低展弦比矩形翼。二维域可以照 NACA0015 的方法按 5c/10c 缩放也可以用笔记给的矩形域长 700mm、高 190mm翼型前端距进口 130mm、高度居中、攻角 0顺流平放。这个域的阻塞比是 70/190≈0.37比 NACA0015 的域大不少做零攻角基准工况够用但做 9° 攻角时上边界离吸力面太近阻塞会抬高空化数阈值。参数NACA0015Clark-Y弦长 c0.1 m70 mm域长×高—700×190 mm前缘距进口5c130 mm攻角6°~10°0°平放展向二维68 mm三维三维计算要注意展向端面。若不做周期性或对称处理端面回流会从翼尖附近倒灌进来污染空化区。常见做法是端面设成对称面或自由滑移再把展向取到足够长让端面影响远离主流空化区。4.2 NACA66(MOD) 参数含义与建模注意笔记给的 NACA66(MOD) 参数是 f/c0.020、a0.8、t/c0.9。这里要小心三个量的常规含义分别是最大弯度比 f/c、中弧线相对位置 a、厚度比 t/c。但 t/c0.9 对水翼来说是极厚的不太可能一般水翼 t/c 在 0.05~0.15怀疑笔记里 t/c 记的是别的量或笔误。实际建模前应该回到原始翼型坐标库核对NACA66 mod 属于 NACA 6 系列修正型常用中弧线 a0.8 指中弧线在 80% 弦长附近达到最大弯度。# NACA 6 系列中弧线近似a0.8, f/c0.020示例 import numpy as np def naca6_meanline(xc, a0.8, fc0.020): # xc: 弦向坐标 0~1 yc np.zeros_like(xc) # a0.8 时最大弯度位置约在 xca for i, x in enumerate(xc): if x a: yc[i] fc / a**2 * (2*a*x - x**2) else: yc[i] fc / (1-a)**2 * ((1-2*a) 2*a*x - x**2) return yc x np.linspace(0, 1, 101)逻辑说明这段是给建模脚本用的中弧线生成把 f/c 和 a 代进去就能得到弯度分布再叠加厚度分布得到完整翼型。参数说明a 越大最大弯度越靠后前缘载荷越低、空化起始越靠后跟 NACA0015 对比正好能看出弯度位置对空化起始的影响。所有坐标生成完建议先画出来跟翼型库叠一下确认没跑偏再进网格。4.3 多工况扫描的批处理思路做 0°、6°、9° 攻角加多个空化数的矩阵扫描一算就是十几个算例。手动改文件不现实常规做法是脚本化生成 casefor aoa in 0 3 6 9; do for sig in 0.5 0.9 1.5; do case${aoa}deg_sigma${sig} cp -r template_case $case # 改入口速度方向为攻角 sed -i s/ALPHA/$aoa/ $case/0/U # 改出口压力为对应空化数反算值 python set_sigma.py --case $case --sigma $sig done done参数说明aoa 控制入口速度矢量方向注意改的是矢量方向而不是网格网格在攻角扫描时保持不变最省事sig 通过 set_sigma.py 按 2.2 的公式反算出口压力写进 p 文件。每个 case 跑完统一后处理把吸力面压力系数分布和汽相体积分数时程抽出来对比才能看出攻角和空化数哪个对脱落频率影响更大。提示多工况扫描时保持时间步、网格、湍流模型完全一致只变攻角和空化数否则结论会互相污染没法归因。5. 云空化非定常结果的验证与两个易踩的坑云空化算完不能只看几张漂亮的汽相云图得做定量的非定常验证。第一步是提取监测点的压力脉动或某截面汽相体积分数时程做 FFT 拿脱落频率再折算斯特劳哈尔数 St f·c/v∞。NACA0015 在 σ≈0.9、攻角 8° 附近文献 St 通常在 0.1~0.2 量级如果算出来差三五倍八成不是物理问题是数值问题。import numpy as np # signal: 某监测点汽相体积分数时程, dt: 采样间隔 freq np.fft.rfftfreq(len(signal), ddt) amp np.abs(np.fft.rfft(signal - signal.mean())) peak freq[np.argmax(amp[1:])1] St peak * c / v_inf print(f脱落频率 {peak:.2f} Hz, St {St:.3f})逻辑说明先减去均值去掉直流再取实FFT看主峰。参数说明dt 是写盘间隔不是时间步要保证一个脱落周期内至少 50 个采样点c 和 v∞ 用各自算例的值别套错翼型。第一个坑是回射流没捕捉到就下结论说云空化机理。回射流是云空化脱落的关键需要足够的流向分辨率和合适的输运系数Cv 给太小回射流会被抹平看云图会觉得空泡直接溃灭。把吸力面后半段也加密、Cv 提到 1.5 左右再跑一轮对比一般能明显看到反向流动。第二个坑是空化数用了进口压力而不是出口压力做参考。进口给速度、出口给压力才闭合若两头都给压力进口速度会自己找平衡实际 v∞ 就不是 12 m/s 了空化数跟着漂。验证方法很简单跑稳定后从后处理里读进口截面的平均速度看是不是 12 m/s不是就回到边界条件重设。最后一招是把翼型吸力面最小压力系数 Cp_min 读出来跟空化数比较Cp_min 的绝对值小于 σ 时理论上不该有明显空化用这个做算例成立性自检比看云图靠谱。本文还有配套的精品资源点击获取