郭春雨， 薛嶸， 胡健， 王戀舟

(1.哈爾濱工程大學(xué) 船舶工程學(xué)院，黑龍江 哈爾濱 150001； 2.西南交通大學(xué) 力學(xué)與工程學(xué)院，四川 成都 611756)

梢渦在有限展長的翼型繞流現(xiàn)象中普遍存在。當(dāng)流動(dòng)流過有攻角的有限長水翼的端部，水翼上下表面的壓力差會(huì)驅(qū)使流動(dòng)從下表面的高壓側(cè)流向上表面的低壓側(cè)，從而形成三維的旋渦，我們把它稱為梢渦[1]。梢渦具有強(qiáng)烈的旋轉(zhuǎn)效應(yīng)，使得渦核中心的壓力值顯著降低。在水動(dòng)力學(xué)領(lǐng)域，梢渦引起的壓降會(huì)驅(qū)使水中的氣核卷入渦核內(nèi)部，產(chǎn)生空化。梢渦空化會(huì)降低水翼結(jié)構(gòu)的水動(dòng)力性能，并會(huì)伴隨出現(xiàn)強(qiáng)烈的噪聲。為了避免梢渦空化帶來的危害，需要對(duì)它的形成機(jī)理進(jìn)行分析。

學(xué)者們進(jìn)行了大量的試驗(yàn)來研究梢渦的細(xì)微結(jié)構(gòu)。Giuni等[2-3]通過煙跡試驗(yàn)，得出梢渦是由起始于壓力面的主渦和2個(gè)分別起始于吸力面和壓力面的次級(jí)渦演化形成的。徐良浩等[4]使用2D-PIV和3D-PIV對(duì)梢渦的二維和三維速度場(chǎng)進(jìn)行了測(cè)量，獲得了不同工況下的渦核位置、尺寸及速度分布，得出了渦核在三維空間內(nèi)的擺動(dòng)特性。Maine等[5]通過對(duì)3種不同翼型剖面的實(shí)驗(yàn)觀測(cè)，研究了水翼的縱傾和側(cè)斜對(duì)梢渦形成的影響。梢渦試驗(yàn)往往伴隨著空化初生現(xiàn)象。Arndt等[6]通過LDV試驗(yàn)指出，梢渦空化發(fā)生的必要條件是分布在水中的氣核進(jìn)入渦核中心并發(fā)展為足夠大的氣泡。在這個(gè)基礎(chǔ)上，Arndt等[7]進(jìn)一步研究了壓力變化和水中含氣量對(duì)梢渦空化的影響。他們得出在含氣量多的水中，梢渦空化多發(fā)生于水翼梢部；而在含氣量較少的水中，梢渦空化的起始位置在梢部流動(dòng)后端。

隨著計(jì)算機(jī)運(yùn)算能力的逐步發(fā)展，采用各種湍流模型來模擬梢渦逐步成為可能。相比各種試驗(yàn)手段，數(shù)值模擬可以觀察到更加細(xì)微的渦結(jié)構(gòu)，從而更好地去觀察梢渦的演化，進(jìn)而了解梢渦空化初生的規(guī)律[8]。盡管如此，由于梢渦尺寸很小，壓力梯度很大，模擬梢渦是一件具有挑戰(zhàn)性的工作[9]。Decaix等[10]采用RANS模型對(duì)截面形狀為NACA 0009的水翼梢渦進(jìn)行了數(shù)值模擬，數(shù)值結(jié)果表明RANS方法可以準(zhǔn)確模擬出梢渦的渦核位置、平均速度場(chǎng)等的時(shí)均特性。蒲汲君等[11]使用k-w、DES和LES模型分別對(duì)水翼的梢渦流場(chǎng)進(jìn)行了數(shù)值模擬，并使用氣泡靜力平衡方程計(jì)算了空化初生現(xiàn)象。結(jié)果表明LES模型可以很好地模擬梢渦處流場(chǎng)，此外應(yīng)用氣泡靜力方程可以準(zhǔn)確計(jì)算初始梢渦空泡數(shù)。當(dāng)涉及空化初生問題時(shí)，可以耦合球形空泡模型來模擬氣泡運(yùn)動(dòng)[12]。Rayleigh[13]建立了球形氣泡模型，用于描述球形氣泡在變化的壓力場(chǎng)中的半徑變化。許多研究致力于分析氣泡及顆粒在水中的運(yùn)動(dòng)，Maxey等[14]給出了一種描述剛性球體的運(yùn)動(dòng)方程，可以很好地描述球形空泡的運(yùn)動(dòng)規(guī)律。Park等[15]利用他的計(jì)算方程，計(jì)算了氣泡在水翼梢部的運(yùn)動(dòng)。

