Re-visiting the head-on collision problem between two solitary waves in shallow water

Upon discovering the wrongness of the statement “although this term does not cause any secularity for this order it will cause secularity at higher order expansion, therefore, that term must vanish” by Su and Mirie [4], in the present work, we studied the head-on collision of two solitary waves propagating in shallow water by introducing a set of stretched coordinates in which the trajectory functions are of order of ϵ², where ϵ is the smallness parameter measuring non-linearity. Expanding the field variables and trajectory functions into power series in ϵ, we obtained a set of differential equations governing various terms in the perturbation expansion. By solving them under non-secularity condition we obtained the evolution equations and also the expressions for phase functions. By seeking a progressive wave solution to these evolution equations we have determined the speed correction terms and the phase shifts. As opposed to the result of Su and Mirie [4] and similar works, our calculations show that the phase shifts depend on both amplitudes of the colliding waves.     

Strongly oblique interactions between internal solitary waves with the same mode

I.IntroductionZabuskyandKruskal(l965)foundthattwoKdVsolitarywavesofthesamemodekeeptheiroriginalshapesandspeedsafterstronginteractions,andcalledthesewavessolitons.However,solitaryx"avessometimestravelintwodimensionalspace,otherthaninonedimensionalspace.Mil…     

General solution for interaction of solitary waves including head-on collisions

A corrected version of the Boussinesq equation for long water waves is derived and its general solution for interaction of any number of solitary waves, including head-on collisions, is given. For two solitary waves in head-on collision (which includes the case of normal reflection) the results agree with the experiments known.     

Analysis of current induced by long internal solitary waves in stratified ocean

An approximate theoretical expression for the current induced by long internal solitary waves is presented when the ocean is continuously or two-layer stratified. Particular attention is paid to characterizing velocity fields in terms of magnitude, flow components, and their temporal evolution/spatial distribution. For the two-layer case, the effects of the upper/lower layer depths and the relative layer density difference upon the induced current are further studied. The results show that the horizontal components are basically uniform in each layer with a shear at the interface. In contrast, the vertical counterparts vary monotonically in the direction of the water depth in each layer while they change sign across the interface or when the wave peak passes through. In addition, though the vertical components are generally one order of magnitude smaller than the horizontal ones, they can never be neglected in predicting the heave response of floating platforms in gravitationally neutral balance. Comparisons are made between the partial theoretical results and the observational field data. Future research directions regarding the internal wave induced flow field are also indicated.     

A numerical study of the load on cylindrical piles exerted by internal solitary waves

A continuously stratified nonlinear model is employed to simulate the generation of internal solitary waves (ISWs) over a sill by tidal flows, and it is shown that the simulated ISW-induced current field basically agrees with that observed. Then the force and torque on a supposed small-diameter vertical cylindrical pile exerted by the simulated ISW packet are calculated. According to the calculation, it is found that, no matter whether the direction of the ISW-induced current is the same as that of the tidal current or not, the force exerted by the ISW would be much larger than that by only the tidal current; if the direction of the ISW-induced current is the same as that of the tidal current, then the torque exerted by the ISW would also be much larger than that by only the tidal current; whilst if the direction of the ISW-induced current is against that of the tidal current, then the torque exerted by the ISW has the same order as that exerted by only the tidal current. It is shown that, under the same conditions, the maximum force on the cylindrical pile is 6.58×10² kN, which is larger than that by the modal separation method of Cai et al., whilst the maximum torque is 2.46×10⁵ kN m, which is less than that given by Cai et al. During the passage of the ISW, the time series of the force and torque on the cylindrical pile can also be shown. Finally, the effect of the characteristics of the Gaussian sill on the force is studied, and the resulted empirical formulas on the force with the wave amplitude and the non-dimensional variable of the sill parameters are put forward.     

Head-on collision between two mKdV solitary waves in a two-layer fluid system

In this paper, based on the equations presented in [2], the head-on collision between two solitary waves described by the modified KdV equation (the mKdV equation, for short) is investigated by using the reductive perturbation method combined with the PLK method. These waves propagate at the interface of a two-fluid system, in which the density ratio of the two fluids equals the square of the depth ratio of the fluids. The second order perturbation solution is obtained. It is found that in the case of disregarding the nonuniform phase shift, the solitary waves preserve their original profiles after collision, which agrees with Fornberg and Whitham's numerical result of overtaking collision161 whereas after considering the nonuniform phase shift, the wave profiles may deform after collision.     

On the Navier–Stokes equations simulation of the head-on collision between two surface solitary waves

The scope of this Note is to show the results obtained for simulating the two-dimensional head-on collision of two solitary waves by solving the Navier–Stokes equations in air and water. The work is dedicated to the numerical investigation of the hydrodynamics associated to this highly nonlinear flow configuration, the first numerical results being analyzed. The original numerical model is proved to be efficient and accurate in predicting the main features described in experiments found in the literature. This Note also outlines the interest of this configuration to be considered as a test-case for numerical models dedicated to computational fluid mechanics. To cite this article: P. Lubin et al., C. R. Mecanique 333 (2005).     

