Formation of solitary waves on gas-sheared liquid layers

The origin of solitary waves on gas-liquid sheared layers is studied by comparing the behavior of the wave field at sufficiently low liquid Reynolds number, R_L, where solitary waves are observed to form, to measurements at higher R_L where solitary waves do not occur. Observations of the wave field with high-speed video imaging suggest that solitary waves, which appear as a secondary transition of the stratified gas-liquid interface, emanate from existing dominant waves, but that not all dominant waves are transformed. From measurements of interface tracings it is found that for low R_L, waves which have amplitude/substrate depth (a/h) ratios of 0.5–1 occur while for higher R_L, no such waves are observed. A comparison of amplitude/wavelength ratios shows no distinction for different R_L. Consequently, it is conjectured that solitary waves originate from waves with sufficiently large a/h ratios; this change of form being similar to wave breaking. The dimensionless wavenumber is found to be smaller at low R_L, where solitary waves are observed. This suggests that perhaps, larger precursor (to solitary wave) waves are possible because the degree of dispersion, which acts to break waves into separate modes, is lower.     

Study on the head-on collisions of two-dimensional trains of solitons

The head on collisions of trains of solitons induced by a two-dimensional submerged elliptical cylinder at critical speed in shallow water are studied based on velocity potential theory. The boundary value problems are solved through boundary element method (BEM). The nonlinear free surface boundary conditions are satisfied. The mixed Euler–Lagrangian method is adopted to track the free surface through a time stepping scheme. The effects of thickness and velocity of the elliptical cylinder on the evolution of solitary waves have been investigated. Two sets of solitons are truncated from these trains of solitary waves. The head-on collisions of these solitons have been simulated. The wave profiles and velocity fields during collision have been analysed. The propagation of solitary waves is the transmissions of kinetic energy and the collision processes are the results of the dynamic balance of potential energy and kinematic energy.     

基于势流理论的内孤立波与立管作用数值研究

内孤立波是一种发生在水面以下的在世界各个海域广泛存在的大幅波浪, 其剧烈的波面起伏所携带的巨大能量对以海洋立管为代表的海洋结构物产生严重威胁, 分析其传播演化过程的流场特征及立管在内孤立波作用下的动力响应规律对于海洋立管的设计具有重要意义. 本文基于分层流体的非线性势流理论, 采用高效率的多域边界单元法, 建立了内孤立波流场分析计算的数值模型, 可以实时获得内孤立波的流场特征. 根据获得的流场信息, 采用莫里森方程计算内孤立波对海洋立管作用的载荷分布. 将内孤立波流场非线性势流计算模型与动力学有限元模型结合来求解内孤立波作用下海洋立管的动力响应特征, 讨论了内孤立波参数、顶张力大小以及内部流体密度对立管动力响应的影响. 发现随着内孤立波波幅的增大, 海洋立管的流向位移和应力明显增大. 由于上层流体速度明显大于下层, 且在所研究问题中拖曳力远大于惯性力, 因此管道顺流向的最大位移发生在上层区域. 顶张力通过改变几何刚度阵的值进而对立管的响应产生明显影响. 对于弱约束立管, 内部流体的密度对管道的流向位移影响较小.      

Correlation between theoretical and experimental solitary waves

Experimental data on surface solitary waves generated by five methods are given. These data and literature information show that at amplitudes 0.2<a/h<0.6 (h is the initial depth of the liquid), experimental solitary waves are in good agreement with their theoretical analogs obtained using the complete model of liquid potential flow. Some discrepancy is observed in the range of small amplitudes. The reasons why free solitary waves of theoretically limiting amplitude have not been realized in experiments are discussed, and an example of a forced wave of nearly limiting amplitude is given. The previously established fact that during evolution from the state of rest, undular waves break when the propagation speed of their leading front reaches the limiting speed of propagation of a solitary wave is confirmed. Lavrent’ev Institute of Hydrodynamics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 3, pp. 44–52, May–June, 1999.     

Internal wave bore in the shelf zone of the sea

This paper presents the observation results for the internal wave bore in the coastal region of the Sea of Japan with the use of vertical thermistor chains. The data obtained is interpreted by the mathematical models of shallow water in which the effect of nonlinearity and dispersion on the propagation of internal wave trains is taken into account. Within the framework of the theory of multilayer shallow water, the problem of transformation of a solitary wave into an internal bore is solved, and the possibility of recovery of a space-time picture of the flow during the passage of an internal bore in the section between adjacent experimental bottom stations is demonstrated.     

Directional Wave Propagation in a Highly Nonlinear Square Packing of Spheres

