• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    一種改進(jìn)誘導(dǎo)因子計(jì)算方法的葉素動(dòng)量理論體積力模型

    2024-01-19 06:57:18王明哲王建華萬德成
    船舶力學(xué) 2024年1期
    關(guān)鍵詞:葉素動(dòng)量螺旋槳

    王明哲,王建華,萬德成

    (上海交通大學(xué)a.船海計(jì)算水動(dòng)力學(xué)研究中心(CMHL);b.船舶海洋與建筑工程學(xué)院,上海 200240)

    0 引 言

    在進(jìn)行螺旋槳與其他結(jié)構(gòu)物相互作用的數(shù)值模擬中,螺旋槳模型因?yàn)閺?fù)雜的旋轉(zhuǎn)運(yùn)動(dòng)以及復(fù)雜的幾何結(jié)構(gòu),往往會(huì)需要較長的計(jì)算時(shí)間和較大的存儲(chǔ)誤差成本。體積力模型通過將螺旋槳與流體的相互作用力?;癁榉植荚诹鲌?chǎng)中的體積力,作為源項(xiàng)直接參與動(dòng)量方程的求解,是代替真實(shí)螺旋槳模型的減小計(jì)算成本的有效手段,在船槳相互作用、槳舵相互作用以及船舶操縱性能模擬等方面得到了大量的應(yīng)用。

    作為一種簡(jiǎn)化模型,體積力模型的復(fù)雜性與準(zhǔn)確性是影響其應(yīng)用性的關(guān)鍵因素。葉素動(dòng)量理論模型作為一種復(fù)雜性與準(zhǔn)確性適中的體積力模型,具有較大的應(yīng)用潛力[1-3]。

    使用葉素動(dòng)量理論計(jì)算得到的螺旋槳性能與真實(shí)螺旋槳模型性能接近,但難以在不同工況下保證相同的準(zhǔn)確性。Benini[4]研究了使用葉素動(dòng)量理論模擬船用螺旋槳的適用性,其中葉元參數(shù)由基于勢(shì)流理論的二維葉元性能計(jì)算軟件XFOIL 計(jì)算,結(jié)果表明使用葉素動(dòng)量理論計(jì)算螺旋槳性能的最小誤差為2%,但誤差對(duì)進(jìn)速系數(shù)的變化較為敏感;Phillips[1]使用基于葉素動(dòng)量理論的體積力模型進(jìn)行了KVLCC 的PMM 試驗(yàn)數(shù)值模擬,忽略了螺旋槳載荷分布的周向不均勻性,最后發(fā)現(xiàn)在小舵角工況下船體側(cè)向力與艏搖力矩預(yù)測(cè)誤差在(2~3)%以內(nèi),而當(dāng)舵角增加時(shí)預(yù)報(bào)誤差增大。種種研究表明,可以通過修正手段提高葉素動(dòng)量理論體積力模型的準(zhǔn)確性。

    葉元參數(shù)準(zhǔn)確性與葉元水動(dòng)力計(jì)算方法是影響葉素動(dòng)量理論體積力模型的關(guān)鍵因素。對(duì)于葉元參數(shù),已存在部分經(jīng)驗(yàn)公式方法用于修正非設(shè)計(jì)工況下的螺旋槳葉元參數(shù),但無法保證此修正方法對(duì)于不同槳型及工況的適用性[5-6]。葉元水動(dòng)力計(jì)算方法包括直接使用當(dāng)?shù)厮俣纫约暗?jì)算誘導(dǎo)因子。馮大奎等[7-10]基于KP505 真實(shí)螺旋槳模型的敞水試驗(yàn)?zāi)M結(jié)果建立了葉元參數(shù)與當(dāng)?shù)厮俣鹊年P(guān)系,分析了不同葉數(shù)下葉元的性能變化,進(jìn)行了KCS 模型尺度與實(shí)尺度在設(shè)計(jì)航速下的自航試驗(yàn)?zāi)M,模擬結(jié)果表明體積力模型與真實(shí)螺旋槳得到的自航因子結(jié)果吻合;Tokgoz[11]通過葉元參數(shù)的一般表達(dá)式與直接使用當(dāng)?shù)厮俣确椒ㄟM(jìn)行了AU螺旋槳敞水與斜流試驗(yàn)?zāi)M,結(jié)果表明敞水條件下得到的螺旋槳扭矩誤差較大,且進(jìn)速系數(shù)越低準(zhǔn)確性越低。

    傳統(tǒng)的誘導(dǎo)因子迭代方法基于理想流體的動(dòng)量理論推導(dǎo)得出,未考慮粘性求解器中的粘流環(huán)境[2]。Villa[12]擬合了單個(gè)進(jìn)速工況下槳葉周向平均載荷的徑向分布函數(shù),其中槳葉載荷由敞水實(shí)驗(yàn)得到,并將此載荷分布作為體積力分布與RANS 求解器結(jié)合,結(jié)果表明該體積力模型與真實(shí)螺旋槳模型作用下的尾流速度分布吻合程度高,證明了按照真實(shí)螺旋槳載荷分布的體積力模型可以得到預(yù)期的誘導(dǎo)速度場(chǎng)。但此種方法需要已知工況,無法處理來流速度變化的非定常及非均勻流場(chǎng)。傅慧萍[13]使用HO模型模擬了KCS船自航,發(fā)現(xiàn)進(jìn)速與槳盤面當(dāng)?shù)厮俣冉咏€性關(guān)系。此研究結(jié)果表明在得到當(dāng)?shù)厮俣鹊那闆r下,可以通過此對(duì)應(yīng)關(guān)系根據(jù)槳盤面當(dāng)?shù)厮俣确赐七M(jìn)速,進(jìn)而得到粘流條件下槳盤面各處誘導(dǎo)因子大小。

    本文通過不同進(jìn)速下螺旋槳敞水試驗(yàn)的CFD 數(shù)值模擬,得到在粘性作用與葉間干擾作用下的螺旋槳各葉元參數(shù),并依據(jù)此葉元參數(shù)得到符合真實(shí)螺旋槳載荷徑向分布的投影盤,進(jìn)而得到此投影盤作用下各進(jìn)速的敞水誘導(dǎo)速度場(chǎng)以及槳盤面處各位置的當(dāng)?shù)厮俣?,最終通過當(dāng)?shù)厮俣扰c進(jìn)速計(jì)算誘導(dǎo)因子分布。此誘導(dǎo)因子分布將植入葉素動(dòng)量理論體積力模型,避免了誘導(dǎo)因子的迭代求解,并通過KP505螺旋槳真實(shí)模型與體積力模型的敞水試驗(yàn)驗(yàn)證本文的體積力改進(jìn)方法。

    1 數(shù)值方法

    1.1 控制方程

    本文使用改進(jìn)的葉素動(dòng)量理論體積力模型與真實(shí)螺旋槳進(jìn)行螺旋槳敞水試驗(yàn)的數(shù)值模擬,使用的基礎(chǔ)數(shù)值求解器為開源平臺(tái)OpenFOAM 中的pimpleFoam 與pimpleDyMFoam,其中pimpleDyMFoam用于真實(shí)螺旋槳模型數(shù)值模擬,pimpleFoam 則用于實(shí)現(xiàn)體積力模型代碼在粘流求解器中的植入。對(duì)于粘性不可壓縮流體,控制方程為

    式中,p為壓力,ρ為流體密度,μ為粘性系數(shù),ui為速度分量為雷諾應(yīng)力。本文使用湍流模型封閉含有雷諾應(yīng)力項(xiàng)的方程。對(duì)于最適用于螺旋槳敞水試驗(yàn)數(shù)值模型的湍流模型仍無定論[14],文中所使用的為k-ωSST 湍流模型,這種湍流模型針對(duì)根據(jù)各區(qū)域到壁面距離的變化在對(duì)應(yīng)位置計(jì)算時(shí)激活不同的湍流變量控制方程形式,是一種應(yīng)用性較廣的方法[15]。(fε)i為體積力源項(xiàng),由單獨(dú)的葉素動(dòng)量理論體積力代碼計(jì)算得到[16],在使用體積力模型進(jìn)行數(shù)值模擬時(shí)激活此項(xiàng)。體積力模型與真實(shí)螺旋槳模型均使用非定常模擬,用以驗(yàn)證體積力模型在處理非定常問題時(shí)的計(jì)算能力。速度場(chǎng)與壓力場(chǎng)的解耦使用PIMPLE 方法,體積力源項(xiàng)的計(jì)算采用顯式方法,即首先利用上一時(shí)間步的速度場(chǎng)計(jì)算體積力分布,再使用得到的體積力分布計(jì)算當(dāng)前時(shí)間步的速度場(chǎng)。

    1.2 葉素動(dòng)量理論

    每次完成控制方程的求解后,將進(jìn)行一次葉素動(dòng)量理論程序的計(jì)算以求取體積力分布。文中各葉元體作用于流場(chǎng)中稱為致動(dòng)點(diǎn)的位置,致動(dòng)點(diǎn)的位置不隨流場(chǎng)的迭代而發(fā)生改變,即螺旋槳轉(zhuǎn)速僅為計(jì)算葉元體水動(dòng)力性能的參量。致動(dòng)點(diǎn)處的葉元體水動(dòng)力在每個(gè)時(shí)間步根據(jù)瞬時(shí)當(dāng)?shù)厮俣扔?jì)算,且各致動(dòng)點(diǎn)互不影響,因此對(duì)于船槳配合的自航問題中螺旋槳載荷徑向與周向分布不均勻的問題,此體積力模型也可以進(jìn)行有效處理,使用致動(dòng)點(diǎn)的軸向當(dāng)?shù)厮俣萔x與周向當(dāng)?shù)厮俣萔θ依照下式計(jì)算葉元體的水動(dòng)力,圖1展示了葉元體水動(dòng)力計(jì)算各相關(guān)物理量的矢量關(guān)系:

    圖1 葉素動(dòng)量理論Fig.1 Blade element momentum theory

    式中:a與b分別表示軸向誘導(dǎo)因子和切向誘導(dǎo)因子,表征螺旋槳引起的誘導(dǎo)速度與來流速度的比值;n為螺旋槳轉(zhuǎn)速;r為葉元體徑向位置,;Vx與Vθ分別表示當(dāng)?shù)剌S向速度和當(dāng)?shù)厍邢蛩俣龋绰菪龢獢_動(dòng)產(chǎn)生的誘導(dǎo)速度與來流速度的和;Va與Vt分別表示大地坐標(biāo)系下的軸向來流速度和切向來流速度,表征無窮遠(yuǎn)處流體相對(duì)于葉元體的速度;β為進(jìn)角;φ為葉元體螺距角;α為幾何攻角,;CL與CD分別表示葉元體的升力系數(shù)和阻力系數(shù);cr為對(duì)應(yīng)徑向位置處的葉元體弦長;dr為葉元體展向長度;dFL與dFD分別代表葉元體的升力和阻力。

    1.3 葉元參數(shù)

    葉元參數(shù)為葉元體升阻力系數(shù)與攻角的對(duì)應(yīng)關(guān)系,在使用體積力模型前存儲(chǔ)到特定文件中,并由體積力計(jì)算程序?qū)崟r(shí)讀取。本文通過對(duì)敞水試驗(yàn)CFD 模擬結(jié)果中的葉元壁面進(jìn)行壓力和粘性力積分,可以分別得到葉元體的升力和阻力系數(shù),進(jìn)而確定各葉元體升阻力系數(shù)與攻角的對(duì)應(yīng)關(guān)系。在執(zhí)行體積力程序時(shí),葉元性能由分段線性插值得到對(duì)應(yīng)攻角下葉元體的升阻力系數(shù)。具體的實(shí)施流程如下:首先進(jìn)行螺旋槳敞水試驗(yàn)的CFD 數(shù)值模擬,并將螺旋槳壁面沿徑向等距分為多個(gè)葉元,各葉元體對(duì)推力與扭矩的貢獻(xiàn)dT與dQ通過表面壓力積分得到(圖2),進(jìn)而計(jì)算葉元體的升阻力系數(shù),并將所得力系數(shù)與攻角的對(duì)應(yīng)關(guān)系保存為數(shù)據(jù)向量,由體積力計(jì)算程序通過分段線性插值方法讀取。在體積力程序的執(zhí)行過程中,首先通過誘導(dǎo)因子與當(dāng)?shù)厮俣鹊玫饺~元來流,進(jìn)而計(jì)算葉元攻角α,最終通過攻角數(shù)據(jù)向量α與升阻力系數(shù)向量CL、CD經(jīng)分段線性插值得到式(8)與式(9)中的升阻力系數(shù)。分段線性插值格式如式(10)、式(11)所示,此種插值方式可盡量減小數(shù)據(jù)處理造成的誤差。

    圖2 螺旋槳葉元體Fig.2 Blade element of propeller

    1.4 真實(shí)螺旋槳誘導(dǎo)因子分布

    誘導(dǎo)因子體現(xiàn)了來流速度與當(dāng)?shù)厮俣鹊年P(guān)系,如圖3 所示。傳統(tǒng)葉素動(dòng)量理論中誘導(dǎo)因子的迭代求解基于理想流體理論[2,4],不適用于粘性求解器中的體積力模型,因此需要一種方法來研究粘性流場(chǎng)螺旋槳載荷作用下各葉元體來流與當(dāng)?shù)厮俣鹊年P(guān)系。為得到這種關(guān)系,本文利用1.3節(jié)得到的葉元參數(shù)得到了多組來流速度下與真實(shí)螺旋槳模型載荷分布相符的體積力分布,并將此體積力分布作用于敞水流場(chǎng)中得到誘導(dǎo)速度場(chǎng),如圖4所示。在得到誘導(dǎo)速度場(chǎng)的過程中,葉元體的水動(dòng)力大小不會(huì)隨流場(chǎng)的變化而迭代。在此種情況下,最終得到的誘導(dǎo)速度場(chǎng)與真實(shí)螺旋槳的誘導(dǎo)速度場(chǎng)吻合程度較高。各葉元體的當(dāng)?shù)厮俣扔捎?jì)算穩(wěn)定后的速度場(chǎng)在致動(dòng)點(diǎn)位置處插值得到,并由同一徑向位置處各致動(dòng)點(diǎn)當(dāng)?shù)厮俣鹊钠骄Y(jié)果代表此徑向位置處的葉元體當(dāng)?shù)厮俣?。最終由預(yù)設(shè)來流速度V與各徑向位置處當(dāng)?shù)厮俣萔x計(jì)算各徑向位置處的誘導(dǎo)因子大小。由于敞水試驗(yàn)中橫向流動(dòng)較弱,因此本研究中尚未考慮切向誘導(dǎo)因子。

    圖3 理想流體理論[2]Fig.3 Ideal fluid theory[2]

    圖4 基于真實(shí)螺旋槳載荷分布的誘導(dǎo)速度場(chǎng)Fig.4 Induced velocity field based on real propeller load distribution

    經(jīng)由上述過程,可以建立單個(gè)工況下誘導(dǎo)因子與徑向位置的關(guān)系,也可建立同一徑向位置處誘導(dǎo)因子與當(dāng)?shù)厮俣鹊年P(guān)系。將各徑向位置處誘導(dǎo)因子數(shù)據(jù)記為向量a,當(dāng)?shù)厮俣扔洖橄蛄縑x并存儲(chǔ)到特定文件中。在體積力模型計(jì)算過程中,通過插值得到致動(dòng)點(diǎn)位置處的當(dāng)?shù)厮俣群?,通過在此數(shù)據(jù)向量中進(jìn)行分段線性插值可得到各徑向位置處的誘導(dǎo)因子,如式(12)所示。由于建立了分段線性插值關(guān)系,在流場(chǎng)迭代求解中插值數(shù)據(jù)也將經(jīng)過反復(fù)迭代,體積力模型能否通過迭代收斂到螺旋槳CFD數(shù)值模擬的載荷分布是本文提出的體積力模型是否具有應(yīng)用性的關(guān)鍵。

    2 算例設(shè)置

    2.1 模型與工況

    本文進(jìn)行了KP505 螺旋槳真實(shí)模型及體積力模型的敞水試驗(yàn)數(shù)值模擬。圖5 為KP505 螺旋槳幾何模型。真實(shí)螺旋槳模型敞水試驗(yàn)轉(zhuǎn)速為10 r/s,包含0.4、0.5、0.6、0.7、0.8 五個(gè)進(jìn)速系數(shù)。擬合誘導(dǎo)因子的體積力模型工況參數(shù)與之相同。KP505螺旋槳主要尺度及擬合誘導(dǎo)因子工況參數(shù)如表1所示,體積力模型直徑及轂徑與之相同。

    表1 KP505螺旋槳主要參數(shù)及敞水試驗(yàn)工況參數(shù)Tab.1 Main parameters of KP505 propeller and open-water test

    圖5 KP505螺旋槳Fig.5 Geometry of KP505 propeller

    2.2 計(jì)算域及網(wǎng)格劃分

    圖6為真實(shí)螺旋槳模型與體積力模型敞水試驗(yàn)?zāi)M計(jì)算域大小。真實(shí)螺旋槳模型及體積力模型的計(jì)算域大小相同,均為圓柱形計(jì)算域,直徑為5D,槳盤面至入口為3D、至出口為5D。圖7 為螺旋槳敞水試驗(yàn)網(wǎng)格劃分,真實(shí)螺旋槳模型算例網(wǎng)格量為251 萬,時(shí)間步長取0.0005 s。體積力模型算例網(wǎng)格量為46 萬,時(shí)間步長取0.001 s。真實(shí)螺旋槳模型采用滑移網(wǎng)格方法,滑移面附近及螺旋槳壁面均進(jìn)行逐級(jí)加密。圖7(b)中圓柱區(qū)域?yàn)轶w積力作用區(qū)域,徑向網(wǎng)格數(shù)量為38 個(gè),周向網(wǎng)格數(shù)量為104個(gè)。網(wǎng)格無關(guān)性驗(yàn)證與時(shí)間步長無關(guān)性驗(yàn)證已由之前相關(guān)工作完成[17]。

    圖6 KP505螺旋槳敞水試驗(yàn)?zāi)M計(jì)算域Fig.6 Computational domain of KP505 propeller open-water test simulation

    圖7 KP505螺旋槳敞水試驗(yàn)?zāi)M網(wǎng)格Fig.7 Grids of KP505 propeller open-water test simulation

    3 結(jié)果分析

    圖8為改進(jìn)的體積力模型與真實(shí)螺旋槳模型敞水曲線模擬結(jié)果,其中real、bf與exp分別代表真實(shí)螺旋槳模型和體積力模型的計(jì)算結(jié)果以及螺旋槳敞水性能的物理實(shí)驗(yàn)結(jié)果??梢园l(fā)現(xiàn)三者有較好的吻合程度,表明真實(shí)螺旋槳的CFD具有準(zhǔn)確性。同時(shí)表明在應(yīng)用了符合真實(shí)螺旋槳模型的葉元性能與誘導(dǎo)因子分布后,在不使用翼型理論經(jīng)驗(yàn)公式的情況下,體積力模型的準(zhǔn)確度得到了有效提高。圖9 展示了各進(jìn)速條件下體積力模型相對(duì)于真實(shí)螺旋槳敞水曲線計(jì)算結(jié)果的誤差,除J=0.4 時(shí)推力系數(shù)誤差為1.05%外,各進(jìn)速條件下推力系數(shù)誤差與扭矩系數(shù)誤差均在1%以內(nèi),其中扭矩系數(shù)誤差均在0.8%以內(nèi),且誤差水平未隨進(jìn)速系數(shù)出現(xiàn)較大變化,表明改進(jìn)后的體積力模型誤差水平不再與進(jìn)速系數(shù)敏感[4,7,11],其準(zhǔn)確性可適用于不同工況。兩種模型模擬條件下槳載荷在其中三種進(jìn)速條件下的徑向分布如圖10 所示。其中axial 與tangential 分別代表軸向載荷與切向載荷,并進(jìn)行了無因次化處理。兩種模型數(shù)據(jù)的吻合程度較高,表明改進(jìn)后的葉素動(dòng)量理論模型可以通過不同徑向位置處的當(dāng)?shù)厮俣鹊玫椒险鎸?shí)螺旋槳載荷分布的葉元水動(dòng)力穩(wěn)定值。

    圖8 敞水曲線模擬結(jié)果Fig.8 Open-water-curve simulation results

    圖9 體積力模型敞水曲線誤差分布Fig.9 Error distribution of open-water-curve of body-force model

    圖10 螺旋槳載荷分布模擬結(jié)果Fig.10 Simulation results of propeller load distribution

    圖11 與圖12 分別為J=0.5 與J=0.7 條件下兩種模型的尾流場(chǎng)軸向速度對(duì)比,其中圖12 為槳盤面后0.5D、1D與1.5D處的位置數(shù)據(jù),此分布范圍囊括了一般槳舵干擾工況下舵模型的位置[17]。真實(shí)螺旋槳模型的尾流場(chǎng)數(shù)據(jù)為一個(gè)旋轉(zhuǎn)周期內(nèi)的時(shí)均數(shù)據(jù)。誘導(dǎo)速度峰值分布在0.6R~0.8R范圍內(nèi),與圖10中最大槳載荷位置相同。可以發(fā)現(xiàn)誘導(dǎo)速度分布峰值附近的較大范圍內(nèi)兩種模型的速度分布在兩種工況下的吻合程度均較高,表明體積力模型可以模擬真實(shí)螺旋槳尾流不同截面中的動(dòng)量輸運(yùn)。在輪轂處兩組數(shù)據(jù)存在一定差距的原因?yàn)樵谶M(jìn)行體積力模型模擬時(shí)未考慮槳軸及輪轂壁面的建模,因此忽略了真實(shí)螺旋槳模型槳軸與輪轂的阻塞效應(yīng),導(dǎo)致體積力模型誘導(dǎo)速度場(chǎng)輪轂附近速度偏高。此差異隨流場(chǎng)的向后發(fā)展逐漸變小,且對(duì)于螺旋槳尾流動(dòng)量輸運(yùn)不起主要影響作用。同時(shí)從圖10中可以發(fā)現(xiàn),槳轂附近的葉元體載荷為整個(gè)槳葉中的最低水平,對(duì)螺旋槳性能不起主要作用,因此槳轂存在對(duì)螺旋槳性能的影響可以忽略。兩種模型尾流分布存在細(xì)微差異的另一原因是真實(shí)螺旋槳壁面的擾動(dòng)作用產(chǎn)生了大量復(fù)雜渦結(jié)構(gòu),而體積力模型由于忽略了葉片旋轉(zhuǎn)效應(yīng),渦結(jié)構(gòu)更加平滑有序。圖13為兩種模型計(jì)算得到的渦結(jié)構(gòu)(Q=50等值面)對(duì)比。體積力模型無法還原較為精細(xì)的流場(chǎng)信息,但可以通過替代真實(shí)螺旋槳模型有效減少螺旋槳尾流動(dòng)量輸運(yùn)模擬的計(jì)算成本。

    圖11 尾流場(chǎng)模擬結(jié)果Fig.11 Simulation results of wake flow field

    圖13 體積力模型與真實(shí)螺旋槳模型的渦結(jié)構(gòu)Fig.13 Vortex structure of body-force model and real propeller model

    4 結(jié) 語

    本文通過擬合誘導(dǎo)因子的方式,將誘導(dǎo)因子與當(dāng)?shù)厮俣鹊膶?duì)應(yīng)數(shù)據(jù)存儲(chǔ)到葉素動(dòng)量理論的執(zhí)行過程中,使用螺旋槳敞水CFD數(shù)值模擬計(jì)算葉元性能,并將改進(jìn)后的體積力模型應(yīng)用于五組進(jìn)速下的KP505螺旋槳敞水試驗(yàn)數(shù)值模擬?;诟倪M(jìn)葉素動(dòng)量理論體積力模型模擬的KP505螺旋槳敞水試驗(yàn)結(jié)果表明,體積力模型可以從初始流場(chǎng)迭代到螺旋槳CFD 數(shù)值模擬的工作狀態(tài),其螺旋槳性能、載荷分布以及尾流分布均與螺旋槳直接數(shù)值模擬的結(jié)果相近,且不同進(jìn)速下均有較好的預(yù)報(bào)精度,敞水誤差在五種進(jìn)速下均未超過1.2%。

    以上研究表明,基于符合真實(shí)螺旋槳誘導(dǎo)因子分布改進(jìn)后的葉素動(dòng)量理論模型,適用于不同進(jìn)速下的螺旋槳敞水性能預(yù)報(bào),克服了以往葉素動(dòng)量理論體積力模型的準(zhǔn)確性隨模擬工況變化敏感度高的問題。在不均勻來流問題如船槳干擾問題中,槳盤面各處的來流速度不同,特別是在船舶操縱條件下螺旋槳的入流條件更加復(fù)雜,因此在較大進(jìn)速范圍內(nèi)保證葉素動(dòng)量理論體積力模型的計(jì)算準(zhǔn)確性對(duì)于解決此類問題具有較大的意義。本文提出的改進(jìn)體積力模型預(yù)報(bào)得到的螺旋槳的誘導(dǎo)速度場(chǎng)與真實(shí)螺旋槳吻合程度較好,因此可以較為準(zhǔn)確地代替真實(shí)螺旋槳進(jìn)行槳與其他結(jié)構(gòu)物的相互作用模擬,在降低計(jì)算時(shí)間與存儲(chǔ)空間需求基礎(chǔ)上還能較好地提升計(jì)算精度,為未來船槳舵干擾下的水動(dòng)力性能高效評(píng)估提供了必要的技術(shù)支撐。

    猜你喜歡
    葉素動(dòng)量螺旋槳
    動(dòng)量守恒定律在三個(gè)物體系中的應(yīng)用
    應(yīng)用動(dòng)量守恒定律解題之秘訣
    基于CFD的螺旋槳拉力確定方法
    動(dòng)量相關(guān)知識(shí)的理解和應(yīng)用
    不能說的秘密
    精心設(shè)計(jì)老婆出軌釀滅門慘案
    新傳奇(2016年29期)2016-09-18 08:30:12
    那個(gè)不能說的秘密
    中外文摘(2016年12期)2016-08-29 08:14:40
    瘋狂的秘密:離婚5天后滅門岳父家
    3800DWT加油船螺旋槳諧鳴分析及消除方法
    廣東造船(2015年6期)2015-02-27 10:52:46
    螺旋槳轂帽鰭節(jié)能性能的數(shù)值模擬
    在线天堂中文资源库| 国产熟女xx| 欧美日韩国产mv在线观看视频| 日日摸夜夜添夜夜添小说| 99国产精品99久久久久| 久久久久久久久中文| 亚洲 欧美一区二区三区| 老司机午夜福利在线观看视频| 1024香蕉在线观看| 欧美日韩黄片免| 欧美日韩福利视频一区二区| 欧美日韩福利视频一区二区| 欧美一级毛片孕妇| 精品一区二区三区av网在线观看| 一级黄色大片毛片| 好男人电影高清在线观看| 日本欧美视频一区| 亚洲精品美女久久久久99蜜臀| 中文字幕另类日韩欧美亚洲嫩草| 日本vs欧美在线观看视频| 欧美另类亚洲清纯唯美| 久久青草综合色| 一本大道久久a久久精品| 又紧又爽又黄一区二区| 日韩 欧美 亚洲 中文字幕| 99在线人妻在线中文字幕| 91av网站免费观看| 亚洲五月婷婷丁香| 国产精品亚洲av一区麻豆| 久久久久久久久中文| av免费在线观看网站| 少妇粗大呻吟视频| 一a级毛片在线观看| www.999成人在线观看| 亚洲视频免费观看视频| 久久精品91无色码中文字幕| 国产麻豆69| 色综合欧美亚洲国产小说| 日本vs欧美在线观看视频| 午夜成年电影在线免费观看| 久久青草综合色| 女人爽到高潮嗷嗷叫在线视频| 老司机亚洲免费影院| 欧美日本亚洲视频在线播放| 婷婷六月久久综合丁香| 欧美黄色片欧美黄色片| 久久人妻熟女aⅴ| 国产av精品麻豆| 久久精品国产亚洲av高清一级| 亚洲国产精品一区二区三区在线| 欧美日韩乱码在线| 香蕉丝袜av| 男女下面插进去视频免费观看| 两性夫妻黄色片| 制服诱惑二区| 精品久久蜜臀av无| 欧美大码av| 咕卡用的链子| 又紧又爽又黄一区二区| 午夜福利一区二区在线看| 国产免费现黄频在线看| 久久久久国内视频| 在线观看日韩欧美| 一级作爱视频免费观看| 男女高潮啪啪啪动态图| 精品久久久精品久久久| 一边摸一边做爽爽视频免费| 久久国产乱子伦精品免费另类| 亚洲,欧美精品.| 一夜夜www| 99香蕉大伊视频| 老司机福利观看| 脱女人内裤的视频| 国产精品国产高清国产av| 999久久久国产精品视频| 一区福利在线观看| 日韩视频一区二区在线观看| 国产精华一区二区三区| 免费av中文字幕在线| 亚洲欧美一区二区三区久久| 久久精品亚洲精品国产色婷小说| 12—13女人毛片做爰片一| 99re在线观看精品视频| 国产免费现黄频在线看| 首页视频小说图片口味搜索| 亚洲精品av麻豆狂野| 老司机在亚洲福利影院| 亚洲专区中文字幕在线| 视频区欧美日本亚洲| 免费人成视频x8x8入口观看| 亚洲人成伊人成综合网2020| 黄片播放在线免费| 色婷婷久久久亚洲欧美| 日韩精品免费视频一区二区三区| 香蕉丝袜av| 亚洲人成电影免费在线| 伦理电影免费视频| 色老头精品视频在线观看| 脱女人内裤的视频| 亚洲午夜理论影院| 日本免费一区二区三区高清不卡 | 亚洲第一av免费看| 亚洲av成人av| 18美女黄网站色大片免费观看| 一进一出抽搐动态| 亚洲在线自拍视频| 午夜福利影视在线免费观看| 亚洲精品在线美女| 在线永久观看黄色视频| 亚洲五月婷婷丁香| 色婷婷av一区二区三区视频| 99精品在免费线老司机午夜| 精品一区二区三区视频在线观看免费 | 欧美一级毛片孕妇| 日日摸夜夜添夜夜添小说| 女同久久另类99精品国产91| www.自偷自拍.com| 99久久久亚洲精品蜜臀av| 欧美一区二区精品小视频在线| 91麻豆精品激情在线观看国产 | 亚洲自偷自拍图片 自拍| 麻豆成人av在线观看| 国产精品久久久av美女十八| 亚洲人成网站在线播放欧美日韩| 在线十欧美十亚洲十日本专区| 亚洲熟女毛片儿| 亚洲一区二区三区色噜噜 | 精品国产乱码久久久久久男人| 久久久久久大精品| 天天添夜夜摸| 精品一区二区三区四区五区乱码| 天堂影院成人在线观看| 99久久久亚洲精品蜜臀av| 国产成年人精品一区二区 | 每晚都被弄得嗷嗷叫到高潮| 久久天堂一区二区三区四区| 1024视频免费在线观看| av天堂久久9| 水蜜桃什么品种好| 国产精品国产高清国产av| 亚洲国产毛片av蜜桃av| 嫩草影院精品99| 视频在线观看一区二区三区| 精品久久久久久电影网| 亚洲av成人一区二区三| 99国产极品粉嫩在线观看| 国产成年人精品一区二区 | 国产成人免费无遮挡视频| 欧美+亚洲+日韩+国产| 国产精品一区二区三区四区久久 | 不卡av一区二区三区| 老司机在亚洲福利影院| 美女高潮喷水抽搐中文字幕| 中国美女看黄片| 校园春色视频在线观看| 成年版毛片免费区| 日本五十路高清| 啦啦啦 在线观看视频| 亚洲七黄色美女视频| 免费在线观看黄色视频的| 最近最新中文字幕大全电影3 | 日本a在线网址| 一进一出抽搐动态| 脱女人内裤的视频| 纯流量卡能插随身wifi吗| 亚洲成av片中文字幕在线观看| 日韩大码丰满熟妇| 久久人人爽av亚洲精品天堂| 欧美成人午夜精品| 啦啦啦免费观看视频1| 99国产综合亚洲精品| 99在线视频只有这里精品首页| 国产成人啪精品午夜网站| 亚洲情色 制服丝袜| 精品卡一卡二卡四卡免费| 欧美在线黄色| 一区在线观看完整版| 高潮久久久久久久久久久不卡| 乱人伦中国视频| 国产一区二区三区在线臀色熟女 | 变态另类成人亚洲欧美熟女 | 久久久水蜜桃国产精品网| 国产色视频综合| 免费观看精品视频网站| 男人舔女人的私密视频| 亚洲av日韩精品久久久久久密| 99久久国产精品久久久| 久久精品国产综合久久久| 欧美+亚洲+日韩+国产| 90打野战视频偷拍视频| 正在播放国产对白刺激| 亚洲成人免费电影在线观看| 午夜免费观看网址| 99久久精品国产亚洲精品| 搡老熟女国产l中国老女人| 国产欧美日韩一区二区三区在线| 欧美在线一区亚洲| 国产成人啪精品午夜网站| 九色亚洲精品在线播放| 久久久国产精品麻豆| 一进一出抽搐动态| 亚洲av成人一区二区三| 丁香欧美五月| 国产精品久久久久久人妻精品电影| 国产精品一区二区免费欧美| 一区二区三区激情视频| 欧美激情高清一区二区三区| 国产成人欧美在线观看| 在线观看66精品国产| 黄色女人牲交| 欧美日韩视频精品一区| 欧美日韩av久久| 亚洲精品久久成人aⅴ小说| 欧美日韩福利视频一区二区| 咕卡用的链子| avwww免费| 99香蕉大伊视频| 看黄色毛片网站| 欧美黑人精品巨大| 国产成人影院久久av| 真人做人爱边吃奶动态| 亚洲国产精品999在线| 岛国视频午夜一区免费看| 99久久精品国产亚洲精品| 中文字幕精品免费在线观看视频| 亚洲成av片中文字幕在线观看| 一区二区三区精品91| 国产精品一区二区免费欧美| 色综合婷婷激情| 国产三级在线视频| 欧美黑人欧美精品刺激| 极品人妻少妇av视频| 老司机亚洲免费影院| 国产主播在线观看一区二区| 亚洲人成伊人成综合网2020| 国产亚洲精品综合一区在线观看 | cao死你这个sao货| 国产高清国产精品国产三级| 老司机靠b影院| 高清黄色对白视频在线免费看| 国内久久婷婷六月综合欲色啪| 午夜两性在线视频| 欧美日本中文国产一区发布| 夫妻午夜视频| 黄色 视频免费看| 亚洲久久久国产精品| 一级,二级,三级黄色视频| 性色av乱码一区二区三区2| 欧美成人午夜精品| 亚洲国产精品999在线| videosex国产| 久久九九热精品免费| 亚洲国产中文字幕在线视频| netflix在线观看网站| 丝袜在线中文字幕| 国产午夜精品久久久久久| 国产精品偷伦视频观看了| 美女福利国产在线| 日韩欧美国产一区二区入口| 午夜a级毛片| 天堂中文最新版在线下载| av片东京热男人的天堂| 亚洲色图综合在线观看| 国产伦一二天堂av在线观看| 国产成+人综合+亚洲专区| 黑人巨大精品欧美一区二区蜜桃| 日韩精品免费视频一区二区三区| 90打野战视频偷拍视频| 欧美丝袜亚洲另类 | 真人做人爱边吃奶动态| 不卡av一区二区三区| 久久久久亚洲av毛片大全| 日本精品一区二区三区蜜桃| 免费不卡黄色视频| 女性被躁到高潮视频| 亚洲中文av在线| 91麻豆av在线| 欧美亚洲日本最大视频资源| 嫁个100分男人电影在线观看| 精品国内亚洲2022精品成人| 啦啦啦免费观看视频1| 一进一出抽搐动态| 欧美精品啪啪一区二区三区| 欧美激情 高清一区二区三区| 麻豆国产av国片精品| www日本在线高清视频| 国产精品久久久av美女十八| 亚洲人成伊人成综合网2020| 变态另类成人亚洲欧美熟女 | 日本免费a在线| 99久久精品国产亚洲精品| 国产精品久久久久成人av| 无遮挡黄片免费观看| 亚洲五月婷婷丁香| 老司机深夜福利视频在线观看| 男女做爰动态图高潮gif福利片 | 亚洲av日韩精品久久久久久密| 成人影院久久| 99国产极品粉嫩在线观看| 久久天躁狠狠躁夜夜2o2o| 欧美一级毛片孕妇| 又黄又粗又硬又大视频| 999久久久国产精品视频| 久久天堂一区二区三区四区| 亚洲一区二区三区欧美精品| 丰满人妻熟妇乱又伦精品不卡| 日韩大尺度精品在线看网址 | 日本五十路高清| 日韩精品免费视频一区二区三区| 宅男免费午夜| 成在线人永久免费视频| 一区在线观看完整版| 女人精品久久久久毛片| 亚洲在线自拍视频| 国产精品99久久99久久久不卡| 午夜视频精品福利| 丰满的人妻完整版| 精品免费久久久久久久清纯| 男人的好看免费观看在线视频 | 交换朋友夫妻互换小说| 极品人妻少妇av视频| 黄片播放在线免费| 色婷婷久久久亚洲欧美| 97碰自拍视频| 18禁美女被吸乳视频| 国产97色在线日韩免费| 不卡一级毛片| 18禁美女被吸乳视频| 国产熟女xx| a级片在线免费高清观看视频| 午夜老司机福利片| 日本免费一区二区三区高清不卡 | 欧美乱妇无乱码| 91在线观看av| 国产成人精品在线电影| 正在播放国产对白刺激| 天天躁狠狠躁夜夜躁狠狠躁| 美女高潮到喷水免费观看| 精品久久久久久成人av| 国产欧美日韩一区二区三区在线| 99在线视频只有这里精品首页| 久久精品亚洲精品国产色婷小说| 老汉色av国产亚洲站长工具| 啪啪无遮挡十八禁网站| 久久久久久久久免费视频了| 欧美黄色片欧美黄色片| 亚洲欧洲精品一区二区精品久久久| 激情视频va一区二区三区| 大型av网站在线播放| 色尼玛亚洲综合影院| 欧美人与性动交α欧美软件| 亚洲熟妇中文字幕五十中出 | 色婷婷久久久亚洲欧美| 天堂中文最新版在线下载| 日韩大尺度精品在线看网址 | 精品国产美女av久久久久小说| 制服诱惑二区| 国产乱人伦免费视频| 宅男免费午夜| 亚洲欧美激情综合另类| 国产男靠女视频免费网站| 嫁个100分男人电影在线观看| 国产熟女xx| av有码第一页| www.精华液| 国产深夜福利视频在线观看| 免费搜索国产男女视频| 国产精品久久久久成人av| 日韩三级视频一区二区三区| 大码成人一级视频| 久久欧美精品欧美久久欧美| 最近最新免费中文字幕在线| 制服人妻中文乱码| 免费在线观看影片大全网站| 日韩 欧美 亚洲 中文字幕| 精品国产美女av久久久久小说| 欧美最黄视频在线播放免费 | www.自偷自拍.com| 在线观看午夜福利视频| 国产片内射在线| 亚洲一区二区三区色噜噜 | 精品福利永久在线观看| 免费久久久久久久精品成人欧美视频| 日韩免费高清中文字幕av| 国产亚洲精品第一综合不卡| 91成年电影在线观看| 久久久精品国产亚洲av高清涩受| 免费在线观看日本一区| 国产成人欧美| 亚洲久久久国产精品| 99riav亚洲国产免费| 久久人人爽av亚洲精品天堂| 亚洲国产欧美日韩在线播放| 国产一区二区三区综合在线观看| 如日韩欧美国产精品一区二区三区| 久久久久久久久免费视频了| 激情视频va一区二区三区| 天堂动漫精品| 黄色片一级片一级黄色片| 国产精品电影一区二区三区| 日本精品一区二区三区蜜桃| 欧洲精品卡2卡3卡4卡5卡区| 搡老岳熟女国产| 国产精品秋霞免费鲁丝片| 国产三级在线视频| 18禁黄网站禁片午夜丰满| 美女国产高潮福利片在线看| 亚洲第一av免费看| 亚洲视频免费观看视频| 无遮挡黄片免费观看| 亚洲九九香蕉| av片东京热男人的天堂| 日韩免费高清中文字幕av| 天天躁夜夜躁狠狠躁躁| 99在线人妻在线中文字幕| 如日韩欧美国产精品一区二区三区| 久久 成人 亚洲| 久久国产精品男人的天堂亚洲| 99国产精品一区二区蜜桃av| 大码成人一级视频| 变态另类成人亚洲欧美熟女 | 欧美日韩一级在线毛片| 精品卡一卡二卡四卡免费| 亚洲第一av免费看| www.www免费av| 一边摸一边抽搐一进一出视频| 岛国视频午夜一区免费看| 午夜成年电影在线免费观看| 夜夜看夜夜爽夜夜摸 | 午夜福利,免费看| 午夜精品久久久久久毛片777| 欧美人与性动交α欧美软件| 欧美日本亚洲视频在线播放| 电影成人av| 人人妻人人添人人爽欧美一区卜| 韩国精品一区二区三区| 亚洲欧美一区二区三区黑人| 天堂影院成人在线观看| 黑人巨大精品欧美一区二区蜜桃| 久久久精品国产亚洲av高清涩受| 看黄色毛片网站| 91成人精品电影| 日韩大码丰满熟妇| 久久国产亚洲av麻豆专区| 国产精品久久久av美女十八| 久久亚洲真实| 91老司机精品| 欧美日本亚洲视频在线播放| 国产精品 国内视频| 免费日韩欧美在线观看| 99精国产麻豆久久婷婷| 日日摸夜夜添夜夜添小说| 99国产精品免费福利视频| 校园春色视频在线观看| 成人三级黄色视频| 久久这里只有精品19| 水蜜桃什么品种好| 国产成人精品无人区| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜免费观看网址| 老司机福利观看| 国产成人系列免费观看| 午夜老司机福利片| 欧美亚洲日本最大视频资源| 桃红色精品国产亚洲av| 久久国产精品男人的天堂亚洲| 久久亚洲真实| 高清在线国产一区| 麻豆成人av在线观看| 日本精品一区二区三区蜜桃| bbb黄色大片| 岛国在线观看网站| 久久精品国产清高在天天线| www.www免费av| 中文字幕最新亚洲高清| ponron亚洲| 国产精品自产拍在线观看55亚洲| 首页视频小说图片口味搜索| 亚洲 欧美 日韩 在线 免费| 亚洲人成77777在线视频| 中文字幕最新亚洲高清| 妹子高潮喷水视频| 亚洲精品中文字幕一二三四区| 熟女少妇亚洲综合色aaa.| 三上悠亚av全集在线观看| 久久精品成人免费网站| 99国产精品99久久久久| 亚洲国产看品久久| 99在线视频只有这里精品首页| 国产免费男女视频| 中文字幕人妻丝袜一区二区| 黄色怎么调成土黄色| 亚洲精品久久成人aⅴ小说| 99在线视频只有这里精品首页| 曰老女人黄片| 日韩视频一区二区在线观看| 国产aⅴ精品一区二区三区波| 黑人巨大精品欧美一区二区蜜桃| 久久久久国内视频| 啦啦啦在线免费观看视频4| av在线播放免费不卡| 18禁黄网站禁片午夜丰满| 一边摸一边做爽爽视频免费| 久久中文字幕人妻熟女| 欧美大码av| 国产av在哪里看| 桃红色精品国产亚洲av| 久久久久久久久中文| 国产精品一区二区三区四区久久 | 精品国产乱子伦一区二区三区| 中文字幕精品免费在线观看视频| 在线永久观看黄色视频| 操美女的视频在线观看| 精品一区二区三区视频在线观看免费 | 欧美日本亚洲视频在线播放| 1024视频免费在线观看| 夜夜躁狠狠躁天天躁| 午夜影院日韩av| 久久久久国内视频| 精品国产国语对白av| av天堂在线播放| www.www免费av| 欧美 亚洲 国产 日韩一| 欧美黄色片欧美黄色片| 性少妇av在线| 亚洲欧美日韩另类电影网站| 免费观看人在逋| av国产精品久久久久影院| 精品电影一区二区在线| 免费人成视频x8x8入口观看| 91在线观看av| 亚洲激情在线av| 欧美中文日本在线观看视频| 国产亚洲av高清不卡| 国产有黄有色有爽视频| 人人妻人人添人人爽欧美一区卜| 精品国产一区二区三区四区第35| 久久香蕉激情| 一级a爱视频在线免费观看| 五月开心婷婷网| 精品日产1卡2卡| 国产高清视频在线播放一区| www.www免费av| 欧美日韩国产mv在线观看视频| 国产欧美日韩精品亚洲av| 中文字幕最新亚洲高清| 夜夜看夜夜爽夜夜摸 | 亚洲精品成人av观看孕妇| 淫秽高清视频在线观看| 婷婷丁香在线五月| 纯流量卡能插随身wifi吗| 国产精品久久久久久人妻精品电影| 99国产精品99久久久久| 国产精品国产av在线观看| 97超级碰碰碰精品色视频在线观看| 精品福利观看| 成人亚洲精品一区在线观看| svipshipincom国产片| 色哟哟哟哟哟哟| 视频在线观看一区二区三区| 高清欧美精品videossex| 深夜精品福利| 国产亚洲欧美精品永久| 久久精品亚洲熟妇少妇任你| 久久久水蜜桃国产精品网| 多毛熟女@视频| 亚洲国产精品999在线| 后天国语完整版免费观看| 曰老女人黄片| 黑丝袜美女国产一区| 9色porny在线观看| 亚洲精品久久成人aⅴ小说| 国产精品秋霞免费鲁丝片| 黄网站色视频无遮挡免费观看| 欧美日韩乱码在线| 久久久久国产一级毛片高清牌| 欧美一区二区精品小视频在线| 电影成人av| 夜夜躁狠狠躁天天躁| 一级a爱片免费观看的视频| 欧美乱妇无乱码| 激情在线观看视频在线高清| 国产精品国产av在线观看| 久久久国产成人精品二区 | 成人国产一区最新在线观看| 国产免费现黄频在线看| 久久午夜综合久久蜜桃| 啦啦啦免费观看视频1| 国产极品粉嫩免费观看在线| 看黄色毛片网站| 在线免费观看的www视频| 精品欧美一区二区三区在线| 亚洲全国av大片| 国产精品成人在线| 国产精品国产高清国产av| 国产91精品成人一区二区三区| 女人被狂操c到高潮| 男人舔女人的私密视频| 无人区码免费观看不卡| 村上凉子中文字幕在线| 久久久国产一区二区| 宅男免费午夜| 精品国内亚洲2022精品成人| 成人亚洲精品一区在线观看| 久久精品亚洲精品国产色婷小说| 99国产精品一区二区蜜桃av| 亚洲精品成人av观看孕妇| 国产av又大| 国产av在哪里看| 欧美成人免费av一区二区三区| 少妇被粗大的猛进出69影院| 亚洲精品久久午夜乱码| 亚洲久久久国产精品| 中文字幕高清在线视频| 欧美成人性av电影在线观看|