An efficient algorithm for strapdown accelerometer-based motion measurement

The paper is concerned with the measurement of motion of free floating bodies using accelerometers strapped to a cross-arm mounted on the body. The measurement system has been studied with the help of a formulation involving a set of coupled nonlinear initial-value equations involving the angular acceleration components. A PC-based software using the Bulirsch-Stoer technique has been developed to solve the initial-value problem so as to deduce the body motions from the measured accelerations. Suitable filtering strategy has been employed at every stage of numerical integration. The reliability of the strapdown accelerometer system together with software developed has been validated using a “dry” test. Typical motion measurements have been done in all the six degrees of freedom of a tug model in a wave flume. The method is deemed to be an efficient and cost-effective technique suitable for free floating bodies and for large motions.     

Complex mathematical model of the WIG motion including the take-off mode

This paper describes an effort to develop a predictive tool for the design of a new promising marine transport — the WIG craft. The presented mathematical model of the WIG craft is capable of modeling the aerodynamics of a WIG system including the ground effect, the hydrodynamics of a stepped planing hull with a hydrofoil and also the simulation of motion for the craft. Based on extensive experience using the model, it is shown that the most important and necessary features of WIG aero- and hydrodynamics are taken into account. The results of simulations have been validated through comparison with other theoretical approaches and also with model experiments. The mathematical model is applied to investigate the dynamics of the small manually piloted WIG craft: “Hydrowing VT01”. The numerical study resulted in recommendations allowing the pilot to overcome the pitch-up tendency and also to perform the take-off manoeuvre smoothly. The stability of the WIG with a hydrofoil and also the dynamic properties of anti-collision manoeuvres have been studied and are presented.     

Study on the dynamic characteristic of a U-type tuned liquid damper

By means of Lagrange's equation, the “coupled” equations of motion for a horizontal plate carrying a U-type tuned liquid damper (TLD) are derived. The “uncoupled” equations of motion for the liquid (in the TLD) and the structural system are then obtained by decoupling the “coupled” ones. Unlike the existing literature to indirectly determine the natural frequencies of a damped vibrating system by using the resonant method, the “complex” eigenvalues of the coupled damped system are obtained directly from the associated eigenvalue equations. Besides, the pressure intensities in the two air chambers and the sizes of the two vertical tanks together with the horizontal conduit are arbitrary in the formulation of this paper. The influence of some key parameters of the TLD on the dynamic responses of the structural system is studied.     

The differences between a wing and a planing plate in two-dimensional flow

