A nonlinear wave profile correction of the diffraction of a wave by a long breakwater: Fixed point approach

The authors of the present paper have suggested an iterative scheme to calculate the nonlinear wave profiles [Jang and Kwon, 2005. Application of nonlinear iteration scheme to the nonlinear water wave problem: Stokes wave. Ocean Engineering 32, 1862–1872]. The scheme was shown to be good for estimating nonlinear wave profiles. In the study, the iterative scheme is applied to the wave-diffraction problem by a long breakwater to calculate a diffracted wave by the breakwater. The iterative solution of diffraction was compared with the linear solution of Sommerfeld, 1896. [Mathematische Theoried der Diffraction. Mathematical Annals 47, 317–374]. For a small wave slope, the two solutions were in good agreement. However, the scheme enabled us to observe the nonlinear behaviors of a beating phenomenon and of wave profile such as Stokes’ wave for a relatively large wave slope: as the wave slope becomes larger, we can examine the nonlinear wave characteristics of the actual shapes of waves, i.e., the crests are steeper and the troughs are flatter.     

波浪的非线性显式弥散关系

波浪的非线性弥散关系在应用于求解波浪的变形问题时很不方便，需要与含非线性效应的缓坡方程一起进行迭代运算，往往导致数值计算的计算量太大，计算过于复杂。采用显式形式表达非线性弥散关系，可以克服上述缺点，大为简化波浪变形数值计算的计算量。本文通过将现有的非线性弥散关系进行分析比较，给出了一个更为一般的非线性弥散关系及其显式表达式，经比较可知，该显式弥散关系与相对应非线性弥散关系吻合的很好。本文最后用该显式结合含弱非线性效应的缓坡方程，对复式浅滩地形上的波浪折射绕射进行了计算。结果表明，考虑弱非线性可以得出与实验数据更为相符的结果，而采用显式弥散关系可以有效提高计算效率,在波浪的非线性计算中不失为一种切实有效的方法。     

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.     

Deterministic and stochastic evolution equations for fully dispersive and weakly nonlinear waves

This paper presents a new and more accurate set of deterministic evolution equations for the propagation of fully dispersive, weakly nonlinear, irregular, multidirectional waves. The equations are derived directly from the Laplace equation with leading order nonlinearity in the surface boundary conditions. It is demonstrated that previous fully dispersive formulations from the literature have used an inconsistent linear relation between the velocity potential and the surface elevation. As a consequence these formulations are accurate only in shallow water, while nonlinear transfer of energy is significantly underestimated for larger wave numbers. In the present work we correct this inconsistency. In addition to the improved deterministic formulation, we present improved stochastic evolution equations in terms of the energy spectrum and the bispectrum for multidirectional waves. The deterministic and stochastic formulations are solved numerically for the case of cross shore motion of unidirectional waves and the results are verified against laboratory data for wave propagation over submerged bars and over a plane slope. Outside the surf zone the two model predictions are generally in good agreement with the measurements, and it is found that the accuracy of e.g., the energy spectrum and of the third-order statistics is considerably improved by the new formulations, particularly outside the shallow-water range.     

Nonlinear Effect of Wave Propagation in Shallow Water

—In this paper,a nonlinear model is presented to describe wave transformation in shallow wat-er with the zero-vorticity equation of wave-number vector and energy conservation equation.Thenonlinear effect due to an empirical dispersion relation(by Hedges)is compared with that of Dalrymple＇sdispersion relation.The model is tested against the laboratory measurements for the case of a submergedelliptical shoal on a slope beach,where both refraction and diffraction are significant.The computation re-sults,compared with those obtained through linear dispersion relation.show that the nonlinear effect ofwave transformation in shallow water is important.And the empirical dispersion relation is suitable for re-searching the nonlinearity of wave in shallow water.     

An explicit finite difference model for simulating weakly nonlinear and weakly dispersive waves over slowly varying water depth

In this paper, a modified leap-frog finite difference (FD) scheme is developed to solve Non linear Shallow Water Equations (NSWE). By adjusting the FD mesh system and modifying the leap-frog algorithm, numerical dispersion is manipulated to mimic physical frequency dispersion for water wave propagation. The resulting numerical scheme is suitable for weakly nonlinear and weakly dispersive waves propagating over a slowly varying water depth. Numerical studies demonstrate that the results of the new numerical scheme agree well with those obtained by directly solving Boussinesq-type models for both long distance propagation, shoaling and re-fraction over a slowly varying bathymetry. Most importantly, the new algorithm is much more computationally efficient than existing Boussinesq-type models, making it an excellent alternative tool for simulating tsunami waves when the frequency dispersion needs to be considered.     

