黃曉婷, 孫鵬楠, 呂鴻冠, 殷曉瑞, 董嘉徐

(1.中山大學(xué) 海洋工程與技術(shù)學(xué)院, 廣東 珠海 519082;2.南方海洋科學(xué)與工程廣東省實(shí)驗(yàn)室(珠海), 廣東 珠海 519082)

近年來,微型飛行器或潛水器的發(fā)展受到廣泛關(guān)注,飛行器或潛水器性能的改進(jìn)離不開對(duì)翼型繞流問題的精確模擬。翼型表面流動(dòng)如分離、附著等現(xiàn)象對(duì)翼型的動(dòng)力性能具有重要影響。為保證在不同工況下飛行器或潛水器運(yùn)行的可靠性,需要了解翼型繞流的流動(dòng)特性和機(jī)理。這就要求對(duì)翼型繞流現(xiàn)象進(jìn)行深入研究和分析。

翼型的幾何形狀對(duì)改善力學(xué)性能具有重要影響[1]?？梢园l(fā)現(xiàn),目前對(duì)于非對(duì)稱翼型的繞流現(xiàn)象的研究較少。其中,Anyoji等[2]經(jīng)過研究發(fā)現(xiàn),非對(duì)稱Ishii弧形翼型由于后緣附近的弧度,其壓力阻力要遠(yuǎn)遠(yuǎn)小于其他對(duì)稱翼型,在同等工況下,可獲得比對(duì)稱翼型更高的升力。因此,展開對(duì)Ishii非對(duì)稱翼型更深入研究,準(zhǔn)確獲取流場(chǎng)信息及其運(yùn)動(dòng)特性,可為未來設(shè)計(jì)出良好動(dòng)力性能的新型翼型提供新思路。

現(xiàn)有翼型繞流研究中,主流研究方法中主要有試驗(yàn)法及數(shù)值模擬法。針對(duì)水翼繞流,季斌等[3]從試驗(yàn)和數(shù)值模擬角度綜述水翼水動(dòng)力特性。在飛行氣翼繞流探究中,Anyoji等[2]采用風(fēng)洞試驗(yàn)探究及驗(yàn)證Ishii翼型氣動(dòng)性能。然而在試驗(yàn)法中,試驗(yàn)條件復(fù)現(xiàn)存在一定挑戰(zhàn),且邊界層湍流細(xì)節(jié)捕捉對(duì)儀器的精密性要求苛刻。此外,試驗(yàn)法的成本高昂,且對(duì)多工況無法同時(shí)進(jìn)行研究。因此現(xiàn)有研究大多采取數(shù)值模擬法展開研究。數(shù)值模擬法中，發(fā)展較為成熟的歐拉網(wǎng)格法較為常見。但對(duì)于網(wǎng)格法而言,翼型后緣狹窄尖角邊界層較薄,網(wǎng)格劃分邊界層存在一定難度[4],在網(wǎng)格離散中需要較高的網(wǎng)格構(gòu)建技術(shù)來滿足復(fù)雜的邊界條件。值得注意的是,近年來發(fā)展的無網(wǎng)格法對(duì)邊界條件、自由液面及大變形等問題[5]處理具有固有優(yōu)勢(shì)。

SPH方法作為無網(wǎng)格法中應(yīng)用廣泛的方法,在流體或固體力學(xué)問題處理中逐漸發(fā)揮重要作用[6]且對(duì)于自由液面變形動(dòng)邊界問題的處理上SPH方法具有顯著優(yōu)勢(shì)[7]。目前,在水動(dòng)力學(xué)領(lǐng)域上得到較多應(yīng)用[8-9]。Sun等[10]采用δ+-SPH模型對(duì)NACA0012翼型固定及自由運(yùn)動(dòng)狀態(tài)進(jìn)行模擬,分析其水動(dòng)力特性。Huang等[11]采用SPH方法對(duì)翼型繞流的運(yùn)動(dòng)性能進(jìn)行探究;在空氣動(dòng)力學(xué)領(lǐng)域,Shadloo等[12]驗(yàn)證了WCSPH可以實(shí)現(xiàn)低雷諾數(shù)下不同攻角翼型繞流計(jì)算;Huang等[13]采用KGF-SPH模型,在低雷諾數(shù)下對(duì)氣翼的升阻力系數(shù)進(jìn)行監(jiān)測(cè),驗(yàn)證該模型的準(zhǔn)確性?？傻贸?SPH方法在翼型繞流問題研究上具有一定的優(yōu)勢(shì)及特色。但由于其本身的精度及穩(wěn)定性問題,SPH方法在翼型氣動(dòng)性上的研究受到一定限制,翼型繞流中邊界層內(nèi)部流場(chǎng)特征獲取仍需進(jìn)一步的方法改進(jìn)和技術(shù)討論。

