Coupled SPHS–BEM method for transient fluid–structure interaction and applications in underwater impacts

Coupled SPHS–BEM method is proposed for transient fluid–structure interaction problems: SPH shell (SPHS) is selected to discretize shell structures, the second-order doubly asymptotic approximations (DAA₂) of boundary element method (BEM) is chosen to analyze flow-field. BEM can remedy the expensive costs for three-dimensional SPH (smoothed particle hydrodynamics), yet SPHS provides a structural solver for BEM. The coupled method is attractive, since only a layer of SPHS particles and a piece of flow-field boundary elements are needed to be modeled; the compatibility conditions of the coupled surface are performed with moving least square (MLS) function. The final two benchmarks on underwater impacts prove the feasibility, stability and accuracy of the proposed method.     

A modified scaled boundary finite-element method for problems with parallel side-faces. Part II. Application and evaluation

The modified scaled boundary finite-element method (SBFEM), keeping the advantages of the original SBFEM, eliminates the restriction of the scaling center location so that this approach can solve two-dimensional problems with parallel side-faces. In this paper, the modified SBFEM is applied to solutions of two types of problems—wave diffraction by a single and twin surface rectangular obstacles and wave radiation induced by an oscillating mono-hull and twin-hull structures in a finite depth of water. For wave diffraction problems, numerical results agree extremely well with the analytic solution for the single obstacle case and other numerical results of a different approach for the twin obstacle case. For wave radiation problems, the particular solutions to the scaled boundary finite-element equation are presented for cases of heave, sway and roll motions. The added mass and damping coefficients for heave, sway and roll motions of a two-dimensional rectangular container are computed and the numerical results are compared with those from independent analytical solution and numerical solution using the boundary element method (BEM). It is found that the SBFEM method achieves equivalent accuracy to the conventional BEM with only a few degrees of freedom. In the last example, wave radiation by a two-dimensional twin-hull structure is analyzed. Comparisons of the results with those obtained using conventional Green's function method (GFM) demonstrate that the method presented in this paper is free from the irregular frequency problems.     

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

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

非线性波浪波面追踪的一种新模式

基于Laplace方程的Green积分表达式和波面BemouUi方程所建立的非线性波动数学模型,是一个时域上具有初始值的边值问题,而精确地追踪自由表面的波动位置,给出波面运动瞬时的波面高度和波面势函数,是建立时域内非线性波浪数值模式的基础。本文采用0-1混合型边界元剖分计算域边界并离散Laplace方程的Green积分表达式,采用有限元剖分自由水面并推导满足自由表面非线性边界条件的波面有限元方程,联立计算域内以节点波势函数和波面位置高度的时间增量为未知量的线性方程组,通过时步内的循环迭代,给出每个时步上的波面位置和波面势函数,从而建立了一种新的非线性波浪波面追踪模式。数值造波水槽内的波浪试验表明,其数值模拟结果具有良好的计算精度。     

Wave loading on axisymmetric bodies using axisymmetric hybrid integral equation method

The wave diffraction problem on axisymmetric structures are solved by treating the fluid field as two separate domains. The velocity potential in the inner domain is represented by a 1/r type Green's function whilst that of the outer domain is represented by an eigenfunction expansion. The simple form of the Green's function in the inner domain reduces significantly the computational effort whilst the eigenfunction expansion in the outer domain is able to satisfy the radiation boundary condition completely. The method requires to have elements cover the entire containing boundary. Results for a number of typical structural geometries are presented and discussions are made on the effect of various parameters.     

An iterative scheme for two-dimensional supercavitating flow

In the present research, supercavitating potential flow is studied numerically by the boundary element method (BEM). Using the advantages of BEM, an iterative algorithm has been introduced to capture cavity boundary in two-dimensional symmetric flows. In this algorithm, the cavity length is known and used to find the related cavitation number and cavity profile. In order to obtain finite length cavities, a cusped cavity closure model has been employed. Applying this cavity closure model, it is possible to change the cavity closure profile and its specified length. By comparing the results of the present analysis with previous analytical and numerical solutions as well as the experimental data, it can be concluded that the present iterative numerical algorithm is reliable and can be applied with BEM or other numerical methods to predict the characteristics of a supercavitating flow. Moreover, the feasibility of the cavity capturing in a flow field with low cavitation number is especially attractive.     

Floating pontoon breakwaters

The hydrodynamic properties of a pair of long floating pontoon breakwaters of rectangular section are investigated theoretically. The structures are partially restrained by linear symmetric moorings fore and aft. The fluid motion is idealized as linearized, two-dimensional potential flow. The breakwater motions are assumed to be two-dimensional, in surge, heave and pitch. The solution for the fluid motion is obtained by the boundary integral equation method using an appropriate Green's function. Numerical results are presented that illustrate the effects of the various wave and structural parameters on the efficiency of the breakwaters as barriers to wave action. It is found that the wave reflection properties of the structures depend strongly on their width, draft and spacing and the mooring line stiffnesses, while their excess buoyancy is of lesser importance.     

