Generation and Propagation of Long-Crested Finite-Depth Surface Waves: Comparison Between the Corresponding Results From Numerical Models and a Wave Tank

Author(s):  
M. Hasanat Zaman ◽  
Wade Parsons ◽  
Okey Nwogu ◽  
Wooyoung Choi ◽  
R. Emile Baddour ◽  
...  

The evolution of long-crested surface waves subject to side-band perturbations is investigated with two different numerical models: a direct solver for the Euler equations using a non-orthogonal boundary-fitted curvilinear coordinate system and an FFT-accelerated boundary integral method. The numerical solutions are then validated with laboratory experiments performed in the NRC-IOT Ocean Engineering Basin with a segmented wave-maker operating in piston mode. The numerical models are forced by a point measurement of the free surface elevation at a wave probe close to the wave-maker and the numerical solutions are compared with the measured time-series of the surface elevation at a few wave probe locations downstream.

2001 ◽  
Vol 42 (3) ◽  
pp. 354-371
Author(s):  
S. W. McCue ◽  
L. K. Forbes

AbstractThe free surface flow of a finite depth fluid past a semi-infinite body is considered. The fluid is assumed to have constant vorticity throughout and the free surface is assumed to attach smoothly to the front face of the body. Numerical solutions are found using a boundary integral method in the physical plane and it is shown that solutions exist for all supercritical Froude numbers. The related problem of the cusp-like flow due to a submerged sink in a comer is also considered. Vorticity is included in the flow and it is shown that the behaviour of the solutions is qualitatively the same as that found in the problem described above.


1998 ◽  
Vol 356 ◽  
pp. 93-124 ◽  
Author(s):  
HARRIS WONG ◽  
DAVID RUMSCHITZKI ◽  
CHARLES MALDARELLI

The expansion and contraction of a bubble pinned at a submerged tube tip and driven by constant gas flow rate Q are studied both theoretically and experimentally for Reynolds number Re[Lt ]1. Bubble shape, gas pressure, surface velocities, and extrapolated detached bubble volume are determined by a boundary integral method for various Bond (Bo=ρga2/σ) and capillary (Ca=μQ/σa2) numbers, where a is the capillary radius, ρ and μ are the liquid density and viscosity, σ is the surface tension, and g is the gravitational acceleration.Bubble expansion from a flat interface to near detachment is simulated for a full range of Ca (0.01–100) and Bo (0.01–0.5). The maximum gas pressure is found to vary almost linearly with Ca for 0.01[les ]Ca[les ]100. This correlation allows the maximum bubble pressure method for measuring dynamic surface tension to be extended to viscous liquids. Simulated detached bubble volumes approach static values for Ca[Lt ]1, and asymptote as Q3/4 for Ca[Gt ]1, in agreement with analytic predictions. In the limit Ca→0, two singular time domains are identified near the beginning and the end of bubble growth during which viscous and capillary forces become comparable.Expansion and contraction experiments were conducted using a viscous silicone oil. Digitized video images of deforming bubbles compare well with numerical solutions. It is observed that a bubble contracting at high Ca snaps off.


Author(s):  
Sylvain Barbot

ABSTRACT Numerical models of rupture dynamics provide great insights into the physics of fault failure. However, resolving stress interactions among multiple faults remains challenging numerically. Here, we derive the elastostatic Green’s functions for stress and displacement caused by arbitrary slip distributions along multiple parallel faults. The equations are derived in the Fourier domain, providing an efficient means to calculate stress interactions with the fast Fourier transform. We demonstrate the relevance of the method for a wide range of applications, by simulating the rupture dynamics of single and multiple parallel faults controlled by a rate- and state-dependent frictional contact, using the spectral boundary integral method and the radiation-damping approximation. Within the antiplane strain approximation, we show seismic cycle simulations with a power-law distribution of rupture sizes and, in a different parameter regime, sequences of seismogenic slow-slip events. Using the in-plane strain approximation, we simulate the rupture dynamics of a restraining stepover. Finally, we describe cycles of large earthquakes along several parallel strike-slip faults in three dimensions. The approach is useful to explore the dynamics of interacting or isolated faults with many degrees of freedom.


2002 ◽  
Vol 454 ◽  
pp. 403-417 ◽  
Author(s):  
A. R. CHAMPNEYS ◽  
J.-M. VANDEN-BROECK ◽  
G. J. LORD

This paper extends the numerical results of Hunter & Vanden-Broeck (1983) and Vanden-Broeck (1991) which were concerned with studies of solitary waves on the surface of fluids of finite depth under the action of gravity and surface tension. The aim of this paper is to answer the question of whether small-amplitude elevation solitary waves exist. Several analytical results have proved that bifurcating from Froude number F = 1, for Bond number τ between 0 and 1/3, there are families of ‘generalized’ solitary waves with periodic tails whose minimum amplitude is an exponentially small function of F−1. An open problem (which, for τ sufficiently close to 1/3, was recently proved by Sun 1999 to be false) is whether this amplitude can ever be zero, which would give a truly localized solitary wave.The problem is first addressed in terms of model equations taking the form of generalized fifth-order KdV equations, where it is demonstrated that if such a zero-tail-amplitude solution occurs, it does so along codimension-one lines in the parameter plane. Moreover, along solution paths of generalized solitary waves a topological distinction is found between cases where the tail does vanish and those where it does not. This motivates a new set of numerical results for the full problem, formulated using a boundary integral method, namely to probe the size of the tail amplitude as τ varies for fixed F > 1. The strong conclusion from the numerical results is that true solitary waves of elevation do not exist for the steady gravity–capillary water wave problem, at least for 9/50 < τ < 1=3. This finding confirms and explains previous asymptotic results by Yang & Akylas.


