本文在非結構性之三角形與四邊形混合網格上,採用一網格頂點有限體積 上風法,來求解二維非穩態可壓縮尤拉方程式;研究重點在了解不同的時 間積分方法,在穩態流場的殘值收斂性,與非穩態流場的時間準確性。以 找出求解穩態與非穩態流場時,所適用的時間積分方法。本文的空間處理 ,以 Roe的近似里曼解子,使用特徵內差變數的 MUSCL差分法,及Van Albada限制函數等,來得到高精確度之數值解。而在時間積分上,使用了 包含顯式的Runge-Kutta四階時間積分法,隱式ALU法,隱式 Symmetric Sweep Point Gauss-Seidel 法,以及改良型的隱式Symmetric Sweep Point Gauss-Seidel法等;並在求解穩態流場問題時,加入局部時間間距 ,殘值平滑法,Three Sweeps Point Gauss-Seidel及改良型的Three Sweeps Point Gauss-Seidel等法來加速殘值收斂。最後,本文以穩態之 斜震波流場,穩態之壁上的斜震波反射流場,及非穩態之移動正震波流場 等物理問題,來測試本文所用的時間積分方式,在穩態與非穩態流場上的 正確性與有效性,並已由所得之高解析解獲得驗證。
|