Yoshimi Goda

ESTIMATION OF INCIDENT AND REFLECTED WAVES IN RANDOM WAVE EXPERIMENTS

A significant portion of the damage by hurricanes is the storm surges. The National Weather Service has developed a dynamical-numerical model to forecast hurricane storm surges. The model is used operationally for prediction, warning, and planning purposes. The model requires fixed oceanographic and real time meteorological input data. The oceanographic data were prepared for the Gulf and East coasts of the U.S. and are stored as an essential part of the program. Meteorological data for any tropical storm are supplied by the forecasters or planners using the model. The model was applied to hurricane Camille 1969. Comparison between the observed and computed surges for Camille was satisfactory for prediction purposes.The main consideration in harbor master planning is to maximize the amount of time that the harbor can be used. The potential level of harbor utilization can be evaluated by analyzing vessel performance during harbor operations in terms of the range of imposed environmental conditions. The harbor utilization level is expressed statistically as the probable amount of time that the harbor can be used as planned.The design and construction of a major ocean outfall and diffuser system for disposal of wastewater effluents is a complex process involving an interplay of requirements originating from various disciplines. These include, among others, considerations of physical oceanography, mixing and dispersion, treatment processes, regulatory requirements, marine geology, economics and construction. The recently completed Sand Island Outfall and the newly designed Barbers Point Outfall are both on the southern coast of the island of Oahu, Hawaii, and are designed for treated sewage effluents from the densely populated portion of the City and County of Honolulu. In this paper, some design considerations of these outfalls will be examined. The emphasis in this paper is on the hydrodynamics, although other design aspects are also discussed briefly.The height of dikes and other coastal structures can only be calculated after determination of the wave run-up. Several formulas for the calculation of wave run-up are developed after model tests as a rule. But the influences of scale effects and natural wind conditions are practically unknown. To clear these questions further investigations and especially field measurements should be carried out. By measuring the markerline of floating trash on the slope of the seadikes the maximum wave run-up could be found out after four storm surges in 1967 and 1973In two graphs it will be shown that on the tidal flats the run-up depends on the waterdepth. The run-up was higher than it could be expected after model tests of 1954. With a newly developed special echo sounder the run-up could be measured in January 1976. The waves and the run-up could be registrated synchronously during two severe storm surges. As shown in Fig. 9 it was found a logarithmic distribution of the wave height, wave period and the higher part of the wave run-up. The found wave run-up is considerably higher than estimated before. The measured 98 % run-up is found about twice the computed value. That is an interesting and important result of the first synchronous recording of wave run-up on sea dikes.In March 1972 the authors firm in association with two Portuguese firms of consulting engineers, Consulmar and Lusotecna, were appointed by the Portuguese Government agency Gabinete da Area de Sines to prepare designs for the construction of a new harbour at Sines on the west coast of Portugal. The location is shown in Figure 1. The main breakwater, which is the subject of this paper, is probably the largest breakwater yet built, being 2 km long and in depths of water of up to 50 m. It is exposed to the North Atlantic and has been designed for a significant wave height of 11 m. Dolos units invented by Merrifield (ref. 1) form the main armour. The project programme required that studies be first made of a wide range of alternative layouts for the harbour. After the client had decided on the layout to be adopted, documents were to be prepared to enable tenders for construction to be invited in January 1973. This allowed little time for the design to be developed and only one series of flume tests, using regular waves, was completed during this period. Further tests in the regular flume were completed during the tender period and a thorough programme of testing with irregular waves was commenced later in the year, continuing until August 1974 when the root of the breakwater was complete and the construction of the main cross-section was about to start. The model tests, which were carried out at the Laboratorio Nacional de Engenharia Civil in Lisbon, were reported by Morals in a paper presented to the 14th International Coastal Engineering Conference in 1974. (ref. 2)Estuaries may be sequentially classified into highly stratified, moderately mixed and vertically homogeneous. An important difference between moderately stratified or vertically homogeneous estuaries, and highly stratified estuaries (salt wedges) is that, in the former, tidal currents are sufficient to cause turbulent mixing of fresh water and sea water over the full depth of the estuary. In the latter, a distinct interface or interfacial layer exists which separates the two nearly homogeneous layers. The vertical advectlon of salt in this two-layer flow is the dominant process in maintaining the salt balance. This paper presents an analytical model describing this process. Experiments have been conducted in the laboratory to compare with the developed theory.

Irregular Wave Deformation in the Surf Zone

A theoretical model of irregular wave breaking and dissipation is developed. The setup of irregular waves is computed by means of radiation stress in association with the above wave breaking model....

NEW WAVE PRESSURE FORMULAE FOR COMPOSITE BREAKWATERS

