Water Of Infinite Depth Research Articles

Linear stability of finite-amplitude capillary waves on water of infinite depth

AbstractWe study the linear stability of the exact deep-water capillary wave solution of Crapper (J. Fluid Mech., vol. 2, 1957, pp. 532–540) subject to two-dimensional perturbations (both subharmonic and superharmonic). By linearizing a set of exact one-dimensional non-local evolution equations, a stability analysis is performed with the aid of Floquet theory. To validate our results, the exact evolution equations are integrated numerically in time and the numerical solutions are compared with the time evolution of linear normal modes. For superharmonic perturbations, contrary to Hogan (J. Fluid Mech., vol. 190, 1988, pp. 165–177), who detected two bubbles of instability for intermediate amplitudes, our results indicate that Crapper’s capillary waves are linearly stable to superharmonic disturbances for all wave amplitudes. For subharmonic perturbations, it is found that Crapper’s capillary waves are unstable, and our results generalize to the highly nonlinear regime the analysis for small amplitudes presented by Chen & Saffman (Stud. Appl. Maths, vol. 72, 1985, pp. 125–147).

A Rigorous Justification of the Modulation Approximation to the 2D Full Water Wave Problem

We consider the 2D inviscid incompressible irrotational infinite depth water wave problem neglecting surface tension. Given wave packet initial data, we show that the modulation of the solution is a profile traveling at group velocity and governed by a focusing cubic nonlinear Schrodinger equation, with rigorous error estimates in Sobolev spaces. As a consequence, we establish existence of solutions of the water wave problem in Sobolev spaces for times in the NLS regime provided the initial data is suitably close to a wave packet of sufficiently small amplitude in Sobolev spaces.

Hydrodynamic exciting forces on a submerged oblate spheroid in regular waves

It is the purpose of this study to provide the analytic solution for the hydrodynamic diffraction problem by stationary, submerged oblate spheroidal bodies subjected to harmonic incident waves in deep water. The analytical process employs the multipole expansion terms derived by Thorne [1] which describe the velocity potential at singular points within a fluid domain with free upper surface and infinite water depth. The multipole potentials are used to analytically formulate the diffraction component of the velocity potential which is initially described by relations involving both spherical and polar coordinates. The goal is to transform the constituent terms of the multipole potentials as well as the incident wave component in oblate spheroidal coordinates. To this end, the appropriate addition theorems are derived which recast Thorne’s [1] formulas into infinite series of associated Legendre functions.

A numerical study of the second-order wave excitation of ship springing in infinite water depth

Assuming the wave steepness of the incident waves and the ship motions are small, the second-order weakly-nonlinear hydrodynamic problem of a ship moving with constant forward speed is studied numerically in a consistent way. The boundary value problem is formulated in a body-fixed coordinate system and the perturbation scheme is used. This formulation does not include any derivatives of the velocity potential on the right-hand side of the body-boundary conditions, and thus avoid the difficulties associated with the terms similar to the so-called mj-terms and their derivatives. The second-order sum-frequency wave excitation of ship springing is studied in both monochromatic and bichromatic head-sea waves. Different Froude numbers are considered. A time-domain Higher-Order Boundary Element Method based on cubic shape function is used as a numerical tool. An upstream finite difference scheme is used for longitudinal derivative terms in the free-surface conditions. For a modified Wigley hull in head-sea waves, it is found that the second-order velocity potential gives dominant contribution to second-order wave excitation of ship springing in the wave frequency region where sum-frequency springing occurs. Quadratic velocity terms in the Bernoulli equation have a relatively small contribution. The numerical results also demonstrate strong dependency of the second-order wave excitation of ship springing on the Froude numbers for small wave lengths. The effect of beam and draft is investigated.

A precorrected-FFT higher-order boundary element method for wave–body problems

The boundary element method (BEM) has been widely applied in the field of wave interaction with offshore structures, but it is still not easy to use in resolving large-scale problems because of computing costs and computer storage being increased by O( N 2) for the traditional BEM. In this paper a precorrected Fast Fourier Transform (pFFT) higher-order boundary element method (HOBEM) is proposed for reducing the computational time and computer memory by O( N). By using a free-surface Green function for infinite water-depth, the disadvantage of the Fast Multipole Boundary Element Method (FMBEM)—i.e. unable to solve infinite deep-water wave problems—can be overcome. Numerical results from the problems of wave interaction with single- and multi-bodies show that the present method evidently has more advantages in saving memory and computing time, especially for large-scale problems, than the traditional HOBEM. In addition, the optimal variable of pFFT mesh is recommended to minimize time cost.

A 3D time-domain method for predicting the wave-induced forces and motions of a floating body

