Propagation of Rayleigh waves on curved surfaces

The propagation and properties of Rayleigh waves on curved surfaces are investigated theoretically. The Rayleigh wave dispersion equation for propagation on a curved surface is derived as a parabolic equation, and its penetration depth is analyzed using the curved surface boundary. Reciprocity is introduced to model the diffracted Rayleigh wave beams. Simulations of Rayleigh waves on some canonical curved surfaces are carried out, and the results are used to quantify the influence of curvature. It is found that the velocity of the surface wave increases with greater concave surface curvature, and a Rayleigh wave no longer exists once the surface wave velocity exceeds the bulk shear wave velocity. Moreover, the predicted wave penetration depth indicates that the energy in the Rayleigh wave is transferred to other modes and cannot propagate on convex surfaces with large curvature. A strong directional dependence is observed for the propagation of Rayleigh waves in different directions on surfaces with complex curvatures. Thus, it is important to include dispersion effects when considering Rayleigh wave propagation on curved surfaces.     

Mutual effect of long waves and turbulence on the surface of a liquid

The interaction of a long coherent wave with the turbulence on the surface of a liquid is investigated within the framework of the theory of weak turbulence. A closed system of equations is obtained which consists of the dynamic equation for the coherent wave and equations of kinetic type describing the turbulent subsystem. It is shown that because of the interaction with the turbulent subsystem, coherent waves with wave vectors identical in magnitude but opposite in direction are coupled. The additional attenuation of the coherent wave because of the interaction is estimated; this attenuation may be considerably greater than that caused by molecular viscosity. A change in the spectrum of height correlators of the liquid surface is seen in the presence of a coherent wave.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 1, pp. 100–109, January–February, 1973.     

An experimental study on turbulent coherent structures near a sheared air-water interface

The turbulence structures near a sheared air-water interface were experimentally investigated with the hydrogen bubble visualization technique. Surface shear was imposed by an airflow over the water flow which was kept free from surface waves. Results show that the wind shear has the main influence on coherent structures under air-water interfaces. Low- and high- speed streaks form in the region close to the interface as a result of the imposed shear stress. When a certain airflow velocity is reached, “turbulent spots” appear randomly at low-speed streaks with some characteristics of hairpin vortices. At even higher shear rates, the flow near the interface is dominated primarily by intermittent bursting events. The coherent structures observed near sheared air-water interfaces show qualitative similarities with those occurring in near-wall turbulence. However, a few distinctive phenomena were also observed, including the fluctuating thickness of the instantaneous boundary layer and vertical vortices in bursting processes, which appear to be associated with the characteristics of air-water interfaces. The project supported by the National Natural Science Foundation of China (Grant No.19672070)     

Improvement of the acoustic black hole effect by using energy transfer due to geometric nonlinearity

Acoustic Black Hole effect (ABH) is a passive vibration damping technique without added mass based on flexural waves properties in thin structures with variable thickness. A common implementation is a plate edge where the thickness is locally reduced with a power law profile and covered with a viscoelastic layer. The plate displacement in the small thickness region is large and easily exceeds the plate thickness. This is the origin of geometric nonlinearity which can generate couplings between linear eigenmodes of the structure and induce energy transfer between low and high frequency regimes. This phenomenon may be used to increase the efficiency of the ABH treatment in the low frequency regime where it is usually inefficient. An experimental investigation evidenced that usual ABH implementation gives rise to measurable geometric nonlinearity and typical nonlinear phenomena. In particular, strongly nonlinear regime and wave turbulence are reported. The nonlinear ABH beam is then modeled as a von Kármán plate with variable thickness. The model is solved numerically by using a modal method combined with an energy-conserving time integration scheme. The effects of both the thickness profile and the damping layer are then investigated in order to improve the damping properties of an ABH beam. It is found that a compromise between the two effects can lead to an important gain of efficiency in the low frequency range.     

Modeling non-collinear mixing by distributions of clapping microcracks

Nonlinear scattering by distributions of clapping cracks in a non-collinear wave mixing setting is modeled. Features of the nonlinear response discriminating distributions of clapping cracks from quadratic nonlinear damage are investigated for distributions of cracks that are parallel to each other or randomly oriented. The effective properties of these distributions are recovered extending an existing model that applies to open cracks. The equation of motion is solved using a perturbation approach, and its solutions are evaluated numerically. Their dependence on the amplitude of the incident field is found to be linear, in contrast with the quadratic dependence characterizing quadratic nonlinearity. The spectrum of the scattered field is shown to contain an infinite number of higher harmonics already at the first order of perturbation. Grating-like structures due to the opening and closing of cracks are responsible for adding diffraction peaks to the directivity functions of waves scattered by open cracks. The locations of the most prominent peaks of these functions do not satisfy the selection rules controlling nonlinear scattering by quadratic nonlinearity. Examples of these are given, together with others showing the possibility of using at least one of several discrimination modalities offered by non-collinear wave mixing.     

Timoshenko beam theory: A perspective based on the wave-mechanics approach