因此,本文克服以往SPH方法難以模擬黏性邊界層流動(dòng)的不足,采用δ+-SPH模型,實(shí)現(xiàn)黏性流動(dòng)下翼型繞流的模擬,展開翼型特性分析。在SPH計(jì)算中,高雷諾數(shù)下翼型繞流模擬存在三大難點(diǎn),主要應(yīng)用以下方法進(jìn)行解決:①針對(duì)流體繞過翼型時(shí),由于負(fù)壓產(chǎn)生的張力不穩(wěn)定現(xiàn)象,采用張力不穩(wěn)定性控制(tensile instability control,TIC)技術(shù),實(shí)現(xiàn)Ishii翼型黏性繞流的壓力場(chǎng)和速度場(chǎng)等穩(wěn)定數(shù)值模擬。②為了避免粒子的拉格朗日結(jié)構(gòu)引起的計(jì)算精度下降現(xiàn)象,在每個(gè)時(shí)間步初始時(shí)刻,施加粒子修正技術(shù)(particle shifting technique,PST),確保粒子的均勻分布。③為實(shí)現(xiàn)邊界層黏性力的精確計(jì)算,采用自適應(yīng)粒子細(xì)化(adaptive particle refinement,APR)技術(shù),減少粒子數(shù)量,提高計(jì)算效率,實(shí)現(xiàn)粒子高分辨率同時(shí)確保計(jì)算精度。

1 δ+-SPH模型

在δ+-SPH模型中,考慮流體的真實(shí)物理黏性,在動(dòng)量方程中施加黏性項(xiàng),以精確模擬邊界層的流動(dòng)。再有,在連續(xù)方程中添加密度耗散項(xiàng)[14],以避免高頻的壓力振蕩。此外,為了避免傳統(tǒng)SPH模型中由于粒子的拉格朗日結(jié)構(gòu)引起的計(jì)算精度降低現(xiàn)象,δ+-SPH模型針對(duì)粒子進(jìn)行位移修正,引入粒子位移修正技術(shù)(particle shifting technique,PST),確保在計(jì)算過程粒子分布的均勻性。因此,基于流體弱可壓假設(shè)下,δ+-SPH模型[15]最終寫為

(1)

式中：i和j表示支持域內(nèi)目標(biāo)粒子及其對(duì)應(yīng)的鄰近粒子；fi是質(zhì)量力,本文計(jì)算均取fi=0；ρi,mi,Vi,ui,ri分別代表粒子i的密度、質(zhì)量、體積、速度和位移；光滑核函數(shù)Wij采用Wendland的C2核函數(shù)[16],光滑長(zhǎng)度h=2Δx,其中Δx是初始粒子間距；耗散項(xiàng)Di和πi的具體參數(shù)可參考文獻(xiàn)[15]。Fij為壓力項(xiàng),在本文中采用張力不穩(wěn)定控制技術(shù)進(jìn)行處理。

為確保粒子的均勻分布,在δ+-SPH模型中,動(dòng)量方程計(jì)算得出的粒子位移ri需要施加人工位移修正量δri,δri具體定義為

(2)

2 張力不穩(wěn)定性控制技術(shù)

