Water wave radiation by a sphere submerged in water with an ice-cover

Using the multipoles method, we formulate the problems of radiation (both heave and sway) of water waves by a submerged sphere in deep as well as in uniform finite depth water with an ice-cover, with the ice-cover being modelled as an elastic plate of very small thickness. In each case this leads to an infinite system of linear equations which are solved numerically by standard techniques. The added-mass and damping coefficients for a heaving and swaying sphere are obtained and depicted graphically against the wave number for various values of the radius of the submerged sphere and flexural rigidity of the ice-cover to show the effect of the presence of ice-cover on these quantities. When the flexural rigidity is taken to be zero, the numerical results for the added-mass and damping coefficient for water with a free surface are recovered.     

Fully nonlinear modeling of radiated waves generated by floating flared structures

The nonlinear radiated waves generated by a structure in forced motion, are simulated numerically based on the potential theory. A fully nonlinear numerical model is developed by using a higher-order boundary element method （HOBEM）. In this model, the instantaneous body position and the transient free surface are updated at each time step. A Lagrangian technique is employed as the time marching scheme on the free surface. The mesh regridding and interpolation methods are adopted to deal with the possible numerical instability. Several auxiliary functions are proposed to calculate the wave loads indirectly, instead of directly predicting the temporal derivative of the velocity potential. Numerical experiments are carried out to simulate the heave motions of a submerged sphere in infinite water depth, the heave and pitch motions of a truncated flared cylinder in finite depth. The results are verified against the published numerical results to ensure the effectiveness of the proposed model. Moreover, a series of higher harmonic waves and force components are obtained by the Fourier transformation to investigate the nonlinear effect of oscillation frequency. The difference among fully nonlinear, body-nonlinear and linear results is analyzed. It is found that the nonlinearity due to free surface and body surface has significant influences on the numerical results of the radiated waves and forces.     

Oscillations of a Circular Cylinder in a Linearly Stratified Fluid

The plane problem of the motions of a three-layer fluid initiated by the oscillations of a circular cylinder is solved in the linear formulation in the Boussinesq approximation. The cylinder is completely immersed in the linearly stratified middle layer, and the upper and lower layers are homogeneous and bounded by rigid horizontal walls. The fluid is assumed to be ideal and incompressible. The added mass and damping coefficients are calculated as functions of the oscillation frequency of the cylinder and the layer thicknesses.     

Experiments to study interactions between baroclinic lower flows and a stably stratified upper layer

Experiments were conducted on a rotating fluid annulus to study the basic interactions between baroclinic lower flows and a stably stratified upper layer. Sufficiently stable stratification is necessary for steady flows to emerge in the lower layer. Upward fluid motions make the baroclinic flows permeate into the upper layer. The stable stratification, however, suppresses upward motions so that zonal fluid velocities decrease with height. In fact, their maximum appears at the top level of the baroclinic lower layer and the sign of the radial temperature gradient changes there; namely, it is warmer on the inner side of the annulus in the upper layer. This temperature profile is reflected in a meridional fluid circulation mixing both layers. In the upper layer of the wave flow, there exists a critical level below and above which the zonal fluid velocities have opposite directions for the wave to have a phase shift of half a wavelength in appearance. The experimental results correspond to real atmospheric phenomena.     

An implicit three‐dimensional fully non‐hydrostatic model for free‐surface flows

An implicit method is developed for solving the complete three‐dimensional (3D) Navier–Stokes equations. The algorithm is based upon a staggered finite difference Crank‐Nicholson scheme on a Cartesian grid. A new top‐layer pressure treatment and a partial cell bottom treatment are introduced so that the 3D model is fully non‐hydrostatic and is free of any hydrostatic assumption. A domain decomposition method is used to segregate the resulting 3D matrix system into a series of two‐dimensional vertical plane problems, for each of which a block tri‐diagonal system can be directly solved for the unknown horizontal velocity. Numerical tests including linear standing waves, nonlinear sloshing motions, and progressive wave interactions with uneven bottoms are performed. It is found that the model is capable to simulate accurately a range of free‐surface flow problems using a very small number of vertical layers (e.g. two–four layers). The developed model is second‐order accuracy in time and space and is unconditionally stable; and it can be effectively used to model 3D surface wave motions. Copyright © 2004 John Wiley & Sons, Ltd.     

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.     

Evolution of the diffusion-induced flow over a sphere submerged in a continuously stratified fluid

The problem of the stabilization of the diffusion-induced flow over a sphere submerged in a continuously stratified fluid is solved using both asymptotic and numerical methods. The analytical solution describes the structure of the main convective cells, including thin meridional jets flowing along the surface and plumes spreading from the flow convergence regions above the upper and lower poles of the sphere which gradually return the fluid particles to the neutral buoyancy horizon. The total width of the flows adjacent to the surface exceeds the thickness of the salinity deficit layer or the density boundary layer. The numerical solution of the complete problem in the nonlinear formulation describes the main convective cells and two systems of unsteady integral waves formed in the vicinity of the sphere poles. At large times, out of the entire system of internal waves only those nearest to the neighborhood of their horizon of formation remain clearly defined. The calculated flow patterns are in agreement with each other and the data of shadow visualization of the stratified fluid structure near a submerged obstacle at rest.     

