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

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

高超声速飞行器加速度导数数值模拟

在Etkin的非定常气动力模型基础上,采用非定常Navier-Stokes方程描述飞行器的非定常强迫沉浮振动流场,采用空间二阶精度的NND格式和"双时间"时间离散法进行数值求解,研究飞行器加速度导数的数值计算方法。并给出了返回舱外形和HBS(Hyper Ballistic Shape)标模外形的加速度导数、旋转导数和俯仰阻尼导数的计算结果及分析。     

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

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

返回舱自由振动跨声速非定常流场数值模拟

本文利用NND差分格式对返回舱俯仰运动和跨声速绕流流场进行了耦合求解，结果表明来流M     

非结构网格、混合网格下计算方法研究

为了提高计算精度，本文研究了计算格式和网格技术。通过比较非结构网格下的NND格式、MUSCL格式、ENO格式及最小二乘格式的计算精度和波动，选择利用NND方法中的守恒变量法来建立非结构网格下的计算格式，通过改进NND方法中的限制器来提高计算精度和分辨率，计算结果比较显示出该格式具有很高的分辨率和计算精度。本文也研究了非定常网格自适应技术，计算给出了双马赫反射问题的新旧限制器及网格自适应计算结果比较，新限制器网格自适应后计算得到的滑移线出现了剪切涡现象，这表明，数值计算中取足够密的网格和构造高精度格式是非常必要的。     

声速再入体表面热流数值模拟研究

本文采用计算效率较高的标量对角化隐式NND格式，通过求解Navier-Stokes方程对影响再入体表面热流计算准确的诸因素进行了综合分析。研究了Steger-Warming通矢量分裂、VanLeer通矢量分裂和通量差分分裂方法及相应熵修正方法对热流的分辨能力，并阐明了在物面边界上采用二阶中心格式、二阶中心和二阶迎风混合格式、以及一阶迎风格式等不同边界格式对热流计算的影响。在此基础上，采用通量差分形式NND格式对钝锥和钝双锥高超声速粘性绕流进行了数值模拟，计算给出了与试验结果相吻合的热流分布。     

尖拱弹身大攻角N-S方程数值计算

本文采用有限体积的NND格式,数值模拟了弹身超音速有攻角绕流,分析了周向网格和法向网格对弹身背风区分离流态及气动性能的影响,计算与实验结果比较,两者相当一致.     

无波动,无自由参数的耗散差分格式（NND）在喷流计算中的应用

本文用文献[1]建立的无波动、无自由参数的耗散差分格式(NND格式),对无粘喷流流动进行了数值模拟,并就一典型的自由射流流动,与国外高精度、高分辨率的PLM差分格式(Piecewise-Linear Method)的计算结果和实验结果作了比较,得到了很一致的结果,且所用的计算时间少于PLM格式。     

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

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

烧蚀外形方程差分计算方法研究

分析了烧蚀外形方程的性质,研究了人工粘性项对烧蚀外形计算结果的影响,发现人工粘性在抹平非物理波动的同时,会严重影响计算精度.借鉴于CFD中NND格式的思想,构造了求解烧蚀外形方程的无波动、无自由参数的耗散差分格式,得到了满意的计算结果.     

时间模式下二维可压缩气固两相混合层流动研究

采用Euler-Lagrange颗粒-轨道双向耦合模型对时间模式下含有固粒的二维可压缩混合层流场进行了研究。气相流场采用具有空间三阶精度的WNND格式进行直接数值模拟,固相方程采用单边三点差分离散。在考虑流场对固粒作用的同时,也计及固粒对流场的反作用。在对流马赫数为0.5时,研究了颗粒相对密度、颗粒尺寸以及Stokes数等对粒子运动和流场结构的影响。     

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

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

超音速来流中侧向喷流干扰流场的数值模拟

通过求解三维可压缩N S方程组,对尖锥/圆柱组合体模型在超音速绕流中的侧向喷流干扰流场进行了模拟,数值方法采用有限体积法,差分格式采用NND格式。分析了干扰流场的结构,研究了干扰效应对气动性能的影响,得到来流Ma∞=3 3和4 5时力放大因子和干扰力矩系数随攻角变化的规律,计算结果与实验一致。     

