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

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

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

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

双曲守恒型方程的二阶摄动有限差分格式

对双曲守恒型方程，将其一阶迎风格式空间差商的常系数摄动展开为时间步长和空间步长的幂级数，通过确定幂级数系数而获得二阶精度的摄动有限差分(PFD)格式。进而从双曲守恒型方程的通量分裂型一阶迎风格式出发，通过类似的摄动展开方法，获得空间精度为二阶的通量分裂形式的摄动有限差分(FPFD)格式。这两类格式保留了一阶守恒迎风格式的简洁结构形式，使用三节点即可达到二阶精度，又避免了三点二阶格式的非物理数值振荡。并将这两类格式推广应用到双曲守恒型方程组，最后通过模型方程和一维激波管流动的数值算例验证了格式的高精度、高分辨率性质。     

中层及低热层大气对局地时变热源的响应

采用二维交替方向隐格式(ADI),数值模拟了中高层大气对一个位于对流层高度范围的局地时变热源的响应．分析了这一过程中大气响应的能量密度,水平波数,垂直波数以及频率．模拟结果表明:响应波动携带能量向上传播,集中在较为狭窄的水平区域,该区域关于源的水平位置对称;相对于动量源驱动,热源的水平尺度对响应的水平波长决定性作用更加明显;对于高频振荡的热源,其响应的频率和源一致;响应的垂直波长(41 km)比源的垂直尺度(6 km)大得多．     

广义拉格朗日乘子MHD方程组的两类逆风分裂格式

开发高性能的磁流体力学数值模拟方法是提高空间天气数值预报研究的一个重要方面.有限体积法的逆风分裂格式具有良好的间断捕获能力,Steger-Warming和AUSM（Advection Upstream Splitting Method）是逆风分裂格式FVS（Flux Vector Splitting）方法中具有代表性的两种格式.采用这两种格式求解具有伽利略不变性的扩展型广义拉格朗日乘子磁流体力学（EGLM-MHD）方程组,对Orszag-Tang涡流问题和三维爆炸波问题进行数值模拟,结果表明两种格式均能得到稳定精确的数值结果.与Steger-Warming格式相比,AUSM格式产生的磁场散度误差更小,计算速度更快.     

跨音速机翼设计方法的改进

以Takanashi提出的"正反迭代、余量修正"的原理为基础,针对其当机翼上表面存在较强的超音速流动时,解的收敛性不好的问题,对Takanashi方法进行改进.具体方法就是在控制方程的双曲部分引入基于迎风格式的附加项.通过对某干线飞机机翼的设计,证明了改进后的方法是成功的,设计结果在超音区有很好的收敛性.      

面向混合网格高精度阻力预测的熵修正方法

Roe通量差分分裂格式具有耗散小、分辨率高等特点,在亚跨超声速流场模拟中得到广泛应用,但空间离散时需要对Roe平均矩阵的特征值进行熵修正,使用熵修正会增大耗散,影响阻力的求解精度。本文针对非结构混合网格中Roe格式熵修正的特点,通过改进传统的Roe格式Harten-Yee熵修正方法,提出了一种可提高非结构混合网格黏性计算精度的Harten-Yee熵修正改进方法。利用改进后的方法,完成了DLR-F4翼身组合体算例的计算和对比分析,改进后的熵修正方法残差收敛特性与原始Harten-Yee熵修正一致,计算精度提高,计算结果和无熵修正时基本一致,说明改进后的熵修正方法既保留了使用熵修正带来的程序鲁棒性等优点,同时把熵修正对阻力预测精度的影响降到了最低,验证了方法的有效性。采用改进后方法对第3届AIAA阻力会议的计算模型DLR-F6翼身组合体进行了详细的模拟,分析了网格收敛性和雷诺数的影响。结果表明：改进的Harten-Yee熵修正,更加适用于非结构混合网格的黏性计算,计算精准度达到国际同类CFD软件水平,进一步验证了改进方法的可靠性。     

Kirchhoff型抛物方程的Galerkin有限元法的超收敛误差分析

文章研究了后向Euler全离散Galerkin格式下的Kirchhoff型抛物方程的超收敛误差分析。首先,讨论了数值解的先验误差估计,并证明了数值解的存在唯一性。其次,使用双线性元的高精度误差估计以及Ritz投影算子与插值算子相结合的技术,通过技巧性地处理非线性项得到了超逼近的误差估计结果。再次,通过插值后处理方法获得了整体的超收敛结果。最后,通过数值试验验证了理论分析的正确性。     

用解析离散法构造 WENO-FCT 格式

本文利用解析离散法的导数计算公式构造高精度ＷＥＮＯ格式,并引入ＦＣＴ通量修正技术以消除非物理解振荡。通过激波－边界层干扰,激波－旋涡相干结构的数值模拟证明方法既能够有效地捕捉激波间断,又对边界层分离及旋涡结构有很好的分辨率,是一类有效的高精度计算格式。     

矩形/非结构混合网格技术及在二维/三维复杂无粘流场数值模拟中的应用

建立了一套模拟复杂无粘流场的矩形／非结构混合网格技术,其中非结构网格仅限于物面附近,发挥非结构网格的几何灵活性,以少量的网格模拟复杂外形;同时在以外的区域采用矩形结构网格,发挥矩形网格计算简单快速的优势,有效地克服全非结构网格计算方法需要较大内存量和较长ＣＰＵ时间的不足。混合网格系统由修正的四分树（２Ｄ）／八分树（３Ｄ）方法生成。本文将ＮＮＤ有限差分与ＮＮＤ有限体积格式有机地融合于混合网格计算,消除了全矩形网格模拟曲面边界的台阶效应,同时保证了网格间的通量守恒。该算法被成功地运用于ＮＡＣＡ００１２翼型的跨、低超声速绕流,二维弯管内运动激波的非定常流,以及三维捆绑火箭的超声速绕流的数值模拟。计算结果表明本方法在模拟复杂无粘流场方面的灵活性和高效性。