8.2 CFD 基础：从网格到解

8.2 CFD 基础：从网格到解

本节摘要：CFD 把守恒律搬进计算机：用网格切分流场，在每个单元上把控制方程写成进出账（有限体积法），迭代求解直到账目闭合。三个决定成败的环节是网格、离散与湍流模型——本节逐一过堂，并给出"一套计算做下来"的完整工序。

风洞的数据买不起太多，计算机用电力换数据。CFD 的核心承诺：同一套守恒律（第 3 章）在离散世界里依然成立——这正是 3.4 节坚持守恒形式的回报。

有限体积法：网格上的控制体

把流场切分成百万级的小单元（网格），每个单元就是一个 3.1 节的控制体。对每个单元写质量、动量、能量三笔账：单元内守计量随时间的变化 = 各面通量的净流入。通量用相邻单元的变量插值估算（离散格式），时间推进用显式或隐式格式（边界层这类刚性问题必须隐式）。守恒形式 + 有限体积 = 网格间通量严格抵消，全域质量、动量自动守恒，激波跳跃关系天然满足——3.4 节的理论承诺在此兑现。

解不可能一步到位：从初始猜想的流场出发，反复迭代，每轮用方程残差（账目的不平衡度）修正场变量，直到残差下降几个量级、监视量（升阻力、流量）稳定——这就是"算收敛了"。

网格：CFD 的地基工程

网格质量直接决定结果可信度，业内有"一半时间画网格"的说法。三种类型各有领地：

边界层内要有足够的网格层（壁面第一层高度按目标 y+ 设定，RANS 常取 30 到 300 的对数律区或小于 1 的低雷诺数解析），否则摩擦阻力与分离位置从第一步就偏。几何再复杂，贴壁那几层不能省——这是 5.1 节"薄却关键"在数值世界的镜像。

湍流模型：全谱系的成本阶梯

网格要细到什么程度才能"直接"算湍流？粗算：雷诺数百万级的流动，直接数值模拟（DNS）的网格量随 Re 的 2 上下次幂增长，工程问题需要数十亿到万亿单元——只有极少数科研算例玩得起。于是有了阶梯：

• DNS（直接数值模拟）：不解任何湍流模型，网格解析到最小湍流尺度。最保真、最昂贵，湍流研究用它生产"数值实验数据"。
• LES（大涡模拟）：大涡直接解析（几何决定的大结构），小涡用亚格子模型。分离流、非定常载荷的当前主力，成本在两者之间。
• RANS（雷诺平均）：把湍流脉动全部平均掉，用湍流模型（如两方程模型）补上雷诺应力。成本最低，工程设计的绝对主力——但模型是经验性的，分离区、强曲率区的误差可达两位数百分比。

选择逻辑是任务驱动：巡航升阻力用 RANS 够用；抖振、涡脱落、进场非定常载荷上 LES；湍流机理研究才请 DNS。"模型误差"不是 CFD 的缺陷而是它的生存方式——知道当前模型在你的流动上偏多少，比换更贵的模型更重要。

图：CFD 工作流——从几何到签发

动手算：给一次 RANS 计算设边界条件

某翼型 RANS 工况：巡航 230 m/s 真速、10 公里高度、弦长 5 m。确定输入参数：

import math
p, T = 26436.0, 223.15          # 10 km 标准大气
R, gma = 287.0, 1.4
rho = p / (R * T)                # 约 0.413 kg/m³
a = math.sqrt(gma * R * T)       # 约 299 m/s
mu = 1.42e-5                     # 223 K 下的粘度
Ma, Re = 230 / a, rho * 230 * 5.0 / mu
print(rho, a, Ma, Re)            # Ma 约 0.77，Re 约 3.3e7
​

边界条件应填：远场来流按 Ma 0.77 与 Re 3.3×10⁷ 组合（或直接给压强密度速度），壁面无滑移绝热，出口远场自由。注意雷诺数三千万——比海平面试验高一个量级，这份 CFD 的雷诺数优势恰恰是风洞最难补齐的部分。CFD 与风洞各守自己的长板：CFD 强在任意状态任意参数组合、全场细节；风洞强在真实物理（真实湍流、真实转捩）。下一节的质检体系就是让两份证据互相咬合。

DNS 到底贵到什么程度

给成本阶梯一个可感的数字。湍流里最小涡的尺度随雷诺数的负四分之三次方缩小，DNS 网格量约随雷诺数的二又四分之一次方增长：雷诺数一万的小算例，数百万网格、一台工作站过夜；雷诺数十万的槽道流，数十亿网格、超算数千核跑数周；雷诺数百万级的完整机翼 DNS，网格量进入万亿量级，世界顶级超算也只做过个别简化构型。这条标律把 DNS 钉死在"湍流机理研究的显微镜"位置上——工程直觉、实验对比、模型标定都靠它供弹药，但它自己永远不直接上生产线。

问题：商用求解器与开源求解器怎么选？

核心方程大同小异，差别在打磨程度与生态：商用工具的网格生成、界面、文档与验证手册厚实，适合以项目交付为节奏的工程团队；开源工具胜在开放（算法可改、可审计）与科研耦合，适合需要定制格式或研究新模型的人。比选型更要紧的是两件事：一是网格与工况的纪律（本章讲的三件套不因工具而打折），二是基准算例的回归测试——换工具、升版本都要先跑一遍已知答案的算例，确认"考卷还考得住"。工具是变量，纪律是常量。

问题：CFD 能完全替代风洞吗？

短期不能，长期是分工重构。CFD 的软肋恰是物理模型的硬伤区：转捩、大分离、激波/边界层干扰、气动声学，这些地方的预测误差仍是两位数百分比或根本不可信；风洞的软肋是雷诺数、成本与周期。行业现状是"CFD 海选、风洞终审、飞行试验盖章"，未来趋势是风洞点数大幅减少但精度要求更高（专门补 CFD 的不确定区）。替代论的对立面不是维持现状，而是把两者的互验流程（8.3 节）做得更便宜——裁决体系在进化，不是在换人。

本节要点

• 有限体积法 = 网格化的控制体记账，守恒形式保证全域守恒与激波正确性。
• 网格是地基：类型按外形选，贴壁边界层加密按 y+ 纪律，几何复杂不能牺牲贴壁层。
• 湍流处理三档：DNS 保真最贵、LES 折中、RANS 工程主力；模型误差要随身携带。
• 七道工序里走最多的是"不合格回炉"的红圈；网格无关性检查是收敛之外的另一道必需门。
• CFD 长在参数自由度与全场细节，风洞长在真实物理——两者是互补证据，不是替代关系。

数据生产出来了，敢不敢用？最后一节上质检体系。