Hydrodynamic instability in a porous layer saturated with a heat generating fluid

Critical Rayleigh numbers determined by linear stabiliy theory are presented for porous-fluid layers of infinite horizontal extent heated internally by a uniform volumetric energy source in the fluid. The thermal coupling between the layer and its environment is represented by a general mixed boundary condition for both the conduction state and the disturbance temperature. Rigid-rigid, rigid-constant pressure, and constant pressure-rigid boundaries are considered in the computations. For a fixed ratio of upper surface Biot number to that at the lower surface, decreasing the Biot number is strictly destabilizing for values of this ratio greater than or equal to one. A layer with a rigid upper surface is generally the most stable; however, a layer with a rigid upper surface and a constant pressure lower surface exhibits the largest values of critical Rayleigh numbers for large values of Biot number.     

Free-field wave solutions in a half-plane exhibiting a special-type of continuous inhomogeneity

This work presents closed-form solutions for free-field motions in a continuously inhomogeneous half-plane that include contributions of incident waves as well as of waves reflected from the traction-free horizontal surface. Both pressure and vertically polarized shear waves are considered. Furthermore, two special types of material inhomogeneity are studied, namely (a) a shear modulus that varies quadratically with respect to the depth coordinate and (b) one that varies exponentially with the said coordinate. In all cases, Poisson’s ratio is fixed at one-quarter, while both shear modulus and material density profiles vary proportionally. Next, a series of numerical results serve to validate the aforementioned models, and to show the differences in the wave motion patterns developing in media that are inhomogeneous as compared to a reference homogeneous background. These results clearly show the influence of inhomogeneity, as summarized by a single material parameter, on the free-field motions that develop in the half-plane. It is believed that this type of information is useful within the context of wave propagation studies in non-homogeneous continua, which in turn find applications in fields as diverse as laminated composites, geophysical prospecting, oil exploration and earthquake engineering.     

Wave regimes of ascending flow of a thin layer of a viscous liquid in contact with a gas

The flow of a liquid in thin layers is one of the hydrodynamic problems of chemistry and heat engineering. The large surface area of films and their small thickness make it possible to accelerate thermal, diffusive, and chemical processes at the gas-liquid boundary.Theoretical studies of liquid flow in a vertical descending thin layer are presented in [1–4]. In this paper we study ascending wave flows of a liquid in a thin vertical layer in contact with a gas, i.e., flows in the direction opposite the action of the force due to gravity, with account for the action of the gas on the liquid surface. Such motions are encountered when oil is extracted from strata that are saturated with gas. At some distance from the stratum the oil and gas separate: the gas travels at high velocity inside the pipe, occupying a considerable portion of the pipe, and the liquid is displaced toward the pipe walls, forming a thin film. In certain cases a wave-like interface develops between the oil and gas that travels with a velocity greater than that of the liquid but less than the average gas velocity. Similar phenomena are observed in high velocity mass exchangers.We examine the effect of the gas for both laminar and turbulent flow.Studies that neglect the effect of the gas flow on the liquid show that for waves on the film surface whose lengths are considerably longer than the average thickness of the layer, the liquid motion in the film is described by boundary layer equations in which account is taken of the mass force, i.e., the force due to gravity. With some approximation, we can assume that in accounting for the effect of the gas on the liquid the liquid flow is described by these same equations.     

Localization of wave motions in a fluid-filled cylinder

The properties of harmonic surface waves in a fluid-filled cylinder made of a compliant material are studied. The wave motions are described by a complete system of dynamic equations of elasticity and the equation of motion of a perfect compressible fluid. An asymptotic analysis of the dispersion equation for large wave numbers and a qualitative analysis of the dispersion spectrum show that there are two surface waves in this waveguide system. The first normal wave forms a Stoneley wave on the inside surface with increase in the wave number. The second normal wave forms a Rayleigh wave on the outside surface. The phase velocities of all the other waves tend to the velocity of the shear wave in the cylinder material. The dispersion, kinematic, and energy characteristics of surface waves are analyzed. It is established how the wave localization processes differ in hard and compliant materials of the cylinder __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 4, pp. 72–86, April 2008.     

Water waves in a suspended container

Analysis and experiments are carried out on a horizontal rectangular wave tank which swings at the lower end of a pendulum. The walls of the tank generate waves which affect the motion of the pendulum. For small displacements of the tank, linearised shallow water equations are used to model the motion, and there exist time-periodic solutions for the system whose periods are governed by a transcendental relation. Numerical and analytic solutions of this relation show that the fundamental period is greater than both the period of the empty tank (moving like a simple pendulum) and the fundamental period of the standing wave which occurs when the tank is removed from its supports and held fixed. For a rectangular tank the theory compares well with some experimental measurements. Qualitative observations are also made of the effect of breaking waves on the tank motion: for a tank which has a mass small compared with its load the energy dissipated by breaking waves can rapidly reduce the amplitude of swing of the tank. Potential flow theory is used with linearised free-surface boundary conditions to find time periodic motions for a tank with a hyperbolic cross section.     

