求解PNS方程的流向通量分裂法论文.pdfVIP

求解PNS方程的流向通量分裂法 求解PNS方程的流向通量分裂法 贺旭照乐嘉陵张若凌杨顺华赵慧勇 中国空气动力研究与发展中心吸气式高超声速技术实验室，四川，绵阳，621000 摘要：依据超音速流场中物理信息传播的规律以及超音速流场中求解点信息依赖域的特点，在空间推进求解 PNS方程中，沿推进求解方向(流动主流方向)对对流项进行正负通量分裂。对负项通量进行合理的改造， 就可以使控制方程沿主流方向具有抛物一双曲性质。结合伪时间迭代法求解PNS方程的思想，构造了适用于现 有求解NS方程数值方法和计算软件的求解PNS方程的方法。对采用这种方法设计的程序进行了数值验证。这 种方法物理概念清晰，数学处理简单，易于从现有的求解Ns方程的计算程序改造实现。所得到的数值模拟结 果以及计算速度和对计算机内存的低消耗都达到了令人满意的程度。 0引言 PNS方程适用于求解超音速流动占主导的沿主流方向无大分离流动的流场，新一代的高超 音速飞行器其流场的很大部分都具有这样的流动特性。采用PNS方程求解这类流场时，具有 显见的优越性。在不失求解流场准确性和精度的情况下，空间推进求解PNS方程方法较全场 迭代求解NS方程方法，具有快速和消耗内存少的优点。目前国外把PNS和NS方法相结合， 完成了对新一代超音速飞行器分区求解，PNS方法在工程设计中得到了极大的应用。 为了空间推进求解，必须对原NS方程的流向对流项和粘性耗散项做一些处理，以使控制 方程在流向具有抛物一双曲性质。以往的对流向对流项的处理都是分裂动量方程中的压力项。 在流场的超音速区域压力项保留不变，在流场的亚音速区域(边界层内)对压力项进行分裂， 部分的舍去压力项或对其作相应的修改u、2、33，以减小在亚音速区域内逆向压力梯度的影响， 使控制方程在主流方向上达到抛物一双曲性质。 本文采用流向通量分裂的思想，对流向的对流项进行通量分裂。若以求解点及推进面为 参照，以主流方向为流场的X轴正方向(下图)，正项通量的信息依赖域位于推进面的左端， 负项通量的信息依赖域位于推进面的右端。在超音速区域，流向对流项的负项通量为零，求 解点的信息完全依赖于推进面左端的流场信息，这时在流向上控制方程具有抛物一双曲性质； 在亚音速区域，求解面的信息依赖域 位于推进面的左右两侧，控制方程沿 流向为椭圆型，若对流向对流项的负 项通量作适当的处理(舍去或部分的 鑫 人工模拟)，那么改造过的控制方程在 喜 亚音速流动部分也会呈现抛物一双曲 藉 性质。对粘性耗散项通常舍去其沿流 丧 向的倒数项，这时粘性耗散项的信息 整 依赖域也位于推进面以左。这样流动 的控制方程沿流向就具有整体的抛物 一双曲性状，空间推进求解超音速占主 导的沿流向无大分离流动的流场变得 可行。 目前有两种方法对PNS方程空间推进求解。一种是空间推进求解定场PNS方程，这种方 法要求空间推进的步长小，求解同一问题要更多地推进步数，但每个推进面上勿须迭代。另 Schemefor PNS SolvingEquations)，这种方法是在每个推进面上迭代求解非定常的PNS方程，直至在 此推进面上收敛，然后进入下一个推进面。这种方法对推进步长有所放宽，求解相同的问题 34l 第十一届全国激波与激波管学术会议 只需要较少的推进步数，但在每一个推进面上需要多次迭代。 本文借鉴伪时间法推进求解PNS方程的思想，结合本文提出的流向通量分裂思想，针对 现有的求解NS方程的数值方法和求解NS方程的程序，构造了使用现有NS程序的求解P