理想流體力學(xué)大作業(yè)Hess-Smith方法求附加質(zhì)量

1、理想流體大作業(yè)應(yīng)用hesssmith方法求附加質(zhì)量物理背景當(dāng)船舶在海洋中運動或潛艇在深海中運動時,可以看成是物體在無界流中的運動,其所受到的力和力矩可以用物體的速度,加速度和附加質(zhì)量決定,因此,附加質(zhì)量是物體最重要的水動力特性.對于附加質(zhì)量的求解就顯得非常重要.由于附加慣性力的作用,物體在理想流體中的變速運動相當(dāng)于物體自身質(zhì)量上增加了一個附加質(zhì)量而在真空中運動,換句話說,理想流體增大了物體的慣性,是物體很難加速也難減速.附加質(zhì)量有物體上的單位運動絕對速度勢來決定,與流體密度,物體形狀和運動方向有關(guān),而與時間無關(guān).二理論依據(jù)物體的附加質(zhì)量可由下時表示(2.1(2.1) )(2.2)(2.2)其中

2、,i i 分別是總速度勢的分量,%表示物法線方向在各個方向分量,山 表示物體在 i i 方向運動引起其在 j j 方向的附加質(zhì)量, ,那么,只要求出在物體外表速度勢在法線方向的各個分量與法向量的分量,物體的附加質(zhì)量便可求出.在無界流中,物體的總速度勢可以分解為定常速度.是和擾動速度勢中.其在物面上滿足定解條件ar-ar- -應(yīng)用分布源來模擬物面上的擾動勢,那么有=4詈所以,在物體外表上有速度沿法線方向的分量表達式- -= =2 2 一?一E口由劣,聞由物面條件得到2 23T3To oP2P2+北一*.與=-%11所以只要求出分布源密度,就可以求出物體外表的速度勢.為了計算附加質(zhì)量,要求單位速度

3、勢巾 L L 在物面上的值,其滿足式2.62.6.從而得到分布源的單位速度勢密度滿足的積分方程.只要求解出單位速度勢密度,便可以求得附加質(zhì)量.三數(shù)值模型主要通過在物面上的分布源來模擬物體在無界流中所產(chǎn)生的擾動勢,應(yīng)用數(shù)值離散的方法列代數(shù)方程組求解積分方程,即用離散量來表示連續(xù)變量的方法.首先,將物體外表離散,劃分成由 n n 個四邊形單元近似組成.n ns s= = 3 3i i=1=1 物面上各個單元的法向可由節(jié)點求出.(2.(3)(2.(3)(2.(4)(2.(4)(2.(5)(2.(5)(2.(6)(2.(6)從而,速度勢密度 0 0可解,單位速度勢可求,附加質(zhì)量可以計算得到.四幾何模型

4、對球體,橢球和圓柱進行網(wǎng)格的劃分,應(yīng)用不同數(shù)量的四邊形網(wǎng)格來近似物體的外表,通過編寫程序?qū)ζ溥M行網(wǎng)格的劃分,并按順時針輸出個節(jié)點的坐標(biāo)值,相應(yīng)的程序如圖五計算參數(shù)主要對半徑 r=10r=10 的球殼,短軸長度為 10,10,長短軸比為 5 5 的橢球,直徑為 20,20,柱體長和截面直徑比為 5 5 的圓柱進行計算,對其進行網(wǎng)格的劃分并求其附加質(zhì)量系數(shù).六結(jié)果討論(1)(1)球殼對球殼進行網(wǎng)格劃分,畫出不同網(wǎng)格數(shù)對應(yīng)的計算得到的附加質(zhì)量系數(shù)曲線如圖其次,可以列出求解(j/=(j/=i),i),有圖可以看出,隨著網(wǎng)格數(shù)的增加,附加質(zhì)量系數(shù)趨于常數(shù),且 milmil 和 m33m33 的值是相同的,這主要是由于球殼是高度對稱的原因(2)(2)旋轉(zhuǎn)橢球?qū)πD(zhuǎn)橢球進行網(wǎng)格劃分,畫出不