基于大涡模拟的并列双方柱扰流特性及影响因素研究

基于大涡模拟的并列双方柱扰流特性及影响因素研究一、引言1.1研究背景与意义在众多工程领域中，并列双方柱绕流现象广泛存在，对其进行深入研究具有至关重要的意义。以建筑工程为例，城市中的高层建筑群，当强风来袭时，风绕过多栋并列的方形建筑，形成复杂的扰流。这种扰流不仅会对建筑结构产生额外的荷载，影响建筑的稳定性和安全性，还可能在建筑周围形成局部强风区域，对行人的安全和舒适性造成威胁。若建筑设计初期没有充分考虑这种绕流效应，随着时间推移，建筑结构可能因长期承受不合理的荷载而出现损坏，如墙体裂缝、结构变形等，严重时甚至可能导致建筑坍塌，造成巨大的生命财产损失。在海洋工程领域，海上的风力发电装置、石油钻井平台等，通常包含多个并列的柱状结构。海水绕流这些结构时，会产生复杂的流动现象，如漩涡脱落、尾流干扰等。这些现象会引起结构的振动，长期的振动可能导致结构疲劳，降低结构的使用寿命。若振动幅度过大，还可能引发结构的共振，直接威胁到海洋工程设施的安全运行，导致石油泄漏、电力中断等严重事故，对海洋生态环境和能源供应产生巨大的负面影响。大涡模拟（LES）方法作为一种重要的湍流数值模拟方法，在研究并列双方柱扰流问题上展现出独特的优势。相较于传统的雷诺平均纳维-斯托克斯（RANS）方法，大涡模拟能够直接求解大尺度涡旋的运动，通过亚网格模型对小尺度涡旋进行模拟，从而更准确地捕捉流场中的复杂流动结构和非定常特性。在研究并列双方柱绕流时，大涡模拟可以清晰地呈现出不同间距比下双方柱之间的相互干扰机制，如涡街的形成、发展和相互作用过程。同时，大涡模拟对复杂几何形状和边界条件具有良好的适应性，这使得它在处理实际工程中的并列双方柱绕流问题时，能够更加贴近真实情况。随着计算机技术的飞速发展，计算资源的不断提升，大涡模拟方法在工程领域的应用潜力日益凸显，有望为建筑、海洋等工程领域的设计和优化提供更可靠的理论依据和技术支持。1.2国内外研究现状在并列双方柱扰流的研究领域，国内外学者通过实验与数值模拟两种途径展开了深入探究。在实验研究方面，国外起步较早，积累了丰富的成果。早在20世纪中叶，一些学者就开始通过风洞实验研究钝体绕流问题，其中包括对并列双方柱绕流的初步探索。他们利用热线风速仪等设备测量流场中的速度分布，观察尾流中的涡旋结构。随着技术的不断进步，粒子图像测速（PIV）技术、激光多普勒测速（LDV）技术等先进测量手段逐渐应用于实验研究中。例如，[具体文献]通过PIV技术，对不同间距比下并列双方柱绕流的流场进行了可视化测量，清晰地展示了流场中的涡旋分布和演化过程，发现间距比的变化会导致涡旋的产生、发展和相互作用方式发生显著改变。国内在该领域的实验研究也取得了长足的进展。近年来，众多高校和科研机构积极开展相关实验，一些团队利用风洞实验研究了不同雷诺数和间距比下并列双方柱的气动特性。通过测量柱体表面的压力分布，得到了平均压力系数、升力系数和阻力系数等重要参数，分析了这些参数随雷诺数和间距比的变化规律。如[具体文献]在风洞实验中，针对不同间距比的并列双方柱，详细测量了其表面压力分布，研究结果表明，在小间距下，并列双方柱的平均压力系数不仅在内侧面有显著影响，在其他面也有轻微影响，并且干扰效应呈现不对称性。数值模拟方面，雷诺平均纳维-斯托克斯（RANS）方法曾是早期研究并列双方柱绕流的常用手段。通过对控制方程进行时均化处理，结合各种湍流模型（如标准k-ε模型、RNGk-ε模型等）来封闭方程，从而求解流场。这种方法在一定程度上能够预测流场的平均特性，例如阻力系数、升力系数的平均值等。然而，RANS方法对湍流脉动信息的描述能力有限，难以准确捕捉流场中的复杂非定常现象和精细的涡旋结构。随着计算机技术的迅猛发展，大涡模拟（LES）方法逐渐在并列双方柱扰流研究中崭露头角。国外学者率先将大涡模拟应用于这一领域，[具体文献]利用大涡模拟对高雷诺数下并列双方柱绕流进行了数值研究，通过对大尺度涡旋的直接求解和小尺度涡旋的亚网格模型模拟，成功捕捉到了流场中丰富的非定常信息和复杂的涡旋结构，揭示了不同间距比下双方柱之间的涡旋相互作用机制，如涡街的形成、发展和相互干扰过程。国内也有诸多学者运用大涡模拟方法深入研究并列双方柱绕流问题。[具体文献]采用大涡模拟结合有限体积法，对不同间距比下的并列双方柱绕流进行了数值模拟，详细分析了流场的速度分布、涡量分布以及柱体的受力情况，研究发现雷诺数和间距比都会对流场产生影响，但间距比对流场的影响更为显著，当间距小于临界间距时，会出现双稳态偏流现象，流场在一段时间内偏向一个柱体，随后又会改变偏流方向偏向另一个柱体；当间距大于临界间距时，流场则为非偏流状态。总体而言，目前国内外对于并列双方柱扰流的研究已取得了丰硕的成果，但仍存在一些有待进一步深入探讨的问题。例如，在大涡模拟中，亚网格模型的精度和适用性仍需不断改进和完善，以提高对小尺度涡旋的模拟准确性；同时，对于复杂工况下（如高雷诺数、多柱体干扰等）并列双方柱扰流的研究还相对较少，需要开展更多的研究工作来揭示其内在的流动规律和物理机制。1.3研究内容与方法本文旨在通过大涡模拟方法，深入研究并列双方柱扰流问题，具体研究内容如下：模型验证与对比：运用CFD软件FLUENT，采用标准k-ε模型和大涡模拟（LES）模型，对雷诺数为22000时的单方柱绕流进行数值模拟。将模拟得到的阻力系数、Strouhal数等结果，与已有的实验结果以及其他作者的计算结果进行详细对比。通过这一过程，一方面验证FLUENT软件中这两种模型在不可压缩绕流流场计算中的可行性，另一方面凸显LES模型相较于标准k-ε模型在捕捉复杂流动细节方面的优越性，为后续并列双方柱绕流的模拟研究奠定坚实基础。低雷诺数下并列双方柱绕流模拟：在低雷诺数条件下，利用大涡模拟方法对并列双方柱绕流进行数值模拟。重点研究不同间距比时流场的特性，包括但不限于流场中速度分布、涡量分布等。通过分析这些特性，深入探讨间距比变化对流场的影响规律，揭示低雷诺数下并列双方柱之间的干扰机制。高雷诺数下并列双方柱绕流模拟：针对高雷诺数工况，同样运用大涡模拟对并列双方柱绕流展开研究。详细模拟不同间距比下的流场情况，获取流场的各种参数和特征。与低雷诺数情况进行对比，分析雷诺数和间距比这两个因素对流场影响的差异，探究高雷诺数下并列双方柱绕流的独特规律，以及间距比在其中所起的关键作用。流场特性分析：对低雷诺数和高雷诺数下的模拟结果进行全面深入的分析。除了关注流场的速度、涡量分布外，还将对柱体的受力情况进行详细研究，包括阻力系数和升力系数的变化规律。同时，利用频谱分析等方法，研究流场的非定常特性，如涡旋的脱落频率、流场的脉动特性等，从多个角度揭示并列双方柱绕流的复杂物理机制。在研究方法上，选用CFD软件FLUENT作为数值模拟工具。该软件具有强大的功能和广泛的适用性，能够处理各种复杂的流动问题。在数值计算过程中，采用有限体积法对控制方程进行离散。有限体积法的基本思想是将计算区域划分为一系列控制体积，使每个控制体积都包围一个网格节点。通过对每个控制体积应用守恒定律，将偏微分形式的控制方程转化为离散的代数方程。这种方法在保证物理量守恒的同时，具有良好的计算稳定性和精度，能够有效地求解不可压缩流体的Navier-Stokes方程。对于压力-速度耦合问题，采用SIMPLE（Semi-ImplicitMethodforPressure-LinkedEquations）算法进行求解。SIMPLE算法通过引入压力修正方程，实现压力和速度的迭代求解，使得计算过程更加稳定和高效。在大涡模拟中，选用动力亚网格尺度（SGS）模型来模拟小尺度涡旋的影响。动力SGS模型能够根据流场的局部特性动态调整亚网格尺度的参数，相比传统的Smagorinsky模型，具有更高的精度和适应性，能够更准确地模拟小尺度涡旋对大尺度流动的影响。二、大涡模拟基本原理与方法2.1大涡模拟基本概念大涡模拟（LargeEddySimulation，LES）是近几十年发展起来的一种在流体力学领域极为重要的数值模拟研究方法。其核心思想是基于湍流运动由不同尺度涡旋组成这一特性，大尺度涡旋对平均流动影响显著，各种变量的湍流扩散、热量、质量、动量和能量的交换以及雷诺应力的产生主要通过大尺度涡旋实现；而小尺度涡旋主要作用于耗散，通过耗散脉动影响各种变量。并且不同的流场形状和边界条件对大涡旋影响较大，使其具有明显的各向不均匀性，小涡旋则近似于各向同性，受边界条件影响小，具有较大共性，便于建立通用模型。基于此，大涡模拟将湍流中的大涡旋（大尺度量）和小涡旋（小尺度量）分开处理，对大涡旋通过N-S方程直接求解，对小涡旋则通过亚格子尺度模型，建立其与大涡旋的关系来进行模拟，大小涡旋通过滤波函数加以区分。大涡模拟的发展历程见证了其在湍流研究领域的逐步崛起。早期，由于计算机技术的限制，对湍流的数值模拟面临诸多挑战。随着计算机性能的不断提升，大涡模拟方法逐渐崭露头角。从最初的理论探索阶段，研究人员致力于完善其基本原理和方法，到后来不断应用于各种简单流场的模拟，验证其可行性和有效性。如今，大涡模拟已广泛应用于航空航天、海洋工程、建筑工程等众多领域，成为研究湍流现象不可或缺的重要手段。在湍流模拟领域，大涡模拟占据着独特而重要的地位，与直接数值模拟（DNS）和雷诺平均（RANS）方法共同构成了湍流数值模拟的主要方法体系。直接数值模拟通过直接求解完整的三维非定常N-S方程组，能够获得湍流场的精确信息，为研究湍流机理提供了有效手段。然而，其对计算机内存容量和计算时间的要求极高，目前仅能局限于较小雷诺数的湍流模拟。以模拟雷诺数仅为3300的槽流为例，所用的网点数N就约达到了2×10⁶，且在向量计算机上运行了250小时，这使得直接数值模拟在实际工程应用中受到极大限制。雷诺平均方法则是将非稳态的N-S方程对时间作平均，求解工程中需要的时均量。这种方法计算速度相对较快，能够满足一定的设计精度要求，成为目前工程应用中较为常用的方法。但它给出的只是平均运动结果，无法反映流场紊动的细节信息，对于一些需要考虑非定常特性和复杂涡旋结构的流动问题，其模拟能力存在局限性。相比之下，大涡模拟兼具两者的部分优势。它基于湍动能传输机制，直接计算大尺度涡的运动，小尺度涡运动对大尺度涡的影响通过建立模型体现出来。这使得大涡模拟既能够得到比雷诺平均方法更多的诸如大尺度涡结构和性质等的动态信息，又比直接数值模拟节省计算量，在计算成本和模拟精度之间找到了较好的平衡，因而得到了越来越广泛的发展和应用。在研究并列双方柱绕流问题时，大涡模拟能够捕捉到流场中复杂的非定常现象和精细的涡旋结构，为深入理解绕流机理提供了有力支持。2.2控制方程与滤波方法大涡模拟所基于的控制方程是经过滤波处理后的Navier-Stokes方程。Navier-Stokes方程是描述粘性不可压缩流体动量守恒的运动方程，其在直角坐标系下的张量形式为：\frac{\partialu_{i}}{\partialt}+u_{j}\frac{\partialu_{i}}{\partialx_{j}}=-\frac{1}{\rho}\frac{\partialp}{\partialx_{i}}+\nu\frac{\partial^{2}u_{i}}{\partialx_{j}\partialx_{j}}+f_{i}\frac{\partialu_{i}}{\partialx_{i}}=0其中，u_{i}和u_{j}分别为i方向和j方向的速度分量，t为时间，\rho为流体密度，p为压力，\nu为运动粘性系数，f_{i}为i方向的体积力。在大涡模拟中，为了将大尺度涡和小尺度涡分离，需要对Navier-Stokes方程进行滤波处理。滤波过程是大涡模拟的关键步骤，它通过一个滤波函数将流场变量分解为大尺度分量和小尺度分量。常用的滤波方法有多种，其中高斯滤波和盒式滤波是较为典型的两种。高斯滤波函数的数学表达式为：G(x,x')=\frac{1}{(2\pi\Delta^{2})^{\frac{3}{2}}}e^{-\frac{|x-x'|^{2}}{2\Delta^{2}}}式中，x和x'分别为空间坐标，\Delta为滤波宽度，它决定了能够被滤除的小尺度涡的大小。高斯滤波函数具有良好的平滑特性，能够在滤除小尺度涡的同时，较好地保留大尺度涡的信息。在实际应用中，对于一些对平滑度要求较高的流场模拟，如航空航天领域中飞行器周围流场的模拟，高斯滤波可以有效地去除高频噪声，使模拟结果更加稳定和准确。盒式滤波函数则定义为：G(x,x')=\begin{cases}\frac{1}{\Delta^{3}},&\text{if}|x-x'|\leq\frac{\Delta}{2}\\0,&\text{otherwise}\end{cases}盒式滤波函数相对简单直观，其滤波宽度\Delta明确，在计算上较为便捷。在一些对计算效率要求较高，且对小尺度涡的滤除精度要求相对不那么苛刻的工程应用中，如建筑风工程中对建筑物周围风场的初步模拟，盒式滤波能够快速地将小尺度涡从流场中分离出来，为后续的分析提供基础。通过滤波操作，将速度u_{i}分解为大尺度速度\overline{u}_{i}和小尺度速度u_{i}^{'}，即u_{i}=\overline{u}_{i}+u_{i}^{'}。对Navier-Stokes方程进行滤波后，得到大涡模拟的控制方程：\frac{\partial\overline{u}_{i}}{\partialt}+\overline{u}_{j}\frac{\partial\overline{u}_{i}}{\partialx_{j}}=-\frac{1}{\rho}\frac{\partial\overline{p}}{\partialx_{i}}+\nu\frac{\partial^{2}\overline{u}_{i}}{\partialx_{j}\partialx_{j}}-\frac{\partial\tau_{ij}}{\partialx_{j}}+\overline{f}_{i}\frac{\partial\overline{u}_{i}}{\partialx_{i}}=0其中，\tau_{ij}=\overline{u_{i}u_{j}}-\overline{u}_{i}\overline{u}_{j}为亚网格应力张量，它体现了小尺度涡对大尺度涡的影响。在大涡模拟中，需要通过亚网格模型对亚网格应力张量进行模拟，从而实现对小尺度涡的间接求解。2.3子网格尺度模型在大涡模拟中，子网格尺度（SGS）模型起着至关重要的作用，它用于模拟小尺度涡旋对大尺度运动的影响，以弥补由于滤波操作而被滤掉的小尺度信息。目前，存在多种子网格尺度模型，它们各自具有独特的特点和适用范围。Smagorinsky模型是最早提出的子网格尺度模型之一，具有形式简单、计算效率高的优点。该模型基于涡粘性假设，将亚网格应力与大尺度应变率联系起来，通过一个经验常数来确定涡粘性系数。其涡粘性系数的表达式为：\nu_{t}=(C_{s}\Delta)^{2}\left|\overline{S}_{ij}\right|其中，C_{s}为Smagorinsky常数，\Delta为滤波宽度，通常取为网格尺度，\overline{S}_{ij}为大尺度应变率张量。在一些对精度要求不是特别高，且计算资源有限的情况下，如对简单几何形状的流场进行初步模拟时，Smagorinsky模型能够快速给出较为合理的结果。然而，该模型的常数C_{s}是固定的，无法根据流场的局部特性进行动态调整，这使得它在模拟复杂流场时存在一定的局限性。例如，在模拟高雷诺数下的湍流边界层时，固定的C_{s}值可能导致对壁面附近的湍流特性模拟不准确，因为壁面附近的流场特性与远离壁面的区域有很大差异。动力亚网格尺度（SGS）模型则克服了Smagorinsky模型中常数固定的缺点，能够根据流场的局部信息动态计算模型系数。它通过引入测试滤波，在不同尺度上对大涡模拟方程进行分析，从而确定模型系数。具体来说，动力SGS模型利用大尺度和测试尺度之间的信息传递关系，动态调整涡粘性系数。这种动态调整使得模型能够更好地适应不同的流场条件，在模拟复杂流场时具有更高的精度。例如，在模拟并列双方柱绕流时，流场中存在复杂的涡旋相互作用和非定常特性，动力SGS模型能够根据流场的实时变化，准确地模拟小尺度涡旋对大尺度运动的影响，从而更真实地反映流场的物理过程。WALE（Wall-AdaptiveLocalEddy-viscosity）模型是一种针对壁面附近流动的子网格尺度模型。它通过引入一个与应变率张量和旋转率张量相关的不变量，来定义涡粘性系数。该模型在壁面附近具有良好的适应性，能够准确模拟壁面附近的湍流特性。这是因为WALE模型考虑了壁面的影响，对壁面附近的小尺度涡旋进行了更合理的模拟。在研究并列双方柱绕流问题时，柱体壁面附近的流场特性对整个流场的影响很大，WALE模型能够很好地捕捉这部分信息。然而，在远离壁面的区域，WALE模型的表现可能不如其他一些通用的子网格尺度模型，因为它主要是为壁面附近流动设计的，对远离壁面区域的流场适应性相对较弱。在本文的并列双方柱扰流大涡模拟研究中，选用动力亚网格尺度（SGS）模型。这主要是因为并列双方柱绕流流场十分复杂，存在着强烈的涡旋相互作用、非定常特性以及复杂的边界条件。动力SGS模型能够根据流场的局部特性动态调整亚网格尺度的参数，从而更准确地模拟小尺度涡旋对大尺度流动的影响。相比之下，Smagorinsky模型由于其常数固定，难以适应流场的复杂变化；WALE模型虽然在壁面附近表现出色，但在整个流场的通用性上不如动力SGS模型。选用动力SGS模型能够充分发挥其优势，为准确揭示并列双方柱绕流的复杂物理机制提供有力支持。2.4边界条件处理在大涡模拟中，边界条件的设置对模拟结果的准确性和可靠性起着至关重要的作用。合理的边界条件能够确保流场在边界处的物理特性得到准确的描述，从而使模拟结果更接近实际情况。入口边界条件通常采用速度入口条件，即给定入口处的速度分布。在研究并列双方柱扰流时，根据实际工况，假设入口来流为均匀流，给定入口速度的大小和方向。这种设置方式能够模拟出流体在进入计算区域时的初始状态，为后续流场的发展提供基础。例如，在模拟风绕过高楼群时，根据当地的气象数据，确定入口风速和风向，以此作为入口边界条件。准确设置入口速度，能够真实反映风在进入建筑区域前的状态，对于研究建筑周围的风场分布至关重要。若入口速度设置不合理，如速度过大或过小，将导致模拟得到的流场与实际情况偏差较大，无法准确预测建筑所受的风荷载以及周围的风环境。出口边界条件一般选用压力出口条件，指定出口处的压力值。在并列双方柱绕流模拟中，通常将出口压力设为环境压力，以保证流场在出口处的压力平衡。这是因为在实际情况中，流体流出计算区域后，会与周围环境相互作用，最终达到环境压力状态。若出口压力设置不当，会导致流场在出口处出现回流或压力异常等现象，影响整个流场的稳定性和准确性。例如，在模拟河流中水流绕过桥墩的情况时，将出口压力设为河流下游的实际压力，能够使模拟结果更符合实际水流的流动情况，准确反映桥墩周围的水流特性和桥墩所受的水动力。壁面边界条件对于模拟柱体表面的流动特性十分关键。在大涡模拟中，采用无滑移边界条件，即认为流体在柱体壁面处的速度为零。这是基于实际物理现象，流体与固体壁面之间存在附着力，使得壁面处的流体速度与壁面速度相同，在静止柱体壁面处，流体速度即为零。无滑移边界条件能够准确模拟柱体壁面附近的边界层发展，捕捉到壁面附近的流动细节，如边界层的分离、再附等现象。这些信息对于研究柱体的受力情况和流场的能量损失具有重要意义。例如，在模拟海洋工程中的立管绕流时，准确的无滑移边界条件能够清晰地呈现出立管壁面附近的水流速度分布和涡旋结构，为分析立管的疲劳寿命和安全性提供重要依据。除了上述主要边界条件外，在一些复杂的模拟中，还可能涉及到其他特殊边界条件的设置。例如，在考虑热传递的流场模拟中，需要设置壁面的热边界条件，如给定壁面温度或热流密度。在研究具有自由表面的流动时，如波浪与柱体的相互作用，需要采用自由表面边界条件来准确描述自由表面的运动。这些特殊边界条件的合理设置，能够进一步拓展大涡模拟的应用范围，使其能够更全面地模拟各种复杂的流动现象。在研究海洋中波浪对海上平台立柱的作用时，自由表面边界条件能够真实地模拟波浪的起伏和传播，以及波浪与立柱之间的相互作用过程，为海上平台的设计和安全评估提供关键信息。三、并列双方柱扰流数值模拟方案3.1计算模型建立本文针对并列双方柱绕流问题构建了三维几何模型，该模型旨在模拟风或水流等流体绕过两个并列放置的方形柱体时的流动情况。模型中，单个方柱的边长D设定为0.1m，这一尺寸的选取主要基于实际工程应用中的常见尺寸范围，同时考虑到计算资源和模拟精度的平衡。例如，在建筑工程中，一些小型建筑构件或桥梁的桥墩等结构，其截面尺寸常常在这一量级范围内。通过选择具有代表性的尺寸，能够使模拟结果更具实际参考价值。双方柱之间的间距比L/D是本研究的关键参数之一，其中L表示双方柱中心之间的距离。本研究选取了L/D=1.5、2.0、2.5、3.0、4.0这几个不同的间距比进行模拟。选择这些间距比的依据在于，它们能够涵盖并列双方柱绕流中常见的间距范围，并且在以往的相关研究中也被广泛采用。较小的间距比如L/D=1.5，双方柱之间的相互干扰作用强烈，流场呈现出高度的复杂性；而较大的间距比如L/D=4.0，双方柱之间的干扰相对较弱，流场特性逐渐接近单方柱绕流。通过对这一系列不同间距比的模拟，可以全面深入地探究间距比变化对并列双方柱绕流流场的影响规律。在建立模型时，充分考虑了实际工程背景和以往研究成果。在实际工程中，如城市建筑群中的建筑布局、海洋工程中的海上平台群等，并列的柱状结构之间的间距各不相同，通过研究不同间距比下的绕流特性，能够为工程设计提供更全面的理论支持。以往的研究表明，间距比的变化会导致流场中涡旋的产生、发展和相互作用方式发生显著改变，进而影响柱体的受力情况和周围流场的分布。因此，对不同间距比进行系统研究具有重要的理论和实际意义。3.2计算域与网格划分在进行并列双方柱绕流的数值模拟时，合理确定计算域的范围和形状至关重要。计算域的设置直接影响到模拟结果的准确性和计算效率。经过综合考虑，本文将计算域设定为一个长方体，其尺寸的选取基于对流体流动特性和边界条件的分析。在入口方向（即流体流入的方向），计算域长度设置为10D，这一长度能够确保入口来流在进入计算域后有足够的发展空间，使流场在到达方柱之前能够充分稳定，避免入口边界条件对柱体附近流场的干扰。例如，在实际的风洞实验中，通常也会在模型前方设置一定长度的稳定段，以保证来流的均匀性和稳定性。在本文的模拟中，10D的入口段长度能够达到类似的效果，为后续流场的模拟提供稳定的初始条件。在出口方向，计算域长度取为20D。这样的长度可以保证流体在绕过方柱后，其尾流能够在计算域内充分发展和衰减，避免出口边界对柱体尾流的反射影响，从而更准确地模拟真实的流场情况。以河流中水流绕过桥墩的情况为例，桥墩下游的水流需要一段较长的距离才能恢复到正常状态，在数值模拟中，足够长的出口段就能够模拟这一过程，使模拟结果更符合实际水流的流动特性。在垂直于来流方向上，计算域的宽度和高度均设定为10D。这一尺寸能够有效避免边界对柱体周围流场的影响，确保流场在横向和纵向都有足够的空间进行发展。当流体绕过并列双方柱时，会在柱体周围形成复杂的三维流场，包括涡旋的产生、发展和脱落等现象。足够大的计算域宽度和高度能够容纳这些复杂的流动结构，使模拟结果能够真实反映流场的全貌。网格划分是数值模拟中的关键环节，它直接影响到计算结果的精度和计算效率。本文采用结构化网格对计算域进行划分，结构化网格具有网格质量高、计算精度好、数据存储量小等优点。在结构化网格中，网格节点按照一定的规律排列，形成整齐的网格结构，这使得在进行数值计算时，计算格式更加简单，计算精度更容易保证。为了提高模拟精度，在柱体附近和尾流区域进行了网格加密。柱体附近是流场变化最为剧烈的区域，存在着边界层分离、涡旋的产生等复杂现象，需要精细的网格来准确捕捉这些流动细节。在尾流区域，涡旋的发展和相互作用也需要较高的网格分辨率来进行模拟。通过在这些关键区域加密网格，可以更好地描述流场的物理特性，提高模拟结果的准确性。例如，在模拟风绕过高楼的情况时，高楼表面和周围的流场变化非常复杂，通过在高楼表面和尾流区域加密网格，可以清晰地呈现出风的流动路径、涡旋的形成和传播等现象，为建筑的抗风设计提供重要依据。为了确保模拟结果的准确性与可靠性，进行了网格无关性分析。分别采用了粗、中、细三种不同密度的网格进行模拟。粗网格的单元数量相对较少，计算成本较低，但可能无法准确捕捉流场的细微变化；细网格的单元数量较多，能够更精确地描述流场，但计算时间和计算资源的消耗也会相应增加；中等网格则介于两者之间。在粗网格划分中，整体网格较为稀疏，柱体附近和尾流区域的网格分辨率相对较低，单元总数约为[X1]个。中等网格在粗网格的基础上，适当增加了柱体附近和尾流区域的网格密度，单元总数约为[X2]个。细网格则进一步加密，尤其是在柱体表面和尾流核心区域，单元总数约为[X3]个。通过对不同网格密度下的模拟结果进行对比分析，重点关注阻力系数、升力系数以及流场中的涡量分布等关键参数。当网格从粗到细逐渐加密时，观察这些参数的变化情况。若随着网格加密，参数的变化趋于稳定，即相邻两种网格密度下的参数差异在可接受范围内（一般认为相对误差小于5%），则认为此时的模拟结果与网格密度无关，即达到了网格无关性。例如，在比较不同网格密度下的阻力系数时，发现粗网格与中等网格的阻力系数相对误差为[Y1]%，中等网格与细网格的阻力系数相对误差为[Y2]%，当[Y2]%小于5%时，表明从中等网格到细网格，阻力系数已基本稳定，达到了网格无关性。经过对比，发现中等网格既能保证计算精度，又能在合理的计算时间内完成模拟，因此最终选择中等网格进行后续的模拟研究。3.3模拟参数设置在本次并列双方柱扰流的大涡模拟研究中，雷诺数（Re）是一个关键的无量纲参数，它表征了流体惯性力与黏性力的比值。根据研究目的和实际工程应用场景，选取了两个具有代表性的雷诺数进行模拟，分别为低雷诺数Re=1000和高雷诺数Re=100000。低雷诺数Re=1000时，黏性力在流场中占据主导地位，流体的流动相对较为平稳，涡旋的生成和发展相对较弱。在这种情况下，流场中的流动特性主要受黏性力的影响，流体质点的运动轨迹相对规则，流场中的速度分布和压力分布相对较为均匀。通过模拟这一低雷诺数下的并列双方柱绕流，可以深入研究黏性力主导下双方柱之间的干扰机制以及流场的基本特性。例如，在微流控芯片中的流体流动，雷诺数通常较低，研究低雷诺数下的绕流特性对于芯片的设计和优化具有重要意义。高雷诺数Re=100000时，惯性力成为主导因素，流场中会产生强烈的湍流现象，涡旋的生成、发展和相互作用十分复杂。在高雷诺数下，流体的流动变得不稳定，速度和压力呈现出剧烈的脉动，流场中存在各种尺度的涡旋结构。模拟这一高雷诺数工况，能够揭示惯性力主导下并列双方柱绕流的复杂流动规律，如涡街的形成、发展和脱落过程，以及柱体所受的非定常气动力等。以航空航天领域中飞行器的绕流问题为例，飞行器在高速飞行时，周围流场的雷诺数很高，研究高雷诺数下的绕流特性对于飞行器的气动性能优化和结构设计至关重要。时间步长的选择对于模拟结果的准确性和计算效率有着重要影响。经过多次试算和分析，确定低雷诺数模拟时的时间步长\Deltat=0.001s，高雷诺数模拟时的时间步长\Deltat=0.0001s。时间步长过小会导致计算量大幅增加，计算时间过长；而时间步长过大则可能会使模拟结果出现较大误差，无法准确捕捉流场的瞬态变化。在低雷诺数下，由于流场变化相对缓慢，较大的时间步长0.001s能够在保证一定精度的前提下，提高计算效率。例如，在模拟低速水流绕过桥墩的情况时，低雷诺数下的水流速度较低，流场变化较为平稳，0.001s的时间步长可以较好地跟踪流场的变化。而在高雷诺数下，流场变化迅速，需要较小的时间步长0.0001s来准确捕捉流场中快速变化的信息，如涡旋的快速脱落和相互作用。比如在模拟高速气流绕过飞机机翼的情况时，高雷诺数下的气流速度快，流场变化剧烈，0.0001s的时间步长能够更精确地模拟机翼周围的复杂流场。迭代次数也是一个重要的模拟参数。在本次模拟中，设定每次计算的迭代次数为2000次。迭代次数的确定需要综合考虑计算精度和计算成本。迭代次数过少，计算结果可能无法收敛，导致模拟结果不准确；而迭代次数过多，则会浪费计算资源和时间。通过大量的模拟实验和分析，发现2000次的迭代次数能够使计算结果在合理的时间内达到收敛，满足模拟精度的要求。在模拟过程中，通过监测流场中的关键物理量，如柱体表面的压力系数、速度分布等，来判断计算是否收敛。当这些物理量在多次迭代后变化很小，趋于稳定时，认为计算已经收敛。例如，在模拟风绕过高楼的过程中，通过监测高楼表面的压力系数，当压力系数在2000次迭代后基本不再变化时，说明此时的模拟结果已经收敛，能够反映真实的流场情况。四、模拟结果与分析4.1流场特性分析4.1.1速度场分布通过大涡模拟，得到了不同时刻并列双方柱周围的速度矢量图和速度等值线图，这些图像为深入分析流场的速度分布规律和变化特征提供了直观而丰富的信息。在低雷诺数Re=1000、间距比L/D=1.5的工况下，从速度矢量图中可以清晰地观察到，流体在接近方柱时，速度逐渐发生变化。在方柱的迎流面，流体受到柱体的阻挡，流速急剧减小，形成明显的滞止区域。这是因为柱体的存在改变了流体的流动路径，使得流体在柱体前方堆积，动能转化为压力能，从而导致流速降低。在方柱的侧面，由于流体的粘性作用，靠近柱体壁面的流体速度迅速减小，形成边界层。边界层内的流体速度梯度较大，这是由于壁面的无滑移条件，使得壁面处的流体速度为零，而随着离壁面距离的增加，流体速度逐渐恢复到主流速度。在双方柱之间的区域，流体的流动受到双方柱的强烈干扰，形成复杂的流动结构。可以看到流体在双方柱之间来回振荡，速度方向不断变化，形成了多个小尺度的回流区域。这是因为双方柱之间的间距较小，流体在通过该区域时受到双方柱的双重阻挡和边界层的相互作用，导致流动变得不稳定。在方柱的尾流区域，流体速度逐渐恢复，但仍然存在明显的速度亏损。尾流中的速度分布呈现出不对称性，靠近双方柱一侧的速度亏损较大，这是由于双方柱之间的干扰使得尾流中的涡旋结构更加复杂，能量耗散更大。从速度等值线图上能够更直观地了解速度的分布情况。在柱体附近，速度等值线密集，表明速度变化剧烈。这是因为柱体表面的边界层和尾流区域存在较大的速度梯度，使得速度等值线在这些区域紧密聚集。在双方柱之间，速度等值线呈现出不规则的形状，反映了该区域复杂的流动特性。随着离柱体距离的增加，速度等值线逐渐稀疏，说明速度变化逐渐减小，流场逐渐趋于均匀。当间距比增大到L/D=4.0时，双方柱之间的干扰作用明显减弱。在速度矢量图中，双方柱之间的流体流动相对较为顺畅，回流区域明显减小。这是因为间距的增大使得双方柱之间的相互影响范围减小，流体在通过双方柱之间时受到的阻挡和干扰减弱，流动更加接近单方柱绕流的情况。在速度等值线图上，双方柱之间的速度等值线分布更加均匀，不再像小间距比时那样呈现出明显的不规则形状。在高雷诺数Re=100000下，流场的速度分布与低雷诺数时有明显的差异。由于惯性力的主导作用，流场中出现了强烈的湍流现象。在速度矢量图中，可以看到流体的速度方向和大小在不断快速变化，呈现出高度的随机性。在柱体表面和尾流区域，存在大量的小尺度涡旋，这些涡旋的旋转和运动使得流体的速度分布更加复杂。在速度等值线图上，高雷诺数下的速度等值线更加密集且分布更加不规则，反映了流场中速度的剧烈变化和高度的不均匀性。与低雷诺数相比，高雷诺数下尾流中的速度亏损范围更大，这是由于湍流的发展导致能量耗散加剧，尾流中的流体需要更长的距离才能恢复到主流速度。4.1.2涡量场分布涡量场是研究流体流动特性的重要物理量，它能够直观地反映涡旋的生成、发展、脱落和相互作用过程。通过大涡模拟得到的涡量云图，为深入分析涡量分布与流场特性的关系提供了关键依据。在低雷诺数Re=1000、间距比L/D=1.5时，从涡量云图中可以清晰地观察到涡旋的生成过程。在方柱的拐角处，由于流体的流速急剧变化，剪切应力增大，导致涡量迅速产生。随着流体的流动，这些初始生成的小涡旋逐渐合并、发展，形成较大的涡旋。在双方柱之间的区域，由于双方柱的相互干扰，涡旋的发展和相互作用变得极为复杂。可以看到多个涡旋相互缠绕、合并，形成了复杂的涡旋结构。在方柱的尾流中，涡旋呈现出周期性的脱落现象。当一个涡旋脱落后，在柱体表面又会重新生成新的涡旋，如此循环往复，形成了稳定的涡街。通过对涡量云图的分析，可以计算出涡旋的脱落频率，发现其与理论预测值基本相符。当间距比增大到L/D=4.0时，双方柱之间的涡旋相互作用明显减弱。在涡量云图中，双方柱之间的涡旋数量减少，且相互之间的干扰程度降低。这是因为间距的增大使得双方柱之间的相互影响范围减小，涡旋在各自的柱体周围发展，相互之间的作用相对较弱。在尾流区域，涡街的形态更加规则，涡旋的脱落频率也相对稳定。在高雷诺数Re=100000下，涡量场的分布更加复杂。由于湍流的强烈发展，流场中存在各种尺度的涡旋。在柱体表面附近，小尺度涡旋不断生成，这些小尺度涡旋在粘性力和惯性力的作用下，迅速扩散和相互作用。在尾流区域，大尺度涡旋和小尺度涡旋相互交织，形成了复杂的湍流结构。高雷诺数下涡旋的脱落频率更高，且涡旋的脱落过程更加不规则。这是由于湍流的脉动特性使得柱体表面的压力分布更加不均匀，从而导致涡旋的脱落更加随机。通过对不同雷诺数和间距比下涡量场的分析，可以发现涡量分布与流场特性之间存在密切的关系。涡旋的生成、发展和脱落过程直接影响着流场的速度分布、压力分布以及柱体的受力情况。在小间距比下，双方柱之间的涡旋相互作用强烈，导致流场的复杂性增加，柱体所受的气动力也更加不稳定。而在大间距比下，双方柱之间的干扰减弱，流场相对较为稳定。高雷诺数下，湍流的发展使得涡量场更加复杂，流场的能量耗散加剧，柱体所受的气动力也呈现出更强的脉动特性。4.2气动力特性分析4.2.1阻力系数与升力系数通过大涡模拟，精确计算得到了不同间距比下双方柱的阻力系数C_D和升力系数C_L随时间的变化曲线，这些曲线为深入分析气动力特性提供了关键数据。阻力系数C_D和升力系数C_L的计算公式如下：C_D=\frac{F_D}{\frac{1}{2}\rhov^2A}C_L=\frac{F_L}{\frac{1}{2}\rhov^2A}其中，F_D为阻力，F_L为升力，\rho为流体密度，v为来流速度，A为柱体的迎风面积。在低雷诺数Re=1000时，当间距比L/D=1.5，从阻力系数随时间变化曲线可以看出，阻力系数呈现出明显的周期性波动。这是因为在小间距比下，双方柱之间的干扰强烈，涡旋的生成、发展和脱落过程呈现出周期性，导致柱体所受的阻力也随之周期性变化。阻力系数的波动幅值较大，平均值约为[X1]。随着间距比逐渐增大到L/D=4.0，阻力系数的波动幅值逐渐减小，平均值约为[X2]。这是由于间距的增大使得双方柱之间的干扰减弱，涡旋的相互作用相对较弱，流场更加稳定，从而导致阻力系数的波动减小。升力系数随时间变化曲线同样显示出周期性波动。在小间距比L/D=1.5时，升力系数的波动频率较高，幅值也较大，这是由于双方柱之间复杂的涡旋相互作用导致柱体表面的压力分布不均匀，从而产生较大的升力波动。随着间距比的增大，升力系数的波动频率逐渐降低，幅值也逐渐减小，表明双方柱之间的干扰对升力的影响逐渐减弱。在高雷诺数Re=100000下，阻力系数和升力系数的变化规律与低雷诺数时有所不同。由于湍流的强烈发展，阻力系数和升力系数的波动更加剧烈，呈现出明显的非周期性。在小间距比L/D=1.5时，阻力系数和升力系数的波动幅值都很大，且波动中包含了丰富的高频成分。这是因为高雷诺数下，流场中的小尺度涡旋大量存在，这些小尺度涡旋的快速生成、发展和相互作用使得柱体所受的气动力更加不稳定。随着间距比的增大，虽然阻力系数和升力系数的波动仍然较为剧烈，但幅值有所减小，高频成分也相对减少。这说明间距的增大在一定程度上能够缓解高雷诺数下湍流对气动力的影响，使流场相对趋于稳定。综合来看，间距比和雷诺数对阻力系数和升力系数都有着显著的影响。小间距比会增强双方柱之间的干扰，导致气动力系数的波动增大；而大间距比则使双方柱之间的干扰减弱，气动力系数的波动减小。高雷诺数下，湍流的作用使得气动力系数的波动更加复杂和剧烈，而低雷诺数下，气动力系数的波动相对较为规则和稳定。这些变化规律对于深入理解并列双方柱绕流的气动力特性，以及在工程应用中准确评估柱体所受的气动力具有重要的指导意义。4.2.2气动力频谱分析为了深入研究气动力的频率特性，对不同间距比下双方柱的气动力进行了频谱分析。通过快速傅里叶变换（FFT），将气动力随时间的变化信号转换到频域，得到了气动力的频谱图。在频谱分析中，斯特劳哈尔数（St）是一个重要的无量纲参数，它与涡脱落频率密切相关，其定义为：St=\frac{fD}{v}其中，f为涡脱落频率，D为柱体边长，v为来流速度。在低雷诺数Re=1000时，对于间距比L/D=1.5的情况，从气动力频谱图中可以清晰地看到，升力系数频谱在频率f_1处出现明显的峰值，对应的斯特劳哈尔数St_1约为[X3]。这表明在该工况下，涡旋以频率f_1周期性地从柱体表面脱落，形成稳定的涡街，而升力系数的波动主要是由涡脱落引起的。阻力系数频谱在频率f_2处也有峰值，且f_2约为f_1的两倍。这是因为每次涡脱落会导致柱体表面压力分布的两次变化，从而使得阻力系数的变化频率是升力系数变化频率的两倍。随着间距比增大到L/D=4.0，升力系数频谱的峰值频率略有降低，对应的斯特劳哈尔数St_2约为[X4]。这说明间距的增大使得涡脱落频率降低，涡街的稳定性增强。阻力系数频谱的峰值频率同样约为升力系数频谱峰值频率的两倍。在高雷诺数Re=100000下，气动力频谱呈现出更为复杂的特征。由于湍流的影响，频谱中出现了多个峰值，且频率分布范围更广。在小间距比L/D=1.5时，升力系数频谱除了在低频段有明显的峰值外，在高频段也存在多个较小的峰值。低频段的峰值对应的涡脱落频率是流场中的主导频率，而高频段的峰值则是由小尺度涡旋的快速运动和相互作用产生的。阻力系数频谱同样在多个频率处出现峰值，且与升力系数频谱的频率关系较为复杂，不再像低雷诺数时那样呈现简单的两倍关系。这是因为高雷诺数下湍流的复杂性导致气动力的变化更加不规则，阻力系数的变化不仅与涡脱落有关，还受到小尺度涡旋的强烈影响。随着间距比的增大，气动力频谱中的高频成分相对减少，主导频率的峰值更加突出。这表明间距的增大有助于抑制小尺度涡旋的产生和发展，使流场中的主导涡脱落频率更加明显，气动力的频率特性相对简化。通过对不同雷诺数和间距比下的气动力频谱分析，可以得出结论：斯特劳哈尔数与涡脱落频率密切相关，能够准确反映流场中涡旋的脱落特性。在低雷诺数下，气动力的频率特性相对简单，主要由涡脱落频率主导；而在高雷诺数下，由于湍流的作用，气动力的频率特性变得复杂，包含了多个频率成分。间距比的变化会影响涡脱落频率和气动力的频谱特性，小间距比下，气动力频谱更为复杂，高频成分较多；大间距比下，气动力频谱相对简单，主导频率更加突出。这些研究结果对于进一步理解并列双方柱绕流的非定常特性，以及在工程应用中预测柱体的振动和噪声等问题具有重要的参考价值。4.3间距比的影响4.3.1对流场结构的影响通过大涡模拟得到的不同间距比下的流场图，能够清晰地观察到间距比变化对流场结构产生的显著影响。在低雷诺数Re=1000时，当间距比L/D=1.5，流场呈现出高度的复杂性。从速度矢量图中可以看到，在双方柱之间，流体的流动受到强烈的干扰，形成了复杂的回流和涡旋结构。由于双方柱之间的间距较小，流体在通过该区域时受到双方柱的双重阻挡和边界层的相互作用，导致流体的流动方向不断改变，形成了多个小尺度的回流区域。在涡量云图中，双方柱之间的涡旋相互作用强烈，多个涡旋相互缠绕、合并，形成了复杂的涡旋结构。这些涡旋的存在进一步加剧了流场的复杂性，使得流场中的速度分布和压力分布更加不均匀。随着间距比增大到L/D=4.0，双方柱之间的干扰作用明显减弱。速度矢量图显示，双方柱之间的流体流动相对较为顺畅，回流区域明显减小。这是因为间距的增大使得双方柱之间的相互影响范围减小，流体在通过双方柱之间时受到的阻挡和干扰减弱，流动更加接近单方柱绕流的情况。在涡量云图中，双方柱之间的涡旋数量减少，且相互之间的干扰程度降低。涡旋主要集中在各自的柱体周围发展，相互之间的作用相对较弱。这表明间距比的增大能够使流场结构更加稳定，减少流场中的复杂性。在高雷诺数Re=100000下，流场结构同样受到间距比的显著影响。当间距比L/D=1.5时，由于湍流的强烈发展，流场中存在各种尺度的涡旋，且涡旋的相互作用更加复杂。在柱体表面附近，小尺度涡旋不断生成，这些小尺度涡旋在粘性力和惯性力的作用下，迅速扩散和相互作用。在尾流区域，大尺度涡旋和小尺度涡旋相互交织，形成了复杂的湍流结构。随着间距比增大到L/D=4.0，虽然流场中仍然存在大量的涡旋，但涡旋的相互作用相对减弱，流场的复杂性有所降低。这是因为间距的增大有助于抑制小尺度涡旋的产生和发展，使流场中的主导涡脱落频率更加明显，流场结构相对简化。综上所述，间距比的变化对并列双方柱绕流的流场结构有着重要影响。小间距比会增强双方柱之间的干扰，导致流场结构更加复杂；而大间距比则使双方柱之间的干扰减弱，流场结构相对稳定。这种影响在低雷诺数和高雷诺数下均有体现，但高雷诺数下由于湍流的作用，流场结构的变化更加复杂。4.3.2对气动力特性的影响间距比对并列双方柱的气动力特性有着显著的影响，通过对不同间距比下阻力系数、升力系数及其频谱特性的分析，可以深入了解这种影响规律。在低雷诺数Re=1000时，随着间距比的变化，阻力系数和升力系数呈现出明显的变化规律。当间距比L/D=1.5时，阻力系数呈现出较大的波动，平均值约为[X1]。这是因为小间距比下双方柱之间的干扰强烈，涡旋的生成、发展和脱落过程呈现出周期性，导致柱体所受的阻力也随之周期性变化。升力系数同样波动较大，其波动频率较高，幅值也较大。这是由于双方柱之间复杂的涡旋相互作用导致柱体表面的压力分布不均匀，从而产生较大的升力波动。随着间距比逐渐增大，阻力系数和升力系数的波动幅值逐渐减小。当间距比增大到L/D=4.0时，阻力系数的平均值约为[X2]，波动幅值明显减小。这是因为间距的增大使得双方柱之间的干扰减弱，涡旋的相互作用相对较弱，流场更加稳定，从而导致阻力系数的波动减小。升力系数的波动频率和幅值也显著降低，表明双方柱之间的干扰对升力的影响逐渐减弱。从气动力频谱分析来看，在低雷诺数下，间距比的变化会影响涡脱落频率。当间距比L/D=1.5时，升力系数频谱在频率f_1处出现明显的峰值，对应的斯特劳哈尔数St_1约为[X3]。随着间距比增大到L/D=4.0，升力系数频谱的峰值频率略有降低，对应的斯特劳哈尔数St_2约为[X4]。这说明间距的增大使得涡脱落频率降低，涡街的稳定性增强。在高雷诺数Re=100000下，间距比对气动力特性的影响更为复杂。由于湍流的强烈发展，阻力系数和升力系数的波动更加剧烈，呈现出明显的非周期性。在小间距比L/D=1.5时，阻力系数和升力系数的波动幅值都很大，且波动中包含了丰富的高频成分。这是因为高雷诺数下，流场中的小尺度涡旋大量存在，这些小尺度涡旋的快速生成、发展和相互作用使得柱体所受的气动力更加不稳定。随着间距比的增大，虽然阻力系数和升力系数的波动仍然较为剧烈，但幅值有所减小，高频成分也相对减少。当间距比增大到L/D=4.0时，气动力频谱中的高频成分相对减少，主导频率的峰值更加突出。这表明间距的增大有助于抑制小尺度涡旋的产生和发展，使流场中的主导涡脱落频率更加明显，气动力的频率特性相对简化。通过分析可以确定临界间距比，当间距比小于临界间距比时，双方柱之间的干扰强烈，气动力特性的波动较大；当间距比大于临界间距比时，双方柱之间的干扰减弱，气动力特性相对稳定。在本研究中，通过对不同间距比下的气动力特性进行分析，确定临界间距比约为[X5]。临界间距比的确定对于工程应用具有重要意义。在建筑设计中，对于高层建筑群的布局，如果间距小于临界间距比，可能会导致建筑所受的风荷载过大，影响建筑的稳定性和安全性。因此，在设计时需要合理规划建筑之间的间距，使其大于临界间距比，以降低风荷载的影响。在海洋工程中，对于海上风力发电装置等柱状结构的布置，也需要考虑临界间距比，以减少结构之间的相互干扰，提高装置的运行效率和安全性。4.4雷诺数的影响4.4.1不同雷诺数下的流场特性为深入探究雷诺数对并列双方柱扰流的影响，进行了不同雷诺数下的大涡模拟。在低雷诺数Re=1000和高雷诺数Re=100000两种典型工况下，对不同间距比的并列双方柱绕流流场特性展开研究，重点对比流场结构和涡量分布。在低雷诺数Re=1000时，流场中的黏性力起主导作用，流体的流动相对较为平稳。从速度场分布来看，流体在接近方柱时，速度变化较为缓慢，边界层发展较为稳定。在双方柱之间，由于间距的不同，流场结构呈现出明显的差异。当间距比L/D=1.5时，双方柱之间的干扰强烈，流体在该区域形成复杂的回流和涡旋结构，速度矢量图显示流体在双方柱之间来回振荡，速度方向不断变化。随着间距比增大到L/D=4.0，双方柱之间的干扰减弱，流体流动相对顺畅，回流区域明显减小。涡量场分布在低雷诺数下也表现出一定的规律。在方柱的拐角处，由于流速的急剧变化，涡量迅速产生。随着流体的流动，涡旋逐渐合并、发展，在双方柱之间形成复杂的涡旋结构。在小间距比下，双方柱之间的涡旋相互作用强烈，多个涡旋相互缠绕、合并；而在大间距比下，涡旋主要集中在各自的柱体周围发展，相互之间的干扰相对较弱。当雷诺数增大到Re=100000时，惯性力成为主导因素，流场中出现强烈的湍流现象。速度场分布变得更加复杂，流体的速度方向和大小在不断快速变化，呈现出高度的随机性。在柱体表面和尾流区域，存在大量的小尺度涡旋，这些涡旋的旋转和运动使得流体的速度分布更加不均匀。在双方柱之间，尽管间距比的变化仍然会影响流场结构，但由于湍流的作用，这种影响相对低雷诺数时更加复杂。高雷诺数下的涡量场分布也与低雷诺数时有显著差异。由于湍流的强烈发展，流场中存在各种尺度的涡旋，小尺度涡旋不断生成、扩散和相互作用。在尾流区域，大尺度涡旋和小尺度涡旋相互交织，形成了复杂的湍流结构。涡旋的脱落频率更高，且脱落过程更加不规则。通过对比不同雷诺数下的流场结构和涡量分布，可以发现雷诺数对流场的影响机制主要体现在以下几个方面：随着雷诺数的增大，惯性力逐渐增强，流体的流动变得更加不稳定，容易产生湍流现象。湍流的出现使得流场中的涡旋结构更加复杂，小尺度涡旋的生成和相互作用加剧，从而导致流场的速度分布和压力分布更加不均匀。雷诺数的变化还会影响双方柱之间的干扰程度。在低雷诺数下，黏性力的作用使得双方柱之间的干扰相对较弱，流场结构相对稳定；而在高雷诺数下，惯性力的主导作用使得双方柱之间的干扰更加复杂，流场结构的变化更加剧烈。4.4.2对气动力特性的影响雷诺数的变化对并列双方柱的气动力特性有着显著的影响，包括气动力系数和频谱特性。通过大涡模拟计算得到的不同雷诺数下双方柱的阻力系数C_D和升力系数C_L随时间的变化曲线，以及气动力的频谱图，能够深入分析这种影响。在低雷诺数Re=1000时，阻力系数和升力系数的变化相对较为规则。阻力系数呈现出明显的周期性波动，这是由于涡旋的周期性脱落导致柱体所受的阻力也随之周期性变化。升力系数同样显示出周期性波动，其波动频率与涡脱落频率相关。随着间距比的变化，阻力系数和升力系数的波动幅值和平均值也会发生改变。在小间距比下，双方柱之间的干扰强烈，气动力系数的波动幅值较大；而在大间距比下，干扰减弱，气动力系数的波动幅值较小。从气动力频谱分析来看，低雷诺数下升力系数频谱在涡脱落频率处出现明显的峰值，阻力系数频谱的峰值频率约为升力系数频谱峰值频率的两倍。这是因为每次涡脱落会导致柱体表面压力分布的两次变化，从而使得阻力系数的变化频率是升力系数变化频率的两倍。当雷诺数增大到Re=100000时，由于湍流的强烈发展，阻力系数和升力系数的波动更加剧烈，呈现出明显的非周期性。在小间距比下，阻力系数和升力系数的波动幅值都很大，且波动中包含了丰富的高频成分。这是因为高雷诺数下，流场中的小尺度涡旋大量存在，这些小尺度涡旋的快速生成、发展和相互作用使得柱体所受的气动力更加不稳定。随着间距比的增大，虽然阻力系数和升力系数的波动仍然较为剧烈，但幅值有所减小，高频成分也相对减少。高雷诺数下的气动力频谱呈现出更为复杂的特征。频谱中出现了多个峰值，且频率分布范围更广。除了在低频段有明显的峰值对应于主导的涡脱落频率外，在高频段也存在多个较小的峰值，这些高频峰值是由小尺度涡旋的快速运动和相互作用产生的。阻力系数频谱与升力系数频谱的频率关系也不再像低雷诺数时那样呈现简单的两倍关系，这是由于高雷诺数下湍流的复杂性导致气动力的变化更加不规则。综合不同雷诺数下的模拟结果，可以总结出雷诺数与气动力之间的定量关系。随着雷诺数的增大，气动力系数的波动幅值增大，波动的非周期性增强，频谱中的高频成分增多。在相同的间距比下，高雷诺数时柱体所受的气动力更加不稳定，能量耗散也更大。这种定量关系对于工程应用中准确评估柱体在不同雷诺数下所受的气动力具有重要的参考价值。例如，在建筑结构设计中，需要考虑不同风速（对应不同雷诺数）下建筑所受的风荷载，通过了解雷诺数与气动力的关系，可以更准确地预测建筑的受力情况，从而优化建筑结构设计，提高建筑的安全性和稳定性。在海洋工程中，对于海上平台的设计，也需要考虑海水流速（对应不同雷诺数）对平台结构的作用力，以确保平台在复杂的海洋环境中能够安全运行。五、结果验证与讨论5.1与实验结果对比为了验证本文大涡模拟结果的准确性和可靠性，将模拟得到的流场特性和气动力特性与已有的实验数据进行了详细对比。在对比过程中，重点关注了速度分布、涡量分布、阻力系数和升力系数等关键参数。在速度分布方面，选取了[具体实验文献]中的实验数据进行对比。该实验采用粒子图像测速（PIV）技术，对不同间距比下并列双方柱绕流的流场速度进行了精确测量。将模拟得到的速度分布云图与实验测量结果进行对比，发现在低雷诺数和高雷诺数下，模拟结果与实验数据在整体趋势上具有较好的一致性。在低雷诺数Re=1000、间距比L/D=1.5时，模拟得到的双方柱之间和尾流区域的速度分布与实验测量结果相似，都呈现出复杂的流动结构和速度变化。在高雷诺数Re=100000、间距比L/D=3.0时，模拟结果同样能够准确地反映出实验中观察到的速度分布特征，如柱体表面附近的边界层速度变化和尾流中的速度亏损情况。然而，在一些局部区域，模拟结果与实验数据仍存在一定的差异。在柱体拐角处，由于模拟中对湍流的处理和实验中的实际情况存在细微差别，导致模拟得到的速度值与实验测量值略有偏差。对于涡量分布，与[另一实验文献]的实验结果进行了对比。该实验利用激光多普勒测速（LDV）技术结合流动显示方法，获得了流场中的涡量分布信息。从对比结果来看，模拟得到的涡量云图与实验图像在涡旋的生成位置、发展过程和脱落特性等方面具有较高的相似度。在低雷诺数下，模拟结果能够准确地捕捉到涡旋在方柱拐角处的生成以及在双方柱之间的相互作用过程，与实验中观察到的现象一致。在高雷诺数下，虽然模拟结果能够反映出流场中复杂的涡量分布特征，但由于湍流的复杂性和实验测量误差的存在，在一些小尺度涡旋的细节上，模拟结果与实验数据存在一定的差异。在气动力特性方面，将模拟得到的阻力系数和升力系数与[相关实验研究]的实验数据进行了对比。该实验通过力传感器精确测量了不同雷诺数和间距比下双方柱所受的气动力，并计算得到了阻力系数和升力系数。对比结果显示，在低雷诺数和高雷诺数下，模拟得到的阻力系数和升力系数随时间的变化曲线与实验数据在趋势上基本一致。在低雷诺数Re=1000时，模拟得到的阻力系数和升力系数的波动周期和幅值与实验数据较为接近，能够准确地反映出涡旋脱落对气动力的影响。在高雷诺数Re=100000下，虽然由于湍流的非定常性和复杂性，模拟结果与实验数据在气动力系数的波动细节上存在一定差异，但整体的变化趋势和量级是相符的。通过对速度分布、涡量分布、阻力系数和升力系数等参数的对比分析，可以得出结论：本文采用的大涡模拟方法在研究并列双方柱绕流问题时，能够较为准确地预测流场特性和气动力特性。模拟结果与实验数据在整体趋势上的一致性，验证了模拟方法和参数设置的合理性和有效性。然而，在一些局部区域和细节方面，仍存在一定的差异，这可能是由于模拟中对湍流模型的简化、网格分辨率的限制以及实验测量误差等多种因素共同作用的结果。在未来的研究中，可以进一步优化湍流模型、提高网格分辨率，并结合更精确的实验测量技术，以减小模拟结果与实际情况的差异，提高模拟的准确性和可靠性。5.2模拟结果的不确定性分析在大涡模拟过程中，存在多种因素可能导致模拟结果产生不确定性。其中，子网格模型的选择对模拟结果有着至关重要的影响。不同的子网格模型基于不同的假设和理论，对小尺度涡旋的模拟方式和精度各不相同。如前文所述，Smagorinsky模型形式简单，计算效率高，但由于其常数固定，在模拟复杂流场时难以准确反映小尺度涡旋的动态变化。在模拟并列双方柱绕流这种复杂流场时，固定的常数可能导致对柱体附近和尾流区域小尺度涡旋的模拟偏差，进而影响对整个流场特性的准确描述。而动力亚网格尺度（SGS）模型虽然能够根据流场的局部特性动态调整亚网格尺度的参数，在一定程度上提高了模拟精度，但它也并非完美无缺。在一些极端工况下，如高雷诺数且强湍流的情况下，动力SGS模型对小尺度涡旋的模拟仍可能存在一定的误差。这是因为在这种复杂工况下，小尺度涡旋的生成、发展和相互作用更加剧烈和复杂，模型的假设和理论可能无法完全准确地描述这些物理过程。网格精度同样是影响模拟结果不确定性的关键因素。在大涡模拟中，网格分辨率直接决定了对大尺度涡旋的捕捉能力。若网格过粗，一些重要的大尺度涡旋结构可能无法被准确解析，导致模拟结果丢失关键的流动信息。在模拟高雷诺数下并列双方柱绕流时，流场中存在各种尺度的涡旋，若网格分辨率不足，可能无法准确捕捉到大尺度涡旋的脱落和相互作用过程，从而使模拟得到的阻力系数、升力系数等气动力特性与实际情况产生较大偏差。虽然进行了网格无关性分析并选择了中等网格，但中等网格在某些局部区域可能仍然无法满足高精度的要求。在柱体拐角等流场变化剧烈的区域，即使是中等网格，其分辨率也可能不足以精确捕捉到流场的细微变化，从而引入一定的不确定性。此外，模拟参数的设置也会对结果的不确定性产生影响。时间步长的选择直接关系到模拟结果的时间精度。时间步长过大，可能无法准确捕捉流场的瞬态变化，导致模拟结果出现偏差。在高雷诺数下，流场变化迅速，若时间步长设置不合理，可能会错过一些重要的涡旋脱落和相互作用事件，使得模拟得到的气动力系数的波动特性与实际情况不符。迭代次数也会影响模拟结果。迭代次数不足，计算结果可能无法收敛，导致模拟结果不准确。在实际模拟中，虽然设定了2000次的迭代次数，但对于一些复杂的流场情况，这一迭代次数可能并不足以使所有物理量都完全收敛，从而增加了模拟结果的不确定性。边界条件的处理也可能引入不确定性。在实际应用中，边界条件的设定往往是基于一定的假设和简化，与实际情况可能存在差异。在设置入口边界条件时，假设入口来流为均匀流，但在实际工程中，入口来流可能存在一定的湍流度和不均匀性，这可能导致模拟结果与实际情况产生偏差。壁面边界条件的处理也可能存在不确定性。虽然采用了无滑移边界条件，但在实际的流体-固体相互作用中，壁面的微观结构和表面性质可能会对流体的流动产生影响，而这些因素在模拟中往往难以完全准确地考虑。为了减小模拟结果的不确定性，可以采取一系列措施。在子网格模型方面，可以进一步研究和改进模型，结合多种模型的优点，开发更精确、更具适应性的子网格模型。在网格划分上，可以采用自适应网格技术，根据流场的变化动态调整网格分辨率，在流场变化剧烈的区域自动加密网格，以提高对流动细节的捕捉能力。对于模拟参数的设置，可以通过多次试算和敏感性分析，确定最优的时间步长和迭代次数。在边界条件处理上，应尽可能收集更准确的实际工况数据，对边界条件进行更合理的设定，同时考虑采用更复杂的边界条件模型，以更准确地描述边界处的物理现象。5.3结果讨论与实际应用启示通过大涡模拟对并列双方柱绕流问题的研究，获得了丰富的模拟结果，这些结果具有深刻的物理意义。在流场特性方面，速度场和涡量场的分布清晰地展示了流体在绕流过程中的复杂运动。速度场中的滞止区域、边界层、回流区域以及速度亏损等现象，反映了流体与柱体之间的相互作用以及能量的传递和耗散。涡量场中涡旋的生成、发展、脱落和相互作用过程，揭示了流场的非定常特性和湍流的产生机制。在低雷诺数下，黏性力主