错列角度对双圆柱涡激振动影响的数值模拟研究1)

引言

涡激振动问题发生在很多实际工程中，比如海洋立管、海底管线、高耸建筑等在受到风或海流的作用时，会因涡旋脱落而导致结构振动.和静止圆柱绕流不同，当涡旋脱落频率接近柱体自振频率时，柱体会发生“锁定”现象[1].“锁定”现象的发生会使结构物发生疲劳破坏，进而缩短结构物的寿命.双圆柱涡激振动作为多柱体涡激振动的基础，除柱体和流体本身的特性外，柱体布置方式也会影响柱体的振动响应和流场状态.如图1所示，仅在双圆柱涡激振动中就存在3种布置方式：串列(α=0°)、并列(α =90°)和错列(0°＜ α ＜ 90°).其中，α是两柱体中心连线和来流方向之间的夹角(锐角)，下文中将其简称为“错列角度”.

  

图1 双柱体布置示意图Fig.1 Sketch of layout for two circular cylinders

涡激振动过程涉及湍流、流动分离、涡旋拓扑、物体振动和流固耦合等复杂问题，其数值分析的难点在于动边界的处理和流场的精确捕捉.Tonmurt[2]利用任意拉欧数值方法对雷诺数Re=100～140范围的单圆柱涡激振动问题进行了系统研究，发现物体在运动时柱体表面网格的正交性不易保持,容易出现网格畸变问题.近些年，对单圆柱涡激振动的研究集中在来流特性对振动结果的影响.其中，剪切流对单圆柱涡激振动影响的相关研究可参考文献[3-5]，周期性振荡流对单圆柱涡激振动的影响可参考文献[6-7].多柱体方面，Mittal等[8]利用重叠网格技术研究了雷诺数Re=100和1000串列及错列布置双圆柱涡激振动问题，发现在同一雷诺数下，错列布置下游圆柱受到的阻力Cd比串列布置大.Prasanth和Mittal[9]使用贴体动网格技术研究了雷诺数Re=100，错列角度α≈7°双圆柱涡激振动下游圆柱的运动轨迹和尾涡模态.研究发现，只有当折合速度7.1＜Ur＜8.3时，下游圆柱运动轨迹才是“8”字型.Borazjani和Sotiropoulos[10]利用浸入边界法处理流固耦合问题，数值模拟了串列双圆柱涡激振动问题，发现串列布置柱体的振幅和升力要比单柱体大，且在不同的折合速度下分别呈现出“2S”和“2P”两种涡旋脱落模态.吴晓笛等[11]基于浸入边界--多松弛时间格子波尔兹曼通量求解法的弱耦合算法对串列双圆柱双自由度涡激振动问题进行了研究.郭晓玲等[12]采用任意拉格朗日--欧拉动网格方法对雷诺数Re=150串列双圆柱涡激振动问题进行了数值模拟.Williamson[13]通过模型试验对并列双圆柱涡激振动进行了研究，发现尾涡脱落存在两种方式，即同相同步.Chen等[14]对雷诺数Re=100不同间距比下并列双圆柱涡激振动进行了数值模拟，发现当间距比大于4.0时，两圆柱之间的相互作用非常弱.可见，关于双圆柱涡激振动的研究大多是在固定两柱体的布置方式下研究不同折合速度时的柱体响应.而实际工程中，比如，两个近距离布置的高耸烟囱、海洋立管等，在不同的来流方向下就会呈现出不同错列角度的布置方式.如何设计常风向或者常流向下最安全的错列角度α是实际工程中必须考虑的问题.基于此，本文拟对不同折合速度Ur和错列角度0°≤α≤90°下双圆柱涡激振动进行数值模拟，系统分析柱体升阻力、位移响应和柱体后方尾涡模态.

本文基于浸入边界方法的流固耦合技术，采用固定笛卡尔正交网格，在计算过程中无需更新网格，避免了重叠网格技术中的大量信息交换和任意拉欧方法中的网格畸变问题，同时可保证模型的高效率，使得模型在处理大位移的流固耦合计算中具有强大优越性.为精确捕捉流场信息，在以前的工作中已将自主研发的紧致插值曲线方法模型应用于圆柱绕流问题[15-16]和强迫振荡圆柱绕流问题[17]的研究中，初步验证了该模型在处理复杂流动问题的准确性.本文拟进一步拓展该模型的应用范围，对不同折合速度Ur和错列角度α下双圆柱涡激振动进行数值模拟，系统分析特定折合速度Ur下柱体升阻力、位移响应和柱体后方尾涡模态随错列角度的变化规律，以期为相关海洋工程设计提供参考数据.

