基于cfd的雙槳式吊艙驅(qū)動器水動力性能數(shù)值預(yù)報

基于cfd的雙槳式吊艙驅(qū)動器水動力性能數(shù)值預(yù)報

1船舶首、船側(cè)、船尾布置水動力作為一種新型的能源集成鼓勵系統(tǒng)，該裝置已廣泛應(yīng)用于許多真實船只。吊艙推進器集成電機和螺旋槳于一體實現(xiàn)模塊化設(shè)計,其在船舶上的布置也不再像傳統(tǒng)螺旋槳那樣受到軸系、主機的嚴(yán)格限制,可以布置在船首、船側(cè)或者船尾等,目前其作為主推進器一般是吊掛在船尾下方,進流比傳統(tǒng)螺旋槳要均勻,這在一定程度上有望改善其水動力、空泡、振動和噪聲等性能。在吊艙推進器研究領(lǐng)域,如何準(zhǔn)確預(yù)報其水動力性能一直是研究的重點之一,因為這關(guān)系到吊艙推進器設(shè)計及性能的提高,也為工程實踐提供重要的依據(jù)。雙槳式吊艙推進器由于在吊艙包前后各布置一個螺旋槳,所以在雙槳及吊艙包的相互影響下,其水動力性能也更為復(fù)雜。本文采用粘流CFD方法對雙槳式吊艙推進器進行了在均流和非均流條件下的水動力性能預(yù)報,并與試驗結(jié)果進行了比較分析,驗證了該方法的可靠性。2cfd理論2.1陣風(fēng)下雷達(dá)平均法模型在工程實際中,流動一般為湍流。而對于湍流的描述,從工程應(yīng)用角度來看,重要的是湍流所引起的平均流場的變化,是整體的效果。由于在湍流流動中,物理量可以分解成一個平均量與一個脈動量的和,因此在模擬湍流時就只計算平均量,脈動量的影響通過某種模型在時均化的方程中體現(xiàn)出來,該方法稱作雷諾平均法,也稱作RANS方法。相應(yīng)的雷諾平均N-S方程為:方程(2)中的-ρui′uj′項,稱之為雷諾應(yīng)力項。雷諾應(yīng)力項屬于新的未知量,因此,要使方程封閉,必須對該應(yīng)力項作某種假設(shè),即建立應(yīng)力的表達(dá)式(或者引進新的湍流模型方程),通過表達(dá)式或湍流模型,把應(yīng)力項中的脈動值與時均值聯(lián)系起來。2.2湍流粘度的定義本文選取理論上發(fā)展較為完善的SSTk-ω模型(shear-stress-transport)來封閉RANS方程。SSTk-ω兩方程湍流模型如下:式中:表示由平均速度梯度引起的湍動能的產(chǎn)生項。Gω表示ω的產(chǎn)生項。Γk和Γω表示k和ω的有效擴散。Yk和Yω分別表示k和ω的耗散。Dω表示交叉擴散項。所有以上這些參數(shù)都定義如下:在這里,σk和σω分別是關(guān)于k和ω的湍流普朗特數(shù)。湍流粘度定義如下:式中S表征了渦量的大小。其中Gk同標(biāo)準(zhǔn)k-ω模型中所定義的湍動能產(chǎn)生項相同。ω的產(chǎn)生項定義為:以上各式中,F1和F2均為混合函數(shù);各經(jīng)驗系數(shù)定義及取值可參考文獻(xiàn)。2.3計算值的方法2.3.1子域控制方程多參考系模型(MultipleReferenceFrame(MRF)model)是一種常用的存在相對運動時的定常計算模型,特別是運動區(qū)域與靜止區(qū)域間的相互作用比較微弱時。MRF模型的另外一個用途就是用來為非定常計算提供初始流場。在MRF方法中,計算區(qū)域劃分為多個子域,每個子域的控制方程是針對子域參考系而寫的。在兩子域間的邊界,子域控制方程的擴散項以及其他項需要鄰近子域的速度值。使用相對速度公式,對每個子域的速度相對于子域的運動進行計算。從而使得速度和速度梯度從移動參考系轉(zhuǎn)換到絕對慣性參考系。2.3.2滑移網(wǎng)格技術(shù)保證了數(shù)值模擬的準(zhǔn)確性滑移網(wǎng)格模型(SlidingMeshmodel)假定流場是非定常的,因此數(shù)值模擬更精確地反應(yīng)流場中旋轉(zhuǎn)物體間的相互作用,從而保證了數(shù)值模擬的精確度?；凭W(wǎng)格技術(shù)的基本原理是將計算模型劃分成多個部分,每個部分網(wǎng)格有自己獨立的網(wǎng)格形式和網(wǎng)格邊界面。兩部分之間邊界面的組合稱為交界面,交界面兩側(cè)網(wǎng)格相互滑動,兩側(cè)的網(wǎng)格結(jié)點不要求相互重合,通過計算兩側(cè)的通量,使其相等。為了計算交界面的通量,首先在每一個新的時間步確定出交界面兩邊交界區(qū)的重合面,通過網(wǎng)格重合面的通量由交界面兩邊交界區(qū)的重合面計算。3計算模型的建立及網(wǎng)格劃分本文所研究的雙槳式吊艙推進器在上海交通大學(xué)空泡水筒中進行了均勻來流狀態(tài)下的水動力性能試驗,試驗按ITTC推薦的程序進行,結(jié)果已進行了筒壁修正。表1為吊艙推進器前、后螺旋槳的主要幾何參數(shù)。圖1為計算對象雙槳式吊艙推進器模型在空泡水筒中試驗的照片。圖2為計算模型的表面網(wǎng)格及計算域縱剖面局部網(wǎng)格。圖3為計算域網(wǎng)格的劃分,將計算域一共劃分為三個區(qū)域:包圍前槳葉的圓柱體旋轉(zhuǎn)域1、包圍后槳葉的圓柱體旋轉(zhuǎn)域2以及一個靜止的計算域3。在兩個旋轉(zhuǎn)域,域1、域2內(nèi)劃分網(wǎng)格時進行局部加密,以保證近壁面網(wǎng)格Y+值的取值范圍(11.63<Y+<500),提高計算結(jié)果的準(zhǔn)確度;域3為吊艙推進器所在的控制域,采用結(jié)構(gòu)化網(wǎng)格,從而達(dá)到合理的減少計算網(wǎng)格的目的。4在均勻流條件下對水動力性能進行計算4.1多參考系和滑移網(wǎng)格相結(jié)合的數(shù)值預(yù)報入口設(shè)置為來流速度入口,出口定義為壓力出口,物面為無滑移表面。特別地,在模擬雙槳吊艙推進器各部件相互干擾情況時,考慮到后槳在前槳的尾流中運轉(zhuǎn),受前槳的影響很大,而同時后槳旋轉(zhuǎn)產(chǎn)生的誘導(dǎo)速度場也會對前槳的性能有所影響。因此這里沒有僅僅使用數(shù)值預(yù)報單槳定常水動力性能中經(jīng)常使用的多參考系方法,而是使用了多參考系和滑移網(wǎng)格相結(jié)合的方法,即前槳采用多參考系方法,后槳采用滑移網(wǎng)格方法。基于壓力變量的SIMPLE耦合求解器,采用二階迎風(fēng)差分離散格式,結(jié)合低雷諾修正SSTk-ω兩方程湍流模型進行數(shù)值模擬。一般迭代計算2000步就可以收斂。4.2試驗結(jié)果比較采用上述數(shù)值方法,對雙槳式吊艙推進器水動力性能進行了預(yù)報,螺旋槳的推力系數(shù)KT、轉(zhuǎn)矩系數(shù)10KQ以及敞水效率η計算結(jié)果(N=1372rpm)見表2,與試驗值的比較如圖4所示。從計算結(jié)果來看,對于雙槳式吊艙推進器的推力系數(shù)KT、轉(zhuǎn)矩系數(shù)10KQ的預(yù)報都較為準(zhǔn)確,推力的誤差不大于4.51%,扭矩誤差不大于3.35%。5在非均質(zhì)流條件下對水動力性能進行計算5.1吊艙流場分布本文為了使非定常特性更為明顯,便于驗證,選取某一散貨船的船后伴流場為給定計算的流場,其分布由圖5給出,而真實的吊艙流場品質(zhì)要比船后的好很多。在FLUENT軟件中,以Profile的形式導(dǎo)入伴流數(shù)據(jù)文件。數(shù)值計算取轉(zhuǎn)速N=1493rpm,此時Vmean=3.602。計算中前、后槳均采用滑移網(wǎng)格方法。5.2結(jié)論分析5.2.1旋轉(zhuǎn)周期內(nèi)總推力分析如圖6至圖9所示,雙槳式吊艙推進器前后槳推力系數(shù)、扭矩系數(shù),及總推力系數(shù)、總扭矩系數(shù)的時域和頻域結(jié)果可以看出,螺旋槳(三葉槳)非定常力的葉頻特征明顯,且以一倍葉頻為主。但對于推進器總推力系數(shù),呈現(xiàn)出六倍軸頻幅值最大的特點,且該幅值小于后槳的一倍葉頻(三倍軸頻)幅值,如圖7所示,說明雙槳式吊艙推進器在前后槳非定常力相互耦合影響下,總的非定常力可能減小,且主要特征頻率為前后槳葉數(shù)總和的一倍葉頻。如圖10至圖11所示,前(后)螺旋槳的側(cè)向力、垂向力的非定常力的一倍葉頻特征明顯,且前槳與后槳在時域上相鄰峰值的角度差值與前后槳之間安裝相位角度差相近。從預(yù)報精度看,一個旋轉(zhuǎn)周期內(nèi)總推力系數(shù)KT-pro的平均值為0.2622,螺旋槳總轉(zhuǎn)矩系數(shù)10KQ-pro的平均值為0.4692,與對應(yīng)工況J=0.6885時,推進器敞水推力系數(shù)的試驗值為0.2491,轉(zhuǎn)矩系數(shù)試驗值為0.4550進行比較,總推力系數(shù)偏差為5.26%,扭矩系數(shù)偏差為3.12%。因此,非均勻伴流情況下雙槳式吊艙推進器非定常水動力預(yù)報精度可滿足要求。5.2.2支架軸向力及側(cè)向力支架在推進器運轉(zhuǎn)過程中,也受到螺旋槳非定常力的作用,由圖12至圖15可知,雙槳式吊艙推進器支架受到的軸向力、側(cè)向力和垂向力的非定常力的成分均以三葉螺旋槳的一倍葉頻為主。支架軸向力系數(shù)KTx-str的數(shù)值為負(fù),說明在吊艙推進器工作時支架提供阻力;支架側(cè)向力系數(shù)KTz-str與螺旋槳側(cè)向力系數(shù)KTz-pro1.0和KTz-pro1.2處于同一數(shù)值量級上,表明支架側(cè)向力與螺旋槳側(cè)向力相關(guān)性較強,且支架側(cè)向力對推進器操縱性也有較大影響,需在設(shè)計中關(guān)注;支架垂向力系數(shù)KTy-str的數(shù)值為負(fù),且脈動幅值最小。5.2.3垂向力及支架軸向力艙體在推進器運轉(zhuǎn)過程中,也受到螺旋槳非定常力的作用,由圖16至圖19可知,雙槳式吊艙推進器艙體受到的軸向力、側(cè)向力和垂向力的非定常力的成分均以三葉螺旋槳的一倍葉頻為主。艙體軸向力系數(shù)KTx-pod數(shù)值為負(fù),說明在吊艙推進器工作時艙體提供阻力,量級和支架軸向力系數(shù)相當(dāng),兩者之和約占螺旋槳總軸向力系數(shù)的6%,不能忽略;艙體側(cè)向力系數(shù)KTz-pod與垂向力系數(shù)KTy-pod處于同一數(shù)值量級上,且與支架側(cè)向力系數(shù)KTz-str相當(dāng),這些在吊艙推進器結(jié)構(gòu)設(shè)計和操縱性研究中也需予以重視。6數(shù)值預(yù)報驗證研究并建立了基于粘流理論,在均流條件下的雙槳式吊艙推進器水動力性能的數(shù)值預(yù)報方法。通過將數(shù)值計算結(jié)果與實驗結(jié)果進行了比較、分析,在設(shè)計工況J=0.9附近,推力系數(shù)的誤差為2.68%,扭矩系數(shù)的誤差為2.14%,驗證了該數(shù)值求解方法對于雙槳式吊艙推進器在均流條件下的水動力性能的數(shù)值預(yù)報