Froude数 — CAE用语解说

分类: 术语集 | 2026-01-15

CAE visualization for froude number 2 - technical simulation diagram

Froude数

🧑‍🎓

老师，我听说Froude数在船舶设计中很重要，请问这是什么数呢？

Froude数的理论基础

Froude数的定义和物理意义

🧑‍🎓

教科书说Froude数是"惯性力与重力比值的平方根"，但具体在哪些情况下很重要呢？

🎓

在具有自由表面的流动中，例如船舶造波阻力或大坝泄放流，它起着本质作用。具体来说，当长度尺度为

$$ L $$

，代表速度为

$$ U $$

时，Froude数

$$ Fr = U / \sqrt{gL} $$

在接近1时，流动状态会发生戏剧性变化。例如，在明渠流中，

$$ Fr < 1 $$

为常流（亚临界流），

$$ Fr > 1 $$

为射流（超临界流）。

🧑‍🎓

你说"接近1"，那么Fr=0.9和Fr=1.1具体有什么区别呢？

🎓

最大的区别在于下游的影响是否能传回上游。在常流Fr=0.9中，水面上的小扰动（例如扔石头）会作为波向上游和下游传播。而在射流Fr=1.1中，这种扰动只向下游传播。这是因为信息传播速度低于流速，是产生冲击波（水力跳跃）的条件。

🧑‍🎓

经常听到的船舶阻力中的"造波阻力"与Froude数有什么关系？

🎓

船体产生的波长

$$ \lambda $$

与

$$ \lambda = 2\pi U^2 / g $$

成正比。当Froude数

$$ Fr = U / \sqrt{gL_{pp}} $$

（

$$ L_{pp} $$

为垂线间长）超过约0.4时，造波阻力急剧增加。例如，全长100m的集装箱船以24节（约12.3 m/s）航行时，Fr约为0.25，但高速渡轮的Fr可能超过0.6，此时造波阻力占总阻力的过半。

Froude数的数值计算方法

CFD中Froude数的处理

🧑‍🎓

用CFD模拟自由表面流时，怎样设置Froude数？是直接作为输入参数指定吗？

🎓

不同求解器有不同做法。多数情况下，分别输入重力加速度

$$ g $$

、流入速度

$$ U $$

和代表长度，求解器内部计算Froude数。例如用Ansys Fluent的VOF法模拟溃坝，需在"Operating Conditions"中设置重力矢量，在入口边界条件中定义速度。Fr是得到的无量纲数。

🧑‍🎓

在模型缩放实验中使用Froude数相似准则，但CFD不需要这样做对吗？

🎓

CFD的优势是可以按实际尺度计算，但为了验证，Froude数相似准则很重要。例如，从1/25缩尺模型的实验中获得数据时，CFD也应按同样缩尺计算，使用相似准则确定的速度（

$$ U_{model} = U_{real} / \sqrt{25} $$

），对比自由表面形状和阻力值。Star-CCM+的"Scalar Parameter"功能可以参数化代表长度和速度，保持Fr恒定进行灵敏度分析。

🧑‍🎓

在非定常计算中，Fr随时间变化的现象（如加速的船）模拟时，需要注意什么？

🎓

关键是时间步长

$$ \Delta t $$

的设置。要捕捉自由表面的变化，需考虑表面重力波的传播时间尺度。一般经验是

$$ \Delta t < 0.01 \times L / \sqrt{gL} $$

（L为代表长度）。同时，使用动态网格或重叠网格时，网格变形速度需能跟随Fr的变化，特别要注意临界Fr=1附近的网格质量退化问题。

Froude数的实务应用

自由表面流分析工作流程

🧑‍🎓

用VOF法分析船体周围流场的具体步骤中，在哪些地方需要考虑Froude数？

🎓

主要有3个环节。1. **前处理**：计算域应在船体后方至少延伸船长的3～5倍，水面上下也要有充分的空间。Fr越高，波延伸越远。2. **初始条件**：通常从静水中的船体开始计算"加速法"，包括达到目标Fr过程中的过渡状态。3. **后处理**：根据Fr评估造波阻力贡献，对船体表面压力积分和用动量理论计算阻力进行交叉验证。

🧑‍🎓

网格如何划分？特别是水面附近是否需要更细？

🎓

必须要细。在预期自由表面存在的区域，至少要有10～20层细网格，以提高界面捕捉精度。用Ansys Meshing的"Inflation"或Siemens Star-CCM+的"Prism Layer Mesher"生成。网格尺寸目标通常为预期波长

