Localization of wave motions in a fluid-filled cylinder

The properties of harmonic surface waves in a fluid-filled cylinder made of a compliant material are studied. The wave motions are described by a complete system of dynamic equations of elasticity and the equation of motion of a perfect compressible fluid. An asymptotic analysis of the dispersion equation for large wave numbers and a qualitative analysis of the dispersion spectrum show that there are two surface waves in this waveguide system. The first normal wave forms a Stoneley wave on the inside surface with increase in the wave number. The second normal wave forms a Rayleigh wave on the outside surface. The phase velocities of all the other waves tend to the velocity of the shear wave in the cylinder material. The dispersion, kinematic, and energy characteristics of surface waves are analyzed. It is established how the wave localization processes differ in hard and compliant materials of the cylinder __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 4, pp. 72–86, April 2008.     

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

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

被均匀流缓慢调制的有限水深毛细重力波

非线性的存在会产生高次谐波,这些谐波又反作用于原来的低次谐波,使波幅发生缓慢变化,从而产生缓慢调制现象．这里从考虑均匀流作用下的毛细重力水波基本方程出发,在不可压缩、无旋、无黏条件假设下,使用多重尺度分析方法推导出了在均匀流影响下有限深水毛细重力波振幅所满足的非线性Schr?dinger方程(NLSE)．分析了NLSE解的调制不稳定性．给出了毛细重力波调制不稳定的条件和钟型孤立波产生的条件．分析了无量纲最大不稳定增长率随无量纲水深和表面张力的变化趋势．同时给出了无量纲不稳定增长率随无量纲微扰动波数变化的曲线,呈现出了先增大后减小的趋势．最后指出均匀顺流减小了无量纲不稳定增长率及最大增长率,逆流增大了它们．由表面张力作用产生的毛细波及重力与表面张力共同作用产生的毛细重力波,与流的相互作用可以改变海表粗糙度和海洋上层流场结构,进而影响海气界面动量、热量及水汽的交换．了解海表这些短波动力机制,对卫星遥感的精确测量、海气相互作用的研究及海气耦合模式的改进等有重要意义．      

Problem of Low-frequency Localized Oscillations in a Thin Film with Growing Islands

Influence of localization waves to islands growing on a thin film is investigated. The film is modelled as a fluid layer covered by an inertial surface with the variable density of mass and surface tension. Mathematically, the problem is reduces to analysis of a system of non-linear equations describing the growth of island nuclei and wave propagation in the films. The existence of trapped modes for the corresponding frequency-domain problem is established. We show that for large time wave localization near islands gives some contribution in the increase of the velocity of island growth.     

The Regularity and Local Bifurcation of¶Steady Periodic Water Waves

Steady periodic water waves on infinite depth, satisfying exactly the kinematic and dynamic boundary conditions on the free surface, with or without surface tension, are given by solutions of a rather tidy nonlinear pseudo-differential operator equation for a 2π-periodic function of a real variable. Being an Euler-Lagrange equation, this formulation has the advantage of gradient structure, but is complicated by the fact that it involves a non-local operator, namely the Hilbert transform, and is quasi-linear. This paper is a mathematical study of the equation in question. First it is shown that its W ^1,2 solutions are real analytic. Then bifurcation theory for gradient operators is used to prove the existence of (non-zero) small amplitude waves near every eigenvalue (irrespective of multiplicity) of the linearised problem. Finally it is shown that when surface tension is absent there are no sub-harmonic bifurcations or turning points at the outset of the branches of Stokes waves which bifurcate from the trivial solution. Accepted: December 6, 1999     

Numerical simulation and experimental verification of silicone oil flow over magnetic fluid under applied magnetic field

Two-layer flow of magnetic fluid and non-magnetic silicone oil was simulated numerically. The continuity equation, momentum equations, kinematic equation, and magnetic potential equation were solved in two-dimensional Cartesian coordinate. PLIC （piecewise linear integration calculation） VOF （volume of fluid） scheme was employed to track the free interface. Surface tension was treated via a continuous surface force （CSF） model that ensures robustness and accuracy. The influences of applied magnetic field, inlet velocity profile, initial surface disturbance of interface and surface tension were analyzed. The computed interface shapes at different conditions were compared with experimental observation.     

