A two-parameter model of wave regimes for viscous liquid film flows

Wave regimes of viscous liquid film flows are considered when the viscosity coefficients vary in a wide range. An approximate model system of differential equations with two external governing parameters for the film layer thickness and the local flow rate is derived. The viscous dissipation of a film layer is taken into account in this system more accurately than in the well-known one-parameter Shkadov model. New properties of linear and nonlinear waves caused by the hydrodynamic instability of high-viscous liquid flows under gravity and surface tension are found.     

Wave motion in a gravity field on the free surface and stratification interface of a fluid with stratified inhomogeneity. Nonlinear analysis

Regularities of the nonlinear gravitational wave motion in a two-layer density-stratified fluid are investigated for a finite thickness of the upper, lighter, layer. The characteristics of the nonlinear internal resonant interaction of the gravity waves generated by the free surface of the upper layer and the medium interface are considered. It is shown that in second-order calculations both degenerate (two-wave) and secondary combined (three-wave) resonant interactions may be realized.     

矩形弹性壳液耦合系统中的重力波分析

根据非线性动力学理论，建立了矩形壳液耦合系统的非线性振动方程组，通过数值求解，发现当激振频率为壳体固有频率与重力波频率之和，且激振力足够大时，会产生大幅低频重力波，通过实验验证，发现了壳液耦合系统中存在的大幅低频重力波现象，实验结果与理论结果基本吻合。     

The analysis of a mechanism of liquid replenishment and draining in horizontal two-phase flow

There is a regime of two-phase flow in which large waves or surges pass rapidly along a horizontal tube accompanied by splashing, wave-breaking and entrainment with the result that water is thrown to the upper surface of the tube. Between surges the film on the top surface is depleted by draining under gravity and by evaporation if the tube is heated. If the interval between surges is sufficiently long a dry patch may begin to form. In this paper, theory is given for the calculation of the film thickness left behind on the top surface and for the calculation of the time to dryout. The theory includes both the effect of the boundary layer developement during replenishment of the film and also the effect of the axial deceleration of the film at the point where the liquid replenishment ceases. Finally, the predicted variation of film thickness is compared with experimental film thickness traces obtained in this type of horizontal two-phase flow. The agreement is found to be very satisfactory. This analysis is of interest in connection with the prevention of permanent and intermittent dryout at low qualities in nuclear power station evaporators.     

Nonlinear vibrations of a vapor film in the presence of intense heat fluxes

Waves propagating along the interface between a thin vapor film and a liquid layer in the presence of a heat flux are investigated. The boundary conditions on the vapor-liquid phase surface take into account the temperature dependence of the pressure and the possibilities of formation of the metastable state of the superheated liquid and mass flow. Variations in the saturation pressure as functions of the temperature and mass flux lead to generation of weakly damped periodic waves of low amplitude whose velocity can be much higher than the velocity of the gravity waves. The waves ensure stability of the vapor film beneath the liquid layer in the gravity field. The finite-amplitude waves on the surface of the vapor film differ from the Stokes surface waves on the free surface of isothermal fluid. Instability regimes related with superheating of the liquid ant its explosive boiling when the amplitude of an initially small wave increases to infinity in a finite time can develop in a certain working-parameter regime.     

Stability of a Liquid Film Flowing Down an Oscillating Inclined Surface

The stability of flow of a liquid film along an inclined plate subject to periodic oscillations under the action of the gravity force is investigated with allowance for the surface tension. An equation of the Orr-Sommerfeld type with time-periodic coefficients is used. A method for determining the eigenvalues of the linear stability problem is developed on the basis of Floquet theory, spectral representation of the variables, and multistep methods of integration of ordinary differential equations. The bifurcation spectrum of the resonance modes is investigated, and the amplification coefficients and phase velocities are calculated for the surface waves, Tollmien-Schlichting waves, and resonance waves. The influence of external parameters, namely, the inclination, the surface tension, and the layer thickness, on the resonance modes and the steady-state flow modes is studied.     

Wave-flow theory for a thin viscous liquid layer

