A kind of extended Korteweg-de Vries equation and solitary wave solutions for interfacial waves in a two-fluid system

This paper considers interfacial waves propagating along the interface between a two-dimensional two-fluid with a flat bottom and a rigid upper boundary. There is a light fluid layer overlying a heavier one in the system, and a small density difference exists between the two layers. It just focuses on the weakly non-linear small amplitude waves by introducing two small independent parameters: the nonlinearity ratio epsilon, represented by the ratio of amplitude to depth, and the dispersion ratio mu, represented by the square of the ratio of depth to wave length, which quantify the relative importance of nonlinearity and dispersion. It derives an extended KdV equation of the interfacial waves using the method adopted by Dullin et al in the study of the surface waves when considering the order up to O(mu(2)). As expected, the equation derived from the present work includes, as special cases, those obtained by Dullin et al for surface waves when the surface tension is neglected. The equation derived using an alternative method here is the same as the equation presented by Choi and Camassa. Also it solves the equation by borrowing the method presented by Marchant used for surface waves, and obtains its asymptotic solitary wave solutions when the weakly nonlinear and weakly dispersive terms are balanced in the extended KdV equation.

Validity ranges of interfacial wave theories in a two-layer fluid system

In the present paper, we endeavor to accomplish a diagram, which demarcates the validity ranges for interfacial wave theories in a two-layer system, to meet the needs of design in ocean engineering. On the basis of the available solutions of periodic and solitary waves, we propose a guideline as principle to identify the validity regions of the interfacial wave theories in terms of wave period T, wave height H, upper layer thickness d(1), and lower layer thickness d(2), instead of only one parameter-water depth d as in the water surface wave circumstance. The diagram proposed here happens to be Le Mehautes plot for free surface waves if water depth ratio r = d(1)/d(2) approaches to infinity and the upper layer water density rho(1) to zero. On the contrary, the diagram for water surface waves can be used for two-layer interfacial waves if gravity acceleration g in it is replaced by the reduced gravity defined in this study under the condition of sigma = (rho(2) - rho(1))/rho(2) -> 1.0 and r > 1.0. In the end, several figures of the validity ranges for various interfacial wave theories in the two-layer fluid are given and compared with the results for surface waves.

Wave generation by the falling rock in the two-dimensional wave tank

Wave generation by the falling rock in the two-dimensional wave tank is experimentally and numerically studied, where the numerical model utilizes the boundary element method to solve the fully nonlinear potential flow theory. The wave profiles at different times are measured in the laboratory, which are also used to test the numerical model. Comparisons show that the experimental and numerical results are in good agreement, and the numerical model can be used to simulate the wave generation due to the submarine rock falling. Further numerical tests on the influences of the rock size, density, initial position and the falling angle on the wave elevation of the generated waves are performed, respectively. The results show that the size and density of the rock have strong effects on the maximum elevation of the generated wave, while the effects of the initial position and the falling angle of the rock are also significant. When the size or the density of the rock increases, the maximum elevation of the generated wave increases. The same effect on the generated wave would be produced if the initial position of the rock becomes closer to the surface, or the falling angle between the falling route and the vertical direction turns larger. In addition, the present numerical tests reveal that the submarine rock falling provides a new generation method for the breaking wave in the wave tank.

Influence of the earth rotation on the interfacial wave solutions and wave-induced tangential stress in a two-layer fluid system

Interfacial waves and wave-induced tangential stress are studied for geostrophic small amplitude waves of two-layer fluid with a top free surface and a flat bottom. The solutions were deduced from the general form of linear fluid dynamic equations of two-layer fluid under the f-plane approximation, and wave-induced tangential stress were estimated based on the solutions obtained. As expected; the solutions derived from the present work include as special cases those obtained by Sun et al. (2004. Science in China, Set. D, 47(12): 1147-1154) for geostrophic small amplitude surface wave solutions and wave-induced tangential stress if tire density of the upper layer is much smaller than that of the lower layer. The results show that the interface and the surface will oscillate synchronously, and the influence of the earth's rotation both on the surface wave solutions and the interfacial wave solutions should be considered.

A new method for the estimation of oceanic mixed-layer depth using shipboard X-band radar images

Shipboard X-band radar images acquired on 24 June 2009 are used to study nonlinear internal wave characteristics in the northeastern South China Sea. The studied images show three nonlinear internal waves in a packet. A method based on the Radon Transform technique is introduced to calculate internal wave parameters such as the direction of propagation and internal wave velocity from backscatter images. Assuming that the ocean is a two-layer finite depth system, we can derive the mixed-layer depth by applying the internal wave velocity to the mixed-layer depth formula. Results show reasonably good agreement with in-situ thermistor chain and conductivity-temperature-depth data sets.

