等截面矩形隔离段内流场的三维数值模拟

针对三维矩形等截面隔离段中由激波／附面层干扰诱导的复杂流场，采用有限体积方法，OC－TVD格式结合B－L代数湍流模型隐式求解了三维雷诺平均N－S方程，数值模拟了三维流场内复杂的激波串和旋涡结构，计算结果表明所采用的数值方法能够很好地捕捉隔离段内流场的复杂流动现象。     

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

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

复杂外形气动力计算的阻力修正方法

基于气动力计算的位流线化理论,发展了一种亚音速附面层阻力修正方法。采用联合流场求解方法获得复杂外形的表面压力分布,在此基础上,将绕复杂外形物体沿流向面元划分为片条,在片条上,把三维附面层转换为二维附面层,对二维附面层进行数值求解,得到附面层粘性阻力修正。算例表明,该方法计算量小,与实验和N-S数值解结果吻合较好。为复杂气动外形飞行器初期设计、方案选型及评估提供了较为理想的工程设计手段。     

高速三维边界层稳定性问题的数值方法

本文把四阶精度紧致差分格式用于求解声速、高超声速三维边界层稳定性问题。提出了适用地求解高速边界层特征值问题的Ｃｈｅｂｙｓｈｅｖ多项式配置点方法,并将各种数值方法进行了比较。数值结果表明上述两种方法是求解高速三维边界层稳定性问题的有效数值方法。     

跨声速三维Euler方程的显式／隐式格式有限元解法

本文采用有限体积Ｇａｌｅｒｋｉｎ法和非结构网格数值求解了跨声速三维Ｅｕｌｅｒ方程,通量的求解采用ＶａｎＬｅｅｒ－阶逆风矢能量发裂方式,并通过外插使其上升为二阶精度,分别采用了显式及隐式两种格式求解。     

基于Euler方程结合附面层修正方法的应用研究

采用欧拉加附面层修正的方法在结构网格上对民用飞机翼身组合体流场进行了数值模拟研究.采用一种适合翼身组合体模型的的网格生成技术来生成计算用的结构网格,该技术是先采用一种基于Higenstock源项修正的椭圆型偏微分方程生成二维网格,然后沿机身轴向合理的光滑性调整成三维网格.在求解Euler方程时采用中心离散格式和显式时间推进,并利用附面层修正的方法来考虑粘性的影响.为了验证此方法的正确性,本文对DLR-F4翼身组合体进行了数值模拟,结果表明此方法和试验结果能够很好的吻合.     

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

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

RANS方程和附面层方程耦合求解转捩位置的方法

研究了基于线性稳定性分析的RANS（Reynolds-Averaged Navier-Stokes）方程和附面层方程耦合求解转捩位置的方法。RANS方程计算得附面层外边界的速度分布作为附面层方程计算的输入,求解附面层方程得到稳定性分析需要的附面层参数。由于附面层的解在分离点处的奇异性,用交互式附面层方程求解方法。用基于线性稳定性分析的e^N方法计算转捩位置,并考虑了Tollmien-Schliching波和层流附面层分离造成的转捩。迭代求解RANS方程和附面层方程,直到转捩位置收敛。RANS求解中使用转捩过渡区模型,避免了点转捩模型引入的数值扰动。通过对NA-CA0012翼型的计算并和实验结果及Xfoil计算的结果进行比较,吻合较好,有较好的工程实用价值。     

高阶精度格式WCNS在三角翼大攻角模拟中的应用研究

采用5阶精度的加权紧致非线性格式(WCNS-E-5)数值模拟了65°三角翼的大攻角绕流流场,主要目的是考核高阶精度格式WCNS在大攻角旋涡流动方面以及跨声速流场的激波附面层干扰、涡破裂位置的模拟能力,重点研究不同网格规模和湍流模型对尖前缘三角翼涡系之间的相互作用的影响。通过求解任意坐标系下的雷诺平均N-S方程,采用5阶精度的加权紧致非线性格式(WCNS-E-5)和多块对接结构网格技术,两种湍流模型分别是一方程SA和两方程SST湍流模型,在与相应试验结果对比的基础上,详细研究了WCNS-E-5格式在跨声速大攻角旋涡流动中的表现,以及不同网格规模、两种湍流模型对主涡二次涡相互作用、涡破裂位置和表面压力分布的影响。本文的研究结果表明,高阶精度格式WCNS-E-5能成功应用于三角翼的跨声速大攻角流动,网格规模的增加进一步提高流场分辨率,SST湍流模型相对SA湍流模型在三角翼大攻角流动中具有更好的适用性。     

三维方柱不可压缩绕流的大涡模拟计算

用非定常不可压缩流动求解方法对三维方柱绕流问题作了大涡模拟计算研究.求解程序采用虚拟压缩方法和双时间步的时间推进方案,对流项的差分采用高阶精度的数值差分格式.通过与实验和文献结果比较分析,验证了基于虚拟压缩方法的湍流大涡模拟计算方案对于复杂的三维钝体绕流计算是可行的.     