本文以剖面為NACA 0015的三維水翼為研究對(duì)象，采用大渦模擬(large eddy simulation, LES)方法來模擬水翼梢渦。將計(jì)算結(jié)果與試驗(yàn)值[16]進(jìn)行對(duì)比，滿足精度要求之后，在這個(gè)基礎(chǔ)上開展基于離散元 (discrete element method, DEM)模型的氣核在水中運(yùn)動(dòng)研究，分析氣核在梢渦作用下的演化規(guī)律。整個(gè)計(jì)算工作在數(shù)值模擬軟件STAR-CCM+下完成。

1 幾何模型與數(shù)值模型及網(wǎng)格劃分

1.1 幾何模型

幾何模型與文獻(xiàn)[17]中試驗(yàn)所用翼型保持一致。如圖1所示，翼型截面形狀為NACA 0015，弦長C=0.15 m，展弦比6.6，攻角12°，端部裝有圓形導(dǎo)流罩。計(jì)算域布置以及邊界條件如圖2所示。為了節(jié)省計(jì)算資源，將水翼的展長選為3.3倍弦長，遠(yuǎn)離梢部的一端貼近壁面，并將壁面設(shè)置為對(duì)稱邊界條件，從而節(jié)省一半的網(wǎng)格數(shù)量。由此，該計(jì)算域的寬為4.5C，高為4.5C，水翼前端距速度入口距離為3.3C，后端距壓力出口距離9C，其余壁面均設(shè)置為無滑移壁面條件。計(jì)算環(huán)境為25 ℃的水，此時(shí)密度998 kg/m3，動(dòng)力粘度0.000 889 9 kg/(m·s)，流速10 m/s，從而雷諾數(shù)Re=1.5×106。

1.2 湍流模型

本研究采用大渦模擬(large eddy simulation, LES)求解流體域中任何位置的空間湍流結(jié)構(gòu)。網(wǎng)格尺寸超出選用尺寸的大尺度旋渦進(jìn)行直接求解，小于選用網(wǎng)格尺寸的細(xì)微旋渦結(jié)構(gòu)采用亞格子進(jìn)行模擬。從而在復(fù)雜的流動(dòng)模擬中得到細(xì)微旋渦結(jié)構(gòu)的流動(dòng)圖像。

圖1 水翼模型及導(dǎo)流罩Fig.1 Hydrofoil model with round end cap

圖2 計(jì)算域布置條件Fig.2 Compute region conditions

(1)

式中φ表示速度分量、壓力、能量或組分濃度。

將分解的求解變量插入納維-斯托克斯方程可得到已濾波物理量的方程。濾波后的質(zhì)量和動(dòng)量的傳輸方程可以寫為：

(2)

(3)

湍流應(yīng)力張量現(xiàn)在表示為亞網(wǎng)格尺度應(yīng)力。這些應(yīng)力產(chǎn)生自較大的已求解渦流與較小的未解析渦流之間的相互作用，使用Boussinesq近似進(jìn)行建模，為：

(4)

(5)

亞網(wǎng)格尺度湍流粘度μt必須由能夠解釋小渦流對(duì)已求解流體的作用的亞網(wǎng)格尺度模型來描述。采用WALE亞網(wǎng)格尺度模型，使用混合長度假設(shè)來對(duì)亞網(wǎng)格尺度應(yīng)力建模：

μt=ρΔ2Sw

(6)

式中：Δ為長度尺度或網(wǎng)格濾波器寬度，使用網(wǎng)格單元體積定義如下：

Δ=min(κd,CwV1/3)

(7)

式中：Cw為模型系數(shù)，典型值介于0.5(對(duì)于均勻各向同性衰變湍流)～0.325(對(duì)于通道流)之間;κ為馮卡門常數(shù)；Sw的定義見文獻(xiàn)[17]。

1.3 DEM模型