基于PSE的单股剪切混合流稳定性分析

将抛物化稳定性方程(PSE)方法应用到可压缩单股剪切混合流的稳定性研究中.采用并发展了适用于自由剪切流的高精度数值方法,包括六阶紧致格式、坐标变换以及渐近边界条件等,对PSE进行有效求解.通过求解相似边界层方程得到更准确的剪切层内基本流;求解线性稳定性理论(LST)控制方程得到扰动的初始条件,并通过流向空间推进方法对扰...     

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

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

DCS5在计算气动声学中的应用研究

对流项采用五阶线性耗散紧致格式(DCS5)离散,粘性项采用6阶中心紧致格式离散,时间项采用5级5阶龙格-库塔时间推进格式求解流体动力学方程组,并由此建立近声场的直接模拟方法。为了应用结构网格在复杂几何形状条件下对流场进行计算,发展了高精度的多块特征对接边界技术。完成了低速无粘振荡圆柱和粘性圆柱绕流等标准计算声学算例,获得了令人满意的结果。     

紧致格式数值模拟超音速粘性绕流问题

通过数值求解可压Navier-Stokes方程的方法计算了球锥粘性绕流问题。方程中粘性项按通常办法逼近,为改善收敛速度利用了文[1,2]中之算子附加修正方法。方程中之流向导数利用了四阶精度的紧致差分。这是由于该方向上流动参量光滑和边界条件易于处理的缘故。在法向方向上利用了中心差分。对差分方程之隐式部分利用了近似因式分解法。由于选用了文[3]中之特殊Jacobiall系数矩阵分裂法而使计算工作量大为减少。用这种方法计算了超音速粘性绕流问题。计算结果与实验结果相符很好。     

先进直升机旋翼悬停状态气动性能计算

 为更好地模拟先进直升机旋翼流场,同时准确计算其悬停状态的气动性能,建立了一套基于Navier-Stokes/Euler方程的混合CFD方法。该方法的求解域由两部分组成：一是围绕旋翼桨叶周围的黏性区域,采用可压缩Navier-Stokes方程来模拟旋翼附近的黏性流动和近场尾涡的捕捉;二是离桨叶较远、黏性可以忽略的远场区域,用Euler方程来描述其流动。在该方法中,将三阶逆风格式（MUSCL）与通量差分分裂方法相结合,无需添加人工黏性,因而可有效地减少旋翼尾迹数值耗散。?阌诹鞒》智蠼庖约爸芷谛员呓缣跫氖凳?采用了嵌套网格方法,并给出旋翼网格与背景网格交界面的信息传递方式。应用所建立的计算方法,首先对二维翼型、三维M6机翼的流场进行了数值模拟,以验证计算方法;然后,着重计算了有实验结果可供对比的具有先进气动外形的HELISHAPE 7A模型旋翼和UH-60A直升机旋翼,通过计算得到旋翼表面压力分布、桨叶展向拉力系数分布、桨叶表面细节流动以及气动性能等,进一步验证了该方法的有效性。     

高温高压超音速气液两相流场数值模拟

采用欧拉-拉格朗日方法对燃气-蒸汽发射动力装置内高温高压超音速燃气中的横向喷雾的气液两相流进行了数值模拟研究，气相采用RK-AUSM+格式求解N-S方程，液相应用颗粒轨道模型，两相之间的耦合通过在气相各守恒方程中添加源项来实现。针对不同的喷水孔径、喷水压差以及水汽质量比进行了数值模拟计算，分析了各种条件下的液滴蒸发情况以及对气相流场和能量输出的影响，为燃气-蒸汽发射动力装置的设计优化和能量调节技术提供了理论依据     

绕非对称头锥的分离系数矩阵差分法

一、亚跨区SCM方法 （一）基本方程和边界条件 建立在球坐标下无因次化的拟非定常欧拉方程组为