有限平板绕流Monte—Carlo方法仿真

本文采用直接仿真Monte-Carlo方法求解有限长平板绕流等问题。本方法是通过计算机跟踪仿真分子的运动来实现数值模拟的。仿真分子间的碰撞计算由统计抽样确定。碰撞模型分别选用了硬球分子模型和负幂律分子模型。仿真分子与固壁作用采用由完全扩散反射和镜面反射按比例混合组成的模型。为了检验方法的可靠性,还计算了激波结构和Rayleigh问题等一维流动,二维计算采用同步并行程序。数值结果表明直接仿真Monte-Carlo方法能够较好地模拟稀薄气体力学中的一些问题。对于二维计算,所花费的机时和所需的内存均在国内机器所允许的范围内。     

改进跨音速位势方法的讨论

 简要地回顾并讨论了传统位势方法之不足。对改进传统位势方法的两种必要的修正(熵修正和涡修正)的量级在小扰动条件下作了分析和估值。结果表明涡修正是更高阶小量,因而作为第一步修正可采用非等熵位势流来近似。采用熵修正激波算子来计及非等熵的激波突跃条件,并考虑了由熵修正带来的尾迹条件的变化。算例表明非等熵位势方法能得到预期的改进位势方法的效果。     

高超声速无粘流场的推进计算

本文将文[1]中用于时间相关法计算的NND格式推广到定常超声速流动的空间推进计算,采用二步的预测、校正方法保证了推进方向的二阶精度,可以证明,这种二阶精度的NND格式具有TVD性质,是MacCormack二步显式格式的推广。本文首先将格式应用于二维平板上斜激波反射流场的推进计算,以检验格式捕捉激波的能力,同时研究了不同的通量分裂方法对格式捕捉激波能力的影响,得到了相当满意的结果。在此基础上,计算了航天飞机简化外形的身部超声速流场,给出了M_∞=10,α=0°,和M_∞=5,α=5°两种状态的部分结果,计算结果清楚地描绘了由于气流在机翼附近受到强烈压缩而产生的内嵌激波与外激波相交的复杂流场结构,与文[7]相比,流场结构更为清晰。     

固体发动机喷管内两相流的有限元数值模拟

本文给出了固体发动机喷管内两相流的数值模似。从非守恒型两相流甚本方程出发,采用加权积分技术,导出了离散点的流动物理量的有限元数值计算公式。数值模拟分别与JPL喷管的两相流壁面静压分布测量结果以及与差分法计算结果作了比较,比较结果表明:二者是吻合的。本文还以中型喷管为例,计算了粒子的轨迹,粒子的速度滞后和燃气的温度滞后,计算结果表明;入口粒子冲击角对粒子轨迹有明显的影响。     

改进格网型中的跨音速分析和设计计算

 <正> 近年来采用不同的主控方程计算跨音速绕流的方法有了不断的发展和改进,越来越多的算例表明了计算流体力学的迅速发展,其中解全速势方程的跨音速计算方法也有了很大进展,并已为国外航空工业界广泛地应用于设计和分析计算中。其中有限体积方法提高了处理复杂外形的能力,若采用多重网格技术则可改进计算效率。但是为了发展一种普遍适用和可靠的格网系统以耦合有限体积算法仍有许多工作要做