Laguerre函数在时域匹配边界元中的应用研究

李志富，任慧龙，石玉云，李辉

(哈尔滨工程大学 船舶工程学院，黑龙江 哈尔滨 150001)

Laguerre函数在时域匹配边界元中的应用研究

李志富，任慧龙，石玉云，李辉

(哈尔滨工程大学 船舶工程学院，黑龙江 哈尔滨 150001)

摘要：利用满足自由表面条件的时域格林函数求解波浪与浮体相互作用问题时，对于外飘型浮体计算结果时常发散。针对此问题，建立了内外场时域快速匹配边界元模型，通过与解析解对比，验证了数值模型的有效性。其次，利用所建立的匹配边界元模型，研究了流场记忆速度势随时间的变化特点。针对长时间模拟，时间卷积积分项导致模拟失败问题，提出了利用Laguerre函数对记忆速度势时间变量进行级数逼近，以对时间卷积积分进行解析计算，并对级数收敛速度进行了研究。最后，为加快级数收敛速度，引入了放缩系数，并对其特性进行了探讨。数值结果表明，通过引入放缩系数，可以利用较少的级数项来逼近速度势的时间变量，能够很好地描述记忆速度势的脉冲特性。

关键词：Laguerre函数；时域匹配；时间卷积；多项式逼近；放缩系数

网络出版地址：http://www.cnki.net/kcms/detail/23.1390.U.20160411.0845.014.html

波浪与浮体相互作用的时域模拟具有重要意义[1]。在势流理论框架下，根据所选取格林函数的不同，主要分为时域Rankine源法和时域自由面格林函数法。对于时域自由面格林函数法，只需要在船体表面进行积分，国内外许多学者都在此方面做了研究[2-3]。此方法对于不同浮体进行模拟时，都需要重新计算时域格林函数在每一时刻的值，而时域格林函数的计算又需要占用较长的CPU时间，所以该方法计算效率较低。此外，格林函数在近水面单元的高频振荡特性[4]，导致外飘型船舶的计算，结果时常发散。对于时域Rankine源法，其形式简单，且易与高阶边界元法相结合，以进行非线性问题的模拟，所以现在仍然有很多学者对此进行研究[5-6]。但是，时域Rankine法需要对计算域进行截断，且需要布置相应的人工远方辐射条件。

通过对时域Rankine源法和时域自由面格林函数法各自优缺点的分析，有学者提出了匹配方法[7-8]，此种方法不仅避免了不同浮体都需对格林函数进行计算的问题，同时适用于外飘型船舶的计算。但是此方法的准确性和有效性仍需进一步验证。无论是时域Rankine源法还是匹配方法，在计算与时间相关的记忆速度势时，时间卷积项的存在，常常导致长时间的模拟无法进行。

针对上述问题，本文首先建立了三维内外场时域快速匹配边界元模型，通过与解析解的对比，验证了数值模型的有效性。其次利用该模型研究了流场中不同点的记忆速度势随时间的变化特点。最后根据记忆速度势随时间的变化特性，提出了利用Laguerre函数[9]对记忆速度势时间变量进行级数逼近，并通过引入一放缩系数加快了级数的收敛速度。

1基础理论

1.1内外场匹配模型

根据文献[10]，假定流场为不可压缩的理想流体，则流场中的物理量可以用流体速度势Φ来进行描述。为建立三维时域匹配模型，引入一虚拟控制面，如图1所示。

图1　内外域匹配求解示意图Fig.1　Coordinate system and boundary surface

流体速度势在内域满足以Rankine源所建立的边界积分方程：

(1)

流体速度势在外域满足以时域自由面格林函数所建立的边界积分方程：

(2)

(3)

(4)

同理，速度势ΦII在外域满足的边界积分方程(2)也可以离散为相应的矩阵形式

(5)

式中：rC表示外域所满足的边界积分方程的时间卷积项。

在虚拟控制面SC上，流体速度势及其法向导数应满足如下连续条件：

(6)

