基于旋转运动的试验气流场产生方法数值研究 Numerical study on experimental flow field generation method based on rotating motion first-author corresp 蒋兴良，教授，从事输配电装备复杂环境灾害防御等研究，E-mail：xljiang@cqu.edu.cn 蒋兴良，教授，从事输配电装备复杂环境灾害防御等研究，E-mail：xljiang@cqu.edu.cn xljiang@cqu.edu.cn 重庆大学 输配电装备及系统安全与新技术国家重点实验室，重庆，400044 重庆大学 输配电装备及系统安全与新技术国家重点实验室，重庆，400044 State Key Laboratory of Power Transmission Equipment & System Security and New Technology, Chongqing University, Chongqing 400044, China State Key Laboratory of Power Transmission Equipment & System Security and New Technology, Chongqing University, Chongqing 400044, China 赵俊 为了丰富气流场试验手段以及满足复杂环境气流场试验的要求，提出一种基于旋转运动的试验气流场产生方法。建立基于旋转参考系以及滑移网格的数值求解模型并对有限空间内的三维旋转流场分布特性进行研究。研究结果表明：涵道结构可以有效提升试验段内的平均流速，其加速效果随着转速和旋转半径的增加而增加。试验段内动压系数及气流偏角则随着转速增加而增加，随着旋转半径的增加而减小。通过增加旋转半径来提升试验流速和流场品质是最合理的方式。涵道内的流场的径向偏移是影响试验段内的试验流场品质的重要原因，后续研究需对其进行校正以增加基于旋转运动的试验气流场产生方法的适用性。 In order to enrich the means of airflow field test and meet the requirements of airflow field test in complex environment, a method of generating experimental airflow field based on rotating motion was proposed. A numerical solution model based on rotating reference frame and slipping mesh was established to study the distribution characteristics of three-dimensional rotating flow field in test space. The results show that the shrinkage duct can effectively improve the average velocity in the test section, and the acceleration effect increases with the increase of rotational speed and radius. Meanwhile, the dynamic pressure coefficient and airflow deflection angle in the test section increase with the increase of rotation speed, but decrease with the increase of rotation radius. Increasing the rotation radius is the most reasonable way to improve the test velocity and flow field quality. The radial direction offset of the flow field in the shrinkage duct is an important factor affecting the quality of the experimental flow field in the test section, which should be overcome in the follow-up study to increase the applicability of this method. 流场试验 数值模拟 高速运动 旋转运动 flow field test numerical simulation high-speed motion rotary motion 气流场试验是空气动力学研究的重要手段，是航空航天、能源、交通运输、环境以及建筑等领域的重要技术 [ 1 1 - 2 2 1-2 ] 。空气动力学来源于流体力学，伴随着航空航天工业的发展而迅速发展，自20世纪60年代以来，随着建筑、能源、运输以及交通等工业的发展，空气动力学朝着多目标多场景发展。目前，空气动力学的主要研究手段包括理论分析、数值模拟以及气流场试验，三大研究手段在空气动力学的发展中都发挥了重大作用 [ 3 3 - 5 5 3-5 ] 。尤其是以风洞试验为代表的气流场试验，不仅为飞行器等设备提供可靠的设计依据，促进了航空航天工业的发展，而且有力地促进了空气动力学的发展 [ 6 6 - 8 8 6-8 ] 。迄今为止，气流场试验仍是研究空气动力学特性的重要手段。 现有的不同形式的风洞主要用于飞行器整体、机翼等设备的气动特性的试验研究 [ 9 9 - 11 11 9-11 ] 。首先，高精度风洞设备的设计与建设成本高昂，占地面积及运行功率极大。其次，对复杂大气环境工况下的模拟试验，如高速动车组牵引机车车顶高压设备在大雾、沙尘环境下的高速运行，风力发电机叶片在高海拔、冰雪环境下的高速旋转等，难以在传统风洞中进行模拟试验 [ 12 12 - 13 13 12-13 ] 。因此，亟需一种新的试验气流场产生方法，以较低成本实现设备在复杂环境高速气流中运行性能的检测，推动气流场试验技术的发展。 由于旋转流场现场实验以及测量手段匮乏，数值模拟仍是研究的主要手段 [ 14 14 - 16 16 14-16 ] ，而旋转参考系方法和搭接网格方法在旋转流场的数值求解中发挥了重要作用，被广泛应用于泵类机械 [ 17 17 - 18 18 17-18 ] 、直升机旋翼 [ 19 19 - 20 20 19-20 ] 、风力发电机 [ 21 21 - 22 22 21-22 ] 、生物搅拌器 [ 23 23 ] 以及车轮制动盘 [ 24 24 ] 等研究。 因此，根据相对运动原理并参考风洞和传统的旋转臂设备，本文提出一种基于旋转运动的试验气流场产生方法。该方法流速的产生方式与传统风洞的不同，是在有限试验空间中通过圆周驱动的方法将加速涵道结构加速至高速，在涵道试验段内形成高速试验气流，以低成本的方式实现高速气流及复杂大气环境甚至高压电场的综合气流场试验条件。建立基于旋转参考系以及滑移的数值求解模型对试验空间内的三维旋转流场进行模拟计算，评估试验段内的流场分布特性，从流速和流场均匀性的角度验证基于旋转运动的试验气流场产生方法的可行性。 1 基于旋转的均质气流产生方法 1 图1 图1 所示为基于旋转运动的试验气流场产生方法的示意图。试验室内部形成密闭试验空间，高速旋转驱动系统放置在试验室内。试验室采用密封绝热设计，以满足模拟低温低压试验环境的要求，可以模拟各类复杂大气环境。 高速旋转驱动系统主要由自平衡导轨与滑台系统、高速驱动与传动系统组成。自平衡导轨与滑台系统是高速运动平台的核心，主要起支承均流加速涵道高速运转的作用。导轨系统由上下2个大尺寸圆形导轨和安装在其上的2个滑台组成，导轨与滑台系统配合限定旋转部件只沿导轨做高速圆周运动。为了保证系统高速旋转时的稳定性，流场控制涵道对称侧安装有补偿质量自动分布型自动平衡装置，以抵消旋转系统本身的不平衡。均流加速涵道由收缩段、试验段以及扩散段组成，通过滑台固定在导轨之间，通过旋转主轴和旋转支撑杆的带动在导轨上做高速圆周运动。 理论上，基于旋转运动的试验气流场产生方法可以在涵道试验段内产生高速气流，并且可以实现高速气流与复杂大气环境的综合气流场试验条件。本文将从试验段内流速和流场品质的角度对其流场分布进行深入分析。 2 有限空间三维旋转流场数值计算模型 1 为了方便对比，将有限空间内的三维旋转流动当作可压缩湍流流动处理。本文采用多参考系计算模型结合滑移网格方法对其进行求解。将流体域划分为外场和内场2部分。内场计算域模型采用旋转参考系，将时变流场求解转变为非时变流场进行求解；而外场计算域模型采用常规静止参考系，通过插值的方法进行内外场计算域数据传递。 2.1 　多区域计算域建模 2 流场计算域如 图2 图2 所示，其中，内场计算域主要为旋转的加速涵道，外场计算域主要为密闭空间的静止壁面， r 为涵道的旋转半径， R 1 为内场计算域外径， R 2 为外场计算域外径， H 1 为内场计算域高度， H 2 为外场计算域高度， D 1 为涵道收缩段入口半径， D 2 为涵道试验段半径， D 3 为涵道扩散段出口半径， L 1 为涵道收缩段长度， L 2 为涵道试验段长度， L 3 为涵道扩散段长度。将外场计算域分为3个部分，分别为上下的圆柱体以及周围的圆环体。计算域以及涵道的具体参数如 表1 表1 所示。 r /m R 1 /m R 2 /m H 1 /m H 2 /m D 1 /m D 2 /m D 3 /m L 1 /m L 2 /m L 3 /m 6 8 10 2 4 0.5 0.25 0.4 0.9 0.5 0.5 2.2 　局部控制混合网格划分策略 2 通过分区划分、混合处理、局部控制等方法形成多区域混合网格。其中，外场计算域分为3个部分后，采用扫略方法生成六面体结构化网格；内场网格由于包含涵道结构，是不规则的形状，则生成以四面体为主的非结构化网格。 为了确保计算精度，采取局部控制的方法对涵道周围网格进行局部加密处理，涵道周围球形区域内的网格尺寸约束为整体网格控制尺寸的40%。其中，在试验段最主要的流动观察区域，将试验段内的网格尺寸约束为整体网格控制尺寸的15%。根据本文使用的 k - ε 模型以及增强型壁面函数的要求，根据壁面距离(量纲一的量) y + = 30确定近壁区的边界层网格的首层高度；边界层的长宽(即涵道结构的面网格尺寸)则按照边界层网格首层高度的100倍确定，然后按照扩散比1.2生成5层边界层网格。数值计算网格如 图3 图3 所示。 为了确保计算的准确性，分别采用RNG模型和Realizable模型对数值网格的独立性进行研究，得到涵道内平均流速随整体网格控制尺寸的变化情况如 图4 图4 所示。从 图4 图4 可以看出：在当整体网格控制尺寸从0.20 mm降到0.14 mm后，网格控制尺寸对涵道内流场分布的影响可以忽略。后续研究中均采用上述的网格划分策略和网格控制尺寸进行数值计算。 2.3 　外流计算域的流动控制方程 2 在静止参考系下求解静态的外流计算域流动控制方程。在连续介质的前提下，计算域外流场流动应该满足质量守恒方程、动量守恒方程，能量守恒方程。其中，质量守恒方程为 <math display="block" id="M1"><mfrac><mrow><mo>∂</mo><mi>ρ</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mo>∇</mo><mo>·</mo><mfenced separators="|"><mrow><mi>ρ</mi><mi mathvariant="bold-italic">v</mi></mrow></mfenced><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705776&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705750&type= 25.14599991 7.28133297 式中： t 为时间； ρ 为流体密度； v 为速度矢量。 动量守恒方程为 <math display="block" id="M2"><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo stretchy="false">(</mo><mi>ρ</mi><mi mathvariant="bold-italic">v</mi><mo stretchy="false">)</mo><mo>+</mo><mo>∇</mo><mo>·</mo><mfenced separators="|"><mrow><mi>ρ</mi><mi mathvariant="bold-italic">v</mi><mi mathvariant="bold-italic">v</mi></mrow></mfenced><mo>=</mo><mo>-</mo><mo>∇</mo><mi>p</mi><mo>+</mo><mo>∇</mo><msup><mrow><mi mathvariant="bold-italic">τ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msup><mo>+</mo><mi>ρ</mi><mi mathvariant="bold-italic">g</mi><mo>+</mo><mi mathvariant="bold-italic">F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705752&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705763&type= 58.50466156 7.28133297 式中： <math id="M3"><msup><mrow><mi mathvariant="bold-italic">τ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705780&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705778&type= 3.13266683 2.96333337 为黏性应力张量； g 为重力加速度； F 为外部体积力，如多相流相互作用力等。 能量方程为 <math display="block" id="M4"><mfrac><mrow><mo>∂</mo><mi>ρ</mi><mi>E</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mo>∇</mo><mfenced separators="|"><mrow><mi>ρ</mi><mi mathvariant="bold-italic">v</mi><mi>H</mi></mrow></mfenced><mo>=</mo><mo>∇</mo><mo>⋅</mo><mfenced separators="|"><mrow><mi>κ</mi><mfenced separators="|"><mrow><mo>∇</mo><mi>T</mi></mrow></mfenced><mo>+</mo><msub><mrow><mi>v</mi></mrow><mrow><mi>i</mi></mrow></msub><msup><mrow><mi>τ</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msup></mrow></mfenced><mo>+</mo><mi>ρ</mi><mi mathvariant="bold-italic">g</mi><mi mathvariant="bold-italic">v</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705766&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705765&type= 65.61666107 7.28133297 式中： H 为流体总熵； T 为流体温度； κ 为流体的热导率； E 为流体的物质能量。 2.4 　内流场计算域的流动控制方程 2 在旋转参考系下求解内流计算域流动控制方程。静止参考系和旋转参考系的速度变换方程为 <math display="block" id="M5"><mi mathvariant="bold-italic">v</mi><mo>=</mo><mi>ω</mi><mo>×</mo><mi mathvariant="bold-italic">r</mi><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705784&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705783&type= 18.11866570 3.72533321 式中： v r 为旋转参考系下的速度矢量。 加速度变换方程为 <math display="block" id="M6"><mi mathvariant="bold-italic">a</mi><mo>=</mo><mfrac><mrow><mi mathvariant="normal">d</mi><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="normal">d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mn mathvariant="normal">2</mn><mi>ω</mi><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub><mo>+</mo><mi>ω</mi><mo>×</mo><mo stretchy="false">(</mo><mi>ω</mi><mo>×</mo><mi mathvariant="bold-italic">r</mi><mo stretchy="false">)</mo><mo>+</mo><mi>α</mi><mo>×</mo><mi mathvariant="bold-italic">r</mi><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705768&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705754&type= 58.75866318 7.28133297 式中： ω 为旋转参考系的旋转角度， 为运动的角加速度； a 为静止参考系下的加速度； a r 为旋转参考系下的加速度。 将 式(4) 式(4) 及 式(5) 式(5) 代入 式(1) 式(1) ，得到旋转参考系下质量守恒方程为 <math display="block" id="M7"><mfrac><mrow><mo>∂</mo><mi>ρ</mi></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mo>∇</mo><mfenced separators="|"><mrow><mi>ρ</mi><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub></mrow></mfenced><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705786&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705769&type= 25.14599991 7.28133297 将 式(4) 式(4) 以及 式(5) 式(5) 代入 式(2) 式(2) ，得到旋转坐标系下的动量守恒方程为 <math display="block" id="M8"><mfrac><mrow><mo>∂</mo><mi>ρ</mi><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>t</mi></mrow></mfrac><mo>+</mo><mo>∇</mo><mfenced separators="|"><mrow><mi>ρ</mi><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub></mrow></mfenced><mo>+</mo><mi>ρ</mi><mo stretchy="false">(</mo><mn mathvariant="normal">2</mn><mi>ω</mi><mo>×</mo><msub><mrow><mi mathvariant="bold-italic">v</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub><mo>+</mo><mi>ω</mi><mo>×</mo><mo stretchy="false">(</mo><mi>ω</mi><mo>×</mo><mi mathvariant="bold-italic">r</mi><mo stretchy="false">)</mo><mo>+</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705771&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705756&type= 62.56866837 7.36600018 <math display="block" id="M9"><mi>α</mi><mo>×</mo><mi mathvariant="bold-italic">r</mi><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">a</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mo>-</mo><mo>∇</mo><mi>p</mi><mo>+</mo><mo>∇</mo><msubsup><mrow><mi>τ</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msubsup><mo>+</mo><mi>ρ</mi><mi mathvariant="bold-italic">g</mi><mo>+</mo><mi mathvariant="bold-italic">F</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705803&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705757&type= 46.05867004 4.82600021 式中： <math id="M10"><msubsup><mrow><mi>τ</mi></mrow><mrow><mi mathvariant="normal">r</mi></mrow><mrow><mi>i</mi><mi>j</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705818&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705789&type= 3.04800010 4.06400013 为运动参考系下的黏性应力张量。 转换坐标系后，能量方程的形式不变。 2.5 　模型设置与求解工况 2 在建立静止参考系以及旋转参考系下的流动控制方程后，采用有限体积法对密闭空间内的三维旋转流场进行求解。计算过程中发现，Coupled算法相较于分离算法具有更好的收敛性，因此，采用Coupled算法进行迭代计算。为了保证计算的精确性，迭代计算的时间步长为涵道旋转周期的1/1 000，即当涵道转速为60 r/min时，对应迭代时间步长为0.001 s。具体的模型求解设置如 表2 表2 所示。 湍流模型 亚松弛 因子 迭代时间步长/s 离散格式 迭代算法 壁面函数 压力插值方法 梯度插值方法 RNG k - ε 模型 0.2~0.5 0.001 二阶迎风 Coupled算法 增强型壁 面函数 PRESTO!方法 基于网格中心的最小二乘法 3 密闭空间旋转流场分布特性 1 3.1 　空间流场的径向分布特性 2 对密闭空间内的旋转流场进行求解，得到某时刻计算域 XOY 截面速度分布如 图5 图5 所示。从 图5 图5 可以看出：涵道的旋转运动将带动密闭空间内的空气随之运动形成旋转尾流，旋转尾流大致沿旋转运动的轨迹分布，并且逐渐向径向外侧扩散。大体上，尾流的速度在涵道周围以及后部较大，然后沿涵道的运动轨迹逐渐减小。 为了进一步分析空间尾流的分布情况，在 XOY 截面内以旋转中心为起点，建立数条放射性的径向监测线( 图5 图5 )，其中，每条相邻的监测线之间相差45°，规定每条监测线的方向均由旋转中心指向外侧。沿不同角度径向监测线的速度分布如 图6 图6 所示。从 图6 图6 可以看出，沿各监测线，速度大致呈现出先增加后减小的趋势。在靠近圆心的区域，速度很小；在旋转尾迹部分，速度达到最大；在外壁面处，速度减小为0 m/s。由于离涵道太近，沿45°监测线，速度分布顶部呈现双峰形式，除此之外其他角度的沿线速度都呈现出单峰的形式。 3.2 　空间流场的轴向分布特性 2 计算域 XOZ 截面的速度分布如 图7 图7 所示。从 图7 图7 可以看到：在涵道旋转轨迹附近，密闭空间尾流的速度最大，在朝四周扩散的过程中已经到达了径向和轴向的壁面处；旋转径向的尾流扩散效果要比旋转轴向的扩散效果好。 为了进一步研究轴向的尾流分布情况，根据径向监测线建立轴向监测线。由于计算域上下的流场对称分布，因此，轴向监测线只位于上半计算域。各轴向监测线起始于 XOY 中心截面内对应径向监测线上尾流最大值出现的位置，一直延伸至轴向边界处，其位置如 图5 图5 中“×”号所示。沿不同角度轴向监测线的速度分布情况如 图8 图8 所示。从 图8 图8 可以看出，不同角度监测线速度分布呈现相似的趋势。当轴向坐标较小时，尾流速度较大，然后随着轴向坐标的增大而逐渐减小。同时，当轴向位置较小时，不同角度的速度分布差异较大，但是随着轴向位置的增加，不同角度的尾流速度的差异逐渐变小。当轴向坐标超过2.0 m之后，各角度的速度线基本上重合。 综上所述，旋转带动的空间尾流将对流场试验产生影响，主要体现在空间尾流的尾迹会和涵道开口相接(转速很小、旋转半径很大时则不会相接)，虽然此时空间尾流的速度已经衰减至较小，但是还是会减小涵道内的相对试验速度，并且转速越高、旋转半径越小，这种减速效果会越明显。 3.3 　试验段内的流场分布 2 涵道内不同截面的速度分布如 图9 图9 所示。从 图9 图9 可以看出：涵道结构具有明显加速效果，流体经过收缩段的收缩加速在试验段内形成高速气流，然后再通过扩散段降速流出涵道。此外，由于旋转时内外圈相对速度的差异导致试验段内 XOY 平面(旋转径向)的速度分布出现偏移，该现象由旋转运动的特性导致，而在 XOZ 平面内并没有类似情况的出现。 (a) XOZ 平面；(b) XOY 平面 在涵道试验段中心处建立检测线1以及检测线2对试验段内的速度分布进行进一步分析。沿监测线1和检测线2的速度分布如 图10 图10 所示。从 图10 图10 可以看出，除近壁面速度迅速减小的区域外，沿监测线1速度则呈现出平稳的趋势，而沿监测线2的速度随着坐标的增大而逐渐增加；模型区内沿监测线2速度的最大值与试验段平均气流速度的相对误差为7.25%，速度最小值与试验段平均气流速度的相对误差为4.97%。 4 试验段内的流速和流场品质特性 1 试验段为气流场试验的核心区域，试验段内的流速和流场品质将直接影响气流场试验效果以及所获数据的准确性。本节将对涵道试验段内的平均流速和流场品质进行分析，试验段内的平均流速以及其他流场品质指标值均为试验段内1 000个均匀取点相应参数的算术平均值。 4.1 　收缩比对试验段流速的影响 2 收缩段进出口面积之比称为收缩比，收缩比与涵道加速效果直接相关。保持旋转半径为6 m，得到不同收缩比试验段内的平均流速的变化如 图11 图11 所示， 图11 图11 中，运动线速度即为涵道中心处的旋转线速度。从 图11 图11 可以看出，涵道结构可以明显提高试验段内的平均流速，当转速相同时，涵道试验段的平均流速随着收缩比的增大而增大。 定义涵道试验段平均流速相对于涵道中心处旋转线速度的提高比例为涵道的加速效果 k 1 ，即 <math display="block" id="M11"><msub><mrow><mi>k</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><msub><mrow><mi>v</mi></mrow><mrow><mi mathvariant="normal">a</mi></mrow></msub><mo>/</mo><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><mi>r</mi><mo>·</mo><mfrac><mrow><mi>n</mi></mrow><mrow><mn mathvariant="normal">60</mn></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705832&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705834&type= 21.59000015 7.28133297 式中： v a 为试验段平均流速； r 为涵道旋转半径； n 为涵道转速，r/min。 保持旋转半径为6 m，计算得到不同转速下 k 1 与收缩比的关系，如 图12 图12 所示。从 图12 图12 可以看出，随着收缩比的增加， k 1 逐渐增加，但是该增加趋势呈现出饱和效应。这是因为当收缩比过大时，收缩段壁面的阻滞效应也会增大，这将会影响涵道内气流的各向同性，甚至产生气壁分离，导致加速效果下降，并且影响试验流场品质 [ 25 25 ] 。因此，不能无限制的增加收缩比来增加试验段内的流速。从 图12 图12 还可以看出，当收缩比较小、转速也较小时， k 1 也较小，说明此时并没有完全发挥出涵道的加速性能，可以适当增加涵道转速以获得更合理的加速效果。采用加速涵道结构后，试验段内平均气流速度相较于旋转线速度都有了明显提升。可见采取涵道加速结构之后，可以以较低的转速和较小的旋转半径获得较大的气流速度，极大地提高了基于旋转运动气流场试验方法的试验效率和可行性。 4.2 　旋转半径和转速对试验段流速的影响 2 保持旋转半径为6 m， k 1 与转速的关系如 图 13 图 13 (a)所示。从 图13 图13 (a)可以看出，随着转速的增加， k 1 呈现先增大后减小的趋势。当收缩比为4时， k 1 在转速60 r/min出现最大值；当收缩比为9和16时， k 1 在90 r/min出现极大值，随后均随着转速的增加而逐渐减小。该趋势是旋转尾流的影响所致；当转速较小时，旋转尾流也较小，不会对涵道入口的相对气流速度造成影响，随着转速的增加，旋转尾流逐渐增大，造成涵道入口处相对气流速度的减小，进而导致试验段内的平均速度减小，削弱涵道加速效果。 (a) 转速的影响；(b) 旋转半径的影响 保持转速为60 r/min，旋转半径对涵道加速效果的影响如 图13 图13 (b)所示。从 图13 图13 (b)可以看出：随着旋转的增加， k 1 呈现单调增加的趋势。旋转半径的增加会减弱旋转尾流对相对试验速度的影响，因此，转速越大涵道加速效果越明显。同时，涵道旋转半径的增加还可以减小旋转运动所导致的速度偏移，因此，从原理上来说，基于旋转运动的气流场产生方法中通过增加旋转半径的来提升试验流速要比增加转速更合理。 4.3 　试验段流场品质评价指标 2 基于旋转运动的气流场产生方法的最大误差源自旋转运动带来的流场的不均匀性。在气流场试验中，与流场的均匀性评价有关的指标主要是动压系数和气流偏角，因此，本研究选择对试验段内的平均动压系数以及平均气流偏角进行计算评估 [ 26 26 - 27 27 26-27 ] 。 动压系数反映的是流速的均匀性，本文将流场中某点的动压系数定义为该点的动压与模型区平均动压的偏差。 <math display="block" id="M12"><msub><mrow><mi>μ</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mfenced open="|" close="|" separators="|"><mrow><mfrac><mrow><msub><mrow><mi>q</mi></mrow><mrow><mi>i</mi></mrow></msub></mrow><mrow><mi>q</mi></mrow></mfrac><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705841&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705851&type= 17.94933319 9.22866726 式中： q i 为某点的动压； q 为模型区的平均动压。 <math display="block" id="M13"><msub><mrow><mi>q</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mi>ρ</mi><msup><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705853&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705867&type= 14.98600006 7.19666624 气流偏角反映的是方向场的均匀性，定义为该位置处的气流的方向偏移主流方向的角度，即： <math display="block" id="M14"><msub><mrow><mi>β</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>=</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfenced separators="|"><mrow><mfrac><mrow><mroot><mrow><msup><mrow><msub><mrow><mi>v</mi></mrow><mrow><mi>y</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><msub><mrow><mi>v</mi></mrow><mrow><mi>z</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot></mrow><mrow><msub><mrow><mi>v</mi></mrow><mrow><mi>x</mi></mrow></msub></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705854&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44705843&type= 35.56000137 13.71599960 式中： β i 为某点的气流偏角， v x 、 v y 和 v z 分别为流速沿 x 方向(主流方向)、 y 方向和 z 方向的分量。 理想情况下，动压系数越小越好，其越接近于0，说明试验段内速度场越均匀；试验段内的平均动压系数也是越小越好，其越接近于0，说明试验段内气流方向越均匀。 4.4 　转速对试验段流场品质的影响 2 转速对试验段动压系数以及气流偏角的影响规律如 图14 图14 所示。从 图14 图14 可以看出：不同旋转半径时，动压系数和气流偏角随着转速的增加而增加，说明随着转速的增加，试验段内的试验流场品质逐渐降低。当转速较小时，动压系数和气流偏角随转速的增加较为缓慢，说明此时转速的增加对试验流场品质的影响并不显著；反之，转速的增加将严重影响试验段内的试验流场品质。一方面，转速的增加加剧了旋转径向相对试验速度的差异，导致试验段内的流场更加不均匀；另一方面，转速的增加也加剧了密闭空间内旋转尾流，进一步影响了试验段内的流场品质，并且随着转速的增加，这种影响越来越明显。因此，为了提升试验流场品质，应该尽量减小转速。 (a) 对动压系数的影响；(b) 对气流偏角的影响 4.5 　旋转半径对试验段流场品质的影响 2 旋转半径对试验段动压系数及气流偏角的影响如 图15 图15 所示。从 图15 图15 可以看出：不同转速时，试验段内的动压系数和气流偏角均随着旋转半径的增加而减小，说明旋转半径的增加将有效提升试验段内的试验流场品质。一方面，旋转半径的增加减少了旋转径向相对试验速度不均匀的情况；另一方面，旋转半径的增加也减弱的旋转尾流的影响。因此，在成本允许的前提下，旋转半径越大越好，在极限情况下，旋转半径无限大时，试验流场不均匀性的问题也不复存在了。 (a) 对动压系数的影响；(b) 对气流偏角的影响 以上对于转速和旋转半径的讨论均忽略了其对流速的影响。试验段内的平均流速与转速和半径成正比，转速太小则涵道内试验流速达不到试验要求，转速太大则涵道内的试验流场品质下降。而旋转半径的增大对于增加试验流速和提升试验流场品质均有利，但是会急剧增加试验装置体积和造价。因此，在试验成本允许的条件下，通过增加旋转半径来提升试验流速和试验流场品质是最合理的方式。 5 结论 1 1) 有限空间内的旋转尾流大致沿旋转运动的轨迹分布，其速度逐渐减小的同时也在朝四周扩散，旋转径向的尾流扩散效果要明显比旋转轴向的好。 2) 涵道结构可以显著提升试验段内的平均流速。涵道的加速效果随着涵道收缩比以及旋转半径的增加而增加，随着转速的增加先增加后减小。 3) 试验段内的平均动压系数以及气流偏角随着转速的增加而增加，随着旋转半径的增加而减小。通过增加旋转半径来提升试验流速和试验流场品质是最合理的方式。 4) 在基于旋转运动的试验气流场产生方法中，旋转运动线速度的差异将使涵道内的流场分布出现径向偏移，严重影响试验段内的试验流场品质。后续研究中需要探索相应校正方法以进一步提升试验内的流场均匀性，使其符合更高精度流场试验的要求。 参考文献 聂万盛 , 车科学 , 何浩波 , 等 . 空气动力学与试验技术 [M]. 北京 : 国防工业出版社 , 2018 : 1 - 4 . 空气动力学与试验技术 NIE Wansheng , CHE Kexue , HE Haobo , et al . Aerodynamics and test techniques [M]. Beijing : National Defense Industry Press , 2018 : 1 - 4 . Aerodynamics and test techniques 王钢林 , 郑遂 . 从战斗机的发展历程看空气动力学的贡献 [J]. 物理 , 2016 , 45 ( 10 ): 645 - 650 . 从战斗机的发展历程看空气动力学的贡献 WANG Ganglin , ZHENG Sui . The development of fighter planes and the contribution of aerodynamics [J]. Physics , 2016 , 45 ( 10 ): 645 - 650 . The development of fighter planes and the contribution of aerodynamics 李正农 , 余世斌 , 吴红华 , 等 . 强风作用下树木周围流场的数值模拟研究 [J]. 中南大学学报(自然科学版) , 2021 , 52 ( 11 ): 3970 - 3980 . 强风作用下树木周围流场的数值模拟研究 LI Zhengnong , YU Shibin , WU Honghua , et al . Numerical simulation of flow field around the tree in strong wind [J]. Journal of Central South University(Science and Technology) , 2021 , 52 ( 11 ): 3970 - 3980 . Numerical simulation of flow field around the tree in strong wind 刘锦生 , 徐胜金 , 王庆洋 , 等 . 基于物理模型的汽车空气动力学研究综述 [J]. 实验流体力学 , 2020 , 34 ( 1 ): 38 - 48 . 基于物理模型的汽车空气动力学研究综述 LIU Jinsheng , XU Shengjin , WANG Qingyang , et al . Review of automotive aerodynamics research based on physical models [J]. Journal of Experiments in Fluid Mechanics , 2020 , 34 ( 1 ): 38 - 48 . Review of automotive aerodynamics research based on physical models 张静 , 刘志刚 , 鲁小兵 , 等 . 高速弓网空气动力学研究进展 [J]. 铁道学报 , 2015 , 37 ( 1 ): 7 - 15 . 高速弓网空气动力学研究进展 ZHANG Jing , LIU Zhigang , LU Xiaobing , et al . Study on aerodynamics development of high-speed pantograph and catenary [J]. Journal of the China Railway Society , 2015 , 37 ( 1 ): 7 - 15 . Study on aerodynamics development of high-speed pantograph and catenary 井国庆 , 杜文博 , 丁东 , 等 . 基于风洞试验和CFD仿真的350km/h高速铁路飞砟研究 [J]. 中南大学学报(自然科学版) , 2020 , 51 ( 12 ): 3546 - 3553 . 基于风洞试验和CFD仿真的350km/h高速铁路飞砟研究 JING Guoqing , DU Wenbo , DING Dong , et al . 350 km/h high speed railway ballast flight analysis based on wind tunnel test and CFD simulation [J]. Journal of Central South University(Science and Technology) , 2020 , 51 ( 12 ): 3546 - 3553 . 350 km/h high speed railway ballast flight analysis based on wind tunnel test and CFD simulation ZHANG Jianhua , JIANG Shoushan . Rigid-flexible coupling model and dynamic analysis of rocket sled [J]. Advanced Materials Research , 2011 , 346 ( 9 ): 447 - 454 . Rigid-flexible coupling model and dynamic analysis of rocket sled GILI P A , BATTIPEDE M . Experimental validation of the wing dihedral effect using a whirling arm equipment [J]. Journal of Aircraft , 2001 , 38 ( 6 ): 1069 - 1075 . Experimental validation of the wing dihedral effect using a whirling arm equipment 罗月培 , 孙宗祥 , 孙杭义 , 等 . 美国高超声速风洞试验能力发展综述 [J]. 飞航导弹 , 2021 , 51 ( 6 ): 33 - 41 . 美国高超声速风洞试验能力发展综述 LUO Yuepei , SUN Zongxiang , SUN Hangyi , et al . Summary of the development of hypersonic wind tunnel test capability in the United States [J]. Aeroclynamic Missiles Journal , 2021 , 51 ( 6 ): 33 - 41 . Summary of the development of hypersonic wind tunnel test capability in the United States 战培国 , 曾慧 , 钟萍 , 等 . 2017年国外风洞试验发展动态综述 [J]. 飞航导弹 , 2018 , 48 ( 8 ): 16 - 20 . 2017年国外风洞试验发展动态综述 ZHAN Peiguo , ZENG Hui , ZHONG Ping , et al . Review on the development of foreign wind tunnel tests in 2017 [J]. Aeroclynamic Missiles Journal , 2018 , 48 ( 8 ): 16 - 20 . Review on the development of foreign wind tunnel tests in 2017 战培国 . 国外先进风洞测试技术综述 [J]. 飞航导弹 , 2017 , 47 ( 11 ): 75 - 79 . 国外先进风洞测试技术综述 ZHAN Peiguo . Summary of foreign advanced wind tunnel testing technology [J]. Aeroclynamic Missiles Journal , 2017 , 47 ( 11 ): 75 - 79 . Summary of foreign advanced wind tunnel testing technology 吴广宁 , 石超群 , 张血琴 . 高速列车车顶绝缘子的研究进展 [J]. 南方电网技术 , 2016 , 10 ( 3 ): 71 - 77 . 高速列车车顶绝缘子的研究进展 WU Guangning , SHI Chaoqun , ZHANG Xueqin . Research progress of roof insulator used for high-speed train [J]. Southern Power System Technology , 2016 , 10 ( 3 ): 71 - 77 . Research progress of roof insulator used for high-speed train SHU Lichun , LIANG Jian , HU Qin , et al . Study on small wind turbine icing and its performance [J]. Cold Regions Science and Technology , 2017 , 134 : 11 - 19 . Study on small wind turbine icing and its performance 王江云 , 毛羽 , 孟文 , 等 . 旋风分离器内非轴对称旋转流场的测量 [J]. 石油学报(石油加工) , 2015 , 31 ( 4 ): 920 - 929 . 旋风分离器内非轴对称旋转流场的测量 WANG Jiangyun , MAO Yu , MENG Wen , et al . Experimental measurement of non-axisymmetric rotating flow field in cyclone separator [J]. Acta Petrolei Sinica(Petroleum Processing Section) , 2015 , 31 ( 4 ): 920 - 929 . Experimental measurement of non-axisymmetric rotating flow field in cyclone separator 张涛 . 基于旋转坐标系的高精度流场数值模拟研究 [D]. 南京 : 南京航空航天大学 , 2018 : 8 - 16 . 基于旋转坐标系的高精度流场数值模拟研究 ZHANG Tao . Study on high-precision numerical simulation of flow field based on rotating coordinate system [D]. Nanjing : Nanjing University of Aeronautics and Astronautics , 2018 : 8 - 16 . Study on high-precision numerical simulation of flow field based on rotating coordinate system 王雷 . 旋转坐标系下旋转圆柱绕流的数值模拟 [J]. 低温建筑技术 , 2014 , 36 ( 2 ): 77 - 79 . 旋转坐标系下旋转圆柱绕流的数值模拟 WANG Lei . Numerical simulation of rotating cylinders streaming in rotating coordinates [J]. Low Temperature Architecture Technology , 2014 , 36 ( 2 ): 77 - 79 . Numerical simulation of rotating cylinders streaming in rotating coordinates CHEW J W , HILLS N J . Computational fluid dynamics for turbomachinery internal air systems [J]. Philosophical Transactions Series A, Mathematical, Physical, and Engineering Sciences , 2007 , 365 ( 5 ): 2587 - 2611 . Computational fluid dynamics for turbomachinery internal air systems DAWES W N . Turbomachinery computational fluid dynamics: asymptotes and paradigm shifts [J]. Philosophical Transactions Series A: Mathematical, Physical, and Engineering Sciences , 2007 , 365 ( 10 ): 2553 - 2585 . Turbomachinery computational fluid dynamics: asymptotes and paradigm shifts STRAWN R C , CARADONNA F X , DUQUE E P N . 30 years of rotorcraft computational fluid dynamics research and development [J]. Journal of the American Helicopter Society , 2006 , 51 ( 1 ): 5 - 21 . 30 years of rotorcraft computational fluid dynamics research and development 徐国华 , 招启军 . 直升机旋翼计算流体力学的研究进展 [J]. 南京航空航天大学学报 , 2003 , 35 ( 3 ): 338 - 344 . 直升机旋翼计算流体力学的研究进展 XU Guohua , ZHAO Qijun . Advances in computational fluid dynamics of helicopter rotor [J]. Journal of Nanjing University of Aeronautics & Astronautics , 2003 , 35 ( 3 ): 338 - 344 . Advances in computational fluid dynamics of helicopter rotor 张立军 , 朱怀宝 , 顾嘉伟 , 等 . 涡流发生器对垂直轴风力机翼型气动性能的影响 [J]. 中南大学学报(自然科学版) , 2020 , 51 ( 2 ): 540 - 550 . 涡流发生器对垂直轴风力机翼型气动性能的影响 ZHANG Lijun , ZHU Huaibao , GU Jiawei , et al . Influence of vortex generator on aerodynamic performance of airfoil of the vertical axis wind turbine [J]. Journal of Central South University(Science and Technology) , 2020 , 51 ( 2 ): 540 - 550 . Influence of vortex generator on aerodynamic performance of airfoil of the vertical axis wind turbine 张立军 , 顾嘉伟 , 朱怀宝 , 等 . 垂直轴风力机尾缘开裂襟翼气动性能及其偏转角调节规律 [J]. 中南大学学报(自然科学版) , 2020 , 51 ( 4 ): 931 - 943 . 垂直轴风力机尾缘开裂襟翼气动性能及其偏转角调节规律 ZHANG Lijun , GU Jiawei , ZHU Huaibao , et al . Aerodynamic performance on the trailing edge split flap and regulation law of deflection angle of VAWT [J]. Journal of Central South University(Science and Technology) , 2020 , 51 ( 4 ): 931 - 943 . Aerodynamic performance on the trailing edge split flap and regulation law of deflection angle of VAWT TAMBURINI A , CIPOLLINA A , MICALE G , et al . CFD simulations of dense solid-liquid suspensions in baffled stirred tanks: prediction of suspension curves [J]. Chemical Engineering Journal , 2011 , 178 ( 1 ): 324 - 341 . CFD simulations of dense solid-liquid suspensions in baffled stirred tanks: prediction of suspension curves JUNG S P , KIM Y G , PARK T W . A study on thermal characteristic analysis and shape optimization of a ventilated disc [J]. International Journal of Precision Engineering and Manufacturing , 2012 , 13 ( 1 ): 57 - 63 . A study on thermal characteristic analysis and shape optimization of a ventilated disc 车科学 , 聂万胜 , 何浩波 , 等 . 低湍流度风洞试验 [M]. 北京 : 国防工业出版社 , 2018 : 4 - 7 . 低湍流度风洞试验 CHE Kexue , NIE Wansheng , HE Haobo , et al . Low turbulence wind tunnel test [M]. Beijing : National Defense Industry Press , 2018 : 4 - 7 . Low turbulence wind tunnel test GJB 1179A—2012 . 低速风洞和高速风洞流场品质要求 [S]. GJB 1179A—2012 . Flow field quality requirements for low speed and high speed wind tunnels [S]. 高丽敏 , 刘哲 , 蔡明 , 等 . 四种风洞收缩段流场特性对比 [J]. 航空动力学报 , 2020 , 35 ( 8 ): 1695 - 1705 . 四种风洞收缩段流场特性对比 GAO Limin , LIU Zhe , CAI Ming , et al . Comparison on flow field characteristics of four wind tunnel contraction sections [J]. Journal of Aerospace Power , 2020 , 35 ( 8 ): 1695 - 1705 . Comparison on flow field characteristics of four wind tunnel contraction sections 10.11817/j.issn.1672-7207.2023.05.011 TB472 A 1672-7207（2023）05-1758-11 黄亚飞, 蒋兴良, 郑华龙, 等. 基于旋转运动的试验气流场产生方法数值研究[J]. 中南大学学报(自然科学版), 2023, 54(5): 1758-1768. 黄亚飞, 蒋兴良, 郑华龙, 等. 基于旋转运动的试验气流场产生方法数值研究[J]. 中南大学学报(自然科学版), 2023, 54(5): 1758-1768. HUANG Yafei, JIANG Xingliang, ZHENG Hualong, et al. Numerical study on experimental flow field generation method based on rotating motion[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1758-1768. HUANG Yafei, JIANG Xingliang, ZHENG Hualong, et al. Numerical study on experimental flow field generation method based on rotating motion[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1758-1768. 国家自然科学基金资助项目(51637002) Project(51637002) supported by the National Natural Science Foundation of China 2022-01-20 2022-03-17 2023-05-26 2023-06-20 38910308 中南大学学报(自然科学版) 5 54