南海文昌海域内波特征研究

2005年中科院海洋所对琼东文昌海域进行了精细内波观测。本文对此次观测所获得的资料进行分析，研究了观测海域内孤立波、内潮、内波波致流和内波传播方向的统计特征，得到观测海域内波的基本特征。首先，观测海域内孤立波频繁发生。在149天的观测当中，一共出现140个孤立波列，出现频率达到了0.94次/天。出现频率最高的振幅是32米，最常见的内孤立波持续时间时20分钟。80%的内孤立波都在以下范围之内：持续时间从13分钟到32分钟，振幅从16米到46米。在这一区间范围之类，内孤立波的持续时间和振幅并没有明显的函数关系，表现出随机分布的特征。同时内孤立波的振幅有着明显的季节变化特征。从4月中旬到6月中旬内孤立波频繁出现且振幅较大，振幅超过35米的内孤立波相当常见；而6月中旬到9月中旬，内孤立波的发生频率明显下降，振幅也偏低，振幅超过35米的内孤立波属于少数情况。其次，观测地的内潮以日频率的振动为主，间隔有半日频率的内潮出现，与当地天文潮情况相似。其振幅小于内孤立波的振幅，在每个观测日都能观测到。内潮的振幅有着明显的季节变化，7月是内潮振幅最小的一个月。再次，观测地内波波致流的大小和剪切层位置有着明显的季节变化。在4月和5月内波波致流流速较大而且靠近海底。而在6月以后，内波波致流较弱，并且剪切层也更为靠近海面。最后，观测地内波可能在局地生成，由于主要呈现全日频率可以排除在吕宋海峡生成的可能。同时传播方向同样随季节变化。4月和5月观测到的内波传播方向主要与等深线线平行，极少被卫星遥感图片观测到。6月、7月、8月和9月传播方向与等深线垂直。

The effect of variable currents on internal solitary waves

The effect of variable currents on internal solitary waves is described within the context of a variable coefficient Korteweg-de Vries (KdV) equation, and the approximate slowly varying, solitary-wave solution of this equation. The general theory which leads to the variable coefficient KdV equation is described; a derivation for the special case when the solitary wave and the current are aligned in the same direction is given in the Appendix. Using further simplifications and approximations, a number of analytical expressions are obtained for the variation in the solitary wave amplitude resulting from variable shear in the basic current or from when the basic current is a depth-independent flow which is a simple representation of a geostrophic current, tidal flow or inertial wave.

Wave loading on floating platforms by internal solitary waves

Morison's equation is used for estimating internal solitary wave-induced forces exerted on SPAR and semi-submersible platforms. And the results we got have also been compared to ocean surface wave loading. It is shown that Morison's equation is an appropriate approach to estimate internal wave loading even for SPAR and semi-submersible platforms, and the internal solitary wave load on floating platforms is comparable to surface wave counterpart. Moreover, the effects of the layers with different thickness on internal solitary wave force are investigated.

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.

Spatial structure of internal tidal waves generated on weak topographies

The generation of internal gravity waves by barotropic tidal flow passing over a two-dimensional topography is investigated. Rather than calculating the conversion of tidal energy, this study focuses on delineating the geometric characteristics of the spatial structure of the resulting internal wave fields (i.e., the configurations of the internal beams and their horizontal projections) which have usually been ignored. it is found that the various possible wave types can be demarcated by three characteristic frequencies: the tidal frequency, wo; the buoyancy frequency, N; and the vertical component of the Coriolis vector or earth's rotation.f. When different possibilities arising from the sequence of these frequencies are considered, there occur 12 kinds of wave structures in the full 3D space in contrast to the 5 kinds identified by the 2D theory. The constant wave phase lines may form as ellipses or hyperbolic lines on the horizontal plane, provided the buoyancy frequency is greater or less than the tidal frequency. The effect that stems from the consideration of the basic flow is also found, which not only serves as the reason for the occurrence of higtter harmonics but also increases the wave strength in the direction of basic flow. (C) 2009 Elsevier B.V. All rights reserved.

Generation of internal and surface waves by seafloor movement in a two layer fluid system

