Regular wave integral approach to numerical simulation of radiation and diffraction of surface waves

A regular wave integral method is developed in the discretisation of a linear hydrodynamic problem on radiation and diffraction of surface waves by a floating or submerged body. The velocity potential of the problem is expressed as a solution of a body boundary integral equation involving the pulsating free surface Green function or pulsating free surface sources distributed on the body surface. With the use of a discretisation on the regular wave integral rather than discretisations on the singular wave integral of the Green function as in earlier investigations, the singular wave integral is approximated as an expansion of regular (or nonirregular) wave potentials. Influence coefficients between pulsating free surface source points are computed by the approximate expansion together with Hess–Smith panel integral formulas. Thus the velocity potential solution is evaluated by a boundary element algorithm. The numerical results produced from the proposed method agree well with semi-analytic solution results.     

Structure of the drift flow initiated by periodic waves traveling over a viscous fluid surface

An analytical procedure used in calculating the Stokes drift velocity (the drift motion initiated by the propagation of a capillary-gravity wave over an ideal fluid surface) is applied to the problem of the calculation of an analogous drift flow in a viscous fluid. An expression for the velocity of the Stokes drift modified with allowance for viscosity is constructed. The properties and the role of the modified Stokes drift in the general pattern of the drift in a viscous fluid are analyzed.     

Stability of the wave motion on the charged interface between two immiscible fluids in the presence of a tangential velocity field discontinuity

In the second-order approximation in the dimensionless wave amplitude, the problem of nonlinear periodic capillary-gravity wave motion of the uniformly charged interface between two immiscible ideal incompressible fluids, the lower of which is perfectly electroconductive and the upper, dielectric, moves translationally at a constant velocity parallel to the interface, is solved analytically. It is shown that on the uniformly charged surface of an electroconductive ideal incompressible fluid the positions of internal nonlinear degenerate resonances depend of the medium density ratio but are independent of the upper medium velocity and the surface charge density on the interface. All resonances are realized at densities of the upper medium smaller than the density of the lower medium. In the region of Rayleigh-Taylor instability with respect to density there is no resonant wave interaction.     

Scattering of water waves by a submerged thin vertical elastic plate

The problem of water wave scattering by a thin vertical elastic plate submerged in infinitely deep water is investigated here assuming linear theory. The boundary condition on the elastic plate is derived from the Bernoulli–Euler equation of motion satisfied by the plate. This is converted into the condition that the normal velocity of the plate is prescribed in terms of an integral involving the difference in velocity potentials (unknown) across the plate multiplied by an appropriate Green’s function. The reflection and transmission coefficients are obtained in terms of integrals involving combinations of the unknown velocity potential on the two sides of the plate and its normal derivative on the plate, which satisfy three simultaneous integral equations, solved numerically. These coefficients are computed numerically for various values of different parameters and are depicted graphically against the wave number for different situations. The energy identity relating these coefficients is also derived analytically by employing Green’s integral theorem. Results for a rigid plate are recovered when the parameters characterizing the elastic plate are chosen negligibly small.     

三层流体中斜入射波作用下半无限板的水弹性响应

解析地研究了在三层流体中斜入射波浪与半无限弹性板的相互作用引起的波散射和板的水弹性响应. 三层流体在界面处的密度发生阶跃, 各层为一常数. 假设流体不可压缩、无黏、流体运动无旋. 在线性势流理论框架下, 使用本征函数展开法和内积式给出波板相互作用的半解析解. 根据色散关系分析, 得到了表面波模态和界面波模态入射时的临界入射角. 随着物理参数的变化, 临界角将随之发生变化. 临界角决定了当由开阔水域向板覆盖水域传播的表面波或界面波的存在性: (1)板覆盖水域入射界面上, 透射波能否存在; (2)入射界面之上界面中, 板覆盖水域中的透射波以及开阔水域中的反射波能否存在. 当下界面波入射时并且入射角足够大时, 开阔水域中的下界面波模态是整个流体域中唯一存在的模态.      

Flow control over an undulating membrane

