论文部分内容阅读
复杂边界下和多孔介质内的流动和传热问题广泛存在于自然界和工程应用中,相关的研究在机械工程、能源工程、化学工程等领域处于十分重要的地位。由于这两类问题涉及复杂的热-流-固相互作用,理论研究和实验方法难以有效地刻画其流动规律和传热特性,而数值模拟以其特有的优势逐渐成为一种被广泛接受的研究方法。近年来流行的格子Boltzmann方法(lattice Boltzmann method,LBM),由于其粒子方法、介观尺度背景和计算方面的优势,非常适合处理复杂边界下和多孔介质内的流动和传热问题。尽管LBM已在这两类问题的研究中取得了一些成果,但仍存在一些亟待解决的问题。本文的主要工作和创新点如下:首先,在模拟复杂边界下的流动和传热问题的算法方面:(1)基于体积格子Boltzmann模型,提出了一种求解流固耦合问题的修正的动量交换法。该方法可较为自由地设置边界点,且计算过程中满足几何守恒律,因而有效地避免了流固相互作用力的非物理振荡。此外,提出了一种求解颗粒运动方程的时间离散格式,可用于流固密度比接近于1或小于1的情形;(2)提出了一种处理带复杂Dirichlet边界条件的流动和传热问题的迭代IB-LBM。通过选择合适的迭代次数和松弛参数,速度和温度Dirichlet边界条件能够得到更好的满足;(3)基于界面跳条件和非光滑区域上的插值公式,提出了一种处理带复杂Neumann边界条件的流动和传热问题的IB-LBM,并采用该方法发现了同心环形腔内的自然对流的一个新的稳态解;(4)推导了流体/固体区域的带源项的统一温度场控制方程,界面温度共轭条件可自动得到满足。基于该统一控制方程,构造了一种求解流固共轭传热问题的纯欧拉格子Boltzmann模型。其次,在模拟多孔介质内的流动和传热问题的算法方面:(1)构造了 REV尺度各向异性多孔介质流动和传热问题的多松弛时间格子Boltzmann模型,并采用Chapman-Enskog方法证明了该模型对应的宏观方程为各向异性Darcy-Brinkman-Forchheimer方程和能量方程。此外,使用该模型模拟了壁面附着各向异性多孔介质层的自然对流问题,并分析了渗透率之比和热传导系数之比对流场结构以及传热性能的影响;(2)基于格子Boltzmann通量格式,提出了一种高效的求解Darcy-Brikmann-Forchhermer方程的有限体积法(FVM)。该方法采用重构基于通用多孔介质流模型的LB方程的局部解的方式来获得单元界面上的数值通量。流场控制方程的时间离散采用中点积分公式,放松了时间步长的限制。