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

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

数值求解二维定常跨音速流场多重解现象探讨

本文用时间推进法求解欧拉方程组计算二维平面叶栅绕流问题数值实验是否有多重解现象。计算结果表明,求解欧拉方程组进行流场数值模拟最终收敛与初场无关,不存在解不唯一问题,进而计算初场熵值大小,也发现初场熵值大小不影响终场结果。     

大扰动跨音速流积分方程的导出及亚临界流的数值解法

本文对一方向上大扰动的平面跨音速势流微分方程提出了新的近似简化方案,从数学上完成了相应积分方程的推导;其中关于边界及激波条件作了小扰动近似处理;最后提供了关于亚临界流的算例。     

关于跨声速流动问题的多重解

本文讨论了跨声速流小扰动方程Riemann问题的多重解及物理上可实现的解的判别。     

新的解跨音速Euler方程的隐式杂交方法

 本文提出一个新的杂交格式,它由Jameson-Turkel LU分解法和Harten TVD（total variation diminishing）格式的修正数值通量耦合而成。LU分解用于Euler方程隐式部分的离散化而Harten的数值通量控制所谓的稳态解残差。精心构造LU算子,可使矩阵算子L和U很容易逐点推进求逆;以Harten的数值通量计算杂交格式右边的通量,这样做在a=1/2的条件下,可保证收敛的稳态解具有二阶精度并能高分辨率地捕获激波。数值实验表明,本文格式对于稳态计算是有生命力的并能保证较高的激波分辨率。     

跨音速非定常流动的Euler方程隐式解

在随时间变化的贴本坐标系中给出求解非定常Euler方程的连续通量分裂法。在此基础上建立了可用于跨音速非定常流动的Euler方程隐式求解法。采用特征向量变换,可在保证原方程组离散化精度的条件下使计算大为简化。针对振动翼绕流特点建立了固连于物体的动坐标与固定坐标间的关系。数值计算在动坐标中进行,既简化网格生成又保证在物面上满足边晃条件。对NACA64A-10冀型绕1/4弦点做简谐俯仰振动的非定常气动力进行了计算,给出了与实验结果基本相符的计算结果。此外,还给出翼型做沉浮及同时进行沉浮与俯仰二自由度振动的非定常气动力的计算结果。     

叶栅跨音速流场的Euler解

应用Jameson的有限体积法求解二维Euler方程，模拟叶栅跨音速流场，并运用当地时间步长和隐式残差平均技术加快收敛速度，计算结果与理论解和试验数据吻合良好，本文比较了采用两种不同的通量计算公式的计算结果。     

跨声速翼型绕流的Euler／边界层方程干扰数值解

本文利用Euler方程和可压缩湍流边界层积分方程研究绕跨声速翼型的有粘与无粘强干扰流动。应用有限差分法在贴体的网格上求解时间相关的Euler方程,以剪功积分方法求解翼面贴附和分离湍流边界层流动,并引入一个松弛方程描述剪应力对上游湍流历程的延迟响应。有粘/无粘干扰采用表面源模型。计算结果表明,对翼面存在强干扰流动情况,获得了与实验值基本吻合的结果。     

平面叶栅跨音流动的数值计算及变工况性能的研究

提出一种隐式矢通量分裂差分格式并用来直接求解Ｒｅｙｎｏｌｄｓ平均Ｎ－Ｓ方程组。该方法避开了近似因子分解及矩阵运算,具有精度高、稳定性好、计算量少等优点。在平面叶栅跨音流场的计算中,较好地捕获了激波,与实验比较,结果令人满意。     

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

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

跨音速单转子压气机三维湍流流场的数值计算

用有限体积表示的三维雷诺平均N-S方程显式格式的时间推进法, 对一个单转子压气机的三维湍流流场进行了计算。方法中应用了Baldwin-Lomax湍流模型。这个单转子压气机是第39届国际燃气轮机会议作为考核三维流场数值方法的盲题。计算结果与实验实测的结果进行了比较, 证明本文提出的方法可以用来验证叶轮机械气动设计。      

喷管二维跨声速两相湍流流场的数值模拟

采用ＳＩＭＰＬＥ算法,结合分散颗粒轨道（ＰＳＩＣ）模型,在同位网格上实现了喷管中跨声速两相湍流流场的数值模拟。根据Ｒｈｉｅ,Ｃｈｏｗ提出的动量插值概念,以协变物理速度分量作为求解变量,完成跨声速气相流场的数值计算。湍流模式采用ｋ－ε双方程模型,计算结果与提供的实验数据相符合。     

跨音叶栅流动的多级有限元分析

将各向异性张力有限单元和Taylor-Galerkin有限元法在多级有限元法的思想下结合起来, 构成了收敛速度和稳定性均较好的适用于跨音流动计算的各向异性多级有限元法。利用这一方法, 基于Euler方程, 研究了平面跨音速叶栅流动, 并与实验结果作了比较。结果表明, 本方法是跨音速计算的强有力手段。      

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

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

三维C-H网格下跨音压气机转子Euler流场数值模拟

介绍了三维C-H网格的生成和全三维转子流场的数值解。以法国宇航院设计的一个高负荷跨音转子为算例, 并与他们的计算结果作了比较。应用这套数值模拟技术得以计算带进口预旋的某跨音压气机转子流场, 并对进口预旋对其下游的转子流场的影响作了数值模拟, 得出了若干有意义的结论。      

跨声速轴流压气机级全工况三维流场数值分析

以LV－SGS－GE隐式格式和MUSCL TVD迎风格式为基础,结合壁面函数方法和简单的混合长度湍流模型,对三维可压缩雷诺平均N－S方程进行求解。叶列间参数的传递采用混合平面方法并应用了微机网络并行计算技术。计算得到了NASA 35号低展弦比、跨声速轴流压气机组70％和90％设计转速下的全工况性能曲线,并重点分析了其中一些典型工况下的内部流场。计算与实验结果的对比表明,此方法能快速得到三维粘性流场的流动特性且计算精度较高,可用来模拟跨声速轴流压气机组内的全工况三维粘性流动。     

考虑顶隙的压气机单转子内三维跨音速粘性流场结构的数值模拟

借助于考虑叶片顶部间隙的压气机转子内跨音速湍流流场的三维数值解,得到了该转子流场中的复杂涡系与波系结构。给出了叶片表面及轮毂面上的壁面流型,分析并揭示了转子叶道内的马蹄涡、通道涡、角涡、尾涡等复杂涡系及其与激波的相互干涉;特别是,分析并揭示了转子叶道内顶隙射流的特点,得到了顶隙附近激波与泄漏流动的干涉形态及泄漏流动的发展,并着重给出了其与激波的干涉形态。     

求解跨声速压气机叶栅粘性流动反问题的数值解

阐述了一种以有限体积的时间推进方法为基础求解跨声速压气机叶栅设计的反问题方法。应用这一技术设计叶型,规定叶片表面无量纲目标速度分布,通过比较目标速度分布和计算获得的速度分布来修改叶片压力面和吸力面的坐标,最终获得要求的叶片形状。计算程序使用分布体力方法模拟粘性,用局部时间步长和多重网格方法加速收敛。还介绍了这种方法求解的基本方程系统,并给出了采用这种方法设计的无激波超临界跨声速压气机叶栅的设计结果。