An interfacial Gerstner-type trapped wave

Coastally trapped rotational interfacial waves are studied theoretically by using a Lagrangian formulation of fluid motion in a rotating ocean. The waves propagate along the interface between two immiscible inviscid incompressible fluid layers of finite depths and different densities, and are trapped at a straight wall due to the Coriolis force. For layers of finite depth, solutions are sought as series expansions after a small parameter. Comparison is made with the irrotational interfacial Kelvin wave. Both types of waves are identical to first order, having zero vorticity. The second order solution yields a relation between the vorticity and the velocity shear in the wave motion. Requiring that the mean motion in both layers is irrotational, then follows the well-known Stokes drift for interfacial Kelvin waves. On the other hand, if the mean forward drift is identically zero, we obtain the second order vorticity in the Gerstner-type wave. The solutions in both layers for the Gerstner-type interfacial wave are given analytically to second order. It is shown that small density differences and thin upper layers both act to yield a shape of the material interfacial with broader crests and sharper troughs. These effects also tend to make the particle trajectories at the interface in both layers become distorted ellipses which are flatter on the upper side. It is concluded that the effect of air excludes the possibility of observing the exact Gerstner wave in deep water.     

Particle trajectory and mass transport of finite-amplitude waves in water of uniform depth

A set of governing equations in Lagrangian form is derived for propagating gravity waves in water of uniform depth. The Lindstedt–Poincaré perturbation method is used to obtain approximations up to fifth order. Recognizing the Lagrangian frequency to be a position function for all particles is a key to find these higher-order approximations. The present solution has zero pressure at the free surface and satisfies exactly the dynamic boundary condition. Under the present approximations, the Lagrangian frequency is composed of two parts. The first part is constant for all particles and equivalent to the term in the fifth-order Stokes' wave theory [J.D. Fenton, A fifth-order Stokes theory for steady waves, J. Waterway, Port, Coastal Ocean Eng. 111 (1985) 216–234]. The second part is a function of the depth. All the particles move as open (nonclosed) loops and have mean drift displacements that decrease exponentially with the water depth. Thus, a new fourth-order mass transport velocity is found.     

CFD simulation of two‐ and three‐dimensional free‐surface flow

The paper describes the implementation of moving‐mesh and free‐surface capabilities within a 3‐d finite‐volume Reynolds‐averaged‐Navier–Stokes solver, using surface‐conforming multi‐block structured meshes. The free‐surface kinematic condition can be applied in two ways: enforcing zero net mass flux or solving the kinematic equation by a finite‐difference method. The free surface is best defined by intermediate control points rather than the mesh vertices. Application of the dynamic boundary condition to the piezometric pressure at these points provides a hydrostatic restoring force which helps to eliminate any unnatural free‐surface undulations. The implementation of time‐marching methods on moving grids are described in some detail and it is shown that a second‐order scheme must be applied in both scalar‐transport and free‐surface equations if flows driven by free‐surface height variations are to be computed without significant wave attenuation using a modest number of time steps. Computations of five flows of theoretical and practical interest—forced motion in a pump, linear waves in a tank, quasi‐1d flow over a ramp, solitary wave interaction with a submerged obstacle and 3‐d flow about a surface‐penetrating cylinder—are described to illustrate the capabilities of our code and methods. Copyright © 2003 John Wiley & Sons, Ltd.     

Bleustein–Gulyaev waves in 6 mm piezoelectric materials loaded with a viscous liquid layer of finite thickness

In this paper, we analyze the propagation of Bleustein–Gulyaev waves in an unbounded piezoelectric half-space loaded with a viscous liquid layer of finite thickness within the linear elastic theories. Exact solutions of the phase velocity equations are obtained in the cases of both electrically open circuit and short circuit by solving the equilibrium equations of piezoelectric materials and the diffusion equation of viscous liquid. A PZT-5H/Glycerin system is selected to perform the numerical calculation. The results show that the mass density and the viscous coefficient have different effects on the propagation attenuation and phase velocity under different electrical boundary conditions. In particular, the penetration depth of the waves is of the same order as the wavelength in the case of electrically short circuit. These effects can be used to manipulate the behavior of the waves and have implications in the application of acoustic wave devices.     

Nonlinear vibrations of a vapor film in the presence of intense heat fluxes