有限深水中骑行波的显式Hamilton描述

小尺度波（扰动波）迭加在大尺度波（未受扰动波）上形成的波动一般之为“骑行波”。研究了有限可变深度的理想不可压缩流体中的骑行波的显式Hamliltn表示，考虑了自由面上流体与空气之间的表面张力。采用自由面高度和自由面上速度势构成的Hamilton正则变量表示骑行波的动能密度，并在未受扰动波的自由面上作一阶展开。运用复变函数论方法处理了二维流动。先用保角变换将物理平面上的流动区域变换到复势平面上的无限长带形区域，然后在复势平面上用Fourier变换解出Laplace方程，最后经Fourier逆变换求出了扰动波速度热所满足的积分方程。作为特例考虑了平坦底部的流动，导出了动能密度的显式表达式。这里给出的积分方程可以替代相当难解的Hamilton正则方程。通过求解积分方程可得出agrange密度的显式表达式。本文提出的方法约研究骑行波的Hamilton描述以及波的相互作用问题提供了新的途径，有助于了解海面的小尺度波的精细结构。     

Wave propagation through mangrove forests in the presence of a viscoelastic bed

The influence of viscoelastic ocean beds on the characteristics of surface waves passing through mangrove forests is analyzed under the assumption of linearized water wave theory in two dimensions. The trunks of the mangroves are assumed to be in the upper-layer inviscid fluid domain, whilst the roots are inside the viscoelastic bed. The associated equation of motion is obtained by coupling the Voigt’s model for flow within the viscoelastic medium with the equation of motion in the presence of mangroves. The modified dynamic conditions are coupled with the kinematic conditions to obtain the boundary condition at the free surface and the interface of the two fluids consisting of the upper layer inviscid fluid and the viscoelastic fluid bed. To understand the effects of bed viscosity as well as elasticity on energy dissipation, the complex dispersion relation associated with the plane progressive wave is derived and analyzed. Effect of physical parameters associated with mangroves and viscoelastic bed on wave motion in surface and internal modes are computed and analyzed to understand their roles in attenuating wave effects. The present model will be useful in the better understanding of wave propagation through mangroves in the coastal zone having muddy seabed.     

Wave regimes on power-law fluid film flowing down a porous plane

Non-linear waves on the surface of a falling film of power-law fluid on a vertical porous plane are investigated. The waves are described by evolution equations generalising equations previously derived in the case of solid plane. It is shown that the slip condition on the interface between pure liquid and the porous substrate drastically changes structure of the steady waves travelling in the film.     

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.     

滚波演化中聚合过程的数值模拟研究

滚波是一种重力作用下自由液面失稳诱发的水面波动现象, 通常可分为具有相对稳定波形和波速的周期性滚波与波形和波速不断变化的不规则滚波(自然滚波). 不规则滚波的相互作用和发展演化过程十分复杂, 至今对其认识尚不成熟. 本文采用基于雷诺平均Navier-Stokes方程的立面二维数值模型, 对不规则滚波发展过程中的吸收聚合和追赶聚合现象进行了数值模拟研究. 分析了两种聚合模式的演化过程, 给出了滚波聚合过程中完整的波形、波速、速度剖面以及湍流黏性等重要信息. 结果表明滚波的聚合过程是不规则滚波演化和增长的重要机制, 在特定条件下滚波增长由自然增长模式转变为以吸收聚合和追赶聚合为主的增长模式. 滚波聚合过程中, 依次经历后波追赶、爬升、与前波合并、内部流场调制等多个步骤, 最终形成一个具有更大波长和波高的滚波. 本文发现了在3个滚波间距较近的情况下, 会发生二重聚合现象, 即后两个滚波首先聚合, 然后与前波进一步聚合形成一个新的滚波.      

Dynamics of a Fluid in a Rotating Horizontal Cylinder

The behavior of a low-viscosity fluid in a rotating horizontal circular cylinder is investigated experimentally. The stability of the centrifuged layer, the motion of the fluid with respect to the cavity, the excitation of inertial waves on the fluid surface, and the effect of the waves on the stability and flow structure are studied over a wide region of relative occupancy of the cavity. The results are analyzed from the viewpoint of vibrational mechanics in which the motion is generated by the oscillations of the fluid with respect to the cavity and the gravity force plays the role of the force oscillating in the cavity reference system.     

