8.3 航空船舶与工业流动实战
8.3 航空船舶与工业流动实战
本节摘要：方法与模型都要到现场兑现。本节做三个完整案例：给机翼做攻角扫描，看升力线斜率、失速与阻力极曲线怎么从 CFD 数据读出来；按拖曳水池流程做船模-实船相关外插，体会第 7 章"Re 与 Fr 无法同时满足"的工程解法；对街区布局做风环境评估，用行人高度风速比判断舒适度。三个现场共用一条工作流：明确问题与判据 → 建模与网格 → 验证 → 参数扫描 → 解读成设计语言。
前两节搭好了求解器与模型选型的骨架，但工程交付从不以"残差收敛"结束——客户要的是：这副机翼在哪个攻角失速、这型船实船油耗多少、这块街区冬天会不会刮穿堂风。把 CFD 输出翻译成这些决策语言，才是最后一公里。
案例一：机翼攻角扫描与极曲线
任务：某低速无人机翼型，攻角从 -4 度扫到 16 度，求升力曲线与失速点。流程是 RANS（SST 模型）二维定常起步，大攻角分离失稳后切换非定常取时均。翼型数据的典型读法：
import numpy as np alpha = np.array([-4, 0, 4, 8, 12, 14, 16]) # 攻角 cl = np.array([-0.38, 0.12, 0.62, 1.08, 1.42, 1.38, 1.05]) cd = np.array([0.010, 0.008, 0.009, 0.014, 0.028, 0.042, 0.075]) slope = np.polyfit(np.radians(alpha[:4]), cl[:4], 1)[0] # 线性段升力线斜率 cl_cd = cl / cd imax = int(np.argmax(cl_cd)) print(f"升力线斜率约 {slope:.2f} 每弧度 二维薄翼理论值约 6.28") print(f"最大升阻比 {cl_cd[imax]:.1f} 出现在攻角 {alpha[imax]} 度") print(f"升力峰在 12 度之后回落 分离启动 失速攻角约 13 到 14 度")
三条判读纪律：线性段斜率对不对（二维理论约 2π，偏低说明网格或湍流模型有系统误差）；失速攻角偏保守还是乐观（RANS 普遍晚报失速，实飞要留余量）；极曲线峰值附近的数据密度（设计点就在最大升阻比附近，那里点要密）。报告附上与风洞数据的对标误差，升力 3% 以内、失速角误差 1 度以内才算立住了 credibility。
案例二：船模拖曳与尺度外插
第 7 章留下的问题在这里清算：模型水池里 Re 与 Fr 不可能同时等于实船值，习惯做法是保 Fr（波浪阻力由重力主导，弗劳德相似是命门），粘性部分用经验外插。标准流程（ITTC-57 方法骨架）：
# 船模 6 m 实船 180 m 同一弗劳德数下阻力换算 Lm, Ls = 6.0, 180.0 lmd = Lm / Ls Vm = 1.5 # 模型拖速 Vs = Vm / np.sqrt(lmd) # Fr 相等 速度按尺度根号放大 print(f"实船对应航速 {Vs:.1f} m/s 约 {Vs*1.944:.1f} 节") Cfm = 0.075 / (np.log10(Vm*Lm/1.14e-6))**2 # 模型摩擦阻力系数 经验式 Cfs = 0.075 / (np.log10(Vs*Ls/1.14e-6))**2 # 实船摩擦阻力系数 print(f"模型 Cf = {Cfm:.5f} 实船 Cf = {Cfs:.5f} 相差 {Cfm/Cfs:.2f} 倍") print("总阻力 = 摩擦差外插加波形阻力 摩擦不能直接放大 正是尺度效应的核心")
模型 Re 比实船小几十倍，摩擦系数大得多——直接按几何比例放大阻力会严重高估实船油耗。外插公式的形状（平板公式加粗糙度补贴加形状因子）各家船研所略有差异，这正是实船试航验收与模型预报总要对账、留裕度的地方。
现场解读：三个案例共用的一条工作流
明确判据 攻角谱 油耗 风速比 建模 简化维度 定湍流模型 验证 对标理论斜率 或 水池数据 扫描 一个变量网格采密 解读 翻译成设计余量与决策
案例三：街区风环境评估
任务：新住宅片区布局报批，需给出行人高度（1.5 m）风速比 R = U_local/U_ref 的分布，判断舒适度等级（步行通勤一般要求 R 均值低于 0.5 量级、无局部超过 1.2 的狭管效应点）。建模取三维 RANS，来流按风速廓线（幂律指数 0.22 的城郊地形），十六个风向扫两个代表性方位：
import numpy as np # 采样点风速比示意 由后处理脚本从场数据提取 R_north = np.array([0.32, 0.41, 0.38, 0.55, 0.29]) # 北风五测点 R_west = np.array([0.61, 1.28, 0.44, 0.37, 0.33]) # 西风 楼间狭管点超标 print(f"北风 R 均值 {R_north.mean():.2f} 合格") print(f"西风最大 R = {R_west.max():.2f} 位于两塔楼之间 需加密绿化或退台")
西风向下两栋塔楼形成狭管，R 达 1.28——比来流还快，行人会有明显推搡感。改造选项按代价排序：底层裙房退台、塔楼错位布置、加防风林带。重跑后 R 降到 0.7 以下，规划部门按此出条件。这个案例的启示是：CFD 输出不是等值线图，是可以改规划的红笔。
⚠️ 风环境结论对来流廓线与湍流度敏感，报批件必须注明气象站数据来源与统计年限，否则数字没有法律效力。
变式思考
先用一张图把三个现场的输出形态放在一起，方便形成条件反射：
三个案例换汤不换药：机翼是"性能曲线扫描"，船是"尺度相关外插"，街区是"判据图叠加"。同一套方法骨架还能搬去数据中心热管理（判据换成机柜进风温度）、隧道通风（换成污染物浓度）、水泵叶轮（换成效率与空化余量）。学会把工程问题翻译成"判据加扫描加验证"的流程，比记住任何软件操作都值钱。
本节要点回顾
- 机翼案例：升力线斜率对标理论、失速要留余量、极曲线峰值是设计点；
- 船舶案例：保 Fr 舍 Re，摩擦系数用经验式外插，尺度效应是模型试验的宿命；
- 风环境案例：判据是行人高度风速比，输出要能改规划，气象依据决定法律效力；
- 共用工作流：判据 → 建模 → 验证 → 扫描 → 决策语言，五步闭环；
- 全书到此合流：物性、守恒、边界层、相似准则、数值方法在同一份工程交付里各就各位。