一种求解涡轮平面叶栅的跨音松驰法

一、引　　言跨音松驰法是七十年代初发展的一种新方法。文献 [1 ]是一篇关于跨音松驰法的奠基性文章 ,文章提出了求解小扰动势方程的混合差分格式。文献 [2 ]发展了求解全位势方程的旋转混合差分格式。近年来 ,求解跨音全位势方程的方法有了更全面的发展 [4、5]。作者于 80年发展了一种求解大弯度二维叶栅的跨音松驰法[3] 。本文将该法推广到求解涡轮平面叶栅跨音带激波流场 ,直接求解全位势方程 ,进一步提高边界区差分格式的精度 ,并注意初场的给法。本法具有节约机时和内存的优点。数值算例表明 ,对于头部不大的涡轮平面叶栅 ,应用本方法…     

直升机旋翼悬停流场的欧拉方程计算

描述了三维Ｅｕｌｅｒ方程数值模拟旋翼悬停流场的过程,欧拉方程的求解采用了风格中心有限体积法、五步Ｒｕｎｇｅ－Ｋｕｔｔａ时间推进格式。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.     

涡轮叶栅叶型损失的数值模拟方法

给出了一个计算亚、跨音涡轮叶栅叶型损失的数值计算方法。主流采用时间推进有限体积法求解积分型欧拉方程,并采用了局部网格修正技术;附面层采用全隐格式求解有限差分方程;在叶栅出口与远后方均匀流之间进行了叶片尾迹与主流的掺混损失计算。算例表明本文的数值方法可准确地预测涡轮叶栅的叶型损失。     

导弹级间热分离流场的数值预示

针对多级导弹级间热分离这一复杂环境条件,用数值模拟的方法对级间段流场进行了预示。计算采用Ｇｏｄｕｎｏｖ有限体积格式,对非定常轴对称欧拉方程进行求解。计算结果为级间分离段和发动机热防护的设计提供了理论依据。     

激波运动与反射流场的ENO算法

在ＴＶＤ格式的基础上，以修正限制函数的方法，将通常的ＴＶＤ格式改造为性能更为优越的ＥＮＯ格式，并利用这种有限差分格式求解欧拉方程和薄层Ｎ－Ｓ方程。分别数值模拟了激波管中运动激波与三种楔形物体相互作用流场；给出了激波在不同物体上的反射及波涡干扰过程中非定常流场的变化情况；算出的波系结构在与ＣＡＲＤＣ小型激波管中做的相关实验结果一致。     

跨音速压气机转子内三维流场数值模拟

采用多步Runge-Kutta与有限差分相结合的格式,离散叶轮机械内三维欧拉方程,应用代数法生成了压气机转子三维网格。用此方法模拟了NASA Rotor37跨音速压气机转子流场,计算得到的流场与试验结果进行了对比。结果表明,该方法收敛速度快,实用可行,计算简单省时。并且,该方法为研究叶轮机械内部流场计算提供了很好的基础。在此基础上进行更深层次的研究会变得很容易。     

二元喷管三维粘性流场的数值模拟

本文用数值方法对二元喷管的三维粘性流场进行模拟,将 Mac Cormack有限差分格式拓广到全三维 N-S方程的求解。在数值迭代中使用加快收敛速度的当地时间步长,并根据格式压缩计算机存储量和计算量。进行了两个算例的数值模拟,计算结果与实验相符。      

跨音速喷管化学反应流数值模拟

用时间相关的半隐格式有限差分数值方法求解了化学非平衡反应跨音速喷管流场，在喷管收敛段，流动接近化学平衡状态，控制方程的刚性问题严重，数值积分困难。通过对时间差分项隐式离散、对空间差分项显式离散，流场边界采用参考平面上的特征线法计算等，成功地解决了由于化学反应有限速率带来的数值解不稳定问题。该格式简单、需要计算机存贮空间少。本文完成了一维和轴对称非平衡化学反应喷管流动计算，并与化学平衡流和冻结流的计算结果做了比较。     

GMRES算法在悬停旋翼数值模拟中的应用

发展了一种基于广义极小残差(GMRES)算法的悬停旋翼数值模拟方法,并对GMRES算法中矩阵与向量乘积的两种计算方法进行了分析和讨论。应用该方法在旋转坐标系中采用非结构混合网格和格点格式有限体积法对以绝对速度为变量的欧拉方程进行了直接求解,其中对流项的离散应用了基于Roe的Riemann近似解的迎风格式。对Caradonna-Tung旋翼跨声速悬停流场进行了数值模拟,计算结果与相关实验数据吻合较好,并与LU-SGS方法进行了对比,表明GMRES算法可以有效地加速流场的收敛,提高计算效率。     

超声速主流中逆向喷流流场的数值模拟

本文在差分格式NND-2基础上^[1]，提出了一个修正格式，形式简单，应用方便。从薄层近似的NS方程出发，用该格式计算了超声速主流中存在逆向喷流的钝体绕流的粘性流场，成功地捕捉到流场内的各种波系和涡系结构。计算的结果与文献[2]中的实验结果比较表明二者符合很好。     

一种求解透平叶栅三维流场的高精度TVD格式