1992 ◽  
Vol 36 (01) ◽  
pp. 38-54
Author(s):  
Zhi Guo ◽  
Allen T. Chwang

The general planar translation of two bodies of revolution through an inviscid and incompressible fluid is considered. The moving trajectories and the hydrodynamic interactions are computed based on the generalized Lagrange equations of motion, including the effects of solid constraints, external forces in the plane of motion, and a uniform stream in any direction parallel to the plane of motion. In a relative coordinate system moving with the stream, the kinetic energy of the fluid is expressed as a function of six added masses due to motions parallel and perpendicular to the line joining the centers of two bodies. Analytical solutions of added masses in series form are obtained for the motion of two spheres. A new iterative formula based on the analysis of velocity potentials around each body is developed for added masses and their derivatives with respect to the separation distance due to the transverse motion. The method of successive images and Taylor's added-mass formula are applied to determine the added masses and their derivatives due to the centroidal motion. These results are compared with the numerical solution of added masses computed by the boundary-integral method and the generalized Taylor added-mass formula. The integral equations, in terms of surface-source distributions on both surfaces, are modified for obtaining accurate numerical solutions. Numerical results are given for several practical engineering problems.


1986 ◽  
Vol 170 ◽  
pp. 479-497 ◽  
Author(s):  
J. R. Blake ◽  
B. B. Taib ◽  
G. Doherty

The growth and collapse of transient vapour cavities near a rigid boundary in the presence of buoyancy forces and an incident stagnation-point flow are modelled via a boundary-integral method. Bubble shapes, particle pathlines and pressure contours are used to illustrate the results of the numerical solutions. Migration of the collapsing bubble, and subsequent jet formation, may be directed either towards or away from the rigid boundary, depending on the relative magnitude of the physical parameters. For appropriate parameter ranges in stagnation-point flow, unusual ‘hour-glass’ shaped bubbles are formed towards the end of the collapse of the bubble. It is postulated that the final collapsed state of the bubble may be two toroidal bubbles/ring vortices of opposite circulation. For buoyant vapour cavities the Kelvin impulse is used to obtain criteria which determine the direction of migration and subsequent jet formation in the collapsing bubble.


1999 ◽  
Vol 399 ◽  
pp. 277-300 ◽  
Author(s):  
SCOTT W. McCUE ◽  
LAWRENCE K. FORBES

Free surface flows of a rotational fluid past a two-dimensional semi-infinite body are considered. The fluid is assumed to be inviscid, incompressible, and of finite depth. A boundary integral method is used to solve the problem for the case where the free surface meets the body at a stagnation point. Supercritical solutions which satisfy the radiation condition are found for various values of the Froude number and the dimensionless vorticity. Subcritical solutions are also found; however these solutions violate the radiation condition and are characterized by a train of waves upstream. It is shown numerically that the amplitude of these waves increases as each of the Froude number, vorticity and height of the body above the bottom increases.


Author(s):  
P. F. Rhodes-Robinson

AbstractThe classical wave-maker problem to determine the forced two-dimensional wave motion with outgoing surface waves at infinity generated by a harmonically oscillating vertical plane wave-maker immersed in water was solved long ago by Sir Thomas Havelock. In this paper we reinvestigate the problem, making allowance for the presence of surface tension which was excluded before, and obtain a solution of the boundary-value problem for the velocity potential which is made unique by prescribing the free surface slope at the wave-maker. The cases of both infinite and finite constant depth are treated, and it is essential to employ a method which is new to this problem since the theory of Havelock cannot be extended in the latter case of finite depth. The solution of the corresponding problem concerning the axisymmetric wave motion due to a vertical cylindrical wave-maker is deduced in conclusion.


2021 ◽  
Vol 929 ◽  
Author(s):  
Anatoliy Khait ◽  
Zhihua Ma

We present an investigation of the fundamental physical processes involved in deep-water gravity wave breaking. Our motivation is to identify the underlying reason causing the deficiency of the eddy viscosity breaking model (EVBM) in predicting surface elevation for strongly nonlinear waves. Owing to the limitation of experimental methods in the provision of high-resolution flow information, we propose a numerical methodology by developing an EVBM enclosed standalone fully nonlinear quasi-potential (FNP) flow model and a coupled FNP plus Navier–Stokes flow model. The numerical models were firstly verified with a wave train subject to modulational instability, then used to simulate a series of broad-banded focusing wave trains under non-, moderate- and strong-breaking conditions. A systematic analysis was carried out to investigate the discrepancies of numerical solutions produced by the two models in surface elevation and other important physical properties. It is found that EVBM predicts accurately the energy dissipated by breaking and the amplitude spectrum of free waves in terms of magnitude, but fails to capture accurately breaking induced phase shifting. The shift of phase grows with breaking intensity and is especially strong for high-wavenumber components. This is identified as a cause of the upshift of the wave dispersion relation, which increases the frequencies of large-wavenumber components. Such a variation drives large-wavenumber components to propagate at nearly the same speed, which is significantly higher than the linear dispersion levels. This suppresses the instant dispersive spreading of harmonics after the focal point, prolonging the lifespan of focused waves and expanding their propagation space.


Sign in / Sign up

Export Citation Format

Share Document