On the basis of a simplified system of equations we study the process of development and stability of wave flows in a thin layer of a viscous liquid. Any unstable disturbance of the laminar flow grows and leads to the establishment of the wave regime. The time to establish the flow changes little for large flow rates, but increases sharply with reduction of the flow rate. Given the same amplitudes of the initial disturbances, the optimum regimes which provide the greatest flow rate in a layer of given average thickness develop more rapidly than the other regimes. All the wave regimes are unstable to disturbances in the form of traveling waves. With moderate flow rates, the optimum regimes will be most stable to near-by disturbances.Strictly periodic wave flows in a thin layer of a viscous liquid under the influence of the gravity force were calculated in [1], Various flow wave regimes which differ in wavelength can theoretically be established for a given liquid flow rate. In particular, there is a wavelength for which the flowing layer exhibits minimum average thickness (and maximum flow rate for a given average thickness). These optimum regimes correspond closely to the experimental data [2]; however, with regard to calculation technique these regimes are no different from the others. In the following we make a comparison of the wave regimes on the basis of the nature of their development and stability.     

Effect of permeability on the instability of a non-Newtonian film down a porous inclined plane

A thin film of a power–law fluid flowing down a porous inclined plane is considered. It is assumed that the flow through the porous medium is governed by the modified Darcy’s law together with Beavers–Joseph boundary condition for a general power–law fluid. Under the assumption of small permeability relative to the thickness of the overlying fluid layer, the flow is decoupled from the filtration flow through the porous medium and a slip condition at the bottom is used to incorporate the effects of the permeability of the porous substrate. Applying the long-wave theory, a nonlinear evolution equation for the thickness of the film is obtained. A linear stability analysis of the base flow is performed and the critical condition for the onset of instability is obtained. The results show that the substrate porosity in general destabilizes the film flow system and the shear-thinning rheology enhances this destabilizing effect. A weakly nonlinear stability analysis reveals the existence of supercritical stable and subcritical unstable regions in the wave number versus Reynolds number parameter space. The numerical solution of the nonlinear evolution equation in a periodic domain shows that the fully developed nonlinear solutions are either time-dependent modes that oscillate slightly in the amplitude or time independent stable two-dimensional nonlinear waves with large amplitude referred to as ‘permanent waves’. The results show that the shape and the amplitude of the nonlinear waves are strongly influenced by the permeability of the porous medium and the shear-thinning rheology.     

Breaking of waves of limiting amplitude over an obstacle

Homogeneous heavy fluid flows over an uneven bottom are studied in a long-wave approximation. A mathematical model is proposed that takes into account both the dispersion effects and the formation of a turbulent upper layer due to the breaking of surface gravity waves. The asymptotic behavior of nonlinear perturbations at the wave front is studied, and the conditions of transition from smooth flows to breaking waves are obtained for steady-state supercritical flow over a local obstacle. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 47, No. 3, pp. 3–11, May–June, 2006.     

TE modes of a layered nonlinear optical waveguide

The general theory is given of the propagation of electromagnetic TE modes in a layered structure consisting of a homogeneous dielectric layer between two homogeneous optically nonlinear semi-infinite dielectrics. Allowance is made for the nonlinear coefficients of the nonlinear media to be different and a detailed analysis is given, together with the appropriate formulae for the power flow down the guide, as a function of guide wave number, and layer thickness. Detailed consideration is given to the existence conditions for surface waves.     

Long time interaction of envelope solitons and freak wave formations