This paper reviews Timoshenko beam theory from the point of view of wave mechanics. Vibration of beam structures can be studied in terms of either normal modes or propagating waves. The latter wave approach has two distinct features: first, it gives rise to clear physical understanding of beam vibration; second, it leads to exact methods for vibration analysis of beam structures, especially in the mid-frequency range. In this paper, the work on wave solutions of an infinite Timoshenko beam is first discussed. The work on the splitting effect of spinning on wave solutions is also reviewed. The wave is treated as constitutive components of standing waves (i.e. normal modes), and a discussion on how the wave components formulate various standing waves is presented. Finally, several numerical examples are presented to illustrate the pros and cons of using different wave approaches to tackle vibration analysis of finite-length Timoshenko beams.     

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.     

Velocity field after wave breaking

The velocity field in breaking water waves is considered in this paper. A numerical simulation describes in detail the transition from a primary overturning and consequent rebounding jets into a bore front, where the vorticity in the coherent large‐scale eddy structures devolves into turbulence. Spatial changes in the frequency spectra of the kinetic energy and the enstrophy are associated with the production, transport and dissipation of the Reynolds stress and the various wave and turbulent mixing length scales. Mean velocity fields and the wave and kinetic energy in a surf zone are evaluated. Fourier and wavelet spectral analysis is applied to study both the surface elevation and energy changes, and the distinction that must be made between spilling and plunging breakers is clarified in this paper. Copyright © 2002 John Wiley & Sons, Ltd.     

Weak turbulence of capillary waves

In recent years the theory of weak turbulence, i.e. the stochastic theory of nonlinear waves [I, 9], has been intensively developed. In the theory of weak turbulence nonlinearity of waves is assumed to be small; this enables us, using the hypothesis of the random nature of the phases of individual waves, to obtain the kinetic equation for the mean squares of the wave aplitudes.In many cases of weak turbulence a situation arises where damping is considerable in the region of large wave numbers and is separated from the region where the basic energy of the waves is concentrated (as a result either of pumping or of the initial conditions) with a wide region of transparency. In [3,4] the hypothesis was stated that weak turbulence in these eases is completely analogous to hydrodynamic turbulence for large Reynolds numbers in the sense that in the region of transparency a univelsal spectrum is established which is determined onIy by the flow of energy into the region of large wave numbers. The spectrum of hydrodynamic turbulence Sk - k 5/s was obtained by A. N. Kolmogorov and A. M. Obukhov [5,6] from dimensional considerations. In the case of weak turbulence the spectrum - obtained as an exact solution of the stationary kinetic equation.Below the ease of weak turbulence of capillary waves on the surface of a liquid is considered.A kinetic equation is obtained for capillary waves. It is significant that in this case the basic contribution to interaction is provided by the process of the decomposition of a wave into two and by the process of two waves merging into one.It is shown that the collision term of the kinetic equation vanishes with the solution ek - k 7/4. Arguments are advanced in favor of the fact that this solution can be interpreted as a universal spectrum in the region of transparency.     

湍流池湍流特征研究

在用于激光束传输和控制模拟实验的湍流池中温度脉动测量结果表明,原始温度变化曲线上有明显的斜坡（ｒａｍｐ）结构,温度三阶及五阶结构函数和理论预计一致和大气对流边界层的测量相似,斜坡结构的出现对二阶结构函数的影响较小,加热面和冷却附近有明显的猝发现象出现,温度脉动的峰度系数大于３,表现了较大的间歇性。     

关于海洋结构物水动力学的研究

建造在海洋中的结构物承受波浪、海流、风、地震等作用力,其中尤以波浪力最为重要且复杂。迄今,在小构件海洋结构物波浪力设计计算中,最普遍使用的仍然是Morison方程。但是,这里所说的Morison方程已经不是Morison等的原始方程,而是各种扩展情况下的Morison方程。当初Morison等把由粘性流体得到的阻力项和由理想流体得到的惯性力项相加即得到总波浪力,这种把同一流体既看成是粘性流体又看成是非粘性不可压缩的理想流体的做法实属自相矛盾。但是,建立Morison方程30多年以来,其改进只限于形式稍有      

An Analysis of Turbulent Motions In and Around a Differentially Forced Density Interface

The Rapid-Distortion-Theory-based analysis proposed by Fernando and Hunt [1] is extended to study the nature of turbulence in and around a density interface sandwiched between turbulent layers with dissimilar properties. It is shown that interfacial motions consist of low-frequency, resonantly excited, nonlinear internal waves and high-frequency, linear internal waves driven by background turbulence. Based on the assumptions that (i) all resonant waves and some nonresonant waves having frequencies close to the resonant frequencies grow rapidly, break, and cause interfacial mixing, (ii) the spectral amplitude of the vertical velocity in the wave-breaking regime is constant, and (iii) kinetic energy is equipartitioned between linear and nonlinear breaking wave regimes, the r.m.s. vertical velocity at the interface and the turbulent kinetic energy flux into the interface are calculated. The migration velocity of the interface is calculated using the additional assumption that the buoyancy flux into a given turbulent layer is a fixed fraction of the turbulent kinetic energy flux supplied to the interface by the same layer. The calculations are found to be in good agreement with the entrainment data obtained in previous laboratory experiments in the parameter regime where the interface is dominated by internal wave dynamics. Received 23 July 1997 and accepted 8 January 1999     

