FINITE ELEMENT ANALYSIS OF A VERTICAL RECTANGULAR PLATE COUPLED WITH AN UNBOUNDED FLUID DOMAIN ON ONE SIDE USING A TRUNCATED FAR BOUNDARY

The dynamic pressure distribution on a rectangular plate attached to a rigid wall and supporting an infinitely large extent of fluid subjected to a harmonic ground excitation is evaluated in the time domain. Governing equations for the fluid domain are set considering the compressibility of the fluid with negligibly small change in density and a linearized free surface. A far boundary condition for the three-dimensional fluid domain is developed so that the far boundary is truncated at a closer proximity to the structure. The coupled problem is solved independently for the structure and the fluid domain by transferring the acceleration of the plate to the fluid and pressure of the fluid to the plate in sequence. Helmholtz equation for the three-dimensional fluid domain and Mindlin's theory for the two-dimensional plate are used for the solution of the interacting domains. Finite element technique is adopted for the solution of this problem with pressure as nodal variable for the fluid domain and displacement for the plate. The time dependent equations are solved in each of the interacting domain using Newmark-b method. The effectiveness of the technique is demonstrated and the influences of surface wave, exciting frequency and flexibility of the plate on dynamic pressure are investigated.     

THREE-DIMENSIONAL HYDROELASTIC ANALYSIS IN TIME DOMAIN WITH APPLICATION TO AN ELASTIC SHIP MODEL

According to the Price-Wu condition at the interface between flexible marine structure and the surrounding fluid flow, a solution to the three-dimensional potential flow in time domain around a flexible structure traveling in waves is expressed with a boundary integral equation in the present paper, where the Green function, satisfying the linearized free surface condition for the time-dependent problem, is employed. A hydroe-lastic analysis directly in time domain to predict loads, motions, and structural responses of ships at steady forward speed in a seaway is formulated. As a example of application, the comparisons are made among the numerical results given by present method, the experimental measurements and the prediction provided by the three-dimensional hydroelas-ticity theory in frequency domain.     

用时域法求解二维二阶非线性水波的绕射问题

本文用二阶理论在时期围内研究非线性水波的绕射问题。一、二阶问题分别满足各自的自由表面条件和物面条件。用一阶Ａｄａｍｓ格式并配合时步处理来计算自由表面上的波高和速度势，而流场内的速度势通过求解有限元方法得到。采用人工阻尼区来吸收反射波，对水平半圆柱在水面的绕射问题进行了计算，得到了自由表面波高时间历程和圆柱所受到的波浪力，计算结果和有关文献进行了比较，其符合程度是令人满意的。     

应用时域格林函数方法模拟有限水深中波浪对结构物的作用

该文采用时域内满足自由水面条件的格林函数方法对无限和有限水深中波浪与任意形状三维物体相互作用所产生的绕射和辐射问题进行了理论研究和数值计算。建立了不包含自由水面积分的积分方程,消除了自由项系数。文中计算了各种形状三维物体在无限水深和有限水深情况下所受到的绕射和辐射波浪力,并与解析解以及频域理论进行了对比,结果符合良好,验证了本文所提理论的正确性。     

NUMERICAL SIMULATIONS OF WAVE-INDUCED SHIP MOTIONS IN TIME DOMAIN BY A RANKINE PANEL METHOD

A time domain prediction of wave-induced ship motions by a Rankine panel method is investigated. Linear boundary conditions on free surface and mean wetted body surface are adopted, while the numerical damping method is used for the radiation conditions. The motions of two ships in regular head waves are computed by the present method. The related numerical results are compared with the experiment data and those from linear strip theory. The comparison shows satisfactory agreements for pitch and heave transfer functions.     

UNSTEADY FREE-SURFACE WAVES GENERATED BY BODIES IN A VISCOUS FLUID

BRIEFINTRODUCTIONOFTHEPAPER :　Inthefirstpartofthisthesis ,thefundamentalsolutionsforthesingularStokesandOseenflowsinanun boundedfluidarederivedinauniversalformwhichinvolvestheHamiltonian ,Hessian ,andLaplacianoperators ,andelementaryfunctions .Inthesecondpart,theinteractionofunsteadylow Reynolds numberflowswithafreesurfaceisinvestigatedanalytically .Thedisturbedflows ,generatedbysubmergedbodiesmovingverticallydownwardsawayfromthesurfaceofthefluid ,aregovernedbytheunsteadyStokeseq…     

Finite element analysis of second order wave radiation by a group of cylinders in the time domain

A finite element based numerical method is employed to analyze the wave radiation by multiple or a group of cylinders in the time domain. The nonlinear free surface and body surface boundary conditions are satisfied based on the perturbation method up to the second order. The first- and second-order velocity potential problems at each time step are solved through a Finite Element Method (FEM). The matrix equation of the FEM is solved through iteration and the initial solution is obtained from the result at the previous time step. The three-dimensional (3-D) mesh required is generated based on a two-dimensional (2-D) hybrid mesh on a horizontal plane and its extension in the vertical direction. The hybrid mesh is generated by combining an unstructured grid away from cylinders and two structured grids near the cylinder and the artificial boundary. The fluid velocity on the free surface and the cylinder surface are calculated by using a differential method. Results for various configurations including the cases of two cylinders and four cylinders and a group of eighteen cylinders are obtained to show the joint influences of cylinders on the first- and second-order waves and forces, including the effects of spacing ratios and wave frequency on the second order waves and the mean force, in particular.     