通过速度势在内域满足的边界积分方程式(4)和外域满足的边界积分方程式(5)，便可以得到内外域匹配求解方程组：

(7)

式中矩阵的系数如下

(8)

(9)

(10)

(11)

由式(7)，在求解流场速度势时，只有虚拟控制面SC上涉及了时域自由面格林函数波动项的计算。

1.2流场速度势分解

由于流体速度势所满足的场方程、边界条件和初始条件均是线性化的，故可以对流体速度势进行线性化分解[11]

(12)

式中：ΦI、ΦD和Φk分别为入射波速度势、绕射波速度势和辐射波速度势，其中ΦI为已知量，ΦD和Φk为待求解的未知量。

对于辐射波速度势Φk，根据其所代表的物理意义，可以将其分解为与时间无关的瞬时速度势Ψk和与时间相关的记忆速度势χk:

(13)

根据瞬时速度势和记忆速度势在内外域满足的边界条件和初始条件[3]，利用式(7)便可求出相应的值。根据线性化的Bernoulli方程[10]便可求出对应的辐射波力

(14)

式中：μjk和Kjk分别为与频率无关的附加质量和时域延迟函数。其可以相应的转化为频域附加质量Ajk和阻尼系数Bjk:

(15)

(16)

1.3时间卷积积分项

(17)

(18)

根据记忆速度势随时间的变化特性，为加快级数的收敛速度，引入一放缩系数a:

(19)

(20)

(21)

2数值结果

2.1匹配边界元有效性验证

利用时域快速匹配边界元法求解浮体和波浪相互作用问题时，影响计算精度的主要因素有时间步长、边界网格数目、自由面网格形式、虚拟控制面和物体之间的距离。

为讨论自由面网格形式的影响，针对相同工况，分别选择了矩形网格、辐射状网格和矩形三角形混合网格进行了计算，通过与解析解[12]的对比验证发现辐射状的自由面网格能够最好的描述辐射波和绕射波的传播，如图2所示，图中半球半径为1 m，控制面半球半径为5 m。此与辐射和绕射波代表外传柱面波的物理意义相一致。

图2　匹配边界元网格划分示意图Fig.2　Grid pattern for Rankine-Green method

为讨论时间步长对计算结果的影响，针对相同工况，分别选取了△t=0.02、0.01、0.005 s进行了计算。计算得到的无量纲时域延迟函数如图3所示，3条曲线几乎重合。此外，针对边界网格数目和虚拟控制面与物体之间的距离也做了同样的对比分析。最后，考虑数值计算精度和计算效率，在文章后续计算中，自由面划分为辐射式网格，计算时间步长选为0.01 s，虚拟控制面和物面间距为5倍的物体特征长度，网格数目为物面130、自由面200、控制面150。

图3　时间步长对计算结果的影响Fig.3　Influence of time step on the numerical results interval

对于半球在静水中做强迫垂荡运动，分别计算了其附加质量系数和阻尼系数，如图4、5所示。由图可知，利用本文建立的快速时域匹配边界元法计算的结果与文献[12]提供的解析解符合良好，说明本文提供的数值模型能够准确模拟波浪和浮体相互作用的流场。

图4　无因次附加质量系数Fig.4　Non-dimensional added mass coefficient

图5　无因次兴波阻尼系数Fig.5　Non-dimensional damping coefficient

2.2Laguerre函数有效性验证

为验证Laguerre函数的有效性，首先对记忆速度势随时间的变化特性进行了研究。记忆速度势代表浮体干扰形成表面波浪，开始波浪传播之后的那部分速度势，并且随时间逐渐衰减，如图6所示，图中点1位于虚拟控制面和自由表面相交处(5,0,0)，点2位于虚拟控制面的中部(4.33,0.00,2.50)，点3位于虚拟控制面的底部(0,0,5)。图6表明随着场点远离自由面，记忆速度势的变化逐渐趋于平缓，此与记忆速度势代表物体兴波部分速度势的物理意义相一致；且记忆速度势随时间出现一个较大的脉冲峰以后，便开始迅速衰减，此与修正Laguerre函数中的指数函数部分相符合，初步说明了利用Laguerre函数进行记忆速度势时间变量级数逼近的可行性。

