为了更适用于复杂边界区域波浪传播的求解,本书基于Beji和Nadaoka改进的Boussinesq方程,采用非结构三角形网格,建立了可用于近岸波浪传播计算的有限元数值计算模型。模型的时间积分采用四阶Adams-Bashforth-Moulton(ABM)预报-校正法,节点上变量的空间一阶导数采用节点周围单元内变量的一阶导数的加权平均值来确定,有效地处理了方程中变量的三阶空间导数项,使模型可以采用线性单元。模型编程采用了索引存储法,只存储方程中二维数组中的非零元素,可节省计算机的内存。同时,针对与坐标轴斜交的全反射边界,通过引入坐标转换的方法,使得在全反射边界处通过求解切向的动量方程来准确求解该类边界处的速度值。一些算例的结果表明,所建立的基于非结构三角形网格的有限元模型能够适用于非规则区域内的波浪传播计算,并可给出很好的模拟结果。
展开