The normal force coefficient on a flat planing surface having arbitrary heave and pitch motion in two-dimensional flow is compared with the lift coefficient of a thin wing in an infinite fluid. Despite the totally different derivations, they are found to be identical (at large Froude numbers and low trim angles and allowing for the wing's interaction with twice as much fluid) at low reduced frequencies. For higher frequency motions, the wing's angle of attack induced lift and its pitch and heave damping are less than those of a planing surface, but the acceleration terms remain identical. The differences at the higher reduced frequencies are due to the fact that, in invisad irrotational flow, the planning plate cannot leave a vortex wake, whereas a wing does.It seems to follow that the “virtual mass” planing hull analysis can be applied to “quasi-static” problems involving wings and bodies in an infinite fluid without the slenderness restriction originally imposed by Jones (1946). Certainly, it is remarkable that the so called “quasi-steady” forces on a two-dimensional wing can be obtained in a few lines of elementary analysis. On the other hand, the method fails entirely when used to compute the pitching moment on a two-dimensional plate, even though it has been found to give good results for the three-dimensional case (Payne, 1981c).This work is offered as a very incomplete study of an intriguing relationship between two very different bodies of analysis. Much more work will need to be done before the relationship between the two approaches will be fully understood.     

Simulation of directional waves in a numerical basin by a desingularized integral equation approach

A numerical solution is developed to investigate the generation and propagation of small-amplitude water waves in a semi-infinite rectangular wave basin. The three-dimensional wave field is produced by the prescribed “snake-like” motion of an array of segmented wave generators located along the wall at one end of the tank. The solution technique is based on the boundary element approach and uses an appropriate three-dimensional Green function which explicitly satisfies the tank-wall boundary conditions. The Green function and its derivatives which appear in the integral equation formulation can be shown to be slowly convergent when the source and field points are in close proximity. Therefore, when computing the velocity potentials on the wave generators, the source points are chosen outside the fluid domain, thereby ensuring the rapid convergence of these functions and rendering the integral equations non-singular. Numerical results are shown which illustrate the influence of the various wavemaker and basin parameters on the generated wave field. Finally, the complete wave field produced by the diffraction of oblique waves by a vertical circular cylinder in a basin is presented.     

A new model for surface wave propagation over undulating topography

A new model is presented for the propagation of monochromatic surface waves over a region of arbitrary, one-dimensional bottom topography. The smoothly varying bed profile is divided into a series of shelves separated by abrupt steps. The wave fields on either side of each step are related by a “transfer matrix”, and the propagation of waves along the shelf between adjacent steps is described by a “rotation matrix”. Starting from a point where the surface wavefield is known, the step by step application of the appropriate combination of these matrices allows computation of the wavefield over the region of interest. If the individual steps are small then the transfer matrix reduces to a simpler plane-wave form, with considerable savings in computational effort. Comparisons are made with an exact potential solution for single and double steps in order to investigate the accuracy and validity of the matrix method. Finally, the model is applied to the case of wave reflection by fixed sinusoidal bottom undulations, and good agreement is found between its predictions and existing laboratory data.     

Estimation of the leeway drift of small craft

Small craft (<6·4 m) leeway is determined as a function of the wind speed in the range of 5–20 knots (3·6–10·3 m/sec). Leeway is calculated relative to the surface current by measurement of the separation distance of the small craft from a dyed patch of surface water at sea, using time-sequenced aerial photography. Leeway increases linearly with wind speed for small craft equipped with or without a sea anchor in the wind range studied. Leeway for small craft without sea anchor can be calculated from the equation U_L = 0.07 U_W + 0.04 where U_W is the wind speed at 2 m elevation. Leeway for small craft drifted off the be calculated from the equation U_LD = 0·05 U_W − 0·12. The small craft drifted off the downwind direction in about 80% of the experiments. The drift angle is variable and difficult to predict.     

The motion of a deep-sea remotely operated vehicle system: Part 1: Motion observations

Rapid and high-resolution motion and tension measurements were made of a caged deep-sea remotely operated vehicle (ROV) system. Simultaneous measurements were made of all six components of motion at the cage and ship A-frame and of the tension in the tether at the ship. Data were collected for cage depths of 0–1765 m. The most significant forcing was in the wave-frequency band (0.1–0.25 Hz) and accounted for over 90% of the variance of vertical acceleration. The vertical acceleration of the cage lagged the acceleration of the A-frame by up to 1.9 s and its variance was larger by up to a factor 2.2. For moderate displacements of the A-frame (≤2 m), the system is only weakly non-linear because the harmonics (3rd and 5th) of the vertical acceleration of the cage account for less than 2% of the total variance. The system is essentially one-dimensional because only the vertical motion of the cage and the vertical motion of the A-frame were coherent, while horizontal motions of the cage were weak and incoherent with any component of motion of the A-frame. The natural frequency of the system is 0.22 Hz at 1730 m, and we estimate that it is within the waveband for depths between 1450 m and the full operating depth of 5000 m.Large vertical excursions of the A-frame produce momentary slack in the tether near the cage. Retensioning results in snap loads with vertical accelerations of 0.5 gravity. Large rates of change of tension and vertical acceleration first occur at the cage during its downward motion and propagate to the surface with the characteristic speed (3870 m s⁻¹) of tensile waves for the tether. Six echoes are clearly detectable at both ends of the tether, and their pattern is extremely repeatable in different snap loads. Due to misalignment of the tether termination with the centres of mass and buoyancy, the cage pitches by up 14° during a snap. The resulting small radius of curvature poses the greatest stress on the tether.     

A coupled numerical model for simulation of wave breaking and hydraulic performances of a composite seawall

A numerical model is developed by combining a porous flow model and a two-phase flow model to simulate wave transformation in porous structure and hydraulic performances of a composite type low-crest seawall. The structure consists of a wide submerged reef, a porous terrace at the top and an impermeable rear wall. The porous flow model is based on the extended Navier-Stokes equations for wave motion in porous media and k−ε turbulence equations. The two-phase flow model combines the water domain with the air zone of finite thickness above water surface. A unique solution domain is established by satisfying kinematic boundary condition at the interface of air and water. The free surface advection of water wave is modeled by the volume of fluid method with newly developed fluid advection algorithm. Comparison of computed and measured wave properties shows reasonably good agreement. The influence of terrace width and structure porosity is investigated based on numerical results. It is concluded that there exist optimum value of terrace width and porosity that can maximize hydraulic performances. The velocity distributions inside and in front of the structure are also investigated.     

A rapid method for recovery of data from expendable surface current probes

A sighting device, called a “probe sight”, is described which simplifies the measurement, taken, from a helicopter, of the float separations from a deployed expendable surface current probe. The separation distance between the floats (along with a known time release) is a measure of the surface current speed. The sighting device can be used at any altitude. Comparisons of the readings using the probe sight and the “normal” measurements taken from aerial photographs for 180 stations in the Chukchi Sea, show that the difference between the two methods is less than 10% when the current speed is greater than 10 cm sec⁻¹. The largest difference (15%) occurs at those stations where the current speed is less than 10 cm sec⁻¹.     

Linear analysis of air cushion vehicle seakeeping response

A linearized analysis of the response of an air cushion vehicle running in waves is described. The analysis uses the linear systems approach where the vehicle is considered to be a “black box”, i.e. the response characteristics are determined experimentally from input-output relationships. The wave forces and moments are expressed in a form that produces the proper limiting behavior for infinite wavelength. Predicted motion response is shown to compare well with experimental data.     

Internal Resonances for Heave,Roll and Pitch Modes of A Spar Platform Considering Wave and Vortex-Induced Loads in the Main Roll Resonance

We present a study of the nonlinear coupling internal resonance for the heave roll and pitch performance of a spar platform under the wave and vortex-induced loads when the ratio of the frequencies of heave, roll and pitch are approximately 2:1:1. In consideration of varying wet surface, the three DOFs nonlinear coupled equations are established for the spar platform under the effect of the first-order wave loads in the heave and pitch, and vortexinduced loads in the roll. By utilizing the method of multi-scales when the vortex-induced frequency is close to the natural roll frequency, the first-order perturbation solution is obtained analytically and further validated by the numerical integration. Sensitivity analysis is performed to understand the influence of the damping and the internal detuning parameter. Two cases with internal resonance are shown. The first case is that no saturation phenomenon exists under small vortex-induced loads. The first order perturbation solution illustrates that only the vortex-induced frequency motion in roll and the super-harmonic frequency motion in heave are excited. The second case is that the vortex-induced loads are large enough to excite the pitch and a saturation phenomenon in the heave mode follows.The results show that there is no steady response occurrence for some cases. For these cases chaos occurs and large amplitudes response can be induced by the vortex-induced excitation.     

Internal Resonances for Heave, Roll and Pitch Modes of A Spar Platform Considering Wave and Vortex-Induced Loads in the Main Roll Resonance

We present a study of the nonlinear coupling internal resonance for the heave roll and pitch performance of a spar platform under the wave and vortex-induced loads when the ratio of the frequencies of heave, roll and pitch are approximately 2:1:1. In consideration of varying wet surface, the three DOFs nonlinear coupled equations are established for the spar platform under the effect of the first-order wave loads in the heave and pitch, and vortex-induced loads in the roll. By utilizing the method of multi-scales when the vortex-induced frequency is close to the natural roll frequency, the first-order perturbation solution is obtained analytically and further validated by the numerical integration. Sensitivity analysis is performed to understand the influence of the damping and the internal detuning parameter. Two cases with internal resonance are shown. The first case is that no saturation phenomenon exists under small vortex-induced loads. The first order perturbation solution illustrates that only the vortex-induced frequency motion in roll and the super-harmonic frequency motion in heave are excited. The second case is that the vortex-induced loads are large enough to excite the pitch and a saturation phenomenon in the heave mode follows. The results show that there is no steady response occurrence for some cases. For these cases chaos occurs and large amplitudes response can be induced by the vortex-induced excitation.     

A new view of near-shore dynamics based on observations from HF Radar

Since 1984 the OSCR HF Radar system has been used in over 50 deployments to measure near-shore surface currents for both scientific and engineering applications. The enhanced scope, resolution and accuracy of these measurements have yielded new insights into the tidal, wind and density driven dynamics of the near-shore zone.Tidal current ellipses obtained from these radar measurements have been shown to be in good aggrement with values calculated by numerical models both for the predominant constituents and also for higher harmonics. Coherent patterns of wind-forced currents ahve been determined with strong evidence of a “slab-like” surface response. In one deployment, with offshore winds blowing over relatively deep water, this “slab” rotated clockwise at near-inertial frequency. Strong (up to 20cm s⁻¹), persistent surface residual currents are commonly observed, these are almost certainly generated by (small) horizontal density gradients. These observed surface residuals provide ideal data for rigorous testing of 3-D numerical models.With a threatened rise in sea level, HF Radar is well-suited for observing the expected changes in the dynamics of near-shore regions. Continuing development of these radar systems offers exciting prospects of remote sensing of both surface waves and currents. Future applications may extend beyond the near-shore region to measurements along the shelf-edge, in oceanic gyres and for “beach-processes”.     

Chromium(VI) distributions in the Arctic and the Atlantic Oceans and a reassessment of the oceanic Cr cycle

Chromium(VI) concentrations ranging between 3.0 and 6.1 nmol l⁻¹ and 3.1 and 7.3 nmol l⁻¹ were found in the Arctic and Atlantic Oceans, respectively. The vertical profiles show modest depletion of chromium(VI) in surface waters, but poor overall correlations between Cr(VI) and nutrient profiles. Given that Cr(VI) is the dominant oxidation state of chromium in open-ocean waters, these data are combined with literature data to reassess the distribution of Cr in oceanic waters. It is concluded that while Cr shows some characteristics of both “recycled” and “accumulated” vertical profiles, it does not fall clearly within either group.     

Recent developments in “added-mass” planing theory

The results of several recent isolated investigations in planing theory are consolidated in this paper, together with new insights generated by a recent numerical solution of the vertically impacting wedge problem by Zhao and Faltinsen [(1992), Water entry of two-dimensional bodies. J. Fluid Mech. 246, 593–612]. As a result, in contrast to some earlier studies, it is found that the “wetted width” associated with the added mass is not that of the intersection of the wedge with the undisturbed water surface, but the wetted width of the splashed-up water, as originally proposed by Wagner [(1932), Uber Stoss-und Gleitvorgange an der Oberflache von Flussig-Keiten, Zeitschrift für Angewandte Mathematik und Mechanik, Band 12, Heft 4 (August)]. However, the splash-up ratio is not the value of (π/2–1) which he proposed, but a value which decreases with increasing deadrise, originally proposed in the late-1940s by Pierson (“Pierson's hypothesis” in the paper). For 30° deadrise, for example, Pierson's splash-up ratio is two-thirds that of Wagner's.The new equations are employed to determine the increase in the “added mass” of prismatic hull sections due to chine immersion, using experimental data. If m_o^′ is the added amss of the hull section whose chines are just wetted, Payne [(1988), Design of High-speed Boats. Volume 1: Planing. Fishergate, Inc., Annapolis, Maryland, U.S.A.] postulated that the increase in added mass due to a chine submergence (z_c) would be

where b is the chine beam and k is a constant which Payne [(1988), Design of High-speed Boats. Volume 1: Planing. Fishergate, Inc., Annapolis, Maryland, U.S.A.] gave as .The present analysis includes the “one-sided flow” correction introduced in Payne [(1990), Planing and impacting forces at large trim angels. Ocean Engng 17, 201–234]. Partly for that reason and partly because of the more precise analysis of the experimental data, the present paper revises the value to k = 2 for wetted length to beam ratios normally employed. For deadrise angles in excess of 40° and wetted keel to beam ratios in excess of 2.0, there is some evidence that k < 2.0.The revised theoretical formulation is compared with eight different sets of experimental data for flat plate and prismatic hull forms and is found to be in excellent agreement when the speed is high enough for “dynamic suction” (a loss of buoyancy at low speeds and low wetted lenghts) to be unimportant. This is true for “chines-dry” operation with deadrise angles up to 50° and chines-wet operation at length to beam ratios far in excess of the most extreme conventional practice.The research involved in performing this analysis led to the realization that different towing tanks measure different wetted chine lengths for the same hulls and test conditions. Some consistently measure more splash-up than “theory” (based on Pierson's splash-up hypothesis) predicts and others measure somewhat less than the theory. Some examples are given in Appendix B. The reason for this is not understood.     

The BCI criterion for the initiation of breaking process in Boussinesq-type equations wave models

The Breaking Celerity Index (BCI) is proposed as a new wave breaking criterion for Boussinesq-type equations wave propagation models (BTE).The BCI effectiveness in determining the breaking initiation location has been verified against data from different experimental investigations conducted with incident regular and irregular waves propagating along uniform slope [Utku, M. (1999). “The Relative Trough Froude Number. A New Criteria for Wave Breaking”. Ph.D. Dissertation, Dept. of Civil and Enviromental Engineering, Old Dominion University, Norfolk, VA; Gonsalves Veloso dos Reis, M.T.L. (1992). “Characteristics of waves in the surf zone”. MS Thesis, Department of Civil Engineering, University of Liverpool., Liverpool; Lara, J.L., Losada, I.J., and Liu, P.L.-F. (2006). “Breaking waves over a mild gravel slope: experimental and numerical analysis”. Journal of Geophysical Research, VOL 111, C11019] and barred beaches [Tomasicchio, G.R., and Sancho, F. (2002). “On wave induced undertow at a barred beach”. Proceedings of 28th International Conference on Coastal Engineering, ASCE, New York, 557–569]. The considered experiments were carried out in small-scale and large-scale facilities. In addition, one set of data has been obtained by the use of the COBRAS model based upon the Reynolds Averaged Navier Stokes (RANS) equations [Liu, P.L.-F., Lin, P., Hsu, T., Chang, K., Losada, I.J., Vidal, C., and Sakakiyama, T. (2000). “A Reynolds averaged Navier–Stokes equation model for nonlinear water wave and structure interactions”. Proceedings of Coastal Structures ‘99, Balkema, Rotterdam, 169–174; Losada, I.J., Lara, J.L., and Liu, P.L.-F. (2005). “Numerical simulation based on a RANS model of wave groups on an impermeable slope”. Proceedings of Fifth International Symposium WAVES 2005, Madrid].Numerical simulations have been performed with the 1D-FUNWAVE model [Kirby, J.T., Wei, G., Chen, Q., Kennedy, A.B., and Dalrymple, R.A. (1998). “FUNWAVE 1.0 Fully Nonlinear Boussinesq Wave Model Documentation and User's Manual”. Research Report No CACR-98-06, Center for Applied Coastal Research, University of Delaware, Newark]. With regard to the adopted experimental conditions, the breaking location has been calculated for different trigger mechanisms [Zelt, J.A. (1991). “The run-up of nonbreaking and breaking solitary waves”. Coastal Engineering, 15, 205–246; Kennedy, A.B., Chen, Q., Kirby, J.T., and Dalrymple, R.A. (2000). “Boussinesq modeling of wave transformation, breaking and run-up. I: 1D”. Journal of Waterway, Port, Coastal and Ocean Engineering, 126, 39–47; Utku, M., and Basco, D.R. (2002). “A new criteria for wave breaking based on the Relative Trough Froude Number”. Proceedings of 28th International Conference on Coastal Engineering, ASCE, New York, 258–268] including the proposed BCI.The calculations have shown that BCI gives a better agreement with the physical data with respect to the other trigger criteria, both for spilling and plunging breaking events, with a not negligible reduction of the calculation time.     

水下物体在波浪力作用下的运动计算

本文通过对水下物体在波浪力作用下所受的力的分析，运用流体力学和波浪理论的基本知识，建立了物体的运动方程，并借助于数值计算的方法对运动方程进行求解，从而得出物体在波浪力作用下的运动特性，为水下机器人的吊放回收提供了理论依据。     

孤立波与升沉水平板相互作用数值模拟研究

合理的刚度和潜深设计可以使升沉水平板获得优异的消浪性能。基于考虑流体黏性的二维不可压缩Navier-Stokes方程,以高阶紧致插值CIP(constrained interpolation profile)方法求解方程对流项,采用VOF(volume of fluid)方法重构自由液面,构建二维数值波浪水槽。采用试验数据验证模型后,研究孤立波与升沉水平板相互作用,分析相对刚度K*、相对潜深d/h、相对波高H/h对于升沉板的消浪性能和运动响应的影响,揭示升沉板对孤立波的消浪机理。研究表明:在孤立波通过时,升沉板会经历一个先上升后下降的运动,随后非线性自由振动,板下方水体近似均匀流动,且水流的垂向流动与板的垂荡方向一致;升沉板主要通过不对称涡旋脱落、浅水变形、波浪反射与辐射波转化等方式消耗孤立波能量;一定条件下,采用最优相对刚度K*=4.0和最优相对潜深d/h=0.52可以取得良好的消浪效果,此时透射系数最小,同时升沉板的运动响应在合理的范围内。