1 数值方法

涡激振动是一种流固耦合问题，需要考虑流场、结构物和流固耦合三个过程.其中，流场模型考虑二维不可压缩流体，控制方程为N-S方程和质量守恒方程(连续性方程)，其张量形式如下

 

其中，i,j是笛卡尔张量符号，t,uj,p,ρ分别表示时间、流速、密度和动水压力，fi是力源项，Sij是黏性项，由下式定义

 

式中，µ表示黏性.

柱体结构可简化为一个质量--弹簧--阻尼系统，考虑横流向和顺流向的运动，其控制方程为

 

Description about acupoints and meridians were standardized according to the WHO Standard Acupuncture Point Locations in the Western Pacific Region.13

在直角坐标系下建立模型，采用多相流理论将计算域分为固相、液相两部分.定义体积函数Φn来区分固相和液相(n=1，代表液相；n=2，代表固相)并捕捉固--液界面.Φn满足如下公式

 

式中，Φ12=Φ1+Φ2，并在一个网格内满足Φ1+Φ2=1.则网格内的流体特性可用下式表示

 

从图6(c)可看出，在3个特征工况下，错列角度较小(0°＜α≤45°)时，下游圆柱的顺流向振动响应都比上游圆柱剧烈；错列角度较大(α＞45°)时两柱体的顺流向振动响应相差较小，这与Prasanth等[26]研究的结论一致.

本文所采用的CIP方法模型，核心部分在于对流体控制方程(1)和(2)进行离散，过程分为3步：(1)采用CIP方法[18]求解对流项；(2)采用中心差分法求解压力项外的非对流项；(3)采用逐次超松驰迭代(successive over relaxation，SOR)方法进行压力速度耦合.流固耦合边界采用浸入边界法(immersed boundary method，IBM)[19]处理.

  

图2 计算域及网格划分示意图Fig.2 Sketch of computational domain and mesh

通过对不同阻流比B(B=d/H，H为计算域宽度，单圆柱时d=D，双圆柱时d=2D−h，h柱为两圆柱在横流向投影的重叠长度，如图1所示)单双圆柱涡激振动问题数值模拟，得知当阻流比小于10%时，本数值模型中阻流比对涡激振动特征参数的影响可忽略.为保证结果可靠性，本文取阻流比B=5%.计算区域及网格划分如图2所示.上下边界的圆柱直径D=1,入口边界距上游圆柱圆心为10D，上边界、下边界距上游圆柱圆心均为20D，出口边界距上游圆柱圆心25D.对于双圆柱涡激振动而言，两柱体并列布置时阻流比达到最大值5%，可保证所有错列角度下阻流比的影响都可忽略.数值模拟参数Re=100,折合速度Ur=2.0～15.0，间隔 1.0，圆柱的质量比 m∗=4m/(ρfπD2)=2.0,阻尼比ζ∗=0.为更好地捕捉柱体的运动和其周围的复杂流动现象，对柱体近壁面网格进行加密，最小网格设为Δx=Δy=D/40.计算流场的边界条件和初始条件设置如下：初始速度场设置为从左向右并保持均匀来流速度u=U，v=0(u,v分别为x,y方向的速度)；初始压力场设定为零；出口边界为开边界(自由出流)；柱体表面采用无滑移边界；两侧壁面采用自由滑移边界.通过对不同柱间距T下错列双圆柱涡激振动进行数值模拟发现，当T＜3D时，两圆柱会出现碰撞现象.为保证两柱体在振动过程中相互影响较大且不出现碰撞现象，综合选取中心间距T=3D.

式中，m,c和k分别表示柱体的质量、阻尼参数和弹簧刚度，Fy表示垂直于流向方向柱体的合力，Fx表示流向方向柱体的合力.

2 模型验证

2.1 网格收敛性验证

为验证网格收敛性，本节选取雷诺数Re=200的单圆柱静止绕流作为算例，设置四套不同尺寸的网格(mesh1,mesh2,mesh3,mesh4)，其对应的最小网格尺寸分别为：Δx= Δy=D/50，D/40，D/25和D/10.将计算出的柱体阻力系数均值Cd.mean和无量纲涡旋脱落频率St和文献[20-22]的结果进行对比，如表1所示.可看出，mesh3和mesh4的阻力系数均值和无量纲涡旋脱落频率较其他结果偏小，mesh1和mesh2与其他文献的结果吻合较好，且mesh2的结果和文献[21]的结果相差最小.表中RE分别是阻力系数均值和无量纲涡脱频率相对mesh1结果的偏差.可看出，mesh4的相对偏差较大，都在10%以上.而mesh2结果的相对偏差都在2%以下，表明mesh2网格已具有收敛性.综合考虑精度及计算效率，本文采用mesh2对接下来的工况进行计算.

 

