Wedge impact with the influence of ice

This work presents a canonical study on a wedge entering water near a single piece of ice using computational-fluid-dynamics (CFD) and a Wagner-type theoretical model with corrections for non-linear effects. Calculations for a series of conditions with ice of different sizes and locations relative to the wedge are conducted. The hydrodynamic force due to impact, the pressure distribution on the wedge surface, and the pile-up phenomenon are examined to study the role of ice in the impact process. The theoretical model is shown to be accurate and can serve as a useful method to assess slamming loads under the influence of ice. It is shown that even for the case of a small piece of ice, the slamming force on the wedge can increase by 30%.     

Decadal and multi-decadal variability of Labrador Sea Water in the north-western North Atlantic Ocean derived from tracer distributions: Heat budget, ventilation, and advection

Time series of profiles of potential temperature, salinity, dissolved oxygen, and planetary potential vorticity at intermediate depths in the Labrador Sea, the Irminger Sea, and the Iceland Basin have been constructed by combining the hydrographic sections crossing the sub-arctic gyre of the North Atlantic Ocean from the coast of Labrador to Europe, occupied nearly annually since 1990, and historic hydrographic data from the preceding years since 1950. The temperature data of the last 60 years mainly reflect a multi-decadal variability, with a characteristic time scale of about 50 years. With the use of a highly simplified heat budget model it was shown that this long-term temperature variability in the Labrador Sea mainly reflects the long-term variation of the net heat flux to the atmosphere. However, the analysis of the data on dissolved oxygen and planetary potential vorticity show that convective ventilation events, during which successive classes of Labrador Sea Water (LSW) are formed, occurring on decadal or shorter time scales. These convective ventilation events have performed the role of vertical mixing in the heat budget model, homogenising the properties of the intermediate layers (e.g. temperature) for significant periods of time. Both the long-term and the near-decadal temperature signals at a pressure of 1500 dbar are connected with successive deep LSW classes, emphasising the leading role of Labrador Sea convection in running the variability of the intermediate depth layers of the North Atlantic. These signals are advected to the neighbouring Irminger Sea and Iceland Basin. Advection time scales, estimated from the 60 year time series, are slightly shorter or of the same order as most earlier estimates, which were mainly based on the feature tracking of the spreading of the LSW₉₄ class formed in the period 1989-1994 in the Labrador Sea.     

On the Development of a Second-Order Bistatic Radar Cross Section of the Ocean Surface: A High-Frequency Result for a Finite Scattering Patch

The development of a model for the second-order bistatic high-frequency (HF) radar cross section on an ocean surface patch remote from the transmitter and receiver is addressed. A new approach is taken that allows a direct comparison with existing monostatic cross sections for finite regions of the ocean surface. The derivation starts with a general expression for the bistatically received second-order electric field in which the scattering surface is assumed to be of small height and slope. The source field is taken to be that of a vertically polarized dipole, and it is assumed that the ocean surface can be described, as is usually done, by a Fourier series in which the coefficients are zero-mean Gaussian random variables. Subsequently, a bistatic cross section of the surface, normalized to patch area, is derived. The result is verified by the following two means: 1) the complete form of the bistatic HF radar cross section in backscattering case is shown to contain an earlier monostatic result that has, itself, been used extensively in radio oceanography applications; and 2) the bistatic electromagnetic coupling coefficient is shown to reduce exactly to the monostatic result when backscattering geometry is imposed. The model is also depicted and discussed based on simulated data     

Impulse modelling of wave impact pressures on vertical wall

Breaking waves on coastal structures cause high magnitude impact pressures which may be important for the structural stability. In estimating the impact pressure distribution on the wall, there have been a lot of theoretical and experimental work. The present study is concerned with a theoretical approach which is based on the pressure impulse, to find the impact pressures on vertical wall. The numerical solution of the governing equation is carried out using the boundary element method. The theoretical impact pressures are determined using the experimental values of impact pressure rising time. The computational results of the impact pressures from the pressure impulse model are found to agree well with the experimental data of an earlier study.     

Predicting tsunami arrivals: Estimates and policy implications

