超声速进气道可压及不可压流动数值模拟

采用Chio-Merkle预处理矩阵对可压NS(Navier-Stokes)方程组进行时间预处理,分析了预处理方法的物理和数学背景.用有限体积方法,结合LU-SGS(Lower-Upper-Symmetric-Gauss-Seidel)隐式时间积分和AUSM+(P)(Precondition Advection Upstream Splitting Method)格式、中心差分2种空间离散格式,求解该预处理NS方程组.通过圆弧凸包流动和二维方腔顶盖驱动流动的数值试验表明,该方法克服了传统时间迭代方法模拟低速流动时的刚性问题,加速了收敛过程,可以同时有效地模拟可压及不可压流动.用此方法与块结构化网格技术相结合,进行了在不同飞行马赫数、攻角和出口反压条件下,混压式轴对称超声速进气道三维流场数值模拟.计算结果表明:该方法能准确地捕获复杂激波系,清楚地揭示进气道三维流动现象;进气道性能随工作条件的变化规律,与理论分析一致.      

复杂外形下的无网格算法研究

研究了二维、三维复杂外形下的无网格算法。在一种布置点云方法的基础上,发展一种曲面拟合的重构方式构造流场物理量,并应用于AUSM+_up格式计算欧拉方程的数值通量;计算采用了隐式时间推进,并引入当地时间步长和残值光顺等加速收敛措施。通过对某三段翼型低速流动、M6机翼跨音速流动、某全机跨音速流动进行了数值模拟,表明本文发展的无网格算法能有效地模拟复杂外形的无粘绕流。     

激波绕楔形物体的流场分析

使用数值计算方法和数值干涉技术并结合实验干涉图片,分析了激波过拐角绕楔形物体运动的流场.并采用数值彩色干涉技术,将数值计算结果转换为彩色条纹编码的数值干涉图,并与实验干涉图比较,可以辨识过激波发生跳跃后的干涉条纹的位移数,从而确定实验干涉照片的密度场.数值计算采用显式二阶迎风TVD格式.采用数值干涉技术将计算密度场转化成数值干涉图与实验干涉图进行了比较,检验了数值结果的可靠性.数值计算结果显示了包括激波、接触间断和旋涡在内的详细的流场结构,与实验相比,更为细致地刻画了流场演变过程的细节.通过数值计算结果观察到激波传播到旋涡中心附近并没有发生明显耗散,旋涡上半部的激波得到减速,下半部分的激波得到加速.     

非结构网格二维理想MHD绕流逆风格式解法

在非结构网格上发展了针对二维理想磁流体方程组的逆风格式求解方法.控制方程中的对流项采用AUSM格式处理,时间推进采用显式5步龙格-库塔方法.为了消除计算中产生的磁场散度的影响,引入了双曲型散度清除方法.通过对磁流体激波管问题的求解验证了该方法对激波的捕捉能力,对有均匀磁场干扰下的喷管流动情况进行了数值模拟,并与文献中结果进行了对比.计算结果显示了磁场对磁流体流动的干扰效应,该结果与参考文献中的数值模拟结果相吻合.     

几种典型迎风格式的分析与比较

重点分析了AUSM 格式的构造方法,同时将其与Van Leer, Roe格式等构造方法进行比较和分析,通过求解一维激波管问题,以及二维和三维跨音速流场的计算比较了AUSM 格式,VanLeer格式,Roe格式和中心格式的数值模拟结果和网格适应性,以及在经典四步Runge-Kutta显式时间推进和LU-SGS隐式时间推进下上述四种格式的计算效率.获取了对于基础研究和工程应用很有意义的结论.     

自适应的高精度中心﹣WE NO混合格式

为对可压缩流动进行高精度高分辨率的数值模拟,基于模板的光滑度量函数,构造了新型的加权因子,在此基础上将7阶精度的WENO格式与8阶中心格式进行加权组合,构造了一种自适应低耗散的中心-WENO混合格式(H-WENO7-CD8),并采用Fourier方法对离散格式数学特性进行了理论分析。该自适应的高精度中心-WENO混合格式相比于七阶的WENO格式具有更小的耗散误差。通过对激波-密度干扰问题、Rayleigh-Taylor不稳定性问题和双Mach反射问题的数值模拟,并与WENO-JS格式的计算结果对比,结果表明:该格式结合了WENO格式和中心格式的特点,能更好地捕捉激波,对物理脉动也具有较高的分辨率,适合于复杂流体流动的数值模拟。     

面向多CAD平台的航空紧固件数据库型式标准研究与实现

分析了当前零件库的局限性和紧固件数据库型式标准的需求,研究了事物特性表设计、中性零件模型表达和零件簇构造等关键技术,开发了中性模型生成工具、紧固件管理系统、CATIA/UG/STEP接口等软件,实现了"数据即标准","一次建模、多次使用","所见即所得","过程特征完整保留"等应用特点。     

采用改进的CE/SE方法模拟方管中氢氧爆轰波的稳定传播结构

采用改进的时空守恒元/解元（CE/SE）格式和新型二步化学反应模型，对小扰动情况下方管通道中的氢氧气相爆轰传播过程进行了数值模拟，对不同截面尺寸和两种扰动模式下的爆轰波三波线结构进行了详细计算和分析。数值模拟花费了相对较少的计算代价，清晰地得到了方形截面通道中3种相对稳定的爆轰波传播结构，即直角模式、对角模式和螺旋模式。随着截面尺寸的减小，爆轰波胞格图样与波阵面结构随之产生适应性的变形，乃至出现胞格数目的减少和部分三波线的消失；最终，直角模式和对角模式都会转化为顺时针或逆时针的螺旋模式；临界尺寸下，得到了转换过程中新的过渡结构和压力脉冲振幅的变化，并以此初步讨论了直角和对角传播模式的稳定性。     

二维Euler方程组的高分辨率隐式差分格式的进行计算

本文首先对标量双曲型守恒律构造了一类新的二阶隐式ＴＶＤ差分格式，再利用Ｒｏｅ方法和隐式近似因子分解法推广格式到一维、二维守恒率方程组。为减少计算工作量，推导了二Ｅｕｌｅｒ方程组隐格式的对角形式。最后，针对江南－Ⅱ并行机设计了适合盱多处理机结构特征的并行计算程序，计算翼型的跨声速流场，结果令人满意。