基于Steger-Warming矢通分裂的NND格式

将Steger-Warming矢通分裂方法和NND格式相结合构成一种基于Steger-Warming矢通分裂且能自动捕捉到激波的差分算法.通过对激波管和对可燃气体DDT过程的数值分析,证明该算法解决了Steger-Warming矢通分裂方法在零特征点附近存在数值解震荡的问题,有着良好的计算精度和计算效率.文中采用的计算方法对各种有关矢通分裂计算方法的实现有一定的参考价值.     

NND格式在多维理想磁流体方程组中的应用

采用修正的四步Runge-Kutta方法求解三维一般曲线坐标系下的理想磁流体方程组,为克服数值振荡,加特征型NND格式进行后处理.特征型NND格式推广到求解三维磁流体(MHD)问题需要知道雅可比通量的左右特征矩阵,在具体计算时需要克服矩阵的奇性.本文用三维程序采用推广的特征NND格式计算了一维MHD激波管和二维(MHD)喷管流动,计算结果表明,特征NND格式保持了TVD格式高精度的优点,又具有计算简单的特点,在包括强弱间断等复杂波系的定常和非定常MHD流场数值模拟中是成功的.     

粘流与无粘流的相互作用计算

本文总结了粘流/无粘流的各种计算方法和结果。重点在于介绍定常流动中的弱相互作用。首先叙述了弱相互作用的数学模型。给出了不可压流动和跨音速流动中粘流/无粘流相互作用的某些正耦合的计算结果。讨论了在分离区附近边界层正方法失效的原因。然后介绍了边界层反方法和适用于带分离的流动中半反方法耦合的粘流/无粘流的相互作用方法。文中也简单地总结了三维情况的应用和强相互作用。     

非结构网格上新型的NND有限元格式

张涵信发展的ＮＮＤ差分格式是由中心差分格式,二阶迎风格式和一阶迎风格式混合组成的杂交型格式。众所周知,和中心差分格式相对应的是Ｇａｌｅｒｋｉｎ有限元格式。通过对中心型有限元格式加修正项的方法本文成功地构造出二阶迎风型有限元格式和一阶迎风型有限元格式。     

关于NND格式的两点注记

对差分格式在激波上计算准确性的标准提出了新的看法,并证明NND格式完全符合这个标准.     

三维坑湍流分离流动的粘流／无粘流干扰计算

本文对坑的三维分离流动做了低速粘流与无粘流的相互作用计算。对三维边界层反方法进行了分析和讨论。用数值试验的方法验证了在H和α作为已知量的情况下,三维边界层反方法的积分方程是双曲型的,并提出了一种近似数值特征线法进行求解。无粘流采用低速位流面元法。计算表明所用方法可计算出三维效应很强(即横向变化很大)的三维边界层分离流动。     

高精度差分格式在机翼分离流计算中的应用

用一个模型方程分析了边界层方程的数值稳定性 ,指出稳定性问题会随计算雷诺数减小而变得严重 ,因而对较低的Re数 ,计算分离边界层流动采用高阶精度差分格式十分必要。本文同时给出用四阶精度差分格式求解在流动分离情况下边界层方程的技术方法 ;对一个大展弦比后掠翼在攻角等于 1 4°时计算了机翼上表面的分离线 ,并与文献上公布的用别的测算方法测出结果做了比较 ,两者总体上符合良好     

跨音速粘流的计算

 概述了一种采用粘流/无粘流相互作用原理计算跨音速粘流的方法。采用熵修正激波算子计及跨越激波时的熵增,形成的非等熵位势方法可比传统的位势方法更准确地计算无粘流动。提出了一种流向速度型,结合其它辅助关系式导出了三维湍流边界层积分方程反方法。用此方法可求得粘流解。利用半反耦合方式耦合了无粘流和粘流解。数值算例表明,计算结果与实验结果吻合;且对计算机的要求较低。     

一类中心型三阶格式及应用

通过分析格式的耗散关系和重线性化技术,发展相应格式构造技术,提出了一类中心型的三阶非线性格式,格式具有与NND格式相当的网格模板.由于格式的线性形式具有很好的耗散特性,同时重线性化技术的使用使格式在光滑区尽量采用线性形式,从而理论上提升了计算对旋涡等结构的分辨和保持能力;另一方面,格式的非线性机制使格式在非光滑区退化为类NND格式,实现耗散和稳定机制.验证计算表明,格式能够开展粘性和旋涡及脉动流场的计算应用,算法形式简单,便于工程应用,容易进行有限体积推广.     

