论文部分内容阅读
矩形柱绕流是流体力学领域重要的流体问题，广泛存在于工程实际中。研究矩形柱绕流的涡旋形态和力学性能对设计矩形截面结构具有重要的指导意义。然而，虽然针对无限流域中的矩形柱绕流已有广泛研究，但是，针对靠近自由水面的矩形柱绕流的研究较少。特别的，水深比(d/l)和宽长比(w/l)对近水面矩形柱绕流的影响、涡旋形态分类和雷诺数(Re)的影响等问题尚未得到解决。
本文采用基于雷诺平均的纳维尔-斯托克斯(Reynolds-averagedNavier-Stokes，简称RANS)方程的方法对近水面矩形柱绕流进行了二维两相流数值模拟。其中，RANS湍流模型采用Mentorshearstresstransport(SST)k-ω两方程模型，两相流模型采用流体体积法(Volumeoffluid，简称VOF)。本文从以下4个方面对靠近自由水面的矩形柱绕流进行了研究：
(1)近水面矩形柱绕流的两相流数值模拟。为验证二维RANSMentorSSTk-ω湍流模型和VOF两相流模型对近水面矩形柱绕流数值模拟的适用性，本文对d/l为4.5时方柱周围的速度和压强分布、不同w/l的矩形柱力系数和不同w/l的矩形柱涡旋流动进行了研究，并与已有的无限流域中矩形柱绕流的实验和数值模拟结果进行了对比。根据数值模拟验证，本文的二维两相流数值模型在0.7≤w/l≤5.0范围内能够较准确地模拟近水面矩形柱绕流。
(2)d/l和w/l对近水面矩形柱绕流的影响。为研究d/l和w/l对近水面矩形柱绕流影响的关系，本文首先对w/l为1.0和3.0的矩形柱绕流随d/l的变化进行了研究。然后，研究了在不同d/l时w/l对矩形柱绕流的影响，以探讨在自由水面的约束作用下w/l对矩形柱绕流影响的变化。研究发现，d/l对近水面矩形柱绕流的影响与w/l有关，并且，在不同d/l时w/l对矩形柱绕流的影响不同。
(3)近水面矩形柱绕流的涡旋形态分类。为研究不同d/l和w/l时近水面矩形柱绕流的涡旋形态，本文考虑了范围广阔的d/l和w/l。d/l和w/l的范围分别是0.3-4.5和0.7-5.0。根据对涡旋流动的观察，本文区分了3大类(反对称、非对称和涡旋周期脱落停止)和6小类(分离、重附着、非对称分离、单侧重附着、非对称重附着和涡旋周期脱落停止)涡旋形态。为确定各涡旋形态的划分边界，本文对速度和力系数随d/l和w/l的变化进行了分析，并根据速度和力系数的变化特点指定了划分边界的规则。涡旋形态划分图表明，当d/l大于1.5时，自由水面与矩形柱绕流之间的相互作用很小。涡旋周期脱落停止的临界d/l在w/l处于0.7≤w/l＜1.5、1.5≤w/l＜2.8和2.8＜w/l≤5.0范围内时分别位于0.5-0.6、0.7-1.0和0.8-0.9。
(4)Re对近水面方柱绕流的影响。为研究Re对近水面方柱绕流的影响，本文对2.0×104≤Re≤9.9×104范围内4个d/l时(0.3、0.7、1.0和4.5)的方柱绕流进行了数值模拟。研究发现，当d/l为4.5时，方柱绕流随Re变化很小，而当d/l为0.3、0.7和1.0时，方柱的涡旋形态和力系数随Re变化显著。另外，针对近水面方柱绕流中的2个重要的流动现象，即水面喷气式流动和涡旋周期脱落停止，本文分别从局部傅汝德数(Fr)和方柱绕流上下方剪切层的相互关系两个角度进行了深入分析。
本文采用基于雷诺平均的纳维尔-斯托克斯(Reynolds-averagedNavier-Stokes，简称RANS)方程的方法对近水面矩形柱绕流进行了二维两相流数值模拟。其中，RANS湍流模型采用Mentorshearstresstransport(SST)k-ω两方程模型，两相流模型采用流体体积法(Volumeoffluid，简称VOF)。本文从以下4个方面对靠近自由水面的矩形柱绕流进行了研究：
(1)近水面矩形柱绕流的两相流数值模拟。为验证二维RANSMentorSSTk-ω湍流模型和VOF两相流模型对近水面矩形柱绕流数值模拟的适用性，本文对d/l为4.5时方柱周围的速度和压强分布、不同w/l的矩形柱力系数和不同w/l的矩形柱涡旋流动进行了研究，并与已有的无限流域中矩形柱绕流的实验和数值模拟结果进行了对比。根据数值模拟验证，本文的二维两相流数值模型在0.7≤w/l≤5.0范围内能够较准确地模拟近水面矩形柱绕流。
(2)d/l和w/l对近水面矩形柱绕流的影响。为研究d/l和w/l对近水面矩形柱绕流影响的关系，本文首先对w/l为1.0和3.0的矩形柱绕流随d/l的变化进行了研究。然后，研究了在不同d/l时w/l对矩形柱绕流的影响，以探讨在自由水面的约束作用下w/l对矩形柱绕流影响的变化。研究发现，d/l对近水面矩形柱绕流的影响与w/l有关，并且，在不同d/l时w/l对矩形柱绕流的影响不同。
(3)近水面矩形柱绕流的涡旋形态分类。为研究不同d/l和w/l时近水面矩形柱绕流的涡旋形态，本文考虑了范围广阔的d/l和w/l。d/l和w/l的范围分别是0.3-4.5和0.7-5.0。根据对涡旋流动的观察，本文区分了3大类(反对称、非对称和涡旋周期脱落停止)和6小类(分离、重附着、非对称分离、单侧重附着、非对称重附着和涡旋周期脱落停止)涡旋形态。为确定各涡旋形态的划分边界，本文对速度和力系数随d/l和w/l的变化进行了分析，并根据速度和力系数的变化特点指定了划分边界的规则。涡旋形态划分图表明，当d/l大于1.5时，自由水面与矩形柱绕流之间的相互作用很小。涡旋周期脱落停止的临界d/l在w/l处于0.7≤w/l＜1.5、1.5≤w/l＜2.8和2.8＜w/l≤5.0范围内时分别位于0.5-0.6、0.7-1.0和0.8-0.9。
(4)Re对近水面方柱绕流的影响。为研究Re对近水面方柱绕流的影响，本文对2.0×104≤Re≤9.9×104范围内4个d/l时(0.3、0.7、1.0和4.5)的方柱绕流进行了数值模拟。研究发现，当d/l为4.5时，方柱绕流随Re变化很小，而当d/l为0.3、0.7和1.0时，方柱的涡旋形态和力系数随Re变化显著。另外，针对近水面方柱绕流中的2个重要的流动现象，即水面喷气式流动和涡旋周期脱落停止，本文分别从局部傅汝德数(Fr)和方柱绕流上下方剪切层的相互关系两个角度进行了深入分析。