$$ \lambda $$

的1/20以下。例如Fr=0.3、船长200m时，

$$ \lambda \approx 70m $$

，水面附近网格应≤3.5m。

🧑‍🎓

计算容易发散的条件，比如Fr接近1的情况，有稳定的技巧吗？

🎓

有的。首先，建议从定常计算的伪时间步长法转向非定常计算开始。其次，速度和压力的耦合求解用"PISO"方案比"Coupled Algorithm"更稳定。另外，通过降低界面锐化系数（从默认1.0降至0.5左右）可缓和界面陡峭性。OpenFOAM的`interFoam`求解器中，增加`nAlphaCorr`（界面重构迭代次数）也有效。

Froude数的软件对比

主要CAE软件中的自由表面流分析功能

🧑‍🎓

Ansys Fluent、Star-CCM+和OpenFOAM在Froude数相关的自由表面流分析中有差异吗？

🎓

核心VOF法思路相同，但实现和易用性有差异。**Ansys Fluent**在"Multiphase Model"中可选显式或隐式VOF，并提供自定义函数定义"Froude Number"用于造波阻力监测。**Siemens Star-CCM+**有专用"VOF Waves"物理模型，可直接设置规则波或不规则波作为入射边界条件，内置Fr相似准则缩放功能。**OpenFOAM**的`interFoam`可在代码级精细控制界面补正算法（MULES）参数，但初始配置难度较高。

🧑‍🎓

专用船舶分析软件（如SHIPFlow或Fine/Marine）相比通用CFD软件有什么优势？

🎓

最大优势是船体阻力成分分解的自动化。SHIPFLOW（Flowtech）采用"分割法"，分别用RANS法计算粘性阻力、用势流理论计算造波阻力，高效求得造波阻力系数

$$ C_w(Fr) $$

。此外，内置ITTC（国际船舶试验会议）推荐的模型尺度补正（1978年性能预报法），能实现符合工程实践的阻力和推进性能预报。通用CFD需要用户自己用脚本构建这些前后处理。

🧑‍🎓

与结构分析耦合（流体-结构耦合：FSI）中，Fr有影响的案例吗？

🎓

高速艇和帆船的帆颤振分析是典型例子。使用Ansys Workbench或SIMULIA（Abaqus/CFD）耦合平台。重点是流体荷载（压力）的变动频率是否与结构固有频率一致。Fr变化会改变尾流涡脱频率（与Strouhal数相关）和造波模式，进而改变激励力谱。特别是Fr在0.4～0.6范围时，船首波与船尾波干涉会产生强激励振动，需用耦合分析评估。

Froude数的故障排除

自由表面流分析的常见错误和对策

🧑‍🎓

计算过程中自由表面越来越扩散，看不出明确的界面，原因是什么？

🎓

主要有3个原因。1. **网格太粗**：如前所述，需相对波长足够细的网格。2. **界面锐化不足**：调整Fluent的"Sharpening/Damping Factors"或Star-CCM+的"Interface Compression"（通常以1.0为基准增减）。3. **时间步太大**：界面追踪需CFL数（Courant数）≤1，理想≤0.5。特别是Fr较大（高流速）的计算，用自动调整时间步的功能效果更好。

🧑‍🎓

造波阻力值与实验值或经验公式对比明显偏大（或偏小），应该查哪里？

🎓

首先查**计算域出口边界**。如果出口离船体太近，或用静压条件，波会在出口反射回船体，导致阻力过大估计。应设置**波吸收区**（Wave Absorption）或用Orlanski型等放射边界条件。其次查**湍流模型**。船体附近边界层解析不足，粘性阻力评估不准，造波阻力分离精度下降。最后查**收敛性**。非定常计算需至少10个波长时间的计算，对阻力值取时间平均。

🧑‍🎓

以Fr为参数做一系列计算，某个点（如Fr=0.45）计算发散了，为什么会这样？

🎓

很可能在该Fr下发生了某种**共振现象**或**不稳定流动模式**。例如，船体产生的波干涉模式在某些Fr（"Hollow干涉条件"）下发生剧烈变化，流场变得极其敏感。对策包括：1) 仅对发散的Fr用更小时间步从头计算，2) 从收敛的邻近Fr结果继承初条件（重启计算），3) 临时增强数值耗散（从高阶格式改为1阶迎风差分）抑制不稳定。这类奇异点本身往往是重要物理发现。