表1 单圆柱绕流结果Table 1 Result of flw past a singular circular cylinder

  

Grids Cd.mean(RE) St(RE)mesh1 407800 1.3135 0.1934 mesh2 92800 1.3046(0.67%) 0.1904(1.5%)mesh3 32175 1.2911(1.7%) 0.1855(4.08%)mesh4 12500 1.1715(10.8%) 0.1697(12.2%)Ref.[20] — 1.336 0.194 Ref.[21] — 1.3 0.196 Ref.[22] — 1.32 0.1922

2.2 单圆柱涡激振动

为便于阐述上文中3种情况下上下游柱体的振动响应、位移响应随错列角度变化的情况，本节选取Ur=3.0,5.0和13.0的工况分别作为上述3种情况的特征工况进行分析.图6给出了3个特征工况下两柱体横流向位移均方根(图6(a))、顺流向位移均值(图6(b))以及顺流向位移均方根(图6(c))随错列角度的变化情况.其中，横、顺流向位移均方根大小分别反映了柱体在横流向和顺流向振动的剧烈程度，顺流向位移均值大小反映了柱体在顺流向上的平均偏移程度.

  

图3 单圆柱涡激振动振幅比随折合速度变化情况Fig.3 Amplitude ratio versus reduced velocity in singular circular cylinder condition

3 结果与讨论

考虑到折合速度Ur对柱体的振动响应有较大影响，本节首先对折合速度Ur=2.0～13.0(间隔1.0)下不同错列角度的双圆柱涡激振动进行数值模拟，重点关注两柱体的阻力系数均值Cd.mean和升力系数均方根Cl.rms.为便于观察，按照图3中单圆柱振幅随折合速度的变化情况将折合速度分成(a)、(b)和(c)三个区域.其中，区域(a)内的折合速度小于锁定区间对应的折合速度，即Ur=2.0～3.0；区域(b)内的折合速度处于锁定区间内，即Ur=4.0～8.0；区域(c)内的折合速度大于锁定区间对应的折合速度，即Ur=9.0～13.0.当折合速度Ur=14.0和15.0时两柱体出现了碰撞，本文暂不予考虑.

3.1 柱体升阻力系数

图4和图5分别是折合速度Ur=2.0～13.0不同错列角度上下游圆柱的阻力系数均值和升力系数均方根随错列角度的变化情况.其中，(a)、(b)、(c)分别对应上文中的3个区域，“UC”代表上游圆柱，“DC”代表下游圆柱.

从图4和图5可看出，除了过渡区的折合速度(Ur=4.0和9.0)外，其他折合速度下柱体升阻力随错列角度的变化规律在同一区域内具有相似性.对于Ur=2.0～3.0，如图4(a)所示，α=0°时下游圆柱阻力几乎为零，说明上游圆柱尾流对下游圆柱的“屏蔽”作用[25]较强，即下游圆柱处在上游圆柱的尾流中，来流速度较小导致下游圆柱阻力较小.

原辅料出库主要面向制丝、卷包和咀棒车间，卷包车间由于辅料中转库的存在，大库到中转库属于移库操作，从中转库到机台的出库为生产出库，而到制丝和咀棒车间的出库直接为生产出库。

  

图4 折合速度Ur=2.0～13.0上下游圆柱的阻力系数均值Cd.mean随错列角度的变化情况Fig.4 Mean value of drag coefficient versus stagger angle for different reduced velocities

  

图5 折合速度Ur=2.0～13.0上下游圆柱的升力系数均方根Cl.rms随错列角度的变化情况Fig.5 Root mean square lift coefficients versus stagger angle for different reduced velocities

图7统计了3个特征工况中两柱体升力系数的主导频率，反映两柱体涡脱的主导频率.图8～图10分别是3个特征工况中不同错列角度下的瞬时流场涡量等值线，反映不同的尾涡脱落模态.

伊秉绶：阴胜国[29]从清初隶书家审美观念的对比中探讨伊秉绶的审美思想，认为伊秉绥的书学思想有3个要点：①取法雄强宽博、静穆方正风格的汉碑；②用笔的中实简省与结体的方正更易；③追求静穆而不失古拙的气息。总之，伊秉绥在审美观念的指引下，凭借自足的勇气和气魄，锐意进取，将隶书革新到底，成功的塑造了崇尚古朴和趋拙避巧的艺术形象，对后来碑派书法的发展产生了深远的影响。