本研究采用DEM模型對(duì)氣核的運(yùn)動(dòng)特性進(jìn)行模擬。離散元方法(DEM)是一種基于拉格朗日方法的數(shù)值方法，通常用來模擬大量相互作用的離散對(duì)象(通常是顆粒)如何運(yùn)動(dòng)。根據(jù)格雷斯等[18]的研究，氣核直徑通常在30～300 μm，在這一尺度下，氣核在運(yùn)動(dòng)過程中可以維持較為穩(wěn)定的球形外形，從而在本研究中，將氣核簡(jiǎn)化為不可形變的球形顆粒，應(yīng)用DEM方法對(duì)氣核進(jìn)行建模。

DEM方法以牛頓第二定律作為控制方程：

(8)

(9)

式中：mi、vi和ωi分別表示顆粒i的質(zhì)量、速度和角速度；Ii為顆粒i的轉(zhuǎn)動(dòng)慣量，對(duì)于球形粒子，Ii=2miRi2/5，Ri為顆粒i的半徑；Fg=mig為顆粒i的重力；,Ffluid為流體對(duì)顆粒i的作用力(包括阻力、升力、附加質(zhì)量力和浮力等)；Fij為顆粒i與顆粒j或墻壁之間的碰撞力以及其他作用在粒子上的非接觸力；Tij為接觸力矩，表示除顆粒重力以外的接觸力在粒子上產(chǎn)生的扭矩。

STAR-CCM+中的DEM模型在拉格朗日架構(gòu)內(nèi)實(shí)現(xiàn)，并使用拉格朗日相模型對(duì)粒子進(jìn)行定義。本文使用DEM方法對(duì)氣核進(jìn)行模擬，并與CFD結(jié)合使用，從而求解由顆粒組成的離散相的氣核在流體介質(zhì)中的運(yùn)動(dòng)。其中，流場(chǎng)的求解由大渦模擬完成；氣核的運(yùn)動(dòng)基于拉格朗日方法,由式(8)、(9)進(jìn)行求解。STAR-CCM+將對(duì)DEM粒子和流體之間的相互作用進(jìn)行建模，并將LES計(jì)算得到的流場(chǎng)中各網(wǎng)格處的力代入到式(8)、(9)中，對(duì)氣核的運(yùn)動(dòng)進(jìn)行求解。

1.4 網(wǎng)格劃分

采用Star-CCM+對(duì)水翼模型進(jìn)行面網(wǎng)格重構(gòu)，形成表面三角化良好的高質(zhì)量面網(wǎng)格，以面網(wǎng)格為基礎(chǔ)生成帶有邊界層和切割體網(wǎng)格的非結(jié)構(gòu)體網(wǎng)格。邊界層共劃分12層棱柱層網(wǎng)格，并將y+值控制在1以下，從而滿足LES模擬的計(jì)算要求。對(duì)水翼表面和梢渦區(qū)域的體網(wǎng)格進(jìn)行加密處理，最終確認(rèn)計(jì)算域網(wǎng)格總數(shù)約為780萬，網(wǎng)格細(xì)節(jié)如圖3所示。

圖3 全局網(wǎng)格、加密區(qū)和邊界層Fig.3 Overall mesh condition, refining region and boundary layer mesh

2 計(jì)算結(jié)果分析

2.1 數(shù)值驗(yàn)證

為驗(yàn)證所采用的數(shù)值模型是否滿足計(jì)算要求，以數(shù)值模擬結(jié)果與文獻(xiàn)[16]中的相關(guān)試驗(yàn)工況的結(jié)果進(jìn)行對(duì)比，來驗(yàn)證模擬結(jié)果的準(zhǔn)確性。監(jiān)測(cè)水翼表面受到的力在x和y方向上的分解值，即阻力和升力，并將升力換算為力系數(shù)并與試驗(yàn)值進(jìn)行對(duì)比?？梢缘贸霰敬螖?shù)值模擬得到的升力系數(shù)為1.18，對(duì)比文獻(xiàn)[16]中對(duì)應(yīng)的試驗(yàn)結(jié)果的升力系數(shù)為1.24，相對(duì)誤差在4.8%以內(nèi)，可以認(rèn)為本文的計(jì)算結(jié)果與文獻(xiàn)中的結(jié)果較為吻合，從而本文采用的數(shù)值模型是可靠的。

