[论文解读] Global Subsonic and Subsonic-Sonic Flows through Infinitely Long Axially Symmetric Nozzles

Chunjing Xie, Zhouping Xin|ArXiv.org|Jul 19, 2009

Navier-Stokes equation solutions参考文献 4被引用 3

一句话总结

该论文利用变分法、椭圆估计和补偿紧致性方法，建立了无限长轴对称喷管中全局亚音速流及亚音速-声速流的存在性。当进气质量通量低于临界值时，存在均匀亚音速流，且在无穷远处趋于均匀流，流场角度一致有界；当通量达到临界值时，通过补偿紧致性方法全局存在弱亚音速-声速流。

ABSTRACT

In this paper, we establish existence of global subsonic and subsonic-sonic flows through infinitely long axially symmetric nozzles by combining variational method, various elliptic estimates and a compensated compactness method. More precisely, it is shown that there exist global subsonic flows in nozzles for incoming mass flux less than a critical value; moreover, uniformly subsonic flows always approach to uniform flows at far fields when nozzle boundaries tend to be flat at far fields, and flow angles for axially symmetric flows are uniformly bounded away from $π/2$; finally, when the incoming mass flux tends to the critical value, subsonic-sonic flows exist globally in nozzles in the weak sense by using angle estimate in conjunction with a compensated compactness framework.

研究动机与目标

• 将亚音速流与亚音速-声速流的存在性理论从二维推广至三维轴对称喷管。
• 解决尽管已知具有单调性性质，但三维轴对称亚音速流仍缺乏存在性结果的问题。
• 利用补偿紧致性方法，在质量通量达到临界值时建立弱亚音速-声速流的全局存在性。
• 证明当喷管边界趋于平坦时，亚音速流的流场角度保持一致有界，并在无穷远处收敛于均匀流。
• 通过流函数形式与椭圆正则性，为长喷管中三维可压缩欧拉流提供严格的理论框架。

提出的方法

• 将三维等熵可压缩欧拉方程用速度势Φ表示，将问题约化为非线性椭圆方程：div(g(|∇Φ|²)∇Φ) = 0。
• 施加边界条件：在喷管壁面上∂Φ/∂n = 0（不可渗透边界），并通过横截面∫_S g(|∇Φ|²) ∂Φ/∂l dS = m₀ 固定质量通量。
• 应用变分法构造质量通量低于临界值时的亚音速流近似解。
• 利用加权Sobolev与Hölder估计控制速度势及其梯度的正则性与一致有界性。
• 采用补偿紧致性框架，在质量通量趋近临界值时取极限，确保弱亚音速-声速解的存在性。
• 引入截断函数与测试函数以处理非线性与边界行为，特别是在声速状态附近。

实验结果

研究问题

• RQ1 在质量通量低于临界阈值时，能否在三维轴对称喷管中构造出全局亚音速流？
• RQ2 当喷管边界趋于平坦时，此类喷管中的均匀亚音速流是否在无穷远处收敛于均匀流？
• RQ3 轴对称亚音速流中的流场角度是否一致有界远离π/2？
• RQ4 当进气质量通量达到临界值时，能否在同一设定下全局构造出弱亚音速-声速流？
• RQ5 补偿紧致性方法在处理三维流中声速状态处的退化问题时起到何种作用？

主要发现

• 在喷管边界为C^{1,α}且远离零的假设下，对于所有低于临界值的进气质量通量，全局光滑亚音速流存在。
• 当喷管边界在无穷远处趋于平坦时，均匀亚音速流在无穷远处收敛于均匀流，且流场角度一致有界远离π/2。
• 当进气质量通量趋近临界值时，弱亚音速-声速解在喷管中全局存在，其意义为分布意义下的解。
• 补偿紧致性方法结合角度估计与椭圆正则性，确保近似解在临界通量下收敛于弱解。
• 密度ρ通过ρ = g(q²)表示为q²的递减函数，(ρq)²表示为G(q²)，且H(1) = 1确保声速状态对应ρ = 1。
• 弱解在发散测度场意义下满足连续性方程与边界条件，符合Anzellotti的定义，确认了物理解释的一致性。