轴对称跨音速进气道外罩绕流的高效差分算法

 本文考虑轴对称跨音速进气道外流场的有限差分计算方法。采用守恒型位流方程、贴体坐标网格和精确边界条件,根据最佳收敛准则对轴对称情形设计出新的近似因式分解迭代格式,并将所设计的格式应用到计算皮托式进气口的轴对称跨音速流场。对几种典型进气口的计算表明,本文格式收敛快,计算结果与实验符合很好。     

关于跨音速叶栅物面边界条件的数值试验

本文讨论了跨音速叶栅物面边界的数值处理方法,目的是寻求与我们所用内点格式相匹配的边界上流动参数的算法,以提高计算精度。本文首先推导了适合于物面边界计算的特征关系式,其次总结了五种边界处理方法,进行了数值试验。计算结果表明,配合我们的内点计算格式,以y向抛物外推的Denton修正法计算效果较好,能够较精确地反映跨音速叶栅带激波流场的物理特点。     

高分辨率TVD格式的应用

在Jameson格式的基础上，利用修正通量法，构筑TVD（Total Variation Diminishing)格式。通过三个典型的一维和二维算例的计算结果表明，TVD格式对激波具有较高的分辨率，它为具有复杂激波系的喷管、叶栅等跨音速流场计算提供了一种有力的工具。     

绕翼型跨音速大扰动对称势流的插值混合差分法

 本文将跨音速定常小扰动势流混合差分法推广到跨音速大扰动定常势流,提出了在局部速度坐标系中求解跨音速精确势流方程的插值混合差分法。作为算例,计算了双圆弧翼型和NACA0015翼型对称问题压力分布,并与已知实验值和双圆弧翼型小扰动混合差分法计算值进行比较,结果接近。试算表明,本文提出的插值混合差分格式是稳定和收敛的。本文解决了M_∞趋近于1的计算难点。     

多级风扇三维粘性流场的数值模拟

 进行了多级跨音速风扇粘性流动的数值模拟研究。求解全三维N-S方程,通过高分辨率格式提高跨音速流场计算的准确性,利用LU-SGS隐式解法加快计算收敛速度。给出了一双级高负荷跨音速风扇设计和非设计工况下的详细计算结果,各工况下的计算与实验数据均吻合较好。     

LTRAN2的一种改进方案

本文提出了LTRAN2的一种改进方案BTRAN2用来分析翼型的各种频率的非定常运动。用E-O调转换格式和ADI方法求解了完全的二维非定常跨音速小扰动位流方程,并用单调转换的AF2格式计算了定常跨音速小扰动方程,以此作为非定常计算的初场。本文给出了带后缘正弦振荡襟翼的NACA 64 A006翼型的绕流和做正弦俯仰振荡的NACA 64 A010翼型的绕流计算结果,它们与Euler方程解或实验数据很吻合。     

先进上风格式跨声速CFD计算研究

介绍了空气动力学跨音速计算中应用广泛的几种上风格式,同时引入了JAMSON中心差分格式作为比较,针对激波分辨率,格式精度,计算效率等方面的问题对ONERAM6机翼进行了对比性计算,讨论了ROE的FDS格式,AUSM+格式和AUSMDV格式的计算性能,得出了很有意义的结论。     

矢通量分裂-TVD杂交新格式及其应用

在贴体曲线坐标系中,本文提出了一个高分辨率的新格式,它是Steger-Warming的矢通量分裂与Harten的TVD(Total Variation Diminishing)格式的杂交,用于计算跨音速定常流动和捕获激波。典型的跨音速内流和叶栅流的算例表明:在60×15网格下捕获的激波过渡区不超过两个网格,证明其分辨率较高;分析激波附近的数值结果,没出现“低亏、过跳”、伪振荡现象;与实验结果的比较表明,对定常流在无人为附加耗散项的条件下它能给出较高质量无数值波动的激波流场解。     

跳点有限体积法及其在跨音速流场计算中的应用

[1]和[2]曾用几种差分格式求解一维N-S方程和欧拉方程,计算表明,Hopscotch格式(跳点格式)计算速度最快。为探讨这种格式在多维流场特别是跨音速流场计算中应用的效果,本文将此差分格式应用于有限体上,成为跳点有限体法,并用该方法求解二维欧拉方程组。对若干个喷管和叶栅所做的计算得到了满意的结果,计算时间较目前常用的差分格式少得多。     

单级跨音速压气机三维非定常粘性流动计算

