重力式网箱数值模拟时波浪理论选择

针对重力式网箱数值模拟中波浪理论选择问题进行了专门的分析与讨论。在数值模型中分别采用线性波浪理论和斯托克斯五阶波浪理论,对重力式网箱的主要部件浮架系统和网衣系统进行模拟。通过将模拟结果与试验结果比较分析选取两种不同波浪理论对模拟结果的影响。分析结果认为在相对水深大于0.18情况下,应用线性波浪理论基本可以准确地模拟重力式网箱波浪作用下的运动和受力。     

Interactions between nonlinear water waves and non-wall-sided 3D structures

Interactions between water waves and non-wall-sided cylinders are analyzed based on velocity potential theory with fully nonlinear boundary conditions on the free surface and the body surface. The finite element method (FEM) is adopted together with a 3D mesh generated through an extension of a 2D Delaunay grid on a horizontal plane along the depth. The linear matrix equation for the velocity potential is constructed by imposing the governing equation and boundary conditions through the Galerkin method and is solved through an iterative method. By imposing the gradient of the potential equal to the velocity, the Galerkin method is used again to obtain the velocity field in the fluid domain. Simulations are made for bottom mounted and truncated cylinders with flare in a numerical tank. Periodic waves and wave groups are generated by a piston type wave maker mounted on one end of the tank. Results are obtained for forces, wave profiles and wave runups. Further simulations are made for a cylinder with flare subjected to forced motion in otherwise still open water. Results are provided for surge and heave motion in different amplitudes, and for a body moving in a circular path in the horizontal plane. Comparisons are made in several cases with the results obtained from the second order solution in the time domain.     

GPS测高技术在无验潮水深测量中的应用

应用双频GPS动态后处理高精度定位技术，建立了一套完整的GPS无验潮海洋深度测量作业模式，通过海上试验与传统作业模式作了数值分析比较，结果表明，该作业模式不仅无需验潮，而且能够有效消除传统作业模式中船只动态吃水和涌浪等因素对测量成果的影响，显著提高水深测量成果的精度。     

The Role of GNSS Vertical Velocities to Correct Estimates of Sea Level Rise from Tide Gauge Measurements in Greece

In this study, we show how the Global Navigation Satellite System (GNSS)-derived vertical velocities contribute to the correction of tide gauge (TG) measurements used for the sea level rise estimation in Greece. Twelve sites with records of local sea level heights are processed in order to estimate their trend. Certain error sources related to TGs, e.g. equipment changes, data noise, may lead to biased or erroneous estimations of the sea level height. Therefore, it would be preferred to follow a robust estimation technique in order to detect and reduce outlier effects. The geocentric sea level rise is estimated by taking into account the land vertical motion of co-located GNSS permanent stations at the Hellenic area. TGs measure the height of the water relative to a monitored geodetic benchmark on land. On the other hand, using GNSS-based methods the vertical land motion can be derived. By means of extended models fitted to the GNSS time-series position, obtained from seven years of continuous data analysis, periodic signals are well described. The synergy of the two co-located techniques results in the correction of TG relative sea level heights taking into account the GNSS vertical velocities and consequently obtaining the conversion to absolute (geocentric) sea level trend.     

Wave reflection and transmission by curtainwall–pile breakwaters using circular piles

This paper presents a mathematical model which computes the hydrodynamic characteristics of a curtainwall–pile breakwater (CPB) using circular piles, by modifying the model developed for rectangular piles by Suh et al. [2006. Hydrodynamic characteristics of pile-supported vertical wall breakwaters. Journal of Waterway, Port, Coastal and Ocean Engineering 132(2), 83–96]. To examine the validity of the model, laboratory experiments have been conducted for CPB with various values of draft of curtain wall, spacing between piles, and wave height and period. Comparisons between measurement and prediction show that the mathematical model adequately reproduces most of the important features of the experimental results. The mathematical model based on linear wave theory tends to over-predict the reflection coefficient as the wave height increases. As the draft of the curtain wall increases and the porosity between piles decreases, the reflection and transmission coefficient increases and decreases, respectively, as expected. As the relative water depth increases, however, the effect of porosity disappears because the wave motion is minimal in the lower part of a water column for short waves.     