A Dimension-Breaking Phenomenon for Water Waves with Weak Surface Tension

It is well known that the water-wave problem with weak surface tension has small-amplitude line solitary-wave solutions which to leading order are described by the nonlinear Schrödinger equation. The present paper contains an existence theory for three-dimensional periodically modulated solitary-wave solutions which have a solitary-wave profile in the direction of propagation and are periodic in the transverse direction; they emanate from the line solitary waves in a dimension-breaking bifurcation. In addition, it is shown that the line solitary waves are linearly unstable to long-wavelength transverse perturbations. The key to these results is a formulation of the water wave problem as an evolutionary system in which the transverse horizontal variable plays the role of time, a careful study of the purely imaginary spectrum of the operator obtained by linearising the evolutionary system at a line solitary wave, and an application of an infinite-dimensional version of the classical Lyapunov centre theorem.     

A film flow model for analysing gravity-driven,thin wavy fluid films

To analyse the physics underlying gravity-driven runoff of thin wavy films, a film flow model is developed, and is solved with computational fluid dynamics. This model is based on the lubrication theory, and takes into account the gravitational, wall shear and surface tension forces. A key characteristic of the model is that it assumes only one computational cell over the film height, which enables studying film flow on larger computational domains. A main aim of this study is to perform a detailed validation of the numerical model. The film flow model is validated against several experiments of gravity-driven, thin fluid films on smooth surfaces. The time-averaged film thickness and the fluid speed profiles predicted by the model show very good agreement with experimental results. Similarly, the film flow model is able to predict the wave speeds with sufficient accuracy. The energy spectra of the waves, where higher frequency waves are present in film flows at higher Reynolds numbers, show an exponentially decaying trend at these high frequencies. The model performs better than the Nusselt equation for film flows, which under-predicts the time-averaged film thickness and over-predicts the time-averaged fluid speeds, even for flows at low Reynolds numbers. The film flow model is compared qualitatively for fingering behaviour. This model also allows to investigate film flows on large surfaces, which can be rough, curved and of complex geometrical shape.     

Features of how surface waves form in a cylinder made of a hard material and filled with a fluid

The properties of harmonic surface waves in an elastic cylinder made of a rigid material and filled with a fluid are studied. The problem is solved using the dynamic equations of elasticity and the equations of motion of a perfect compressible fluid. It is shown that two surface (Stoneley and Rayleigh) waves exist in this waveguide system. The first normal wave generates a Stoneley wave on the inner surface of the cylinder. If the material is rigid, no normal wave exists to transform into a Rayleigh wave. The Rayleigh wave on the outer surface forms on certain sections of different dispersion curves. The kinematic and energy characteristics of surface waves are analyzed. As the wave number increases, the phase velocities of all normal waves, except the first one, tend to the sonic velocity in the fluid from above __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 9, pp. 48–62, September 2007.     

Small Reynolds number instabilities in two-layer Couette flow

Instabilities in two-layer Couette flow are investigated from a small Reynolds number expansion of the eigenvalue problem governing linear stability. The wave velocity and growth rate are given explicitly, and previous results for long waves and short waves are retrieved as special cases. In addition to the inertial instability due to viscous stratification, the flow may be subject to the Rayleigh–Taylor instability. As a result of the competition of these two instabilities, inertia may completely stabilise a gravity-unstable flow above some finite critical Froude number, or conversely, for a gravity-stable flow, inertia may give rise to finite wavenumber instability above some finite critical Weber number. General conditions for these phenomena are given, as well as exact or approximate values of the critical numbers. The validity domain of the many asymptotic expansions is then investigated from comparison with the numerical solution. It appears that the small-Re expansion gives good results beyond Re = 1, with an error less that 1%. For Reynolds numbers of a few hundred, we show from the eigenfunctions and the energy equation that the nature of the instability changes: instability still arises from the interfacial mode (there is no mode crossing), but this mode takes all the features of a shear mode. The other modes correspond to the stable eigenmodes of the single-layer Couette flow, which are recovered when one fluid is rigidified by increasing its viscosity or surface tension.     

CFD simulation of two‐ and three‐dimensional free‐surface flow

