波浪问题中唯一可解的高阶边界元方法

本文就波浪与结构物相互作用问题，提出了一个适用于高阶边界元应用的新的积分方程，并利用修改积分区域的方法得了适用于本积分方程的不规则频率消除方法。最后，通过数值计算对附加区域的选择、单元的离散做了研究     

Numerical Investigation on Submerged Horizontal Plate

Hydrodynamic characters on a horizontal, thin, rigid plate located beneath the free surface are numerically investigated. Assuming a linear, time-harmonic potential flow and utilizing Green identity, the governing Laplace equation can be simplified into Fredholm integral equation ofthe second kind. Supposing linear-order discontinuous elements along intersecting vertical boundaries, and by use of the boundary element method, numerical solution about source strength distribution on the plate can be changed into a series of algebraic equations. The 3D Green function is introduced to set up the integral equations, and the GMRES solver is performed for solving the large dense linear system of equations. The added-mass, damping force and exciting force are evaluated directly from the equations. It is found that the added-mass coefficient becomes negative for a range of frequencies when the plate is sufficiently close to the free surface.     

Removing irregular frequencies by a partial discontinuous higher order boundary element method

The phenomenon of irregular frequencies is a puzzle in the course of calculating the interaction of waves and structures by the boundary element method. To remove the irregular frequencies, the modified integral domain method is adopted, and continuous higher order elements and partial discontinuous higher order elements are used for discretization. By these means, the effects of the irregular frequencies are effectively removed. Effective strategies have been adopted to deal with singular integrals and nearly singular integrals at different situations. The numerical results of the horizontal wave force on a uniform cylinder in the first order and second order diffraction problems show that the present method has a good validity. At the same time, the influence of collocation parameter on accuracy of numerical results is examined in detail.     

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

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

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

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

A higher-order σ-coordinate non-hydrostatic model for nonlinear surface waves

A higher-order non-hydrostatic model in a σ-coordinate system is developed. The model uses an implicit finite difference scheme on a staggered grid to simultaneously solve the unsteady Navier-Stokes equations (NSE) with the free-surface boundary conditions. An integral method is applied to resolve the top-layer non-hydrostatic pressure, allowing for accurately resolving free-surface wave propagation. In contrast to the previous work, a higher-order spatial discretization is utilized to approximate the large horizontal pressure gradient due to steep surface waves or rapidly varying topographies. An efficient direct solver is developed to solve the resulting block hepta-diagonal matrix system. Accuracy of the new model is validated by linear and nonlinear standing waves and progressive waves. The model is then used to examine freak (extreme) waves. Features of downshifting focusing location and wave asymmetry characteristics are predicted on the temporal and spatial domains of a freak wave.     

Numerical Analysis of Hydrodynamic Pressure Induced by Fluid-Solid Impact

—As a further development of the authors'work(Huang and Qian,1993),in this paper a newnumerical method based on the time domain boundary element technique is proposed for solving fluid-sol-id coupling problems,in which a rigid body impacts normally on the calm surface of a half-space fluid.Afundamental solution to the half-space potential flow problem is first derived with the method of images.Then,an equivalent boundary integral equation in the Laplace transform domain is established by meansof Green's second identity.Through the inverse Laplace transform and discretization in both time andboundary of the fluid region,the numerical calculation for the problem under consideration has been car-ried out.Several examples demonstrate that the present method is more efficient than existing ones,fromwhich it is also seen that the shape of the impacting body has a considerable effect on the total impactforce.     

Fully Nonlinear Boussinesq-Type Equations with Optimized Parameters for Water Wave Propagation

For simulating water wave propagation in coastal areas, various Boussinesq-type equations with improved properties in intermediate or deep water have been presented in the past several decades. How to choose proper Boussinesq-type equations has been a practical problem for engineers. In this paper, approaches of improving the characteristics of the equations, i.e. linear dispersion, shoaling gradient and nonlinearity, are reviewed and the advantages and disadvantages of several different Boussinesq-type equations are compared for the applications of these Boussinesq-type equations in coastal engineering with relatively large sea areas. Then for improving the properties of Boussinesq-type equations, a new set of fully nonlinear Boussinseq-type equations with modified representative velocity are derived, which can be used for better linear dispersion and nonlinearity. Based on the method of minimizing the overall error in different ranges of applications, sets of parameters are determined with optimized linear dispersion, linear shoaling and nonlinearity, respectively. Finally, a test example is given for validating the results of this study. Both results show that the equations with optimized parameters display better characteristics than the ones obtained by matching with padé approximation.     

Multi-resolution MPS for incompressible fluid-elastic structure interactions in ocean engineering