We studied the dynamic response of a two-dimensional square packing of uncompressed stainless steel spheres excited by impulsive loadings. We developed a new experimental measurement technique, employing miniature tri-axial accelerometers, to determine the stress wave properties in the array resulting from both an in-plane and out-of-plane impact. Results from our numerical simulations, based on a discrete particle model, were in good agreement with the experimental results. We observed that the impulsive excitations were resolved into solitary waves traveling only through initially excited chains. The observed solitary waves were determined to have similar (Hertzian) properties to the extensively studied solitary waves supported by an uncompressed, uniform, one-dimensional chain of spheres. The highly directional response of this system could be used as a basis to design granular crystals with predetermined wave propagation paths capable of mitigating stress wave energy.     

Generation of acoustic solitary waves in a lattice of Helmholtz resonators

This paper addresses the propagation of high amplitude acoustic pulses through a 1D lattice of Helmholtz resonators connected to a waveguide. Based on the model proposed by Sugimoto (1992), a new numerical method is developed to take into account both the nonlinear wave propagation and the different mechanisms of dissipation: the volume attenuation, the linear viscothermal losses at the walls, and the nonlinear absorption due to the acoustic jet formation in the resonator necks. Good agreement between numerical and experimental results is obtained, highlighting the crucial role of the nonlinear losses. Different kinds of solitary waves are observed experimentally with characteristics depending on the dispersion properties of the lattice.     

Symmetry and Uniqueness of the Solitary-Wave Solution for the Ostrovsky Equation

We consider herein the Ostrovsky equation which arises in modeling the propagation of the surface and internal solitary waves in shallow water, or the capillary waves in a plasma with the effects of rotation. Using the modified sliding method, we prove that the solitary wave moving to the left to the Ostrovsky equation is symmetric about the origin and unique up to translations. We also establish the regularity and decay properties of solitary waves and obtain some results of the nonexistence of solitary wave solutions depending on the wave speed, weak rotation, and dispersive parameter.     

INTERNAL SOLITARY WAVE LOADS EXPERIMENTS AND ITS THEORETICAL MODEL FOR A CYLINDRICAL STRUCTURE

在大型重力式密度分层水槽中, 对内孤立波与圆柱型结构的相互作用特性开展了系列实验. 基于两层流体中 内孤立波的KdV,eKdV和MCC理论, 建立了圆柱型结构内孤立波载荷的理论预报模型, 给出了该载荷理论预报模型中3类内孤立波理论的适用性条件.研究表明, 圆柱型结构内孤立波水平载荷包括水平Froude-Krylov力、附加质量力和拖曳力3个部分, 可以由Morison公式计算, 而内孤立波垂向载荷主要为垂向Froude-Krylov力, 可以由内孤立波诱导动压力计算.系列实验结果表明, 附加质量系数可以取为常数1.0, 拖曳力系数与内孤立波诱导速度场的雷诺数之间为指数函数关系, 而且基于理论预报模型的数值结果与系列实验结果吻合.     

On the influence of material properties on the wave propagation in Mindlin-type microstructured solids

A mathematical model describing 1D wave propagation in Mindlin-type microstructured solids with nonlinearities in the macro- and microscale is used for studying propagation of solitary waves in such media. The results could be used for the stress analysis as well as for the nondestructive testing of material properties. The model equations are solved numerically under the localized initial conditions and periodic boundary conditions by the pseudospectral method. It is demonstrated how the values of the model parameters influence the wave propagation, the evolution and the interaction of waves under the framework of considered models. For this reason the solutions of the model equations are compared under different parameter combinations against one fixed combination of material parameters which is called ‘the reference case’.     

柱形汇聚激波冲击球形重气体界面的实验研究

采用高速纹影法实验研究了柱形汇聚激波与球形重气体界面相互作用的 Richtmyer-Meshkov不稳定性问题. 激波管实验段基于激波动力学理论设计, 将马赫数为1.2 的平面激波转化为柱形汇聚激波, 气体界面由肥皂膜分隔六氟化硫(内)和空气(外)得到. 采用高速摄影机在单次实验中拍摄激波运动的全过程, 对柱形激波的形成进行了实验验证, 并进一步观测了汇聚激波与球形气体界面相互作用过程中的波系发展和气体界面变形以及反射激波同已变形界面二次作用的流场演化. 结果表明: 当柱形汇聚激波穿过气泡界面以后, 气泡左侧界面极点沿激波传播方向保持匀速运动, 气泡右侧界面发展成为射流结构, 气泡主体发展成为涡环结构; 在反射激波的二次作用下, 流场中无序运动显著增强并很快进入湍流混合阶段.     