在傳統(tǒng)SPH模型中,由于負(fù)壓會(huì)導(dǎo)致張力不穩(wěn)定現(xiàn)象[18],流場(chǎng)信息將無法很好呈現(xiàn)?；赟un等[15]的研究結(jié)果,為更好實(shí)現(xiàn)漩渦附著、分離、脫落細(xì)節(jié)的捕捉,展開對(duì)Ishii翼型流場(chǎng)規(guī)律分析,本文采用TIC技術(shù)改進(jìn)壓力梯度項(xiàng)。對(duì)壓力項(xiàng)Fij進(jìn)行改進(jìn)[15]，如(3)式所示

(3)

DF表示液面區(qū)域自由液面粒子及相鄰粒子的集合。

3 拉格朗日擬序結(jié)構(gòu)可視化

在歐拉算法的模擬結(jié)果中,流場(chǎng)特征量(如渦量)通?；跉W拉框架下的瞬態(tài)速度場(chǎng)計(jì)算得出,難以給出流體介質(zhì)在某個(gè)時(shí)間段內(nèi)的輸運(yùn)特征,而本文采用拉格朗日擬序結(jié)構(gòu)對(duì)流場(chǎng)中漩渦結(jié)構(gòu)進(jìn)行可視化,可更加清晰地獲取流場(chǎng)細(xì)節(jié)。拉格朗日擬序結(jié)構(gòu)通過計(jì)算有限時(shí)間李亞普諾夫指數(shù)(finite-time Lyapunov exponent,FTLE)實(shí)現(xiàn)。在FTLE云圖中呈現(xiàn)的脊線為拉格朗日擬序結(jié)構(gòu)[19]。

在SPH計(jì)算程序中采用逆向FTLE展開對(duì)FTLE計(jì)算,以便能清晰反映漩渦結(jié)構(gòu)邊界特征[20]。逆向FTLE計(jì)算方程為

(4)

4 多級(jí)粒子分辨率技術(shù)

翼型的首尾部往往較為狹窄,翼型繞流中表面邊界層內(nèi)湍流流動(dòng)的捕捉對(duì)于網(wǎng)格法和無網(wǎng)格法均具有一定挑戰(zhàn)性。因此,本文將利用多級(jí)粒子分辨率技術(shù)[21],將翼型表面離散成分辨率高的固定虛粒子,實(shí)現(xiàn)精確捕捉湍流流動(dòng)細(xì)節(jié),提高邊界黏性力的計(jì)算精度。多級(jí)粒子分辨率技術(shù)原理如圖1所示。

圖1 多級(jí)粒子分辨率技術(shù)原理示意圖[21]

主要原理闡述如下:不同分辨率的粒子定義為不同的粒子集合。首先定義翼型附近區(qū)域?yàn)榧?xì)化區(qū),在細(xì)化區(qū)邊緣定義過渡區(qū),如圖1b)所示,充當(dāng)其邊界條件,獲取進(jìn)入細(xì)化區(qū)域母粒子的信息。過渡區(qū)域的寬度設(shè)置要大于核函數(shù)的半徑,防止其邊緣處核函數(shù)截?cái)喽斐捎?jì)算誤差。在計(jì)算中,參與流場(chǎng)SPH計(jì)算的粒子稱為SPH粒子;活性關(guān)閉,不參與SPH計(jì)算的粒子為非活性粒子。在細(xì)化區(qū),母粒子分裂,如圖1a)所示,在正方形4個(gè)頂點(diǎn)處形成更小粒子即子粒子,子粒子活性激活,參與流場(chǎng)的SPH計(jì)算,即為SPH粒子。而母粒子在撕裂之后,作為非活性粒子不再參與SPH計(jì)算,但仍然隨著流場(chǎng)運(yùn)動(dòng),通過從周圍細(xì)化的子粒子插值獲得流場(chǎng)信息如速度、壓力等。相反的,在過渡區(qū)域,母粒子為SPH粒子,子粒子屬性從周圍SPH粒子即母粒子插值得出。采用Shepard插值方法,非活性粒子的物理量φ從周圍SPH粒子完成插值,插值公式為

(5)