The hydrodynamic coefficients for an oscillating rectangular structure on a free surface with sidewall

Based on a two dimensional linear water wave theory, the boundary element method (BEM) is developed and applied to study the heave and the sway problem of a floating rectangular structure in water to finite depth with one side of the boundary is a vertical sidewall and the other boundary is an open boundary. Numerical results for the added mass and radiation damping coefficients are presented. These coefficients are not only depend on the submergence and the width of the structure, but also depend on the clearance between structure and sidewall. Negative added mass and sharp peaks in the damping and added mass coefficients have been found when the clearance with a value close to integral times of half wave length of wave generated by oscillation structure. The important effect of the clearance on the added mass and radiation damping coefficients are discussed in detail. An analytical solution method is also presented. The BEM solution is compared with the analytical solution, and the comparison shows good agreement.     

Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

两层流体中多个圆柱浮体波浪力的解析方法

采用解析方法研究了线性入射波作用下两层流体中多个圆柱形淹没浮体的渡浪力特性.首先基于多极子展开方法,建立了散射势函数的解析表达式,并进一步得到浮体散射渡浪力的计算公式,然后利用边界元方法验证了本文的解析解,最后分析了不同参数的变化对双圆柱形浮体结构波浪力的特有影响.     

Scaled Boundary Finite Element Analysis of Wave Passing A Submerged Breakwater

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique combining the advantage of the finite element method （FEM） and the boundary element method （BEM） with its unique properties. In this paper, the SBFEM is used for computing wave passing submerged breakwaters, and the reflection coeffcient and transmission coefficient are given for the case of wave passing by a rectangular submerged breakwater, a rigid submerged barrier breakwater and a trapezium submerged breakwater in a constant water depth. The results are compared with the analytical solution and experimental results. Good agreement is obtained. Through comparison with the results using the dual boundary element method （DBEM）, it is found that the SBFEM can obtain higher accuracy with fewer elements. Many submerged breakwaters with different dimensions are computed by the SBFEM, and the changing character of the reflection coeffcient and the transmission coefficient are given in the current study.     

A numerical modeling of hydrodynamic characteristics of various planing hull forms

A numerical algorithm based on the boundary element method (BEM) is presented for predicting the hydrodynamic characteristics of the various planing hull forms. The boundary integral equation is derived using Green's theorem on the wetted body surface and the free surface. The ventilation function at the transom is estimated with Doctor's empirical formula. This function is defined as the transom zone free surface boundary condition. The combined boundary integral equation and modified free surface boundary condition are simultaneously solved to determine the dipole on the wetted hull surface and the source on the free surface. The method is applied to investigate three examples of planing hulls, which include flat-plates, as well as wedge-shaped and variable deadrise planing hulls. Their hydrodynamic characteristics are calculated for different speeds. Computational results are presented and compared with existing theories and experiments. On the whole, the agreement between the present method and the selected experimental and numerical data is satisfactory.     

A simple two-way coupling method of BEM and VOF model for random wave calculations

A numerical method, which combines the boundary element method (BEM) and the volume of the fluid method (VOF method), has been presented to solve wave–structure interactions; the intense wave motion at the proximity of the structure is modeled by the VOF method and the rest of the fluid region is modeled by the BEM. The combined method can considerably reduce the time-consuming VOF domain, and thus practically makes it possible to apply the VOF method for random wave calculations, in which long time computations are usually required to obtain statistically meaningful results, and therefore the use of the single-VOF model often becomes prohibitive in terms of computational time and storage memories. A VOF model CADMAS-SURF, which is based on SMAC scheme and had been constructed by a number of VOF researchers in coastal engineering in Japan, is used in the combined BEM–VOF model. The two-way coupling treatment, which enables us to deal with bidirectional wave propagations, which was originally given for the SOLA-VOF model by Yan et al. (2003a) and later improved by Kim et al. (2007), was modified for the SMAC scheme. The coupling treatments are described in detail in the paper. The validity of the combined BEM–VOF model was investigated by comparing the numerical results with the theoretical results for the propagations of Stokes 5th order waves and random waves.     

海底管道冲刷势流模型的改进

对Li&Cheng的势流模型在数值方法、床面平衡条件及冲刷床面调整技术做出了3方面的改进。采用边界元法代替差分法求解Laplace方程,前者可以准确地拟合地形与管道边界,因而可以准确反映固壁边界对流态的影响,此外,还具有数据准备简单,降低计算维数,计算速度快等优点。Li&Cheng模型以床面水流切应力等于泥沙起动切应力,τb=τc,作为床面平衡条件,这只适用于清水冲刷。以沿程输沙相同作为平衡剖面条件,理论上,将模型推广到了动床冲刷。此外,为了提高模型的收敛性,提出了最速下降法与牛顿迭代法相结合的床面调整技术。采用实验资料对模型进行了检验,计算的冲刷深度与实验结果符合良好。     

应用比例边界有限元法求解狭缝对双箱水动力的影响