Combined effects of topography and bottom friction on shoaling internal solitary waves in the South China Sea

A numerical study to a generalized Korteweg-de Vries (KdV) equation is adopted to model the propagation and disintegration of large-amplitude internal solitary waves (ISWs) in the South China Sea (SCS). Based on theoretical analysis and in situ measurements, the drag coefficient of the Chezy friction is regarded as inversely proportional to the initial amplitude of an ISW, rather than a constant as assumed in the previous studies. Numerical simulations of ISWs propagating from a deep basin to a continental shelf are performed with the generalized KdV model. It is found that the depression waves are disintegrated into several solitons on the continental shelf due to the variable topography. It turns out that the amplitude of the leading ISW reaches a maximum at the shelf break, which is consistent with the field observation in the SCS.Moreover, a dimensionless parameter defining the relative importance of the variable topography and friction is presented.     

Frontal collision of internal solitary waves of first mode

The dynamics and energetics of a frontal collision of internal solitary waves (ISW) of first mode in a fluid with two homogeneous layers separated by a thin interfacial layer are studied numerically within the framework of the Navier–Stokes equations for stratified fluid. It was shown that the head-on collision of internal solitary waves of small and moderate amplitude results in a small phase shift and in the generation of dispersive wave train travelling behind the transmitted solitary wave. The phase shift grows as amplitudes of the interacting waves increase. The maximum run-up amplitude during the wave collision reaches a value larger than the sum of the amplitudes of the incident solitary waves. The excess of the maximum run-up amplitude over the sum of the amplitudes of the colliding waves grows with the increasing amplitude of interacting waves of small and moderate amplitudes whereas it decreases for colliding waves of large amplitude. Unlike the waves of small and moderate amplitudes collision of ISWs of large amplitude was accompanied by shear instability and the formation of Kelvin–Helmholtz (KH) vortices in the interface layer, however, subsequently waves again become stable. The loss of energy due to the KH instability does not exceed 5%–6%. An interaction of large amplitude ISW with even small amplitude ISW can trigger instability of larger wave and development of KH billows in larger wave. When smaller wave amplitude increases the wave interaction was accompanied by KH instability of both waves.     

Interaction of rarefaction waves with wet foam

The interaction of rarefaction waves of different shapes with wet water foams is studied experimentally. It is found that the observed values of the pressure are greater, while the surface velocity is lower than the corresponding values predicted by the pseudogas model. The foam breakdown starts as the pressure decreases by 0.3 atm relative to the initial pressure. During downstream propagation of the rarefaction-wave leading edge the propagation velocity decreases.Using of water-based foams as effective screens for damping blast waves in different technological processes has caused considerable interest in studying wave propagation in such systems. The pressure wave dynamics in a foam have been investigated in much detail, both experimentally and theoretically [1–3]. However, the interaction of rarefaction waves with foam has practically never been studied, although it was mentioned in [4] that the unloading phase following the compression wave phase is one of the factors defining the damaging action of blast waves. Besides blast-wave damping, rarefaction wave propagation takes place if such waves are used to breakup foam in oil-producing wells [5].Below, the interaction of rarefaction waves of different shapes with wet water foams is studied. The vertical shock tube described in detail in [3] was used in these experiments.Brest. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 2, pp. 76–82, March–April, 1995.     

Interaction of interfacial and internal waves *

Reflection of internal waves at an interface between the mixed layer and a stratified fluid under the layer is investigated experimentally. Reflection coefficients of the internal waves and transmission coefficients of interfacial waves are obtained for the case in which the interfacial waves propagate along the interface. In order to reveal the energy budget of the two kinds of waves, the above reflection and transmission coefficients are employed in consideration of dissipation of the waves. These coefficients are investigated for different values of a ratio of the frequencies of the interfacial and internal waves. When the ratio is , 1 or 2 at the same horizontal wavenumber, overreflection takes place (i.e., reflection coefficient > 1), while under the same condition the transmission coefficient of the interfacial wave is less than unity. This indicates that the energy source of the overreflection of the internal wave is the interfacial wave. Overreflection also takes place when the two waves have the same frequencies and are in phase. On the other hand, the interfacial wave overtransmits the energy (i.e., transmission coefficient > 1) when the two waves are out of phase.     

Numerical analysis of macrocrack propagation along a bimaterial interface under dynamic loading processes