COMPUTATIONAL RESEARCH ON WAVE MAKING OF MOVING WIGLEY HULL IN TIME DOMAIN

Based on Green＇s theorem, a time domain numerical model was constructed to simulate wave making phenomenon caused by a moving ship. In this article, the Rankine sources and dipoles were placed on boundary surfaces （i.e., the ship surface and free surface）, and a time-stepping scheme was employed. Its unique characteristic is that steady state can be realized from initial value by employing the time-stepping scheme and unsteady free surface conditions. In time domain, if the results of unsteady flow problem tend to data stabilization after many time steps of computation, they could be regarded as the data of steady ones. This model could be employed to steady or unsteady problems. Theoretical reasoning and computational process of this method was described in detail The linear and nonlinear boundary conditions on body surface were studied, and the relative means to realize these boundary conditions in iterative computation were also discussed. Some proper parameters about the model of the Wigley hull were determined by many numerical tests, and their influences on wave making resistance and wave pattern were discussed. According to the comparison between numerical results and data available in relative references, the method used in this work is proven to be a reliable method in time domain. And the lattice reorganization in every time step computation is a feasible numerical approach.     

浮体的非线性辐射效应

本文在空间上应用面源法，在时间域上应用步进法，研究了浮体的非线性辐射效应，本文探讨了：（1）不同形式的辐射条件对非线性辐射问题解的影响；（2）不同距离的辐射面对开边界条件的使用所产生的影响；（3）大振幅的浮体受近运动所产和的非线性效应。     

WAVE MAKING COMPUTATION IN TIME DOMAIN FOR MULTI-HULL SHIPS

1. INTRODUCTION For high-speed displacement ships, the wave making at the free surface is an important hydro- dynamic performance. The wave making reduction has attracted more attention for ship designers. The ship wave resistance and ship hull form improvement have been widely investigated. The thin ship theory adopting the Havelock source distribution on the ship centerline plane can be used to analyze the ship wave making profile and the ship wave resistance may be obtained by ship wav…     

THE WAVE-MAKING CHARACTERISTICS OF A MOVING BODY IN A TWO-LAYER FLUID

The Wave-making characteristics of a moving body in a two-layer fluid with free surface is investigated numerically and experimentally. The numerical analysis is based on the modified layered boundary integral equation system. The wave characteristics on the free surface and interface generated by a moving sphere and an ellipsoid is numerically simulated in both finite depth and infinite depth of lower layer model. The numerical results of the sphere are compared with the analytical results for a dipole with the same velocity in a two-layer fluid of finite depth. The dependence of the wave systems and structures on the characteristic quantities is discussed. Three kinds of measurement techniques are used in model experiments on the internal waves generated by a sphere advancing in a two-layer fluid. The effects of the varying velocity and stratification on the wavelength, wave amplitudes and the maximum half angles of internal waves are analyzed qualitatively.     

Iterative Rankine HOBEM analysis of hull-form effects in forward-speed diffraction problem

A time-domain numerical algorithm based on the higher-order boundary element method and the iterative time-marching scheme is proposed for seakeeping analysis. The ship waves generated by a hull advancing at a constant forward speed in incident waves and the resultant diffraction forces acting on the hull are computed to investigate the hull-form effects on the hydrodynamic forces. A rectangular computational domain travelling at ship's speed is considered. An artificial damping beach for satisfying the radiation condition is installed at the outer portion of the free surface except the downstream side. An iterative time-marching scheme is employed for updating both kinematic and dynamic free-surface boundary conditions for numerical accuracy and stability. The boundary integral equation is solved by distributing higher-order boundary elements over the wetted body surface and the free surface. The hull-form effects on the naval hydrodynamics are investigated by comparing three different Wigley models. Finally, the corresponding unsteady wave patterns and the wave profiles around the hulls are illustrated and discussed.     

强震下流体对渡槽槽身的作用

本文应用边界元法计算了强震下流体的非线性晃动及其对槽身的水平力及翻转力矩，将所得数值结果与线性解析方法的计算结果作了比较，分析了两类结果的异同点。计算表明：强震时，流体对槽身可能会产生很大的横向水平力，建议在抗震设计时考虑这一作用因素。     

An iterative Rankine BEM for wave-making analysis of submerged and surface- piercing bodies in finite water depth