The flow along a flexible membrane forced to undulate in the form of a streamwise travelling wave pattern is studied experimentally in detail. Flow field and force measurements confirm that the form drag of the wavy wall is significantly reduced when starting the undulatory motion. A mechanical model of an undulating membrane was built, based on previous investigations described in literature, and placed in an open water channel. The motion pattern of the membrane was prescribed in such a way to achieve a downstream travelling wave with increasing amplitude. The exploratory focus laid on the identification of hydrodynamic mechanisms of drag reduction due to undulatory motion. The wave-speed c of the travelling wave was set proportional to the incoming flow velocity U, according to an optimum ratio identified by previous numerical and experimental investigations. Poisson’s equation for the pressure was used to calculate the 2D pressure field from the experimental data of the unsteady flow field. In addition, the integral drag force of the membrane, as a function of c/U, was measured with a force balance to compare with previous published numerical findings. Furthermore, the velocities close to the surface of the membrane were measured, and the boundary layer profiles were determined. The resulting normalised velocity profiles affirm an oscillation between laminar and turbulent flow over one period of the motion. The results are in good agreement with previous experimental and numerical findings. Additionally, the characteristics of the flow along a travelling wave with increasing amplitude are discussed in more detail.     

考虑表面张力的非传播孤立波

本文考虑了表面张力,用多重尺度法导出了与立方 Schrodinger 方程相类似的非传播孤立波的基本方程,得到了非传播孤立波解。用毛细重力波理论解释了非传播孤立波横向谐振中波峰尖、波谷平的原因。在σ～kh 平面上首次给出了可产生非传播孤立波的二个参数区,但现有的实验点都在区域(1)中。     

Measurements of velocity,temperature and velocity fluctuation distributions in falling liquid films

The velocity, temperature and velocity fluctuation distributions within falling spindle oil films in an inclined rectangular channel were measured using hot-wire techniques and thin thermocouples. The interfacial shear was caused by cocurrent air flow.The results indicate that the liquid films are as a whole much more laminar-like than turbulent in a range of Reynolds numbers (4γ/μ) up to the experimental limit of 6000. Mixing motion occurs in the vicinity of the interface; however, the flow near the wall surface exhibits no sign of such eddy motions, as predicted by the wall law for single phase turbulent flow. Although velocity fluctuation is observed within films with interfacial shear, mean velocity profiles are approximately the same as those obtained by the laminar film prediction.     

On the propagation of energy in linear conservative waves

Summary This paper is concerned with the question when and why the rate of energy propagation in a system of waves equals the group velocity. It is shown by the method of stationary phase that this equality holds, for travelling waves without dissipation, whenever this method applies. The reason why this result can be obtained by this kinematical method is investigated by a discussion of simple harmonic waves. It is shown that the choice of an expression for the energy density to be used in connection with a given wave equation is restricted by the conservation of energy in such a way that the average rate of work done divided by the average energy density always equals the group velocity. Finally some examples of wave motion are discussed to illustrate the derived formulae.     

Linear resonance in viscous films on inclined wavy planes

We study viscous gravity-driven films flowing over periodically undulated substrates. Linear analysis describes steady flow along small amplitude corrugations for films of arbitrary thickness. Solving the resulting system numerically, we demonstrate resonance (or, possibly, near resonance) and identify different behaviours for thin, intermediate and thick films. Approximating the leading-order velocity profile by the free surface value allows for an analytic solution, which – in the limit of high Reynolds numbers – recovers the different regimes and reveals the relevant physical mechanisms. Our results support the view that the resonance is associated with an interaction of the undulated film with capillary-gravity waves travelling against the mean flow direction. As a consequence, the resonance peak is attained under conditions that render the wave phase velocity equal to zero in the laboratory reference frame, and thus permit direct exchange of energy between the steadily deformed film and the free surface.     

In-plane forced vibration of curved pipe conveying fluid by Green function method

The Green function method(GFM) is utilized to analyze the in-plane forced vibration of curved pipe conveying fluid, where the randomicity and distribution of the external excitation and the added mass and damping ratio are considered. The Laplace transform is used, and the Green functions with various boundary conditions are obtained subsequently. Numerical calculations are performed to validate the present solutions, and the effects of some key parameters on both tangential and radial displacements are further investigated. The forced vibration problems with linear and nonlinear motion constraints are also discussed briefly. The method can be radiated to study other forms of forced vibration problems related with pipes or more extensive issues.     