This paper concerns long time interaction of envelope solitary gravity waves propagating at the surface of a two-dimensional deep fluid in potential flow. Fully nonlinear numerical simulations show how an initially long wave group slowly splits into a number of solitary wave groups. In the example presented, three large wave events are formed during the evolution. They occur during a time scale that is beyond the time range of validity of simplified equations like the nonlinear Schrödinger (NLS) equation or modifications of this equation. A Fourier analysis shows that these large wave events are caused by significant transfer to side-band modes of the carrier waves. Temporary downshiftings of the dominant wavenumber of the spectrum coincide with the formation large wave events. The wave slope at maximal amplifications is about three times higher than the initial wave slope. The results show how interacting solitary wave groups that emerge from a long wave packet can produce freak wave events.Our reference numerical simulation are performed with the fully nonlinear model of Clamond and Grue [D. Clamond, J. Grue, A fast method for fully nonlinear water wave computations, J. Fluid Mech. 447 (2001) 337–355]. The results of this model are compared with that of two weakly nonlinear models, the NLS equation and its higher-order extension derived by Trulsen et al. [K. Trulsen, I. Kliakhandler, K.B. Dysthe, M.G. Velarde, On weakly nonlinear modulation of waves on deep water, Phys. Fluids 12 (10) (2000) 2432–2437]. They are also compared with the results obtained with a high-order spectral method (HOSM) based on the formulation of West et al. [B.J. West, K.A. Brueckner, R.S. Janda, A method of studying nonlinear random field of surface gravity waves by direct numerical simulation, J. Geophys. Res. 92 (C11) (1987) 11 803–11 824]. An important issue concerning the representation and the treatment of the vertical velocity in the HOSM formulation is highlighted here for the study of long-time evolutions.     

Turbulent breaking of overturning gravity waves below a critical level

The interaction of an internal gravity wave with its evolving critical layer and the subsequent generation of turbulence by overturning waves are studied by three-dimensional numerical simulations. The simulation describes the flow of a stably stratified Boussinesq fluid between a bottom wavy surface and a top flat surface, both without friction and adiabatic. The amplitude of the surface wave amounts to about 0.03 of the layer depth. The horizontal flow velocity is negative near the lower surface, positive near the top surface with uniform shear and zero mean value. The bulk Richardson number is one. The flow over the wavy surface induces a standing gravity wave causing a critical layer at mid altitude. After a successful comparison of a two-dimensional version of the model with experimental observations (Thorpe [21]), results obtained with two different models of viscosity are discussed: a direct numerical simulation (DNS) with constant viscosity and a large-eddy simulation (LES) where the subgrid scales are modelled by a stability-dependent first-order closure. Both simulations are similar in the build-up of a primary overturning roll and show the expected early stage of the interaction between wave and critical level. Afterwards, the flows become nonlinear and evolve differently in both cases: the flow structure in the DNS consists of coherent smaller-scale secondary rolls with increasing vertical depth. On the other hand, in the LES the convectively unstable primary roll collapses into three-dimensional turbulence. The results show that convectively overturning regions are always formed but the details of breaking and the resulting structure of the mixed layer depend on the effective Reynolds number of the flow. With sufficient viscous damping, three-dimensional turbulent convective instabilities are more easily suppressed than two-dimensional laminar overturning.     

Wave regimes of ascending flow of a thin layer of a viscous liquid in contact with a gas

The flow of a liquid in thin layers is one of the hydrodynamic problems of chemistry and heat engineering. The large surface area of films and their small thickness make it possible to accelerate thermal, diffusive, and chemical processes at the gas-liquid boundary.Theoretical studies of liquid flow in a vertical descending thin layer are presented in [1–4]. In this paper we study ascending wave flows of a liquid in a thin vertical layer in contact with a gas, i.e., flows in the direction opposite the action of the force due to gravity, with account for the action of the gas on the liquid surface. Such motions are encountered when oil is extracted from strata that are saturated with gas. At some distance from the stratum the oil and gas separate: the gas travels at high velocity inside the pipe, occupying a considerable portion of the pipe, and the liquid is displaced toward the pipe walls, forming a thin film. In certain cases a wave-like interface develops between the oil and gas that travels with a velocity greater than that of the liquid but less than the average gas velocity. Similar phenomena are observed in high velocity mass exchangers.We examine the effect of the gas for both laminar and turbulent flow.Studies that neglect the effect of the gas flow on the liquid show that for waves on the film surface whose lengths are considerably longer than the average thickness of the layer, the liquid motion in the film is described by boundary layer equations in which account is taken of the mass force, i.e., the force due to gravity. With some approximation, we can assume that in accounting for the effect of the gas on the liquid the liquid flow is described by these same equations.     

