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

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

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