Modeling and computation of cavitation in vortical flow

The paper presents an investigation of Euler–Lagrangian methods for cavitating two-phase flows. The Euler–Euler methods, widely used for simulations of cavitating flows in ship technology, perform well in regions of moderate flow changes but fail in zones of strong, vortical flow. Reasons are the strong approximations of cavitation models in the Euler concept. Alternatively, Euler–Lagrangian concepts enable more detailed formulations for transport, dynamics and acoustic of discrete vapor bubbles. Test calculations are performed to study the influence of different parameters in the equations of motion and in the Rayleigh–Plesset equation for bubble dynamics. Results confirm that only Lagrangian models are able to describe correctly the bubble behavior in vortices, while Eulerian results deviate strongly. Lagrangian formulations enable additionally the determination of acoustic pressure of cavitation noise. Two-way coupling between the phases is required for large regions of the vapor phase. A new coupling concept between continuous fluid flow and discrete bubble phase is developed and demonstrated for flow through a nozzle. However, the iterative coupling between the phases via volume fractions is computationally expensive and should therefore be applied only in regions where Eulerian treatment fails. A corresponding local concept for combination with an Euler–Euler method is outlined and is in progress.     

Shock structure in bubbly liquids: comparison of direct numerical simulations and model equations

Shock wave structure in a bubbly mixture composed of a cluster of gas bubbles in a quiescent liquid with initial void fractions around 10% inside a 3D rectangular domain excited by a sudden increase in the pressure at one boundary is investigated using the front tracking/finite volume method. The effects of bubble/bubble interactions and bubble deformations are, therefore, investigated for further modeling. The liquid is taken to be incompressible while the bubbles are assumed to be compressible. The gas pressure inside the bubbles is taken uniform and is assumed to vary isothermally. Results obtained for the pressure distribution at different locations along the direction of propagation show the characteristics of one-dimensional unsteady shock propagation evolving towards steady-state. The steady-state shock structures obtained by the present direct numerical simulations, which show a transition from A-type to C-type steady-state shock structures, are compared with those obtained by the classical Rayleigh–Plesset equation and by a modified Rayleigh–Plesset equation accounting for bubble/bubble interactions in the mean-field theory.      

Free nonlinear oscillations of a thermally relaxing spherical gas bubble in an incompressible liquid

The quasi-adiabatic regime of free oscillation of a bubble in the presence of irreversible interphase heat transfer between the bubble and the ambient liquid is studied. On the basis of simplified model equations of a rarefield bubble mixture, a nonlinear-oscillation equation of the relaxation type is obtained. In constructing an exact particular solution of this equation, the heat transfer law associated with bubble compression is established. For studying the harmonic oscillations, the Krylov-Bogolyubov-Mitropol’skii asymptotic method is used. It is shown that, for a small bubble, the viscosity and heat transfer effects are of the same order. For a small bubble, the influence of these effects on the formation of the natural-oscillation frequency, which is small in the linear approximation, may be significant in the nonlinear formulation. For a large bubble, the influence of these effects is negligible in both approximations. For the approximate solution of the nonlinear equation, a uniformly valid second-order expansion is constructed.     

Bubble dynamics: a review and some recent results

Several aspects of small-amplitude oscillations of bubbles containing gas, vapor, or a gas-vapor mixture are discussed. An application to pressure-wave propagation in a bubbly liquid is described. Nonlinear forced oscillations are considered in the light of recent research on forced oscillations of nonlinear systems. The growth of vapor bubbles, an extension of the Rayleigh-Plesset equation to non-Newtonian liquids and appreciable mass transfer at the interface, and a boundary integral numerical method for nonspherical cavitation bubble dynamics are also briefly discussed.     

The Response of a Rayleigh–Liénard Oscillator to a Fundamental Resonance

A Rayleigh–Liénard oscillator excited by a fundamentalresonance is investigated by using an asymptotic perturbation method based on Fourier expansion and time rescaling. Two first-order nonlinear ordinarydifferential equations governing the modulation of the amplitude andthe phase of solutions are derived. These equations are used todetermine steady-state responses and their stability. Excitationamplitude-response and frequency-response curves are shown and checkedby numerical integration. Dulac's criterion, the Poincaré–Bendixsontheorem, and energy considerations are used in order to study the existenceand characteristics of limit cycles of the two modulation equations. Alimit cycle corresponds to a modulated motion for the Rayleigh–Liénardoscillator. For small excitation amplitude, the analytical results arein excellent agreement with the numerical solutions. In certain caseswhen the excitation amplitude is very low, an approximate analyticsolution corresponding to a modulated motion can be obtained andnumerically checked. Moreover, if the excitation amplitude is increased,an infinite-period bifurcation occurs because the modulation periodlengthens and becomes infinite, while the modulation amplitude remainsfinite and suddenly the attractor settles down into a periodic motion.     