經(jīng)過此過程,2層分辨率粒子完成了相互之間的信息交換。多級(jí)粒子分辨率技術(shù)對(duì)局部流場(chǎng)粒子加密,在不同的分辨率區(qū)域之間,過渡區(qū)域及細(xì)化區(qū)域并不隨著粒子的移動(dòng)而移動(dòng)。在計(jì)算中,當(dāng)母粒子撕裂而來的子粒子流出細(xì)化區(qū)及過渡區(qū)域時(shí),不發(fā)生粒子合并而將直接被刪除,而母粒子流出細(xì)化區(qū)域時(shí),將正常參與SPH流場(chǎng)計(jì)算?；谝陨狭Ｗ蛹?xì)化思想,在多層分辨率粒子間,保持粒子尺度一樣,每層粒子的過渡區(qū)域流場(chǎng)信息從下一層插值取得,從而實(shí)現(xiàn)不同層次間粒子信息交換,實(shí)現(xiàn)多級(jí)分辨率的SPH數(shù)值計(jì)算,提高流場(chǎng)計(jì)算精度。在下文將首先通過圓柱繞流基準(zhǔn)算例驗(yàn)證在δ+-SPH中添加多級(jí)粒子分辨率技術(shù)的有效性,從而應(yīng)用于翼型繞流問題的研究中。

5 基準(zhǔn)算例驗(yàn)證

圓柱繞流問題是流體力學(xué)中的經(jīng)典問題,是檢驗(yàn)數(shù)值算法的經(jīng)典算例。選取雷諾數(shù)Re=3 000,人工聲速c0=15 m/s時(shí)的圓柱繞流模擬算例,與試驗(yàn)結(jié)果對(duì)比以檢驗(yàn)δ+-SPH模型的數(shù)值精度。

在該算例中,施加多級(jí)粒子分辨率技術(shù),采用5層分辨率不同的粒子。圓柱表面的粒子間距為D/Δx=400,最遠(yuǎn)處粒子間距為D/Δx=25。在計(jì)算中,粒子總數(shù)為1.45×106。假設(shè)在計(jì)算域中粒子分辨率均取為D/Δx=400,粒子總數(shù)將達(dá)到7.68×107?？梢?多級(jí)粒子分辨率技術(shù)對(duì)重點(diǎn)區(qū)域流場(chǎng)進(jìn)行加密,在保證計(jì)算精度條件下,有效減少了粒子總數(shù),從而減少計(jì)算時(shí)間,提高計(jì)算效率,對(duì)實(shí)現(xiàn)SPH方法在工程上應(yīng)用具有重要意義。

圖2 雷諾數(shù)Re=9 500圓柱繞流流場(chǎng)特征

圖2中選取圓柱上表面的δ+-SPH計(jì)算圖展示圓柱繞流的漩渦發(fā)展過程。通過漩渦流線圖與試驗(yàn)結(jié)果[22]對(duì)比,可以發(fā)現(xiàn)試驗(yàn)結(jié)果與δ+-SPH結(jié)果吻合良好。δ+-SPH模型與多級(jí)粒子分辨率技術(shù)相結(jié)合,漩渦發(fā)展細(xì)節(jié)能較好被模擬,δ+-SPH結(jié)合多級(jí)粒子分辨率技術(shù)可以較為準(zhǔn)確模擬繞流現(xiàn)象中分離、附著現(xiàn)象。SPH中拉格朗日擬序結(jié)構(gòu)能清晰展示漩渦演化過程,呈現(xiàn)清晰的內(nèi)部細(xì)節(jié),與試驗(yàn)結(jié)果基本保持一致,說明拉格朗日擬序結(jié)構(gòu)獲取流場(chǎng)特征的高精度?；谝陨嫌懻?可見SPH方法在漩渦結(jié)構(gòu)捕捉等方面具有獨(dú)特優(yōu)勢(shì)。在下節(jié)中,本文將采用多級(jí)粒子分辨率技術(shù)對(duì)計(jì)算域粒子進(jìn)行處理,結(jié)合δ+-SPH計(jì)算結(jié)果中的拉格朗日擬序結(jié)構(gòu)展開對(duì)翼型繞流現(xiàn)象研究。

6 Ishii翼型繞流

6.1 Ishii翼型繞流的初始工況和邊界條件設(shè)置