本文以数值求解雷诺平均非定常N-S方程方法完成了单级跨音速压气机动静叶相干三维非定常粘性流动计算,给出了高分辨率NND格式的时间推进有限差分解。对某跨音速单级压气机进行了定常和非定常流的详细数值计算,设计转速下定常粘性流场计算所得到的性能参数与实验结果吻合较好。在此基础上,对最高效率状态进行了动/静叶相互干扰的非定常流动计算,并对结果进行了分析。      

LU-AUSMPW混合格式在可压缩无粘和粘性流动数值分析中的应用

本文首先将 AUSMPW格式与三阶 MUSCL格式融合, 给出了其在任意曲线坐标下的三维形式, 并将其与 LU-SGS格式结合, 应用于可压缩 Euler和 Navier-Stokes方程的求解。其次, 构造了一种新的通量限制器。最后, 为了验证 LU -AUSMPW混合格式的性能, 将平面叶栅跨音速无粘流动以及喷管超音速粘性流动作为算例。本文计算结果与文献计算结果和实验数据相符很好, 表明采用 LU-AUSMPW混合格式数值模拟可压缩流场, 具有较高的计算精度、较快的收敛速度和良好的稳定性。      

跨音速小扰动方程的多重解

本文利用Engquist-Osher格式,离散跨音速小扰动方程(TSD)和边界条件,选择不同初场计算NACA0012翼型二维定常位势流场,得到了与试验值较符合的结果,同时用Murman-Cole非守恒和守恒两格式计算并进行了比较。着重用数值实验进一步探讨多重解问题,发现并总结出了几点很有参考价值的规律和现象。     

直接的矩阵分裂计算跨音速喷管流场

给出了分离系数矩阵求解欧拉方程组时系数矩阵的直接分离公式，减少和简化了应用该方法时的矩阵运算，参照Ｂｅａｍ－Ｗａｒｍｉｎｇ显式格式构造与ＳＣＭ方法一致的内点差分格式，并分析了它的数值特性，对有斜激波的跨音速内流场进行了数值仿真。计算中采用变系数的当地时间步长，加快了流场达到的收敛速度，计算无需引入任何人工参数，计算值与实验结果吻合。     

跨音速压气机级三维紊流流场数值模拟

在对跨音速压气机级的动／静双叶排中的三维定常紊流流场数值的模拟中，利用三阶高分辨率ＮＮＤ格式和ＬＵ－ＳＧＳ隐式推进迭代法，既保证了流场中激波的模拟质量和粘性流动特征的正确预估，也实现了求解过程的高效率。对某单级跨音速压气机的计算及其与实验数据的对比证实了本文方法的有效性。将计算推广至多级轴流叶轮机械的情况将会是十分简便的     

进气道物面贴体坐标生成的几何自适应控制方法

针对轴对称进气道跨音速绕流问题的特点，用几何自适应的方法和椭圆型方程构筑了适合流场特征的贴体坐标系，并用AF3格式计算了A－27进气道物面上的压力分布。计算结果表明：几何自适应方法能够控制进气道物面上网点的疏密，给出进气进物面上网点的良好分布，从而为进气道跨音速绕流问题的准确计算奠定了基础。     

几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

实际机身跨音速流场分析计算

本文应用修正的跨音速小扰动(MSD)方程和疏密网格交替迭代法,将Jameson格式在非等步长网格中的推广形式与Engquist-Oslier格式相结合,在机身侧边局部横向网格线上使用横向线松弛法,经数值实验得到了一个收敛速度比较快的计算方案。许多算例均证明本方案迭代次数少,计算结果与实验结果符合得也很好。     

皮托管式轴对称进气道跨音速内,外流场的AF—2方法计算

本文构造了轴向大扰动,径向小扰动轴对称速势方程的AF-2分解格式,计算了IC型皮托管式进气道的跨音速内外流场,并对这种分解格式进行了线化的稳定性分析。分析表明:AF-2方法的稳定性优于SLOR方法,其收敛速度至少比SLOR方法快10倍。     

跨音速定常小扰动压力方程的混合差分法

本文提出直接求解跨音速定常小扰动压力方程的数值方法。对于研究某些洞壁干扰问题,与传统的速势方程相比,用压力方程作为求解跨音速流场的主管方程,边界条件为Dirichlet形式,易于处理,且待求变量为压力,可减少积累误差,提高计算精度。 本文采用混合差分法求解压力方程,通过数值试验,确定合适的差分格式及迭代线化方法。其收敛解与相应的以速势方程为主管方程求到的解相比吻合得比较好,从而证实了本文方法的可行性。 最后给出应用本文方法计算鉴定跨音速翼型风洞壁干扰以及由给定的压力分布计算翼型外形的典型算例。     

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

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