A Unique Solvable Higher Order BEM for Wave Diffraction and Radiation

- For the discretization of higher order elements, this paper presents a modified integral domain method to remove the irregular frequencies inherited in the integral equation of wave diffraction and radiation from a surface-piercing body. The set of over-determined linear equations obtained from the method is modified into a normal set of linear equations by superposing a set of linear equations with zero solutions. Numerical experiments have also been carried out to find the optimum choice of the size of the auxiliary domain and the discretization on it.     

波浪在三维圆形岛地形上的绕射研究

波浪在传播过程中遇到岛屿就会发生绕射。本文使用混合元方法对修正型缓坡方程进行了数值求解,并与KUO et al的解析解进行了比较验证。在此基础上研究了工程尺度背景下,波浪在三维圆形岛地形上的绕射,计算了不同入射波浪周期、浅滩形状参数和岛屿尺寸情况下,沿波浪传播方向断面上和岛屿岸线上的相对波高大小。计算结果表明：随着入射波周期的减小、浅滩形状参数的增大和岛屿尺寸的减小,圆形岛迎浪侧的相对波高振荡幅度、圆形岛背浪侧的相对波高大小以及岛屿岸线上的相对波高振幅和大小均随之增大。不同情况下,岛屿岸线上的相对波高最大值大多数发生在迎浪点,个别发生在迎浪点两侧20°~25°处;最小值发生在背浪点两侧30°附近。     

A composite numerical model for wave diffraction in a harbor with varying water depth

A composite numerical model is presented for computing the wave field in a harbor. The mild slope equation is discretized by a finite element method in the domain concerned. Out of the computational domain, the water depth is assumed to be constant. The boundary element method is applied to the outer boundary for dealing with the infinite boundary condition. Because the model satisfies strictly the infinite boundary condition, more accurate results can be obtained. The model is firstly applied to compute the wave diffraction in a narrow rectangular bay and the wave diffraction from a porous cylinder. The numerical results are compared with the analytical solution, experimental data and other numerical results. Good agreements are obtained. Then the model is applied to computing the wave diffraction in a square harbor with varying water depth. The effects of the water depth in the harbor and the incoming wave direction on the wave height distribution are discussed.     

A composite numerical model for wave diffraction in a harbor with varying water depth

A composite numerical model is presented for computing the wave field in a harbor. The mild slope equation is discretized by a finite element method in the domain concerned. Out of the computational domain, the water depth is assumed to be constant. The boundary element method is applied to the outer boundary for dealing with the infinite boundary condition. Because the model satisfies strictly the infinite boundary condition, more accurate results can be obtained. The model is firstly applied to compute the wave diffraction in a narrow rectangular bay and the wave diffraction from a porous cylinder. The numerical results are compared with the analytical solution, experimental data and other numerical results. Good agreements are obtained. Then the model is applied to computing the wave diffraction in a square harbor with varying water depth. The effects of the water depth in the harbor and the incoming wave direction on the wave height distribution are discussed.     

A 3-D Numerical Model for the Calculation of Water Wave Transformation in Large Area

Based on the integral equation transformed from three dimensional Laplace equation and by the adoption of the division manner of sub- region boundary element method, the numerical computations of the velocity potential of each sub-region are given considering the continuity conditions of potential and normal derivatives at the interface of sub-regions. Therefore, computation of wave deformation in offshore flow field is realized. The present numerical model provides a good solution for the application of boundary element method to the calculation of wave deformation in large areas.     

波浪与结构物作用分析的一种高阶边界元方法--自由项和柯西主值积分的直接数值计算

采用直接数值计算方法计算了势流问题高阶边界元方法中的自由项系数和柯西主值积分,建立了波浪与结构物作用的一种高阶边界元方法.通过算例研究了物体表面上固角系数的计算精度和不同网格剖分、不同阶高斯积分点对柯西主值积分的影响.对截断圆柱上的波浪作用力与解析解做了对比,发现本方法具有很高的计算精度,随网格的加密迅速收敛于解析解.     

An extension of general identities for 3D water-wave diffraction with application to the Diffraction Transfer Matrix

