Diffraction of internal waves by a circular cylinder near the pycnocline

Propagation of internal waves over a circular cylinder under the conditions of a continuous stratification characterized by the presence of a high-gradient density layer (the pycnocline) of finite thickness is studied. The dependences of the coefficent of wave propagation on the wavelength of the first-mode incident wave for various thicknesses of the pycnocline are obtained. In the diffraction of internal waves, substantial nonlinear effects are shown to occur, which result in the appearance of waves of double oscillation frequency compared to the frequency of the incident waves. The generation coefficient for these waves is found. Lavrent’ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 2, pp. 79–85, March–April, 1999.     

Experimental Study of the Dynamic Effect of an Internal Solitary Wave on a Submerged Circular Cylinder

The hydrodynamic loads due to the interaction of a small-amplitude internal solitary wave with a submerged circular cylinder in a two-layer system of miscible fluids were studied experimentally. The dependence of the internal-wave transmission coefficient on the position of the center of the cylinder relative to the pycnocline and on the ratio of the cylinder diameter to the fluid-layer thickness was obtained. The effects of the pycnocline thickness and the depth of the center of the cylinder on the value of the hydrodynamic loads were studied. Visualization of the flow structure was performed. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 46, No. 6, pp. 36–44, November–December, 2005.     

Laboratory experiments on internal wave interactions with a pycnocline

Laboratory experiments have been performed to investigate the interaction of internal waves with a pycnocline. An oscillating cylinder generated internal wave beams, which were observed using the synthetic schlieren technique. Internal waves incident on the pycnocline layer excited higher-frequency modes. In the absence of shear, a discrete spectrum of harmonic modes was generated due to nonlinear effects. These harmonic modes might play a role in the formation of internal solitary waves which have been observed in ocean pycnoclines. With shear, a continuous spectrum of excited modes was found.     

内孤立波作用下潜深对潜体运动响应和载荷特性的影响研究

内孤立波常发生于海洋密度跃层, 因其峰高谷深、携带能量巨大, 在传播过程中会导致跃层上下的海水流动呈现剪切状态, 并引起突发性的强流. 潜体在水下悬停时极有可能会遭遇内孤立波, 由于内孤立波的流场特性, 置于跃层上下的悬浮潜体所产生运动响应和水动力载荷变化不尽相同, 甚者会出现掉深现象. 为探究潜深对波体耦合作用的影响, 基于不可压缩N-S方程和mKdV理论, 采用速度入口造波, 结合重叠网格技术和流固耦合方法, 建立了分层流中内孤立波耦合水下潜体多自由度运动的数值模型, 通过该模型分析了不同潜深下悬浮潜体的运动响应和载荷特性. 结果表明: 在内孤立波作用下, 位于密度跃层上方和跃层中的潜体顺着波的前进方向运动, 先下沉后抬升, 位于跃层下方的潜体则会逆流持续下沉; 潜体与波面的垂向距离越小, 对其纵荡、垂荡和速度的影响越显著, 而位于密度跃层中的潜体在分界面处沿着波形运动, 其运动响应和载荷变化受影响较小; 潜体在跃层上、下流体中所受水平力的方向相反, 水平力峰值小于垂向力峰值, 且位于跃层下方的潜体一直受到低头力矩, 最终导致掉深.      

Flow pattern transition in liquid-liquid flows with a transverse cylinder

The effect of a cylindrical bluff body on the interface characteristics of stratified two-phase, oil-water, pipe flows is experimentally investigated with high speed Particle Image Velocimetry (PIV). The motivation was to study the feasibility of flow pattern map actuation by using a transverse cylinder immersed in water in the stratified pattern, and particularly the transition from separated to dispersed flows. The cylinder has a diameter of 5 mm and is located at 6.75 mm from the bottom of the pipe in a 37 mm ID acrylic test section. Velocity profiles were obtained in the middle plane of the pipe. For reference, single phase flows were also investigated for Reynolds numbers from 1550 to 3488. It was found that the flow behind the cylinder was similar to the two dimensional cases, while the presence of the lower pipe wall diverted the vorticity layers towards the top. In two-phase flows, the Froude number (from 1.4 to 1.8) and the depth of the cylinder submergence below the interface affected the generation of waves. For high Froude numbers and low depths of submergence the counter rotating von Karman vortices generated by the cylinder interacted with the interface. In this case, the vorticity clusters from the top of the cylinder were seen to attach at the wave crests. At high depths of submergence, a jet like flow appeared between the top of the cylinder and the interface. High speed imaging revealed that the presence of the cylinder reduced to lower mixture velocities the transition from separated to dual continuous flows where drops of one phase appear into the other.     