This study presents a 3D time-domain Green function method for predicting the wave-induced forces and motions of a floating body. The problem is considered in the cases of both infinite and finite water depths. Based on convolution theory, the convolution of the time-domain Green function and the velocity potential is replaced by the product of their Fourier transformations. A recursive formula for the convolution is proposed in terms of the frequency-domain Green function. This formula provides an efficient way to perform a long-time numerical simulation because the associated computational cost is constant over time. Furthermore, higher-order panels are used to give precise and continuous representations of the body geometry and the velocity potential over the body hull. Computed results for a hemisphere and a Wigley hull are compared with published results, and the comparisons indicate that the various results are in good agreement.

전산유체역학을 이용한 제한수로에서의 선박 침하 해석

When a ship proceeds in confined water, like canal, the water ahead of ship is pushed by hull. This pushed water returns to the side and under the hull, and this returned water will make fluid velocity higher at the side and under the hull, compared to the case in the infinite water depth. Due to the higher velocity, the pressure under the hull will decrease, resulting in the ship drop. This phenomenon is called "ship squat" and ship squat will result in various marine accidents. In this paper, for predicting ship squat, numerical calculation was carried out using commercial CFD code, FLUENT. To confirm wave pattern profile around the ship, VOF(Volume of Fluid) method was applied. The calculated results were compared with other paper's results and empirical methods.

On capillary gravity-wave motion in two-layer fluids

Generalized expansion formulae for the velocity potentials associated with plane gravity-wave problems in the presence of surface tension and interfacial tension are derived in both the cases of finite and infinite water depths in two-layer fluids. As a part of the expansion formulae, orthogonal mode-coupling relations, associated with the eigenfunctions of the velocity potential, are derived. The dispersion relations are analyzed to determine the characteristics of the two propagating modes in the presence of surface and interfacial tension in both the cases of deep-water and shallow-water waves. The expansion formulae are then generalized to deal with boundary-value problems satisfying higher-order boundary conditions at the free surface and interface. As applications of the expansion formulae, the solutions associated with the source potential, forced oscillation and reflection of capillary–gravity waves in the presence of interfacial tension are derived.

On a nonlocal reaction–diffusion–advection equation modelling phytoplankton dynamics

We investigate a reaction–diffusion–advection equation that models the dynamics of a single phytoplankton species in a eutrophic vertical water column. First, we extend the results of Du and Hsu (2010 SIAM J. Math. Anal.42 1305–33) to show that even with variable diffusion and sinking rates, the global dynamics of the model is completely determined by its unique steady-state solution. This implies that the bistable behaviour observed through numerical simulation in Ryabov et al (2010 J. Theor. Biol. 263 120–33) for the phytoplankton dynamics can only occur when one assumes limitation of nutrients in the model. Second, we examine the asymptotic profiles of the positive steady-state solution for small diffusion, large diffusion and deep water column, respectively. Our results reveal that for small diffusion, the phytoplankton population concentrates at the bottom of the water column, while for large diffusion, the population tends to distribute evenly in the water column, and when all the other factors are the same, in a water column with positive background turbidity, the total biomass is bigger in the large diffusion case than in the small diffusion case, and in a water column with zero (or negligible) background turbidity, the total biomass tends to the same limit in both cases; when the water column depth goes to infinity, the population distribution approaches that obtained in Ishii and Takagi (1982 J. Math. Biol. 16 1–24) with infinite water depth, and it reaches a unique maximum at a certain finite water level. We also give a complete answer to a question left open in Hsu and Lou (2010 SIAM J. Appl. Math. 70 245–54) regarding the behaviour of the critical death rate for deep water column, which plays a key role in determining whether a critical water depth exists.

On the Existence and Conditional Energetic Stability of Solitary Gravity-Capillary Surface Waves on Deep Water

This paper presents an existence and stability theory for gravity-capillary solitary waves on the surface of a body of water of infinite depth. Exploiting a classical variational principle, we prove the existence of a minimiser of the wave energy epsilon subject to the constraint I = root 2 mu, where I is the wave momentum and 0 < mu << 1. Since epsilon and I are both conserved quantities a standard argument asserts the stability of the set D-mu of minimisers: solutions starting near D-mu remain close to D-mu in a suitably defined energy space over their interval of existence. In the applied mathematics literature solitary water waves of the present kind are modelled as solutions of the nonlinear Schrodinger equation with cubic focussing nonlinearity. We show that the waves detected by our variational method converge (after an appropriate rescaling) to solutions of this model equation as mu down arrow 0.

A three-dimensional surface wave–ocean circulation coupled model and its initial testing