Dynamic Behavior of Beam–Pile Structure Under Vertical Transient Excitation

The dynamic response of beam–pile–soil system under vertical transient excitation is investigated. Both piles and beam are assumed to be one-dimensional rods and subjected to vertical exciting forces. The uniformly distributed Voigt models are introduced to simulate the pile tip resistances, and the dynamic interactions between piles and beam are simplified as a set of concentrated point loads. Then, the plane strain model, the theory of longitudinal vibration of one-dimensional rod, and the Timoshenko beam theory are used to establish the mathematical models for the motion of soil, piles, and beam, respectively. On this basis, the matrix equation for solving the governing equations is constructed in the Laplace domain and the time-domain response is then obtained by the discrete inverse Fourier transform. Comparisons with numerical simulations and model tests are conducted to evaluate the rationality of the present solution. The results show that the dynamic responses calculated by the proposed solution are generally consistent with simulated curves and experimental data.     

Viscous drift forces in regular and irregular waves

A method to compute wave- and current-induced viscous drift forces and moments on floating platforms in regular and random waves is presented. The relative velocity drag term of Morison's equation is used in conjunction with frequency domain first-order motion transfer functions to compute the drift forces and moments. Mean viscous drift forces and moments in regular waves in all six degrees-of-freedom of a tension leg platform are computed. The relative importance of the free-surface force integration, steady current, wave-current interaction and platform motions on the computed drift forces and moments are discussed. The results from this method, in the frequency domain, are used to compute the drift forces and responses in irregular waves using existing methods developed for potential drift computations. Comparisons with results from time-domain computations are also presented and good agreement between the frequency-domain and time-domain results is found. Some comparisons with experimental data are also made. The frequency-domain method is found to be an efficient and useful tool for the analysis of semi-submersible and tension leg platforms during the preliminary design stage in which extensive parametric studies need to be undertaken.     

Analytical solutions of the diffraction problem of a group of truncated vertical cylinders

An exact analytical method is described to solve the diffraction problem of a group of truncated vertical cylinders. In order to account for the interaction between the cylinders, Kagemoto and Yue's exact algebraic method is utilised. The isolated cylinder diffraction potential due to incident waves is obtained using Garret's solution and evanescent mode solutions are derived in a similar manner.Numerical results are presented for arrays of two and four cylinders. Comparisons between the results obtained from the method presented here and those obtained from numerical methods show excellent agreement.     

Effects of the asymmetric riser and bilge keel arrangements on FPSO green water assessment

This article aims to report important findings on how the asymmetric riser and bilge keel arrangements affect the motion response and green water assessment by using a real FPSO conversion project. Recently, the authors have proposed a practical approach for short-term and long-term green water prediction. In this paper, the method has been further extended to include the effect of truncated bilge keel by using Morrison elements. Numerical studies are conducted focusing on the effect induced by asymmetric riser arrangement and truncated bilge keels. Comparisons of short-term and long-term results between different models indicate that the FPSO’s motion is significantly affected by asymmetrically arranged appendages and attachments in a complicated way. The relative wave elevation is also affected by appendages and attachments, but not the same trend as the motion response. The effect of the asymmetric arrangement of risers and bilge keel on long-term relative wave elevation response has been captured by both traditional contour line approach and response-based analysis, but some discrepancy identified between the results from the two methods indicates the limitation of the traditional contour line approach.     

Axisymmetric liquid block impact on a solid surface

The hydrodynamic problem of collision of an axisymmetric water block and solid boundary has been analyzed based on the velocity potential theory. The three dimensional axisymmetric hydrodynamic problem is converted to a quasi two dimensional one and is solved in axisymmetric coordinate system using the boundary element method. Numerical simulations are made for collision of a solid cone with a liquid cone or a spheroidal liquid block. The results in the former are compared with those obtained from similarity solution and excellent agreement is found. Numerical results are then presented for the collision between different solid cones and various spheroidal liquid blocks.     

Action of Cnoidal Waves on Vertical Walls

