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

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

基于N-S方程和离散共轭方法的气动设计方法研究

对于基于梯度信息的优化设计方法,很重要的一点是快速准确获得目标函数对设计变量的梯度.本文采用离散共轭方法计算目标函数关于设计变量的梯度,流动控制方程为N-S方程.对于离散共轭方程和流动控制方程均采用LU-SGS方法求解.算例表明,对于亚声速和跨声速两种情况,该方法都能快速准确地获得升力和阻力关于设计变量的梯度.本文采用该方法进行了翼型优化设计,成功地减弱了激波,降低了总阻力.算例证明了本文方法可靠性好、收敛快, 特别适合工程实际.     

基于特征值理论的轴流跨声速压气机失稳预测

流动失稳问题是叶轮机械所面临的严峻挑战之一,而旋转失速是轴流压气机及压缩系统中普遍存在的一种典型流动不稳定状态。基于特征值理论发展了一个轴流压气机流动失稳开始点预测模型,应用谱方法和奇异值分解(SVD)算法求解所建立的特征值方程,通过特征频率的虚部判断压缩系统的流动稳定性。并通过一个高速单转子在高低两个转速时的算例,定量分析了流动压缩性对跨声速压气机流动失稳开始点预测的重要影响。该模型可在合理的精度范围内,以工程应用可接受的计算资源消耗来预测跨声速单级压气机旋转失速的开始点流量,同时不应用大量经验数据,为压气机设计阶段的稳定性评估提供指导。     

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

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

跨声速非等熵势函数方程的数值求解

基于动量方程的特点,本文将只与温度有关的“等熵密度”与熵的变化分离,而由散度型动量方程直接计算出气体经过激波的熵增。与经典势函数方程迭代,可以方便地求解非等熵跨声速流动。从对任意回转面叶栅的跨声速流动计算来看,这种非等熵势函数得到的激波位置可比等熵势函数前移1～2个网格,强度有所减弱,和试验结果较接近。由于各流线的熵不相等,对Kutta条件和出口边界条件作了修正。本文建议的方法也可用于计算流函数的熵增。     

轴流压气机孤立转子内流场数值模拟

针对跨声速轴流压气机转子叶排中的流动，利用高分辨率的三阶NND格式和LU-SGS隐式推进迭代法，建立了孤立转子内三维N-S方程的高分辨率和高效率的数值分析程序。用该程序分别计算了NASA-37跨声速转子，计算结果表明：所建立的孤立转子程序能够较好的模拟轴流压气机叶排内流动情况。数值分析结果与实验吻合较好，从而证明了该程序的有效性。     

分区并行跨声速流的计算

本文采用求解Ｎ－Ｓ方程方法计算了轴对称钝头倒锥体外形的跨声速大迎角层流流动。     

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

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

跨声速叶栅粘流计算的多级Taylor—Galerkin有限元法

本文对Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法在收敛速度和稳定性两方面作了改进，与多级有限元法相结合，构成了多级Ｔａｙｌｏｒ－Ｇａｌｅｒｋｉｎ有限元法。基于Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，研究了平面跨声速叶栅流动，并与实验结果作了比较。计算结果表明新方法具有较好的稳定性及较快的收敛速度。     

含摩擦和加热的跨声速流动数值模拟研究

喷管中跨声速流动的反应流体动力学是实验和理论研究的一个重要领域.本文主要由两部分构成,在第一部分中对文献[2]中的Murman格式[1]的一般形式进行了简化,得到了守恒形式的欧拉方程,并对其求解方法进行了讨论,同时针对Chattot等的工作展开了分析,并给出了带激波的流动计算结果;而文章的第二部分则研究了包含摩擦、加热以及耗散的跨声速流动,并重点关注了一步Arrhenius化学反应模型下的爆震现象.本文还探究了摩擦力和带燃烧的跨声速流动之间的相互作用,同时详细给出了新的欧拉方程数值求解算法并同已有文献算例进行了对比验证.     

超燃冲压发动机流量匹配的临界流量法

为了使用基于积分方程的超燃冲压发动机性能计算模型模拟跨燃、亚燃工作模态,用临界流量法建立用于表征燃烧室热壅塞程度的残差方程,并以隔离段出口静压为独立变量,通过牛顿迭代法迭代求解残差方程,模拟了燃烧室内存在临界截面（即喉道）但不出现热壅塞的实际跨燃、亚燃工作状态。结果表明：该方法能够确定有喉道的超燃冲压发动机燃烧室一维参数分布,并具有计算精度高、计算速度快、收敛性好等优点,可以为基于积分方程的超燃冲压发动机性能计算模型求解提供一定的参考。     

