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.     

A modal approach to second-order analysis of sloshing using boundary element method

This paper aims at developing a modal approach for the non-linear analysis of sloshing in an arbitrary-shape tank under both horizontal and vertical excitations. For this purpose, the perturbation technique is employed and the potential flow is adopted as the liquid sloshing model. The first- and second-order kinematic and dynamic boundary conditions of the liquid-free surface are used along with a boundary element model which is formulated in terms of the velocity potential of the liquid-free surface. The boundary element model is used to determine the natural mode shapes of sloshing and their corresponding frequencies. Using the modal analysis technique, a non-linear model is presented for the calculation of the first- and second-order potential which can be used to obtain a reduced-order model for the sloshing dynamics. The results of the presented model are verified with the analytical solution for the second-order analysis of sloshing in a rectangular tank and very good results were obtained. Also, the second-order sloshing in some other example tanks with complex bed shapes is studied. The second-order resonance conditions of liquid sloshing in the example tanks are investigated and some conclusions are drawn.     

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.     

预修正快速傅里叶变换高阶边界元方法在波浪对物体作用问题中的应用

边界元方法被广泛应用于波浪对海上婕筑物作用领域,但由于传统边界元方法的存储量和计算量均为未知量的平方量级,很难满足大范围多未知量问题的计算需要.本文基于高阶边界元方法,应用预修正快速傅里叶变换方法,使计算量与存储量分别降低至O(NlogN)和O(N)量级,并得到一个连续的压强分布以适应结构设计的要求,同时可以通过使用满...     

基于高阶边界元的三维数值波浪港池--波浪破碎的模拟

在势流理论的框架内,采用高阶边界元方法和混合欧拉-拉格朗日法,实现了对三维波浪破碎过程的数值模拟.数值模型使用可调节时间步长的基于二阶显式泰勒展开的混合欧拉-拉格郎日时间步进来求解自由表面的演化过程.在所使用的边界元方法中,采用16节点三次滑移四边形单元来表示,这种单元在单元内具有高阶的精度同时在单元之间具有良好的连续性.给出了孤立波的传播和周期性非线性波浪沿缓坡传播的计算结果,表明数值模型具有良好的稳定性.     

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

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

A novel decomposition of the quadratic transfer function (QTF) for the time-domain simulation of non-linear wave forces on floating bodies

In the present study, a novel method is proposed for the separation of the second-order sum- and difference-frequency wave forces—that is, quadratic transfer functions (QTFs)—on a floating body into three components due to wave–wave, wave–motion, and motion–motion action. By applying the new QTF components, the second-order wave forces on a floating body can be strictly computed in the time domain. In this work, the boundary value problems (BVPs) corresponding to the three kinds of QTF components were derived, and non-homogeneous boundary conditions on the free surface and the body surface were obtained. The second-order diffraction potentials were determined using the boundary integral equation method. In the solution procedure, the highly oscillatory and slowly converging integral on the free surface was evaluated in an accurate and effective manner. Furthermore, the application of the QTF components in the time domain was demonstrated. The second-order exciting forces in the time domain were divided into three parts. Each part of these forces was computed via a two-term Volterra series model based on the incident waves, the first-order motion response, and the QTF components. This method was applied to several numerical examples. The results demonstrated that this decomposition yields satisfactory results.     

Dynamic modeling of cable towed body using nodal position finite element method

This paper analyses nonlinear dynamics of cable towed body system. The cable has been modeled and analyzed using a new nodal position finite element method, which calculates the position of the cable directly instead of the displacement by the existing finite element method. The newly derived nodal position finite element method eliminates the need of decoupling the rigid body motion from the total motion, where numerical errors arise in the existing nonlinear finite element method, and the limitation of small rotation in each time step in the existing nonlinear finite element method. The towed body is modeled as a rigid body with six degrees of freedom while the tow ship motion is treated as a moving boundary to the system. A special procedure has been developed to couple the cable element with the towed body. The current approach can be used as design tool for achieving improved directional stability, maneuverability, safety and control characteristics with the cable towed body. The analysis results show the elegance and robustness of the proposed approach by comparing with the sea trial data.     

应用比例边界有限元法求解多种二维浮体水动力特性

比例边界有限元法(SBFEM)是一种半解析的数值方法,比完全数值方法具有更高的精度,该方法结合了有限元和边界元的优点,采用相对少的剖分单元就可以得到较高精度的模拟结果。通过改变有限子域内部比例中心的位置,使这种方法可以应用到多种形式浮体在波浪作用下的水动力特性的计算中。同时还给出了各种形式浮体的波浪力及反射、透射系数的数值结果,并与边界元方法(BEM)计算结果和特征函数展开方法得到解析解进行了比较,均吻合良好。研究表明比例边界有限元不仅可以计算矩形的浮体结构,而且对于多种结构形式的浮体都可以计算,这为多种结构形式浮体的水动力分析提供了一个可行的方法。     

应用贝叶斯方法对有效信息较少的渔业资源进行资源评估-以印度洋长鳍金枪鱼的资源评估为例

众所周知,对有效信息较少的渔业资源进行资源评估面临很大的挑战,而贝叶斯方法在数据数量较少、质量较差的情况下能利用其它种群高质量的数据或已知的先验信息提高资源评估结果的可靠性。由于印度洋长鳍金枪鱼的数据质量较差而数据量有限,长鳍金枪鱼的资源评估结果存在很大的不确定性,为此,本文以印度洋长鳍金枪鱼的资源评估为例,以调查贝叶斯方法在有效信息较少的资源评估中的优势。本文根据不同的先验假设与捕捞数据系列,共构建了8个贝叶斯动态产量模型,以评估长鳍金枪鱼资源。结果表明:(1)分析参数的后验分布能提高捕捞数据系列选择与参数假设的合理性; (2) 利用种群统计学方法为内禀增长率(r)构建有信息先验,能提高资源评估结果的可靠性。与传统方法相比,当基于贝叶斯框架时,能将已知的知识表示为先验信息并能分析参数的后验分布,从而在数据较少或数据质量较差的情况下,能利用各种信息提高参数估计的合理性与资源评估的可靠性。因此,对数据量较少或数据质量较差情况下的渔业资源评估而言,贝叶斯方法非常有效,如本文所示的印度洋长鳍金枪鱼的资源评估。