比例边界有限元法(SBFEM)是一种半解析数值分析的新方法,既融合了有限元法和边界元法的优点,又有其特有的优点。用该方法可求解有限水深下狭缝对双箱水动力作用的影响,为波浪与多浮体超大型结构的相互作用探索一些规律。整个计算域划分成2个无限子域和4个有限子域,并利用加权余量法在各个子域上推导了SBFEM的积分方程;计算了4个数值算例并与边界元等其它数值方法进行了比较,验证了该方法是一种用很少单元便能得到精确结果的高效方法。应用SBFEM对不同箱体宽度、不同狭缝宽度、不同吃水深度条件的双箱作了计算,得出了狭缝对双箱水动力干涉影响的一些规律,对超大型浮体水动力分析和结构设计具有一定的参考价值。     

Radiation and diffraction of linear water waves by an infinitely long submerged rectangular structure parallel to a vertical wall

The radiation and the diffraction of linear water waves by an infinitely long floating rectangular structure submerged in water of finite depth with leeward boundary being a vertical wall are analyzed in this paper by using the method of separation of variables. Analytical expressions for the radiated and diffracted potentials are derived as infinite series with unknown coefficients determined by the eigenfunction expansion matching method. The expressions for wave forces and hydrodynamic coefficients are given. A comparison is made between the results obtained by the present analytical solution and those obtained by the boundary element method. By using the present analytical solution, the hydrodynamic influences of the submergence, the width, the thickness of the structure, and the distance between the structure and the wall on the wave forces and hydrodynamic coefficients are discussed in detail.     

Pressure characteristics of bubble collapse near a rigid wall in compressible fluid

High speed liquid jet and shockwave can be produced when a bubble collapses near a rigid wall, which may cause severe damage to solid structures. A hybrid algorithm was adopted to simulate bubble motion and associated pressures near a wall combining Level Set-Modified Ghost Fluid-Discontinuous Galerkin (LS-MGF-DG) method and boundary element method (BEM). Numerical results were compared with experimental data to validate the presented algorithm. Jet formation was simulated by BEM and the induced pressure on the wall was calculated with auxiliary function. The pressure at the point on the wall where the jet points to reaches its peak value after the jet penetrates the bubble. Bubble collapse and rebounding were simulated by the LS-MGF-DG method. Shock-wave is induced when the bubble collapse toroidally to a minimum volume and the pressure at wall center reaches the maximum due to shockwave superposition. A third pressure peak is found associated with the bubble rebounds and bubble splitting. In the case studied, a higher pressure was found due to collapse shockwave than bubble jet and affects a larger area of the wall. In addition, the three pressure peaks due to jet impact, collapse impact as well as bubble rebounding and splitting decrease with the increase of the standoff distance.     

A fast multipole boundary element method for three-dimensional potential flow problems

A fast multipole methodology (FMM) is developed as a numerical approach to reduce the computational cost and memory requirements in solving large-scale problems. It is applied to the boundary element method (BEM) for three-dimensional potential flow problems. The algorithm based on mixed multipole expansion and numeric, al integration is implemented in combination with an iterative solver. Numerical examinations, on Dirichlet and Neumann problems, are carried out to demonstrate the capability and accuracy of the present method. It has been shown that the method has evident advantages in saving memory and computing time when used to solve huge-scale problems which may be prohibitive for the traditional BEM implementation.     

完全非线性波浪与二维固定结构物作用的IGN-BEM耦合分析模型

为高效准确地对完全非线性波浪与二维固定结构物的相互作用进行模拟分析,建立了二维完全非线性时域耦合模型。耦合模型将计算域划分为靠近结构物的内域和远离结构物的外域,每个区域均采用满足完全非线性自由水面边界条件的波浪模型进行求解。在内域使用Laplace方程描述流体运动并采用高阶边界元法（BEM）对其进行求解;而在没有结构物的外域,波浪运动的控制方程为Irrotational Green-Naghdi（IGN）方程并采用有限元法（FEM）对其进行求解。内域和外域通过一段重叠区域进行耦合,从而实现模型间变量的传递。首先利用耦合模型分别对规则波的传播、直墙前立波的生成以及相关物理模型试验进行模拟,数值结果与精确解和试验结果的良好吻合验证了耦合模型耦合方式的合理性以及处理非线性问题的准确性;然后使用耦合模型模拟分析了波浪与固定结构物间的相互作用,并将结果与线性解析解以及完全非线性BEM模型的结果进行了对比分析,进一步证明了耦合模型的正确性与高效性。     

A B-spline-based method for radiation and diffraction problems

An open uniform B-spline-based panel method is developed for solution of potential flow problems. In this method, both geometry as well as the field variables are represented by the same open uniform B-spline basis function. The method is initially applied for the radiation problem in unbounded fluid. Computed results for a spheroid of different aspect ratio are found to be in excellent agreement with analytical results. The method is then applied for diffraction problem formulated based on the transient (time-domain) Green's function. Computed results for a hemisphere and Wigley hull are compared with published results and the comparison shows good agreement.