图6　记忆速度势随时间变化Fig.6　 Memory potential for different points

根据图6中的3组数据，本文对随时间变化最为剧烈的点1进行了Laguerre函数逼近研究。在不引入放缩系数的情况下，函数逼近随级数项数的变化如图7所示，绝对误差随项数的变化如图8所示。由图可知，随着逼近级数项数的增加，级数解越趋近于原数据，但是收敛效果随项数的增加改善较为缓慢，例如当级数项数增加一倍时，逼近绝对误差缩小不到一倍，而计算求解矩阵的维数却大大增加。

图7　Laguerre级数逼近精度与截断项数的关系(a=1)Fig.7　Relationship between series approximation and 　　truncated series terms (a=1)

图8　Laguerre级数逼近绝对误差与截断项数的关系(a=1)Fig.8　Relationship between absolute error and truncated 　　series terms (a=1)

图9　Laguerre级数逼近与放缩系数的关系(N=15)Fig.9　Relationship between series approximation and　　 flexible parameter (N=15)

图10　Laguerre级数逼近绝对误差与放缩系数的关系(N=15)Fig.10　Relationship between absolute error and flexible 　　parameter (N=15)

为了增加Laguerre级数逼近的收敛速度，同时又不增加计算机求解负担，引入了放缩系数的概念，如式(21)所示。由图10可知，随着放缩系数的减小，相对误差逐渐减小，当放缩系数等于0.2时，级数项取为15的结果要明显优于不取放缩系数的级数项取为20的结果，如图11所示。由此可以确定，通过引入放缩系数，可以利用较少的项数来逼近函数的真解，从而确保了利用Laguerre函数进行记忆速度势时间变量进行级数逼近的经济性和可行性。

图11　放缩系数对级数收敛速度影响Fig.11　 Influence of flexible parameter and series terms

4结论

本文通过引入虚拟控制面，将流场划分为內域和外域，并根据速度势及其法向导数在控制面上连续的条件建立了时域快速匹配边界元模型。此外，提出了利用Laguerre函数对记忆速度势的时间轴进行级数逼近，以解析计算时间卷积积分的思路，并通过引入一放缩系数，大大加快了级数的收敛速度。通过对计算结果的分析，得到以下结论：

1)通过与解析解的对比表明本文建立的内外场时域快速匹配边界元法具有较高的数值计算精度，且由求解矩阵可知，时域快速匹配边界元法与直接边界元法相比具有较高的计算效率，为研究浮体和波浪相互作用下流场的特性提供了一种精确、高效的方法；

2)自由面划分为辐射状网格，控制面距离物面为5倍的物体的特征尺度，时间步长取为0.01 s时，时域快速匹配边界元法便能给出令人满意的计算结果；

3)代表浮体兴波之后流场特性的记忆速度势在出现一个很大的峰之后，随时间迅速衰减，当t=2.5 s后，基本趋于平缓；

4)利用Laguerre函数对记忆速度势的时间变量做级数逼近时，随着级数项数的增加，级数逐渐收敛，但级数收敛速度较为缓慢；

5)分析Laguerre函数及其正交权函数的性质，引入一放缩系数，数值结果显示，放缩系数可大大加快级数的收敛速度，N=15，a=0.2的逼近结果要明

显的优于N=20，a=1时的结果,保证了时域快速解析匹配边界元的可行性与高效性。

参考文献：

[1]SEO M G, KIM Y. Numerical analysis on ship maneuvering coupled with ship motion in waves[J]. Ocean engineering, 2011, 38(17/18): 1934-1945.

[2]BECK R, LIAPIS S. Transient motions of floating bodies at zero forward speed[J]. Journal of ship research, 1987, 31(3):164-176.

[3]KING B. Time-domain analysis of wave exciting forces on ships and bodies[R]. Michigan: University of Michigan, 1987.