Hydrodynamics of flow over a transversely oscillating circular cylinder beneath a free surface

In this paper, hydrodynamic force coefficients and wake vortex structures of uniform flow over a transversely oscillating circular cylinder beneath a free surface were numerically investigated by an adaptive Cartesian cut-cell/level-set method. At a fixed Reynolds number, 100, a series of simulations covering three Froude numbers, two submergence depths, and three oscillation amplitudes were performed over a wide range of oscillation frequency. Results show that, for a deeply submerged cylinder with sufficiently large oscillation amplitudes, both the lift amplitude jump and the lift phase sharp drop exist, not accompanied by significant changes of vortex shedding timing. The near-cylinder vortex structure changes when the lift amplitude jump occurs. For a cylinder oscillating beneath a free surface, larger oscillation amplitude or submergence depth causes higher time-averaged drag for frequency ratio (=oscillation frequency/natural vortex shedding frequency) greater than 1.25. All near-free-surface cases exhibit negative time-averaged lift the magnitude of which increases with decreasing submergence depth. In contrast to a deeply submerged cylinder, occurrences of beating in the temporal variation of lift are fewer for a cylinder oscillating beneath a free surface, especially for small submergence depth. For the highest Froude number investigated, the lift frequency is locked to the cylinder oscillation frequency for frequency ratios higher than one. The vortex shedding mode tends to be double-row for deep and single-row for shallow submergence. Proximity to the free surface would change or destroy the near-cylinder vortex structure characteristic of deep-submergence cases. The lift amplitude jump is smoother for smaller submergence depth. Similar to deep-submergence cases, the vortex shedding frequency is not necessarily the same as the primary-mode frequency of the lift coefficient. The frequency of the induced free surface wave is exactly the cylinder oscillation frequency. The trends of wave length variation with the Froude number and frequency ratio agree with those predicted by the linear theory of small-amplitude free surface waves.     

Internal wave radiation by a turbulent fountain in a stratified fluid

Large eddy simulation is applied to model a fountain in a density-stratified fluid. The fountain is formed, as a vertical turbulent jet penetrates through a pycnocline. The jet flow is initiated by the formulation of a boundary condition in the form of an upward neutral-buoyancy fluid flow with the Gaussian axisymmetric mean-velocity profile and a given fluctuation level. It is shown that at a Froude number Fr higher than a certain critical value the fountain executes self-oscillations accompanied by internal wave generation within the pycnocline. The predominant self-oscillation mode is axisymmetric, when the fountain top periodically breaks down generating internal wave packets traveling toward the periphery of the computation domain. The characteristic frequency of the internal waves coincides with that of the fountain top oscillations and monotonically decreases with increase in Fr. The Fr-dependence of the fountain top oscillation amplitude obtained in the numerical solution is in good agreement with the predictions of the theoretical Landau model for the instability mode in the soft self-excitation regime.     

Recurrence and Reconnection Effect in an Attached Two-Dimensional Internal Wave Field

The shadow pattern of the wake behind a horizontal cylinder uniformly towed in a linearly stratified fluid is analyzed. In the “imbedded turbulent vortex” and “narrow turbulent wake” regimes the structure of the attached internal waves with fairly large amplitudes in the neighborhood of the body satisfactorily agrees with calculations based on the linear theory. As the distance from the source increases, the pattern becomes more complex. The internal waves radiated by the body itself and by the isolated vortices in the wake interfere and nonlinearly interact with each other and with the other flow elements. At large ages τ=t/T_b > 5.8 — time is normalized with respect to the buoyancy period) recurrence (spatial structure restored) and restructuring (phase surfaces reconnected) of a group of longest internal waves (of the attached type) are observed. The recurrent waves penetrate through the wake. These effects are associated with a general restructuring of the stratification modulated by a particular wave scale.     