Dynamics of soluble gas bubbles

The problem of the mass, thermal and dynamic interaction between a bubble containing a soluble gas and a liquid is considered. It is shown that this problem can be reduced to the problem of the behavior of a vapor bubble with phase transitions investigated in detail in [1–3]. Expressions are obtained for the rate of decay of the radially symmetric oscillations of the bubbles due to the solubility of the gas in the liquid. The effective coefficients of mass transfer between the radially pulsating bubbles and the liquid are determined. A numerical solution is obtained for the problem of the radial motion of a bubble created by a sudden change of pressure in the liquid which, in particular, corresponds to the behavior of the bubbles behind the shock front when a shock wave enters a bubble screen.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 52–59, November–December, 1985.     

Effect of viscoelastic properties of a liquid on the dynamics of small oscillations of a gas bubble

A series of papers has been devoted to questions of gas bubble dynamics in viscoeiastic liquids. Of these papers we mention [1–4]. The radial oscillations of a gas bubble in an incompressible viscoeiastic liquid have been studied numerically in [1, 2] using Oldroyd's model [5]. Anexact solution was found in [3], and independently in [4], for the equation of small density oscillations of a cavity in an Oldroyd medium when there is a periodic pressure change at infinity. The analysis of bubble oscillations in a viscoeiastic liquid is complicated by properties of limiting transitions in the rheological equation of the medium. These properties are of particular interest for the problem under investigation. These properties are discussed below, and characteristics of the small oscillations of a bubble in an Oldroyd medium are investigated on the basis of a numerical analysis of the exact solution obtained in [3].Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, No. 3, pp. 82–87, May–June, 1976.The authors are grateful to V. N. Nikolaevskii for useful advice and for discussing the results.     

One-dimensional nonstationary motions of a liquid

An effective method is developed for solving the problem of the nonstationary motion of a liquid with plane, cylindrical, and spherical symmetry [1]. It is based on the idea of dividing the region of disturbed motion into two parts and using matched asymptotic expansions. Solutions are presented to typical problems associated with the motion of a piston, and these make it possible to obtain the solution to problems of an explosion in a liquid, oscillations of a bubble, and so forth. It is also shown that the well-known solutions for such problems given, for example, in the book of Naugol'nykh and Roi on the basis of the acoustic approximation with allowance for nonlinear terms are incorrect.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 3–8, March–April, 1980.     

Pulsating flow in a heat exchanger

The problem of heat transfer in industrial processes, heat exchangers, and combustion chambers is formulated for a case where flow inside the chamber consists of a periodic motion imposed on a fully developed turbulent flow. It is shown that the velocity pulsations induce harmonic oscillations in temperature, thus breaking the temperature field into a steady mean part and a harmonic part. The interaction between the velocity and temperature oscillations introduces an extra term into the energy equation which reflects the effect of pulsations in producing higher heat transfer rates. The analysis shows that when the mean temperature is fully developed with constant heat flux at the wall, there is no effect of the velocity pulsations on the total heat transfer rate along the chamber. For the case where the mean temperature profile is not fully developed, analytical solutions are obtained for asymptotic values of the pulsations frequency. The results show the temperature gradient and its dependence on the frequency. These results are used to evaluate the feasibility of pulsating the flow in a heat exchanger for obtaining higher rates of heat transfer.     

Spherical bubble collapse in viscoelastic fluids