The steady state profile of the longshore current induced by regular, obliquely incident, breaking waves, over a bottom with arbitrary parallel bottom contours, is predicted. A momentum approach is adopted. The wave parameters must be given at a depth outside the surf zone, where the current velocity is very small. The variation of the bottom roughness along the given bottom profile must be prescribed in advance. Depth refraction is included also in the calculation of wave set-down and set-up. Current refraction and rip-currents are excluded. The model includes two new expressions, one for the calculation of the turbulent lateral mixing, and one for the turbulent bottom friction. The term for the bottom friction is non-linear. Rapid convergent numerical algorithms are described for the solution of the governing equations. The predicted current profiles are compared with laboratory experiments and field measurements. For a plane sloping bottom, the influence of different eddy viscosities and constant values of bottom roughness is examined.The calculation of turbulent flow using Naviers equations assumes the introduction of a turbulent viscosity coefficient the value of which is normally constant, conforming with Boussinesqs hypothesis. It was shown that setting aside this hypothesis, a velocity profile quite different to that resulting from the classic theory is obtained in the case of flow induced by wind. This result appears to be confirmed by the tests carried out in the Mediterranean. The advantage of this method is that it gives the vertical turbulent diffusion which is of particular interest to pollution studies.In the numerical method of prediction of wind waves in deep water, Hasselmanns nonlinear interaction theory is applied. This method assumes the energy balance of individual component waves. However, the total energy balance must exist in the transformation of irregular waves in shoaling water. In this investigation, experiments were carried out on the transformations in shoaling water of composite waves having two components and random waves having one or two main peaks. It was found that the elementary component wave height of the composite waves and the elementary peak power of the random waves decrease with decrease in the water depth. This reason can be explained qualitatively by the theory of the elementary component wave height change of finite amplitude waves in shoaling water. The secondary component wave height of the composite waves and the secondary peak power of the random waves increase with decrease in the water depth. This can be explained qualitatively by Hamadas theory of nonlinear interaction in uniform depth.Experiments have been carried out by using non-breaking waves and breaking waves to investigate the wave forces on a vertical circular cell located in the shallow water. Based on the experimental data, the drag coefficient and the inertia coefficient of a circular cylinder and the curling factor of breaking waves are estimated, and the computation methods of wave forces are examined. As a result, it is shown that the phase lag of inertia forces behind the accelerations of water particles should be considered for the estimation of the drag coefficient as well as the inertia coefficient. In addition the previous formula of the maximum breaking wave forces acting on a cell or a pile is revised by introducing the effects of the above-mentioned phase lag and another phase difference, both of which are functions of the ratio of the cell diameter to the wave length. • It is confirmed that the proposed formula is applicable even to the large cell with the diameter comparable to the wave length. INTRODUCTION Many studies have been done on the impulsive pressures acting on a vertical wall, but there has been very little investigation of breaking wave forces on a cell-type structure. The breaking wave forces should be taken into consideration all the same in the design of pile-type or cell-type structures in nearshore area, because breaking waves cause extreme shock pressures on a cell structure asThe air bubble plume induced by the steady release of air into water has been analyzed with an integral technique based on the equations for conservation of mass, momentum and buoyancy. This approach has been widely used to study the behavior of submerged turbulent jets and plumes. The case of air-bubble induced flow, however, includes additional features. In this study the compressibility of the air and the differential velocity between the rising air bubbles ,and the water are introduced as basic propertie s of the air bubble plume in addition to a fundamental coefficient of entrainment and a turbulent Schmidt number characterizing the lateral spreading of the air bubbles. Theoretical solutions for twoand three-dimensional air-bubble systems in homogeneous, stagnant water are presented in both dimensional and normalized form and compared to existing experimental data. The further complication of a stratified environment is briefly discussed since this case is of great practical interest. This paper is to be considered as a progress report, as future experimental verification of various hypotheses is needed.

Recommended practice for extreme wave analysis

Recommended methods for the statistical analysis of extreme waves are presented. Proper data selection is stressed. Data from different seasons and/or of different type or origin should be analysed separately. Use of the peaks over threshold method is recommended. It is advised that a three parameter Weibull distribution is fitted to the observed storm peak data giving due consideration to data censoring.

A COMPARATIVE REVIEW ON THE FUNCTIONAL FORMS OF DIRECTIONAL WAVE SPECTRUM

Standard functional forms of directional wave spectrum are presented and compared. Frequency spectra of the JONSWAP, Wallops, and TMA are presented as the functions of significant wave height and s...

ON THE METHODOLOGY OF SELECTING DESIGN WAVE HEIGHT

