第
34
卷第
11
期
2013
年
11
月
哈尔滨工程大学学报
Vol.34
No
.11
Nov.2013
Journal of
Harbin
Engineering University
预条件
GMRES(
m
)算法在大型浮体
水动力边界元分析中的应用
段文洋
1
,刁峰
2
,陈纪康
1
(
1.
哈尔滨工程大学船舶工程学院,黑龙江哈尔滨
150001
;2.
中国船舶科学研究中心上海分部,上海
200011)
摘
要:针对大型浮体水动力边界元分析产生的复系数线’性方程组结构复杂,直接方法难以求解或者求解费时的问题,
提出一种带有预条件技术的重启动型
GMRES
算法.选取
2
种不同的预条件处理技术对方程组系数矩阵进行预处理,通
过具体算例给出
2
种预条件方法的数值比较.数值试验表明对于大型浮体水动力边界元分析产生的复系数线性方程
组,带有不完全
LU
分解预条件处理技术的
GMRES(m
)算法求解效率最高,优于直接解法.
关键词:大型浮体;边界元分析;重启动
GMRES
;不完全
LU
分解预条件
doi: 10.3969/j.issn.1006-7043.201211097
网络出版地址:
http://wv
仰
.cnki.neι
/
kcms/
detaiV23.1390. U.20131108.1314.006.html
中图分类号:
U662.9
文献标志码:
A
文章编号:
1006-
7043
( 2013) 11-1363-06
Application
of
preconditioned
and
restarted
GMRES to
the boundary element analysis
of
large floating structures
DUAN
Wenyang1,
DIAO
Feng2
,CHEN
Jikang1
(I.College
of
Shipbuilding Engineering, Harbin
Engineering
University,
Harbin
150001, China;
2.
China
Ship
Scientific
Research
Center,
Shanghai
200011, China)
Abstract:
Hydrodynamic boundary element analysis of large :floating bodies involves
linear
equations with complex
coefficients. Sometimes the direct method cannot solve it or will require a large amount of time. In
order
to solve it
efficiently, the restarted
and
preconditioned generalized minimal residual (
GMRES)
algorithm is presented. Two
different kinds of preconditioners are selected to preprocess the coefficient matrix of the
equations,
and
are com-
pared
numerically by detailed examples. Numerical experiments indicate
that,
for
linear
equations with complex co-
efficients generated from hydrodynamic
bounda
巧
r
element analysis of large :floating
bodies,
the
restarted GMRES al-
gorithm preconditioned with incomplete
LU
factorization is the optimal choice for
the
problems considered.
Keywords:
large :floating
structures;
boundary element
method;
restarted GMRES algorithm; incomplete LU factori-
zat10n
在船舶与海洋工程领域使用边界元法分析大
型浮体在规则波中的线性稳态响应是一种常用方
法,该方法使用奇点分布将描述水波问题的拉普拉
斯方程转化为边界积分方程.在大型浮体频域水动
力分析中,通过边界离散所形成的影响系数矩阵为
复系数满秩矩阵,由于面元尺寸的不一致和对物面
法向求导而导致格林函数的对称性的丢失,使得影
响系数矩阵非对称非正定.若采用直接法求解该类
收稿日期:
2012-11-29.
网络出版时间:
2013-11-8
13:
14.
基金项目:国家自然科学基金资助项目(
51079032)
.
作者简介:段文洋(
1967
-),男,教授,博士生导师
通信作者:段文洋,
E-mail:
duanwenyang@hrbeu.edu.cn.
线性方程组,所需要的运算量为
O
(
肘)
[!],这里
N
为问题的面元数.为了减少求解方程组所需的运算
量,更多的是采用迭代方法求解,其运算量约为
。
(
nN2)
,n
为迭代次数.
近几年,国内外许多学者在运用迭代法求解边
界元形成线性方程组方面进行了大量研究,广义极
小残值法(
generalized minimal
residual,
GMRES)
[2]
是一种适合于求解稠密非对称线性方程组的有效算
法.随着所研究问题越来越复杂以及规模的扩大,为
了加快迭代收敛速度并节约内存,往往需要采用预
条件处理及重启动措施重启动
GMRES
算法通常
第
34
卷第
11
期
2013
年
11
月
哈尔滨工程大学学报
Vol.34
No
.11
Nov.2013
Journal of
Harbin
Engineering University
预条件
GMRES(
m
)算法在大型浮体
水动力边界元分析中的应用
段文洋
1
,刁峰
2
,陈纪康
1
(
1.
哈尔滨工程大学船舶工程学院,黑龙江哈尔滨
150001
;2.
中国船舶科学研究中心上海分部,上海
200011)
摘
要:针对大型浮体水动力边界元分析产生的复系数线’性方程组结构复杂,直接方法难以求解或者求解费时的问题,
提出一种带有预条件技术的重启动型
GMRES
算法.选取
2
种不同的预条件处理技术对方程组系数矩阵进行预处理,通
过具体算例给出
2
种预条件方法的数值比较.数值试验表明对于大型浮体水动力边界元分析产生的复系数线性方程
组,带有不完全
LU
分解预条件处理技术的
GMRES(m
)算法求解效率最高,优于直接解法.
关键词:大型浮体;边界元分析;重启动
GMRES
;不完全
LU
分解预条件
doi: 10.3969/j.issn.1006-7043.201211097
网络出版地址:
http://wv
仰
.cnki.neι
/
kcms/
detaiV23.1390. U.20131108.1314.006.html
中图分类号:
U662.9
文献标志码:
A
文章编号:
1006-
7043
( 2013) 11-1363-06
Application
of
preconditioned
and
restarted
GMRES to
the boundary element analysis
of
large floating structures
DUAN
Wenyang1,
DIAO
Feng2
,CHEN
Jikang1
(I.College
of
Shipbuilding Engineering, Harbin
Engineering
University,
Harbin
150001, China;
2.
China
Ship
Scientific
Research
Center,
Shanghai
200011, China)
Abstract:
Hydrodynamic boundary element analysis of large :floating bodies involves
linear
equations with complex
coefficients. Sometimes the direct method cannot solve it or will require a large amount of time. In
order
to solve it
efficiently, the restarted
and
preconditioned generalized minimal residual (
GMRES)
algorithm is presented. Two
different kinds of preconditioners are selected to preprocess the coefficient matrix of the
equations,
and
are com-
pared
numerically by detailed examples. Numerical experiments indicate
that,
for
linear
equations with complex co-
efficients generated from hydrodynamic
bounda
巧
r
element analysis of large :floating
bodies,
the
restarted GMRES al-
gorithm preconditioned with incomplete
LU
factorization is the optimal choice for
the
problems considered.
Keywords:
large :floating
structures;
boundary element
method;
restarted GMRES algorithm; incomplete LU factori-
zat10n
在船舶与海洋工程领域使用边界元法分析大
型浮体在规则波中的线性稳态响应是一种常用方
法,该方法使用奇点分布将描述水波问题的拉普拉
斯方程转化为边界积分方程.在大型浮体频域水动
力分析中,通过边界离散所形成的影响系数矩阵为
复系数满秩矩阵,由于面元尺寸的不一致和对物面
法向求导而导致格林函数的对称性的丢失,使得影
响系数矩阵非对称非正定.若采用直接法求解该类
收稿日期:
2012-11-29.
网络出版时间:
2013-11-8
13:
14.
基金项目:国家自然科学基金资助项目(
51079032)
.
作者简介:段文洋(
1967
-),男,教授,博士生导师
通信作者:段文洋,
E-mail:
duanwenyang@hrbeu.edu.cn.
线性方程组,所需要的运算量为
O
(
肘)
[!],这里
N
为问题的面元数.为了减少求解方程组所需的运算
量,更多的是采用迭代方法求解,其运算量约为
。
(
nN2)
,n
为迭代次数.
近几年,国内外许多学者在运用迭代法求解边
界元形成线性方程组方面进行了大量研究,广义极
小残值法(
generalized minimal
residual,
GMRES)
[2]
是一种适合于求解稠密非对称线性方程组的有效算
法.随着所研究问题越来越复杂以及规模的扩大,为
了加快迭代收敛速度并节约内存,往往需要采用预
条件处理及重启动措施重启动
GMRES
算法通常