圖4為梢渦附近流場(chǎng)無因次垂向速度沿弦長方向分布結(jié)果與試驗(yàn)結(jié)果的對(duì)比。在水翼隨邊后端0.2倍弦長距離的位置布置一排監(jiān)測(cè)點(diǎn)，用于監(jiān)測(cè)梢渦附近的垂向速度值，如圖5所示。本研究采用的大渦模擬方法是一種瞬態(tài)的模擬方法，因此各監(jiān)測(cè)點(diǎn)需要將監(jiān)測(cè)一段時(shí)間后的速度值進(jìn)行平均化，再將平均化后的垂向速度值UZ除以來流速度U∞進(jìn)行無因次化，由于梢渦軌跡近似于X軸方向，從而認(rèn)為UZ/U∞的值作為梢渦的無因次切向速度。從圖5中可以看出，數(shù)值模擬對(duì)渦核位置附近的切向速度值捕捉與試驗(yàn)值較為相符，除了在峰值和谷值處有一定的偏差。這是因?yàn)闊o量綱化速度的極值位于梢渦渦核位置，由于渦核高強(qiáng)度的旋轉(zhuǎn)帶來了切向速度的極值。試驗(yàn)中使用的渦量計(jì)布置位置為一系列等間距點(diǎn)，在峰值附近，布置點(diǎn)較少而漏掉了極值；而模擬中在渦核位置附近布置了更多的監(jiān)測(cè)點(diǎn)，從而監(jiān)測(cè)到渦核的極值。通過對(duì)數(shù)值模擬中速度場(chǎng)的驗(yàn)證，可以認(rèn)為本研究對(duì)于流場(chǎng)的模擬具有較高的精確度。

圖4 水翼后端0.2倍弦長處梢渦附近垂向速度對(duì)比Fig.4 Comparison of vertical velocity across the tip vortex at r=0.2C

圖5 水翼后端0.2倍弦長處垂向速度監(jiān)測(cè)位置Fig.5 The location of vertical velocity monitoring at r=0.2C

2.2 結(jié)果分析

2.2.1 流場(chǎng)結(jié)構(gòu)

本節(jié)從梢渦結(jié)構(gòu)方面來分析梢渦流場(chǎng)數(shù)值模擬的結(jié)果。使用Q判據(jù)來可視化渦結(jié)構(gòu)，Q判據(jù)的計(jì)算公式為：

(10)

式中:Ω為渦張量;S為應(yīng)變率張量。

圖6為Q=10 000 s-2的瞬時(shí)梢渦等值面，并采用渦量幅值進(jìn)行著色。為了展現(xiàn)整個(gè)梢渦位置渦量場(chǎng)的細(xì)節(jié)，分別從水翼側(cè)面和后方視角展示了渦量場(chǎng)的整體與局部。由圖6(a)展現(xiàn)的整體渦結(jié)構(gòu)可知，本文采用的LES模型結(jié)合水翼后端位置網(wǎng)格加密，可以捕捉到精細(xì)的梢渦結(jié)構(gòu)，且梢渦主渦的強(qiáng)度值在1 600 s-1附近可以保持在5倍弦長之上的距離，這同Bailey等[19]的試驗(yàn)結(jié)果相符。從等值面渦量場(chǎng)可得，本研究對(duì)于梢渦區(qū)域的網(wǎng)格細(xì)化程度滿足分析要求。

圖6(b)展示了近梢部端渦量場(chǎng)的渦結(jié)構(gòu)。從圖中可以看出，梢部渦系由2部分組成，序號(hào)為1的渦來源于水翼壓力面梢部的位置邊界層發(fā)生分離，進(jìn)而產(chǎn)生渦強(qiáng)很高的梢渦；序號(hào)為2的渦來源于水翼吸力面靠近梢部的位置，這部分的流動(dòng)在從水翼隨邊分離后，在梢渦的強(qiáng)旋轉(zhuǎn)效應(yīng)下形成了包裹在梢渦外邊的渦系。

圖7為不同截面上的渦量云圖，截面的選取參考圖8，選取水翼隨邊為起始位置，每隔0.5倍弦長距離共作6個(gè)截面，形成圖7的截面。從渦量云圖可以觀察到，梢渦位置具有多重渦結(jié)構(gòu)組成的渦系。壓力面形成一個(gè)主渦和一個(gè)二次渦，在向下游遷移到1.5倍弦長附近發(fā)生融合；吸力面發(fā)生流動(dòng)分離并形成多個(gè)較弱的渦心，在向下游發(fā)展的過程中渦量逐漸耗散。

圖6 Q=10 000 s-2等值面的渦量場(chǎng)及細(xì)節(jié)展示Fig.6 Vortex field and detail display of iso-surface with Q=10 000 s-2

2.2.2 氣核運(yùn)動(dòng)特性分析