跨音速叶栅流场的无粘-边界层迭代计算

二维跨音速叶栅流场的计算是在无粘时间相关法计算程序的基础上采用无粘-边界层迭代而得。马赫数从亚音速到超音速范围的几种数值试验均与测试结果吻合一致。在边界层无强烈脱流时可获得很好的结果。     

粘性高超声速激波—激波碰撞的计算和分析

本文用特片变量型的隐式ＮＮＤ格式计算了五种类型的高超声速激波－激波干扰流场，得到了清晰的流场结构图画。激波和剪切层等分辨得很清楚。计算的物面压力与实验甚为吻合。借助于激波极曲线，对第４类激波干扰结构进行了分析，解释了干扰引起物面压力增高的原因。     

GRP格式推广用于多维复杂流动的数值计算

采用算子分裂技术，将原始求解一线流动问题的高精度、高分辨率GRP格式推广用于包含各种不同间断解及其相互作用的多维复杂流动计算。以轴对称管端出口处激波绕射传播为例，作了数值模拟，得到与实验结果吻合很好的流场精细结构。     

粘性流场三维NS方程有限体积法数值解

提出了有限体积法中的隐式ＮＮＤ格式数值求解三维全ＮＳ方程的薄层三维ＮＳ方程的方法，对绕带座舱的再入飞行器高超声速粘性流动进行数值模拟。     

飞船高超声速粘性绕流的数值模拟

本文采用稳式ＮＮＤ格式，通过求解三维薄层近似Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，数值模拟了飞船高超声速粘性流场。在计算中充分考虑了网格雷诺数和网格分布对计算结果的影响，不但得到了满意的表面压力、热流分布和流场的物理量分布，而且得到与实验观察一致的表面油流图像和分离结构。     

粘性不可压缩流动三维复杂流场分块耦合求解

提出了一个粘性不可压缩流动三维复杂流场的分块耦合求解方法。可解决科学和工程计算中的各种稳态和瞬态粘性不可压缩流动计算问题。     

预估校正的高分辨率格式及其应用

本文从基于通量分裂的一阶迎风格式出发,构造了一类校正型具有高分辨瘁的二阶全变差递减(TVD)格式,并用此格式进行了一维激波管问题及二维问题的计算。结果表明格式的解为二阶精度,激波处无明显的振荡,激波过渡区为1～2个点。     

运动激波绕射过程中的粘性效应

应用 TVD格式求解二维非定常全 N- S方程 ,以研究运动激波绕射现象。考虑粘性后 ,绕凹角流动时 ,压缩拐角附近表面压力分布与实验结果更为符合 ;绕凸起物时 ,出现反射激波因边界层干扰根部加宽 ,前、后角点附近有局部分离与再附 ,并引起局部压力上升 (与无粘流场模拟相比 )等现象。结果表明 ,全 N- S方程可较好地模拟实际粘性流中复杂激波障碍物干扰流场特性     

跨声速轴流压气机三维黏性流场的数值模拟

改进了动静叶级间的掺混面处理方法,提出了用参数分布修正方法来保证参数在掺混面上的周向不均匀性,并使用亚松弛的参数分布系数来解决程序在动静叶轴向间距较小时的计算稳定性和收敛性问题。采用经改进的NAPA软件对NASA37级压气机的设计点和在不同转速下的流场进行了数值模拟,计算结果与试验结果吻合较好。     

含内嵌激波流场的自适应网格求解

为了精确求解含内嵌激波的流场，应采用自适应网格技术。本文用Ｂｒａｃｋｂｉｌｌ和Ｓａｌｔｚｍａｎ的基于变分原理的自适应网格生成方法得到了在“人造斜激波”附近自动加密的网格，然后将该方法用于圆柱体底部流场的ＮＳ方程的求解。计算结果表明，用自适应网格求解比用固定网格求解捕捉的激波宽度较窄，质量较高。     

喷流流场旋涡激波干扰的数值模拟研究

本文采用隐式CSCM格式,通过数值求解NS方程,对平板上二维超声速主流中横向喷流流场进行了数值模拟,并与实验和其它计算结果进行了分析比较,结果令人满意。