Tsunami is one of a few kinds of natural disasters that leave people some time for escape. This escape time, which is essentially the time for the giant wave to propagate from the epicentre to a coast, has to be estimated without delay upon the occurrence of the incident. With the advancement of water wave theories, much work has been done to model the propagation of tsunamis from deep oceans to shallow water. The authors argue that while much emphasis has been put on the expansion of the high-tech early warning system and the development of complicated tsunami models, a simple-to-use yet accurate predictive model is still wanting. This paper presents a handy linear wave model, which is capable of estimating the arrival time of a tsunami with very good accuracy, as has been verified by comparison with past incidents. With the availability of such a simple model, even local communities without access to a high-tech warning system can readily estimate the time left for emergency evacuation.     

Extreme Statistics of Storm Surges in the Baltic Sea

Statistical analysis of the extreme values of the Baltic Sea level has been performed for a series of observations for 15–125 years at 13 tide gauge stations. It is shown that the empirical relation between value of extreme sea level rises or ebbs (caused by storm events) and its return period in the Baltic Sea can be well approximated by the Gumbel probability distribution. The maximum values of extreme floods/ebbs of the 100-year recurrence were observed in the Gulf of Finland and the Gulf of Riga. The two longest data series, observed in Stockholm and Vyborg over 125 years, have shown a significant deviation from the Gumbel distribution for the rarest events. Statistical analysis of the hourly sea level data series reveals some asymmetry in the variability of the Baltic Sea level. The probability of rises proved higher than that of ebbs. As for the magnitude of the 100-year recurrence surge, it considerably exceeded the magnitude of ebbs almost everywhere. This asymmetry effect can be attributed to the influence of low atmospheric pressure during storms. A statistical study of extreme values has also been applied to sea level series for Narva over the period of 1994–2000, which were simulated by the ROMS numerical model. Comparisons of the “simulated” and “observed” extreme sea level distributions show that the model reproduces quite satisfactorily extreme floods of “moderate” magnitude; however, it underestimates sea level changes for the most powerful storm surges.     

Analysis of dynamic effects relevant for the wear damage in hydraulic machines for wave energy conversion

The present paper deals with a mathematical model of a heaving-buoy Wave Energy Converter (WEC) equipped with high-pressure hydraulic power take-off machinery for energy conversion. This model is based on linear hydrodynamic theory, and a hybrid frequency-time domain model is used to study the dynamics of the heaving-body exposed to an irregular incident wave. For the power take-off system, end-stop devices are provided to protect the hydraulic machinery when the buoy is exposed to severe sea states. The model also takes into account the lubricated friction force and pressure drops of orifice flow through the valves in the hydraulic system. All the forces mentioned in the hydraulic power take-off machinery have non-linear features. A complete non-linear state space model for the WEC system is presented in this study.The WEC system was numerically simulated for different cylinder lengths under a fixed volume. The effect of fluid compressibility in the cylinder has been investigated in the mathematical model. High frequency oscillations (HFOs) caused by the compressibility of the fluid are displayed in the time series and in corresponding power spectra, and variation is shown for different cylinder sizes. Piston ring and cylinder bore wear damage is estimated by using Archard’s equation on the basis of the simulation results. A comparison of these results with a performance of an identical WEC system which neglects fluid compressibility has been done in this work. It shows that although the spectral power is small, HFO can make a large contribution to both the ring and cylinder bore wear. For the purpose of wear prediction, oscillations at or below the wave frequency and HFO may be equally important.     

Trajectory model for vertical sphere water-entry in presence of deep-seal cavity

This study presents a piecewise model for determining the vertical distance and velocity evolution with time for a sphere impacting a water surface and submerging to depths beyond deep-seal cavity pinch-off. Experimental data taken with a high-speed camera are presented for varying sphere mass ratios and impact velocities. The semi-empirical model incorporates results from previously published research and is shown to be in good agreement with experiments for heavier spheres but deviates when the sphere is only slightly denser than water. Two causes for the deviation are presented which relate to the dynamics of the cavity pinch-off event and the inception of a trailing vortex ring after the trailing cavity sloughs from the sphere. A model for predicting cavity pinch-off time and sphere position and velocity at the moment of cavity pinch-off is shown to agree well with experimental results for varying sphere mass ratios and impact velocities. The key experimental values are provided for comparison with current and future modeling efforts.     

Prediction of hydrological phenomena using self-organizing mathematical models

