粘性尺度效应对球鼻艏砰击压力及全尺度外推的多尺度CFD研究

《Journal of Marine Science and Engineering》：Multi-Scale CFD Investigation of Viscous Scale Effects on Bulbous Bow Slamming Pressures and Full-Scale Extrapolation

　　

预测球鼻艏波激砰击压力对船舶结构安全至关重要。研究人员采用重叠网格（overset-grid）雷诺平均纳维-斯托克斯-流体体积法（RANS-VOF）框架，研究粘性尺度效应对球鼻艏砰击载荷的影响。在临界纵摇-升沉共振顶浪（λ/LWL=1.2）下，针对1:50、1:20、1:15和1:10（α=50、20、15、10）四种几何尺度比开展多尺度模拟。全局运动响应遵循弗劳德相似（Froude similitude），而局部动态砰击压力呈现显著尺度差异。较小物理模型发育相对更厚的粘性边界层，充当水动力缓冲层，降低峰值压力并拓宽脉冲时长。因此，由小尺度模型直接弗劳德缩放会低估全尺度冲击载荷。为计入此类粘性尺度效应，研究人员提出基于多尺度回归的工程外推方法，在主冲击区所考察范围内表现出合理的线性相关性（R2=0.88–0.99）。该研究为球鼻艏砰击的尺度行为提供物理洞见，并为全尺度载荷估算提供实用参考。

研究背景与动机

船舶砰击是船体与水面的突发强非线性水动力冲击，球鼻艏因三维曲面前端在大幅海况运动中承受最严酷载荷。传统经验公式、切片理论、Wagner模型及线性势流求解器难以刻画自由面破碎、射流与复杂三维分离；RANS-VOF已成为极端波体相互作用的标准手段。工程上常借弗劳德相似由模型尺度外推全尺度，但砰击压力涉及流体-结构-空气耦合物理，对三维凸形球鼻艏，空气垫效应弱，粘性边界层厚度随雷诺数降低而增大，构成主要尺度偏差源。目前系统隔离三维凸形球鼻艏在升沉-纵摇耦合下粘性尺度效应的多尺度定量数据稀缺，直接弗劳德外推可能危险低估全尺度载荷。研究人员以一艘140.0 m设计水线高速瘦长排水型裸船体（带显著球鼻艏）为对象，在STAR-CCM+中建立重叠网格RANS-VOF数值框架，于四种几何尺度比（α=50、20、15、10）下，在λ/L_WL=1.2临界共振顶浪（实尺度波高H=9.75 m，海况8级）中量化偏离，并提出多尺度线性回归外推方法。论文发表于《Journal of Marine Science and Engineering》。

主要关键技术方法

研究人员采用STAR-CCM+求解瞬态不可压多相RANS方程，自由面用流体体积法（VOF）结合高分辨率界面捕获（HRIC）格式，湍流闭合用可实现k-ε两方程模型（Realizable k-ε Two-Layer）配合all-y⁺壁面处理；刚体运动由动态流体-体相互作用（DFBI）解算升沉与纵摇两自由度；计算域用静止背景块与随船重叠块（overset mesh）线性插值传递，按ITTC指南布置五阶斯托克斯波、松弛消波区与对称面；网格用各向异性裁剪单元，自由面纵向Δx:Δz=4:1，球鼻艏处逐级加密并铺棱柱层；验证包括空槽波面校准、30°楔形体入水压力系数对比Zhao和Faltinsen边界元解、三套网格GCI收敛（中网格3.43×10⁶单元作基线）、时间步1.0 ms与0.5 ms敏感性；尺度分析取第三至第五稳定遭遇周期峰值，砰击压力同时取点探头与0.6 m×0.6 m板面积分平均，脉冲用无量纲冲量I*与三角脉冲近似评估。

研究结果

波长影响：在1:20模型扫描λ/L_WL=1.0、1.2、1.4，λ/L_WL=1.2时遭遇频率使艏部相对垂向速度最大，Panel 1峰值C_p,max=10.97，判定为临界结构工况并用于后续尺度研究。

时空压力分布：砰击非空均匀，底中线Panel 1–5载荷最大，侧向因近垂直受载小，入水后期上舷因波爬升二次升压；沿中线由前向后（Panel 1→5、BottomP4→BottomP12）顺序湿润化致压力递减，前20%船长按压力集中；出水阶段出现C_p<0的负压吸力但高于汽化压力无空化。

板与探头压力：面积平均使主冲击区Panel 1–4峰值C_p较点探头降低3.2%–13.5%（1:20下Panel 1降11.2%），Panel 5因分离流略高7.7%；点峰值用于局部凹痕疲劳，面积平均用于板屈服与有限元边界更合理。

砰击冲量：定义I含升压时间τ_rise，CFD非对称脉冲前陡后拖尾，理想三角脉冲涵盖主冲击72.6%–80.6%动量；Probe 1峰值C_p比Probe 3高21%，但I几乎相等（0.0320对0.0318），表明前端点高尖短促、中部低宽长时，总动量补偿，峰值单指标不能代结构风险。

运动相似性：四尺度下无量纲升沉、纵摇曲线重合，相对1:20基线偏差≤9%（1:50升沉?8.6%源于绝对网格粗化），证实全局运动服从弗劳德相似。

压力系数尺度效应：C_p随α增大（模型变小）而降，Panel 1从α=10的12.00降至α=50的9.28，大尺度间非单调微振荡但总体下行，违背经典C_p不变假设。

尺度效应物理机制：弗劳德缩放下Re∝α^?1/2，小尺度Re低、边界层厚，作为水动力垫降低射流速度削峰；小尺度（1:50）脉冲峰低但 Rise/decay 宽，总冲量跨尺度变化小，呈粘性延寿耗能自补偿。

多尺度线性回归外推：对各测点建立弗劳德缩放物理压力p_scale=A·α+B，主冲击区Panel 1–3/Probe 1–3 R²=0.88–0.99（S_e<10.3 kPa），尾部Panel 4–5 R²=0.59–0.83因三维分离散度增大；外推α=1得Panel 1全尺度448.6±29.2 kPa、Probe 1为499.8±38.5 kPa（95%区间），而1:50直接弗劳德仅336.9 kPa，低估约32%。

讨论与结论翻译

讨论指出，线性外推在1:10–1:50区间适用（R²高），但粘性底层理论非线性（∝Re^?1/2），广尺度延拓需谨慎；尾部散射源于分离与尾迹；上舷Panel 11–25属波拍非冲击，次级行Panel 6–10仍保R²=0.93–0.98。结论：（1）λ/L_WL=1.2为临界砰击工况，相位叠加放大艏相对垂速；（2）峰值集中于球鼻艏前腹中线并向后衰减，上舷有波爬升二次压；（3）面积平均使主区峰值降11.2%–13.5%，三角冲量近似涵盖72.6%–80.6%主冲击动量；（4）小尺度厚边界层作缓冲削峰拓宽脉冲，直接弗劳德外推低估全尺度，四尺度线性回归在主区R²=0.88–0.99，外推Probe 1为499.8±38.5 kPa、Panel 1为448.6±29.2 kPa，可作初步工程参考；研究纯数值未获物理水池或实船校验，绝对精度待未来试验完善。