一种双斜切双压缩面气道性能的数值模拟

采用LU－SGS隐式推进方法及ROE的通量差分分裂格式，结合带非线性限量因子的MUSCL插值方法编写了雷诺平均的三维N－S方程的数值求解程序，并供助该流场计算程序研究了一种双斜切双压缩面进气道在未经流场控制情况下的流动结构和气动性能（M＝1.9,H=11Km）。计算结果表明：所编制的流场计算程序能够准确快速地求解三维粘性流场，并对流场间断具有较高的分辨率；所研究的进气道即使在未经流场控制的情况下也具有较高的总压恢复系数，这体现了双斜切双压缩面进气道总压恢复系数高的优点，但由于流场中存在分离流动，应采取抽吸附面层等流场控制措施。     

附面层分离流的逆算法

所谓附面层的逆算法是指给定位移厚度δ~*或壁面剪切应力τ_w,而把外流速度分布u_e或压强分布p做为未知函数通过附面层方程的求解反算出来。这样就使得附面层差分计算顺利地通过分离点而不出现任何奇性,从而开辟了一条求解分离流动的非常有意义的途径。与直接求解N-S方程相比,它的计算量与贮存量要小得多。对实际应用来说,它是求解附面层分离气泡的好的近似方法。     

应用 TVD 格式求解 N-S 方程预测透平叶栅二次流损失

应用ＴＶＤ格式和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数紊流模型求解三维Ｎ－Ｓ方程,详细模拟出了一个直列透平叶栅中二次流及其损失在叶栅通道中的产生和发展过程,并给出了叶栅损失的分布状况,其结果深刻揭示了透平叶栅二次流损失的产生发展机制。本文计算得到的叶栅损失的栅距平均值沿叶高的分布曲线与实验值吻合良好,表明格式求解精度较高。     

不定常差分方法中的二阶三层格式

本文介绍一种求解三维理想可压流问题的二阶三层格式并给出它的稳定性证明。通过大量二维和三维问题的实际计算,表明此方法的优点是精度好,时间省,公式简单,程序量少,便于实际应用。文中对应用此法如何求解间断流场的问题,针对具体算例做了较详细的讨论。对不同外形在不同流动条件下,如何设计差分网格也做了分析。最后给出部分复杂外形的二维和三维计算结果。     

涡轮叶栅冷气掺混数值模拟方法

采用具有三阶精度TVD性质的有限差分格式、自由型曲面网格生成技术以及分区网格算法,对某型涡轮叶栅进行了任意数目冷气孔的薄膜冷却全三维N-S方程数值求解,较为详细的阐述了冷气掺混流场中卵型涡的发展规律。结果表明,本文所开发的求解系统能够较好的实现任意分区的流场求解,对涡轮叶栅中冷气喷射流场的描述是合理的。      

跨声速压气机数值模拟格式效应研究

针对跨声速叶轮机械复杂内部粘性流动,开发了全三维高精度RANS数值求解器。分别采用FVS格式与AUSM+格式结合多种限制器对Rotor 67跨声速压气机转子进行数值模拟,分析和比较了各格式的计算效果并与试验结果对比。结果表明:与FVS格式相比,AUSM+格式的数值粘性更小,边界层的模拟精度更高;Hemker限制器的综合表现最优,Van Albada限制器的粘性分辨率略低于Hemker限制器,Minmod限制器对流动分离现象的捕捉能力较差,Van Leer限制器容易引入色散误差。     

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

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

涡喷发动机涡轮密封舱小气路流场数值计算

运用 NUMECA的三维粘性流动数值计算软件 FINE/ TURBO对涡喷发动机涡轮密封舱小气路中的复杂流场进行研究 ,计算中采用 Jameson的中心差分格式结合 Baldwin- Lomax湍流模型求解雷诺平均 N- S方程 ,得到了三维复杂流场中气体参数分布的详细结构和规律 ,具有较高的精度 ,为小气路中的故障的诊断提供数值依据     

三元Navier-Stokes方程的数值计算——激波与附面层的相互干扰

前言 对于空气动力学所提出的问题,通常所需要的只是其定常解。稳定法被广泛地用来解决这类问题。显式格式更为简单,但稳定性要求很小的时间步长,因而计算机时间题是一个突出问题。为克服这一问题,近年来人们作了很多工作。本文采用了一个单步差分格式来计算三元激波与附面层干扰问题。此方法简单,在解的光滑区,格式具有形式上的二阶精度,并且格式是高度稳定的。计算中另一个问题是,在激波附近的     

机翼附面层层流区的积分计算方法

本文提供了一种精度较高的层流附面层积分计算方法。它改进了国外机分法的计算精度,且不受附面层外流必须是位势流的限制。