Advances in computing as well as measurement instrumentation have recently allowed for the investigation of a wider spectrum of physical phenomena in dynamic failure than previously possible. With increasing demand for specialized lightweight, high strength structures, failure of inhomogeneous solids has been receiving increased attention. Such inhomogeneous solids include structural composites such as bonded and sandwich structures, layered and composite materials as well as functionally graded solids. Many of such solids are composed of brittle constituents possessing substantial mismatch in wave speeds, and are bonded together with weak interfaces, which may serve as sites for catastrophic failure (Rosakis and Ravichandran (2000)).In the present study numerical analysis of macrocrack propagation along a bimaterial interface under dynamic loading processes is presented. A general constitutive model of elasto-viscoplastic damaged polycrystalline solids is developed within the thermodynamic framework of the rate type covariance structure with finite set of the internal state variables. A set of the internal state variables is assumed and interpreted such that the theory developed takes account of the effects as follows: (i) plastic non-normality; (ii) softening generated by microdamage mechanisms; (iii) thermomechanical coupling (thermal plastic softening and thermal expansion); (iv) rate sensitivity.To describe suitably the time and temperature dependent effects observed experimentally during dynamic loading processes the kinetics of microdamage has been modified. The relaxation time is used as a regularization parameter. By assuming that the relaxation time tends to zero, the rate independent elastic–plastic response can be obtained. The identification procedure is developed basing on the experimental observations. The finite difference method for regularized elasto-viscoplastic model is used. The edge-cracked bimaterial specimen is considered. In the initial configuration, the height of the specimen is equal to 30 cm, width is 12.5 cm and the length of the initial crack is equal to 2.5 cm. The length of the boundary over which impact is applied is equal to 5 cm, the rise time is fixed at 0.1 μs and the impact velocity is varied. The impact area is localized symmetrically or asymmetrically to the shorter axis of the specimen (symmetry axis of the cohesive band). Basing on the available data of recent experimental observation Rosakis et al. (1999) that have been carried out for relatively thin specimens both the plane stress and plane strain conditions are considered. The material of the specimen is AISI 4340 steel, while PMMA is the cohesive band, both modelled by thermo-elasto-viscoplastic constitutive equations with effects of isotropic hardening and softening generated by microdamage mechanisms and thermomechanical coupling. Fracture criterion based on the evolution of microdamage is assumed. Both, isothermal and adiabatic processes are considered.Particular attention is focused on the investigation of the interactions and reflections of stress waves and the influence of these waves on the propagation of macrocrack within the interface band. The propagation of the macroscopic crack within the material of the interface band for both symmetrical and asymmetrical impact cases has been investigated. It has been found that macrocrack-tip speeds vary from the shear wave speed to the dilatational wave speed of the material and is higher than the Rayleigh surface wave speed. This result is in accord with the experimental observations performed by Rosakis et al. (1999).     

Propagation of a soliton along a fluid-filled pipe

In [1, 2] a mathematical model of the motion of a fluid in a pipe whose axis is a curve in space was discussed and certain simplifications of the problem were studied. The propagation of linear and nonlinear waves in the framework of the model was studied. In the present paper we consider a simple wave flow in a pipe with elastic walls suing one of the models introduced in [1], which, unlike [2], takes into account axial displacements of the pipe. The basic equations describing the propagation of waves in the pipe are obtained.Translated from Zhurnal Prikladnoi Mekhaniki i Technicheskoi Fiziki, No. 3, pp. 58–63, May–June, 1986.     

Comparisons of CIP,compact and CIP‐CSL2 schemes for reproducing internal solitary waves

Six different models were evaluated for reproducing internal solitary waves which occur and propagate in a stratified flow field with a sharp interface. Three stages were used to compute internal solitary waves in a stratified field: (1) first‐phase computation of momentum equations, (2) second‐phase computation of momentum equations, which corresponds to computing the Poisson's equation, and (3) density computation. The six models discussed in this paper consisted of combinations of four different schemes, a three‐point combined compact difference scheme (CCD), a normal central difference scheme (CDS), a cubic‐polynomial interpolation (CIP), and an exactly conservative semi‐Lagrangian scheme (CIP‐CSL2). The residual cutting method was used to solve the Poisson's equation. Three tests were used to confirm the validity of the computations using KdV theory; i.e. the incremental wave speed and amplitude of internal solitary waves, the maximum horizontal velocity and amplitude, and the wave form. In terms of the shape of an internal solitary wave, using CIP for momentum equations was found to provide better performance than CCD. These results suggest one of the most appropriate scheme for reproducing internal solitary waves may be one in which CIP is used for momentum equations and CCD to solve the Poisson's equation. Copyright © 2005 John Wiley & Sons, Ltd.