Characterizing texture in polycrystalline materials by ultrasonic waves

A method for characterizing texture from measurements of ultrasonic wave velocities is proposed. In polycrystalline aggregates, ultrasonic wave velocities are strongly affected by orientation distribution coefficients (ODCs), which are usually used to describe the degree of preferred grain orientation in textured materials. In this work, velocities of longitudinal and transverse waves propagating into aluminum alloy 6061 were measured under pure shear, simple shear and uniaxial tension. From the measured ultrasonic wave velocities, the ODCs W₄₀₀ and W₄₂₀ were calculated to infer the deformation-induced texture. The predicted pole figures, obtained using ultrasonic velocities, were in good qualitative agreement with the finite element polycrystal model analyzed pole figures.     

Boundary element modeling for defect characterization potential in a wave guide

Wave scattering analysis implemented by boundary element methods (BEM) and the normal mode expansion technique is used to study the sizing potential of two-dimensional shaped defects in a wave guide. Surface breaking half-elliptical shaped defects of three opening lengths (0.3, 6.35 and 12.7 mm) and through-wall depths of 10–90% on a 10 mm thick steel plate were considered. The reflection and transmission coefficients of both Lamb and shear horizontal (SH) waves over a frequency range 0.05–2 MHz were studied. A powerfully practical result was obtained whereby the numerical results for the S₀ mode Lamb wave and n₀ mode SH wave at low frequencies showed a monotonic increase in signal amplitude with an increase in the defect through-wall depth. At high frequency (usually above the cut-off frequency of the A₁ mode for Lamb waves and the n₁ mode for SH waves, respectively), the monotonic trend does not hold in general due to the energy redistribution to the higher order wave modes. Guided waves impinging onto an internal stringer-like an inclusion were also studied. Both the Lamb and SH waves were shown to be insensitive to the stringer internal inclusions at low frequency. Experiments with piezoelectric Lamb wave transducers and non-contact SH wave electro-magnetic acoustic transducers (EMAT) verified some of the theoretical results.     

Internal wave generation by a fountain in a stratified fluid

The purpose of the study is the direct numerical and theoretical modeling of fountain dynamics in a fluid with density stratification in the form of a pycnocline. The fountain is formed as a vertical jet penetrates through the pycnocline. In numerical simulation the jet flow is initiated by means of preassigning a boundary condition in the form of an upward-directed laminar flow of a neutral-buoyancy fluid with an axisymmetric Gaussian velocity profile. The calculations show that at a Froude number Fr greater than a certain critical value the flow becomes unstable and the fountain executes self-oscillations accompanied by internal wave generation in the pycnocline. Depending on Fr, two self-oscillation modes can be distinguished. At fairly low Fr the fountain executes circular motion in the horizontal plane, in the vicinity of the center of jet, its shape remaining almost invariant. In this case, internal waves in the form of unwinding spirals are radiated. At fairly high Fr another mode predominates, when the fountain top chaotically “strays” in the vicinity of the center of the jet and, periodically breaking down, generates wave packets propagating toward the periphery of the computation domain. In both cases, the main peak in the frequency spectrum of the internal waves coincides with the fountain top oscillation frequency which monotonically decreases with increase in Fr. In numerical simulation the Fr-dependence of the fountain top oscillation amplitude is in good agreement with that predicted by the theoretical model of the concurrence of the interacting modes in the soft self-excitation regime.     

Transient thermoelastic waves in an anisotropic hollow cylinder due to localized heating

This paper presents a theoretical study of transient ultrasonic guided waves generated by concentrated heating of the outer surface of an infinite anisotropic hollow circular cylinder. Generalized thermoelastic theory proposed by Lord and Shulman is adopted to model the dynamic thermoelastic behavior of the cylinder. The concentrated heat source model used is to represent heating due to a pulsed laser beam, which is focused on the outer surface of the cylinder. A semi-analytical finite element (SAFE) method is employed to evaluate guided wave modes in the cylinder. Using integral transform techniques, the modal wave forms are obtained in frequency and wave number domains. Time histories of the propagating modes are then calculated by applying inverse Fourier transformation in the time domain. Numerical results showing the dispersion curves for the group velocities of the propagating modes and transient radial displacements are presented. For this purpose it is assumed that the cylinder is made of transversely isotropic silicon nitride (Si₃N₄). Attention is focused on the propagation characteristics of longitudinal and flexural modes separately.     

