rGFM方法在二維流-固問題模擬中的應(yīng)用

◎剛蕾徐爽 唐強

rGFM方法在二維流-固問題模擬中的應(yīng)用

◎剛蕾1徐爽 唐強

本文結(jié)合rGFM（real Ghost Fluid Method），Hydro-elasto-plastic固體模型，推廣應(yīng)用于二維流-固界面的模擬計算。通過在界面處構(gòu)造并求解Riemann問題， 得到界面處流體的準確的流動狀態(tài)，以此狀態(tài)來定義界面邊界條件，從而將多介質(zhì)流動問題轉(zhuǎn)化為單介質(zhì)流動問題來求解。二維數(shù)值結(jié)果表明，該方法可以用于流-固界面的數(shù)值計算，界面和其它間斷計算準確，二維問題與精確解吻合較好。

1　方程

考慮二維可壓縮多介質(zhì)流體力學(xué)方程組

p為壓力， E為單位體積的總能，可取為

這里e為單位質(zhì)量的內(nèi)能.

可壓縮氣體和水的狀態(tài)方程可以表示下列統(tǒng)一形式，

這里γ和B為常數(shù)。對于理想氣體，γ=1.4， B=0。對水（Tait方程）γ=7.15，

B=3309..通常情況下流體和固體會具有完全不同的特性，但當(dāng)固體受到很大的外力作用時，能表現(xiàn)出流體的特性. Hydroelasto-plastic body即為具有這個特性的固體模型，狀態(tài)方程為

ρ為單位質(zhì)量體積，

這里m和β分別表示體積模量和一個正常數(shù)，F(xiàn)和μ分別是楊氏系數(shù)和泊松比。

下標a表示在標準大氣壓下。用Huber-Mises標準和Hooke定律來得到剪切應(yīng)力，

為了跟蹤流體的運動界面，建立如下二維LevelSet方程【3】

2.界面邊界條件

對于二維問題，要定義界面邊界條件，需要給定四個變量，分別取為密度、x方向速度、y方向速度和壓力.首先考慮流體1，構(gòu)造Riemann問題如下：任意選擇流體1中緊臨界面的網(wǎng)格點A（圖1），尋找流體2中緊鄰界面與A相對應(yīng)的點B，使得這兩點的法向所成的角最小，由于A點和B點流體的狀態(tài)已知，分別設(shè)其中和分別為A點和B點的法向速度，這樣在界面法線方向定義的Riemann問題的初始條件可以取求解Riemann問題可以得到界面處流體的流動狀態(tài)為

圖1　界面處Riemann問題得構(gòu)造

更新A點的狀態(tài)，由于在界面處法向速度和壓力連續(xù)， A點的法向速度和壓力分別取為切向速度和密度會出現(xiàn)間斷，切向速度仍然取A點原來的切向速度，密度取。如果流體1中緊鄰界面的所有網(wǎng)格點的狀態(tài)都得到了更新， 求解方程可以定義流體1的界面邊界條件。同理定義流體2的邊界條件【1，6】.

3　數(shù)值試驗

二維問題主要針對激波與氣-固界面和液-固界面的相互作用進行數(shù)值計算，主要目的是驗證二維界面處理方法對強間斷問題的有效性，數(shù)值結(jié)果表明，該算法能準確地捕捉到流場內(nèi)的物理現(xiàn)象。由于固體內(nèi)部激波相對較弱，為了清楚地顯示激波的位置，本文用Schlieren圖像來顯示密度間斷，計算公式為：

算例 氣體中激波與鋼柱相互作用問題

圖（a）-（d）給出了激波與界面作用的完整過程。由圖可以看出，激波與界面作用后分解為一個入射激波在固體中傳播，一個反射激波在氣體中傳播；由圖（a）可以看到氣體中的反射激波到達邊界并被反射，圖（b）和（c）可以看到被反射的激波在氣體中傳播；圖（d）中被邊界反射的激波到達界面并與界面相互作用，產(chǎn)生新的入射和反射波.

（作者單位：1.南京航空航天大學(xué)金城學(xué)院；2.中國直升機設(shè)計研究所；3.上海航天局第八設(shè)計部）