对于Ur=9.0～13.0，如图4(c)所示，错列角度较小时上游圆柱对下游圆柱的“屏蔽”作用较弱，在错列角度α=0°时下游圆柱阻力不是最小值.如图5(c)所示，两柱体的升力随错列角度的增大呈“下凹”趋势，在α=45°左右取得最小值.若固定错列角度，对比不同折合速度下柱体的升阻力，可发现不同折合速度下上下游柱体的升阻力基本相等，说明折合速度对上下游柱体升阻力的影响都较小.

3.2 柱体振动响应、位移响应

本节中选取质量比m∗=2.0，雷诺数Re=100，折合速度Ur=2.0～15.0(间隔1.0)的单圆柱涡激振动工况作为算例，并将振幅比A∗=(Ymax−Ymin)/D与文献[20,23-25]结果相比较，如图3.其中，Ymax和Ymin分别是柱体在一个响应周期内横流向的最大、最小位移.如图3所示，当 Ur=5.0时，振幅比A∗=0.58，达到最大值，意味着圆柱进入频率锁定状态，与及春宁等[25]的结果吻合较好.随折合速度的增大，振幅比逐渐变小，在Ur=9.0时，振幅比A∗=0.1178，远小于锁定区间的振幅比，柱体脱离频率锁定区.通过对比可看出，柱体振幅比随折合速度变化的规律和现有文献结果吻合较好.

对Ur=3.0，从图8可看出错列角度较小时的“屏蔽”作用.α=0°和15°时，柱间流态呈现单钝体模式和剪切层重附着模式[27]，上游圆柱没有出现涡旋脱落.随错列角度的增大，上游圆柱开始出现涡旋脱落，柱后涡街数量增大，涡街形态逐渐变得规则，导致图7(a)中下游圆柱涡脱的两个主导频率逐渐接近，使得柱体升力和横流向位移均方根增大，解释了上下游圆柱升力和横流向位移均方根随错列角度增大呈单调增大的现象.

  

图6 不同折合速度上下游柱体位移响应随错列角度的变化情况Fig.6 Displacement response of both cylinders with stagger angle under different reduced velocity Ur

  

图7 不同折合速度上下游圆柱涡旋脱落主导频率随错列角度的变化情况Fig.7 The main vortex shedding frequency of both cylinders versus stagger angle under different reduced velocities

Ur=5.0时，上下游柱体的横流向位移均方根及顺流向位移均值随错列角度变化的规律分别和升、阻力的变化规律一致.其中，上游圆柱的横流向位移均方根和柱体升力在α=60°左右同时取得最小值，下游圆柱的横流向位移均方根和柱体升力在α=30°左右同时取得最大值，而后随错列角度的增大而减小.对于顺流向位移均值，下游圆柱在α=0°时最小，并随错列角度的增大先增大后减小，在α=60°左右取得最大值.而上游圆柱的顺流向位移均值随错列角度的变化波动较小，在阻力最小即α=60°时取得最小值.Ur=13.0时，上下游圆柱顺流向位移均值和所受阻力随错列角度变化规律一致，但柱体的横流向位移均方根不再和升力变化规律一致.下游圆柱升力在α=45°左右取得最小值时，其横流向位移均方根却不是最小值.总体来说，下游圆柱在错列角度较小时(α＜45°)其横流向位移均方根普遍较大，而错列角度较大时(α＞45°)横流向位移均方根普遍较小.在错列角度α＞45°后，随错列角度的增大下游圆柱所受升力逐渐增大，而横流向位移均方根却逐渐减小.这是由于本文折合速度Ur的改变是通过改变柱体的刚度k实现的.折合速度较大时柱体的刚度较小，此时折合速度并不在锁定区间，相应单柱体的横流向振动不剧烈.但是，错列布置方式使得下游圆柱在错列角度较小时容易受到上游圆柱尾涡的影响而在横流向产生较为剧烈的振动，使得下游圆柱横流向位移均方根在错列角度较小时普遍较大，下游圆柱的横流向振动主要由错列角度控制.上游圆柱的横流向位移均方根随错列角度的变化规律和柱体所受升力保持一致.

赛珍珠女士(Pearl S.Buck)将中国名著《水浒传》的书名译为《All Men Are Brothers》，单从字面看觉得该翻译没有表达出原著的意图，但是译者考虑到无法用“water”直译“水”，“浒”指“水边，离水稍远的岸上陆地”(waterside)，“传”可译为“novel”，如果只是将中文对应的英文串联起来则毫无意义。因此，译者有意识地误译让外国读者从书名也可以猜测出该著作有关侠义之事。