The paper considers the application of self-organizing models, specifically, the method of grouped arguments consideration (MGAC), to forecast short and non-stationary time series of observations in the ocean. A sequence of operations for the treatment of observational series is suggested. To assess its efficiency, we have used mean monthly oxygen concentration data collected in the surface and near-bottom layers of the Taganrog Bay. It is shown that the application of the MGAC model allows one to reduce by two times the root-mean-square error of that of the series prediction by five points, in comparison with the Jenkins-Box regressional model. It has been concluded that the predictors' non-linear functions may be effectively used in the treatment of short samplings. Translated by Vladimir A. Puchkin.     

Investigation of ship hydroelasticity

The importance of hydroelastic analysis of large and flexible container ships of today is pointed out. A methodology for investigation of this challenging phenomenon is drawn up and a mathematical model is worked out. It includes definition of ship geometry, mass distribution, structure stiffness, and combines ship hydrostatics, hydrodynamics, wave load, ship motion and vibrations. Based on the presented theory, a computer program is developed and applied for hydroelastic analysis of a flexible segmented barge for which model test results of motion and distortion in waves have been available. A correlation analysis of numerical simulation and measured response shows quite good agreement of the transfer functions for heave, pitch, roll, vertical and horizontal bending and torsion. The tool checked in such a way can be further used for reliable hydroelastic analysis of ship-like structures.     

Analytical convolution model for shipping water evolution on a fixed structure

The shipping of water is a problem that affects naval and offshore structures. Estimating its propagation on the decks of these structures by using analytical methods has been a main concern of projects. However, classical approaches disregard the decay tendency of water elevation time series and tend to overestimate the resultant water on deck. This paper is concerned with estimating the evolution of water along the deck of a fixed structure due to shipping water events. An analytical convolution model is proposed to estimate water elevations. The model considers the freeboard exceedance time series and the mean shipping flow velocity as inputs and the frictional effects of the bottom by resistance coefficients, which enables an approximated representation of the water elevation time series over the deck. It was validated with experiments of isolated shipping water events that were generated with the wet dam-break approach. The results obtained with the proposed model captured the experimental results, approximating the peak values and the decay trend of time series. Improvement of the proposed approach over classical models to represent shipping water elevations was demonstrated by comparing the results obtained with those of the dam-break model of Stoker.     

Wave interactions with multiple-row curtainwall-pile breakwaters

In this study, a mathematical model has been developed that can compute various hydrodynamic characteristics of a multiple-row curtainwall-pile breakwater. To examine the validity of the developed model, laboratory experiments have been conducted for double- and triple-row breakwaters with various combinations of drafts of curtain walls, porosities between piles, and distances between rows. Comparisons between measurement and prediction show that the mathematical model adequately reproduces most of the important features of the experimental results. As a whole, the transmission coefficient decreases with an increase in relative water depth, whereas the reflection coefficient, normalized run-up and force exhibit an opposite trend in their variations. With fixed values of the draft of the curtain wall and the porosity of lower perforated part of the first row of a double-row breakwater, as these values of the second row increase and decrease, respectively, the transmission coefficient decreases, as expected. On the other hand, their effects on wave reflection, run-up, and wave force change with the relative depth. As for the distance between the rows, the transmission coefficient becomes a maximum when it is about one half of the wave length, suggesting that this condition should be avoided to achieve the advantage of the breakwater in reducing wave transmission. It is shown that for prototype breakwaters, on an average, the transmission coefficient would be smaller than 0.3 for wave periods less than 6.0 s, and it would be about 0.45 even for the wave period of 9.0 s, although there would be a variation depending on the geometry of the breakwater. It is also shown that wave transmission is significantly reduced by multiple-row breakwaters compared with a single-row breakwater, while the difference between double-row and triple-row breakwaters is marginal. Finally, engineering monograms are provided for double-row breakwaters to be used in practical engineering applications of the breakwaters.     

Marine riser emergency disconnection analysis using scalar elements for tensioner modelling

This work proposes nonlinear parallel spring-damper scalar elements in series with a rigid beam element to represent drilling riser tensioner behavior during an emergency disconnection scenario in global riser analysis. The software package used was Orcaflex. This is a simpler approach than the more complex, full hydraulic-pneumatic-mechanic dynamic simulation modelling of the tensioning system, commonly used in literature, which can be time consuming and numerically more sensitive. The model can also be used to represent top tension variation in global riser analysis, allowing for more complex and precise simulations, especially for fatigue life and overpull estimations. The results show adequate agreement with the researched bibliography.     

A general ecosystem model for applications to primary productivity and carbon cycle studies in the global oceans