Ishii翼型具有翼厚比低、前緣半徑小等特征,翼型形狀[23]如圖3所示。本文模擬將Ishii翼型固定在矩形流場(chǎng)中,流場(chǎng)左側(cè)為流入邊界,來流速度設(shè)為U0=1 m/s,右側(cè)為流出邊界,上下邊界為自由滑移壁面邊界。翼型的上表面施加無滑移邊界條件。

圖3 攻角α=6°時(shí)的Ishii翼型形狀

對(duì)Ishii翼型表面的無量綱量升阻力系數(shù)測(cè)量計(jì)算,獲取翼型的氣動(dòng)力特性。對(duì)阻力系數(shù)CD,升力系數(shù)CL定義為：

(6)

式中，fD和fL分別定義為繞流Ishii翼型的阻力和升力,在SPH程序計(jì)算中,這兩者的計(jì)算基于流體粒子與結(jié)構(gòu)粒子之間相互作用力的平衡算得。在本文翼型模擬中,弦長(zhǎng)L=1 m。

6.2 多級(jí)粒子分辨率技術(shù)對(duì)計(jì)算效率的提升

為提高計(jì)算效率,采用6層不同的粒子分辨率,在翼型表面的最高粒子分辨率為L(zhǎng)/Δx=800,最低粒子分辨率即距翼型最遠(yuǎn)處取L/Δx=25。在本計(jì)算中粒子總數(shù)為6.40×105,如果采用L/Δx=800對(duì)全流場(chǎng)進(jìn)行計(jì)算,粒子總數(shù)將達(dá)到3.07×107,可見采用多級(jí)粒子分辨率技術(shù),粒子總數(shù)減少將近43倍,除去不同尺度粒子之間的插值耗時(shí),多級(jí)粒子分辨率技術(shù)可將計(jì)算效率提高約20倍。

6.3 雷諾數(shù)Re=3 000下繞流模擬結(jié)果和驗(yàn)證

首先,在雷諾數(shù)Re=3 000、攻角α=0°條件下,模擬了Ishii翼型繞流特性。阻力系數(shù)的計(jì)算結(jié)果如圖4所示。FVM結(jié)果與SPH結(jié)果基本吻合,驗(yàn)證了SPH方法的計(jì)算精度。

圖4 雷諾數(shù)Re=3 000, 攻角α=0°Ishii翼型阻力曲線

6.4 雷諾數(shù)Re=30 000下繞流模擬結(jié)果和分析

6.4.1 升阻力系數(shù)

隨后,在雷諾數(shù)Re=30 000、攻角α=0°,3°,6°條件下,模擬了Ishii翼型繞流的特性。α=6°時(shí),翼型阻力及升力系數(shù)的SPH計(jì)算結(jié)果與FVM結(jié)果對(duì)比如圖5所示。在粒子初始階段到歷經(jīng)劇烈抖動(dòng)之后,最終CL在0.60～0.85之間動(dòng)態(tài)波動(dòng)。同樣地,曲線可發(fā)現(xiàn)阻力CD系數(shù)在0.044上下波動(dòng)。

圖5 雷諾數(shù)Re=30 000下的升阻力系數(shù)時(shí)歷曲線

表1給出了本文計(jì)算結(jié)果及文獻(xiàn)中試驗(yàn)結(jié)果[23]與其他數(shù)值結(jié)果[24](2-D Lam、2-D RANS)的對(duì)比?？梢钥闯?多種算法的升阻力系數(shù)值基本吻合,SPH結(jié)果與2-D Lam計(jì)算結(jié)果吻合最佳。通過比較圖6中攻角α=0°,3°,6°的升力系數(shù)可以看出,攻角變大、升力系數(shù)增加。此外,攻角增大引起黏性漩渦的非周期性脫落,使得升力系數(shù)隨時(shí)間產(chǎn)生振蕩。

表1 升阻力結(jié)果驗(yàn)證

圖6 不同攻角α下升力系數(shù)時(shí)歷曲線

6.4.2 SPH計(jì)算黏性繞流問題的速度場(chǎng)與壓力場(chǎng)

圖7展示了雷諾數(shù)Re=30 000、攻角α=6°條件下,Ishii翼型繞流的速度場(chǎng)和粒子分布情況。