式中，γ表示密度ρ或黏滞系数µ.

对于全民学习共享平台而言，教学质量是平台社会服务能力的体现，对教学质量提出高要求，才能规范教学活动的各个环节，促进教学效果的优化。教学质量管理是全民学习共享平台的品质保证：

3.3 涡脱频率及尾涡模态

随错列角度的增大，“屏蔽”作用逐渐减弱，下游圆柱阻力明显增加，在α=90°时两柱体阻力相等.和下游圆柱相比，上游圆柱的阻力随错列角度的变化较小.如图5(a)所示，两柱体升力随错列角度的增大几乎都呈现单调增加的趋势.另外，若固定错列角度，对比不同折合速度下两柱体的升阻力，可发现上下游柱体的阻力基本相等，而升力出现差异，说明折合速度对上下游柱体的阻力影响较小，对升力的影响较大.对于Ur=4.0～8.0，如图4(b)所示，柱体阻力随错列角度变化规律上呈现的主要差异是“屏蔽”作用的弱化.如图5(b)所示，两柱体的升力随错列角度的变大不再单调增大，而是呈“上凸”趋势，在α=15°～30°取得最大值.若固定错列角度，对比不同折合速度下柱体的升阻力，可发现随折合速度的增大，上下游圆柱所受升阻力逐渐减小，说明折合速度对上下游柱体升阻力的影响都较大.

如图6(a)和图6(b)所示，Ur=3.0时，上下游圆柱的横流向位移均方根在α=0°时最小，随错列角度的增大两柱体在横流向的振动程度都在增大，在α=75°附近取得最大值，和柱体升力大小变化规律一致.由于上游圆柱“屏蔽”作用的影响，下游圆柱的顺流向位移均值在α=0°时几乎为零.随错列角度的增大，下游圆柱的顺流向偏移程度逐渐变大，并逐渐接近上游圆柱，总体上和柱体阻力大小变化规律一致.

  

图8 折合速度Ur=3.0不同错列角度下尾涡模态Fig.8 Instantaneous vorticity contour versus different stagger angle at Ur=3.0

  

图9 折合速度Ur=5.0不同错列角度下的尾涡模态Fig.9 Instantaneous vorticity contour versus different staggered angle at Ur=5.0

对 Ur=5.0，从图9可看出，在 α=0°～30°时，由于上下游圆柱横流向振动都比较剧烈，两柱体在横流向的位移差使得上游圆柱在错列角度较小时也有足够的空间产生涡旋脱落，下游圆柱在整个振动周期内会穿过上游圆柱尾涡，产生较强的相互作用，这和Ur=3.0的尾涡模态不同.在α=0°时上游圆柱的尾涡交替作用在下游圆柱两侧并和符号相同的涡融合，对下游圆柱起到一个“屏蔽”作用，但这种“屏蔽”作用要比Ur=3.0上游圆柱剪切层静态附着在下游圆柱时的“屏蔽”作用弱.“屏蔽”作用随错列角度的增大而变弱，上游圆柱脱落的涡旋逐渐只能作用在下游圆柱的下侧，促进了下游圆柱的涡旋脱落，使得下游圆柱升力在α=0°～30°的变化过程中出现了明显增大的趋势.虽然总体上上下游圆柱的涡脱频率随错列角度的增大几乎保持不变，但是柱体的升力却呈现“上凸”的变化规律，在α=30°处出现拐点.从图9可看到，这是因为在α=30°到45°转变时尾涡模态从单一涡街模式转变为双涡街模式，导致升力在此产生突变.α=45°时，上下游圆柱内侧的涡旋出现“撞击”的现象.即在整个振动周期内，上下游柱体横流向距离最小时两柱内侧涡旋会“撞击”，使得内侧将要脱落的两个涡旋“断裂”成小涡后向下游运动，小涡的脱落使得上下游圆柱升力幅值减小，导致升力减小.同时，从α=45°和60°的涡量图可以看到，这种“撞击”限制了上游圆柱的横流向振动，使得上游圆柱涡街变窄，而下游柱体的振动没有出现明显变化，尾涡脱落后形成两条平行涡街.这也解释了图6(a)中α=60°时上游柱体横流向位移均方根突然减小的现象.随错列角度的继续增大(α=75°和90°)，柱体间横流向间距变大，上述“撞击”现象所引起的抑制作用变弱，且两柱体顺流向距离的减小使得两柱尾涡出现“同相同步”和“反相同步”两种模态，导致上下游圆柱升力系数均方根、阻力系数均值及相应的位移响应逐渐相同.值得一提的是，和Ur=3.0并列布置下的柱后尾涡相比，此时“同相同步”尾涡模态持续的时间更长.Chen等[28]在研究并列双圆柱涡激振动时也发现了这种现象，Chen认为这是因为当两柱体横流向振幅较大时，“同相同步”状态可以使圆柱之间保持近似的距离，有效避免碰撞.这种情况下，“同相同步”状态对两柱体的振动响应是更稳定的状态.

