三维green公式_一种求解三维时域Green函数的数值方法

一种求解三维时域

Green

函数的数值方法

孙丽萍，王建伟，戴绍仕

【摘

要】

利用三维时域方法计算船舶水动力问题的关键是如何准确高效地处理

时域

Green

函数及其空间导数．基于

Bessel

函数的性质推导出

Green

函数与

其空间导数的关系方程式；由时域

Green

函数的表达式建立其满足的常微分方

程并采用

Runge-Kutta

法求解．利用关系方程式直接求解出

Green

函数的空

间导数值．研究结果表明，该方法不仅保证了数值的精度，而且

CPU

的计算时

间和文件的存储空间都明显减少．该方法为时域计算方法提供了一种新的求解

思路．

【期刊名称】

天津大学学报

【年

(

卷

),

期】

2015(000)007

【总页数】

6

【关键词】

水动力；时域；

Green

函数；

Bessel

函数；关系方程

准确预报波浪的诱导响应和载荷对船舶和海洋结构物设计是一个关键的环

节．在采用三维时域法求解水动力问题时，会涉及到对时域

Green

函数的求

解．可是，由于时域

Green

函数的公式含有无穷积分项，以及其函数本身具有

高频振荡和随时间增幅的特点，使得对于该函数的计算变得非常困难．

国内外的很多学者对于时域

Green

函数的求解方法做了很多的研究．最初

Newman[1]

、

Lin

等

[2]

主要是将

Green

函数的计算区域按时间大小分成小时

间区域和大时间区域．在小时间区域内，将

Green

函数的波动项表示成级数形

式；在大时间区域内，将其表示成渐近展开的形式

.

这种方法的优点是计算较为

准确，缺点是小时间区域和大时间区域的分界线没有给出理论的依据，只能靠