考虑非线性弥散影响的波浪变形数学模型

提出了逼近Kirby和Dalrymple的非线性弥散关系的显式非线性弥散关系的表达式,该显式表达式与他们的非线性弥散关系的精度几乎完全相同.采用显式非线性弥散关系,结合含弱非线性效应的缓坡方程,得到考虑非线性弥散影响的波浪变形数学模型,并对该数学模型进行了数值验证.结果表明,考虑非线性弥散影响的波浪变形数学模型更为精确.     

An extended time-dependent numerical model of the mild-slope equation with weakly nonlinear amplitude dispersion

In the present paper, by introducing the effective wave elevation, we transform the extended ellip- tic mild-slope equation with bottom friction, wave breaking and steep or rapidly varying bottom topography to the simplest time-dependent hyperbolic equation. Based on this equation and the empirical nonlinear amplitude dispersion relation proposed by Li et al. (2003), the numerical scheme is established. Error analysis by Taylor expansion method shows that the numerical stability of the present model succeeds the merits in Song et al. (2007)’s model because of the introduced dissipation terms. For the purpose of verifying its performance on wave nonlinearity, rapidly vary- ing topography and wave breaking, the present model is applied to study: (1) wave refraction and diffraction over a submerged elliptic shoal on a slope (Berkhoff et al., 1982); (2) Bragg reflection of monochromatic waves from the sinusoidal ripples (Davies and Heathershaw, 1985); (3) wave transformation near a shore attached breakwater (Watanabe and Maruyama, 1986). Comparisons of the numerical solutions with the experimental or theoretical ones or with those of other models (REF/DIF model and FUNWAVE model) show good results, which indicate that the present model is capable of giving favorably predictions of wave refraction, diffraction, reflection, shoaling, bottom friction, breaking energy dissipation and weak nonlinearity in the near shore zone.     

Linear and nonlinear irregular waves and forces in a numerical wave tank

A time-domain higher-order boundary element scheme was utilized to simulate the linear and nonlinear irregular waves and diffractions due to a structure. Upon the second-order irregular waves with four Airy wave components being fed through the inflow boundary, the fully nonlinear boundary problem was solved in a time-marching scheme. The open boundary was modeled by combining an absorbing beach and the stretching technique. The proposed numerical scheme was verified by simulating the linear regular and irregular waves. The scheme was further applied to compute the linear and nonlinear irregular wave diffraction forces acting on a vertical truncated circular cylinder. The nonlinear results were also verified by checking the accuracy of the nonlinear simulation.     

Detection of ocean waves using satellite altimetry: Application to equatorial Kelvin waves

A new method for wave motion detection from satellite altimetric measurements of sea surface height is presented. The essence of the approach is to construct a two‐dimensional traveling‐wave Fourier series representation of the amplitude field within a prespecified oceanic region. The method employs an iterative, nonlinear least‐squares technique based on the Marquardt‐Levenberg algorithm to solve for model parameters describing characteristic features of the evolving wave system. The Marquardt‐Levenberg Fourier series (MLFS) algorithm was applied to Kelvin waves active during the 1986–1987 El Nino event in the equatorial Pacific ocean using GEOSAT Exact Repeat Mission altimetry data. Characteristics of the wave system were found to be in essential agreement with earlier field measurements and the observations of Cheney and Miller (1987) obtained using time series developed from GEOSAT data. The advantage of the present detection scheme lies in its speed and ability to determine a wave system's dispersion relation over a finite range of wavenumbers, and hence the group velocity of that system.     