圖7 Ishii翼型繞流的速度場(chǎng)和粒子分布

對(duì)比t=0.7 s時(shí)刻,施加與不施加PST時(shí)的粒子分布,可見：不施加PST時(shí),翼型周圍粒子分布十分不規(guī)則,且速度場(chǎng)呈現(xiàn)出不穩(wěn)定;而施加PST,翼型表面邊界層內(nèi)流體粒子分布十分均勻,保證了速度場(chǎng)計(jì)算的精度。

在SPH模擬中,實(shí)現(xiàn)穩(wěn)定的壓力場(chǎng)計(jì)算至關(guān)重要。然而在傳統(tǒng)SPH算法中,常常存在壓力振蕩及張力不穩(wěn)定性現(xiàn)象。在δ+-SPH模型中,通過在連續(xù)方程中添加密度耗散項(xiàng)[14],可得到光滑的壓力場(chǎng)。針對(duì)張力不穩(wěn)定現(xiàn)象,本文采用張力不穩(wěn)定性控制技術(shù)進(jìn)行解決。為驗(yàn)證其有效性,圖8中給出Ishii翼型繞流模擬中施加TIC與不施加TIC的壓力云圖?？梢?不施加TIC時(shí),翼型周圍流場(chǎng)中出現(xiàn)數(shù)值空洞,導(dǎo)致計(jì)算結(jié)果失真;施加TIC時(shí),Ishii翼型周圍壓力分布均勻,且沒有數(shù)值噪聲。

圖8 壓力場(chǎng)云圖(t=4.076 s,Re=30 000)

以上對(duì)比說明,δ+-SPH模型結(jié)合TIC技術(shù)可以有效消除壓力噪聲,進(jìn)而較好實(shí)現(xiàn)翼型的壓力數(shù)值變化預(yù)報(bào)。此外,觀察壓力云圖可以發(fā)現(xiàn),翼型首部附近區(qū)域存在巨大負(fù)壓區(qū),這是引起張力不穩(wěn)定現(xiàn)象的主要原因。因此,采取張力不穩(wěn)定性控制技術(shù)是準(zhǔn)確獲取流場(chǎng)特征的重要前提。

6.4.3α=6°時(shí)Ishii翼型繞流中漩渦結(jié)構(gòu)分析

圖9 翼型周圍流場(chǎng)中的漩渦結(jié)構(gòu)(t=12.326 s,Re=30 000)

圖9給出了翼型周圍流場(chǎng)中的拉格朗日擬序結(jié)構(gòu)和渦量場(chǎng)?？梢?翼型表面附近高分辨率粒子的布置使得邊界層內(nèi)湍流細(xì)節(jié)能能夠清晰地捕捉,反映出多級(jí)粒子分辨率技術(shù)在高雷諾數(shù)下黏性繞流問題中的重要作用。相比于渦量場(chǎng)(網(wǎng)格法較常使用該方法捕捉漩渦結(jié)構(gòu)),在SPH框架中捕捉的拉格朗日擬序結(jié)構(gòu)更能清晰重現(xiàn)渦結(jié)構(gòu)的輪廓和內(nèi)部細(xì)節(jié)。此外,相較于渦量場(chǎng)云圖,FTLE云圖受閾值的影響不敏感,更加便于進(jìn)行流場(chǎng)細(xì)節(jié)分析。這也是SPH粒子法處理翼型繞流問題的優(yōu)勢(shì)之一。

從圖9可見,在翼型形狀及壓力梯度影響下,Ishii翼型漩渦的生成、附著、分離等主要發(fā)生在翼型上緣?；诙嗉?jí)粒子分辨率技術(shù),可以發(fā)現(xiàn),尺度很小的漩渦也能清晰捕捉,再次說明了多級(jí)粒子分辨率技術(shù)的有效性。

7 拓展應(yīng)用:拍動(dòng)翼型繞流問題

