流体力学在体育运动中的应用越来越广，在体育运动中的研究和进展正为我国的体育事业提供着借鉴，其研究成果在体育运动的各个方面已发挥了越来越大的作用。下面是搜索整理的流体力学论文8篇，供大家借鉴参考。
流体力学论文第一篇：双层U型管中有无挡板及不同位置的流动特征分析
摘要：针对双层U型管中有无挡板及其安放位置对流动特征的影响展开研究。热水入口流速0.5m/s，冷水入口流速1.0m/s，使用商用计算流体动力学软件ANSYS FLUENT进行流场数值模拟，得到了管道内流体速度、涡量、压力和温度的分布。设置挡板后，管道内湍动剧烈，涡强度增强，平均流速先增大后减小。设置挡板可以显著改善隔板出口处的压力分布情况，而随着挡板到隔板距离的增大，挡板处的压差逐渐增大。无挡板时，出口处流体间换热效果最好，设置挡板后，随着挡板到隔板距离的增大，换热效果变差。
关键词：双层U型管; 数值模拟; 流动特征; 隔板; 挡板;
双层管道由于其结构的特殊性，使得2股流体被中间隔板分开，汇合之后可起到促进流体之间混合和热交换，因而被广泛应用于石油化工、能源动力等行业。国内外对工业弯管和T型管道内的流动特性的研究较多，而对于小尺寸的双层U型管内的流动特性的研究较少。无论从尺寸的层面还是从流动混合的角度，双层U型管内隔板出口处的流动特征都具有重要的研究价值。
挡板常在工业中作为改变流体流动方向的装置，其结构和位置影响管道内液体的流动与换热。Graham等[1]对U型带肋通道中的流场进行了测量，发现肋前后存在明显的循环区域。赵国寿等人[2]研究发现叶片表面障碍物能够增加近壁面流场湍动能从而使转捩提前发生，减少流动分离引起的压差阻力，从而对空化产生了抑制作用。Lim等[3]通过大涡模拟计算方法和试验方法研究了湍流边界层中的贴壁立方体的流场结构和流动情况，分析了立方体周围的涡结构与平均压力和脉动压力的关系。王庆锋等[4]研究了折流板间距与换热器对流传热系数的关系，研究发现折流板间距越大，对流传热系数越小，压降越小。
基于优化双层管道结构以提高换热性能的初衷，针对小尺寸的某双层U型管道结构，在隔板出口处布置挡板，研究挡板到隔板出口的距离对管道内流动特性和换热效果的影响。由于模型尺寸小，在实验的实施上存在一定的难度，所以采用数值模拟的方法对管内的流动特征进行初步研究。本文在数值模拟中采用ANSYS Fluent软件，借助k-ε湍流模型对双层U型管内的流动特征进行研究，重点对挡板到隔板出口距离的大小对隔板出口到挡板之间区域的流动特征和管道出口处流体温度分布的影响进行描述与分析。
1 模型的建立
1.1 几何模型
某双层U型管几何模型如图1所示，流体分别由入口1和入口2流入，在流过180°的弯管后，撤去隔板，2个管道合并成1个管道，入口截面采用边长为5mm的方形布置，出口截面为长35mm、宽15mm的矩形，隔板厚5mm，弯管半径为100mm。分别在距离隔板出口10mm、20mm、30mm和50mm处布置挡板，建立4种模型并建立无挡板模型作为参照，建立三维坐标系（图1），坐标原点位于隔板出口。挡板结构如图2所示，挡板宽15mm、长10mm、高5mm。挡板到隔板出口的距离用字母d表示。
图1 模型几何结构
1.2 边界条件与网格划分
数值计算采用k-ε湍流模型，使用清水作为流体介质。入口1流速0.5m/s，温度60℃（333K）；入口2流速1.0m/s，温度20℃（293K）。壁面采用壁面无滑移条件，设置管壁为绝热壁。出口压力设为环境压力，值为101325Pa，环境温度20℃（293K）。收敛标准是各监测项的残差均小于10-4。
图2 挡板几何结构
网格划分采用六面体结构网格，每个模型计算域共有100万个网格。本文重点研究隔板出口到挡板之间的流动特征与温度分布，所以对隔板厚度的范围内沿x方向加密处理，并对挡板处的网格进行加密，如图3所示。
图3 局部网格划分
2 数值模拟与分析
2.1 速度分布
如图4所示，挡板到隔板距离的大小对速度分布有较大影响。隔板前方来流方向的流速比隔板后方的流速略大。在高流速区域，无挡板时，最大流速等值线分为2段，增加挡板后，两段最大流速等值线逐渐靠近，最终2段等值线合为1条等值线，随着d的增大，挡板对隔板出口处流体流动的扰动作用减弱，最大流速等值线再次逐渐断裂，d=50mm时，断裂为3段等值线，随着d的增大，高速区域平均流速经历了先增大再减小的过程。在低流速区，d=10mm时挡板对隔板出口处流体流动的扰动作用最强，低速区域平均流速最大。d=30mm和50mm与无挡板时在低速区隔板附近的流速等值线形态相似，量值略有差别。随着d的增大，挡板附近来流方向的流速逐渐增大，挡板后尾流的流速也逐渐增大，但是尾流的范围先增大后减小，d=20mm时挡板后尾流的范围最大。挡板对隔板后的尾流有较大影响，当挡板距离隔板较近时，挡板处于隔板后的尾流范围内，流体阻力增大导致尾流内的低速区域范围增大，流速降低；随着d的增大，挡板对隔板尾流的阻力作用减弱，但d=50mm时隔板尾流的范围仍然比无挡板时的隔板尾流范围大。
图4 速度分布
2.2 涡量分布
如图5所示，隔板与挡板靠近高流速的区域比靠近低流速的区域高涡量区域的范围更大，且涡旋尾部延伸较长。2高涡量区之间低涡量的形态沿X方向上下波动。如图5(b）所示，当d=10mm时，挡板处于隔板高涡量区范围内，涡尾沿Z轴方向变宽。隔板附近高涡量区与挡板附近的高涡量区连成一体。如图5(c）所示，当d=20mm时，隔板出口下方高涡量区域尾部上下波动幅度更大，尾部形成一个微小的高涡量区，隔板与挡板附近上方的高涡量区分离，隔板中间的低涡量区域与挡板上方低涡量区域连接成一体。挡板上壁面高涡量区域长度变大。如图5(d）所示，当d=30mm时，挡板上壁面高涡量区域进一步扩大，延伸至挡板圆周前端，挡板后端下方的高涡量区域拉长，范围变大。如图5(e）所示，当d=50mm时，隔板出口处下方的高涡量区尾端出现多个封闭的等值线，涡量大小变化频繁，同时尾部形状较无挡板时波动更剧烈。d=10mm时较无挡板时，隔板出口处的高涡量区范围缩小，随着d的增加，隔板出口处的高涡量区域范围逐渐增大，涡形态上下波动更剧烈。有挡板时，随着d的增大，挡板壁面附近的涡强度增强，挡板后的高涡量区涡强度增强，涡尾的长度延长。
2.3 压力与温度分布
隔板出口处压力沿Z轴分布如图6所示，隔板两侧的压力较高，在隔板沿Z轴的宽度范围内形成一个低压区。无挡板时隔板出口处的低压值最小为负压，与隔板两侧的压差最大；当d=10mm时隔板出口处的低压值最大，在20Pa左右，与隔板两侧的压差最小。有挡板时，随着d的增大，隔板两侧的压力逐渐增大，隔板处低压值逐渐降低，压差增大。因此，d=10mm时，可以显著改善隔板出口处的压力分布情况。有挡板时，在挡板附近出现管道内的压力最大值和最小值，如表1所示。随着d的增加，压力最大值先增大后减小，d=30mm时压力的最大值最大；压力的最小值逐渐减小，且减幅增大；挡板附近的最大压力和最小压力的压差随着d的增大而增大。
图5 涡量分布
图6 隔板出口处压力分布
流体出口截面的温度分布如图7所示，无挡板时如图7(a）所示，出口截面中间区域等温线较为平缓，温度在310K到319K之间，两侧等温线向外鼓起，且鼓起程度逐渐增大。有挡板时，随着d的增大等温线倾斜程度增大，左侧等温线密集，温度梯度大，冷热水之间热量交换更加不均匀；同时随着d的增加，高温区域范围逐渐增大，d=10mm时出口最高温度为331.67K,d=50mm时出口最高温度为332.21K，换热效果变差。与无挡板的情况相比，d=10mm时，出口最高温范围略有扩大，高温区域等值线鼓起程度较小，温度梯度增大，说明挡板的设置不利于2股流体间的热量交换。从图7可以看出，出口截面下方的等温线在形态和位置上几乎无变化，出口截面上方等温线变化显著，因此挡板对于热水区域的热量交换的影响更为显著。
表1 摇挡板附近压力极值
图7 出口温度分布
3 结论
在某双层U型管内隔板出口后方布置挡板，对管道内的流动结构有显著影响，可以得出下面的结论。
设置挡板后，管道内湍动剧烈，随着d的增大，高涡量区域范围扩大，涡强度增强，隔板出口处的涡形态上下波动更加剧烈，挡板后方的高涡量区域的涡尾变长；管道内平均流速先小，流速增大；设置挡板可以显著改善隔板出口处的压力分布情况，d=10mm时隔板出口处与隔板两侧的压差最小，随着d的增大，隔板两侧的压强升高，隔板后的压强降低，压差增大，无挡板时隔板后低压区出现负压，压差最大，而随着d的增大，挡板处的压差逐渐增大；无挡板时，出口处流体间换热效果最好，设置挡板后，对热水区域的换热影响较大，冷水区域的换热情况基本无变化，随着d的增大，换热效果变差。
参考文献
[1] Graham A,Sewall E,Thole K A.Flowfield measurements in a ribbed channel relevant to internal turbine blade cooling[R].ASME Paper2004-GT-53361,2004.
[2] 赵国寿.离心泵内空化流动不稳定性及其控制研究[D].兰州理工大学，2018.
[3] Lim H C,Thomas T G,Castro I P.Flow around a cube in a turbulent boundary layer:LES and experiment[J].Journal of Wind Engineering&Industrial Aerodynamics,2009,97(2):96-109.
[4] 王庆锋，庞鑫，赵双.管壳式换热器传热效率影响因素及数值模拟分析[J].石油机械，2015,43(10):102-107.
[5] 杨龙滨，高振森，吕玉庭.不同挡板开度下静态分离器分离性能数值模拟[J].黑龙江科技学院学报，2009,19(04):253-257.
[6] 张静，刘晓亮，龚斌，等.挡板的相对曲率对分离器入口局域流场流体力学性能的影响[J].化工进展，2017,36(11):3963-3970.
[7] 童传琛，方弘毅，康涌，等.U型管式回热器设计与流动换热性能试验研究[J].燃气涡轮试验与研究，2017,30(02):17-21+31.
[8] 刘喜岳，张靖周，李刚团，等.双U型管束模型换热器的流动和传热特性[J].航空学报，2015,36(12):3832-3842.
[9] 林泽鹏，包芸.隔板对流系统中影响流动和传热的压力分布特性[J].中国科学：物理学力学天文学，2018,48(10):94-102.
[10] 林泽鹏，包芸.隔板对流系统中影响流动和温度特性的关键几何参数[J].中国科学：物理学力学天文学，2018,48(05):60-68.
[11] 李晓敏，王立军，王河.管壳式换热器流体与传热模拟分析[J].真空科学与技术学报，2018,38(10):919-923.
[12] 刘晓亮.结构参数对挡板局域流场特性的影响[D].沈阳化工大学，2018.
[13] 姚甜.肋片上游障碍物扰流特性及其对换热的影响研究[D].南京航空航天大学，2018.
[14] Li Y,Chan C K,Mei B,et al.LES of incompressible heat and fluid flows past a square cylinder at high Reynolds numbers[J].International Journal of Computational Fluid Dynamics,2015,29(3-5):272-285.
[15] 邱庆刚，沈胜强.矩形肋片导流角度对内冷通道流动与换热特性的影响[J].热科学与技术，2008,7(4):331-336.
[16] 王赛.隔离板对振荡柱体尾流控制的研究[D].中国计量学院，2012.
流体力学论文第二篇：旋涡的动力学涡定义方案探析
摘要：针对旋涡定义这个长期未解决的难题，在分析历史上提出的几种思路的基础上，本文提出：定义旋涡的总目标，应当体现其管状运动形态，并与其作为流体运动肌腱的动力学功能有机结合；该定义应能以统一的方式覆盖经典涡动力学业已成熟的结果，并且指引对湍流中复杂涡状结构及其相互作用的识别。据此，本文根据涡量场演化的已知规律，提出了旋涡的一个动力学定义，作为继续深化讨论的参考。同时提出，人们针对湍流结构研究构造的各种涡判据，虽然对湍流涡状结构的可视化起了重要作用，但因其难以遵循管式涡量场的因果演化而无法取代完备的涡定义。相比之下，近十年来形成的涡量面理论，有望在未来发展中使旋涡的定义和检验问题回归涡量动力学的自然规律。
关键词：旋涡; 涡动力学; 涡定义; 涡判据; 湍流结构;
Thoughts on vortex definition
WU Jiezhi YANG Yue
State Key Laboratory for Turbulence and Complex Systems,College of Engineering,Peking University
Abstract：Aiming at the long unsolved issue of vortex definition,this article reviewed critically those historically appeared definitions,and thereby proposes that the general goal of vortex definition should reflect an inherent combination of the solenoidal nature of vortices in kinematics and their role as the sinews and muscles of fluid motion in dynamics.The definition should,in a unified manner,cover all mature achievements of classic vortex dynamics and guide the identification of complex vortical structures and interactions in turbulence.Thus,on the basis of the known evolution laws of vorticity field,we propose a dynamic definition of vortices as a reference of further discussion.Although various “vortex criteria”constructed for visualizing turbulent structures have played important role in one's understanding,they cannot replace a complete vortex definition since they do not follow the causal evolution rules of the solenoidal vorticity field.In contrast,it is expected that the vorticity-surface theory formed in the past decade will,in its future development,be able to bring the definition of vortex and its testing back to the natural laws of vorticity dynamics. 0 引言
根据动量方程的Helmholtz分解可以发现，单一性质的流体运动有两个基本动力学过程：以涡量为特征变量的横过程即剪切过程，和以胀量或热力学变量（压力、密度、焓、熵）为特征变量的纵过程即胀压过程（Wu et al.2015)[1]。纵过程的特征结构很清晰：声波、激波以及势流，它们有明确的数学定义。纵场理论的难点在于其特征量不唯一，对这些结构的数学描述与源的辨识也不唯一。横过程则与此恰恰相反。其特征量和数学描述唯一，但偏偏其最重要的基本结构———旋涡或涡，没有公认的定义。这已成为流体力学的一个老大难问题。本文在回顾经典定义演化过程的基础上，对“什么是涡”这个问题提出自己的初步看法，抛转引玉，以期促成更广泛的讨论。
1 涡的运动学定义及其动力学功能
横场有一个独特的纯运动学性质，即涡量场的无散性。它的直接结果是三维无界流场的总涡量守恒和二维流场的总环量的时间不变性。特别对于涡量管，它导致Helmholtz第一涡量管定理（简称H1）：涡量管的强度或其环量与管面上计算它的位置无关（不同著作对Helmholtz三个涡量管定理的编号与说法不同，这里据Truesdell(1954)[2]）。无散性决定了“涡”（vortex）最根本的普适形态特性，是自Helmholtz以来人们把细涡量管称为涡丝并简称为涡的原初依据，业已指引人们理解了若干重要的物理现象。最著名的经典事例来自空气动力学。对于二维流，Prandtl从翼型突然起动形成的起动涡，根据总环量守恒推断出产生升力的等值反号的翼型环量（图1)[3]，而这个环量只能是翼型上下表面边界层内涡量的代数和。对于三维流，Lanchester最先根据H1从翼尖涡推断出它们必然来自机翼表面边界层里的涡量在翼尖附近的汇聚和折转（图2a)[4]，这启发Prandtl在1918年基于H1提出了升力线理论的涡系模型（图2b）。没有涡量的无散性，这些认识是不可能的。更一般地说，如果某一时刻在流场某处产生了新的涡量场，则必在同一时刻会在流场别的地方产生新的总量相同、方向相反的涡量场，无论这两处涡量场中的一个或两个是否发展成集中的旋涡。
图1 从左向右突然起动的翼型的起动涡和升力涡[3]
图2 三维尾涡
从涡量场的无散性到涡定义，却是一段坎坷的历程。回顾过去，可看到四种不同的思路。
1.1 基于涡量管的运动学定义（Helmholtz-Lamb)
从Helmholtz(1858）把细涡量管称为涡丝（vortex filament），而后Lamb(1932）把涡丝简称涡（vortex）开始，人们就用细涡量管即涡丝（vortex filament）定义旋涡（参见Wu,2018[5]）。这是涡定义的一个主流方向，因为它抓住了H1所体现的涡量场作为无散管式场（solenoidal field）的特性。但这个运动学定义并不完备。首先，Helmholtz和Lamb并未规定该细涡量管之外有没有涡量。设有一族连续分布于空间的涡量线，从空间中任何不是涡量线的闭曲线出发均可生成一个涡量管。由于在这条闭曲线的邻域还可引入无限多与之不相交或相交的闭曲线（如奥运会的五环会标），从这些环又可生成一系列涡量管，或者互不相交如一套涡量管，或者相互嵌套的涡量管。但这样生成的结构显然不唯一，不能把一组涡量管中的任何一个识别为通常意义下的旋涡。
其次，人们通常认为，根据H1，在无界流体中涡量管和涡线一定闭合成环，并据此讨论过大量物理问题。但Fuente(2007)[6]举例证明，涡量线和涡量管并非如人们常说的那样要么自行封闭要么延伸到流体边界。涡量线可能无穷长但局限于有限空间，或者在流体内部终止于零涡量点；而依其生成闭曲线的不同，涡量管可以自行封闭，也可以是无穷长而且自相交无数次，但不封闭。可见，任意选定的涡量线和涡量管的拓扑结构并不简单。
事实上，基于涡量管的纯运动学定义的确是不充分的，因为严格的不可压流模型也有无散性，其流管也满足H1；而且若流动无黏，Helmholtz第二、三涡量管定理（H2,H3）对流管也成立。但除了计算流量之外这些流管并无其他重要特性。
1.2 基于流动有旋性的定义（Kelvin-Saffman)
和上述基于涡量管的主流定义不同，Kelvin(1869,p.225)[7]在关于涡运动的经典论文中提出，“I now define a vortex as a portion of fluid having any motion that it could not acquire by fluid pressure transmitted through itself from its boundary。”意思是，旋涡是一团任意运动的流体，该运动并不来自从其边界传输过来的流体压力。显然，这部分流体只能理解为整个涡量场。这个定义排除了不可压流的流管。Saffman&Baker(1979)[8]给出的涡定义实际上是Kelvin定义的某种直白说法：旋涡是被无旋流包围的一个有旋流体。由于涡量场的管式特性，KelvinSaffman定义实际上指的是最大可能的涡量管，其侧面就是一个有界涡量场的外边界。此定义的好处是，前述涡量管定义的不唯一性问题，以及任选的涡量线与涡量管的复杂拓扑结构问题，就都不存在了。
然而，如果任何一个有界的涡量场都能被整个定义为一个旋涡，这和人们的直观理解相去甚远。于是，为避免过于宽泛，Saffman(1992)[9]又补充说，涡主要指无旋流包围的细涡丝。值得注意的是，当Kelvin提出他的涡定义时，隐含地用到无旋流中联系压力和动能的Bernoulli定理。而这个定理已经属于动力学，但却在Saffman的定义中被忽略了。
1.3 涡的功能性刻画
由上可见，单靠涡量场的无散性或者流体的有旋性来定义涡，都是不充分的。它们无法解释，为什么并不是所有细涡量管或有旋流都具有必须给以特殊关注的特性。例如，在小雷诺数下忽略惯性力的定常Stokes流中，压力的梯度（除以黏性系数）和涡量的旋度互相平衡，在二维流中它们无非是一个复解析函数的实虚部（Wu et al.2006,p.151)[9]，虽然人们总可以画出许多涡量管，却没有一个能被称为涡；而作为有旋流的整个流场也不具有任何旋涡的意义。为了定义涡，必须进入动力学，而且不能像Kelvin那样简单地用是否满足Bernoulli定理作为判据。
在这方面，Küchemann(1965)[10]沿着人们对涡的物理直觉迈出了重要的一步。他全然不顾涡的运动学定义，提出了涡是“流体运动的肌腱”这个著名命题。随后Moffatt,Kida&Ohkitani(1984)[11]又把这句话延伸到湍流，认为被拉伸的涡是湍流的“筋”（sinews）。这就在实际上给出了一个阐明涡的功能的形象化判据，它强调了只有强到足以成为肌腱的细管状结构才够格被称为涡。这个命题符合人们的物理直觉，因而得到广泛的认可。但涡在什么意义下能成为流体运动的肌腱？这个问题也只能到动力学中去寻找答案。
1.4 基于旋转度量的涡判据
在经典涡动力学中，旋涡的理论模型总以孤立的形态出现，或者是无穷长的单个或一对直涡，或者是涡环。这种理论模型在湍流研究中不够用了，而代之以复杂的、相互纠缠和相互作用的各种“相干结构”或“涡状结构”（vortical structures）。虽然人们确信这些结构都来自涡量场，但长期无法给出它们的明确定义，而只能用不同实验手段、在不同阶段和条件下观察到的现象，对这些结构给出瞎子摸象式的形象化名称。这种状况导致了各种“涡判据”的诞生。
早先，Truesdell(1954,Secs.51-62)[2]基于涡量是有量纲的量这个事实，提出了需要一个无量纲的数来度量流体运动旋转程度的问题（“to find a quantitative measure of the amount of rotation in a motion”）。他基于Helmholtz分解（他根据历史的考察称之为Cauchy-Stokes分解），用涡量与应变率张量的两次缩并之比定义了运动学的涡量数（vorticity number），并用Lamb矢量的绝对值和物质加速度扣除Lamb矢量后的绝对值之比定义了一个动力学的涡量数。前者已成了研究湍流相干结构时引入的各种“涡判据”（如见Epps,2017)[13]的滥觞。这些运动学判据实质上都对涡的强度给出了规范。虽因相互有别而导致对同一结构有不同理解，但它们对于可视化湍流中复杂的涡状结构毕竟起了重大作用，做出了不可替代的贡献。在公认的更好且易用的高效涡识别方法出现之前，这些判据还将较长时间发挥作用。但它们和涡定义不属于同一个概念（Wu 2018)[5]。从运动学上看，现有的判据都和涡量场的管式特性无关，识别出来的结构都不满足H1。满足H1的一个必要但远非充分的条件是（Wu et al.2005)[14]，涡不能因轴向拉伸（导致涡量管变细、涡量强度与卷吸的动能增加，但保持环量沿管子不变）而被不合适的判据拉断。根据Koláǐ&ístek(2019)[15]最近的检验，在几个常用的判据中只有称为Rortex的判据（如Gao&Liu 2018)[16]能满足这个必要条件。然而，即使拉不断，Rortex的外边界仍未必能近似保持涡量场无散性导致的管状性质。
从动力学上看，这些判据最根本的缺陷，是它们不能统一地揭示涡状结构演化中不同时刻可识别结构之间的因果关联。它们是时空中孤立的存在，不是涡量场本身演化的运动学-动力学规律的自然结果，也就无法由之导出各种极有价值的经典涡结构模型及其动力学理论，例如边界层、面涡、涡丝、点涡、涡斑（vortex patch）等等。比如Prandtl升力线理论的机翼束缚涡和线化自由尾流面涡这样简单的理论模型（图2b），乃至涡动力学中所有业已证明有效的简化模型，都是用这些判据难以得到的。
可以看到，上述前三个思路虽然都不完备、不充分，却各自抓住了涡的某个方面的特征。一个好的涡定义应当把这些特征有机地整合起来。关键词是管状、强旋转。需要把“肌腱”这种形象化的描述改造成精炼而完整的物理理论概念，即进一步阐明什么样的细涡量管能“升级”成流体运动的肌腱。而为此就得解释它们从何而来、到何处去的动力学来源与演化。同时，作为管和“肌腱”，就只能包括“轴状涡”。有了好的定义之后，发展量化的识别判据就有了明确的目标和方向。
总之，涡量场的无散性和H1是任何旋涡必备的独特运动学形态，但基于涡量管的纯运动学涡定义并不充分，还必须体现旋涡作为流体运动肌腱的独特动力学功能，而只顾及动力学功能而不顾H1的定义也是不充分的。因此，定义旋涡的总体目标应当是：（1）能将体现其管状的运动学形态与肌腱的动力学功能有机统一；（2）能以统一的方式覆盖经典涡动力学业已成熟的结果，并且指引对湍流中复杂涡状结构及其相互作用的识别。
下面是我们沿着这个总体思路的一个初步探索。本文所用的论证材料大多可在Wu et al.[1,9]中找到。
2 有效无黏流：从面涡到旋涡
先从雷诺数趋于无穷的有效无黏流说起。这时，涡量场的基本形态是无限薄面涡（vortex sheet，或称涡片），即切向间断面。可能是Prandtl最早提出了面涡会卷绕成轴状旋涡的思想。Kaden(1931)[17]首先用二维流证明，一个有边界的面涡，如果它在边界上的涡量是非零有限值，就会“诱导”面涡迅速卷绕成很紧并不断伸长的旋进（spiral in）单螺旋结构。随卷绕层数增加、层间距离和每层面涡的强度趋于零，这个结构的横截面形成卷紧的闭圈（图3）。
这种机制已由大量实验证实。据此，Betz[18]令人信服地论证：在外无界的小黏性不可压流体中，旋涡只能来自面涡卷绕。他先假定，如果在静止流体中搓一根细圆柱令其旋转，最后把圆柱取出，那么在无限大雷诺数下旋转圆柱得用无限大的时间才能完全带动流体转动，而且其动能需要无穷大的功率输入。而事实是一块平板一动，边缘就搓出面涡，快速卷绕而成旋涡。Betz的这个论证解释了从飞机到昆虫的定常与非定常空气动力学中遇到的各种旋涡的形成。从Krasny(1987)[19]发明消除奇异性的数值研究开始，这些螺旋状结构已经被陆续用不同的方法证实，数值方法先是针对二维流，后来拓展到轴对称流和一般三维流。
图3 一端有界的面涡卷绕，（Betz,1950[18])
此外，人们还从实验、理论和计算中发现，一个两端无限长的平直面涡，也会因Kelvin-Helmholtz不稳定性而非线性地形成旋进-旋出（spiral in and out）的双螺旋结构，而形成一个涡列（图4）。在一般的情况下，二维面涡的这种自诱导性质服从高度非线性的BirkhoffRott奇异积分-微分方程。值得强调的是，面涡卷绕成集中涡这个机制，是在Helmholtz的工作70多年之后发现的，代表了人们对涡结构认识的重要进展。
图4 受扰无限长面涡因Kelvin-Helmholtz不稳定性和自诱导而发展成卷绕涡列（Krasny,1988[19])
对三维流，一个轴状涡一旦形成，背景流场沿轴向的加速度就会导致涡的拉伸而使其所含的拟涡能指数加强。这是独立于卷绕的另一个强化旋涡的运动学机制。浴盆中携带涡量的流体在出口处由于重力加速度的拉伸作用而形成的浴盆涡就是这种机制的一个简单而常见的例子。更重要的典型例子，是积雨云中由斜压效应产生的强烈翻滚的有旋流体由于上冷下热的强对流而拉出的龙卷风（图5）。在地面上的飞机发动机形成的小型龙卷风状旋涡也源于类似的机制。这种热对流和地球旋转的科氏力与风切变的相互作用，还能形成更大尺度的台风甚至飓风。
图5 积雨云和龙卷风[20]
连续性方程的运动学效应和动量方程不同分量间耦合的动力学效应显示，涡的轴向拉伸往往伴随着向心的径向速度，导致周围流体被卷并到旋涡中而形成动能的集聚，使旋涡成为流动的组织者。这正是涡之所以成为流体运动的肌腱的重要原因。当然，涡的轴状运动学形态也正是筋与肌肉的形态。
在湍流中，涡因拉伸而加强只是个大概率事件，有时也会因相反的轴向加速度导致轴向收缩、旋涡变粗和涡量密度减弱。但沿轴向挤压一根细橡皮棍最容易导致的结果是，一个微小的非轴对称扰动就把皮棍挤弯；同样，在有轴向挤压的背景场中，旋涡发生弯折也是比它单纯变粗的概率大得多的事件。只要把一条直涡丝弯成U状，其两条平行反号的腿所诱导的动能必定大为减小。因此，Chorin(1994)[21]指出：由于能量守恒，涡的拉伸必然伴随着它的弯折。在涡丝动力学中，弯折的涡丝要沿副法线方向运动，主控方程是局部诱导近似的非线性薛定谔方程。而涡丝的拉伸和弯折本身则涉及非局部自诱导和背景流效应，主控方程是更复杂的积分-微分方程（Wu et al.2006,401-407;413-418)[9]。
现在，如果我们把图4那个最初平直的无限长面涡换成复杂的曲面，并把流动换成三维的，则可以想象这样的不稳定曲面也大量存在于图5的积雨云内部的湍流中。这里既有面涡的卷绕效应，又有轴状涡的拉伸效应和折转效应，其复杂程度可想而知，而在低空显得极为强烈的龙卷风，也只是这庞大的积雨云的一个小尾巴，是个有一定概率的事件。
值得注意的是，在有效无黏流中，面涡卷绕的最终结果不会是闭合涡量管，而是无限多层紧密卷绕的螺旋结构，它和涡量管的拓扑性质还是不同的，所以H1用不上。另一方面，这个螺旋结构的轴心是条直径无限小的奇异线，那里的周向速度和动能趋向无穷大。如果奇异线是弯曲的（通常必然如此），它还会有无穷大的自诱导速度。这两个问题只能靠引入黏性才能解决，讨论如下。
3 黏性的作用和Moore-Saffman尾涡模型
真实的流体都有黏性，引入黏性表明问题进入了动力学。黏性有几个方面的作用。首先，面涡变成有限厚度的剪切层，但这个厚度必须足够薄才能使剪切层发生卷绕。剪切层这样的结构只在小黏性流体或大雷诺数流动中才出现。这是旋涡出现的一个必要的动力学条件。其实纵场的结构也是如此。很黏的流体既没有旋涡，也没有长程声波和激波。换言之，纵横场的结构都是由大雷诺数动力学过程决定的。
黏性的第二个作用，是消除前述卷绕面涡不可避免的奇异性。在大雷诺数下，剪切层的厚度为的量级，L是整体长度尺度。涡层卷得越紧，其强度也随之减小，涡量分布和速度分布越趋向光滑，如图3所示。当相邻剪切层间的距离小到每层的厚度这个量级时，卷绕结构就变模糊，退化为环量Γ=O(1）的黏性涡量管，特征雷诺数是。
黏性的第三个作用，是不仅使Helmholtz第2、第3涡量管定理失效，而且在抹光卷绕面涡轴心的奇异性的同时，也使管状涡核失去了清晰的外边界面。所以，这样得到的旋涡，乃至于涡量场本身，本性上就是没有清晰外边界的黏性结构，正如边界层和自由剪切层没有清晰的外边界一样，有的只是指数衰减的“尾巴”（所以边界层厚度的定义不唯一）。顺便提一句，黏性激波层也是这种物理上紧致的结构，但激波层没有自诱导，不会卷成螺旋而已。这样一来，我们不得不承认，涡核的外边界是个有限厚度的“外壳”，大体上对应于从核内切向速度单调增长转变成在核外势流中代数衰减的过渡层（参见图7a）。Saffman在和本文第一作者的一次交谈中（约为1994年）也承认，他把涡定义为无旋流包围的有旋流体，只适用于无黏流，在黏流中涡量场没有清晰的外边界。下面，我们姑且把作为涡核外边界的有限厚度的“外壳”称为准涡量管。这里重要的是，在准涡量管内必存在满足H1的严格涡量管，在每个瞬间其环量沿管不变。
为了具体理解在黏性流体中涡层卷绕成轴状涡并形成光滑涡核的过程，不妨回顾Moore和Saffman(1973)[22]根据Kaden相似律，针对飞机尾涡结构提出的飞机尾涡模型（简译MS模型）。尾流面涡从后缘起的三维定常卷绕过程可按xUt化简为一个横流截面上的二维非定常问题。首先考虑初始时刻在（x,y）平面上的一个位于从机翼左端起y>0的半无穷大平直面涡，强度分布为γ0(x)～y-n,n∈[0,1](n=0.5代表环量有椭圆分布）。则这个面涡会因自诱导而迅速卷起，在极坐标（r，θ）下可用量纲分析证明卷紧的结构服从Kaden相似律（β是个常数）：
Moore和Saffman用渐近匹配展开证明，尾涡有个三层状的核结构：最外层是松散的无黏卷绕面涡，可用Kaden相似律定义一个核半径a0；内层是黏性子核以消除切向速度的奇异性，半径为avis；中间是从卷绕面涡过渡到黏性子核的光滑无黏区，其外边界为asm。见图6。
图6 Moore-Saffman尾涡模型的核结构
MS模型包含了轴向压力梯度引起的拉伸-收缩作用。据其算得的尾涡周向速度V与轴向速度W的剖面取决于机翼载荷参数n（见图7），该参数决定着轴向流的分布和正负号（射流型还是尾流型）。这个模型的预测结果已被一些实验所证实。
当然，MS模型不能直接搬用到其他旋涡上。但是其定性的物理特征具有普遍意义。首先，它体现了Betz的论断，即只有小黏性流体中的涡层卷绕才能形成旋涡；其次，它证实构成旋涡的闭合涡量管不是任何运动学意义下的涡量管，而是涡层卷绕与黏性扩散相抗衡的产物，由于卷绕过程把大片面涡携带的涡量集聚到狭小的管状区域内，这种涡量管的涡量的确有流体运动肌腱的强度。在t=O(1）时，黏性子核半径avis=O(Re-1/2）处的切向速度达v=O(Re1/2）。由于周向速度和轴向速度的动力学耦合（参见Wu et al.2015,Eq.(6.4.8))[1]和质量守恒，这种强涡量管能把周围流体卷吸过来而成为动能高度集聚的结构。第三，H1对光滑核内的流体成立，这对于旋涡定义是十分必要的。对于椭圆载荷分布，在t=O(1）时光滑核半径为asm=O(Re-1/5），其内有闭合涡量管。第四，它揭示了卷绕之前的剪切层强度分布对整个旋涡性态的影响。这些特性是旋涡的大多数理论模型所没能做到的。值得一提的是，在湍流级串过程中，这种卷绕结构也是小尺度间歇性的来源。Lundgren(1982)[23]首次证明，带旋臂的拉伸涡模型的系综平均具有-5/3的能谱。与此相反，由于湍流研究中构造的各种涡判据只能辨识轴状结构，都不可能揭示剪切层卷绕成涡的过程，必然对MS模型和Lundgren模型掐头去尾。
图7 Moore-Saffman尾涡模型的周向与轴向速度分布
4 旋涡的一个动力学定义
作为对以上三节讨论结果的概括，如果要写一个“旋涡”的词条，我们可以试着给旋涡下一个新的动力学定义：
旋涡，或简称涡，是有旋流体在各种运动学效应和大雷诺数动力学效应下，通过剪切层卷绕和涡量拉伸自组织成的管状结构。它们在流体运动中导致动能向其周围的高度集聚，起着肌腱的作用。
如果再多说一点，可以加上：其光滑涡核中的涡量管满足Helmholtz第一定理，在其自身演化和与其他流动结构相互作用的过程中，每个旋涡都经历从生成到消亡的过程。
这个定义排除了Helmholtz-Lamb基于涡量管的运动学定义所过多包含的东西，而添加了作为流体运动肌腱的动力学机制，包括黏性的必要性。定义本身无需列入各种运动学和动理学的具体效应以及它们发生的具体条件和流动参数范围。这些具体化的分析是需要在词条的定义之后阐明的。我们期望，基于这些分析，能够证明这个定义涵盖了过去150年来所有人们称之为旋涡的横场结构的共同特征，包括航空航天、流体机械、风工程、大气、海洋等自然界和工程技术领域遇到的旋涡和湍流中的涡。否则，定义还需修改。
在本文第一作者过去的论著中，曾或隐或显地把旋涡分类成层状涡，即剪切层和轴状涡。现在看来这个分类不妥，因为没有卷绕的层状涡不具有管状结构。除非硬把剪切层中的任何运动学意义下的涡量管都算进去，否则它与Helmholtz三定理无关。上述定义廓清了剪切层和旋涡之间的因果关系，把前者视为旋涡形成前的结构形态看来更为贴切。换言之，上述旋涡的动力学定义只适用于轴状涡。
值得强调的是，这个动力学定义并不要求旋涡必须闭合成环。一个非环状旋涡可以只有有限长度，在其一端或两端打开成非卷绕或弱卷绕的剪切层。一个重要的典型例子是图4的龙卷风，其底部和顶部散开，涡量线在底部折向，平行于地面进入边界层，涡量管呈喇叭状扩大；在顶部拐进猛烈翻滚的复杂云层。正是这个底部和顶部的边界条件决定了大气剪切层为何能卷绕成龙卷风，而这个过程迄今也没彻底澄清，使得龙卷风的预报仍是个难题。另一个例子是机翼前缘涡，它由前缘分离的剪切层卷绕而成，有上游端；在下游遇到后缘附近的逆压梯度会经历高度非线性的失稳过程突然破裂、散开，并转捩成湍流，故也有下游端。这也是一个有限长的旋涡。这种有限长旋涡的两端演化在传统模型中是无法描绘的，上述定义应当有助于推动人们更多地关注旋涡形成或终结的过程。轴状涡的这种有限长度，也使只能识别轴状涡的所有判据不可能分辨涡的来龙去脉。
有件事需要解释：因为允许涡量管打开成非管状的涡量面，在大多数实际情况下几乎无法实际检验构成有限长旋涡的所有运动学涡量管，包括它们在旋涡端部的延伸，是否满足H1。对这个困惑有个简单的回答———不必检验，只要真正抓住了涡量线、涡量面和涡量管，知道H1在理论上必然满足就够了。再去检验H1就有可能不必要地陷入复杂的几何拓扑问题，而这些纯几何的复杂性并不影响本文关于涡定义的讨论。
作为一个一般的概念定义，上述涡的动力学定义应能涵盖所有的旋涡，既包括湍流中复杂的旋涡，又包括空气动力学、地球流体力学等领域观察到的孤立旋涡，所以也需要受到这两类领域典型事例的检验。其实，这个动力学定义本来就是对经典旋涡理论与实验结果的概括，可以认为它已经得到了这类事例的检验（若有差误，定义还应修改）。目前需要进一步做的，是用湍流中的复杂涡状结构来检验它。在这种检验过程中，人们已经遇到经典孤立旋涡案例所观察不到的若干新现象需要解释，例如涡对的重联；也还将会继续发现新的现象，如涡相互作用引起的级串（Yao&Hussain,2020[24]）。这些发现都将导致涡动力学的丰富和继续发展。
5 涡量面理论的意义
为检验第4节的动力学定义，对计算/实验流体力学（CFD/EFD）方法有个基本需求———能跟踪大雷诺数下薄剪切层的拉格朗日演化与扩散，及涡量管的拉伸、折转与相互作用。更一般地说，至少对低马赫数流，需要CFD/EFD具有模拟复杂流场中涡量的产生及其运动学-动力学演化和相互作用的能力。对经典旋涡案例来说，这个需求已可满足，但尚无系统的验证。而对湍流中的复杂结构，还需要CFD/EFD有进一步的重大发展。
前面提到的面涡即切向间断面这个简化概念，在有限的大雷诺数下变成剪切层或涡层，而一个涡层无非是一族无限多个涡量面，它们的动力学（运动学和动理学）演化完整地刻画着涡层的演化，包括其可能的卷绕成旋涡的过程。在CFD中，一个涡量连续分布的剪切层总要离散化为一族强度足够小的面涡，每个这样的面涡就是一个涡量面，其中最需要关注的是代表稳定性理论判定的最早失稳的临界层的那个涡量面。因此，对典型涡量面演化过程的DNS研究，开辟了对各种相干结构给出统一物理解释的重要方向。这个研究方向始于Yang&Pullin(2010)[25]。
这个方向近10年来取得了可观的进展。Zhao等（2016,2018)[26,27]初步算出了边界层中精心选定的初始平直涡量面如何在适当扰动下失稳而拱起来并逐步卷成准周期发卡涡列的过程；槽道流中上下壁面形成的发卡涡如何相互作用而重联；发卡涡的尾部如何形成“条带”结构；多个发卡涡的后期演化如何导致“猝发”和“湍斑”，等等。这就初步实现了原来的各种碎片化认识在理性基础上的统一。的确，由于对剪切层和旋涡相互关系有清晰的表述，上述旋涡的动力学定义在追踪涡量面演化的过程中得到了很好的印证。发卡涡就是开放涡量面卷成的有限长旋涡。它的头和腿是卷紧的旋涡，脚是开放涡量面。
我们说“初步”，因为只跟踪了极少数涡量面，其行为未必能代表该涡量面所属的整个剪切层（剪切层的边缘也是模糊的）；而且黏性扩散使涡量面不能保持为同一物质面，数值上准确追踪这样的涡量面非常困难。尽管现在误差可控制在5%以下，但它在复杂湍流中或许有时足以使追踪偏离到另一个涡量面上去了。由于这些局限，目前涡量面理论的计算结果尚未非常清晰地解读湍流中复杂涡结构的相互作用图景。但无论如何，涡量面理论开辟了一个回归到涡量动力学的正确方向，其今后发展应注重描述代表剪切层的一族涡量面。事实上，Zhao等（2018)[28]已作出一定尝试，即利用体渲染的方法绘制一族涡量面来显示转捩中复杂的湍斑结构，如图8所示。
图8 用一族涡量面算出的发卡涡列的形成与演化[28]
6 结语
旋涡是流体运动横过程的典型结构，其发生和演化过程及其对整个流动的影响，应当根据而且只能根据涡动力学（关于横过程及其与纵过程相互作用的理论）来解释。其定义也必能在涡动力学框架内改进和完善。本文提出的动力学涡定义方案是这个信念下的一个尝试，和CFD中涡量面理论的概念与方向相契合，也应对涡量面理论的未来发展有所启示。
参考文献
[1] WU J Z,MA H Y,ZHOU M D.Vortical flows[M].Springer,2015.
[2] TRUESDELL C.The kinematics of vorticity[M].Bloomington:Indiana University Press,1954.
[3] PRANDTL L, TIETJENS O G. Applied Hydro-and Aeromechanics[M].McGraw-Hill,1934.
[4] LANCHESTER F W.Aerodynamics[M].London:Constable,1907.
[5] WU J Z.Vortex definition and“vortex criteria”[J].Science China Physics,Mechanics&Astronomy,2018,61(2):024731.
[6] FUENTES O U V.On the topology of vortex lines and tubes[J].Journal of Fluid Mechanics,2007,584:147-156.
[7] Kelvin(Sir W.Thomson).On vortex motion[J].Trans R Soc Edinb,1869,25:217-260.
[8] SAFFMAN P G,BAKER G R.Vortex interactions[J].Annual Review of Fluid Mechanics,1979,11(1):95-121.
[9] SAFFMAN P G.Vortex dynamics[M].Cambridge University Press,1992.
[10] WU J Z,MA H Y,ZHOU M D.Vorticity and vortex dynamics[M].Berlim Heidelberg,2006:151.
[11]KCHEMANN D.Report on the I.U.T.A.M.symposium on concentrated vortex motions influids[J].Journal of Fluid Mechanics,1965,21(1):1.
[12] MOFFATT H K,KIDA S,OHKITANI K.Stretched vorticesthe sinews of turbulence;large-Reynolds-number asymptotics[J].J Fluid Mech,1994,259:241-264.
[13]EPPS B.Review of vortex identification methods[R].AIAA2017-0989.
[14] WU J Z,XIONG A K,YANG Y T.Axial stretching and vortex definition[J].Physics of Fluids,2005,17(3):038108.
[15]KOLR V，STEK J.Stretching response of Rortex and other vortex-identification schemes[J].AIP Advances,2019,9(10):105025.
[16]GAO Y S,LIU C Q.Rortex and comparison with eigenvaluebased vortex identification criteria[J].Physics of Fluids,2018,30(8):085107.
[17]KADEN H.Aufwicklung einer unstabilen Unstetigkeitsflche[J].Ingenieur-Archiv,1931,2(2):140-168.
[18]BETZ A.Wie entsteht ein Wirbel in einer wenig zhen Flüssigkeit?[J].Naturwissenschaften,1950,37(9):193-196.
[19]KRASNY R.Computation of vortex sheet roll-up in the Trefftz plane[J].Journal of Fluid Mechanics,1987,184:123-155.
[20] 今井明子，小松研吾，等．魏俊霞（译）．天气预报的科学[J]．科学世界，2019,10:49.
[21]CHORIN A J.Vorticity and turbulence[M].New York,NY:Springer New York,1994.
[22] WU X S.Nonlinear theories for shear flow instabilities:physical insights and practical implications[J].Annual Review of Fluid Mechanics,2019,51(1):451-485.
[23] MOORE D W,SAFFMAN P G.Axial flow in laminar trailing vortices[J].Proceedings of the Royal Society A:Mathematical,Physical and Engineering Sciences,1973,333:491-508.
[24]LUNDGREN T S.Strained spiral vortex model for turbulent fine structure[J].Physics of Fluids,1982,25(12):2193.
[25]YAO J,HUSSAIN F.A physical model of turbulence cascade via vortex reconnection sequence and avalanche[J].Journal of Fluid Mechanics,2020,883:A51.
[26]YNG Y,PULLIN D I.On Lagrangian and vortex-surface fields for flows with Taylor-Green and Kida-Pelz initial conditions[J].Journal of Fluid Mechanics,2010,661:446-481.
[27]ZHAO Y M,YANG Y,CHEN S Y.Vortex reconnection in the late transition in channel flow[J].Journal of Fluid Mechanics,2016,802:R4.
[28]ZHAO Y M,XIONG S Y,YANG Y,et al.Sinuous distortion of vortex surfaces in the lateral growth of turbulent spots[J].Physical Review Fluids,2018,3(7):074701.