Non-propagating solitary waves with the consideration of surface tension

In this paper, the governing equation for the non-propagating solitary waves, similar to the cubic Schrödinger equation, is derived by the multiple scales with the consideration of surface tension. The non-propagating solitary wave solution is given. It is explained by the capillary-gravity wave theory that the crests are sharpened and the troughs are flattened in the transversal harmonic of the non-propagating solitary waves. On σ～kh plane, two parameter regions are obtained in which the non-propagating solitary wave can occur, but all existing experimental parameters are in region 1 (Fig. 1).     

Nonlinear Capillary-Gravity Waves on the Charged Surface of an Ideal Fluid

A second-order asymptotic expression for the profile of a capillary-gravity wave traveling over the charged surface of an ideal incompressible fluid is calculated analytically. Two types of steady-state profiles of nonlinear periodic capillary-gravity waves are found. For a certain fixed dimensionless surface charge the shape of the tops of the nonlinear waves changes: from blunt to pointed for short waves and from pointed to blunt for long waves.     

Numerical Simulation of Three-Dimensional Bubble Oscillations by a Generalized Vortex Method

A numerical method is implemented for simulating the simultaneous three-dimensional volume and shape oscillations of a compressible vapor or gas bubble suspended in an inviscid ambient fluid in the presence of interfacial tension. The flow generated by the bubble expansion, contraction, and deformation is represented by an interfacial distribution of potential dipoles supplemented by a point source situated inside the bubble, accounting for changes in the bubble volume. The mathematical formulation is completed by setting the strength of the point source proportional to the integral of the density of the double-layer potential over the interface. The motion of marker points distributed over the interface is computed using a boundary-element implementation of Baker's generalized vortex method in which the normal component of the interfacial velocity is computed in terms of tangential derivatives of the vector potential associated with the dipoles, whereas the tangential component of the interfacial velocity is computed in terms of the surface gradient of the scalar harmonic potential. The density of the double-layer distribution is computed by solving an integral equation of the second kind using an iterative method, while the evolution of the interfacial distribution of the harmonic potential is computed using Bernoulli's equation for irrotational flow. The onset of interfacial irregularities due to numerical instabilities is prevented by truncating the Fourier–Legendre spectrum of the interfacial distribution of the harmonic potential. With smoothing implemented, the numerical method is capable of describing simultaneous volume and shape oscillations for an indefinite period of time. Received 7 September 2001 and accepted 30 April 2002 Published online 30 October 2002 RID="*" ID="*" This research was supported by a grant provided by NASA. Communicated by J.R. Blake     

Impulse wave run-over: experimental benchmark study for numerical modelling

This research intends to provide a detailed data basis for numerical modelling of impulse waves. Three tests are described involving a rectangular wave channel, in which a trapezoidal ‘breakwater’ was inserted to study wave run-over. In addition, a reference test is also described, in which the breakwater was removed. Two-dimensional impulse waves were generated by means of subaerial granular slides accelerated by a pneumatic landslide generator into the water body. Wave propagation and run-over over the artificial breakwater are documented by a set of high-quality photographs. Water surface profiles were recorded using capacitance wave gages upstream and downstream of the breakwater, and velocity vector fields were determined for the run-over zone by means of Particle Image Velocimetry. The measurements are compared with predictive formulae for wave features and wave non-linearity. The present data set involves both simple channel topography and wave features to allow for numerical simulations under basic laboratory conditions.     

二层流体中波动问题的Hamilton正则方程

