4. 射流入口与边界条件建模

射流CFD计算的精度，往往在入口边界条件设定时就已经决定了大半。喷嘴出口的速度型、湍流参数、伴随流与静压条件、远场处理方式，共同决定了剪切层初期的增长率、势核长度以及下游湍流统计量的可信度。本节围绕这些入口与边界条件，说明其物理含义、建模方式与实现要点。

喷嘴出口速度型的选择

喷嘴出口速度型是射流入口最核心的输入。常见有三类：顶帽型、抛物线型和湍流边界层型。顶帽型假设出口截面上轴向速度均匀，适用于长径比较大、出口已充分发展的收缩喷嘴，或作为简化对比算例。抛物线型对应充分发展的层流管流，速度分布为 u(r)=u_max(1-(r/R)^2)，中心速度约为截面平均速度的2倍。湍流边界层型则更接近真实喷嘴：壁面附近速度迅速下降，中心区域近似平坦，可用1/7次幂律 u/u_c=(1-r/R)^(1/7) 近似，或直接由上游管流计算结果映射得到。

选择速度型时需注意：若喷嘴出口存在明显收缩，实际速度型介于顶帽与湍流边界层之间；若用抛物线型模拟高雷诺数射流，会人为放大中心速度、削弱剪切层，导致势核偏长。工程上建议优先采用湍流边界层型或由上游管流映射，仅在验证性算例中使用顶帽型。

湍流参数的指定

入口湍流参数通常以湍流强度 I 与湍流长度尺度 L（或水力直径 D_h）给出。对圆喷嘴，D_h 取喷嘴直径 D；对矩形喷嘴，D_h=4A/P。湍流强度可按经验估计：充分发展管流出口 I≈5%~10%，收缩喷嘴出口可低至 1%~3%。湍动能 k=1.5(u_avg I)^2，比耗散率 ω=k^0.5/(C_μ^0.25 L)，其中 C_μ=0.09，L=0.07D_h 为常用估计。若采用SST k-ω 模型，还需给出比耗散率；若采用Realizable k-ε，则给出耗散率 ε=C_μ^0.75 k^1.5/L。

需要强调的是，入口湍流参数对射流近场影响显著，对远场逐渐减弱。若研究重点是势核长度和近场剪切层，应尽量使用实测或上游计算得到的 k、ω 剖面，而非单一均匀值。

合成湍流方法

当入口需要非均匀、含脉动的湍流场时，均匀湍流参数不足以描述真实入口。合成湍流方法可在入口生成满足指定湍流统计的脉动速度场，常用两类：涡方法和数字滤波法。

涡方法（如Lamb-Oseen涡叠加）在入口平面布置若干随机涡，每个涡具有随机位置、环量和半径，通过叠加诱导速度生成脉动。其优点是实现简单、可指定湍流长度尺度；缺点是脉动场不完全满足连续性方程，需配合散度修正。数字滤波法则通过对随机白噪声场做空间相关滤波，使生成的脉动速度具有指定两点相关函数，能较好保持湍流统计特性，但计算量略大，且需要给定相关长度。

实际使用中，若仅需提升剪切层初期的湍流发展速度，涡方法已足够；若需精确匹配入口湍流谱或相关函数，应选择数字滤波法。无论哪种方法，都应在入口下游留出足够长的过渡区，避免合成湍流未发展就进入关键测量区。

伴随流与静压边界

伴随流（co-flow）是喷嘴外部同向或横向的环境流动。同向伴随流会降低射流与环境的速度差，从而减弱剪切层增长率、延长势核；横向伴随流则使射流发生偏转，形成弯曲射流。建模时，伴随流速度通常在入口边界的外环区域给定，或通过独立的伴随流入口引入。若伴随流速度远小于射流速度（如 U_co/U_j < 0.1），可近似忽略；若达到 0.2~0.5，则必须显式建模。

静压边界用于出口和远场。对亚声速射流，出口通常设为压力出口，给定静压等于环境压力 p_amb，并允许回流。若出口出现持续回流，说明计算域下游太短或出口位置处于回流区，应延长计算域或改用压力远场边界。对超声速射流，出口需采用压力远场或超声速出口，避免非物理反射。

远场边界处理

远场边界应布置在射流影响区之外。经验上，径向远场距轴线不小于 20D，轴向出口距喷嘴不小于 50D~100D（视研究目标而定）。远场类型优先选择压力远场或无反射边界，给定环境静压、温度和来流速度。若计算域受限，可在远场外侧加阻尼区或采用对流出口条件，减少压力波反射对射流近场的影响。

参数设置示例

以下给出一个圆喷嘴射流入口的典型参数设置示例（示例参数，非实测结果）：喷嘴直径 D=0.02 m，射流出口平均速度 U_j=30 m/s，环境空气密度 ρ=1.225 kg/m³，动力粘度 μ=1.79e-5 Pa·s。则 Re=ρU_jD/μ≈4.1e4，为湍流射流。取湍流强度 I=5%，水力直径 D_h=D=0.02 m，则 k=1.5(30×0.05)²≈3.375 m²/s²，L=0.07×0.02=0.0014 m，ω=k^0.5/(0.09^0.25×0.0014)≈1.84/(0.5477×0.0014)≈2400 s⁻¹。若采用Realizable k-ε，则 ε=0.09^0.75×3.375^1.5/0.0014≈0.164×6.20/0.0014≈726 m²/s³。这些值可直接填入入口边界面板。

# 入口边界示例（圆喷嘴，D=0.02 m）
velocity_type   = turbulent_boundary_layer   # 1/7幂律
U_jet           = 30 m/s
turb_intensity  = 5%
D_h             = 0.02 m
k               = 3.375 m2/s2
omega           = 2400 1/s        # SST k-omega
epsilon         = 726 m2/s3       # Realizable k-epsilon

# 伴随流
U_coflow        = 3 m/s           # 同向，U_co/U_j = 0.1
coflow_region   = annular

# 出口/远场
outlet_type     = pressure_outlet
p_gauge         = 0 Pa
farfield_type   = pressure_farfield
farfield_radius = 20D
domain_length   = 80D

设置完成后，建议先做网格无关性验证：在相同入口条件下，比较两套网格的轴线速度衰减曲线。若势核长度和衰减率差异小于5%，可认为网格足够。随后再开展参数研究或与实验对比。