Interactions of atmospheric solitary waves of different modes

The interactions of atmospheric solitary waves with different modes are investigated by a perturbation method. The model considered in this paper consists of a lower layer with exponential density profile and an infinitely deep upper layer with constant density. The analysis show that the waves obey the Benjamin-Ono equation before and after interaction, and the main effect of the interaction is the phase shifts for each wave. The project supported by the National and Shanghai Education Commission of Science Foundation     

Diffraction of solitary water waves around three-dimensional floating bodies

By using the matched asymptotic expansion method and the idea of edge layer, a mathematic model for describing the interaction between weakly nonlinear shallow-water waves and three-dimensional floating bodies is formed in the paper. As a numerical example, the diffraction of a solitary wave around a vertically floating circular cylinder has been investigated and the results are presented. The present method can further be extended to the study of wave diffraction around floating bodies of general shape. The project is supported by the National Natural Science Foundation of China.     

The stability of finite-amplitude interfacial solitary waves

The linear stability of finite-amplitude interfacial gravity solitary waves propagating in a two-layer fluid is investigated analytically focusing on the occurrence of an exchange of stability. We make an asymptotic analysis for small growth rates of infinitesimal disturbances, and explicitly obtain their growth rates near an exchange of stability. The result indicates that an exchange of stability occurs at every stationary value of the total energy of the solitary waves. It also gives us information whether the number of growing modes increases or decreases after experiencing the exchange of stability. We apply these analytical results to specific interfacial solitary waves, and find various features on their stability that are not seen in the case of surface solitary waves.     

三维海洋内孤立波数值水槽造波研究

海洋内孤立波因其分布广泛和携带巨大能量,对于潜艇安全航行影响很大.本文采用有限体积自适应半结构多重网格法求解Navier-Stokes方程,并用VOF方法追踪两层流体界面,应用双推板造波法进行内孤立波数值造波,建立了两层流体中的内孤立波数值水槽.数值模拟结果证实了该数值水槽数值造波的有效性和可靠性,为后续研究打下了基础...     

Asymptotic models of internal stationary waves

Equations of stationary long waves on the interface between a homogeneous fluid and an exponentially stratified fluid are considered. An equation of the second-order approximation of the shallow water theory inheriting the dispersion properties of the full Euler equations is used as the basic model. A family of asymptotic submodels is constructed, which describe three different types of bifurcation of solitary waves at the boundary points of the continuous spectrum of the linearized problem. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 49, No. 4, pp. 151–161, July–August, 2008.     

Interpretation of wake instability at slip line in rotating detonation

ABSTRACT

In studies on instabilities of flowfield in rotating detonation, one of the most common concerns is the instability at the slip line originating from the conjunction of the detonation wave and oblique shock. Using Euler equations associated with the 7-species-and-8-reaction finite-rate chemical reaction model of hydrogen/air mixtures, further studies are performed to simulate the 2-D rotating detonation, and the flow mechanism of instability at the slip line is investigated in depth. The results show that the distinct wake profile exists at the slip line, which is different from the typical mixing layer. Analysis indicates that the generation of wake is caused by the transition shock between the detonation wave and oblique shock. Because of the wake profile, the vorticity distribution therein appears in a double-layer layout, and different evolutions exist in different vorticity layers. Based on the velocity profile across the slip line, the analysis by the linear stability theory is made, and two main unstable modes which have different shape profiles and phase velocities are found. Discrete Fourier transformation is utilised to analyse the numerical results, and similar shape profiles are obtained. A general coincidence in velocity of vortex movement is also attained between the theoretical predictions and simulations. Investigations show that the wake instability is responsible for the unstable mechanism, and corresponding unstable structures differ from the canonical ones in typical mixing layers.     

Numerical detection and generation of solitary waves for a nonlinear wave equation

This paper presents an automatic algorithm for detecting and generating solitary waves of nonlinear wave equations. With this purpose, dynamic simulations are carried out, the solution of which evolves into a main pulse along with smaller dispersive tails. The solitary waves are detected automatically by the algorithm by checking that they have constant amplitude and are symmetric respect to its maximum value. Once the main wave has been detected, the algorithm cleans the dispersive tails for time enough so that the solitary wave is obtained with the required precision.In order to use our algorithm, we need a spatial discretization with local basis. The numerical experiments are carried out for the BBM equation discretized in space with cubic finite elements along with periodic boundary conditions. Moreover, a geometric integrator in time is used in order to obtain good approximations of the solitary waves.     

A design of a two-dimensional,supersonic KH experiment on OMEGA-EP