[4]CLÉMENT A H. An ordinary differential equation for the Green function of time-domain free-surface hydrodynamics[J]. Journal of engineering mathematics, 1998, 33(2): 201-217.

[5]SONG M J, KIM K H, KIM Y. Numerical analysis and validation of weakly nonlinear ship motions and structural loads on a modern containership[J]. Ocean engineering, 2011, 38(1): 77-87.

[6]SEO M G, YANG K K, PARK D M, et al. Numerical analysis of added resistance on ships in short waves[J]. Ocean engineering, 2014, 87: 97-110.

[7]LIU S, PAPANIKOLAOU A. Application of a nonlinear time domain hybrid method to the study of a semi-submersible in waves[C]//Proceedings of the 22nd international offshore (ocean) and polar engineering conference,Rhodes, Greece, 2012.

[8]韩旭亮, 段文洋. 时域匹配直接边界元方法及其数值特性[J]. 哈尔滨工程大学学报, 2013(7): 837-843.

HAN Xuliang, DUAN Wenyang. Study on the numerical characteristics of time domain matching direct boundary element method[J]. Journal of Harbin engineering university, 2013(7): 837-843.

[9]ABRAMOWITZ M, STEGUN I A. Handbook of mathematical functions: With formulas, graphs, and mathematical tables[M]. (s.l.):Courier Dover Publications, 1972.

[10]WEHAUSEN J V, LAITONE E V. Surface waves[M]. Springer, 1960.

[11]TANG K, ZHU R, MIAO G, et al. Domain decomposition and matching for time-domain analysis of motions of ships advancing in head sea[J]. China ocean engineering, 2014, 28: 433-444.

[12]HULME A. The wave forces acting on a floating hemisphere undergoing forced periodic oscillations[J]. Journal of fluid mechanics, 1982, 121: 443-463.

本文引用格式：

李志富，任慧龙，石玉云，等. Laguerre函数在时域匹配边界元中的应用研究[J]. 哈尔滨工程大学学报， 2016, 37(5): 629-633.

LI Zhifu，REN Huilong，SHI Yuyun， et al. Application of Laguerre function in the time domain Rankine-Green panel method[J]. Journal of Harbin Engineering University, 2016, 37(5): 629-633.

Application of Laguerre function in the time domain Rankine-Green panel method

LI Zhifu，REN Huilong，SHI Yuyun，LI Hui

(College of Shipbuilding Engineering, Harbin Engineering University, Harbin 150001, China)

Abstract：The divergence phenomenon will happen when the transient free surface Green function method is adopted to simulate the motion response of the ships with large flare. To avoid that problem, the Rankine-Green panel method is proposed , and the accuracy of this method is verified by the comparision with the analytical solutions. Then the hybrid method is used to investigate the variation properties of the memory velocity potential wiht respect to time, and the Laguerre function is introduced to express the memory potentials in the Laguerre series form. Therefore, the temporal convolution integral of the memory potentials with the Green function can be computed analytically, which ensures the long time simulations. Besiedes, to accelerate the convergence speed of the Laguerre series, a flexible parameter is introduced and its features are discussed. The numerical results show that the temporal variable of the memory potential can be approximated by less Laguerre series terms with a flexible parameter, which means that the Laguerre series can give a reasonable description of the impulse properties of the memory velocity potential.

Keywords：Laguerre function; Rankine-Green panel method; temporal convolution integral; series approximation; flexible parameter

收稿日期：2015-04-10.

基金项目：国家“973”计划项目(2011CB3703)；国家自然科学基金资助项目(51109046).

作者简介：李志富(1990-), 男, 博士研究生； 通信作者：任慧龙, E-mail：renhuilong@263.net.

DOI:10.11990/jheu.201504023

中图分类号：O352

文献标志码：A

文章编号：1006-7043(2016)05-0629-05

网络出版时间：2016-04-11.

任慧龙(1965-), 男, 教授, 博士生导师.