基于非結(jié)構(gòu)動網(wǎng)格技術(shù)及運動物面粘流場數(shù)值計算方法研究.pdf

1、運動物面粘流場數(shù)值求解技術(shù)已成為現(xiàn)代計算流體力學(xué)研究的熱點之一,而其中的關(guān)鍵技術(shù)就是動網(wǎng)格的生成方法以及基于動網(wǎng)格技術(shù)的數(shù)值求解方法。
　　本文針對非結(jié)構(gòu)運動網(wǎng)格生成技術(shù)研究,并將動網(wǎng)格技術(shù)與非結(jié)構(gòu)同位網(wǎng)格有限體積法結(jié)合對Navier-Stokes方程進行數(shù)值求解,建立了基于非結(jié)構(gòu)動網(wǎng)格的非定常數(shù)值離散計算模型及數(shù)值計算方法,編制了數(shù)值求解程序,并完成了算例驗證。
　　本文發(fā)展了兩種非結(jié)構(gòu)動網(wǎng)格生成技術(shù),一種是在標準彈簧法基

2、礎(chǔ)上進行了改進,提出了一種改進的彈簧近似法,該方法除考慮標準彈簧法中的防止網(wǎng)格節(jié)點相互碰撞,還考慮了網(wǎng)格變形的扭轉(zhuǎn)效應(yīng),空間及邊界效應(yīng),使得網(wǎng)格的變形能力大幅提高。第二種動網(wǎng)格技術(shù)為局部重構(gòu)法,該方法能夠有效的解決彈簧法不能解決的物體大幅值運動問題。本文將Delaunay三角化方法應(yīng)用到運動網(wǎng)格的重構(gòu)中,充分利用Delaunay三角化方法的快速生成效率及接近正三角形的性質(zhì),引入的加密邊界作為“窗口”的方式提高了Delaunay三角化方法

3、的布點能力,并編制了基于改進彈簧法的非結(jié)構(gòu)動網(wǎng)格程序和改進局部重構(gòu)法的動網(wǎng)格生成程序。
　　本文在ALE法的基礎(chǔ)上建立了動網(wǎng)格技術(shù)與SIMPLE算法相結(jié)合的數(shù)值計算模型,給出了適合動網(wǎng)格技術(shù)的,滿足幾何守恒定律的數(shù)值離散方法及邊界條件。動網(wǎng)格技術(shù)中網(wǎng)格信息傳遞方式的精度決定計算精度,本文采用雙線性插值結(jié)合ENO方法計算待插點的流場信息參數(shù),這種插值方法的計算精度達到二階精度。利用FORTRAN語言編寫了基于非結(jié)構(gòu)動網(wǎng)格技術(shù)的數(shù)值計

4、算程序,并研究了圓柱橫向周期性振蕩問題,經(jīng)過與 FLUENT軟件的計算結(jié)果對比,結(jié)果表明本文編制的基于非結(jié)構(gòu)運動網(wǎng)格技術(shù)的數(shù)值計算程序是可靠的,可用于開展相關(guān)運動物體流場的數(shù)值計算模擬。
　　本文利用所編制的程序開展了圓柱柔性變形及圓柱后緣柔性板變形的流場特性研究,研究發(fā)現(xiàn)物體的柔性變形能夠改變原流場的結(jié)構(gòu),使得渦的脫落位置后移,圓柱后緣柔性板的應(yīng)用還可以使渦的脫落變得不再規(guī)律。柔性變形物體的這一特性對于抑制渦激振動有著現(xiàn)實的意義