A 3-D iterative Rankine Boundary Element Method （BEM） for seakeeping problem in time domain is developed in the framework of linear potential theory. Waves generated by both submerged and surface-piercing bodies moving at a constant forward speed in otherwise calm water, and the resultant steady wave pattern, wave profile and resistance are computed to validate this newly-developed code. A rectangular computational domain moving with the same forward speed as the body is introduced, in which an artificial damping beach is installed at an outer portion of the free surface except the downstream side for satisfying the radiation condition. The velocity potential on the ship hull and the normal velocity on the free surface are obtained directly by solving the boundary integral equation, with the Rankine source used as the kernel function. An iterative time-marching scheme is employed for updating both kinematic and dynamic free surface boundary conditions to stabilize the calculation. Extensive results including the wave patterns, wave profiles and wave resistances for a submerged spheroid and a Wigley hull with forward speed are presented to validate the efficiency of the proposed 3-D time-domain higher-order approach. Finally, the sensitivity of ship-generated waves to the water depth is investigated. Computed results show satisfactory agreement with the corresponding experimental data and other numerical solutions.     

Water entry of decelerating spheres simulations using improved ISPH method

In this paper, we simulated the vertical impact of spheres on a water surface using three-dimensional incompressible smoothed particle hydrodynamics(3-D ISPH) method. The sphere motion is taken to be a rigid body motion and it is modeled by ISPH method. The governing equations are discretized and solved numerically using ISPH method. A stabilized incompressible SPH method by relaxing the density invariance condition is adopted. Here, we computed the motions of a rigid body by direct integration of the fluid pressure at the position of each particle on the body surface. The equations of translational and rotational motion were integrated in time domain to update the position of the rigid body at each time step. In this study, we improved the boundary treatment between fluid and fixed solid boundary by using virtual marker technique. In addition, an improved algorithm based on the virtual marker technique for the boundary particles is proposed to treat the moving boundary of the rigid body motion. The force exerted on the moving rigid boundary particles by the surrounding particles, is calculated by the SPH approximation at the virtual marker points. The applicability and efficiency of the current ISPH method are tested by comparison with reference experimental results.     

造波板运动造波的实时模拟

基于势流理论和时域高阶边界元方法,建立了完全非线性数值波浪水槽模型,通过实时模拟造波板运动来产生波浪.采用混合欧拉-拉格朗日方法追踪流体瞬时水面,运用四阶Runga-Kutta方法更新下一时间步的波面和速度势,同时应用镜像格林函数消除了水槽两个侧面和底面上的积分,在水槽末端布置人工阻尼层消除反射波浪,通过加速势的方法准确计算速度势的时间导数.利用所建模型分别对规则波和不规则波进行了模拟研究,得到了稳定的波形,与试验结果、理论结果和已发表结果均吻合良好,同时计算了造波板做单一频率运动产生的非线性波浪力.     

TIME-DOMAIN ANALYSIS OF WAVE FORCES ON SHIPS WITH FORWARD SPEED

Analyses of wave forces on Wigley and Series 60 ships at forward speedare presented in time domain which is based on a free surface transient Green functionwith linearized condition. The main effort focus on the investigation of the numerical sta-bility of time stepping calculation and the waterline intergal contribution to the excitingforces. Numerical study shows that the stable results can be obtained when a wall-sided as-sumption is used at the stern only and the effect of waterline can be neglected. Compar-isons have been made between the present time domain calculations and other works. Nu-merical results for the added resistance on a C_B=0.8 ship in short waves are presentedbased on the assumption that the steady disturbance of ship hull can be neglected.     

Simulation of water entry of an elastic wedge using the FDS scheme and HCIB method

The hydroelasticity of water entry of an elastic wedge is simulated using a code developed using the flux-difference splitting scheme for immiscible and incompressible fluids and the hybrid Cartesian/immersed boundary method. The free surface is regarded as a moving contact discontinuity and is captured without any additional treatment along the interface. Immersed boundary nodes are distributed inside a fluid domain based on the edges that cross an instantaneous body boundary. Dependent variables are reconstructed at each immersed boundary node with the help of an interpolation along a local normal line for providing a boundary condition for a discretized flow problem. A dynamic beam equation is used for modeling the elastic deformation of a wedge. The developed code is validated through comparisons with other experimental and computational results for a free-falling wedge. The effects of the elastic deformation of the wedge on the pressure fields and time histories of the impact force are investigated in relation to the stiffness and density of the structure. Grid independence test is carried out for the computed time history of the force acting on an elastic wedge.     

FULLY NONLINEAR CALCULATION OF 2－D WATER WAVE GENERATED BY ROCKER FLAP WAVEMAKER

ＦＵＬＬＹＮＯＮＬＩＮＥＡＲＣＡＬＣＵＬＡＴＩＯＮＯＦ２－ＤＷＡＴＥＲＷＡＶＥＧＥＮＥＲＡＴＥＤＢＹＲＯＣＫＥＲＦＬＡＰＷＡＶＥＭＡＫＥＲ￥ＨｅＷｕ－ｚｈｏｕ（Ｄｅｐｔ．ｏｆＨｙｄｒａｕｌｉｃＥｎｇ．，ＴｓｉｎｇｈｕａＵｎｉｖｅｒｓｉｔｙ，Ｂｅｉｊｉ...