论文部分内容阅读
钝体绕流广泛存在于航天、建筑、交通、环境等许多研究领域,因流体在钝体表面产生大范围的分离并形成宽阔的伴有旋涡脱落现象的尾流,使得具有简单几何边界的方形钝体周围的湍流场的模拟变得非常复杂。不仅为考察数学模型的计算能力,更为了保证建筑结构安全、减弱热岛现象、强化换热、促进污染物的扩散,对这种流动现象进行深入研究是十分必要和有意义的。本文通过分析国内外研究现状及存在问题基础上,结合理论分析、数值研究和实验方法主要完成了以下工作:首先研究了边界条件对钝体湍流场数值模拟计算结果的影响,通过对不同入口流速与湍流强度的多种工况的数值计算,考察了入流速度变化对流场中方形钝体顶部回流区分离长度及是否会产生附着的影响,分析了影响机理,确定了不同条件下速度函数中指数m的适用值,从而提高模型对绕流湍流场的预测的能力。分析了湍流强度对方形钝体顶面压力分布的影响,当湍流强度逐渐增大时,最大湍流动能kmax会向前移动并开始靠近钝体前缘拐点。通过对标准k-ε模型在研究方形钝体绕流运动中产生的问题,特别是钝体尖角处模拟计算值过高估算的问题深入分析发现,不能准确计算钝体顶部回流涡及湍流动能最大值kmax位置的原因,并不是完全由标准k-ε模型本身造成,更要考虑如何处理钝体边界处的壁面处理方法,通过对尖角处w和v相交过程的详细分析,明确了问题的原因,提出通过限制流体在钝体尖角转弯处来流方向的速度的方法,降低数值离散中控制体内数据传递的误差,消除横向速度在尖角处出现比实际值大的现象,从而达到准确预测钝体湍流场特征参数的目的。采用限定边界条件的方法,使用标准k-ε模型和LES模型对风洞内贴壁方形钝体绕流进行模拟计算,与经典实验结果进行比较,验证了标准k-ε模型的精度和适用性。同时分析了网格精度和计算域敏感度对LES模型模拟的影响,提出网格加密不会自动导致更好的结果,标准网格模拟的结果与实验的吻合度比细网格的更好。本文利用粒子成像PIV测速技术测量了通道内矩形截面贴壁方柱钝体模型的速度场,与模拟结果进行对比分析,除尾部回流区高于实验结果,其他位置结果较吻合,验证了限定迎风尖角壁面条件的基础上的标准k-ε模型的预测能力。选用德克萨斯州技术大学校园内的一栋测试建筑物(TTU)和玻璃屋(GLH)两种全尺寸建筑模型,考察了在大气边界层内标准k-ε模型的适用性,通过模拟风场的数值计算,在确定合适的入口边界后获得与实验数据较吻合的压力场结果,进一步验证了RANS法可用于大型工程流场研究,并研究明确了网格尺寸和分布、计算域尺寸、计算维数等的影响。并研究了LES法适用的入流边界条件。从计算精度和计算成本分析确定了RANS的适用性,并对复杂实体建筑物模型进行了模拟计算,考察了三种流速和两种流动方向的流场情况,分析了风速对建筑周围分离涡的影响及屋面压力的变化。本文的研究结果为加深流体问题物理现象的理解提供了理论依据,同时对钝体绕流场的研究可为相关工程领域设计提供参考,为以后更深入的研究打下基础。