Waves propagating along the interface between a thin vapor film and a liquid layer in the presence of a heat flux are investigated. The boundary conditions on the vapor-liquid phase surface take into account the temperature dependence of the pressure and the possibilities of formation of the metastable state of the superheated liquid and mass flow. Variations in the saturation pressure as functions of the temperature and mass flux lead to generation of weakly damped periodic waves of low amplitude whose velocity can be much higher than the velocity of the gravity waves. The waves ensure stability of the vapor film beneath the liquid layer in the gravity field. The finite-amplitude waves on the surface of the vapor film differ from the Stokes surface waves on the free surface of isothermal fluid. Instability regimes related with superheating of the liquid ant its explosive boiling when the amplitude of an initially small wave increases to infinity in a finite time can develop in a certain working-parameter regime.     

Boundary layer flow of Maxwell fluid due to torsional motion of cylinder: modeling and simulation

This paper investigates the boundary layer flow of the Maxwell fluid around a stretchable horizontal rotating cylinder under the influence of a transverse magnetic field. The constitutive flow equations for the current physical problem are modeled and analyzed for the first time in the literature. The torsional motion of the cylinder is considered with the constant azimuthal velocity E. The partial differential equations (PDEs) governing the torsional motion of the Maxwell fluid together with energy transport are simplified with the boundary layer concept. The current analysis is valid only for a certain range of the positive Reynolds numbers. However, for very large Reynolds numbers, the flow becomes turbulent. Thus, the governing similarity equations are simplified through suitable transformations for the analysis of the large Reynolds numbers. The numerical simulations for the flow, heat, and mass transport phenomena are carried out in view of the bvp4c scheme in MATLAB. The outcomes reveal that the velocity decays exponentially faster and reduces for higher values of the Reynolds numbers and the flow penetrates shallower into the free stream fluid. It is also noted that the phenomenon of stress relaxation, described by the Deborah number, causes to decline the flow fields and enhance the thermal and solutal energy transport during the fluid motion. The penetration depth decreases for the transport of heat and mass in the fluid with the higher Reynolds numbers. An excellent validation of the numerical results is assured through tabular data with the existing literature.     

Discrete particle simulation of mixed sand transport

An Eulerian/Lagrangian numerical simulation is performed on mixed sand transport. Volume averaged Navier–Stokes equations are solved to calculate gas motion, and particle motion is calculated using Newton's equation, involving a hard sphere model to describe particle-to-particle and particle-to-wall collisions. The influence of wall characteristics, size distribution of sand particles and boundary layer depth on vertical distribution of sand mass flux and particle mean horizontal velocity is analyzed, suggesting that all these three factors affect sand transport at different levels. In all cases, for small size groups, sand mass flux first increases with height and then decreases while for large size groups, it decreases exponentially with height and for middle size groups the behavior is in-between. The mean horizontal velocity for all size groups well fits experimental data, that is, increasing logarithmically with height in the middle height region. Wall characteristics greatly affects particle to wall collision and makes the flat bed similar to a Gobi surface and the rough bed similar to a sandy surface. Particle size distribution largely affects the sand mass flux and the highest heights they can reach especially for larger particles.     

Turbulent Dynamics of a Critically Reflecting Internal Gravity Wave

Results from computational fluid dynamics experiments of internal wave reflection from sloping boundaries are presented. In these experiments the incident wave lies in a plane normal to the slope. When the angle of wave energy propagation is close to the bottom slope the reflection causes wave breakdown into a quasi-periodic, turbulent boundary layer. Boundary layer energetics and vorticity dynamics are examined and indicate the importance of the three-dimensional turbulence. The boundary layer exhibits intermittent turbulence: approximately every 1.2 wave periods the boundary layer mixes energetically for a duration of about one-third of a wave period, and then it restratifies until the next mixing event. Throughout the wave cycle a strong thermal front is observed to move upslope at the phase speed of the incident waves. Simulations demonstrate that the net effects of turbulent mixing are not confined to the boundary layer, but are communicated to the interior stratified fluid by motions induced by buoyancy effects and by the wave field, resulting in progressive weakening of the background density gradient. Transition to turbulence is determined to occur at Reynolds numbers of approximately 1500, based upon the wavelength and maximum current velocity of the oncoming wave train. The boundary layer thickness depends on the Reynolds number for low Richardson numbers, with a characteristic depth of approximately one-half of the vertical wavelength of the oncoming wave. Received 21 May 1997 and accepted 14 October 1997     

Composite steady waves in a gravity fluid stream of finite depth