壁面加热湍流相干结构的实验研究

用热线风速仪测量了风洞中壁面加热平板湍流边界层不同法向位置的流向脉动速度,用条件采样方法中的m-u-level法和VITA法研究了壁面加热对壁湍流相干结构猝发条件平均波形的影响,发现壁面另热能够提高猝发强度,缩短喷射时间,使一次猝发中出现多次喷射事件的现象明显增多,增强猝发的间歇性。     

Laboratory study on the evolution of waves parameters due to wave breaking in deep water

Understanding of the occurrence of the wave breaking, the process of the wave breaking and evolution of waves after they break in deep water is crucial to simulate the growth of wind wave in ocean. In this study, deep-water breaking waves with various spectral types, center frequencies and frequency bandwidths are generated in a wave flume based on energy focusing theory. The time series of the wave surface elevation along the flume are obtained by 22 wave probes mounted along the central line of the flume. The characteristics of deep-water wave breaking are analyzed using the spectrum analysis based on the Fast Fourier Transform (FFT). For small center frequency the maximum height of wave surface generated using the Pierson–Moskowitz (P–M) spectrum is produced and the impact of the frequency width is small in wave breaking zone. While the spectral type has a significant impact on the local wave steepness during breaking, the influence of center frequency and frequency width on the local wave steepness is very weak. The significant wave steepness changes significantly after wave breaking, but it remains stable in the upstream or the downstream of wave breaking zone. After wave breaking, the peak frequency remains stable, but the spectrally weighted wave frequency changes significantly. The relationship between the level of downshift and the incident wave steepness is approximately linear. By analyzing the energy spectra, it is found that the energy loses near high frequency of controlling frequencies range and increases near peak frequency during the wave breaking. After wave breaking, the total energy dissipates remarkably with increasing breaking intensity.     

Turbulent breaking of overturning gravity waves below a critical level

The interaction of an internal gravity wave with its evolving critical layer and the subsequent generation of turbulence by overturning waves are studied by three-dimensional numerical simulations. The simulation describes the flow of a stably stratified Boussinesq fluid between a bottom wavy surface and a top flat surface, both without friction and adiabatic. The amplitude of the surface wave amounts to about 0.03 of the layer depth. The horizontal flow velocity is negative near the lower surface, positive near the top surface with uniform shear and zero mean value. The bulk Richardson number is one. The flow over the wavy surface induces a standing gravity wave causing a critical layer at mid altitude. After a successful comparison of a two-dimensional version of the model with experimental observations (Thorpe [21]), results obtained with two different models of viscosity are discussed: a direct numerical simulation (DNS) with constant viscosity and a large-eddy simulation (LES) where the subgrid scales are modelled by a stability-dependent first-order closure. Both simulations are similar in the build-up of a primary overturning roll and show the expected early stage of the interaction between wave and critical level. Afterwards, the flows become nonlinear and evolve differently in both cases: the flow structure in the DNS consists of coherent smaller-scale secondary rolls with increasing vertical depth. On the other hand, in the LES the convectively unstable primary roll collapses into three-dimensional turbulence. The results show that convectively overturning regions are always formed but the details of breaking and the resulting structure of the mixed layer depend on the effective Reynolds number of the flow. With sufficient viscous damping, three-dimensional turbulent convective instabilities are more easily suppressed than two-dimensional laminar overturning.     

Visualisation and modelling studies of churn flow in a vertical pipe

Churn flow is an important flow regime intermediate between slug flow and annular flow. A feature of this regime is the occurrence of very large waves travelling upwards over a liquid film substrate which may intermittently travel downwards. These waves are often formed close to the liquid inlet where their behaviour is usually difficult to observe. This paper describes a series of experiments using a test section with a specially constructed transparent liquid inlet. High-speed video recordings show clearly the process of wave formation and analysis of the recordings gave data on wave frequencies and typical velocities. Also, predictions of velocity and distance travelled by the waves were obtained via the application of a simple mathematical model that takes into account the forces acting on a circumferentially coherent wave.     

On the analysis of 2D nonlinear gravity waves with a fully nonlinear numerical model

A fully nonlinear numerical method, developed on the basis of Euler equations, is used to study the dynamics of nonlinear gravity waves, mainly in the aspects of the propagation of Stokes wave with disturbed sidebands, the evolution of one wave packet and the interaction of two wave groups. These cases have previously been studied with the higher order spectral method, which will be an approximately fully nonlinear scheme if the order of nonlinearity is not large enough, while the present method in the case of the 2D model has an integration scheme that is exact to the computer precision. As expected, in most cases the results are consistent between these two numerical models and it is confirmed again that this fully nonlinear numerical model is also capable of maintaining a high accuracy and good convergence, particularly in the long-term evolutionary process.