三角翼大攻角绕流场数值模拟研究

本文基于N-S方程采用Jameson的显式中心有限体积格式和Roc的上风格式对三角翼大攻角绕流场进行初步数值模拟，空间网格为H-0型网格。计算结果表明两种计算格式均可较好的模拟大攻角绕流场。本文计算结果为今后开展全机大攻角绕流场数值模拟提供了计算经验。     

悬停旋翼无粘流场数值模拟中的多重网格方法

发展了一种基于不同空间离散格式的多重网格算法,并应用于悬停旋翼无粘绕流的Euler方程数值模拟。由于悬停旋翼流场中存在不可压区域,同时旋翼尾涡系统的发展需要较长的时间,使得旋翼流场的收敛速度远低于固定翼流场,因此研究旋翼流场的多重网格算法具有重要意义。空间离散采用了Roe s FDS格式和Jameson中心有限体积格式,时间推进应用了五步Runge-Kutta方法。采用多重网格V循环方式,对一跨声速悬停旋翼无粘流场进行了数值计算。计算结果表明:多重网格算法可以显著加速悬停旋翼无粘流场的数值计算收敛速度;无论在激波分辨率还是在计算精度上,Roe s FDS格式都优于JST格式。     

非结构网格上的TV-van Leer混合迎风格式

为了改进迎风型对流通量格式的激波捕捉计算稳定性,对迎风格式产生数值激波不稳定现象的原因进行分析,并且在分析结果基础上构造新型迎风型格式。为此,首先设计了数值算例,分析非结构网格上迎风格式出现数值激波不稳定现象的主要影响因素。结果表明,在一定条件下,网格和流动状态共同作用导致产生不合理的数值结果,并影响计算的稳定。其中,网格间的动量通量对数值激波不稳定现象影响显著。因此,基于TV格式构造了一种自适应调节动量通量耗散项的混合格式,从而在激波附近引入耗散,抑制扰动的产生和传播。数值计算结果显示新型混合格式在非结构网格离散条件下能够稳定地捕捉激波,同时混合格式保持了TV格式原有的黏性流动计算分辨率。     

激波／边界层干扰数值模拟的格式效应

用不同的数值格式,包括Ｒｏｅ的ＦＤＳ格式和ｖａｎＬｅｅｒ的ＦＶＳ格式并配备不同限制器,求解二维Ｎ－Ｓ方程,模拟研究了激波／边界层干扰的流动问题,分析了模拟出的流动现象,仔细研究了其中的格式效应,得出同前人不尽一致的结论：ｖａｎＬｅｅｒ格式对限制器的依赖性比Ｒｏｅ格式更强;在本文情况下,ｍｉｎｍｏｄｌｉｍｉｔｅｒ的激波、膨胀波的分辨率比ａｌｂａｄａｌｉｍｉｔｅｒ还好     

二维埋入式进气道的数值分析

采用并改进了适用于分离流动的 Ｊｏｈｎｓｏｎ Ｋｉｎｇ 紊流模型,结合 Ｌ ＵＳ Ｇ Ｓ近似隐式格式和基于 Ｍ Ｕ Ｓ Ｃ Ｌ非线性插值的 Ｒｏｅ 格式编写了雷诺平均的二维 Ｎ-Ｓ方程的求解程序。然后借助于该流场计算程序完成了二维埋入式进气道的数值研究,分析了几何设计参数对进气道性能的影响。     

用隐式WENO格式计算悬停旋翼跨声速流场

发展了一种基于高阶迎风格式的计算悬停旋翼跨声速流场的隐式有限体积法。对流项采用Roe通量差分分裂格式,使用五阶WENO格式进行左右状态重构,并与MUSCL插值进行比较;粘性项采用中心有限体积法。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。数值模拟采用了一种能够有效传递网格间流场信息的重叠网格,其中使用了三层内边界和贡献边界的方法以便插值的直接进行。用该方法对一跨声速悬停旋翼粘性流场进行了数值计算,数值结果表明:与MUSCL格式相比,WENO格式对激波位置捕捉得更准确,具有更强的涡捕捉能力,同时还表明了WENO格式在很大程度上能够克服涡耗散问题。     

不同矢量偏角下二元喷管内流特性的数值研究

应用多重网格MacCormack显式格式对不同矢量角下二元喷管的内部流动进行了模拟分析。计算时, 在细网格上采用了雷诺平均N-S方程与k-ε湍流模型方程;在粗网格上则采用Euler方程。数值计算给出了较清晰的流谱, 壁面压力分布的计算也与实验结果相吻合。通过计算结果分析了矢量喷管几何参数对流场的影响, 并对矢量喷管的推力特性进行了数值研究。      

高超声速数值激波失稳的网格依赖性

高超声速强激波的稳定捕捉仍极具挑战性。目前工程计算中普遍应用的数值格式在模拟多维强激波时通常会遭遇明显的激波不稳定现象,且数值格式的激波稳定性对计算网格表现出严重的依赖性。基于矩阵稳定性分析法,对比了具有不同耗散性质的数值格式稳定捕捉激波的能力,分析了空间二阶精度格式的激波稳定性及限制器对激波稳定性的影响。在此基础上,重点探究了计算网格对激波稳定性的影响规律,研究了网格长宽比和畸变角度对激波稳定性的影响。结果显示,在激波附近通过增大网格长宽比或改变网格畸变角度可有效改善激波捕捉的稳定性;相比于增大网格长宽比,改变网格畸变角度提升激波捕捉稳定性的效果更加明显。在此基础上,结合数值耗散分析对高马赫数下数值激波失稳现象的网格依赖性机制进行了探讨。     