The collapse of a spherical bubble in an infinite expanse of viscoelastic fluid is considered. For a range of viscoelastic models, the problem is formulated in terms of a generalized Bernoulli equation for a velocity potential, under the assumptions of incompressibility and irrotationality. The boundary element method is used to determine the velocity potential and viscoelastic effects are incorporated into the model through the normal stress balance across the surface of the bubble. In the case of the Maxwell constitutive equation, the model predicts phenomena such as the damped oscillation of the bubble radius in time, the almost elastic oscillations in the large Deborah number limit and the rebound limit at large values of the Deborah number. A rebound condition in terms of ReDe is derived theoretically for the Maxwell model by solving the Rayleigh–Plesset equation. A range of other viscoelastic models such as the Jeffreys model, the Rouse model and the Doi-Edwards model are amenable to solution using the same technique. Increasing the solvent viscosity in the Jeffreys model is shown to lead to increasingly damped oscillations of the bubble radius.     

Oscillations of a viscous fluid in a thin layer between two coaxial disks rotating together

A study is made of the problem of the motion of an incompressible viscous fluid in the space between two coaxial disks rotating together with constant angular velocity under the assumption that the pressure changes in time in accordance with a harmonic law. The problem is solved using the equations of unsteady motion of an incompressible viscous fluid in a thin layer. It is shown that the velocity field in this case is a superposition on a steady field of damped oscillations with cyclic frequency equal to twice the angular velocity of the disks and forced oscillations with cyclic frequency equal to the cyclic frequency of the oscillations of the pressure field. It is shown that the amplitude of the forced oscillations of the velocity field depends strongly on the ratio of the cyclic frequency of the oscillations of the pressure field to the angular velocity of the disks. It is shown that there is a certain value of the ratio at which the amplitude of the forced oscillations has a maximal value (resonance). It is shown that even for very small amplitudes of the pressure oscillations the amplitude of the oscillations of the relative velocity at resonance may reach values comparable with the mean velocity of the main flow.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 166–169, January–February, 1984.     

Heat exchange between a gas bubble and a liquid

The nonlinear problem of thermal and dynamic interaction between a single gas bubble and surrounding liquid is considered. This problem is met in studies of gas-liquid mixture flows, in particular, in Shockwave propagation in such media. A numerical solution is presented for various modes of bubble surface radial motion. The modes correspond to bubble behavior directly beyond a shock-wave front, where the latter enters the bubble screen, and to the behavior of a bubble located in the depths of the bubble curtain, where the wave becomes diffuse. Analytic solutions of the linearized problem of thermal conductivity for free and constrained small harmonic oscillations of a gas bubble in a liquid were obtained in [1, 2]. Cooling of a hot gas bubble was considered in [3], that study, however, contains inaccuracies. In particular, it was assumed in the solution that the gas density in the bubble was homogeneous. The equation for heat flux in dimensionless variables was written inaccurately. However, in the examples considered in [3] these inaccuracies do not lead to significant errors in the numerical results.     

Weakly nonlinear forced convection patterns in rectangular planform containers

This paper considers the structure of weakly nonlinear steady-state convection patterns in shallow rectangular planform containers heated from below. The lateral dimensions of the container are assumed to be much larger than the characteristic wavelength of convection, and the lateral boundaries are subject to forcing equivalent, for example, to imperfect thermal insulation in the Rayleigh–Benard problem. This has the effect of generating rolls parallel and perpendicular to the lateral boundaries. The resulting patterns are modelled by a coupled pair of nonlinear amplitude equations derived from a phenomenological model of convection introduced by Swift and Hohenberg [Phys. Rev. A15 (1977) 319]. These equations are applicable in the weakly nonlinear limit to a variety of pattern-forming systems such as the Rayleigh–Benard system. Solutions are found using both numerical and asymptotic methods. The boundary imperfection is shown to give rise to some novel effects, including the possibility of patterns containing square cells. More generally, patterns evolve that are dominated by rolls but with transitions to more complex bimodal forms near the edges of the container. The emergence and structure of transition lines, or grain boundaries, is analysed in detail.     

Conjugate heat transfer inside a porous channel

Analytical and numerical analyses have been performed for fully developed forced convection in a fluid-saturated porous medium channel bounded by two parallel plates. The channel walls are assumed to be finite in thickness. Conduction heat transfer inside the channel wall is also accounted and the full problem is treated as a conjugate heat transfer problem. The flow in the porous material is described by the Darcy–Brinkman momentum equation. The outer surfaces of the solid walls are treated as isothermal. A temperature dependent volumetric heat generation is considered inside the solid wall only. Analytical expressions for velocity, temperature, and Nusselt number are obtained after simplifying and solving the governing differential equations with reasonable approximations. Subsequent results obtained by numerical calculations show an excellent agreement with the analytical results.     