对Ur=13.0，从图10可看出，α=0°～30°时，两柱体横流向振幅较小，尾涡模态和Ur=3.0的尾涡模态相似.在α=30°和45°时，刚度较小的特点使得下游圆柱在横流向的振动容易受上游圆柱尾涡脱落的影响，在横流向产生较大的位移.下游圆柱的振动对上游圆柱尾涡产生“抑制”作用，使得上游圆柱尾涡变得“细长”.且错列角度越大抑制作用越明显，进而导致上游圆柱在错列角度α=0°～45°的变化过程中升力出现减小的趋势.随着错列角度的继续增大，柱体间横向间距变大，这种影响逐渐变弱，涡街数量增大，尾涡也变得比较规则，出现了“同相同步”和“反相同步”两种模态.如图7(c)所示，α=0°时，两柱体的涡脱频率都有一个次高频存在，约为单柱体涡脱频率的3倍(f=0.3564).杨立红[29]认为这是由于柱体间相互作用较强烈导致升力中的二次谐波(3f)的振幅较大，不可忽略.α=60°时，涡脱模态发生转变，导致两柱体升力突然增大.横流向间距的增大使得上游圆柱发生明显的涡旋脱落，柱后涡街数量增大.在α=60°～90°的变化过程中，两柱体顺流向距离的减小和横流向间距的增加使得尾涡相互作用变弱，逐渐出现较为规则的涡街，主导频率间的差异逐渐减小.升力系数时程曲线中“拍”的周期因主导频率差异的减小而增大，导致柱体升力系数均方根随错列角度的增大而变大.特别的，α=90°时尾涡脱落只有“反相同步”一种模态.

  

图10 折合速度Ur=13.0不同错列角度下尾涡模态Fig.10 Instantaneous vorticity contour versus different staggered angle at Ur=13.0

4 结论

基于自主研发的CIP方法模型对雷诺数Re=100，不同错列角度及不同折合速度下的双圆柱涡激振动问题进行了研究，两柱体中心间距固定为T=3D，重点分析了上下游柱体的升阻力、位移响应、涡脱频率随错列角度变化的规律，并从涡脱模态变化的角度解释其原因.主要结论如下：

(1)错列布置时，上下游圆柱是否处于锁定状态仍由折合速度决定.不同的折合速度下，上下游柱体升阻力、振动响应、位移响应、涡脱频率及尾涡模态随错列角度的增大所呈现的规律不同，在实际工程设计中必须根据流速大小考虑较为安全的错列角度.

(2)当折合速度较小时(Ur=2.0～3.0)，上游圆柱尾涡对下游圆柱“屏蔽”作用较明显，柱体升阻力随错列角度的增大基本呈单调增大的趋势，且柱体升阻力的大小分别决定了柱体横流向位移均方根和顺流向位移均值的大小.错列角度固定时，折合速度的改变对上下游柱体阻力的影响较小，对升力的影响较大.

染料的广泛使用导致染料废水造成的环境污染日益严重[1,2]。染料废水污染组分复杂、毒性强，是一种难处理的工业废水[3,4]。亚甲基蓝就是一种使用最为广泛的染料。常用的染料废水的处理方法主要有沉淀法[5]、氧化法[6]、膜分离法[7]及吸附法[8]等。吸附法因其易操作、高效等优点而受到重视。

(3)当折合速度Ur=5.0～8.0时，上下游柱体处于“锁定”状态，横流向振动较剧烈，上游圆柱尾涡对下游圆柱“屏蔽”作用较弱，随错列角度的增大上下游柱体阻力变化较小，升力呈“上凸”趋势，在α=15°～30°左右取得最大值，且柱体升阻力的大小分别决定了柱体横流向位移均方根和顺流向位移均值的大小.错列角度固定时，折合速度的改变对上下游柱体升阻力的影响都较大.

