4. CFD求解内流道压降的总体流程
内流道压降CFD分析不是“画完网格点一下求解”这么简单。它是一条从CAD几何到最终压降数值的完整数据链，任何一环处理不当，都会让最后的Δp失真。本节把这条链拆成五个阶段：几何清理、网格划分、物理模型选择、边界条件设定、后处理提取压降，并说明每个阶段真正要做的判断和操作。
一、几何清理：把“能画出来”变成“能算出来”
CAD模型面向制造，CFD模型面向流动。导入的几何常带有对流动无影响、却会破坏网格的细节：螺栓孔、倒角、铭牌凹槽、装配缝隙、微小台阶。这些特征尺寸可能只有0.2 mm，而流道主尺度是20 mm，若保留，局部网格会被迫加密到极小尺度，网格量暴涨且质量恶化。
清理的核心判断是：该特征是否显著改变流动分离、涡系或流通面积。以某液冷板内流道为例，流道截面为8 mm×6 mm，进出口为Φ10 mm圆管。清理时删除M4安装螺纹孔、0.5 mm工艺倒角、表面丝印凹槽；保留流道内所有转折、分流脊和扰流柱，因为它们直接决定压降。清理后还需做“水密性”检查：用填充孔洞、缝合自由边、修复相交面，确保几何是单一封闭流体域。若存在微小缝隙，网格会在此处产生负体积或泄漏，求解直接发散。
一个实用做法是：在CAD中先抽取内流道流体域，而不是把固体零件直接丢给网格工具。抽取后测量关键截面面积，与设计值对比，偏差超过2%就说明清理过程中误删或误补了特征。
二、网格划分：在精度与规模之间找平衡
内流道压降对近壁面速度梯度敏感，因此边界层网格是重点。以雷诺数Re≈8000的湍流工况为例，若采用k-ε模型配合壁面函数，第一层网格中心到壁面的无量纲距离y+宜落在30～300；若采用SST k-ω模型做低雷诺数近壁解析，则要求y+≈1。两者网格策略完全不同，必须在划网格前就确定物理模型方向。
体网格优先选多面体或六面体核心。多面体网格在复杂流道中收敛快、对偏斜容忍度高；六面体在规则直道中精度和数量更优。边界层通常设5～10层，增长率1.2～1.3，总厚度覆盖边界层厚度的1.2倍以上。全局尺寸可先按流道最小特征尺寸的1/5起步，再通过网格无关性验证调整。
网格质量看三个指标：最小正交质量（建议＞0.15）、最大长宽比（边界层内可放宽，但核心区宜＜100）、负体积数量（必须为0）。若出现负体积，优先回到几何清理阶段查缝隙和尖角，而不是反复调网格尺寸。
三、物理模型选择：定常还是瞬态，层流还是湍流
先算雷诺数。以水为介质、入口流速1.5 m/s、水力直径6 mm、运动黏度1×10⁻⁶ m²/s估算，Re≈9000，属于湍流。若实际工况是低流速0.05 m/s，Re≈300，则为层流，此时用层流模型即可，湍流模型反而会引入不必要的耗散。
湍流模型选择上，内流道压降计算常用SST k-ω或Realizable k-ε。SST k-ω对逆压梯度和分离预测更好，适合有分流、转弯、扰流柱的流道；Realizable k-ε计算稳定、资源消耗低，适合长直道为主的流道。若流道内存在明显非定常涡脱落，例如大空间突扩后接弯头，应做瞬态计算，时间步长按最小涡时间尺度的1/10取，例如0.001 s，并监测出口压力随时间波动，取统计平均值作为压降。
介质属性必须按实际工况输入：水在60℃时密度约983 kg/m³、黏度约0.47 mPa·s，与20℃的998 kg/m³、1.0 mPa·s差别明显，直接影响压降约10%～15%。不可默认用常温物性代替。
四、边界条件设定：压降的“输入”和“出口”
内流道压降通常用质量流量入口+压力出口组合，因为流量是实际运行中最容易确定的量。入口给定质量流量，例如0.12 kg/s；出口给定0 Pa表压。壁面采用无滑移、绝热条件。若流道有外部散热，则壁面改为恒热流或对流换热，但压降计算中热边界对动量影响较小，可先按绝热简化。
收敛判断不能只看残差。残差降到1×10⁻⁴以下只是基本要求，更关键的是监测进出口总压差是否稳定。当进出口压差在连续200步内波动小于0.5%，且入口质量流量不平衡小于0.1%，才可认为收敛。若压差持续缓慢漂移，说明网格或模型仍有问题。
五、后处理提取压降：算完了，怎么读
压降定义为入口截面总压平均值减去出口截面总压平均值：Δp = Pin,total − Pout,total。总压包含静压和动压，若进出口截面积相同，动压项近似抵消，也可用静压差近似；但若进出口面积不同，必须用总压，否则会漏掉动压变化带来的压降贡献。
提取时不要只取一个点。应在入口和出口截面上做面积加权平均，避免局部回流区拉偏结果。同时检查进出口质量流量是否守恒，若偏差超过0.5%，说明收敛不充分或存在泄漏。
最后做网格无关性验证：用粗、中、细三套网格分别计算，若中网格与细网格的压降偏差小于2%，可认为中网格已足够。以示例工况为例，粗网格Δp=18.6 kPa，中网格Δp=17.2 kPa，细网格Δp=17.0 kPa，中细偏差1.2%，选中网格即可。若三套网格结果持续单调变化，说明还没进入无关区，需继续加密。
把五个阶段串起来看，压降CFD求解的本质是：用干净的流体域、合适的网格和物理模型，在正确的边界条件下求解动量方程，最后从收敛的流场中提取进出口总压差。每一步的判断依据都指向同一个目标——让Δp可信、可复现、可用于设计决策。