A novel method of solving the exact equations for capillary-gravity waves

A new method is presented for the computation of two-dimensional periodicprogressive surface waves propagating under the combined influence of gravity and surfacetension.The nonlinear surface is expressed by Fourier series with finite number of terms,after the computational domain is transformed into a unit circle.The dynamic boundaryequation is used in its exact nonlinear form and the coefficients of Fourier series are foundby the Nweton-Raphson method successively.This is a neat method,Yielding highprescision with little computational effort.     

Nonstationary motion of a shell on the surface of a heavy fluid

A three-dimensional nonstationary problem of vibrations of a flexible shell moving on the surface of an ideal heavy fluid. The forces due to surface tension are ignored. The problem is formulated in the space of the acceleration potential. The potential of the pulsating source is found by solving the Euler equation and the continuity equation taking into account the free-surface conditions (linear theory of small waves) and the conditions at infinity. The density distribution function of the dipole layer is determined from the boundary conditions on the surface of the shell. Formulas for determining the shape of gravity waves on the fluid surface and the natural frequencies of vibrations of the shell are obtained. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 66–75, July–August, 2009.     

Generation of Nonlinear Waves on a Viscoelastic Coating in a Turbulent Boundary Layer

Self–induced excitation of periodic nonlinear waves on a viscoelastic coating interacting with a turbulent boundary layer of an incompressible flow is studied. The response of the flow to multiwave excitation of the coating surface is determined in the approximation of small slopes. A system of equations is obtained for complex amplitudes of multiple harmonics of a slow (divergent) wave resulting from the development of hydroelastic instability on a coating with large losses. It is shown that three–wave resonant relations between the harmonics lead to the development of explosive instability, which is stabilized due to the deformation of the mean (Sover the wave period) shear flow in the boundary layer. Conditions of soft and hard excitation of divergent waves are determined. Based on the calculations performed, qualitative features of excitation of divergent waves in known experiments are explained.     

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.     

Gravity waves and Rayleigh Taylor instability on a Jeffrey-fluid

A theory for linear surface gravity waves on a semi-infinite layer of viscoelastic fluid described by a Jeffrey model is presented. Results are given for the decay rate and the phase velocity as a function of the parameters of the fluid: a nondimensional time constant, and a ratio of the retardation time to the relaxation time. At small wave numbers the behavior is Newtonian. In other cases depending on the nondimensional parameters, a number of possible other behaviors exist. The influence of the non-dimensional parameters on the growth rate of Rayleigh-Taylor instability is also discussed.     

高度非线性孤立波与弹性大板的耦合作用研究

一维颗粒链的一端受到一个有初速度颗粒的撞击,导致颗粒连中产生稳定传播的应力波——高度非线性孤立波,该应力波的波长、波速以及幅值都能保持很好的稳定性,且遇到边界才会反射. 孤立波是一种良好的信息载体,广泛应用于无损检测技术中. 基于孤立波的特性,研究高度非线性孤立波与弹性大板耦合作用,基于赫兹定律和板的内在非弹性理论,推导出晶体链与大板的耦合微分方程组. 用龙格库塔法求解该微分方程组,得到颗粒链中各颗粒的位移、速度曲线. 通过分析回弹波出现的时间、回弹波所携带的能量以及模量、厚度、重力等对孤立波的影响,发现反射孤立波对大板的弹性模量和厚度尤为敏感,此外,颗粒链的摆放对整个耦合过程也有影响. 研究的结果为孤立波对结构体的无损探伤提供了理论依据,该技术可实现对结构体的快速检查和可控性研究.     

Influence of rotation on the damping of surface waves in a cylindrical layer of viscous fluid

A study is made of the oscillations of the free surface of a cylindrical layer of viscous incompressible fluid attracted by gravity to a solid, uniformly rotating cylinder. Logarithmic decay rates are found for the damping of the surface gravitational waves when large Reynolds numbers are assumed. It is shown that rotation introduces an asymmetry into the damping of the waves traveling in and against the direction of rotation of the fluid.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6. pp. 171–173, November–December, 1984.