用非结构直角网格和欧拉方程计算飞机绕流

 采用八叉树结构 ,生成复杂外形绕流计算的非结构直角网格。物面附近用投影方法 ,使网格贴体。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。对歼击机模型的绕流计算表明 ,网格生成的机时花费很少 ,总体质量好 ,因而欧拉方程解算的收敛质量也好。     

H型叶栅通道流计算网格的改进

提出一种壁面正交、尾迹区网格节点分布逐渐均匀化叶栅通道网格生成方法。为了研究壁面处网格正交性对数值计算结果的影响，推导出壁面附近流动参数一阶导数离散误差与网格线交角的数学关系。采用N-S方程分别对平板紊流附面层流和叶栅通道流进行计算分析。结果表明：壁面处网格正交性对附面层内流动计算预测精度影响较大，而对势流区流动影响很小；尾迹区切向网格节点逐渐均匀化分布能更准确计算尾迹区流动。     

多重网格法求二元跨声速流在一种半自适应网格上的欧拉解

本文采用MacCormack两步显式格式,用有限体积法求解了二元跨声速欧拉流。推导了物面边界条件,采用了特征远场边界条件及远场环量修正。利用保角变换方法生成O型贴体网格,并修改得到了一种在激波处局部加密的半自适应贴体网格。采用多重网格及焓阻尼加速收敛技术计算了NACA0012翼型的跨声速气动特性,得到了十分满意的收敛过程和计算结果。     

基于离散曲率计算的三角网格模型优化调整

采用面积夹角加权的三角网格模型顶点法矢及三角片质心权值对Taubin的三角网格模型离散曲率计算方法进行了改进,在此基础上提出了一种新的三角网格模型优化调整方法。用该方法调整三角网格模型,在模型上曲率变化较平缓的平坦区域及曲率变化较剧烈的特征区域都能取得较好的调整效果。     

战斗机弹射救生系统数值模拟研究

对弹射救生系统这种几何外形复杂、不规则非流线型钝体,以Jameson有限体积法为基础发展了一套基于高质量混合网格的Euler方程解算器,在Euler方程解算器基础上,采用S-A湍流模型发展了一套同样基于混合网格的N-S方程解算器,初步建立了适合弹射救生系统这类几何外形复杂物体的数值模拟技术平台,并采用多种加速收敛措施和合适的耗散项模型使得所发展的解算器在软件实现上鲁棒性很强。通过选取合适的计算起始迎角,对弹射救生系统成功地进行了大迎角和大侧滑角绕流计算,获得了与风洞试验比较吻合的计算结果,并且可以清楚地观察流场中马赫数分布、弹射救生系统各部位压力分布等流场细节,对弹射救生系统空气动力特性分析具有重要的指导意义。该计算平台可作为风洞试验的合理补充、加强和完善。     

求解Euler方程的无网格与有限体积混合算法研究

研究了无网格方法与有限体积法相结合的、求解Euler方程的混合算法.与单一的无网格方法相比较,由于在大部分计算区域采用了有限体积法,使得发展的算法在运算效率上能与有限体积法相当;同时在物体附近嵌入了无网格区,使得在处理几何外形上更具灵活性.有关算法涉及的无网格与网格区域搭接处理,通过在交界面处引入辅助卫星点和网格单元构成边界信息,实现了区域间的流动信息传递.Euler方程的空间导数分别在两区域用有限体积法和无网格法离散近似,时间方向都采用四步显式Runge-Kutta格式推进求解,数值模拟了喷管内流和绕翼型外流,并分别与整体有限体积法和整体无网格方法进行了比较.算例展示出,提出的混合算法能有效捕捉激波间断,且两区域等值线过渡光滑,算法效率如预期与有限体积法相当,表明该方法是可行的.     

三维分区对接网格Euler方程的并行数值模拟

本文根据CFD中分区对接网格的特点，应用基于工作站集群系统的MPI网络并行环境，实现了典型双三角翼模型亚声速绕流Euler方程的并行计算。文中对计算结果的收敛情况、耗时、加速比及效率等并行性能指标进行了讨论。结果表明，多域多机与单域单机的计算结果基本吻合；2节点并行效率达到88．21％，4节点并行效率达到78．99％。     

基于细胞自动控制算法的硅各向异性湿法刻蚀物理仿真

应用细胞自动控制算法对硅湿法刻蚀进行物理仿真,按硅原子晶体结构构建细胞单元,以原子共价键连接状态和相邻细胞的状态来判断刻蚀中某个细胞是否去除或者保留,最后获得与理论分析相符的仿真结果.     

在非结构动网格上追踪多介质界面的ALE方法

在处理带有介质分界面的多介质流动问题时,把介质边界当作一种特殊的边界,并且在该边界上的点都有着双重的定义,分别对应边界两侧的不同介质状态,通过利用HLLC方法求解任意拉格朗日-欧拉方程(ALE),以及运用界面跟踪和动网格方法分别对一维和二维多介质流进行了数值模拟,并得到了满意的结果.     

非结构动网格技术及其在超声速飞行器头罩分离模拟中的应用

针对激波诱导物体运动的二维动力学模型,比较了采用流动定常、结合非定常弹道方程的所谓"准定常"计算方法和非结构动网格计算方法,通过流场和运动特性分析,指出发展有相对运动的非定常流动模拟方法的重要性;另外采用修正弹簧近似和局部重构的非结构动网格技术,耦合求解Euler方程及6DOF弹道方程,数值模拟高度H=25km、飞行马赫数Ma=6.0、迎角α=0°状态下火箭抛罩的动力学过程.研究表明,稠密大气层内头罩受到的气动力比惯性力大几个数量级,导致头罩的运动特性完全由气动力控制.分离过程中,头罩和火箭头部之间的气动干扰在时间和空间上都表现出强烈的非线性特征,头罩姿态变化剧烈.本文还对分离过程中部分位置上的流场进行定常模拟,比较了定常和非定常的气动力,所得结论对于分析地面风洞可控轨迹系统(CTS)试验具有指导意义.