We have developed a general 1-D multi-component ecosystem model that incorporates a skillful upper ocean mixed layer model based on second moment closure of turbulence. The model is intended for eventual incorporation into coupled 3-D physical–biogeochemical ocean models with potential applications to modeling and studying primary productivity and carbon cycling in the global oceans as well as to promote the use of chlorophyll concentrations, in concert with satellite-sensed ocean color, as a diagnostic tool to delineate circulation features in numerical circulation models. The model is nitrogen-based and the design is deliberately general enough and modular to enable many of the existing ecosystem model formulations to be simulated and hence model-to-model comparisons rendered feasible. In its more general form (GEM10), the model solves for nitrate, ammonium, dissolved nitrogen, bacteria and two size categories of phytoplankton, zooplankton and detritus, in addition to solving for dissolved inorganic carbon and total alkalinity to enable estimation of the carbon dioxide flux at the air–sea interface. Dissolved oxygen is another prognostic variable enabling air–sea exchange of oxygen to be calculated. For potential applications to HNLC regions where productivity is constrained by the availability of a trace constituent such as iron, the model carries the trace constituent as an additional prognostic variable. Here we present 1-D model simulations for the Black Sea, Station PAPA and the BATS site. The Black Sea simulations assimilate seasonal monthly SST, SSS and surface chlorophyll, and the seasonal modulations compare favorably with earlier work. Station PAPA simulations for 1975–1977 with GEM5 assimilating observed SST and a plausible seasonal modulation of surface chlorophyll concentration also compare favorably with earlier work and with the limited observations on nitrate and pCO₂ available. Finally, GEM5 simulations at BATS for 1985–1997 are consistent with the available time series. The simulations suggest that while it is generally desirable to employ a comprehensive ecosystem model with a large number of components when accurate depiction of the entire ecosystem is desirable, as is the prevailing practice, a simpler formulation such as GEM5 (N²PZD model) combined with assimilation of remotely sensed SST and chlorophyll concentrations may suffice for incorporation into 3-D prediction models of primary productivity, upper ocean optical clarity and carbon cycling.     

Study on key performance parameters of hydro-pneumatic tensioner for top tensioned riser

In this paper an improved mathematical model of a hydro-pneumatic tensioner (HPT) system for top tensioned riser (TTR) is derived by consideration of friction, mass of piston and piston rod, tension loss in hydraulic oil piping, and compressibility of hydraulic oil. The vertical motion of the riser string is also considered. Subsequently, the proposed detailed model and the conventional simplified model are comparatively studied. Finally, the tension characteristic and performance parameters of the HPT are analyzed based on the proposed model. Results show that the conventional simplified model indeed overestimates the tension of tensioner as it neglects some performance parameters resulting in tension loss. The diameters of piston and piston rod and the initial pressure of high pressure (HP) gas have the most significant influence on the tension of the tensioner. The initial volume of HP gas and the initial pressure of low pressure (LP) gas also have some impact. The influence of the initial volume of LP gas and the inner diameter and length of the piping is relatively small. The research has some reference value for HPT and TTR design.     

长江河口航道拦门沙冲淤变化的数学模拟和预测

本文选用1975—1978年长江口南槽逐月水下地形图以及相应的大通站的流量和含沙量,中浚、横沙和高桥等站的潮位资料,应用数理统计理论和方法,对长江口南槽航道拦门沙的冲淤变化进行数学模拟。结果表明,多元回归数学模式能够客观地反映拦门沙变化的内在规律。且该模式的拟合计算值与实测值基本相符。因此所建模式可以预测三峡工程和上游调水后,因水沙改变而引起的拦门沙的变化趋势。     

安德拉水位计测量中水位值计算的误差分析和估算

在安德拉水位计的实际应用中,除由压力测量引起的水位测量误差外,在水位数据的计算中,也存在着误差因素。本文分析了这些因素的误差分量,并粗略地估算了它们的大小;由此提出了参与水位值计算的几个主要订正参数的选值方法。     

波浪作用下岸坡和海底平衡条件的研究

提出波浪作用下岸坡和海底动态和静态平衡条件的数学模型。在已建立的推移质泥沙体积输沙率基本关系式的基础上,根据连续方程,计算出底坡、泥沙、波浪三要素在动态和静态平衡情况下的关系式,得出反映这种关系的底坡平衡函数曲线图。用实际资料对这一函数曲线进行了验证,并对实际资料相对模型的某些差异作出解释。