An experiment meant to investigate the evolution of single mode Kelvin–Helmholtz (KH) instability in the supersonic regime is presented and theoretically analyzed. This experiment is intended to provide a direct measurement of the two-dimensional vortex evolution so that the high-Mach-number effects can be measured. The proposed design takes advantage of the ability of OMEGA-EP to drive experiments for up to 30 ns to produce steady conditions for KH that endure long enough to observe substantial growth. KH growth for the proposed design has been analyzed using two-dimensional numerical simulations. The results were compared to synthetic temporal KH numerical simulations using non-dimensional scaling in the low and high Mach number regime. The comparisons show that the growth in the high Mach number regime is expected to be suppressed by up to a factor of two. The effects of two-dimensional rarefactions from the lateral boundaries of the experimental system were also investigated. It was found that they introduce no major uncertainties or hazards to the experiment. We produced simulated radiographs, which show that the proposed experimental system will enable observation of the KH structures. An experiment of this kind has not yet been performed, and therefore would serve to validate numerical results and analytical models presented here and in the literature.     

General solution for interaction of solitary waves including head-on collisions

A general solution of the Boussinesq equation is presented which solves the problem of interaction of any number of right-going and left-going solitary waves. The solution relies on the exact solution of Gardner, Greene, Kruskal, and Miura (1967), and has the same degree of accuracy as that solution, but has a wider scope of application. It is much simpler than, but as accurate as, Hirota's exact solution (1973) of the Boussinesq equation, to which the present solution is compared for the simplest case of two solitary waves in head-on collision.     

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.     

Recognition of coherent structures in the boundary layer of a low-pressure-turbine blade for different free-stream turbulence intensity levels

Particle Image Velocimetry (PIV) has been adopted to analyze the instantaneous flow field developing on a high-lift turbine blade profile operating under low and elevated free-stream turbulence conditions (FSTI). Results reported in the paper allow us to analyze the dynamics leading to transition and separation of the suction side boundary layer, looking to generation, propagation and breakdown of coherent structures observed in the two different FSTI cases. To this end, measurements have been performed in two orthogonal planes. Results obtained in the blade-to-blade plane allow the detailed characterization of the propagation of Kelvin–Helmholtz (KH) rolls generating, at low FSTI condition, as a consequence of a non-reattaching separation. Otherwise, data in the wall-parallel plane allow recognizing the presence of three-dimensional disuniformities induced at high FSTI by low and high speed streaks (Klebanoff mode). The sinuous breakdown of boundary layer streaks generates other complex three-dimensional coherent structures such as hairpin or cane-like vortices that induce transition. Proper Orthogonal Decomposition (POD) has been adopted to in depth characterize these structures, thus further explaining the mechanisms through which the free-stream turbulence intensity modify the transition/separation processes of the suction side boundary layer of an highly loaded low pressure turbine blade.     

A solitary group of two-dimensional deep-water waves due to a concentrated force-couple

In 1986, Yih found a solitary group of deep-water waves caused by a travelling and oscillatory concentrated force. In this paper, such a single group caused by a travelling and oscillatory force-couple is presented. The expression of wave elevation is also derived and an intuitive explanation of the results is provided.     

Nonlinear dynamics of hydrostatic internal gravity waves

Stratified hydrostatic fluids have linear internal gravity waves with different phase speeds and vertical profiles. Here a simplified set of partial differential equations (PDE) is derived to represent the nonlinear dynamics of waves with different vertical profiles. The equations are derived by projecting the full nonlinear equations onto the vertical modes of two gravity waves, and the resulting equations are thus referred to here as the two-mode shallow water equations (2MSWE). A key aspect of the nonlinearities of the 2MSWE is that they allow for interactions between a background wind shear and propagating waves. This is important in the tropical atmosphere where horizontally propagating gravity waves interact together with wind shear and have source terms due to convection. It is shown here that the 2MSWE have nonlinear internal bore solutions, and the behavior of the nonlinear waves is investigated for different background wind shears. When a background shear is included, there is an asymmetry between the east- and westward propagating waves. This could be an important effect for the large-scale organization of tropical convection, since the convection is often not isotropic but organized on large scales by waves. An idealized illustration of this asymmetry is given for a background shear from the westerly wind burst phase of the Madden–Julian oscillation; the potential for organized convection is increased to the west of the existing convection by the propagating nonlinear gravity waves, which agrees qualitatively with actual observations. The ideas here should be useful for other physical applications as well. Moreover, the 2MSWE have several interesting mathematical properties: they are a system of nonconservative PDE with a conserved energy, they are conditionally hyperbolic, and they are neither genuinely nonlinear nor linearly degenerate over all of state space. Theory and numerics are developed to illustrate these features, and these features are important in designing the numerical scheme. A numerical method is designed with simplicity and minimal computational cost as the main design principles. Numerical tests demonstrate that no catastrophic effects are introduced when hyperbolicity is lost, and the scheme can represent propagating discontinuities without introducing spurious oscillations.