研究了两种常密度不可压缩理想流体组成的垂直分层的二流体系统的无旋等熵流动,考虑了上层流体与空气及两层流体间的表面张力。流动区域在水平方向无限伸展,上层流体有限深度,下层流体无限深。利用自由面及分界面相对于静止时平衡位置的偏移以及两层流体的速度势构造了Hamilton函数。为导出Hamilton正则方程引用了Euler描述下的流体运动的变分原理。自由面的位移是Hamilton意义下的正则变量,其对偶变量是上层流体在自由面上取值的速度势与密度的乘积。另一个正则变量是分界面的位移,其对偶变量是下层流体的密度与下层流体速度势在分界面上所取值的乘积减去上层流体密度与上层流体速度势在分界面上所取值的相应乘积。导出的Hamilton结构对分析分层流动中表面波与内波的相互作用是重要的。     

Generalized Fourier transforms for the acoustic pressure and velocity in compressible dissipative stratified media

Generalized Fourier transforms are derived for the acoustic pressure and velocity in compressible dissipative plane stratified media. The acoustic source terms are accounted for in the equations of continuity and force. The acoustic pressure and velocity are each expressed as sums of two infinite (branch cut) integrals and a discrete term. In the far field the infinite integrals correspond to the direct and specularly reflected waves and the lateral wave. The discrete term associated with the pole of the reflection coefficient is the surface wave. The transforms provide a suitable basis for the expansion of the acoustic pressure and the velocity when the height of the interface, the adiabatic bulk modulus, the equilibrium density, and the absorption of the medium vary. Both exact boundary conditions and the approximate impedance boundary condition are considered in this work.     

An experimental investigation on the effects of surface gravity waves on the water evaporation rate in different air flow regimes

Estimating rate of evaporation from undisturbed water surfaces to moving and quiet air has been the topic a vast number of research activities. The obvious presence of various shapes of gravity waves on the water body surfaces was the motivation of this experimental investigation. In this investigation experimental measurements have been done to quantify evaporation rate from wavy water surfaces in free, mixed and forced convection regimes. The effects of a wide range of surface gravity waves from low steepness, round shaped crest with slow celerity, to steep and very slight spilling crest waves, on the water evaporation rate have been investigated. A wide range of ${\text{Gr}}/{\text{Re}}^{2} (0.01 \le {\text{Gr}}/{\text{Re}}^{2} \le 100)$ was achieved by applying different air flow velocities on a large heated wave flume equipped with a wind tunnel. Results reveal that wave motion on the water surface increase the rate of evaporation for all air flow regimes. For free convection, due to the effect of wave motion for pumping rotational airflows at the wave troughs and the dominant effect of natural convection for the air flow advection, the maximum evaporation increment percentage from wavy water surface is about 70 %. For mixed and forced convection, water evaporation rate increment is more sensitive to the air flow velocity for the appearance of very slight spilling on the steep wave crests and the leeward air flow structures.     

Wave propagation through mangrove forests in the presence of a viscoelastic bed

The influence of viscoelastic ocean beds on the characteristics of surface waves passing through mangrove forests is analyzed under the assumption of linearized water wave theory in two dimensions. The trunks of the mangroves are assumed to be in the upper-layer inviscid fluid domain, whilst the roots are inside the viscoelastic bed. The associated equation of motion is obtained by coupling the Voigt’s model for flow within the viscoelastic medium with the equation of motion in the presence of mangroves. The modified dynamic conditions are coupled with the kinematic conditions to obtain the boundary condition at the free surface and the interface of the two fluids consisting of the upper layer inviscid fluid and the viscoelastic fluid bed. To understand the effects of bed viscosity as well as elasticity on energy dissipation, the complex dispersion relation associated with the plane progressive wave is derived and analyzed. Effect of physical parameters associated with mangroves and viscoelastic bed on wave motion in surface and internal modes are computed and analyzed to understand their roles in attenuating wave effects. The present model will be useful in the better understanding of wave propagation through mangroves in the coastal zone having muddy seabed.     

Long wavelength approximation to peristaltic motion of a power law fluid

Peristaltic motion of a power law fluid in a two-dimensional channel is studied. Assuming that the wavelength of the peristaltic wave is large in comparison to the mean half-width of the channel, a solution for the stream function is obtained as an asymptotic expansion in terms of slope parameter. Expressions for axial pressure gradient and shear stress are derived. The effect of flow behaviour indexn on the streamline pattern and shear stress is studied and the phenomenon of trapping is discussed.