瞬态热传导问题的光滑粒子流体动力学研究

摘要

温度场的计算是流体机械的一个重要研究领域,在瞬态温度场函数中,温度的求解包括时间域和空间域,因此如何精确、有效地求解温度场的分布情况,在实际应用领域上显得尤为重要。
　　 本文研究了求解瞬态热传导问题的一种新型无网格数值计算方法——光滑粒子流体动力学(Smooth Particle Hydrodynamics,SPH)方法,它是一种无网格自适应拉格朗日粒子法,因其拥有不受网格限制、自适应、拉格朗日等优点,已成为目前科学和工程计算方法的研究热点之一。
　　 本文用两种方法将SPH方法应用于瞬态热传导问题的数值模拟:一是将瞬态热传导问题的二阶导的偏微分方程降阶为两个一阶导的偏微分方程,再依次应用SPH方法；二是将SPH方法直接应用到二阶导的瞬态热传导问题中。对于SPH方法的边界处理,本文采用了镜像粒子法和固定边界粒子法,前者通过镜像原理在边界外生成与内部粒子相似的粒子,后者是通过内部粒子加权求和计算出边界粒子的参数。
　　 本文从SPH方法开始入手,分析如何将其应用于瞬态热传导问题的数值模拟,然后介绍用FORTRAN语言编写SPH方法的计算程序,最后通过大量的数值算例,验证本文方法的可行性和有效性。本文对算例结果、计算精度及SPH方法的边界处理等问题进行对比分析得出:对规则边界和不规则边界的几何模型,两种SPH方法都可进行求解,同时两种方法均能通过粒子数的增加来提高计算精度；二阶导SPH方法在满足()2W(r,h)≥0的情况下,与降阶求解SPH方法对同一算例选用Lucy核函数进行求解对比,发现二阶导SPH方法的计算精度和光滑性略优于降阶求解SPH方法；镜像粒子法对边界粒子的处理要优于固定边界粒子法。