A theoretical framework to include the influences of nonbreaking surface waves in ocean general circulation models is established based on Reynolds stresses and fluxes terms derived from surface wave-induced fluctuation. An expression for the wave-induced viscosity and diffusivity as a function of the wave number spectrum is derived for infinite and finite water depths; this derivation allows the coupling of ocean circulation models with a wave number spectrum numerical model. In the case of monochromatic surface wave, the wave-induced viscosity and diffusivity are functions of the Stokes drift. The influence of the wave-induced mixing scheme on global ocean circulation models was tested with the Princeton Ocean Model, indicating significant improvement in upper ocean thermal structure and mixed layer depth compared with mixing obtained by the Mellor–Yamada scheme without the wave influence. For example, the model–observation correlation coefficient of the upper 100-m temperature along 35° N increases from 0.68 without wave influence to 0.93 with wave influence. The wave-induced Reynolds stress can reach up to about 5% of the wind stress in high latitudes, and drive 2–3 Sv transport in the global ocean in the form of mesoscale eddies with diameter of 500–1,000 km. The surface wave-induced mixing is more pronounced in middle and high latitudes during the summer in the Northern Hemisphere and in middle latitudes in the Southern Hemisphere.

Explicit Representation for the Infinite-Depth Two-Dimensional Free-Surface Green's Function in Linear Water-Wave Theory

In this paper we derive an explicit representation for the two-dimensional free-surface Green's function in water of infinite depth, based on a finite combination of complex-valued exponential integrals and elementary functions. This representation can easily and accurately be evaluated in a numerical manner, and its main advantage over other representations lies in its simplicity and in the fact that it can be extended towards the complementary half-plane in a straightforward manner. It seems that this extension has not been studied rigorously until now, and it is required when boundary integral equations are extended in the same way. For the computation of the Green's function, the limiting absorption principle and a partial Fourier transform along the free surface are used. Some of its properties are also discussed, and an expression for its far field is developed, which allows us to state appropriately the involved radiation condition. This Green's function is then used to solve the two-dimensional infinite-depth water-wave problem by developing a corresponding boundary integral equation, whose solution is determined by means of the boundary element method. To validate the computations, a benchmark problem based on a half-circle is presented and solved numerically.

Wave interaction with multiple articulated floating elastic plates

Flexural gravity wave scattering by multiple articulated floating elastic plates is investigated in the three cases for water of finite depth, infinite depth and shallow water approximation under the assumptions of two-dimensional linearized theory of water waves. The elastic plates are joined through connectors, which act as articulated joints. In the case when two semi-infinite plates are connected through a single articulation, using the symmetric characteristic of the plate geometry and the expansion formulae for wave-structure interaction problem, the velocity potentials are obtained in closed forms in the case of finite and infinite water depths. On the other hand, in the case of shallow water approximation, the continuity of energy and mass flux are used to obtain a system of equations for the determination of the full velocity potentials for wave scattering by multiple articulations. Further, using the results for single articulation and assuming that the articulated joints are wide apart, the wide-spacing approximation method is used to obtain the reflection coefficient for wave scattering due to multiple articulated floating elastic plates. The effects of the stiffness of the connectors, length of the elastic plates and water depth on the propagation of flexural gravity waves are investigated by analysing the reflection coefficient.

Time-dependent motion of a two-dimensional floating elastic plate

Three methods are presented to determine the motion of a two-dimensional finite elastic plate floating on the water surface, which is released from rest and allowed to evolve freely. The first method is based on a generalized eigenfunction expansion and it is valid for all water depths. The second method is based on an integral equation derived from the Fourier transform, and it is valid for all water depths, although computations are made only for water of infinite depth. These two methods are both based on the frequency-domain solution—however no obvious connection exists between the two methods. The third method is valid only for shallow water, and it expresses the solution as the sum over decaying modes. We present a new derivation of the integral equation for a floating plate based on the Fourier transform of the equations of motion in the time domain. The solution obtained by each method is compared in the appropriate regime, and excellent agreement is found, thereby providing benchmark solutions. We also investigate the regime of validity of the infinite and shallow-depth solutions, and show that both give good results for a quite wide range of depths.

The effect of third-order nonlinearity on statistical properties of random directional waves in finite depth

Abstract. It is well established that third-order nonlinearity produces a strong deviation from Gaussian statistics in water of infinite depth, provided the wave field is long crested, narrow banded and sufficiently steep. A reduction of third-order effects is however expected when the wave energy is distributed on a wide range of directions. In water of arbitrary depth, on the other hand, third-order effects tend to be suppressed by finite depth effects if waves are long crested. Numerical simulations of the truncated potential Euler equations are here used to address the combined effect of directionality and finite depth on the statistical properties of surface gravity waves; only relative water depth kh greater than 0.8 are here considered. Results show that random directional wave fields in intermediate water depths, kh=O(1), weakly deviate from Gaussian statistics independently of the degree of directional spreading of the wave energy.

