Parabolic Wave Propagation Modelling in Orthogonal Coordinate Systems

This paper presents a numerical model study of the propagation of water waves using the parabolic approximation of the mild-slope equation in the orthogonal coordinate system. Two types of coordinate systems are studied: (a) a general form of orthogonal coordinate system and (b) the conformal system, a special form of orthogonal coordinate system. Two typical examples, namely, expanded breakwaters and a circular channel, are studied to validate the model. First, the examples are studied by use of the general orthogonal coordinates. Then the same examples are computed by use of the confonnal system. The computational results show that the confonnal coordinate system generally gives better predictions than the general orthogonal system. A numerical technique for generating the conformal grid is combined with the numerical model to improve the practicability of the model. The comparison between the result from the numerical grid system and that from the analytical grid system shows that reliable computational results can be obtained by use of the numerical confonnal grid system.     

Marine tether dynamics: retrieval and deployment from a heaving platform

A tether management system of a winch module in a marine environment is proposed. For the purpose of this study the subsea body is considered as a mass attached at the free end of a long tether which is wrapped around a circular drum controlled by an external torque. The winch is considered mounted on a heaving platform. The effect of the heaving platform on the motions of the drum, tether and attached mass are studied with respect to the longitudinal vibrations of the tether in one dimension. The hydrodynamic effects are considered on the deployed portion of the tether and the tethered mass, both of which are assumed submerged in otherwise still water. The resulting nonlinear system of equations of motion is developed and solved numerically for an example of a marine tether system. The effect of the tether extensibility on the operation threshold of the system is identified. Comparison with an inextensible tether case is provided. The numerical results and their analysis for the retrieval/deployment of the system are presented.     

Regeneration of sand waves after dredging

Sand waves are large bed waves on the seabed, being a few metres high and lying hundreds of metres apart. In some cases, these sand waves occur in navigation channels. If these sand waves reduce the water depth to an unacceptable level and hinder navigation, they need to be dredged. It has been observed in the Bisanseto Channel in Japan that the sand waves tend to regain their shape after dredging. In this paper, we address modelling of this regeneration of sand waves, aiming to predict this process. For this purpose, we combine a very simple, yet effective, amplitude-evolution model based on the Landau equation, with measurements in the Bisanseto Channel. The model parameters are tuned to the measured data using a genetic algorithm, a stochastic optimization routine. The results are good. The tuned model accurately reproduces the measured growth of the sand waves. The differences between the measured weave heights and the model results are smaller than the measurement noise. Furthermore, the resulting parameters are surprisingly consistent, given the large variations in the sediment characteristics, the water depth and the flow field. This approach was tested on its predictive capacity using a synthetic test case. The model was tuned based on constructed predredging data and the amplitude evolution as measured for over 2 years. After tuning, the predictions were accurate for about 10 years. Thus, it is shown that the approach could be a useful tool in the optimization of dredging strategies in case of dredging of sand waves.     

Parameter identification of a compliant nonlinear SDOF system in random ocean waves by reverse MISO method

The determination of the drag and inertia coefficients, which enter into the wave force model given by Morison's equation, is particularly uncertain and difficult when a linear spectral model is used for ocean waves, and the structure is compliant and has nonlinear dynamic response. In this paper, a nonlinear System Identification method, called Reverse Multiple Inputs–Single Output (R–MISO) is applied to identify the hydrodynamic coefficients as well as the nonlinear stiffness parameter for a compliant single-degree-of-freedom system. Four different types of problems have been identified for use in various situations and the R–MISO has been applied to all of them. One of the problems requires iterative solution strategy to identify the parameters. The method has been found to be efficient in predicting the parameters with reasonable accuracy and has the potential for use in the laboratory experiments on compliant nonlinear offshore systems.     

New Numerical Scheme for Simulation of Hyperbolic Mild-Slope Equation