Based on the 2nd order cnoidal wave theory, the characters of shallow water standing waves and their action on vertical walls are studied in this paper. The theoretical expressions of the wave surface elevation in front of and the wave pressure on the vertical wall are obtained. In order to verify the theoretical results, model tests were made in the State Key Laboratory of Coastal and Offshore Engineering at DUT. For the wave surface elevation in front of the wall and the wave forces on the wall at the moment when the wave surface at the wall surface goes down to the bottom of the wave trough, the calculated results coincide quite well with the experimental results. For the wave forces on the wall at the moment when the wave surface at the wall surface goes up to the top of the wave crest, the theoretical expressions are modified by the experimental results. For the convenience of practical use, calculations are made for the wave conditions which usually occur in enginering practice by use of the inves     

Finite Element Analysis to Two-Dimensional Nonlinear Sloshing Problems

A two-dimensional nonlinear sloshing problem is analyzed by means of the fully nonlinear theory and time domain second order theory of water waves. Liquid sloshing in a rectangular container subjected to a horizontal excitation is sim-ulated by the finite element method. Comparisons between the two theories are made based on their numerical results, It is found that good agreement is obtained for the case of small amplitude oscillation and obvious differences occur for large amplitude excitation. Even though, the second order solution can still exhibit typical nonlinear features of nonlinear wave and can be used instead of the fully nonlinear theory.     

Numerical prediction of green water incidents

Green water loads on moored or sailing ships occur when an incoming wave significantly exceeds the freeboard and water runs onto the deck. In this paper, a Navier–Stokes solver with a free surface capturing scheme (i.e. the VOF model; [Hirt and Nichols, 1981]) is used to numerically model green water loads on a moored FPSO exposed to head sea waves. Two cases are investigated: first, green water on a fixed vessel has been analysed, where resulting waterheight on deck, and impact pressure on a deck mounted structure have been computed. These results have been compared to experimental data obtained by [Greco, 2001] and show very favourable agreement. Second, a full green water incident, including vessel motions has been modelled. In these computations, the vertical motion has been modelled by the use of transfer functions for heave and pitch, but the rotational contribution from the pitch motion has been neglected. The computed water height on deck has been compared to the experimental data obtained by [Buchner, 1995a] and it also shows very good agreement. The modelling in the second case was performed in both 2-D and 3-D with very similar results, which indicates that 3-D effects are not dominant.     

Finite vertical oscillation of fluid parcel in a heterogeneous,incompressible sea water—Part I

The solutions of vertical oscillation of finite amplitude of a fluid parcel in a heterogeneous, incompressible sea water were obtained, solving rigorously the equation of motion, retaining small terms to the third order of the amplitude non-dimensionized by a vertical reference length. The motion is asymmetric with respect to its equilibrium position and the period of motion is longer than that of infinitesimal Väisälä oscillation. When the vertical velocity of the parcel exceeds certain critical value, oscillation of large amplitude sets in near the outskirts of the thermocline, its period being more than twice larger than that of Väisälä oscillation. The apprehension whether this type of the increased vertical motion might be influenced or denied by the horizontal motion in such a case as in internal waves is also discussed. From some observations which will be appeared in Part II, we found oscillations of large amplitude with the period which is in good agreement with the theoretical result.     

Diagnosis of vertical velocities with the QG Omega equation: a relocation method to obtain pseudo-synoptic data sets

The quantification of vertical motion and vertical fluxes is essential in our ability to predict and tolerate climate change. However, diagnostic estimations might be affected by the errors arising from the necessary compromise between spatio-temporal resolution and cost of hydrographic surveys. Observations of a numerical ocean model have been made in order to test the accuracy of different sampling strategies and their possible a-posteriori corrections. A simple first-order correction method, computing a pseudo-synoptic data set from a non-synoptic data set and involving a geostrophic relocation of the stations is shown to correct significantly the synopticity error in hydrographic data, derived QG vertical motion and vertical temperature fluxes. Sensitivity analyses also show that the lack of synopticity is more critical than other factors, including the sampling resolution, the level of no-motion and the analysis.     

Wave interaction with tension leg platforms

