非等熵势函数及其在叶轮机械内跨声速流动数值计算的应用

在对近年来有了很大发展的各种非等熵势函数方法进行理论分析的基础上,本文指出,从散度型动量方程直接求解密度并与势函数方程迭代的方法较为完善。这种方法既保持了势函数计算简便的优点,又扩大了势函数的适用范围,可以在跨声速计算中发挥很大的作用。     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

联合使用激波捕获—拟合法求解激波间断的跨声速流场

吸取激波捕获和激波拟合两种方法各自的优点,发展了一种求解有激波的跨声速流场的方法,并对回转叶栅中定常跨声速流场做了计算。计算中,通过捕获法确定初始激波位置,然后经反复修正,拟合出确定的流场通道激波。计算表明,这种方法可自动获得清晰的激波,流场中气流各参数分布合理。这一方法可适用于复杂边界和不同进口M数,计算时间仅比势函数方法多一倍左右。     

考虑NND格式修正的非定常Taylor—Galerkin有限元法

以Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ法为基础，吸收了差分格式中ＮＮＤ格式的思想，并在时间方向上采用四阶Ｒｕｎｇｅ－Ｋｕｔｔａ方法，构成了适用于定常和非定常跨声速流动计算的有限元算法。利用文的方法，分别计算了跨声速叶栅流动和工程中一类复杂非定常轴对称跨声速流动。计算表明，方法克服了Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ法的局限性，对激波及流场的分辨率较好，可用于跨声速流动计算，并可用于非定常流场的计算。     

修正的跨声速小扰动方程以及与其相容的激波方程

本文给出了与修正的跨声速小扰动(MSD)方程相容的激波方程,计算证明它的极曲线与精确的激波极曲线十分接近。     

跨声速轴流压气机特性计算用激波损失模型

基于实验结果,发展了一种适用于跨声速叶型的新激波损失模型,该模型考虑了压气机工作状态从堵塞点向喘振点变化时激波系结构的变化。将新激波损失模型应用于某两级跨声速轴流压气机特性预测中,计算与实验结果的比较表明,在跨声速轴流压气机非设计点性能预测中,该模型能够给出更为准确的结果,对于效率曲线的改进尤为明显,并且激波损失系数在全工况下的变化规律更加合理,这表明新激波损失模型能够准确地反映出跨声速压气机内部的流动机理和激波损失大小。     

一种求解核函数方程的新方法—辅助核法

本文提出了求解核函数方程的辅助核法。此法克服了偶极子栅格法对网格分布的限制,具有更大的适用范围。数值计算表明,此法对求解绕翼型的跨声速激波流动是很有利的。     

时间推进积分方程法求解非定常跨声速流动

本文提出了一种求解非定常跨声速流动的新方法——时间推进积分方程法,此法克服了时间线化积分方程法的限制,能较好地模拟激波的运动。本文首先用一维波(?)问题——模型问题阐明此法的基本思想,然后将它应用于二维非定常跨声速流动中。本文还首次引入拟速度位的概念,使时间推进积分方程式得到简化,尾涡条件和Kutta条件更易处理。数值计算表明时间推进积分方程法是合理可靠的。     

超燃冲压发动机一维稳态跨声速流动奇异初值问题的一个解法

超燃冲压发动机燃烧室一维稳态跨声速流动方程在临界声速点存在奇异初值问题，现有的基于L’Hospital法则的求解方法在原理上存在较大的初值误差，影响一维稳态跨声速计算的精度。为此，本文提出了一种基于变量代换的改进算法，通过定义新的流动变量W=-Ma^2＋2Ma，构造出了非奇异的一维稳态跨声速流动方程。消除了微分方程的奇异性，有效的解决了一维稳态跨声速流动计算过程中存在的奇异初值问题。     

A Structured Mesh Euler and Interactive Boundary Layer Method for Wing/Body Configurations

To compute transonic flows over a complex 3D aircraft configuration, a viscous/inviscid interaction method is developed by coupling an integral boundary-layer solver with an Eluer solver in a ＂semi-inverse＂ manner. For the turbulent boundary-layer, an integral method using Green＇s lag equation is coupled with the outer inviscid flow. A blowing velocity approach is used to simulate the displacement effects of the boundary layer. To predict the aerodynamic drag, it is developed a numerical technique called far-field method that is based on the momentum theorem, in which the total drag is divided into three component drags, i.e. viscous, induced and wave-formed. Consequently, it can provide more physical insight into the drag sources than the often-used surface integral technique. The drag decomposition can be achieved with help of the second law of thermodynamics, which implies that entropy increases and total pressure decreases only across shock wave along a streamline of an inviscid non-isentropic flow. This method has been applied to the DLR-F4 wing/body configuration showing results in good agreement with the wind tunnel data.     

Stagnation temperature effect on the conical shock with application for air

