rGFM方法在二維流-固問題模擬中的應用

◎剛蕾徐爽 唐強

rGFM方法在二維流-固問題模擬中的應用

◎剛蕾1徐爽 唐強

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

1　方程

考慮二維可壓縮多介質流體力學方程組

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

這里e為單位質量的內能.

可壓縮氣體和水的狀態方程可以表示下列統一形式，

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

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

ρ為單位質量體積，

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

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

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

2.界面邊界條件

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

圖1　界面處Riemann問題得構造

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

3　數值試驗

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

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

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

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