The Labrador Sea Extreme Waves Experiment (LEWEX), is an international basic research programme concerned with full-scale measurements, analysis, modelling and simulation in test basins of 3-dimensional seas. The research is carried out in order to assess the significance of 3-dimensional sea states in engineering applications. The first phase of the programme full scale wave measurements in the North Atlantic Ocean was performed at a site and time that had a high probability of encountering severe sea states. The present publication shows examples of measured bi-modal directional sea spectra obtained with the WAVESCAN buoy and directional sea spectra measured with an airborne Synthetic Aperture Radar (SAR). Directional spectra of gravity waves are obtained with the SAR both in open waters and below an ice cover. Further work is needed in order to verify SAR-measurements with in-situ observations. In-situ measured directional spectra are also compared with hindcast spectra from the 3G-WAM model. Hindcast significant wave heights were found to be lower than the in-situ measurements.The littoral drift model developed at DHI and ISVA, see Deigaard et al. (1986b) has been extended to include the effects of the irregularity of the waves, of a coastal current and a wind acting on the surf zone. Further, a mathematical model to simulate the near-shore current pattern along a barred coast with rip channels has been developed. The influence on the littoral drift of the irregularity of waves, wind, coastal current, and rip channels is discussed. It is concluded that irregularity of waves and presence of rip channels must be considered while coastal current and wind action are of minor importance.At present, the Prodhoe Bay oil field in Alaska contributes a substantial amount of the domestic oil production of the United States. Oil is also expected to be present on the continental shelf of Alaska, and it is estimated that approximatedly 28 percent of the total U. S. reserve is located beneath the shallow ice covered seas of the Alaskan continental shelf. To expolre and to exploit these oil rich resources, engineers are confronted with hostile oceanographic conditions such as high tides, waves, strong currents and sea ice. The same area is also rich in fishery resources. Being one of the most productive fishing grounds in the nothern Pacific, the potential ecological impact due to an oil spill is of a major concern. This paper describes the methologies used for the development of a modeling system for the oil risk analysis. The system is designed with generality in mind so it can be used for other coastal areas. The development of three dimensional models used in the modeling system described here have been published in the earlier International Coastal Engineering Conferences (Liu and Leendertse, 1982, 1984, 1986) and a report published recently by RAND (Liu and Leendertse, 1987). In the oil-spill risk analysis, these three dimensional hydrodynamic models are coupled to a two-dimensional stochastic weather model and an oil weathering model.The two projects (LUBIATOWO 79 and LUBIATOWO 86) were aimed at study pore pressure behaviour in natural sand bed in the coastal zone of the Baltic sea under the action of storm waves. During both projects, the wave induced pore pressures at the various levels in the sea bed were measured. The collected data were used to verify the applicability of the various theoretical approaches. In the conclusion, the range of the application of the particular analytical method is given.The rapid recession of the shingle bank of Hurst Beach (up to 3.5m/yr) makes it an excellent natural laboratory for the study of the factors which influence the stability of shingle beaches. Studies have included: the significance of long period, high energy, swell waves the classification and quantification of overwash processes run-up and seepage characteristics the effect of settlement of the underlying strata and the implications for practices in shingle nourishment. The studies have revealed the distinctive character of shingle beaches as compared with the more fully researched sand beaches. More detailed research on shingle beaches is justified particularly in relation to (i) the run-up characteristics including its interaction with swash cusps and (ii) the influence of the subsidiary sand fraction on the beach characteristics.Environmental assessment, engineering studies and designs were completed for a new 26.5 m3/s seawater intake system in the Persian Gulf. The original intake facility consisted of a curved, 60m breakwater with one end attached to the shoreline, a settling basin immediately adjacent to the shoreline and dredged to a maximum depth of approximately 5m, and a pumphouse structure located on shore such that the seaward wall formed one side of the settling basin. The facility located on an island in the Gulf, which served multiple seawater uses, had experienced both structural and operational problems, the latter consisting principally of excessive ingestion of sediment and seaweed. These factors plus the requirement for additional demands for seawater beyond plant capacity caused the owner to initiate a study of alternative intake systems, produce a design for the most effective solution and construct the new intake system.A two-dimensional wave prediction model suitable for use on personal computers is described. The model requires the twodimensional time-dependent wind field as input. Output consists of wave height, wave period, and wave direction estimates at all grid points on a computational grid representing an enclosed or semi-closed basin. Model predictions compare favorably with observations from a wave research tower in Lake Erie. A formula is provided to estimate how long a model simulation would take on a personal computer given the surface area of the computational domain, the grid size, and the computer clock speed.Cullera Bay is a neritic ecosystem placed on the Spanish Mediterranean Littoral largely influenced by the Jucar River, that brings about lower salinities than surrounding waters, and broad variations of its values. An extensive research, with 9 samplings throughout the year, was carried out, measuring both physical and chemical parameters, and the planktonic communities. The trophic status of the ecosystem, the spatial and temporal variations of the nutrients and the planktonic communities were studied, evaluating the influence of the river loads and the littoral dynamics. Some essential basis to allow a suitable emplacement of waste waters disposals along the Valencian littoral are set up in order to minimize the gradual eutrophication of this coast.In the last two years a whole of studies was realized in order to determine precise solutions to the regeneration of Villajoyosas beach, in the Spanish mediterranean coast. Investigations were carried out to the surrounding coastal areas based in field investigations and laboratory analyses of the beaches materials.

