3. 控制方程与湍流建模选择
射流问题的CFD建模，第一步不是打开软件画网格，而是判断“该解什么方程、该用什么湍流模型”。射流的核心特征是：高速剪切层产生强湍流，同时伴随卷吸、掺混、势核衰减和下游自相似发展。若方程与模型选错，即使网格再密，也可能得到错误的势核长度、扩散角或中心线速度衰减率。本节围绕可压缩/不可压缩N-S方程、RANS/LES/DES混合方法，以及k-ε、k-ω SST、SA、Smagorinsky、WALE、动态Smagorinsky等模型在射流中的适用性展开，帮助读者建立“按射流物理选模型”的判断能力。
一、先判断可压缩还是不可压缩
射流是否可压缩，不能只看“气体”二字，而要看马赫数与密度变化。对于低速射流，例如出口速度30 m/s、环境温度300 K的空气射流，马赫数约0.09，密度变化小于5%，此时采用不可压缩N-S方程即可，压力以相对压力求解，能量方程可关闭或仅作被动标量处理。不可压缩求解的优点是压力-速度耦合稳定、时间步可较大，适合研究射流扩散角、势核长度和浓度衰减。
当出口马赫数超过约0.3，或存在强加热、强压缩、激波、欠膨胀/过膨胀等工况时，必须采用可压缩N-S方程，并保留能量方程与状态方程。可压缩射流中，密度不再是常数，湍流与压缩性相互作用，甚至出现激波-湍流干扰。此时若仍用不可压缩假设，会低估射流膨胀、错估激波位置和噪声源区。工程上常用判据是：Ma < 0.3且温差小，用不可压缩；Ma > 0.3或存在激波，用可压缩。
二、RANS、LES与DES的层级选择
RANS对全部湍流尺度做时间平均，计算量最小，适合参数扫描和工程初步设计，但射流中的大尺度涡卷起、剪切层不稳定和噪声源会被平均掉，因此对势核长度和远场扩散角的预测常偏保守。LES直接解析大尺度涡，只模化小尺度，能捕捉射流剪切层的涡配对和卷吸过程，适合研究掺混、燃烧和噪声，但近壁或近喷口区域网格量巨大。DES属于混合方法，在边界层内用RANS降低网格需求，在分离区和射流主体用LES解析大涡，适合高雷诺数射流与分离流动。
选择逻辑可以概括为：若只关心时均速度衰减和扩散角，RANS足够；若关心涡结构、掺混效率和噪声，选LES；若关心高雷诺数下分离与射流共存，选DES。需要注意的是，LES对网格和数值格式敏感，若网格不足以解析剪切层，结果可能比RANS更差。
三、RANS模型在射流中的差异
k-ε模型是射流计算中最常见的两方程模型，鲁棒性好、收敛快，但对圆形射流的扩散角预测偏小，常需要修正Cε1或采用RNG k-ε。标准k-ε在强旋流和强曲率下误差明显。k-ω SST模型结合了k-ω近壁优势和k-ε远场优势，对逆压梯度和剪切层更敏感，在射流中通常比标准k-ε更准确地预测势核长度和速度衰减，但入口湍流参数设置不当会导致剪切层过早或过晚发展。
SA模型是单方程模型，计算量小，在航空外流和边界层中表现优秀，但对自由射流的剪切层增长预测偏弱，因为其模型常数针对壁面剪切流标定。若射流以自由剪切为主，SA通常不是首选；若射流与壁面、分离区强耦合，SA可作为折中。
四、LES亚格子模型的选择
Smagorinsky模型是最经典的亚格子模型，形式简单，但耗散偏大，在射流剪切层中会抑制小尺度涡，导致掺混不足。其系数Cs固定，近壁处需要阻尼函数。WALE模型在近壁处具有正确的渐近行为，不需要壁面阻尼，对射流近喷口区域和转捩预测更友好。动态Smagorinsky模型通过Germano恒等式动态确定系数，能自适应层流和湍流区域，在射流中可减少近场过度耗散，但计算量增加，且需要二次过滤，对数值实现要求较高。
实际选择时，若射流以时均特性为主，优先k-ω SST；若关注涡结构和噪声，优先LES，亚格子模型可选WALE或动态Smagorinsky；若高雷诺数且计算资源有限，选DES，近壁RANS用SST，分离区LES用WALE。入口湍流参数必须与模型匹配：k-ε需要给定k和ε，SST需要给定k和ω，LES需要给定湍流强度和长度尺度，否则射流剪切层发展会失真。
五、一个可复现的示例设置
以圆形空气射流为例：喷口直径D=0.05 m，出口速度Uj=50 m/s，环境静止，温度300 K。马赫数约0.15，选不可压缩N-S。若关注势核长度和中心线衰减，选k-ω SST，入口湍流强度5%，湍流长度尺度0.07D。若关注涡卷起和噪声，选LES，亚格子用WALE，喷口下游2D内网格尺度约D/40，时间步满足CFL<0.5。计算后观察中心线速度Uc/Uj随x/D的变化：SST通常预测势核约6D，LES约5D，若差异过大，检查入口湍流参数和网格分辨率。若出现非物理振荡，检查对流格式是否过耗散或过色散。
模型选择没有绝对最优，只有与射流物理和计算目标匹配。先明确马赫数、关注尺度和计算资源，再按“方程—方法—模型—入口参数”的顺序逐层确定，才能让射流CFD结果既有物理意义又可复现。