Flow pattern behind a submerged body near a pycnocline

A flow past an ellipsoid immersed in a flow of a viscous stratified fluid is studied using a Large Eddy Simulation (LES) method for different locations of the body relative to a density discontinuity. It is shown that, even when the internal Froude number is large, for small angles of attack of the body the stratification affects its drag force. When the body is located above the pycnocline, the presence of the discontinuity results in the ascending of a vortex filament. A spectral analysis showed that the internal waves in the body wake have a multimode structure.     

Two-dimensional unstructured elastic model for acoustic pulse scattering at solid-liquid interfaces

Abstract. A numerical model to simulate elastic waves and acoustic scattering in two spatial dimensions has been developed and thoroughly tested. The model universally includes elastic solids and liquids. The equations of motion are written in terms of stresses, displacements and displacement velocities for control volumes constructed about the nodes of a triangular unstructured grid. The latter conveniently supports various geometries with complex external and internal boundaries separating sub-domains of different elastic properties. Theoretical dispersion for zero mode symmetric () and antisymmetric () waves in a plate has been reproduced numerically with high accuracy, thus verifying the method and code. Comparison of simulated acoustic pulse scattering at water-immersed steel plate with the respective experiments reveals a very good agreement in such delicate features as excitation of the surface (A) wave. The numerical results explain the peculiar location of the surface wave relative to the other ones in experimental registrations. Examples of acoustic pulse interactions with curvilinear metallic shells in water demonstrate flexibility of the method with respect to complex geometries. Potential applications as well as some directions for further improvement to the technique are briefly discussed. Received 5 September 2002 / Accepted 25 November 2002 Published online 4 February 2003 RID="*" ID="*"Permanent address: Ioffe Physical-Technical Institute, 26 Polytekhnicheskaya, 194021 St. Petersburg, Russia Correspondence to: P. Voinovich (e-mail: vpeter@scc.ioffe.ru)     

Vortex shedding around a near-wall circular cylinder induced by a solitary wave

This study developed a two-dimensional generalized vortex method to analyze the shedding of vortices and the hydrodynamic forces resulting from a solitary wave passing over a submerged circular cylinder placed near a flat seabed. Numerical results for validation are compared with other numerical and experimental results, and satisfactory agreement is found. A series of simulations were performed to study the effects of gap-to-diameter ratio and incident wave height on vorticity pattern as well as the forces exerted on the cylinder. The range of the heights of incident waves is from 0.3h to 0.7h, where h is the still water depth. The range of the gap-to-diameter ratios is from 0.1 to 0.8. The results indicate that the flow pattern and the pressure distribution change significantly because of the close proximity of the seabed where the vorticity flux on the seabed-side surface of the cylinder is suppressed. Placing the cylinder nearer the seabed increases the drag and the positive lift on the cylinder. When the gap-to-diameter ratio increases, the pattern of vortices changes because of the interaction between the main recirculation zone and the shear layers separated from the gap. The maxima of drag, lift and total force increase linearly with the height of the incident wave.     

Analysis of the Propagation of Sound Waves in Partially Saturated Soils by Means of a Macroscopic Linear Poroelastic Model

In this paper, the propagation of sound waves in partially saturated soils is investigated. A macroscopic linear model that is based on the two-component model of Biot and on the Simple Mixture Model by Wilmanski is used. For the construction of the model by a micro-macro transition, see Albers, Géotechnique, 2007. We investigate a porous medium consisting of a deformable skeleton and two compressible, chemically non-reacting, pore fluids (liquid and gas). The wave analysis of the poroelastic model reveals the number of acoustic waves and the dependence of velocities and attenuations of these waves on the initial saturation and frequency. There appear four body waves: three longitudinal waves, P1, P2, P3, and one shear wave, S. The P2-wave shows a similar feature as in air–water mixtures: from the theory of suspensions, it is well known that the existence of air bubbles in water reveals a minimum in the sonic velocity. This is also the case for the P2 -speed in the unsaturated porous medium. The P1-velocity increases very abruptly for a certain degree of saturation. This provides the hope for the development of a nondestructive testing method.     