热流CFD计算中格式和网格效应若干问题研究

通过数值实验、理论分析、试验数据和经验公式的综合研究,探讨了CFD热流计算的格式效应及网格效应:研究了Roe的FDS格式、Van Leer的FVS格式、AUSM 格式及中心差分格式的热流计算分辨率,探讨了不同格式对热流计算精度的影响,认为AUSM 格式在热流计算精确性方面具有优势;在热流绝对值的模拟上,网格的影响很大,能否计算出正确的驻点热流值将是获得精确物面热流解的关键;通过大量的数值计算、理论分析及其同试验数据的对比研究,对可能影响热流计算的网格因素进行了细致的分析,揭示了各因素对热流计算结果的影响规律。     

通量限制器与LED性质、空间离散精度

本文从标量形式的ＵＳＬＩＰ格式出发,从理论上分析了通量限制器与ＬＥＤ性质和空间离散精度的关系。利用Ｒｏｅ平均矩阵的Ｕ特性和广义数值Ｒｉｅｍａｎｎ不变量的概念,将此类格式应用于二维翼型跨声流场的数值模拟,针对提供的来流条件,考查了不同的通量限制器对流场计算精度的影响。     

透平跨音速叶栅的优化设计

将遗传算法用于透平跨音速叶栅的优化设计,提出了相应的数学模型和优化方法。优化算例表明,本文所提出的模型及算法能有效的优化叶型方案,避免了一般优化算法的局部最优或维数灾难问题,从而在新的领域内将优化计算在叶轮机械设计中的应用进行了有益的探索。      

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

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

低展弦比高负荷跨音速冷却涡轮特性的计算方法

本文在AMDC叶栅损失公式基础上引进了KQ系统,根据高负荷跨音速的特点,对AMDC超临界波阻公式作了修正,选用了合适的栅后气流角公式。用一元可压缩轴对称理论推导了冷气同主气流掺混公式。最后给出高亚音速和跨音速模型涡轮特性以及低展弦比高亚音速冷却涡轮设计点性能和试验结果的比较,精度为1％。     

一种削弱跨音叶栅流场中激波强度的数值方法

本文提出了削弱跨音叶栅流场中激波强度的两个计算方案,并给出了两套叶栅的计算结果。数值结果表明,削弱激波强度的效果是明显的。文中提出的数值方法可方便地纳入现有压气机气动设计系统付诸实用。本文的工作是初步的,并有待实验验证。     

跨声速涡轮叶栅三维流场的数值模拟

以具有试验测量数据的一个高压涡轮动叶中截面平面叶栅和一个高压涡轮导向器环形叶栅为对象,采用计算流体动力学(CFD)软件CFX-TASCflow,对两个涡轮叶栅的跨声速三维粘性流场进行了数值计算,验证了CFX-TASCflow软件在航空燃气涡轮设计计算中的可靠性,并对所研究的跨声速涡轮叶栅粘性流场进行了分析.      

跨音速平面叶栅全位势无粘-粘性迭代计算

本文发展了一种全位势与附面层、尾迹流正逆模式积分法迭代计算程序,给出了计算结果与国内外四个不同的叶栅风洞试验结果的对比分析。结果表明,对比是较为满意的,本文在网格生成,无粘跨音流求解和无粘流与附面层、尾迹流迭代等方面所作的发展和改进适用,使整个程序系统具有工程实用价值。      

跨音涡轮中三维流动的数值试验

 本文从积分型的基本方程组出发,以我们的关于平面叶栅和任意旋成面叶栅的工作为基础,将文献（3、4）中的计算方案推广到三维流动情况。编制了相应的计算机程序,进行了一些数值试验,得到了初步的计算结果。 数值试验表明,拟三维计算结果可以导致相当大的误差,从而说明了完全三维流动计算的必要性。另一方面,如果取2000个左右的计算点,则在国产的013机上算一个例子,约需100分钟。从实际应用的观点来看也是可行的。     

点隐式龙格-库塔方法的应用研究

为了提高求解器的效率,在显式龙格-库塔时间推进的欧拉方程求解器之上,发展了点隐式龙格-库塔时间推进格式。给出了其推导过程和非结构网格下中心格式和迎风格式（Roe和Van Leer格式）预处理算子的构造方法。NACA0012翼型和RAE2822翼型的跨音速无粘流动模拟表明：与显式龙格-库塔方法相比,方法能提高求解效率且内存需求相当,具有一定的工程应用价值。     

两个涡轮叶栅跨声速流场的数值模拟

采用CFD软件CFXTASCflow时两个涡轮叶栅跨声速流场进行了分析，其中一个为高压涡轮导向器环形叶栅，另一个为高压涡轮动叶中截面平面叶栅，计算结果与试验数据吻合良好。     

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

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