Internal and surface waves generated by the deformations of the solid bed in a two layer fluid system of infinite lateral extent and uniform depth are investigated. An integral solution is developed for an arbitrary bed displacement on the basis of a linear approximation of the complete description of wave motion using a transform method (Laplace in time and Fourier in space) analogous to that used to study the generation of tsunamis by many researchers. The theoretical solutions are presented for three interesting specific deformations of the seafloor; the spatial variation of each seafloor displacement consists of a block section of the seafloor moving vertically either up or down while the time-displacement history of the block section is varied. The generation process and the profiles of the internal and surface waves for the case of the exponential bed movement are numerically illustrated, and the effects of the deformation parameters, densities and depths of the two layers on the solutions are discussed. As expected, the solutions derived from the present work include as special cases that obtained by Kervella et al. [Theor Comput Fluid Dyn 21:245-269, 2007] for tsunamis cased by an instantaneous seabed deformation and those presented by Hammack [J Fluid Mech 60:769-799, 1973] for the exponential and the half-sine bed displacements when the density of the upper fluid is taken as zero.

Micro-shock waves generated inside a fluid jet impinging on plane wall

An Nd:glass laser pulse (18 ns, 1.38 J) is focused in a tiny area of about 100-mum diam under ambient conditions to produce micro-shock waves. The laser is focused above a planar surface with a typical standoff distance of about 4 mm, The laser energy is focused inside a supersonic circular jet of carbon dioxide gas produced by a nozzle with internal diameter of 2.9 mm and external diameter of 8 mm, Nominal value of the Mach number of the jet is around 2 with the corresponding pressure ratio of 7.5 (stagnation pressure/static pressure at the exit of the nozzle), The interaction process of the micro-shock wave generated inside the supersonic jet with the plane wall is investigated using double-pulse holographic interferometry. A strong surface vortex field with subsequent generation of a side jet propagating outward along the plane wail is observed. The interaction of the micro-shock wave with the cellular structure of the supersonic jet does not seem to influence the near surface features of the flowfield. The development of the coherent structures near the nozzle exit due to the upstream propagation of pressure waves seems to be affected by the outward propagating micro-shock wave. Mach reflection is observed when the micro-shock wave interacts with the plane wall at a standoff distance of 4 mm, The Mach stem is slightly deflected, indicating strong boundary-layer and viscous effects near the wall. The interaction process is also simulated numerically using an axisymmetric transient laminar Navier-Stokes solver. Qualitative agreement between experimental and numerical results is good.

The induced flow field by internal solitary wave and its action on cylindrical piles in the stratified ocean

The induced flow fields by internal solitary waves and its actions on cylindrical piles in density stratified ocean with a basic density profile and a basic velocity profile are investigated. Some results, such as the time evolution of flow fields and hydrodynamic forces on the piles are yielded both by theoretical analysis and numerical calculation for general and specific cases. Several kinds of ambient sea conditions of the South China Sea are specified for numerical simulation. Moreover, the effects of relative density difference, depth ratio and wave steepness on maximal total force and total torque are analyzed.

Interaction of linear waves with an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth

The linear water wave scattering and radiation by an array of infinitely long horizontal circular cylinders in a two-layer fluid of infinite depth is investigated by use of the multipole expansion method. The diffracted and radiated potentials are expressed as a linear combination of infinite multipoles placed at the centre of each cylinder with unknown coefficients to be determined by the cylinder boundary conditions. Analytical expressions for wave forces, hydrodynamic coefficients, reflection and transmission coefficients and energies are derived. Comparisons are made between the present analytical results and those obtained by the boundary element method, and some examples are presented to illustrate the hydrodynamic behavior of multiple horizontal circular cylinders in a two-layer fluid. It is found that for two submerged circular cylinders the influence of the fluid density ratio on internal-mode wave forces is more appreciable than surface-mode wave forces, and the periodic oscillations of hydrodynamic results occur with the increase of the distance between two cylinders; for four submerged circular cylinders the influence of adding two cylinders on the wave forces of the former cylinders is small in low and high wave frequencies, but the influence is appreciable in intermediate wave frequencies.

Dynamic Characteristics of Spherically Converging Detonation Waves

The spherically converging detonation wave was numerically investigated by solving the one-dimensional multi-component Euler equations in spherical coordinates with a dispersion-controlled dissipative scheme. Finite rate and detailed chemical reaction models were used and numerical solutions were obtained for both a spherical by converging detonation in a stoichiometric hydrogen-oxygen mixture and a spherically focusing shock in air. The results showed that the post-shock pressure approximately arises to the same amplitude in vicinity of the focal point for the two cases, but the post-shock temperature level mainly depends on chemical reactions and molecular dissociations of a gas mixture. While the chemical reaction heat plays an important role in the early stage of detonation wave propagation, gas dissociations dramatically affect the post-shock flow states near the focal point. The maximum pressure and temperature, non-dimensionalized by their initial value, are approximately scaled to the propagation radius over the initial detonation diameter. The post-shock pressure is proportional to the initial pressure of the detonable mixture, and the post-shock temperature is also increased with the initial pressure, but in a much lower rate than that of the post-shock pressure.