The original hyperbolic mild-slope equation can effectively take into account the combined effects of wave shoaling, refraction, diffraction and reflection, but does not consider the nonlinear effect of waves, and the existing numerical schemes for it show some deficiencies. Based on the original hyperbolic mild-slope equation, a nonlinear dispersion relation is introduced in present paper to effectively take the nonlinear effect of waves into account and a new numerical scheme is proposed. The weakly nonlinear dispersion relation and the improved numerical scheme are applied to the simulation of wave transformation over an elliptic shoal. Numerical tests show that the improvement of the numerical scheme makes efficient the solution to the hyperbolic mild-slope equation. A comparison of numerical results with experimental data indicates that the results obtained by use of the new scheme are satisfactory.     

A new model for the vibration isolation via pile rows consisting of infinite number of piles

In this study, on the basis of the Floquet transform method, a numerical model for the simulation of the vibration isolation via multiple periodic pile rows with infinite number of piles is established. By means of the fictitious pile method due to Muki and Sternberg, the second kind of Fredholm integral equations for the pile rows are developed by using the fundamental solutions for the half‐space and the compatibility conditions between the piles and half‐space. Employing the Floquet transform method, integral equations for the pile rows in the wavenumber domain are then derived. Solution of the integral equations yields the wavenumber domain solution for the pile rows. The space domain solution can then be retrieved by inversion of the Floquet transform. Numerical results show that the proposed model with the Floquet transform method is in a good agreement with those of the conventional direct superposition method. On the basis of the new model, influences of the spacing between neighboring piles, the Young's modulus of the piles, and the pile length on the vibration isolation effect of the pile rows are investigated. Numerical simulations conducted in this study show that compared with the direct superposition method, the efficiency of the proposed model for simulation of the vibration isolation via pile rows is very high. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical algorithm for the dynamic analysis of base-isolated structures with dry friction

This paper deals with the techniques for solving ordinary differential equations with essential nonlinearity arising from the representation of frictional force by the sign function of the relative velocity. This problem is connected with dynamic behavior studies of base-isolated structures with dry friction devices.Mathematical evidence is presented concerning the inaccuracies of the solution which are connected with the representation of the frictional force by the sign function. An approximation of this sign model is made on the basis of the appropriate linearization of a small part of the frictional force function. A numerical algorithm for the dynamic analysis of base-isolated structures with dry friction is proposed. Numerical tests are carried out for establishing the exactness of the proposed techniques. The obtained results prove the accuracy of the techniques.     

On a power series solution to the Boussinesq equation

For certain initial and boundary conditions the Boussinesq equation, a nonlinear partial differential equation describing the flow of water in unconfined aquifers, can be reduced to a boundary value problem for a nonlinear ordinary differential equation. Using Song et al.'s (2007) [7] approach, we show that for zero head initial condition and power-law flux boundary condition at the inlet boundary, the solution in the form of power series can be obtained with Barenblatt's (1990) [2] rescaling procedure applied to the power series solution obtained in Song et al. (2007) [7] for the power-law head boundary condition. Polynomial approximations can then be obtained by taking terms from the power series. Although for a small number of terms the newly obtained approximations may be worse than polynomial approximations obtained by other techniques, any desired accuracy can be achieved by taking more terms from the power series.     

Hydraulics of runoff and loess loam deposition

On many more or less loamy soils, rill erosion is reported to start on slopes that are equal to or steeper than 2–3°; critical Froude numbers for the start of rill wash on these slopes vary between 2·0 and 3·0. This explains why colluvial deposition often occurs on slopes below 2–3°, when water spreads out at the downslope extremities of the rills. The critical hydraulic conditions for loess loam deposition were tested in the laboratory for slopes of 0·5° and 2°, applying unit-discharges (q) up to 10 cm²/s. It appeared from these experiments that for afterflow, without raindrop impact, deposition starts for critical load concentrations (c_cr) varying between several g/1 and about 60 g/l. Under rain c_cr amounts to a minimum value of 100–125 g/l and it increases when the runoff film becomes thinner. Nevertheless, deposition in pluvial runoff is also possible, as was the case during the Weichselian, according to data from quarries in Belgium and in The Netherlands. A modified Kalinske equation is proposed for c_cr prediction, with the introduction of a typical empirical coefficient C_r and considering such factors as shear stress and mean particle size. Massive sedimentation may occur when it stops raining and afterflow starts, since c_cr values are then much lower. It is shown from the Shields' diagram that loamy suspensions are more sensitive to sedimentation than sands in clear water.