经济发展水平(ED)：一般而言，工业化进程中经济发展水平越高，制造业结构越高级，本文使用人均国内生产总值(GDP)度量经济发展水平，其值越大，说明经济发展水平越高。

(4)当折合速度较大时(Ur=10.0～13.0)，上下游柱体不在“锁定”状态，但“刚度小”的特点使得下游圆柱横流向的振动在错列角度α=30°～45°时容易受上游圆柱尾涡脱落的影响，在横流向产生较大的振动，“屏蔽”作用较弱.随错列角度的增大上下游柱体阻力变化较小，升力呈“下凹”趋势，在α=30°～45°左右取得最小值，且柱体横流向位移不再由升力决定.错列角度固定时，折合速度的改变对上下游柱体升阻力的影响都较小.

(5)下游圆柱在顺流向的振动响应受错列角度的影响较大.当0°＜α≤45°时，下游圆柱的顺流向振动都比上游圆柱剧烈；当45°＜α≤90°，上下游圆柱的顺流向振动相差较小.

可视化是三维模型设计的最大特点，也是相比传统设计的优势所在。以变送器为例，设计人员可以非常直观地看到变送器与取源点的位置关系，以及安装场地的周边情况，进而根据气上液下的原则调整变送器与取源点的相对位置关系，旋转、移动变送器到便于安装检修的朝向及位置，同时避开其他的设备、管道和结构。变送器建模示例如图1所示。

学术英语智能写作评估的新研究——评《基于语类的二语研究论文智能写作评估：从设计、评估到优化》………………………………………………………………………………………… 钟家宝 钟兰凤（4.108）

参考文献

1 Bearman PW.Vortex shedding from oscillating bluffbodies.Annual Review of Fluid Mechanics,2012,16(1)∶195-222

2 Nomurt T.Finite element analysis of vortex-induced vibrations of bluffcylinder.Journal of Wind Engineering and Industrial Aerodynamics,1993,46∶587-594

3 及春宁,花阳,许栋等.不同剪切率来流作用下柔性圆柱涡激振动数值模拟.力学学报,2018,50(1)∶21-31(Ji Chunning,Hua Yang,Xu Dong,et al.Numerical simulation of vortex-induced virbration of a flxible cylinder exposed to shear flw at different shear rates.ChineseJournalofTheoreticalandAppliedMechanics,2018,50(1)∶21-31(in Chinese))

4 高云,付世晓,熊友明等.剪切来流下柔性圆柱体涡激振动响应试验研究.振动与冲击,2016,35(20)∶142-148(Gao Yun,Fu Shixiao,Xiong Youming,et al.Experimental study on vortex induced vibration responses of flxible cylinder in sheared current.Journal of Vibration and Shock,2016,35(20)∶142-148(in Chinese))

5 宋磊建,付世晓,于大鹏等.剪切流下发生涡激振动的柔性立管阻力特性研究.力学学报,2016,48(2)∶300-306(Song Leijian,Fu Shixiao,Yu Dapeng,et al.Investigation of drag forces for flxible risers undergoing vortex-induced vibration in sheared flw.Journal of Vibration and Shock,2016,48(2)∶300-306(in Chinese))

6 Zhao M,Cheng L,An HW.Numerical investigation of vortexinduced virbration of circular cylinder in transverse direction in oscillatory flw.Ocean Engineering,2012,41(5)∶39-52

7 王俊高,付世晓,许玉旺等.正弦振荡来流下柔性立管涡激振动发展过程.力学学报,2014,46(2)∶173-182(Wang Jungao,Fu Shixiao,Xu Yuwang,et al.VIV developing process of a flxible cylinder under oscillatory flw.Chinese Journal of Theoretical and Applied Mechanics,2014，46(2)∶173-182(in Chinese))

8 Mittal S,Kumar V,Raghuvanshi A.Unsteady incompressible flw past two cylinders in tandem and staggered arrangements.International Journal of Numerical Methods in Fluids,1997,25∶1315-1344

9 Prasanth TK,Mittal S.Vortex-induced vibration of two circular cylinders at low Reynolds number.Journal of Fluid and Structures,2009,25(4)∶731-741

10 Borazjani I,Sotiropoulos F.Vortex-induced vibrations of two cylinders in tandem arrangement in the proximity-wake interference region.Journal of Fluid Mechanics,2009,621∶321-364

11 吴晓笛,刘华坪,陈浮.基于浸入边界--多松弛时间格子玻尔兹曼通量求解法的流固耦合算法研究.物理学报,2017,66(22)∶252-265(Wu Xiaodi,Liu Huaping,Chen Fu.A method combined immersed boundary with multi-relaxation-time lattice Boltzmann flu solver for fluid-structurinteraction.Acta Physica Sinica,2017,66(22)∶252-265(in Chinese))

