非靜水壓模式模擬自由液面波之發展

本研究之目的為發展一套高精度的三維非靜壓sigma模式模擬自由液面波由深水向淺水傳遞之特性。本模式採用隱式有限差分法同時求解非穩態的納維-史托克斯方程式(Navier-Stokes equation)及自由液面邊界條件。本研究提出高精度的表層非靜水壓處理法使得本模式可精確地模擬延散波(dispersive wave)的傳遞。遇到非線性波時,使用sigma座標轉換在計算壓力梯度上會產生極大的數值誤差。因此,本模式採用四階精度空間離散方法來處理水平壓力梯度以克服數值誤差。直接求解三維的水波問題會產生很高的數值計算量,故本模式將一個三維的問題分解成一系列的垂直二維問題。於每個垂直二維問題中,所有的未知數可以水平速度來代換使得整個系統變成水平速度的塊狀七對角矩陣。此種矩陣系統,不需要迭代求解,只需使用解矩陣的副程式(如:雙掃法)就可以直接有效地求解。本研究首先檢視此非靜壓模式的特性包含延散性與非線性。接著,將此模式應用到一系列自由液面波的問題,例如:波流共存、非線性深水波群與流過不規則底床的沿岸波。水波與水流、水波與水波以及水波與底床間的交互作用亦被仔細地探討。本研究證實非靜壓模式利用少許的垂直格網(如:二至五層)能精確有效率地模擬自由液面波問題。

作者:楊智傑