The problem of plane steady gravitational waves of finite amplitude, caused by a periodically distributed pressure over the surface of an ideal incompressible gravity fluid stream of finite depth, is considered. It is assumed that these waves do not vanish as the pressure becomes constant, but become free waves, which exist at constant pressure and special values of the stream velocity. As in [1], where a stream of finite depth is considered, such waves will be designated composite as contrasted with forced waves which vanish together with the variable part of the pressure. A general method is given for computing the composite wave characteristics. The first three approximations are computed to the end. An approximate equation for the wave profile is found.     

Simulation of particles released near the wall in a turbulent boundary layer

A direct numerical simulation was used along with a Lagrangian particle tracking technique to study particle motion in a horizontal, spatially developing turbulent boundary layer along an upper-wall (with terminal velocity directed away from the wall). The objective of the research was to study particle diffusion, dispersion, reflection, and mean velocity in the context of two parametric studies: one investigated the effect of the drift parameter (the ratio of particle terminal velocity to fluid friction velocity) for a fixed and finite particle inertia, and the second varied the drift parameter and particle inertia by the same amount (i.e. for a constant Froude number). A range of drift parameters from 10⁻⁴ to 10⁰ were considered for both cases. The particles were injected into the simulation at a height of four wall units for several evenly distributed points across the span and a perfectly elastic wall collision was specified at one wall unit.Statistics collected along the particle trajectories demonstrated a transition in particle movement from one that is dominated by diffusion to one that is dominated by gravity. For small and intermediate sized particles (i.e. ones with outer Stokes numbers and drift parameters much less than unity) transverse diffusion away from the wall dominated particle motion. However, preferential concentration is seen near the wall for intermediate-sized particles due to inhomogeneous turbulence effects (turbophoresis), consistent with previous channel flow studies. Particle–wall collision statistics indicated that impact velocities tended to increase with increasing terminal velocity for small and moderate inertias, after which initial conditions become important. Finally, high relative velocity fluctuations (compared to terminal velocity) were found as particle inertia increased, and were well described with a quasi-one-dimensional fluctuation model.     

Gravity-capillary waves in the presence of constant vorticity

Periodic and solitary gravity-capillary waves propagating at a constant velocity at the surface of a fluid of finite depth are considered. The vorticity in the fluid is assumed to be constant. Analytical solutions are presented for waves of small amplitude. For waves of large amplitude, numerical solutions are computed by boundary integral equation methods. The results unify previous findings for irrotational gravity capillary waves and gravity waves with constant vorticity. In particular solitary waves with oscillatory tails and branches of solutions which exist only for waves of large amplitude are found.     

Unsteady RANS method for surface ship boundary layer and wake and wave field

Results are reported of an unsteady Reynolds‐averaged Navier–Stokes (RANS) method for simulation of the boundary layer and wake and wave field for a surface ship advancing in regular head waves, but restrained from body motions. Second‐order finite differences are used for both spatial and temporal discretization and a Poisson equation projection method is used for velocity–pressure coupling. The exact kinematic free‐surface boundary condition is solved for the free‐surface elevation using a body‐fitted/free‐surface conforming grid updated in each time step. The simulations are for the model problem of a Wigley hull advancing in calm water and in regular head waves. Verification and validation procedures are followed, which include careful consideration of both simulation and experimental uncertainties. The steady flow results are comparable to other steady RANS methods in predicting resistance, boundary layer and wake, and free‐surface effects. The unsteady flow results cover a wide range of Froude number, wavelength, and amplitude for which first harmonic amplitude and phase force and moment experimental data are available for validation along with frequency domain, linear potential flow results for comparisons. The present results, which include the effects of turbulent flow and non‐linear interactions, are in good agreement with the data and overall show better capability than the potential flow results. The physics of the unsteady boundary layer and wake and wave field response are explained with regard to frequency of encounter and seakeeping theory. The results of the present study suggest applicability for additional complexities such as practical ship geometry, ship motion, and maneuvering in arbitrary ambient waves. Copyright © 2001 John Wiley & Sons, Ltd.     

Theory of steady waves in horizontal flow with a linear velocity profile

The structure and characteristics of nonlinear steady waves on the surface of horizontal shear flow of an ideal homogeneous incompressible fluid of finite depth with a linear velocity profile are studied using two-dimensional theory and the Euler approach. The wave motion is considered irrotational. A modification of the first Stokes method is proposed that allows algebraic calculations of terms of perturbation series. Nonlinear dispersion relations are obtained and analyzed for both upstream and downstream traveling waves. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 43–48, May–June, 2006.     

