预修正快速傅里叶变换方法在波浪绕射和辐射问题中应用

在波浪对海上建筑物绕射和辐射问题的计算中,面元法被广泛使用,但由于传统面元法的存储量和计算量均为未知量的平方量级,很难满足大范围多未知量问题的计算需要。采用预修正快速傅里叶变换方法(pFFT方法),使计算量与存储量都降低到未知量的线性量级。以淹没圆球与漂浮圆柱两个典型算例为基础,通过不同未知量时pFFT方法与传统面元法的计算量与存储量的对比,以及pFFT方法自身各步骤计算时间的对比,研究了不同网格方案的选取对pFFT方法计算量和存储量的影响,推荐根据未知量个数采用计算时间最小化原则选取pFFT网格参数。     

无限水深下波浪与二维水面物体作用的简单格林函数方法

针对无限水深下波浪与二维水面物体相互作用问题,传统的波浪格林函数形式复杂、计算缓慢,为了提高计算效率和计算精度,将流域分为物体周围的内域及远离物体的外域,内域采用简单格林函数法,外域采用多极子展开方法,通过内外域边界匹配,耦合求解得到流域中任意一点的速度势,并可计算物体在波浪作用下的波浪激振力、附加质量、辐射阻尼及透射和反射系数。应用该方法计算了二维水面漂浮半圆和水面漂浮方箱的算例,数值计算结果表明,该方法可以方便、准确、快速地计算无限水深下波浪与任意漂浮物体的作用问题。     

高阶边界元方法求解规则波在三维局部渗透海床上的传播

基于线性势流理论,利用高阶边界元法研究了规则波在三维局部渗透海床上的传播。根据Darcy渗透定律推导出渗透海床的控制方程,利用渗透海床顶部和海底处法向速度和压强连续条件得到渗透海床顶部满足的边界条件。根据绕射理论,利用满足自由水面条件的格林函数建立了求解渗透海床绕射势的边界积分方程,采用高阶边界元方法求解边界积分方程进而得到自由水面的绕射势和波浪在局部渗透海床上传播过程中幅值的变化情况。通过与已发表的波浪对圆柱形暗礁的时域全绕射结果对比,证明了本文建立的频域方法计算波幅的正确性和有效性。利用这一模型研究了三维矩形渗透海床区域上波浪的传播特性,并分析了入射波波长、海床渗透特性系数等参数对波浪传播的影响。     

An Efficient Model for Transient Surface Waves in Both Finite and Infinite Water Depths

A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropr...     

An analysis of second-order wave forces on floating bodies by using a higher-order boundary element method

The wave diffraction and radiation problem is studied numerically by using a higher-order boundary element method. The convergence of the higher-order boundary element method is tested systematically for bodies of different shapes. For the second-order force, particular attention is paid to the contribution of the second-order potential, following the line of Molin's approach. For numerical evaluation, the free surface is divided into three subregions; inner, intermediate and outer ones. In the inner region, the integral is evaluated numerically by using higher-order boundary elements. In the intermediate region, semi-analytic form is constructed with the help of eigen functions. In the outer domain, the analytic solution is available. This subdivision scheme reduces the numerical burden remarkably.     

An Efficient Model for Transient Surface Waves in Both Finite and Infinite Water Depth

A numerical model is developed to simulate fully nonlinear extreme waves in finite and infinite water-depth wave tanks. A semi-mixed Eulerian-Lagrangian formulation is adopted and a higher-order boundary element method in conjunction with an image Green function is used for the fluid domain. The boundary values on the free surface are updated at each time step by a fourth-order Runga-Kutta time-marching scheme at each time step. Input wave characteristics are specified at the upstream boundary by an appropriate wave theory. At the downstream boundary, an artificial damping zone is used to prevent wave reflection back into the computational domain. Using the image Green function in the whole fluid domain, the integrations on the two lateral walls and bottom are excluded. The simulation results on extreme wave elevations in finite and infinite water-depths are compared with experimental results and second-order analytical solutions respectively. The wave kinematics is also discussed in the present study.     