The aim of this work is to realize a new numerical program based on the development of a mathematical model allowing determining the parameters of the supersonic flow through a conical shock under hypothesis at high temperature, in the context of correcting the perfect gas model. In this case, the specific heat at constant pressure does not remain constant and varies with the increase of temperature. The stagnation temperature becomes an important parameter in the calculation. The mathematical model is presented by the numerical resolution of a system of first-order nonlinear differential equations with three coupled unknowns for initial conditions. The numerical resolution is made by adapting the higher order Runge Kutta method. The parameters through the conical shock can be determined by considering a new model of an oblique shock at high temperature. All isentropic parameters of after the shock flow depend on the deviation of the flow from the transverse direction. The comparison of the results is done with the perfect gas model for low stagnation temperatures, upstream Mach number and cone deviation angle. A calculation of the error is made between our high temperature model and the perfect gas model. The application is made for air.     

湿压缩对压气机转子失稳边界工况性能的影响

针对某型跨声速压气机转子的失速边界工况,采用数值模拟方法系统研究了不同加湿条件对其性能的影响。获得了水滴粒径、速度、温度和喷水量与转子进气流量、压比、效率、比压缩功等性能参数之间的关系曲线。结果表明:在失速边界上对跨声速压气机进气加湿,可以增加空气流量,提高总压比和效率,降低比压缩功。通过对流场的分析,发现进气加湿后,向低能区注入了高动量的水滴颗粒,削弱叶顶阻塞团,将主流与泄漏流交界面和激波位置向下游推进,减小激波强度,降低叶顶载荷,增大主流区载荷,改善了叶顶流场。提出了以水滴的面积流量定量评价不同进气加湿条件对叶顶流场的作用效果。      

基于动网格的涡轮流场仿真方法及特性分析

为提高涡轮流场仿真计算的收敛性和结果的精度,基于动网格的流场仿真技术,探讨了涡轮瞬态流场仿真的 2步方法——先稳态后瞬态计算。计算结果表明:在转子—静子间滑动界面上出现了计算方程守恒引起的流场参数阶跃,流场参数变化趋势与一维等熵流动理论分析结果一致。     

跨声速非匀熵全速位方程的有限元数值解(英文)

本文改进了传统全速位方程的有限元数值解,既保留了全速位方程简单的优点又考虑了出现激波后流场不是匀熵的特点,对密度和后缘库塔条件进行修正,改善了计算结果。文中给出翼型和机翼的算例,其结果与传统全速位方程的结果相比,激波的位置偏往上游,强度较弱,这个趋势使结果较接近于欧拉方程的解。     

渗透壁面竖直平板Blasius流熵产特性

对渗透壁面竖直平板Blasius流层流边界层流动熵产及传热熵产与总熵产的比值(Bejan数)进行了深入研究,将耦合的动量与能量偏微分方程组通过相似变换转换为非线性常微分方程组,用Runge-Kutta法求解获得了无量纲流函数与无量纲温度的相似解,进而研究了无量纲参数对熵产及Bejan数的影响.结果表明:吸入时最大熵产位于壁面处,吸入速度越大最大熵产越大,熵产随相似变量的增加单调下降;喷注速度越大壁面处熵产越小;近壁面处熵产随变黏度参数增大而增大,远离壁面处随喷注/吸入参数或变黏度参数的增大而减小;远离壁面处Bejan数随喷注/吸入参数或变黏度参数的增大而增大,近壁面处与此相反;熵产和Bejan数随着辐射参数和Reynolds数增大而减小,随着Biot数增大而增大;Biot数在0~2范围内对传热熵产的影响很剧烈,Biot数较小时总熵产由流动熵产控制,Biot数较大时总熵产由传热熵产控制;熵产随滑移参数的减小而增大,Bejan数随滑移参数的增大而增大;Grashof数对熵产的影响较大,在较大的Grashof数条件下熵产随相似变量的变化较剧烈;滑移参数越大,流动熵产越小,Bejan数越大.      

超音速气流绕局部透气凹角流动的数值模拟

以具有压力分裂形式的简化Ｎ Ｓ方程为控制方程,数值模拟了超音速来流条件下的激波 边界层干扰被动控制（ｐａｓｓｉｖｅｃｏｎｔｒｏｌｏｆｓｈｏｃｋ ｂｏｕｎｄａｒｙｌａｙｅｒｉｎｔｅｒａｃｔｉｏｎ）。模拟是以预先给定激波前吹气和激波后吸气的流量来实现的。为了定性地确定吹气或吸气对激波 边界层干扰的影响,首先计算了单独吹气和单独吸气两种情况。数值计算时采用了多重扫描法对控制方程差分离散,以反映亚音速区压力对流场的椭圆性影响。     

变截面管道中激波诱导的两相流动

本文研究扩张型与收缩型变截面管道中激波诱导的非定常两相流动。对于稀相气固悬浮体采用了双连续介质模型,对于两相流动按照准一维近似处理。控制方程则是利用算子分裂技术和二阶GRP方法数值求解。文中给出了波后气固两相的流动结构,并讨论了粒子对两相流场的影响。