The paper presents a multi-resolution MPS (Moving Particle Semi-implicit)-based FSI (Fluid-Structure Interaction) solver for efficient and accurate simulations of incompressible fluid flows interacting with elastic structures. The fluid model is founded on the projection-based MPS solution of continuity and Navier-Stokes equations. The structure model is set based on MPS-based discretization of linear and angular momenta corresponding to an isotropic elastic solid. Fluid-structure coupling is conducted in a mathematically-physically consistent manner along with implementation of a multi-resolution scheme comprising of common radius of influence, revised weight function, revised number density and potential number density concepts to enhance i) consistency of particle-based discretizations, ii) imposition of boundary conditions and iii) volume conservation at fluid-structure interface. A set of previously developed enhanced schemes are also adopted for the fluid model. The robustness and efficiency of proposed Enhanced Multi-resolution MPS-MPS FSI solver are investigated through a set of ocean engineering-related benchmark tests. To the best knowledge of authors, this study presents the first multi-resolution particle method for FSI corresponding to incompressible fluid and elastic structures.     

求解反应扩散方程的紧致隐积分因子方法

积分因子(IF)方法是近年来提出的求解刚性常微分方程组的一种有效的数值方法。本文应用改进的紧致隐积分因子(cIIF)方法求解二维反应扩散方程。在空间离散上采用二阶中心差分方法。数值模拟得到各种斑图结构与所引文献结果相当一致。     

Two phase modal analysis of nonlinear sloshing in a rectangular container

Sloshing, or liquid free surface oscillation, in containers has many important applications in a variety of engineering fields. The modal method can be used to solve linear sloshing problems and is the most efficient reduced order method that has been used during the previous decade. In the present article, the modal method is used to solve a nonlinear sloshing problem. The method is based on a potential flow solution that implements a two-phase analysis on sloshing in a rectangular container. According to this method, the solution to the mass conservation equation, with a nonpenetration condition at the tank walls, results in velocity potential expansion; this is similar to the mode shapes used in modal method. The kinematic and dynamic boundary conditions create a set of two-space-dimensional differential equations with respect to time. The numerical solution of this set of differential equations, in the time domain, predicts the time response of interfacial oscillations. Modal method solutions for the time response of container sloshing due to lateral harmonic oscillations show a good agreement with experimental and numerical results reported in the literature.     

Investigation of motions of catamarans in regular waves—II

By extending the linear frequency domain theory, a quasi-non-linear time-domain technique has been developed to investigate the large amplitude motions of catamarans in regular waves. The non-linearity of hydrodynamic forces included in this practical method comes from variations of a ship's submerged portion. These forces are obtained from a database generated by the linear frequency domain method at each time step. The coupled equations, heave and pitch, are solved in the time domain by using the Runge-Kutta method with proper initial values. In order to investigate the non-linear effects of large amplitude motions of the V-1 catamaran in the head-sea condition, numerical results obtained from the linear and non-linear strip methods have been compared with those obtained from a series of experiments carried out in the towing tank of the Hydrodynamics Laboratory at the University of Glasgow. Based on the comparative studies, the numerical results obtained from the time-domain program can provide better predictions for the large amplitude motions of catamarans than the linear frequency domain method. It is concluded that the non-linear effects are significant when the model speeds and wave amplitudes increase. The peak values of large amplitude motions around the resonance frequencies, as obtained from the non-linear time-domain predictions as well as from measurements, are smaller than those obtained from the linear theory.     

Internal wave generation in an improved two-dimensional Boussinesq model

A set of Boussinesq-type equations with improved linear frequency dispersion in deeper water is solved numerically using a fourth order accurate predictor-corrector method. The model can be used to simulate the evolution of relatively long, weakly nonlinear waves in water of constant or variable depth provided the bed slope is of the same order of magnitude as the frequency dispersion parameter. By performing a linearized stability analysis, the phase and amplitude portraits of the numerical schemes are quantified, providing important information on practical grid resolutions in time and space. In contrast to previous models of the same kind, the incident wave field is generated inside the fluid domain by considering the scattered wave field in one part of the fluid domain and the total wave field in the other. Consequently, waves leaving the fluid domain are absorbed almost perfectly in the boundary regions by employment of damping terms in the mass and momentum equations. Additionally, the form of the incident regular wave field is computed by a Fourier approximation method which satisfies the governing equations accurately in water of constant depth. Since the Fourier approximation method requires an Eulerian mean current below wave trough level or a net mass transport velocity to be specified, the method can be used to study the interaction of waves and currents in closed as well as open basins. Several computational examples are given. These illustrate the potential of the wave generation method and the capability of the developed model.     

Internal wave generation in an improved two-dimensional Boussinesq model

