论文部分内容阅读
随着计算流体力学(CFD,Computational Fluid Dynamics)计算方法和计算机技术的快速发展,人们逐渐倾向于以具有快速、高效、经济等优点的数值计算为研究手段,对型号部件或全机(弹)进行气动性能评估,然后进行实验,再进行精细的气动设计。这种研究思路被国内外普遍采用,是流体力学的发展结果,因此,计算流体动力学在飞行器气动设计中会发挥着越来越大的作用。但是,目前在热流、摩阻等方面CFD还有不足,离工程实用有较大的差距。本文旨在提高热流计算精度方面进行探索,开展了有限元以及有限元/有限体积混合方法方面的研究,构造了非正交单元有限元计算格式,提出了有限元/有限体积混合算法,并成功地应用于复杂外形流场计算中,计算结果表明本文建立的计算方法具有较高的精度和分辨率,热流计算结果令人满意。本文共分七章,概述如下:第一章为引言,简单回顾了计算流体力学中的有限元和有限体积计算方法的发展过程,概述了有限元、有限体积方法以及网格生成技术的发展现状和应用情况,介绍了本文的主要工作。第二章为数值方法介绍。在RKDG有限元方法的基础上,通过局部坐标变换,构造了适合粘性计算的二维和三维非正交单元DG有限元计算格式,提出了非正交三棱柱单元限制器方法,发展了适应于间断Galerkin有限元方法的隐式计算方法。鉴于有限元方法内存需求大和计算效率低的不足,在分析有限元和有限体积方法的各自优缺点的基础上,综合二者的优势,提出了有限元/有限体积区域混合计算方法。在该章中,同时给出了本文采用的控制方程和数值方法。第三章进行了RKDG有限元方法验证工作。采用RKDG有限元方法数值模拟了Sod问题、Shu问题、Lax问题、激波碰撞问题、来流马赫数为3的前台阶流动、强运动斜激波的双Mach反射问题,计算结果表明该方法在捕捉激波和接触间断方面具有较高的能力,显示了对流场结构的高分辨率,而且流场中没有明显的非物理波动。第四章进行了非正交单元间断Galerkin有限元计算格式验证工作。利用该方法数值模拟了二维圆柱粘性绕流、第四类激波干扰问题以及三维球头高超声速粘性流动问题。数值模拟结果显示了本文建立的非正交单元间断Galerkin有限元计算格式均具有较高的计算精度和流场分辨率,热流计算精度令人满意。第五章开展了有限元/有限体积混合算法的验证工作。利用该方法数值模拟了高超声速粘性流动问题,算例包括球头和20°攻角钝锥的粘性绕流。计算结果表明混合算法能够得到与完全有限元方法同样清晰的流场结构和较好的热流分布。同时,通过对两种方法占用计算机资源和计算效率的比较,表明了有限元/有限体积混合算法能大幅度减少内存需求和提高计算效率,克服了有限元方法的缺点。第六章进行了有限元/有限体积混合算法的应用研究。在第五章研究工作的基础上,将该混合算法应用于复杂外形的流场计算中,给出了与风洞实验一致的结果,显示了工程应用的潜力。第七章是本文的结束语。