本文提出了一种应用三阶精度Chakravarthy-OsherTVD格式求解环形透平叶栅三维定常流场的隐式近似因子分解方法。控制方程为柱座标速度分量表示的任意曲线座标系下的三维欧拉方程。采用Beam-Warming近似分解, 对隐式项进行迎风差分, 得到了时间精度为一阶、空间精度为5点三阶的隐式格式。算例对5点三阶和5点二阶精度的格式进行了比较。表明空间三阶精度隐格式对计算结果没有明显改进, 但收敛残差有较大降低。      

超燃冲压发动机三维喷流干扰流场的数值研究

用隐式有限体积法求解三维N－S方程数值研究超燃冲压发动机三维平行、垂直喷流干扰流场。计算中采用了具有当地极值递减（LED）性质的JST和SLIP中心差分格式。为了减缓计算网格数很多时对计算机内存的压力以及提高计算效率，运用了通量守恒的分区计算方法。     

ISOMS——计算粘性流场的一种隐式差分格式

用稳定化方法求解可压缩完全Navier-Stokes方程(CCNS方程)是计算复杂流场的一种有效的手段。我们在文献[1]、[2]中利用新的开关函数,构造了一种新的混合格式,称为SOMS格式,用于数值求解CCNS方程。本文是将SOMS格式推广到隐式混合格式,称为ISOMS格式(The implicit second order monotone scheme)。它具有如下几个特点:ⅰ)它是隐式格式,稳定性良好;ⅱ)穿过激波是单调的;ⅲ)对于     

微型飞行器低雷诺数三维流场数值模拟

针对当前所研究的微型飞行器(MAV)流动特征,采用时间相关法和有限体积法求解N-S方程,选择了Van Albada插值限量因子的高分辨率空间差分格式,建立了模拟MAV三维低雷诺数流场的数值计算方法。通过熵校正技术改进和粘性计算准确性与网格密度关系的分析,确保了流场计算的准确可靠。使用该数值算法,对某常规气动布局的MAV进行流场模拟和分析,并比较了其流场与常规尺寸飞行器的不同。     

高超声速化学非平衡欧拉方程数值模拟

本文给出了求解高超声速化学非平衡无粘流场的数值模拟方法。计算中采用守恒变量型的ＮＮＤ差分格式和全流场捕捉技术;建立了（数学）奇性轴上的控制方程,克服了通常在其上插值方法所引起的扰动,改善了收敛性和计算精度;化学反应模型由五组元五反应组成;温度和压力的计算是非迭代的,有别于国内外惯用的迭代算法;此外,还发展了背风过低压区的处理方法,从而提高了欧拉方程的应用范围,使无粘非平衡绕流的计算攻角达到３０°     

多重网格法求解二元跨音速欧拉流

 近年来,采用时间相关欧拉方程或N-S方程求解跨音速流得到了迅速发展。为了获得二元或三元复杂形体正确的表面气动参数和足够高分辨率的流场物理特性,必须采用十分细密的网格,因而导致计算时间成倍增长。为此,除继续发展各种差分格式,探讨各种自适应网格技术外,运用多重网格技术加速收敛过程得到了国外普遍关注并取得了一定的成功。国内也已有采用多重网格法解跨音速欧拉流的文章,然而均局限于两重网格的计算。为探索多重网格法和焓阻尼技术对加速收敛的作用,我们以NACA0012翼型的跨音速流动为例,采用MacCormack两步显式格式进行了数值实验。     

二维Euler方程的自适应网格矢通量算法

给出了二维轴对称欧拉方程的自适应网格隐式矢通量分裂差分格式和模拟结果,隐式矢通量分裂格式不需人工粘性项与人工耗散项,可以计算跨声速流场,也可以较准确地捕获激波,利用流场解的梯度信息对网格进行调整,可以提高激波的捕获精度,算例表明,该方法精度高稳定性好,不依赖于初始条件。     

烧蚀外形三维粘性流场数值模拟

本文通过求解层流Navier-Stokes方程,数值模拟了三维奶头型烧蚀外形粘性流场。差分格式为Yee和Harten的迎风隐式TVD格式,在空间方向具有二阶精度。计算得到了高质量的激波和合理的流场结构。物面压力分布与实验结果进行比较,符合较好。     

叶轮机转子内流动的三维跨声Navier-Stokes方程组数值计算

应用有限体积法及Baldwin-Lomax代数紊流模型求解叶轮机转子内流动的三维可压缩非定常薄层Navier-Stokes方程组,计算一个绕X轴等转速旋转的转子内流场问题,用当地时间步长方法来加速收敛。采用隐式算法对对流通量进行求解,而扩散通量采用显式算法。采用高阶TVD下的Roe格式计算对流通量,采用中心差分格式计算扩散通量。关于守恒变量Q的非线性方程组,采用牛顿迭代方法求解,线化后的方程用对称的Gauss-Seidle方法求其渐近解。应用该算法模拟了转子NASARotor67流场,与试验做了大量的对比研究,包括性能、等马赫线以及总压和总温等等。在最高效率点附近计算与试验比较接近,但在靠近失速点时计算与试验略有差别。