本節(jié)在得到穩(wěn)定發(fā)展的渦量場(chǎng)的基礎(chǔ)上，向流場(chǎng)中加入DEM粒子所表示的氣核，來研究水翼梢部氣核的空間演化過程。

拉格朗日相的DEM粒子需要通過噴射器加入到歐拉相的流場(chǎng)中進(jìn)行求解。如圖9為DEM粒子的噴射點(diǎn)位置，在水翼的端部，沿弦長方向每隔1/10弦長設(shè)置一排噴射點(diǎn)，這些噴射點(diǎn)沿圓形導(dǎo)流罩長度方向間距為1 mm。隨后在拉格朗日相下設(shè)置氣核的物理屬性，在本研究中，氣核的半徑為100 μm，密度為氣體密度。如此設(shè)置后在流場(chǎng)中噴射點(diǎn)的位置各出現(xiàn)1個(gè)氣核，隨后氣核將在DEM方程的控制下，基于受到流場(chǎng)中的各種力隨梢渦向水翼后端遷移。

圖7 水翼隨邊后端渦量云圖截面分布Fig.7 The distribution of vorticity contour section behind hydrofoil′s trailing edge

圖8 不同截面上渦量云圖Fig.8 Vorticity contour in typical sections

圖9 DEM粒子噴射點(diǎn)位置Fig.9 DEM injector positions distribution

如圖10所示，選取噴射點(diǎn)位于端部中軸線并距端部垂直距離為1 mm處的氣核流跡進(jìn)行分析。圖中用速度場(chǎng)對(duì)氣核流跡進(jìn)行著色，深色線條為和氣核初始位置相同的流線。從流線軌跡可知，初始流動(dòng)位置越靠近水翼前緣，在向下游遷移的過程中會(huì)出現(xiàn)遠(yuǎn)離梢渦核心的趨勢(shì)；而初始位置接近水翼后緣的流動(dòng)會(huì)被梢渦捕獲。對(duì)比經(jīng)過各氣核初始位置的流線，可以對(duì)這一現(xiàn)象進(jìn)行解釋。由于梢渦在水翼壓力面后端開始產(chǎn)生，因此水翼端部后半端的流動(dòng)受旋渦影響較大，從而這些位置的流線緊緊跟隨梢渦運(yùn)動(dòng)；水翼端部前半段產(chǎn)生的流線表明，經(jīng)過這部分的流動(dòng)會(huì)沿水翼展向遷移，在隨邊與水翼分離后，受到圖6(b)中渦系2的影響，隨著包裹在渦核外部的渦系旋轉(zhuǎn)運(yùn)動(dòng)。氣核具有一定的流線跟隨性，從而伴隨著2種不同的流線運(yùn)動(dòng)狀態(tài)，氣核也具有2種不同的運(yùn)動(dòng)狀態(tài)。

圖10 初始位置距端部1 mm處氣核流跡及流線對(duì)比Fig.10 Comparison between nuclei track and streamline with initial location from round end=1 mm

氣核的流跡雖然表現(xiàn)出一定的流線跟隨性，但是與流線還是有一些區(qū)別。從圖6(b)細(xì)節(jié)可以看出，氣核的遷移軌跡極易受到渦強(qiáng)的影響。在水翼壓力面梢渦的主渦剛形成的位置產(chǎn)生的氣核，會(huì)迅速遷移到旋渦中心，跟隨渦核向下流旋轉(zhuǎn)運(yùn)動(dòng)；在水翼前端部產(chǎn)生的氣核，在隨流線遷移到與水翼分離后，遷移軌跡位于梢渦外部包裹的渦系內(nèi)?？梢缘贸?，氣核先隨流線向后遷移，當(dāng)軌跡附近出現(xiàn)強(qiáng)渦量區(qū)域后，氣核會(huì)向渦強(qiáng)較強(qiáng)的位置進(jìn)行遷移，從而不再跟隨流線。

