力与力矩系数：结果诊断与可信度验证

力与力矩系数的诊断难点不在 functionObject 能不能跑通，而在参考量定义是否与文献一致。同一组力，换参考面积或参考速度，$C_D$ 可以差 25% 以上。本文给出 forces 的真实配置、系数定义与一次可核对的动压手算，并把参考面积、参考速度、力矩参考点三类口径陷阱逐条落到判定试验上。

forces functionObject 的真实配置

functions
{
    forces
    {
        type            forces;
        libs            ("libforces.so");
        patches         (walls);
        rho             rhoInf;
        rhoInf          1.2;
        CofR            (0 0 0);
        log             true;
        writeControl    timeStep;
        writeInterval   1;
    }
}

rho rhoInf 让工具用常数密度 1.2 kg/m³ 计算动量通量，这对不可压求解器是必需的；写成 rho rho 时工具会去读密度场，不可压算例没有该场会直接报错。CofR 是力矩参考点，力矩对它的敏感度是线性的：把参考点沿 x 平移 0.1 m，$C_M$会多出$0.1 F_L$ 的贡献。

系数定义与参考量

$$C_D = \frac{F_D}{\tfrac12 \rho_\infty U_\infty^2 A_{ref}}, \qquad C_L = \frac{F_L}{\tfrac12 \rho_\infty U_\infty^2 A_{ref}}$$ $$C_M = \frac{M_z}{\tfrac12 \rho_\infty U_\infty^2 A_{ref}\, c}$$

动压项用远场量，$A_{ref}$ 与 $c$ 由几何约定给出。以 $\rho_\infty = 1.2\ \mathrm{kg/m^3}$、$U_\infty = 20\ \mathrm{m/s}$、$A_{ref} = 0.05\ \mathrm{m^2}$ 为例：

$$q_\infty = \tfrac12 \times 1.2 \times 20^2 = 240\ \mathrm{Pa}, \qquad q_\infty A_{ref} = 12.0\ \mathrm{N}$$

若 forces.dat 的阻力列读到 $F_D = 7.8\ \mathrm{N}$，则$C_D = 7.8/12.0 = 0.65$。这个 0.65 只在参考面积等于 0.05 m² 时成立；若文献用的是迎风面积 0.04 m²，同一个$F_D$给出$C_D = 0.8125$，相差 25%。

三种参考量口径的陷阱

• 参考面积：翼型用弦长乘单位展长，圆柱用直径乘单位展长，汽车常用迎风面积。报告必须同时给出 $A_{ref}$ 数值。
• 参考速度：外流用自由流 $U_\infty$，管内流常改用体积平均速度 $\bar U$，两者在充分发展段可差 20%。
• 力矩参考点：$C_M$必须写明 `CofR`；气动中心处的$C_M$近似不随$C_L$ 变化，其他点不具备这一性质。

从压力积分独立复核力

若壁面阻力主要由压差贡献，可用表面积分复核：

$$F_x = -\int_A p\, n_x\, \mathrm{d}A \approx -\sum_f p_f S_{f,x}$$

某算例用 forces 得到 $F_x = 7.8\ \mathrm{N}$，再对同一 patch 做 areaIntegrate(p) 投影得 $-7.6\ \mathrm{N}$，差的 0.2 N 来自黏性剪切与法向投影。两者应在 5% 内一致，否则先怀疑 patch 集合不同——例如 walls 是否漏掉了端板。

postProcess -func "forces(patch=walls,rho=rhoInf,rhoInf=1.2,CofR=(0 0 0))" -time 100:500
postProcess -func "patchIntegrate(patch=walls,fields=(p))" -time 500

第一条命令复用同一组参考量重算系数，第二条给出压力积分作为独立对照。

瞬态算例的统计口径

非定常绕流的 $C_L$ 是周期量，报告应给出时间平均值与脉动幅值：

$$\overline{C_L} = \frac{1}{T}\int_{t_0}^{t_0+T} C_L(t)\,\mathrm{d}t, \qquad C_{L,rms} = \sqrt{\overline{(C_L-\overline{C_L})^2}}$$

某圆柱在 $U_\infty = 1.0\ \mathrm{m/s}$、$D = 0.1\ \mathrm{m}$ 下测得 $St = fD/U = 0.198$，即脱落频率 $f \approx 1.98\ \mathrm{Hz}$，周期 0.505 s。统计窗口至少覆盖 10 个周期，即 $T \ge 5.05\ \mathrm{s}$，否则 $C_{L,rms}$ 的重复性会差 10% 以上。

失败模式

现象	根因	判定试验
日志报找不到 rho 场	rho rhoInf 误写成 rho rho	改回 rhoInf 并保留 rhoInf 1.2，重跑一步
$C_D$ 与文献差近 4 倍	参考面积口径不同（迎风 vs 湿润面积）	用几何算出 $A_{ref}$ 反推 $F_D$，与 forces 输出比对
$C_M$ 随迭代线性漂移	CofR 未设或落在默认原点	打印 CofR，用 $M = M_0 + \mathbf{r}\times\mathbf{F}$ 换算到目标参考点
$C_D$ 出现负值	patch 法向朝内或坐标轴与阻力方向相反	检查 patch 法向，核对 forces 的方向约定
瞬态 $C_L$ 幅值不可重复	统计窗口不足或采样频率低于脱落频率	把窗口扩到 10 个周期、采样改为每步输出

参考文献

1. S. F. Hoerner, Fluid-Dynamic Drag, Hoerner Fluid Dynamics, 1965.
2. R. D. Blevins, Applied Fluid Dynamics Handbook, Van Nostrand Reinhold, 1984.
3. B. R. Munson, T. H. Okiishi, W. W. Huebsch, A. P. Rothmayer, Fundamentals of Fluid Mechanics, 7th ed., Wiley, 2013.
4. OpenFOAM Foundation, OpenFOAM User Guide, §“Forces and force coefficients”.
5. C. Norberg, “Fluctuating lift on a circular cylinder: review and new measurements”, Journal of Fluids and Structures, 17(1), 2003.