The shock-wave structure in a two-velocity mixture of compressible media with different pressures

The problem of the shock-wave structure in a mixture of two compressible media with different velocities and pressures of components is considered. The problem is reduced to solving a boundary-value problem for two ordinary differential equations that describe the velocity relaxation and pressure equalization of the components. Using methods of the qualitative theory of dynamic systems on a plane, the existence and uniqueness of four types of waves are shown: (a) fully dispersed waves; (b) frozen-dispersed waves; (c) dispersed-frozen waves; (d) frozen waves of two-front configuration. A chart of solutions of the corresponding flow types is constructed in the plane of the following parameters: the initial velocity of the mixture and the initial volume concentration of one of the components. The numerical calculations conducted illustrate the obtained analytical structures of the shock wave. It is shown that the results obtained using the suggested mathematical model are in agreement with experimental data on the dependence of the velocity of the dispersed shock wave on the equilibrium pressure behind the shock-wave front for a mixture of silica sand and water. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 2, pp. 10–19, March–April, 1998.     

Three-dimensional green water velocity on a model structure

Flow kinematics of green water due to plunging breaking waves impinging on a simplified, 3D model structure was investigated in the laboratory. Two breaking wave conditions were tested: one with waves impinging on the vertical wall of the model at still water level, and the other with waves impinging on the horizontal deck surface. The bubble image velocimetry (BIV) technique was used to measure flow velocities. Measurements were taken on both vertical and horizontal planes. Evolution of green water flow kinematics in time and space was revealed and was found to be quite different between the two wave conditions, even though the incoming waves are essentially identical. The time history of maximum velocity is demonstrated and compared. In both cases, the maximum velocity occurs near the green water front and beneath the free surface. The maximum horizontal velocity for the deck impinging case is 1.44C with C being the wave phase speed, which is greater than 1.24C for the wall impingement case. The overall turbulence level is about 0.3 of the corresponding maximum velocity in each wave condition. The results were also compared with 2D experimental results to examine the 3D effect. It was found that the magnitude of the maximum vertical velocity during the runup process is 1.7C in the 3D model study and 2.9C in the 2D model study, whereas the maximum horizontal velocity on the deck is similar, 1.2C in both 3D and 2D model studies.     

Free v/s forced vibrations of a cylinder at low Re

Free vibrations of a circular cylinder of low non-dimensional mass are investigated at low Reynolds numbers. Computations are carried out for 5% blockage. Lock-in is observed for a range of Re and is accompanied with hysteresis at both lower as well as higher Re ends of the synchronisation/lock-in region. It is well known that the lock-in regime for free vibrations depends on the non-dimensional mass of the oscillator. The results from the present computations are compared with the data for forced vibrations from Koopmann (Journal of Fluid Mechanics, 28, 501–512, 1967) on a Y _max/D vs. f* plot, where Y _max is the maximum oscillation amplitude and f* is the ratio of cylinder vibration frequency to the vortex shedding frequency for a stationary cylinder. Good agreement is observed for the critical amplitude needed for onset of synchronisation between the forced and free vibrations. The results from the free vibrations are compared to the predictions from the linear oscillator model by assuming that the forces on the cylinder are unaffected as a result of vibrations. It is found that, for low mass oscillators, the modification of vortex shedding frequency and lift coefficient due to cylinder oscillations leads to the enhancement of the lock-in regime.     

Stability of a fluid-saturated porous medium contained in a vertical cylinder heated from below by forced convection

A linear stability analysis determining the critical Rayleigh number R _c for onset of convection in a bounded vertical cylinder containing a fluid-saturated porous medium is performed for insulated sidewalls, isothermal top surface, and bottom surface heated by forced convection. This Newtonian heating of the bottom surface involves a Biot number Bi that allows consideration of the continuum of boundary conditions ranging from constant heat flux, with global minimum R _min=27.096 found as Bi→0, to isothermal, with global minimum R _min=4π² found as Bi→ ∞. In both cases and for most cylinder aspect ratios, incipient convection sets in as an asymmetric mode, though islands of aspect ratio exist where the onset mode is symmetric. Sample three-dimensional renderings of disturbance temperature distributions showing preferred modes at onset of convection for fixed Bi are provided and an analytical fit to R _min as a function of Bi is given.