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

可压缩准一维无粘管道流动 摘 要 本题利用一维欧拉方程求解可压缩一维无粘管道流动，并针对出口不同的条件，出口亚音速和出口超音速两种不同条件，将流道进行网格划分，利用MacCormack进行差分求解，得到管道内的总压、马赫数、总焓、内能的分布，并给出计算过程中残差收敛的过程。 关键词 准一维；欧拉方程；MacCormack 1 问题提出 如Figure 1所示，流动为流经变截面流道内、理想、不可压、定常、平面的流动。设进出口截面速度均匀分布。来流条件为，，。计算出口超音速和亚音速()时，这个准一维流动的压强和马赫数分布，并给出残差收敛过程。 Figure 1准一维管道示意图 本题的分析思路：首先，建立数学模型，之后运用边界条件定解。接下来将模型问题转化为数值求解问题，通过网格划分、差分离散，并运用matlab编程进行求解。 2 模型建立 选定控制体，根据积分型方程导出该一维流动的欧拉方程。控制体Figure 1中红框，、动量和方程如和 以上方程对于任意控制体均成立，因此可得到如下微分方程即 式中， 其中。对于此处的欧拉方程，有三个方程组，但是未知量分别为共5个，因此为了方程组的封闭需要补充两个方程。补充方程如下： 至此，已实现方程组的封闭，可以进行求解。 3 边界条件及初始条件 进口边界条件：进口为题中给