為進(jìn)一步展示SPH在處理運(yùn)動(dòng)翼型繞流問題中的優(yōu)勢(shì),本節(jié)采用δ+-SPH模型對(duì)拍動(dòng)水翼問題進(jìn)行模擬和驗(yàn)證。在網(wǎng)格法中,對(duì)運(yùn)動(dòng)翼型繞流問題的模擬容易引起網(wǎng)格變形而降低計(jì)算精度。即使采用重疊網(wǎng)格法,也常常因涉及不同網(wǎng)格間的插值而使算法變得復(fù)雜。相比而言,得益于SPH的拉格朗日粒子特性,對(duì)于大幅運(yùn)動(dòng)邊界問題的模擬較為方便。

本節(jié)選擇頭部半圓尾部尖角的拍動(dòng)翼型(具體細(xì)節(jié)參考文獻(xiàn)[25])進(jìn)行SPH計(jì)算與驗(yàn)證;采用無量綱參數(shù)包括雷諾數(shù)Re、斯特勞哈爾數(shù)Sr、幅值無量綱系數(shù)AD表征流體黏性、翼型拍動(dòng)頻率、翼型拍動(dòng)振幅對(duì)翼型繞流模擬的影響。拍動(dòng)幅值無量綱系數(shù)AD=2A/D,其中,A為翼型拍動(dòng)最大幅值,D取為翼型頭部半圓直徑。

圖10給出了雷諾數(shù)Re=440,斯特勞哈爾數(shù)Sr=0.035,無量綱幅值系數(shù)AD=1.47條件下,SPH模擬的漩渦結(jié)構(gòu)與試驗(yàn)結(jié)果[25]的對(duì)比。在t1時(shí)刻,SPH結(jié)果與試驗(yàn)結(jié)果[25]吻合較好,體現(xiàn)了δ+-SPH模型在運(yùn)動(dòng)翼型繞流模擬上的精度和獨(dú)特優(yōu)勢(shì)。觀察SPH不同時(shí)刻的結(jié)果,可見拉格朗日擬序結(jié)構(gòu)能呈現(xiàn)拍動(dòng)翼周期拍動(dòng)時(shí)漩渦演化過程,這為未來翼型推進(jìn)性能與漩渦關(guān)系的研究提供了新思路。

圖10 拍動(dòng)水翼繞流尾跡的試驗(yàn)結(jié)果[25]與SPH結(jié)果

8 結(jié) 論

本文采用δ+-SPH方法對(duì)不同雷諾數(shù)和攻角下的翼型繞流問題進(jìn)行了數(shù)值模擬研究,分析了升阻力特性和流場(chǎng)漩渦特征,并將SPH方法從固定翼繞流拓展應(yīng)用到了拍動(dòng)翼繞流的模擬中。研究結(jié)果表明:

1) SPH方法在模擬真實(shí)黏性邊界層和漩渦運(yùn)動(dòng)的流體動(dòng)力學(xué)問題中具有競(jìng)爭(zhēng)力。不同于網(wǎng)格法,SPH方法易于實(shí)現(xiàn)拉格朗日擬序結(jié)構(gòu)的捕捉,動(dòng)態(tài)呈現(xiàn)漩渦結(jié)構(gòu)演化過程,因此在未來探究漩渦與流體動(dòng)力特性關(guān)系上具有重要優(yōu)勢(shì)。

2) 采用TIC技術(shù)及PST技術(shù)能夠有效提高翼型繞流的速度場(chǎng)和壓力場(chǎng)計(jì)算精度,并實(shí)現(xiàn)對(duì)升阻力系數(shù)較高精度的預(yù)報(bào);采用多級(jí)粒子分辨率技術(shù)能夠有效減少粒子總數(shù),縮短計(jì)算時(shí)間,從而便于在工程中應(yīng)用。

3) SPH方法不需要生成網(wǎng)格,能自動(dòng)追蹤運(yùn)動(dòng)和變形的流固界面,而且該方法的拉格朗日和無網(wǎng)格粒子特性使其在模擬運(yùn)動(dòng)翼型繞流問題上具有先天優(yōu)勢(shì)和巨大潛力。在未來研究中,利用該優(yōu)勢(shì),SPH方法在翼型繞流的渦振或顫振問題研究中有望發(fā)揮重要作用。