Comments on `Theory and application of calibration techniques foran NDBC directional wave measurements buoy' by K.E. Steele, et al.:nonlinear effects

For original paper see ibid., vol. OE-10, no.4, p.382-96 (1985). The authors of the above mentioned paper present an extensive set of linear calibration techniques that are applied to National Data Buoy Center wave-buoy sensor spectral output before calculating and disseminating directional wave spectra. The commentators identify and estimate the nonlinear effects that produce biases still present in the output, due both to wave nonlinearities themselves and to constraints on the buoy and mooring system to the driving forces. Simple models show that these nonlinearities can produce spectral energy biases of 5-15% at and above the spectral peak frequency, and even greater errors below it. NDBC presently records wave data from vertically stabilized and fixed accelerometers and slope sensors. Calculations show that these sensors all incur bias due to wave nonlinearities: this is greater for vertically stabilized accelerometers and least for slope sensors. Effects of the resulting inconsistencies between the different sensors are most pronounced below the spectral peak where the nonlinear terms dominate; these effects are illustrated with measured data     

Nonlinear wave profiles of wave–wave interaction in a finite water depth by fixed point approach

In this paper, a superposition of two periodic wave profiles in a finite water depth was investigated. This paper is focused on the improvement of a wave profile on the linear superposition of two waves. This improvement was realized by introducing an iterative method, which was based on a fixed point approach. Application of the fixed point approach to the wave superposition made it possible to obtain a wave profile of wave–wave interaction. The improved result of the wave profile was in good agreement with that of the nonlinear perturbation solution of the second order. It was interesting that the improved result revealed the higher-order nonlinear frequencies for two interacting Stokes waves while Dalzell's solution by a perturbation method could not predict them.     

Applicability of numerical models to nonlinear dispersive waves

The applicability of three different wave-propagation models in nonlinear dispersive wave fields has been investigated. The numerical models tested here are based on three different wave theories: a fully nonlinear potential theory, a Stokes second-order theory, and a Boussinesq-type theory with an improved dispersion relation. Physical experiments and computations were conducted for wave evolutions during passage over a submerged shelf under various wave conditions. As expected, the fully nonlinear solutions agree better with the measurements than do the other solutions. Although the second-order solution has sufficient accuracy for smaller-amplitude wave cases, the truncation after the third harmonics causes significant discrepancies in wave form for larger waves. In addition, the second-order model markedly overestimates the first- and second-harmonic amplitudes in transmitted waves. The Boussinesq model provides excellent predictions of wave profile over the shelf even in larger wave cases. However, this model also overestimates the magnitudes of several higher harmonics in transmitted waves. These facts may indicate that energy transfer from bound components into free waves in these higher harmonics cannot be accurately evaluated by the Boussinesq-type equations.     

Fixed-frequency Stokes wave expansion

The study describes a new fixed-frequency Stokes wave theory that differs from previous Stokes wave theories that fix the wave number. The present wave expansion analytically reveals that the wavelength increases with wave height and exceeds than the wavelength obtained by linear wave theory. A method proposed to comparably transform the wave celerity of Fenton's [Fenton, J.D., 1985. A fifth-order Stokes theory for steady waves. Journal of Waterway, Port, Coastal and Ocean Engineering 111, 216–234.] wave theory to the present one. A direct calculation of the wavelength is introduced for practical solutions, avoiding the need to solve a nonlinear equation using an iterative numerical method.     

Freely floating-body simulation by a 2D fully nonlinear numerical wave tank

Nonlinear interactions between large waves and freely floating bodies are investigated by a 2D fully nonlinear numerical wave tank (NWT). The fully nonlinear 2D NWT is developed based on the potential theory, MEL/material-node time-marching approach, and boundary element method (BEM). A robust and stable 4th-order Runge–Kutta fully updated time-integration scheme is used with regriding (every time step) and smoothing (every five steps). A special φ_n-η type numerical beach on the free surface is developed to minimize wave reflection from end-wall and wave maker. The acceleration-potential formulation and direct mode-decomposition method are used for calculating the time derivative of velocity potential. The indirect mode-decomposition method is also independently developed for cross-checking. The present fully nonlinear simulations for a 2D freely floating barge are compared with the corresponding linear results, Nojiri and Murayama’s (Trans. West-Jpn. Soc. Nav. Archit. 51 (1975)) experimental results, and Tanizawa and Minami’s (Abstract for the 6th Symposium on Nonlinear and Free-surface Flow, 1998) fully nonlinear simulation results. It is shown that the fully nonlinear results converge to the corresponding linear results as incident wave heights decrease. A noticeable discrepancy between linear and fully nonlinear simulations is observed near the resonance area, where the second and third harmonic sway forces are even bigger than the first harmonic component causing highly nonlinear features in sway time series. The surprisingly large second harmonic heave forces in short waves are also successfully reproduced. The fully updated time-marching scheme is found to be much more robust than the frozen-coefficient method in fully nonlinear simulations with floating bodies. To compare the role of free-surface and body-surface nonlinearities, the body-nonlinear-only case with linearized free-surface condition was separately developed and simulated.     

Extreme Wave Simulation with Iterative Adaptive Approach in Numerical Wave Flume

Extreme wave is highly nonlinear and may occur due to diverse reasons unexpectedly.The simulated results of extreme wave based on wave focusing,which were generated using high order spectrum method,are presented.The influences of the steepness,frequency bandwidth as well as frequency spectrum on focusing position shift were examined,showing that they can affect the wave focusing significantly.Hence,controlled accurate generation of extreme wave at a predefined position in wave flume is a difficult but important task.In this paper,an iterative adaptive approach is applied using linear dispersion theory to optimize the control signal of the wavemaker.The performance of the proposed approach is numerically investigated for a wide variety of scenarios.The results demonstrate that this approach can reproduce accurate wave focusing effectively.     

Direct time domain analysis of floating structures with linear and nonlinear mooring stiffness in a 3D numerical wave tank

In this paper, motion response of a moored floating structure interacting with a large amplitude and steep incident wave field is studied using a coupled time domain solution scheme. Solution of the hydrodynamic boundary value problem is achieved using a three-dimensional numerical wave tank (3D NWT) approach based upon a form of Mixed-Eulerian–Lagrangian (MEL) scheme. In the developed method, nonlinearity arising due to incident wave as well as nonlinear hydrostatics is completely captured while the hydrodynamic interactions of radiation and diffraction are determined at every time step based on certain simplifying approximations. Mooring lines are modelled as linear as well as nonlinear springs. The horizontal tension for each individual mooring line is obtained from the nonlinear load-excursion plot of the lines computed using catenary theory, from which the linear and nonlinear line stiffness are determined. Motions of three realistic floating structures with different mooring systems are analyzed considering various combinations of linear and approximate nonlinear hydrodynamic load computations and linear/nonlinear mooring line stiffness. Results are discussed to bring out the influence and need for consideration of nonlinearities in the hydrodynamics and hydrostatics as well as the nonlinear modelling of the line stiffness.     

A comparative linear and nonlinear ship motion study using 3-D time domain methods

This paper reports on 3D time domain seakeeping computations at different level of modeling of nonlinear effects. Four successively improved levels of computations are considered: (i) a linear computation, (ii) Froude-Krylov nonlinear computation, (iii) body nonlinear computation where the perturbation potential is computed based on the instantaneous hull under the mean free-surface, and finally (iv) the body-exact nonlinear computation, where the perturbation potential is determined based on the wetted hull under the incident wave profile after suitable mapping of the hull into a computational domain. Computations are carried out for a Wigley hull (having less ‘geometric’ nonlinearity due to vertical sides), and a S175 hull at different forward speeds. The results are obtained for both regular and irregular waves.     

Dispersion relation of internal waves in the western equatorial Pacific Ocean

Based mainly on TOGA COARE data, that is, the CI''D data from R/V Xiangyanghong No.5 (Pu et al.,1993),the temperature and current data from the Woods Hole mooring and other deep current data, the layered numerical profiles of buoyancy frequency and mean current components are figured out.A numerical method calculating internal wave dispersion relation without background shear current, used by Fliegel and Hunkins (1975),is improved to be fit for the internal wave equation with mean currents and their second derivatives.The dispersion relations and wave functions of the long crested internal wave progressing in any direction can be calculated inveniently by using the improved method.A comparison between the calculated dispersion relation in the paper and the dispersion relation in GM spectral model of ocean internal waves (Garret and Munk, 1972) is performed.It shows that the mean currents are important to the dispersion relation of internal waves in the western equatorial Pacific Ocean and that the currents make the wave progressing co-directional with (against) the currents stretched (shrink).The influence of the mean currents on dispersion relation is much stronger than that of their second derivatives, but that on wave function is less than that of their second derivatives.The influences on wave functions result in the change of vertical wavenumber, that is, making the wave function stretch or shrink.There exists obvious turning depth but no significant critical layer absorption is found.