Fully nonlinear numerical wave tank (NWT) simulations and wave run-up prediction around 3-D structures

A finite-difference scheme and a modified marker-and-cell (MAC) algorithm have been developed to investigate the interactions of fully nonlinear waves with two- or three-dimensional structures of arbitrary shape. The Navier–Stokes (NS) and continuity equations are solved in the computational domain and the boundary values are updated at each time step by the finite-difference time-marching scheme in the framework of a rectangular coordinate system. The fully nonlinear kinematic free-surface condition is implemented by the marker-density function (MDF) technique developed for two fluid layers.To demonstrate the capability and accuracy of the present method, the numerical simulation of backstep flows with free-surface, and the numerical tests of the MDF technique with limit functions are conducted. The 3D program was then applied to nonlinear wave interactions with conical gravity platforms of circular and octagonal cross-sections. The numerical prediction of maximum wave run-up on arctic structures is compared with the prediction of the Shore Protection Manual (SPM) method and those of linear and second-order diffraction analyses based on potential theory and boundary element method (BEM). Through this comparison, the effects of non-linearity and viscosity on wave loading and run-up are discussed.     

The dimensions of sand ripples in full-scale oscillatory flows

New large-scale experiments have been carried out in two oscillatory flow tunnels to study ripple regime sand suspension and net sand transport processes in full-scale oscillatory flows. The paper focuses on ripple dimensions and the new data are combined with existing data to make a large dataset of ripple heights and lengths for flows with field-scale amplitudes and periods. A feature of the new experiments is a focus on the effect of flow irregularity. The combined dataset is analysed to examine the range of hydraulic conditions under which oscillatory flow ripples occur, to examine the effects of flow irregularity and ripple three-dimensionality on ripple dimensions and to test and improve existing methods for predicting ripple dimensions.The following are the main conclusions. (1) The highest velocities in a flow time-series play an important role in determining the type of bedform occurring in oscillatory flow. Bedform regime is well characterised by mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth peak velocities in the case of irregular flow. (2) For field-scale flows, sand size is the primary factor determining whether equilibrium ripples will be 2D or 3D. 2D ripples occur when the sand D₅₀ ≥ 0.30 mm and 3D ripples occur when D₅₀ ≤ 0.22 mm (except when the flow orbital diameter is low). (3) Ripple type (2D or 3D) is the same for regular and irregular flows and ripple dimensions produced by equivalent regular and irregular flows follow a similar functional dependence on mobility number, with mobility number based on maximum velocity in the case of regular flow and based on the mean of the highest one tenth velocities in the case of irregular flow. For much of the ripple regime, ripple dimensions have weak dependency on mobility number and ripple dimensions are similar for regular and irregular flows with the same flow orbital amplitude. However, differences in ripples produced by equivalent regular and irregular flows become significant at the high mobility end of the ripple regime. (4) Ripple dimensions predicted using the Wiberg and Harris formulae are in poor agreement with measured ripple dimensions from the large-scale experiments. Predictions based on the Mogridge et al. and the Nielsen formulae show better overall agreement with the data but also show systematic differences in cases of 3D ripples and ripples generated by irregular flows. (5) Based on the combined large-scale data, modifications to the Nielsen ripple dimension equations are proposed for the heights and lengths of 2D ripples. The same equations apply to regular and irregular flows, but with mobility number appropriately defined. 3D ripples are generally smaller than 2D ripples and estimates of 3D ripple height and length may be obtained by applying multipliers of 0.55 and 0.73 respectively to the 2D formulae. The proposed modified Nielsen formulae provide an improved fit to the large-scale data, accounting for flow irregularity and ripple three-dimensionality.     

Examination of empirical formulas for wave shoaling and breaking on steep slopes

Estimation of the wave height transformation of shoaling and breaking is essential for the nearshore hydrodynamics and the design of coastal structures. Many empirical formulas have been well recognized to the wave height transformation, but most of them were only applicable for gentle slopes. This paper reports the experimental results of wave shoaling and breaking over the steep slopes to examine the applicability of the previous empirical formulas. Two steep bottom slopes of 1/3 and 1/5, and one gentle slope of 1/10 were conducted in the present experiments. Experimental results show that the shoaling distance of steep slopes become short and the surface waves may be partially reflected from the steep bottom, thus the estimation of wave shoaling using the well-known previous formula did not conform completely to the experimental results. The previous empirical formulas for the wave breaking criteria were also examined, and the modified equations to the steep beaches were proposed in this work. A numerical model was finally adopted to calculate the wave height transformation in the surf zone by introducing the modified breaking index.     

The nonlinear hydrodynamic model for simulating a ship steering in waves with autopilot system