12 郭晓玲,唐国强,刘名名等.低雷诺数下串联双圆柱涡激振动机理的数值研究.振动与冲击,2014,33(4)∶60-69(Guo Xiaoling,Tang Guoqiang,Liu Mingming,et al.Numerical investigation on vortex-induced vibration of twin tandem circular cylinders under low Reynolds number.Journal of Vibration and Shock,2014,33(4)∶60-69(in Chinese))

13 Williamson CHK.Evolution of a single wake behind a pair of bluffbodies.Journal of Fluid Mechanics,2006,159(159)∶1-18

14 Chen W,Ji C,Xu W,et al.On the responses and wake patterns of two side-by-side elastically circular cylinders in uniform laminar flw.Journal of Fluids and Structure,2015,55∶218-236

15 赵西增,付英男,张大可.柱体绕流的 CIP方法模拟.哈尔滨工程大学学报,2016,37(3)∶297-305(Zhao Xizeng,Fu Yingnan,Zhang Dake.Numerical simulation of flw past a cylinder using a CIP-Based model.Journal of Harbin Engineering University,2016,37(3)∶297-305(in Chinese))

16 Zhao X,Cheng D,Zhang D,et al. Numerical study of low-Reynolds-number flw past two tandem square cylinders with varying incident angles of the downstream one using a CIP-based model.Ocean Engineering,2016,121∶414-421

17 赵西增，付英男，张大可等.紧致插值曲线方法在流向强迫振荡圆柱绕流中的应用.力学学报,2015,47(3)∶441-450(Zhao Xizeng,Fu Yingnan,Zhang Dake,et al.Application of a CIP-based numerical simulation of flw past an in-line forced oscillating circular cylinder.Chinese Journal of Theoretical and Applied Mechanics,2015,47(3)∶441-450(in Chinese))

18 Yabe T,Xiao F,Utsumi T.The constrained interpolation profil method for multiphase analysis.Journal of Computational Physics,2001,169(2)∶556-593

19 Peskin CS.Flow patterns around heart valves∶a digital computer method for solving the equations of motion.[PhD Thesis].New York∶Albert Einstein College of Medicine,1972

20 Zhao M.Flow induced vibration of two rigidly coupled circular cylinders in tandem and side-by-side arrangements at a low Reynolds number of 150.Physics of Fluids,2013,25(12)∶355-381

21 Meneghini JR,Saltara F,Siqueira CLR,et al.Numerical simulation of flw interference between two circular cylinders in tandem and side-by-side arrangements.Journal of Fluids and Structures,2001,15(2)∶327-350

22 Zhou CY,So R,Lam K.Vortex-induced vibrations of an elastic circular cylinder.Journal of Fluids and Structures,1999,13(2)∶165-189

23 Ahn HT,Kallinderis Y.Strongly coupled flw/structure interactions with a geometrically conservative ALE scheme on general hybrid meshes.Journal of Computational Physics,2006,219(2)∶671-696

24 Bao Y,Zhou D,Tu J.Flow interference between a stationary cylinder and an elastically mounted cylinder arranged in proximity.Journal of Fluids and Structures,2011,27(8)∶1425-1446

25 及春宁,陈威霖,黄继露等.串列双圆柱流致振动的数值模拟及其耦合机制.力学学报,2014,46(6)∶862-870(Ji Chunning,Chen Weilin,Huang Jilu,et al.Numerical investigation on flw-induced vibration of two cylinders in tandem arrangements and its coupling mechanism.ChineseJournalofTheoreticalandAppliedMechanics,2014,46(6)∶862-870(in Chinese))

26 Prasanth TK,Mittal S.Vortex-induced vibration of two circular cylinders at low Reynolds number.Journal of Fluids and Structures,2009,25(4)∶731-741

27 Zravkovich MM.The effects of flw interference between two circular cylinders in various arrangements.Journal of Fluids and Structures,1987,1∶239-261

28 Chen HC,Huang K,Chen CR,et al.CFD simulation of riser VIV.Final Project Report,2006∶15-20

29 杨立红.低雷诺数下串列多圆柱流致振动的数值模拟及其机理研究.[博士论文].天津：天津大学,2014(Yang Lihong.Numerical investigation on flw-induced vibrations of multiple circular cylinders in tandem arrangement at low Reynolds number.[PhD Thesis].Tianjin∶Tianjin University,2014(in Chinese))