可压缩准一维管道无粘流动_欧拉方程_MacCormack解答.docVIP

可压缩准一维无粘管道流动 摘 要 本题利用一维欧拉方程求解可压缩一维无粘管道流动，并针对出口不同的条件，出口亚音速和出口超音速两种不同条件，将流道进行网格划分，利用MacCormack进行差分求解，得到管道内的总压、马赫数、总焓、内能的分布，并给出计算过程中残差收敛的过程。 关键词 准一维；欧拉方程；MacCormack 1 问题提出 如Figure 1所示，流动为流经变截面流道内、理想、不可压、定常、平面的流动。设进出口截面速度均匀分布。来流条件为，，。计算出口超音速和亚音速()时，这个准一维流动的压强和马赫数分布，并给出残差收敛过程。 Figure 1准一维管道示意图 本题的分析思路：首先，建立数学模型，之后运用边界条件定解。接下来将模型问题转化为数值求解问题，通过网格划分、差分离散，并运用matlab编程进行求解。 2 模型建立 选定控制体，根据积分型方程导出该一维流动的欧拉方程。控制体Figure 1中红框，、动量和方程如和 以上方程对于任意控制体均成立，因此可得到如下微分方程即 式中， 其中。对于此处的欧拉方程，有三个方程组，但是未知量分别为共5个，因此为了方程组的封闭需要补充两个方程。补充方程如下： 至此，已实现方程组的封闭，可以进行求解。 3 边界条件及初始条件 进口边界条件：进口为题中给定条件为，，。当地音速按公式，则进口速度为，总焓和总内能按照式(4)和式(5)进行求解。 出口边界条件：出口边界条件分别为超音速和亚音速。超音速则不给出口条件，在。 初始条件：其余节点初始条件假设与进口边界条件相同。 4 网格划分 由于本题所需求解的是一维管道流动问题，在同一截面位置，各个参数均匀分布，因此只需要沿着管道周向进行网格的划分即可，划分500个节点，空间步长为。 5 模型的差分离散 MacCormack格式广泛应用于求解流动方程，是一种显示预估/校正二步格式，对于求解一维欧拉方程其离散格式为： 两个方程中均加入了人工粘性进行修正，取人工粘性系数，本题计算时取。 在求解过程中，方程(7)求出的是在第时间层的临时估算值，而第时间层的最终值由校正方程(8)提供。需要注意到对空间进行差分离散时，预估步中采用的是向前差分，而在校正步中采用的是向后差分。每步迭代过程中，先将连续性方程、动量方程和能量方程进行预估迭代，之后将能量方程进行校正迭代，得到的新值，根据(3)式求出各点的压力新值，带入动量方程，并进行校正迭代，得到新的，再将其带入连续性方程，并进行校正迭代，得到新的，最终得到一组新的值。 为了使算法稳定，时间步长取法为： 定义收敛残差为相邻两次迭代过程所得的压力列向量之差的1范数： 计算时取为收敛条件。 6 计算结果及分析 (1) 出口为超音速条件时，计算结果如Figure 2出口超音速条件下计算结果所示 (a) 压力沿轴向分布 (b) 马赫数沿轴向分布 (c) 沿轴向分布 (d) 残差收敛曲线 Figure 2出口超音速条件下计算结果 从上述图中可以观察得知，当入口超音、出口超音时，流道内的压强以及马赫数变化较为平滑，符合超音在扩张管中降压的物理事实，且为等熵过程中波动在，可以认为）但同时也注意到，和马赫数的是因为出口给定条件和计算的出口不一致上来说就是出现了激波2000次之后开始呈减小趋势，经过4300步迭代，达到判定收敛标准。 根据等熵流动变截面管道流量公式结合条件，可求得理论解 、为总参数，流动中为常数 求出理论，Figure 3所示=10处的突变，其余完全重合 Figure 3 Ma理论和数值解对比，计算结果如所示： (a) 压力沿轴向分布 (b) 马赫数沿轴向分布 (c) 沿轴向分布 (d) 残差收敛曲线 图 1 出口为亚音速条件下计算结果 从上述结果可以看出，当入口超音、出口亚音时，流场内的压力、马赫数均在流场的中间某一位置发生了剧烈的突变现象，即流动过程中产生了激波，位置在左右。在激波之前等熵流动，经过激波之后，突增，之后仍等熵过程。后收敛。 激波，超音流在扩张管中持续加速，实现亚音出口，所以激波进口、出口下的必然产物之后的亚音流减速增压，亚音在扩压中的流动特征。调小出口速度，需要有足够的空间亚音流减速，那么激波必然提前，若增大出口速度，激波会。给出激波位置的理论计算方法。 为激波等熵方程， 、和分别为面积 马赫数马赫数和激波处面积同理，有，、、分别为面积马赫数马赫数