In the paper, a hydrodynamic numerical model including wave effect is developed to simulate ship autopilot systems by using the time domain analysis. The PD controller and the sliding mode controller are adopted as the autopilot systems. The differences of simulation results between two controllers are analyzed by cost function composed of heading angle error and rudder deflection, either in calm water or in waves. The results in calm water show that both controllers are tracking well for the desired route with the similar cost function value by tuning the key design parameters. However, the course tracking ability of the controller using sliding mode in waves is generally better even the cost function value is similar.     

Analysis of seabed instability using element free Galerkin method

Wave-induced seabed instability, either momentary liquefaction or shear failure, is an important topic in ocean and coastal engineering. Many factors, such as seabed properties and wave parameters, affect the seabed instability. A non-dimensional parameter is proposed in this paper to evaluate the occurrence of momentary liquefaction. This parameter includes the properties of the soil and the wave. The determination of the wave-induced liquefaction depth is also suggested based on this non-dimensional parameter. As an example, a two-dimensional seabed with finite thickness is numerically treated with the EFGM meshless method developed early for wave-induced seabed responses. Parametric study is carried out to investigate the effect of wavelength, compressibility of pore fluid, permeability and stiffness of porous media, and variable stiffness with depth on the seabed response with three criteria for liquefaction. It is found that this non-dimensional parameter is a good index for identifying the momentary liquefaction qualitatively, and the criterion of liquefaction with seepage force can be used to predict the deepest liquefaction depth.     

On the fairing of resistance test data: A software-based approach

A fairing method is presented for giving a good representation of ship model resistance test data, which in general tend to exhibit a scattered trend and need to be smoothed appropriately. One component of the test data originates from the frictional resistance and can be formulated smoothly. The second component is the wave resistance which exhibits a specific harmonic character with peaks and troughs. In this method, this component is modelled using suitable harmonic cosine terms for obtaining the best representation of the corresponding resistance data. The performance of the method is then tested on two example cases, using the characteristic parameters of the sample data as input during the smoothing process.     

Diet restriction induced autophagy: a lysosomal protective system against oxidative- and pollutant-stress and cell injury

Nutrient deprivation or dietary restriction (DR) confers protection against ageing and stress in many animals and induced lysosomal autophagy is part of this mechanism. The effects of dietary restriction on the toxicity of copper and the polycyclic aromatic hydrocarbon phenanthrene have been investigated in the common marine mussel Mytilus edulis. The findings show that DR-induced autophagy facilitates the recovery of the digestive gland (i.e., molluscan liver analogue) from cell injury caused by both copper and phenanthrene. It is inferred that DR-induced autophagy and lysosomal proteolysis results in improved cellular "housekeeping" through the more efficient removal of oxidatively and pollutant damaged proteins (e.g., protein carbonyls, protein adducts, etc.) and that this contributes to stress resistance.     

海浪周期关系的研究

本文利用现场观测数据、实验室风浪槽观测数据和由文氏谱及JONSWAP谱生成的模拟数据,研究了有效周期与谱平均周期、谱峰周期与谱平均周期以及有效周期与谱峰周期的关系,通过数据拟合给出了相应的关系式.研究发现,有效周期与由谱的负阶矩计算的平均周期之间的关系更加稳定,并且有效周期与负2阶矩计算的平均周期几乎相等,均代表了海浪主导波对海浪平均周期的贡献.     

青岛近海地区SO2和SO42-干沉降速率的季节变化

用小麦岛和团岛观测站提供的青岛近海气象和水温资料,以适于SO2气体和硫酸盐(SO42-)气溶胶的干沉降模型研究了模型中的输送阻尼对两种污染物质干沉降速率的贡献,并计算了2003年青岛近海地区两种污染物质干沉降速率的季节变化。结果表明:对SO2而言,粘性副层传输系数对沉降速率的贡献大于空气动力学传输系数;对于SO42-,空气动力学传输系数对沉降速率的贡献要大于表面传输系数。SO2和SO42-的干沉降速率变化范围分别为0.187~0.868cm/s和0.188~0.532cm/s。两种污染物质干沉降速率的四季变化有相似的规律,即冬季>秋季>春季>夏季。     

Wave reflection from partially perforated-wall caisson breakwater

In 1995, Suh and Park developed a numerical model that computes the reflection of regular waves from a fully perforated-wall caisson breakwater. This paper describes how to apply this model to a partially perforated-wall caisson and irregular waves. To examine the performance of the model, existing experimental data are used for regular waves, while a laboratory experiment is conducted in this study for irregular waves. The numerical model based on a linear wave theory tends to over-predict the reflection coefficient of regular waves as the wave nonlinearity increases, but such an over-prediction is not observed in the case of irregular waves. For both regular and irregular waves, the numerical model slightly over- and under-predicts the reflection coefficients at larger and smaller values, respectively, because the model neglects the evanescent waves near the breakwater.