大尖削机翼非定常跨音速绕流的数值解法(英文)

本文给出一种计算飞机机翼上定常和非定常跨音速气动力的数值解法,使用了一种特殊设计的坐标变换,是用时间精确交替方向隐式(ADI)有限差分算法来求解非定常跨音速修正三元小扰动位势方程。给出了F5战斗机机翼的数值结果并与XTRAN3S,ATRAN3S及试验结果进行了比较,表明本方法是有效的和经济的。     

超燃冲压发动机二维热环境数值模拟

对超燃冲压发动机热环境进行了研究和计算.以二维N-S方程和一维瑞利加热规律为基础, 建立了超燃冲压发动机内流场的热环境计算模型;用MacCormark预测-校正格式编写了数值计算程序;对飞行马赫数Ma=6和Ma=8两种工况的超燃冲压发动机内部热环境进行了数值模拟, 得到了发动机内部流场的温度分布和壁面热流分布;对计算结果进行了分析.结果表明, 采用二维无反应流体计算和一维加热规律相结合的方法, 在超燃冲压发动机热环境研究中是可行且有效的.      

释热分布对超燃冲压发动机性能的影响及优化

为了改进燃烧室、尾喷管的设计和优化燃料在燃烧室内的释热方案,建立了广义一维流动/释热分析方法,对超燃冲压模型发动机的主要性能参数进行了研究。通过与单凹腔模型发动机实验结果的对比,着重分析了燃料在燃烧室内不同的释热分布对模型发动机推力性能的影响,同时采用模拟退火优化算法对发动机燃烧室内释热分布和内流道作了优化,结果表明,在燃烧室内采用多阶段分散释热方案对模型发动机推力性能有较大的提高。     

预处理法求解定常/非定常混合网格的全速流场

运用时间导数预处理法在混合网格上求解三维定常/非定常N-S方程。预处理后的N-S方程,用二阶迎风格式的有限体积法离散,隐式、LU-SGS迭代求解,双时间步推进求解非定常流。本文以一维预处理欧拉方程为基础进行分析,提出了一种特征变量的边界条件以适应预处理后的特征系统。从极低速到跨声速范围的众多算例表明,本文的边界条件处理稳定有效,预处理法在各种速度下都具备同样快速的收敛性,在较宽的使用范围内均能得到理想的计算结果。     

跨音速三维机翼及翼身组合体非守恒型全位势流AF3算法

 本文对跨音速非守恒型全位势流,用Baker的快速有限差分法,改变其网络拓扑,研制出分析机翼及翼身组合体跨音速绕流的实用程序。其中包括:对“C-H”型贴体坐标网格提出AF3算法的具体实施方案取得了快速收敛效果;采用张量形任意曲线坐标上的非守恒全位势方程,使得程序中跨音速求解部分既能用于分析单独机翼又能用于翼身组合体跨音速绕流计算。     

旋成体轴对称跨声速全位势流的高效差分算法

本文用守恒型全位势方程,贴体坐标网格,对旋成体轴对称跨声速绕流的差分数值计算方法进行了研究;根据最佳收敛准则,提出了轴对称情形的AF2迭代算法,并将此算法应用于半球头柱体、弹体等各种外形的旋成体。与一般方法仅适用于亚声速自由流不同,本计算可从亚声速、跨声速自由流一直到低超声速自由流。计算结果表明,本文方法收敛快,与实验及其它方法的结果符合较好。     

二维非定常跨音速时间线化积分方程的推导

1.引言 求解非定常跨音速流动的主要困难是非线性问题。对于微幅翼型振动问题(同时引起激波的微幅振动)可做时间线化简化处理,得到时间线化微分方程定解问题。时间线化积分方程是由时间线化微分方程推导出来的。     

用隐式近似因子分解法计算喷管跨音速流场

本文用隐式近似因子分解法结合贴体曲线坐标计算了定常、无粘、跨音速喷管流场,选择了参数变化比较剧烈的小喉部曲率半径喷管,陡壁收敛段喷管和潜入喷管三种算例.计算结果表明,该方法收敛快,并且具有良好的精度.Courant数可以取至7.这种隐式法与显式法相比较,花费的机时少得多.本文的方法可望推广到求解N-S方程、PNS方程,和二相有粘喷管流动计算.