The paper describes the implementation of moving‐mesh and free‐surface capabilities within a 3‐d finite‐volume Reynolds‐averaged‐Navier–Stokes solver, using surface‐conforming multi‐block structured meshes. The free‐surface kinematic condition can be applied in two ways: enforcing zero net mass flux or solving the kinematic equation by a finite‐difference method. The free surface is best defined by intermediate control points rather than the mesh vertices. Application of the dynamic boundary condition to the piezometric pressure at these points provides a hydrostatic restoring force which helps to eliminate any unnatural free‐surface undulations. The implementation of time‐marching methods on moving grids are described in some detail and it is shown that a second‐order scheme must be applied in both scalar‐transport and free‐surface equations if flows driven by free‐surface height variations are to be computed without significant wave attenuation using a modest number of time steps. Computations of five flows of theoretical and practical interest—forced motion in a pump, linear waves in a tank, quasi‐1d flow over a ramp, solitary wave interaction with a submerged obstacle and 3‐d flow about a surface‐penetrating cylinder—are described to illustrate the capabilities of our code and methods. Copyright © 2003 John Wiley & Sons, Ltd.     

Modeling of long nonlinear waves on the interface in a horizontal two-layer viscous channel flow

The dynamics of two-dimensional waves of small but finite amplitude are theoretically studied for the case of a two-layer system bounded by a horizontal top and bottom. It is shown that for relatively large steady-state flow velocities and at certain fluid depth ratios the vertical velocity profile is nonlinear. An evolutionary equation governing the fluid interface disturbances and allowing for the long-wave contributions of the layer inertia and surface tension, the weak nonlinearity of the waves, and the unsteady friction on all the boundaries of the system is derived. Steady-state solutions of the cnoidal and solitary wave type for the disturbed flow are determined without regard for dissipation losses. It is found that the magnitude and the direction of the flow can alter not only the lengths of the waves but also their polarity.__________Translated from Izvestiya Rossiiskoi Academii Nauk, Mekhanika Zhidkosti i Gaza, No. 1, 2005, pp. 143–158. Original Russian Text Copyright © 2005 by Arkhipov and Khabakhpashev.     

Fully non‐linear free‐surface simulations by a 3D viscous numerical wave tank

A finite difference scheme using a modified marker‐and‐cell (MAC) method is applied to investigate the characteristics of non‐linear wave motions and their interactions with a stationary three‐dimensional body inside a numerical wave tank (NWT). The Navier–Stokes (NS) equation is solved for two fluid layers, and the boundary values are updated at each time step by a finite difference time marching scheme in the frame of a rectangular co‐ordinate system. The viscous stresses and surface tension are neglected in the dynamic free‐surface condition, and the fully non‐linear kinematic free‐surface condition is satisfied by the density function method developed for two fluid layers. The incident waves are generated from the inflow boundary by prescribing a velocity profile resembling flexible flap wavemaker motions, and the outgoing waves are numerically dissipated inside an artificial damping zone located at the end of the tank. The present NS–MAC NWT simulations for a vertical truncated circular cylinder inside a rectangular wave tank are compared with the experimental results of Mercier and Niedzwecki, an independently developed potential‐based fully non‐linear NWT, and the second‐order diffraction computation. Copyright © 1999 John Wiley & Sons, Ltd.     

Nonlinear steady states of hyperelastic membrane tubes conveying a viscous non-Newtonian fluid

We study possible steady states of an infinitely long tube made of a hyperelastic membrane and conveying either an inviscid, or a viscous fluid with power-law rheology. The tube model is geometrically and physically nonlinear; the fluid model is limited to smooth changes in the tube’s radius. For the inviscid case, we analyse the tube’s stretch and flow velocity range at which standing solitary waves of both the swelling and the necking type exist. For the viscous case, we first analyse the tube’s upstream and downstream limit states that are balanced by infinitely growing upstream (and decreasing downstream) fluid pressure and axial stress caused by fluid viscosity. Then we investigate conditions that can connect these limit states by a single solution. We show that such a solution exists only for sufficiently small flow speeds and that it has a form of a kink wave; solitary waves do not exist. For the case of a semi-infinite tube (infinite either upstream or downstream), there exist both kink and solitary wave solutions. For finite-length tubes, there exist solutions of any kind, i.e. in the form of pieces of kink waves, solitary waves, and periodic waves.