点隐式龙格-库塔方法的应用研究

为了提高求解器的效率,在显式龙格-库塔时间推进的欧拉方程求解器之上,发展了点隐式龙格-库塔时间推进格式。给出了其推导过程和非结构网格下中心格式和迎风格式（Roe和Van Leer格式）预处理算子的构造方法。NACA0012翼型和RAE2822翼型的跨音速无粘流动模拟表明：与显式龙格-库塔方法相比,方法能提高求解效率且内存需求相当,具有一定的工程应用价值。     

结合响应面法与网格法的低速翼型优化设计

本文首先了介绍响应面法和网格法的基本原理,翼型绕流场分析使用有限体积法空间离散和RUNGE-KUTTA时间迭代格式进行N-S方程计算,由翼型表面压力和摩擦力数值积分得到升力和阻力.然后,以NACA0012和GAw-1为算例,结合响应面法和网格法进行低速翼型的减阻优化设计,算例表明:在应用响应面法进行低速翼型优化时,采用...     

翼型表面结冰准定常数值模拟

对翼型表面的结冰过程进行了准定常数值模拟。在每个时间步长内完成网格随着壁面边界的移动而更新、周围流场和水滴撞击特性重新计算、冰形计算及壁面边界的重构工作,如此循环直至所需的结冰计算时间。采用拉格朗日轨迹追踪法获得水滴撞击特性,应用考虑壁面粗糙度影响的边界层积分法计算壁面的对流换热系数,在此基础上求解Messinger结冰热力学模型,冰层始终沿壁面外法线方向生长。对NACA0012翼型在不同环境下生成的三种典型冰形进行了预测,并与实验结果进行了比较,表明本文所述方法的有效性。     

拟压缩性方法数值模拟绕弹体的不可压缩层流脱体涡流动

本文应用在连续生方程中增加压力对时间导数项的拟压缩性方法，数值求解三维定常不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程。数值求解过程采用隐式的Ｂｅａｍ－Ｗａｒｍｉｎｇ近似因式分解格式。计算绕尖头旋成柱体的大迎角脱体涡流动，得到与实验相符的结果。     

二维Euler方程SCB格式的并行实现

本文讨论二维双曲型守恒律方程组的高分辨ＳＣＢ格式，并应用这些格式求解二维Ｅｕｌｅｒ方程组。计算过程中在江南－Ⅱ型并行计算机上进行的，计算过程运用了许多并行计算技术，计算结果令人满意。     

对提高统计线性化精度及高阶矩传播法的研究

提出提高统计线性化精度的新方法,同时研究高阶矩传播法的计算方法,并以此方法做为衡量计算精度的标准,计算结果表明,提出的修正方法改善了计算精度,并且易于推广使用。     

推力级数、推力比和推力作用时间对弹道散布的影响

用详细的计算数据，在保持计算条件和总冲不变的情况下，分析了二级推力和三级推力方案时的方向密集度，从而得出了二级推力方案的方向密集度优于三级的结论，同时，提出了使方向密集度最小的最佳转换时间。     

翼型外形高气动效率/低可探测性的优化

 采用ｖａｎＬｅｅｒ矢通量分裂格式求解Ｅｕｌｅｒ方程的方法计算了绕翼型的气动特性;采用矢通量分裂方法计算了绕翼型的时域电磁散射场特性及雷达散射截面积（ＲＣＳ）;采用一种简单而有效的数值优化方法对流场解和电磁场解进行了翼型外形高气动效率／低可探测性的优化计算。算例结果表明,本方法提供了一种对翼型既可作气动优化设计亦可进行多学科综合设计的有效工具。     

等速上仰翼型绕流结构的数值模拟

通过数值求解涡量-流函数形式的N-S方程研究等速上仰翼型的分离流动结构。用ADI方法解涡量方程;用Poisson方程直接法解流函数方程。在计算得到的流场中除可看到实验中观察到的多涡结构外,还发现一些实验中未观察到的现象。将分离流动结构的计算结果与翼型气动力变化相联系,并结合涡动力学理论,定性地解释了上仰翼型产生高升力的原因。     

关于平面激波在凹两段楔面上的传播

众所周知,当一平面激波同任意形状的障碍物相互作用时,将产生十分复杂的波系,其中激波在曲面壁上的绕射和反射,将产生非定常波系。关于平面激波在凹圆柱壁面上的传播,我们可以用一连串、长度为无限小的楔面所构成的多边形来近似和代替凹圆柱面,因此分析和研究激波在两段楔面上反射和绕射的目的在于得出激波在圆柱壁面以及其它较复杂形状壁面上的传播规律。本文采用马赫-曾德干涉仪在激波管上进行光测模型实验。根据所拍摄的流场干涉图像,以分析其波系的演变过程及传播规律,得到有关流场的密度分布和压力分布,更为重要的是给出流场参数的定量结果。