Fully non‐linear free‐surface simulations by a 3D viscous numerical wave tank

A finite difference scheme using a modified marker‐and‐cell (MAC) method is applied to investigate the characteristics of non‐linear wave motions and their interactions with a stationary three‐dimensional body inside a numerical wave tank (NWT). The Navier–Stokes (NS) equation is solved for two fluid layers, and the boundary values are updated at each time step by a finite difference time marching scheme in the frame of a rectangular co‐ordinate system. The viscous stresses and surface tension are neglected in the dynamic free‐surface condition, and the fully non‐linear kinematic free‐surface condition is satisfied by the density function method developed for two fluid layers. The incident waves are generated from the inflow boundary by prescribing a velocity profile resembling flexible flap wavemaker motions, and the outgoing waves are numerically dissipated inside an artificial damping zone located at the end of the tank. The present NS–MAC NWT simulations for a vertical truncated circular cylinder inside a rectangular wave tank are compared with the experimental results of Mercier and Niedzwecki, an independently developed potential‐based fully non‐linear NWT, and the second‐order diffraction computation. Copyright © 1999 John Wiley & Sons, Ltd.     

Viscous fluid flow induced by rotational-oscillatory motion of a porous sphere

Viscous fluid flow induced by rotational-oscillatorymotion of a porous sphere submerged in the fluid is determined. The Darcy formula for the viscous medium drag is supplementedwith a term that allows for the medium motion. The medium motion is also included in the boundary conditions. Exact analytical solutions are obtained for the time-dependent Brinkman equation in the region inside the sphere and for the Navier–Stokes equations outside the body. The existence of internal transverse waves in the fluid is shown; in these waves the velocity is perpendicular to the wave propagation direction. The waves are standing inside the sphere and traveling outside of it. The particular cases of low and high oscillation frequencies are considered.     

A coupled boundary element–finite difference model of surface wave motion over a wall turbulent flow

An effective numerical technique is presented to model turbulent motion of a standing surface wave in a tank. The equations of motion for turbulent boundary layers at the solid surfaces are coupled with the potential flow in the bulk of the fluid, and a mixed BEM–finite difference technique is used to model the wave motion and the corresponding boundary layer flow. A mixing‐length theory is used for turbulence modelling. The model results are in good agreement with previous physical and numerical experiments. Although the technique is presented for a standing surface wave, it can be easily applied to other free surface problems. Copyright © 2005 John Wiley & Sons, Ltd.     

Oscillatory Motion of a Viscous Fluid in Contact with a Flat Layer of a Porous Medium

Analytical solutions are obtained for two problems of transverse internal waves in a viscous fluid contacting with a flat layer of a fixed porous medium. In the first problem, the waves are considered which are caused by the motion of an infinite flat plate located on the fluid surface and performing harmonic oscillations in its plane. In the second problem, the waves are caused by periodic shear stresses applied to the free surface of the fluid. To describe the fluid motion in the porous medium, the unsteady Brinkman equation is used, and the motion of the fluid outside the porous medium is described by the Navier–Stokes equation. Examples of numerical calculations of the fluid velocity and filtration velocity profiles are presented. The existence of fluid layers with counter-directed velocities is revealed.     

Obtainable accuracy in the solution of practical problems by small-amplitude wave theory

Results obtainable using the theory of progressive waves of small amplitude and certain fundamental solutions relating to waves of finite height are investigated. The theoretical findings are compared with existing experimental data. It is established that the best agreement between the theoretical and experimental profiles of a plane wave is achieved with constructions based on Kozhevnikov's [1] graphs and the equations of motion in the second approximation with respect to the wave height in the form proposed by Mich [2]. The limits within which it is expedient to use the theoretical formulas of the theory of small-amplitude waves and the theory of the second approximation with respect to the wave height are found and proved for the particle velocity, excess pressure, energy flux, and the energy of a single wave.     

Waves due to an oscillatory pressure on a stratifield fluid: A uniformly valid asymptotic solution

The three-dimensional axisymmetrical initial-value problem of waves in a two-layered fluid of finite depth by an oscillatory surface pressure is solved. The exact integral solutions for velocity potentials of each layer and wave elevations at the surface and interface are obtained. The uniform asymptotic analysis of the unsteady state of waves is carried out when lower fluid is of infinite depth.