Interaction theories are used in numerous branches of physics to efficiently evaluate wave scattering by multiple obstacles. An example of these interaction theories is the direct matrix method introduced by Kagemoto and Yue [1], which enables fast computation of three-dimensional water-wave multiple-scattering problems. The building block of interaction theories is a mathematical operator that encapsulates the mapping between incident and scattered waves. This operator is generally referred to as T-matrix and satisfies both reciprocity and energy identities. In some branches of physics, such as acoustics and electromagnetism, these identities are well established; in hydrodynamics, however, they have only been derived for a T-matrix that maps two-dimensional incident and scattered water waves. In three dimensions, water waves can be represented as a series expansion of cylindrical eigenfunctions. In this paper, we use this representation of water waves to derive the reciprocity and energy identities satisfied by the T-matrix of the direct matrix method, known as Diffraction Transfer Matrix (dtm). The identities derived herein represent an extension of existing general relations between two diffraction solutions. We show that this extension can be applied to verify the accuracy of the dtm entries, thereby increasing the reliability of existing schemes for computing the dtm. We present results for the dtm of two geometrically different isolated obstacles, as well as for the dtm of an asymmetric array. Finally, we demonstrate that the results presented herein can be extended to floating bodies found in a wide range of ocean engineering problems.     

新型三维内波及分层流水槽系统

本文提供了作者研制的三维内波及分层流水槽系统的背景技术资料。经源致三维内波的实验证实 ,利用该水槽系统中的双缸制取分层水系统可快速制取所需的分层水 ,不必利用分子扩散而 1次得到预先设置的密度剖面。同时 ,还介绍了利用该水槽系统进行的 1个典型的源致内波的内波增阻实验。从所获得的增阻实验结果得知 ,对于一给定的分层、拖曳速度和水下回转体 ,存在唯一 1个内Froude数使得强迫源受到最大的内波增阻。大量的实验证实 :该水槽系统是研究三维自由和强迫内波的 1个基本实验系统。     

A finite difference method for effective treatment of mild-slope wave equation subject to non-reflecting boundary conditions

A numerical model for coastal water wave motion that includes an effective method for treatment of non-reflecting boundaries is presented. The second-order one-way wave equation to approximate the non-reflecting boundary condition is found to be excellent and it ensures a very low level of reflection for waves approaching the boundary at a fairly wide range of the incidence angle. If the Newman approximation is adopted, the resulting boundary condition has a unique property to allow the free propagation of wave components along the boundary. The study is also based on a newly derived mild-slope wave equation system that can be easily made compatible to the one-way wave equation. The equation system is theoretically more accurate than the previous equations in terms of the mild-slope assumption. The finite difference method defined on a staggered grid is employed to solve the basic equations and to implement the non-reflecting boundary condition. For verification, the numerical model is then applied to three coastal water wave problems including the classical problem of plane wave diffraction by a vertical circular cylinder, the problem of combined wave diffraction and refraction over a submerged hump in the open sea, and the wave deformation around a detached breakwater. In all cases, the numerical results are demonstrated to agree very well with the relevant analytical solutions or with experimental data. It is thus concluded that the numerical model proposed in this study is effective and advantageous.     

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.     

Evaluation of the Risk of Marine Slope Instability: A Pseudo-3D Approach for Application to Large Areas

Abstract

This article presents a methodology developed to evaluate the instability of submarine slopes that extend over a large area. Special attention was paid to (1) the complex geometry (bathymetry) and the expanse of the slope, (2) the heterogeneity of the sediment, and (3) the distribution of the pore pressure. The safety factor was considered as a spatially varying quantity. The General Formulation (GLE, Fredlund and Krahn 1977), which fully satisfies equilibrium conditions, was used for evaluating the stability of the marine slope. The submarine slope failure, which occurred on 16 October 1979 during the construction of the new Nice airport, was studied in order to test the developed model. Geotechnical parameters were taken from experimental tests carried out by IFREMER on 19 cores extracted at different depths (from 27 m to 1300 m) (Cochonat, Bourillet, and Savoye, 1993; Mulder et al., 1994). Many scenarios were proposed in order to explain the cause of the Nice slope failure (Habib, 1994). In this article, two of those scenarios were tested. Simulation results are presented and discussed.     

Evaluation of the Risk of Marine Slope Instability: A Pseudo-3D Approach for Application to Large Areas

This article presents a methodology developed to evaluate the instability of submarine slopes that extend over a large area. Special attention was paid to (1) the complex geometry (bathymetry) and the expanse of the slope, (2) the heterogeneity of the sediment, and (3) the distribution of the pore pressure. The safety factor was considered as a spatially varying quantity. The General Formulation (GLE, Fredlund and Krahn 1977), which fully satisfies equilibrium conditions, was used for evaluating the stability of the marine slope. The submarine slope failure, which occurred on 16 October 1979 during the construction of the new Nice airport, was studied in order to test the developed model. Geotechnical parameters were taken from experimental tests carried out by IFREMER on 19 cores extracted at different depths (from 27 m to 1300 m) (Cochonat, Bourillet, and Savoye, 1993; Mulder et al., 1994). Many scenarios were proposed in order to explain the cause of the Nice slope failure (Habib, 1994). In this article, two of those scenarios were tested. Simulation results are presented and discussed.