Nonlinear interaction of longitudinal—Transverse acoustic waves in a circular cylinder

The limiting amplitudes of acoustic oscillations in a cylindrical volume of a heat releasing medium in which one or several modes are unstable in the linear approximation are determined. One of the mechanisms limiting the amplitudes of unstable acoustic modes is the transfer of energy from them to damped modes by nonlinear interaction. The nonlinear interactions of plane acoustic waves in a long channel have been considered by Artamonov and Vorob'ev [1]; in the present paper, the interaction of mixed longitudinal—transverse acoustic modes in a closed cylindrical volume is considered. The equations describing the interaction of two and three longitudinal—transverse modes are derived and investigated in the quadratic approximation by the method of slowly varying amplitudes and phases of the oscillations [2]. The treatment is applicable to a high-temperature gas, for which general stability conditions in the linear approximation have been formulated by Artamonov [3].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 3–9, September–October, 1982.I should like to express my thanks to K. I. Artamonov (deceased) for suggesting the problem and for scientific supervision and A. P. Vorob'ev for constant interest in the work and helpful advice.     

Dynamics of a Flexible Cantilever Beam Carrying a Moving Mass

The motion of a flexible cantilever beam carrying a moving spring-mass system is investigated. The beam is assumed to be an Euler–Bernouli beam. The motion of the system is described by a set of two nonlinear coupled partial differential equations where the coupling terms have to be evaluated at the position of the mass. The nonlinearities arise due to the coupling between the mass and the beam. Due to the nonlinearities the system exhibits internal resonance which is investigated in this work. The equations of motion are solved numerically using the Rayleigh–Ritz method and an automatic ODE solver. An approximate solution using the perturbation method of multiple scales is also obtained.     

Convection of a viscous incompressible gas in rectangular regions with inlet and outlet channels

A study is made of two-dimensional problems of thermal convection of a viscous incompressible gas in rectangular regions that have gas inlet and outlet channels in the presence of a temperature difference between the bottom and the top (the bottom is heated). In contrast to the well-studied case of natural convection, when no-slip conditions are specified on all boundaries of the region and motion in the region occurs only through the temperature difference [1–4], the heat transfer in the investigated flows is complicated by the additional influence of the forced convection of the gas due to the motion of gas through the inlet and outlet channels. Flows of such type simulate well the processes that take place in many heat transfer devices and in ventilated and air-conditioned industrial premises. Two formulations of the problem are considered. In the first, the gas flow through the inlet and outlet channels is assumed given, and the solution of the problem is determined by the dimensionless Prandtl, Grashof, and Reynolds numbers. In the second case, this flow rate is not given but determined during the solution of the problem. The motion in the region arises from the difference between the temperatures of the bottom and the top of the region, and the motion, in its turn, causes a flow of gas through the inlet and outlet channels. As in the case of natural convection, the solution of the problem in this case is determined by only two dimensionless numbers — the Grashof and Prandtl numbers. By numerical solution of the boundary-value problems for the equations of heat transfer a study is made of the influence of the characteristic dimensionless numbers on the hydrodynamic and temperature fields and the heat fluxes through the boundaries of the region. The solutions of the problems in the two formulations are compared for different positions of the outlet channels.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 126–131, September–October, 1979.We thank G. I. Petrov for discussing the results.     

Nonlinear Oscillations of Viscoelastic Rectangular Plates

Nonlinear oscillations of viscoelastic simply supported rectangular plates are studied by assuming the Voigt–Kelvin constitutive model. Using Hamilton's principle in conjunction with the kinematics associated with Kirchhoff's plate model, the governing equations of motion including the effect of damping are represented in terms of the transversal deflection and a stress function. Utilizing the Bubnov–Galerkin method, the nonlinear partial differential equations are reduced to an ordinary differential equation which is studied geometrically by approximate construction of the Poincaré maps. Explicit expressions are given for periodic solutions.     

Dynamics of vapor-gas bubbles

The nonlinear problem of thermal, mass, and dynamic interaction between a vapor-gas bubble and a liquid is considered. The results of numerical solution of the problem of radial motion of the bubble caused by a sudden pressure change in the liquid, illustrating the behavior of vapor-gas bubbles in compression and rarefaction waves, are presented. The corresponding problem for single-component gas and vapor bubbles was considered in [1, 2].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 56–61, November–December. 1976.