The design of deep water offshore platforms requires the analysis of wave-structure interaction phenomena which have not been as critical for shallower water platform designs. In the case of tension leg platforms (TLPs) interaction phenomena such as wave run-up on the vertical legs and the amplification of the waves beneath the deck are major design considerations. The research investigation reported here focuses on a series of small scale wave tank tests on four column TLP models examining these phenomena. The role of vertical leg spacing and comparative tests of the TLP models with and without pontoons was investigated. As the vertical legs were moved closer an increase in wave run-up and a shifting of the incident wave period corresponding to the maximum wave upwelling were noted. Comparisons with wave measurements for single cylinders from previous experimental studies and the TLP configurations used in this study are presented. A design formula for estimating wave run-up on TLPs is suggested based upon these experiments. The wave run-up on a leg directly in the wake of another leg is presented. A comparison of the wave upwelling measurements with previously published numerical results are discussed. A wave uplift force model which allows for the inclusion of the experimentally obtained wave upwelling measurements is presented and discussed with regard to the design specification of platform deck elevation.     

Drift velocity at the wave breaking point

Wave-tank studies were conducted on the measurement of the drift velocity at the breaking point under different types of breaking waves on a rigid, plane beach. The drift velocity has onshore direction near the surface and close to the bottom; in the main flow column, the drift velocity is always offshore. The offshore drift velocity shows a more uniform vertical distribution than that in the offshore region. The experimental data are compared with theoretical values of three different second-order constant-depth wave theories. Comparisons with data from other sources are also made.     

Experimental study on bubble pulse features under the combined action of horizontal and vertical walls

The pulse features of a bubble have a close connection with the boundary condition. When a bubble moves near a rigid wall, it will be attracted by the Bjerknes force of the wall, and a jet pointing at the wall will be generated. In real application, the bubble may move under the combined action of walls in different directions when it forms at the corner of a pipe or at the bottom of a dam. The motion of the bubble shows complex and nonlinear characteristics under these conditions. In order to investigate the bubble pulse features near complex walls, a horizontal wall and a vertical wall are put into the experimental water tank synchronously, and an electric circuit with 200 voltages is designed to generate discharge bubbles, and then experimental study on the bubble pulse features under the combined action of horizontal and vertical walls is carried out. The influences of the combined action of two walls on the bubble shape, pulse period, moving trace and inside jet are obtained by changing the distances from bubble center to the two walls. It aims at providing references for the relevant theoretical and numerical research.     

Experimental investigation of a fast monohull in forced harmonic motions

The paper presents the results of an experimental investigation of added masses and damping coefficients of a model of a fast monohull. A model of 4.5 m length between perpendiculars was constructed of fiber glass reinforced plastic (FRP) with four segments connected by a backbone. The backbone was instrumented with load cells at the positions of the cuts. This configuration, combined with load cells measuring the force exerted by the forced motion actuators, made it possible to obtain the hydrodynamic coefficients for each of the four hull segments.

The investigation focused on the vertical motions. Thus, the experimental program included forced harmonic heave and pitch motions in calm water (no incident waves). Subtracting inertial and restoring forces from total measured forces, one obtained the hydrodynamic component, which then resulted in the hydrodynamic coefficients. The effects of steady forward speed on the radiation forces were investigated by conducting model tests at four forward speeds. Finally, nonlinear effects were assessed by conducting model tests for three amplitudes of forced heave and forced pitch motions.     

Wave enhancement due to blockage in semi-submersible and TLP structures

An exact analytical method is described to solve the diffraction problem of a group of truncated vertical cylinders. In order to account for the interaction between the cylinders, Kagemoto and Yue's exact algebraic method is utilised (Kagemoto, H., Yue, D.K.P., 1986. Interactions among multiple three-dimensional bodies in water waves: an exact algebraic method. J Fluid Mech, 46, 129–139). The isolated cylinder diffraction potential is obtained using Garret's solution and evanescent mode solutions are derived in a similar manner (Garret, C.J.R., 1971. Wave forces on a circular dock. J Fluid Mech, 46, 129–139).Free surface elevations are calculated for an array of four cylinders and compared with experiments. Comparisons show good agreement.