A set of Boussinesq-type equations with improved linear frequency dispersion in deeper water is solved numerically using a fourth order accurate predictor-corrector method. The model can be used to simulate the evolution of relatively long, weakly nonlinear waves in water of constant or variable depth provided the bed slope is of the same order of magnitude as the frequency dispersion parameter. By performing a linearized stability analysis, the phase and amplitude portraits of the numerical schemes are quantified, providing important information on practical grid resolutions in time and space. In contrast to previous models of the same kind, the incident wave field is generated inside the fluid domain by considering the scattered wave field in one part of the fluid domain and the total wave field in the other. Consequently, waves leaving the fluid domain are absorbed almost perfectly in the boundary regions by employment of damping terms in the mass and momentum equations. Additionally, the form of the incident regular wave field is computed by a Fourier approximation method which satisfies the governing equations accurately in water of constant depth. Since the Fourier approximation method requires an Eulerian mean current below wave trough level or a net mass transport velocity to be specified, the method can be used to study the interaction of waves and currents in closed as well as open basins. Several computational examples are given. These illustrate the potential of the wave generation method and the capability of the developed model.     

A Compact Explicit Difference Scheme of High Accuracy for Extended Boussinesq Equations

Presented here is a compact explicit difference scheme of high accuracy for solving the extended Boussinesq equations.For time discretization,a three-stage explicit Runge-Kutta method with TVD property is used at predicting stage,a cubic spline function is adopted at correcting stage,which made the time discretization accuracy up to fourth order;For spatial discretization,a three-point explicit compact difference scheme with arbitrary order accuracy is employed.The extended Boussinesq equations derived by Beji and Nadaoka are solved by the proposed scheme.The numerical results agree well with the experimental data.At the same time,the comparisons of the two numerical results between the present scheme and low accuracy difference method are made,which further show the necessity of using high accuracy scheme to solve the extended Boussinesq equations.As a valid sample,the wave propagation on the rectangular step is formulated by the present scheme,the modelled results are in better agreement with the experimental data than those of Kittitanasuan.     

A finite volume discretization of the pressure gradient force using analytic integration

Layered ocean models can exhibit spurious thermobaric instability if the compressibility of sea water is not treated accurately enough. We find that previous solutions to this problem are inadequate for simulations of a changing climate. We propose a new discretization of the pressure gradient acceleration using the finite volume method. In this method, the pressure gradient acceleration is exhibited as the difference of the integral “contact” pressure acting on the edges of a finite volume. This integral “contact” pressure can be calculated analytically by choosing a tractable equation of state. The result is a discretization that has zero truncation error for an isothermal and isohaline layer and does not exhibit the spurious thermobaric instability.     

Two sets of higher-order Boussinesq-type equations for water waves

Based on the classical Boussinesq model by Peregrine [Peregrine, D.H., 1967. Long waves on a beach. J. Fluid Mech. 27 (4), 815–827], two parameters are introduced to improve dispersion and linear shoaling characteristics. The higher order non-linear terms are added to the modified Boussinesq equations. The non-linearity of the Boussinesq model is analyzed. A parameter related to h/L₀ is used to improve the quadratic transfer function in relatively deep water. Since the dispersion characteristic of the modified Boussinesq equations with two parameters is only equal to the second-order Padé expansion of the linear dispersion relation, further improvement is done by introducing a new velocity vector to replace the depth-averaged one in the modified Boussinesq equations. The dispersion characteristic of the further modified Boussinesq equations is accurate to the fourth-order Padé approximation of the linear dispersion relation. Compared to the modified Boussinesq equations, the accuracy of quadratic transfer functions is improved and the shoaling characteristic of the equations has higher accuracy from shallow water to deep water.     

An integro-variational method for interior and exterior free surface flow problems

An integro-variational method is used to solve free surface problems of linear potential flow. Results obtained by the proposed method are compared with solution of the finite element formulation and the boundary integral equation. The I.V. method uses isoparametric element distributed on the contour of the fluid domain.     

Wave sloshing in corrugated bottom tanks with two-dimensional flow

The boundary integral element method based on Green's formula is applied to the analysis of transient flow problem in corrugated bottom tanks. The problem is formulated as a two-dimensional linear, initial boundary value problem in terms of a velocity potential. The Laplace equation and the boundary conditions, except the dynamic boundary condition on the free surface, are transformed into an integral equation by the application of Green's formula. Finite Difference discretization is applied timewise. Initially a triangular wave on the free surface is assumed to be formed. The height of the triangular corrugated bottom is varied between 1/10 and 1/5 of the tank depth. The form of the free surface and the equipotential lines for the flow in the tank are presented at different time steps. An accuracy analysis is performed and distortion in time is considered. Proper coefficients for solutions are derived and presented. The results show that utilization of triangular corrugated bottoms may help to regulate the flow in tanks.