Experiments to study baroclinic waves penetrating into a stratified layer by a quasi-geostrophic potential vorticity equation

The upward penetration of the mid-latitude baroclinic waves has been well known to influence the air motion in the stratosphere. To study this in the laboratory, we present a new type of rotating fluid annulus experiment where a radial temperature difference is imposed on water near the flat bottom to produce a baroclinic wave in a lower layer and the water surface is heated to create a stratified upper layer. To test this model, we carried out experiments to study baroclinic waves penetrating into the stratified layer by a quasi-geostrophic potential vorticity equation in the atmospheric dynamics. The baroclinic waves produced in the lower fluid layer were observed to penetrate into the stratified upper layer. As expected from the equation, the vorticities of the flows for wavenumbers 3, 4 and 5 decayed exponentially with height and showed a weaker decay for lower wavenumbers.     

Effect of Internal Heat Generation on the Onset of Marangoni Convection in a Fluid Layer Overlying a Layer of an Anisotropic Porous Medium

Linear stability analysis has been performed to investigate the effect of internal heat generation on the criterion for the onset of Marangoni convection in a two-layer system comprising an incompressible fluid-saturated anisotropic porous layer over which lies a layer of the same fluid. The upper non-deformable free surface and the lower rigid surface are assumed to be insulated to temperature perturbations. The fluid flow in the porous layer is governed by the modified Darcy equation and the Beavers–Joseph empirical slip condition is employed at the interface between the two layers. The resulting eigenvalue problem is solved exactly. Besides, analytical expression for the critical Marangoni number is also obtained by using regular perturbation technique with wave number as a perturbation parameter. The effect of internal heating in the porous layer alone exhibits more stabilizing effect on the system compared to its presence in both fluid and porous layers and the system is least stable if the internal heating is in fluid layer alone. It is found that an increase in the value of mechanical anisotropy parameter is to hasten the onset of Marangoni convection while an opposite trend is noticed with increasing thermal anisotropy parameter. Besides, the possibilities of controlling (suppress or augment) Marangoni convection is discussed in detail.     

Propagation of Love waves in a functionally graded piezoelectric material (FGPM) layered composite system

In this theoretical study, we investigate the propagation of Love waves in a layered structure consisting of two different homogenous piezoelectric materials, an upper layer and a substrate. A functionally graded piezoelectric material (FGPM) buffer layer is in between the upper layer and the substrate. We employ the power series technique to solve the governing differential equations with variable coefficients. The influence of the gradient coefficients of FGPM and the layer thicknesses on the dispersion relations, the electro-mechanical coupling factor, and the stress distributions of Love waves in this structure are investigated. We demonstrate that the low gradient coefficient raises the significant variation of the phase velocity within a certain range of ratios of upper layer thickness to equivalent thickness. The electro-mechanical coupling factor can be increased when the equivalent thickness equals one or two wavelengths, and the discontinuity of the interlaminar stress can be eliminated by the FGPM buffer layer. The theoretical results set guidelines not only for the design of high-performance surface acoustic wave (SAW) devices using the FGPM buffer layer, but also for the measurement of material properties in such FGPM layered structures using Love waves.     

Surface electro-elastic Love waves in a layered structure with a piezoelectric substrate and two isotropic layers

Propagation of electro-elastic surface Love waves in a structure consisting of a piezoelectric half-space substrate of crystal class 6, 4, 6 mm or 4 mm and two layers, one of which (adjacent to the substrate) is a conducting material and the second is either a conducting or a dielectric material, is considered. The mathematical model obtained includes all the above crystal classes i.e. the surface wave problems related to all these classes are presented in a single mathematical model. The dispersion equation for the existence of Love surface waves with respect to phase velocity is obtained. Numerical calculations are carried out for three different layered structures. The effect of the second layer on the propagation behaviour of the surface Love wave in the structure is revealed.     

Wave propagation along a non-principal direction in a compressible pre-stressed elastic layer

The dispersive behavior of small amplitude waves propagating along a non-principal direction in a pre-stressed, compressible elastic layer is considered. One of the principal axes of stretch is normal to the elastic layer and the direction of propagation makes an angle θ with one of the in-plane principal axes. The dispersion relations which relate wave speed and wavenumber are obtained for both symmetric and anti-symmetric motions by formulating the incremental boundary value problem for a general strain energy function. The behavior of the dispersion curves for symmetric waves is for the most part similar to that of the anti-symmetric waves at the low and high wavenumber limits. At the low wavenumber limit, depending on the pre-stress and propagation angle, it may be possible for both the fundamental mode and the next lowest mode to have finite phase speeds, while other higher modes have an infinite phase speed. At the high wavenumber limit, the phase speeds of the fundamental mode and the higher modes tend to the Rayleigh surface wave speed and the limiting wave speeds of the layer, respectively. Numerical results are presented for a Blatz–Ko material and the effect of the propagation angle is clearly illustrated.