手机版

ANSYS气体仿真

发布时间：2026-09-02 来源：科研学术网

小中大

字号：

一、项目背景与科学问题

可压缩气体在拉伐尔（收敛-扩张）喷管中的流动是空气动力学与推进技术的基本问题。当气体在收敛段加速至喉部达到声速（M=1）后，在扩张段继续加速至超声速，实现压力能向动能的充分转化。喷管是火箭发动机、超声速风洞、蒸汽轮机与射流切割设备的核心部件，其流量特性与沿程参数分布直接影响推进性能。

本项目的科学问题是：在给定总压、总温下，气体流经收敛-扩张喷管时，马赫数沿轴向如何变化？喉部是否达到声速？静压与静温沿轴向如何下降？一维等熵理论能否准确描述喷管内的流动？利用 ANSYS CFD 对喷管内可压缩流动进行仿真可以定量回答。

二、计算方法与理论背景

2.1 一维等熵流动

对可压缩无黏等熵流动，面积-马赫数关系为 (A/A*)² = (1/M²)·[(2/(γ+1))·(1+(γ−1)M²/2)]^((γ+1)/(γ−1))，其中 A* 为喉部面积，γ 为比热比（空气取 1.4）。对给定面积比，亚声速与超声速两支分别对应 M<1 与 M>1 的解：收敛段取亚声速支，扩张段取超声速支，喉部 M=1。

2.2 等熵参数关系

静压与总压之比 p/p₀ = (1+(γ−1)M²/2)^(−γ/(γ−1))，静温与总温之比 T/T₀ = (1+(γ−1)M²/2)^(−1)。随马赫数增大，静压与静温单调下降：气体膨胀加速，内能转化为动能。

2.3 数值方法

CFD 求解可压缩 N-S 方程（密度基求解器），采用二阶迎风格式离散对流项，标准 k-ε 或 SST 湍流模型处理边界层。压力远场或总压入口、压力出口作为边界条件。

三、计算设置与模型构建

软件：ANSYS Fluent，二维轴对称密度基稳态求解器。
几何：收敛-扩张喷管，喉部位于中部（x/L=0.5），喉部面积最小，进出口面积比约 1:3。
气体：空气，理想气体，γ=1.4。
边界：入口总压 p₀=5 bar、总温 T₀=300 K；出口静压设置使喷管处于设计工况（扩张段全超声速）。
湍流模型：SST k-ω，近壁 y⁺≈1。
网格：结构化网格，喉部加密，总网格约 8 万。
收敛：残差低于 10⁻⁶，出口流量守恒；监控喉部马赫数趋于 1。
输出：沿轴向的马赫数、静压比 p/p₀、静温比 T/T₀ 分布，与一维等熵理论解对比。

四、结果分析与案例图解读

图 14 的左面板（a）给出马赫数沿喷管轴向的分布。气体在收敛段从约 0.3 缓慢加速，接近喉部时快速上升，在喉部（x/L=0.5）达到 M=1（图中灰色虚线标注），随后在扩张段继续加速至出口约 M=1.8。这一”亚声速收敛-声速喉部-超声速扩张”的经典形态与一维等熵面积-马赫数关系一致。喉部恰好达到声速是壅塞（choking）条件的体现：在给定总压下，喷管流量达到最大，下游扰动无法向上游传播。出口 M=1.8 表明气流在扩张段完成了充分的膨胀加速，对应设计工况。

图 14 的右面板（b）给出静压比 p/p₀ 与静温比 T/T₀ 沿轴向的分布。随马赫数增大，静压与静温单调下降：入口处 p/p₀≈1，喉部降至约 0.53，出口降至约 0.16；静温比相应从 1 降至出口约 0.61。压降与温降体现了气体绝热膨胀过程中压力能与内能向动能的转化。例如出口静温约 183 K，对应气体在喷管中被显著冷却，这是高速气流的典型特征。数值结果（曲线）与一维等熵理论解的吻合验证了仿真模型的正确性，也说明对于无激波的设计工况，喷管流动可用简单的一维理论良好描述。

两面板结合，左面板展示了马赫数的加速过程与壅塞特性，右面板揭示了静压静温的膨胀下降规律，共同完成收敛-扩张喷管气体仿真的核心分析，为喷管设计与推进性能评估提供了可靠依据。