4. 外流场计算域与边界条件

外流场仿真的可信度，很大程度上不取决于网格有多密、湍流模型有多高级，而取决于计算域是否足够大、边界条件是否与真实物理环境一致。计算域是人为截取的有限空间，边界条件则是在这个截断面上告诉求解器“外部世界是什么状态”。如果边界设置错误，即使内部流场算得再精细，气动系数也可能整体偏移几个百分点甚至出现非物理结果。

计算域尺寸：远场边界要“足够远”

外流场没有天然的外边界，必须人为设定远场。远场边界的核心要求是：不能对模型附近的流动产生明显干扰。对于常见的汽车、翼型、建筑绕流，经验做法是：

• 上游入口距模型前端约 3～5 倍特征长度；
• 下游出口距模型尾端约 10～15 倍特征长度，因为尾流恢复很慢；
• 侧向和顶部距模型约 5～10 倍特征长度；
• 若关注升力，顶部距离应更保守，避免出现“虚拟地面”效应。

这里的特征长度通常取模型迎风方向的最大尺寸，例如汽车取车长，翼型取弦长。计算域过小会人为加速绕流、抬高阻塞比，使阻力系数偏大；计算域过大则网格量激增、计算成本上升。因此需要在精度和成本之间取平衡，一般建议阻塞比控制在 1%～3% 以内。

入口、出口与壁面边界

入口边界通常采用速度入口，给定来流速度大小和方向。例如模拟 25 m/s 的侧风工况，可设置入口速度为 (25, 0, 0) m/s，同时给定湍流强度（如 1%～5%）和水力直径。入口湍流参数会直接影响边界层转捩和尾流发展，不能随意留默认值。

出口边界最常用压力出口，给定静压（通常为 0 Pa 表压），让流体自由流出。压力出口对回流有一定容忍度，但若出口离尾流太近，回流会污染内部解。若流动方向明确且尾流已充分发展，也可用出流边界，但它对回流不友好，需谨慎使用。

壁面边界分两类：无滑移壁面用于真实固体表面，速度为零，配合壁面函数或近壁网格处理边界层；滑移壁面用于模拟自由表面或对称面，法向速度为零、切向无剪切。对于风洞地面，若研究静止车辆，地面设为无滑移移动壁面或固定壁面，取决于是否模拟相对运动。

对称边界与周期边界

当几何和流动关于某个平面近似对称时，可用对称边界只计算一半模型，网格量和计算时间直接减半。对称边界要求法向速度为零、法向梯度为零，且不能有穿过对称面的质量交换。注意：只有流动本身对称时才成立，若存在侧风、横摆或非对称尾流，强行用对称边界会得到错误结果。

周期边界用于重复性结构，如叶栅、换热器翅片、阵列建筑。它把一对边界面“粘”在一起，使流出计算域的量从另一侧重新进入。平移周期用于直线重复，旋转周期用于圆周重复。使用周期边界的前提是几何和流动在周期方向上严格重复，且周期长度内包含完整流动特征。

边界条件对气动系数的影响

边界条件对气动系数的影响是系统性的。以下用一组教学示例说明（数值为合理假设，用于说明趋势，非实测结果）：

从表中可以看出，边界条件不是“填完就行”的选项，而是直接参与物理建模的组成部分。设置前应先明确：来流方向、参考坐标系、特征长度和参考面积，并在整个仿真中保持一致。气动系数通常按 Cd = 2Fd / (ρV²A) 计算，其中 A 为参考面积，V 为来流速度，ρ 为密度。若参考量前后不一致，系数对比就失去意义。

实操中的检查要点

完成边界设置后，建议按以下顺序检查：第一，确认计算域各方向尺寸满足经验倍数，阻塞比在合理范围；第二，确认入口速度方向与坐标系一致，参考面积、参考速度、密度与后处理一致；第三，确认出口没有明显回流，必要时把出口继续后移；第四，确认壁面类型与物理问题匹配，地面是否移动、是否考虑风洞壁干扰；第五，若使用对称或周期边界，确认流动本身满足对称或重复条件。

总之，计算域和边界条件共同定义了“外部世界如何与模型相互作用”。远场要足够远，入口要给出真实来流，出口要让尾流自由离开，壁面要符合真实表面状态，对称和周期边界只能在物理条件成立时使用。把这些关系理清，气动系数才有讨论的基础。