以高斯幅狀基底函數之無網格方法數值模擬完全非線性水面波

Numerical simulation of fully nonlinear surface waves by meshless method with Gaussian radial basis functions

吳南靖(Nan-Jing Wu)

指導教授 ： 蔡丁貴

國立臺灣大學/工學院/土木工程學研究所/博士(2008年)

https://doi.org/10.6342/NTU.2008.01061

摘要

本論文利用無網格數值方法，建立一個能模擬完全非線性水面波的數值模式。在空間方面，以拉普拉斯方程式(Laplace Equation)基本解(Fundamental Solution) 之線性組合擬合流體速度勢(Velocity Potential)，將前述基本解之中心(Source Point)置於計算邊界以外的位置，僅需由適當少數邊界點上的邊界條件即可進行求解，不需要積分，也不需要處理任何奇異點(Singularity)。在時間方面，利用二階中央差分，將自由液面邊界條件離散成為顯式的差分式，不需要迭代，也不需在時間上做積分，即可推測下個時間步（Time Step）求解所需之自由液面邊界條件。此種離散方式亦稱為跳蛙法（Leap Frog Scheme）。而推算下個時間步自由液面邊界條件需要自由液面梯度，則透過高斯幅狀基底函數(Gaussian Radial Basis Function)來擬合自由液面的水位高程，進而得到求解所需之條件。
本論文先模擬孤立波在三維渠道內之傳遞，並檢查質量與能量之守?琚A來證明本數值方法之正確性。然後，再分別模擬二維及三維的問題，來說明本數值方法之適用性。各項數值模擬，均與前人研究之解析解、數值解或實驗結果比較，以說明模擬結果之準確可靠。

關鍵字

無網格數值方法 ； 非線性水波 ； 基本解 ； 高斯輻狀基底函數

並列摘要

A numerical model for fully nonlinear free surface waves is developed in present study, by applying a boundary-type meshless approach with leap frog time marching scheme. Adopting Gaussian Radial Basis Functions to fit the free surface, a non-iterative approach to discretize the nonlinear free surface boundary conditions is formulated. Using the fundamental solutions of the Laplace equation as the solution form of the velocity potential, free-surface wave problems can be solved by collocations at only a few boundary points since the governing equation is automatically satisfied.
The accuracy of present method is verified by comparing the simulated propagation of a solitary wave in a three dimensional flume with an exact solution. The applicability of present model is illustrated by employing it to simulate both two dimensional and three dimensional problems. Good agreement can be found in comparing the results of present numerical model with other schemes, as well as with analytic solutions and available experimental data.

並列關鍵字

Meshless Numerical Method ； Nonlinear Water Wave ； Gaussian Radial Basis Function ； Fundamental Solution

參考文獻

Alves C.J.S., Antunes P.R.S., The method of fundamental solutions applied to the calculation of eigenfrequencies and eigenmodes of 2D simply connected shapes, Computers, Materials and Continua, 2005; 2: 251-266.

Ambrosi D., Quartapelle L., A Taylor-Galerkin Method for Simulation Nonlinear Dispersive Waves, Journal of Computational Physics, 1998; 146: 546-569.

Ata R., Soulaimani A., A stabilized SPH method for inviscid shallow water flows, International Journal for Numerical Methods in Fluids, 2005; 47: 139-159.

Berkhoff J.C.W., Computation of combined refraction-diffraction, Proceedings of 13th Conference in Engineering, ASCE, 1972: 471-490.

Beylkin G., Coifman R., Rokhlin V., Fast wavelet transforms and numerical algorithms I, Communications on Pure and Applied Mathematics, 1991; 44: 141-183.