Applicability of the method of fundamental solutions to 3-D wave–body interaction with fully nonlinear free surface

A numerical model for three-dimensional fully nonlinear free-surface waves is developed by applying a boundary-type meshless approach with a leap-frog time-marching scheme. Adopting Gaussian Radial Basis Functions to fit the free surface, a non-iterative approach to discretize the nonlinear free-surface boundary is formulated. Using the fundamental solutions of the Laplace equation as the solution form of the velocity potential, free-surface wave problems can be solved by collocations at only a few boundary points since the governing equation is automatically satisfied. The accuracy of the present method is verified by comparing the simulated propagation of a solitary wave with an exact solution. The applicability of the present model is illustrated by applying it to the problem of a solitary wave running up on a vertical surface-piercing cylinder and the problem of wave generation in infinite water depth by a submerged moving object.

Transport of energy in polychromatic fluid gravity waves

Potential theory is used to derive closed-form expressions for the energy flux in propagating polychromatic surface gravity waves for water of infinite depth. The energy flux vector has been obtained as convolution-type integrals of the wave-elevation time series at one point on the surface. An expression for the integrated flux through a vertical surface is also given. It is shown that the integrated flux cannot be negative at any instant in time.

HF Radar Sea-echo from Shallow Water.

HF radar systems are widely and routinely used for the measurement of ocean surface currents and waves. Analysis methods presently in use are based on the assumption of infinite water depth, and may therefore be inadequate close to shore where the radar echo is strongest. In this paper, we treat the situation when the radar echo is returned from ocean waves that interact with the ocean floor. Simulations are described which demonstrate the effect of shallow water on radar sea-echo. These are used to investigate limits on the existing theory and to define water depths at which shallow-water effects become significant. The second-order spectral energy increases relative to the first-order as the water depth decreases, resulting in spectral saturation when the waveheight exceeds a limit defined by the radar transmit frequency. This effect is particularly marked for lower radar transmit frequencies. The saturation limit on waveheight is less for shallow water. Shallow water affects second-order spectra (which gives wave information) far more than first-order (which gives information on current velocities), the latter being significantly affected only for the lowest radar transmit frequencies for extremely shallow water. We describe analysis of radar echo from shallow water measured by a Rutgers University HF radar system to give ocean wave spectral estimates. Radar-derived wave height, period and direction are compared with simultaneous shallow-water in-situ measurements.

Numerical simulation of non-linear regular and focused waves in an infinite water-depth

Inviscid three-dimensional free surface wave motions are simulated using a novel quadratic higher order boundary element model (HOBEM) based on potential theory for irrotational, incompressible fluid flow in an infinite water-depth. The free surface boundary conditions are fully non-linear. Based on the use of images, a channel Green function is developed and applied to the present model so that two lateral surfaces of an infinite-depth wave tank can be excluded from the calculation domain. In order to generate incident waves and dissipate outgoing waves, a non-reflective wave generator, composed of a series of vertically aligned point sources in the computational domain, is used in conjunction with upstream and downstream damping layers. Numerical experiments are carried out, with linear and fully non-linear, regular and focused waves. It can be seen from the results that the present approach is effective in generating a specified wave profile in an infinite water-depth without reflection at the open boundaries, and fully non-linear numerical simulations compare well with theoretical solutions. The present numerical technique is aimed at efficient modelling of the non-linear wave interactions with ocean structures in deep water.

Expansion formulae for wave structure interaction problems with applications in hydroelasticity

Alternate derivations of the expansion formulae for wave structure interaction problems are obtained in case of water of infinite depth and utilized to analyze the hydroelastic behavior of large floating structures. Considering the boundary value problem associated with Laplace equation having higher order boundary condition on the horizontal boundary and a Dirichlet type boundary condition on the vertical boundary in a quarter plane, Fourier sine transform is applied in the horizontal direction to convert the problem to a Sturm–Liouville type boundary value problem associated with non-homogeneous ordinary differential equation (ODE) in the transformed variable. Finally, inverting the transformed functions and applying the regularity criterion of the transformed function, the required expansion formula is derived. The expansion formula thus derived is extended to deal with similar boundary value problems having Neumann type boundary condition. The expansion formulae are applied to (i) analyze oblique scattering of flexural gravity waves by an articulated floating elastic plate and (ii) study the effect of compression on the oblique scattering of flexural gravity waves by a line discontinuity in a large floating ice sheet in water of infinite depth, which find applications in marine technology and arctic engineering, respectively. The present derivations of the expansion formulae are very simple and straightforward and can be easily used to study a large class of problems in the area of fluids and structures in mathematical physics and engineering.