接下來選取部分氣核，對(duì)它們流跡的特征進(jìn)行分析。如圖11分別為1、2、5、7、9號(hào)的噴射器發(fā)射的氣核流跡，在向下游發(fā)展的過程中垂向速度的變化曲線。對(duì)遷移距離除以弦長得到L/C，垂向速度除以來流速度得到V/U進(jìn)行無量綱化處理。由于水平長度7C以后的網(wǎng)格位于梢渦加密區(qū)外，因此不再取之后的速度進(jìn)行分析。從曲線中可以看出，為1號(hào)和2號(hào)的氣核由于流跡位于旋渦核心外部，在遷移的過程中速度振蕩頻率較小。而其余氣核流跡位于渦核內(nèi)部，可以按照振蕩特性將軌跡分為3個(gè)階段。如圖11(c)所示，第1階段是從氣核剛進(jìn)入流場(chǎng)到遷移到水翼尾端位置，這一階段氣核開始進(jìn)入剛形成的梢渦內(nèi)部，垂向速度的變化表現(xiàn)出較大的隨機(jī)性。這是由于這一階段旋渦剛開始形成，而氣核不再跟隨流線運(yùn)動(dòng)并開始進(jìn)入渦核內(nèi)部，造成了無規(guī)律的垂向速度變化。第2階段是氣核從到水翼尾端位置遷移到距尾端水平距離為4倍弦長處。這一階段渦核具有相對(duì)較大的渦強(qiáng)，氣核隨渦核向下游段遷移過程中具有較為穩(wěn)定的振蕩頻率和幅值。第3階段位于水平距離4倍弦長之后，在這一階段，渦核的強(qiáng)度逐步下降，而氣核的垂向速度變化頻率也隨之變緩慢。

圖11 部分氣核流跡上的垂向速度分布Fig.11 Vertical velocity distribution in typical nuclei track

圖12中，(a)～(e)分別為1、2、5、7、9號(hào)的噴射器發(fā)射出的氣核運(yùn)動(dòng)過程中相對(duì)壓力的變化曲線。對(duì)遷移距離除以弦長得到L/C，對(duì)計(jì)算得到的相對(duì)壓力轉(zhuǎn)化為相對(duì)壓力系數(shù)Cp，完成無量綱化處理。壓力變化曲線與速度變化曲線有一定的相似之處，位于梢渦核心位置處和梢渦核心以外區(qū)域的流跡在壓力變化上具有2種明顯不同的趨勢(shì)。1號(hào)和2號(hào)的氣核在遷移過程中，壓力首先經(jīng)歷了一段明顯的上升，在遷移到水翼后端1倍弦長的位置處相對(duì)壓力達(dá)到峰值；隨后下降到-0.3左右。而5、7和9號(hào)的氣核在遷移到水翼后端1倍弦長的位置處壓力先急劇下降，隨后再急劇上升，最后平穩(wěn)上升。從氣核的遷移軌跡來看，1號(hào)和2號(hào)的氣核在剛開始遷移的過程中隨流線向展長的方向進(jìn)行遷移，因而經(jīng)歷了一次壓力上升；隨后被梢渦外端渦系捕獲后，壓力再次下降并穩(wěn)定。而5、7、9號(hào)的氣核剛開始遷移時(shí)被渦核所捕獲，故而壓力急劇下降；隨后在隨渦核向后遷移的過程中，隨著渦強(qiáng)降低，壓力也隨之上升；而在渦強(qiáng)下降到一定程度后，下降趨勢(shì)逐步變慢，此時(shí)壓力上升趨勢(shì)也開始放緩。

圖12 部分氣核流跡上的相對(duì)壓力分布Fig.12 Relative pressure distribution in typical nuclei track

3 結(jié)論

1)采用LES模型對(duì)帶導(dǎo)流罩的三維水翼受到的升阻力、無因次速度分布進(jìn)行計(jì)算，與文獻(xiàn)中數(shù)據(jù)符合較好，計(jì)算結(jié)果具有良好的可信度。

2)帶導(dǎo)流罩的圓形導(dǎo)邊水翼的梢渦流場(chǎng)有著獨(dú)特的結(jié)構(gòu)。在壓力面出現(xiàn)主渦和二次渦，在吸力面出現(xiàn)一系列小渦。主渦和二次渦在發(fā)展的過程中合并，形成梢渦核心，而吸力面形成的渦系在向后遷移的過程中，包裹在核心外部，形成總體的渦系。

3)氣核在遷移過程中，在水翼首部出現(xiàn)的氣核的遷移軌跡位于梢渦核心外，而水翼首部后端的氣核遷移軌跡位于梢渦核心內(nèi)。此外，渦核會(huì)對(duì)氣核的遷移產(chǎn)生很大的影響，使其偏離流線的方向；氣核被渦核捕獲后，垂向速度會(huì)有較大的振蕩值，受到的壓力也會(huì)產(chǎn)生顯著的變化。