A numerical study of second-order springing of an elastic body using higher-order boundary element method (HOBEM)

Second-order springing on an elastic body with forward speed is analyzed by numerical simulations. The boundary-value problem for the velocity potential is solved by means of the direct time-domain higher-order boundary element method (HOBEM). The free-surface boundary condition in the boundary-value problem is approximated on the mean surface up to second order by use of perturbation and Taylor-series expansion methods. The body boundary condition for an elastic body is derived with various quantities which are redefined in the generalized mode. These variables such as mode shape, normal vector, etc. are obtained by using directional derivative and continuum mechanics, and the same mathematical expressions are used to obtain several second-order generalized forces. To validate the numerical results, the second-order hydrodynamic force on the bottom-mounted rigid/elastic cylinders without forward speed is compared with other semi-analytic results. The property of second-order forces on an elastic ship is studied by changing the flexural rigidity and forward speed with elastic response. It is confirmed that the second-order velocity potential is important for a body with forward speed and investigation should be made more on numerical methods for accurate computation of the second-order velocity-potential force with forward speed.     

面元法的快速算法与B样条高阶方法综述

介绍了面元法的快速算法及B样条高阶方法的发展概况，对快速多极算法及预修正快速傅里叶变换方法做了详细的介绍，指出了这些方法今后的发展方向。     

Extension of the Frequency-Domain pFFT Method for Wave Structure Interaction in Finite Depth

To analyze wave interaction with a large scale body in the frequency domain, a precorrected Fast Fourier Transform (pFFT) method has been proposed for infinite depth problems with the deep water Green function, as it can form a matrix with Toeplitz and Hankel properties. In this paper, a method is proposed to decompose the finite depth Green function into two terms, which can form matrices with the Toeplitz and a Hankel properties respectively. Then, a pFFT method for finite depth problems is developed. Based on the pFFT method, a numerical code pFFT-HOBEM is developed with the discretization of high order elements. The model is validated, and examinations on the computing efficiency and memory requirement of the new method have also been carried out. It shows that the new method has the same advantages as that for infinite depth.     

基于机器学习的脉动源格林函数预报初探

自由面格林函数是边界元法求解海洋工程水动力学问题的基础,如何精确而快速地计算格林函数及其导数是水动力求解的难题。对无因次化表达的脉动点源格林函数计算建立的数据库,采用深度学习函数库Keras,对数据库进行学习,建立了神经网络预报模型,探讨了全局和局部学习及预报精度,研究了模型预报效率。研究表明机器学习模型预报的格林函数能够保证较高的精度,计算效率高于数值积分计算,低于解析函数为主的多项式逼近,为提高水动力问题求解效率,解决传统计算难题提供了新的思路。     

Second-order Taylor expansion boundary element method for the second-order wave radiation problem

A novel Boundary Element Method (BEM) named the second-order Taylor Expansion Boundary Element Method (the 2nd order TEBEM) is developed for the solution of the second-order wave radiation velocity potential and sum-frequency wave loads for floating bodies. The radiation condition is enforced by a hybrid method of the multi-transmitting formula and damping zone. For the interior domain problem of a cube and a sphere, numerical results demonstrate that the 2nd order TEBEM can accurately solve the first and second-order gradients of velocity potential on the no-smoothed and smoothed boundary compared to the low-order BEM. The double frequency forces acting on the truncated cylinder are calculated under finite water depth. The agreement between the 2nd order TEBEM and others' numerical results is good. Moreover, all of the singular integrals in the 2nd order TEBEM can be solved analytically, so its implementation is much easier compared to the high-order BEM.     

消除"不规则频率"的非连续高阶元方法

针对使用边界元法计算波浪与结构物相互作用时所出现的“不规则频率”现象，采用连续高阶元和部分非连续高阶元对通过修改积分区域所获得的边界积分方程进行离散，有效地消除了“不规则频率”现象的发生。波浪作用下的截断圆柱所受到的水平波浪力和垂向波浪力的数值计算结果验证了该方法的有效性，同时考虑了非连续单元配置点的选择及单元划分数目对消除效果的影响。