A RELIABILITY DESIGN METHOD OF CAISSON BREAKWATERS WITH OPTIMAL WAVE HEIGHTS

The paper presents an improved version of the deformation-based reliability design method by means of expected sliding distance proposed by Shimosako and Takahashi (1998 and 1999). The process of economic optimization is added, and the concept of the optimal return period for selection of design wave heights is introduced with a sample diagram for the optimal return period for three values of service lifetime. It proposes a new tolerable limit of expected sliding distance of 0.10 m instead of the value of 0.3 m adopted by Shimosako and Takahashi. The ratio of water depth to the significant wave height is demonstrated to govern the reliability design of caisson breakwaters, because of random wave breaking effects. In shallow water, a breakwater can be designed with the cross section smaller than that by the conventional design, but a larger cross section is needed in deeper water. Discussions are made on the comparison between the partial safety factor systems proposed by European groups and on the future problems to be solved for further refinements of the deformation-based design method.

REANALYSIS OF REGULAR AND RANDOM BREAKING WAVE STATISTICS

Statistics of breaking waves across the surf zone are reanalyzed on the basis of various sets of field and laboratory data so as to provide coastal engineers with reliable information on breaking waves. The breaker index or the ratio of wave height to water depth is to be expressed as a function of the two parameters of beach slope and relative depth, and Godas breaker index formula is slightly to reduce the slope effect. The breaker index for regular waves has inherent variability as expressed with the coefficient of variability, which increases from 6% to 14% as the beach slope becomes steep up to 1/10. The incipient breaking height of the significant wave is about 30% lower than that of regular waves, but the ratio of significant wave height to water depth gradually increases within the surf zone toward the shoreline. The wave height distribution is the narrowest in the middle of the surf zone, but it returns to the Rayleigh distribution near the shoreline owing to wave regeneration after breaking. The nonlinearity of random waves is strongest at the outer edge of the surf zone, but it is destroyed by the wave breaking process inside the surf zone. The ratios of statistical wave heights H 1/10, H 1/3 and H rms to the spectral significant wave height H m0 are shown to increase as the wave nonlinearity parameter increases up to the outer edge of the surf zone.

DIFFRACTION DIAGRAMS FOR DIRECTIONAL RANDOM WAVES

The equations describing conservation of mass, momentum and energy in a turbulent free surface flow are derived for a controle volume extending over the whole depth. The effect of the turbulent surface oscillations are discussed but neglected in the following analysis, where the equations are applied to the energy balance in a surf zone wave motion. This leads to results for the wave height variation and the velocity of propagation. The results cannot be reconciled completely with measurements and the concluding discussion is aimed at revealing how the model can be improved.A three-dimensional morphodynamic model of sequential beach changes Is presented. The model Is based on variations in breaker wave power generating a predictable sequence of beach conditions. The spectrum of beach conditions from fully eroded-dissipatlve to fully accreted reflective is characterised by ten beach-stages. Using the breaker wave power to beach-stage relationship the model Is applied to explain temporal, spatial and global variations In beach morphodynamlcs.The agents of initial damage to the dunes are water, which undermines them, and animals (including man) which damage the protective vegetation by grazing or trampling. Of these, man has recently assumed predominant local importance because of the popularity of sea-side holidays and of the land-falls of certain marine engineering works such as oil and gas pipelines and sewage outfalls. The need is therefore increasing for active dune management programmes to ensure that under these accentuated pressures, the coast retain an equilibrium comparable with that delicately balanced equilibrium which obtains naturally at a particular location.

Dynamic response of upright breakwaters to impulsive breaking wave forces

Abstract A formulation is made for the motion of an upright section of a composite breakwater activated by an impulsive force due to breaking waves. The foundation, including a rubble mound, is represented by a system of mass and dual springs for rotational and horizontal motions. The validity of the model was confirmed by a small scale model test in which a concrete block rested upon a crushed stone mound was hit by a pendulum of known momentum. The model is then combined with the momentum theory of impulsive breaking wave force to yield an estimate of the sliding distance of an upright breakwater. The coefficients of elastic uniform shear and nonuniform compression of the foundation are estimated to be a few hundred tons per cubic meters on the basis of the vibration data of the Alger and Haboro breakwaters. The equivalent pressure of breaking waves effective for breakwater sliding is then estimated to be 3 times or less the hydrostatic head of the wave height.