Small amplitude pressure wave in subsonic and supersonic bubble flows in convergent-divergent nozzles

This paper makes a theoretical analysis of the propagation phenomena of the small amplitude pressure wave in the subsonic and supersonic bubble flow with a velocity slip between bubble and liquid in the convergent-divergent nozzle. From an analysis of the time-mean flow, the nondimensional parameter $\text{m = {u}\text{2}\text{G}\text{·α(1 ? α)ρ}\text{l}\text{β(2 ? 1/S)/P·[αβS + (1 ? α)βS}{}^2\text{+ α(1 ? α)]}}{}^{\mathrm{<mtext>1</mtext><mtext>2</mtext>}}$ corresponds to Mach number is gasdynamics where u_G is the gas velocity, α: the void fraction, ρ_L: the liquid density, P: the pressure, S: the velocity ratio of the gas and liquid flows and β: the proportional constant for the virtual mass. From a theoretical analysis of the small disturbance field, it is clarified that the parameter m also plays an essential and important role as Mach number, although the propagation performance of the disturbance is very complicated compared with that in gasdynamics. It is also shown that the pressure waves are divided into four groups depending on the velocity ratio S. Two of them are rather realistic, but the other two are required of a further investigation in future.     

Experimental and numerical investigation a bubble-driven laminar flow in a axisymmetric vessel

Measurements have been obtained, by laser-Doppler anemometry (LDA), of the axisymetric, recirculating liquid flow caused by a column of air bubbles ($\text{5–6}\text{1}\text{2}\text{mm dia.}$) rising through caster oil in a cylindrical enclosure (100 mm dia.). The liquid velocities correspond to creeping flow. Axial and radial liquid velocity profiles are reported at eight axial stations and, close to within the bubble column, as a function of time. The maximum liquid velocity found outside the bubble column is about 0.5 of that of the bubbles and a very rapid radical decay from this value is noted. The temporal variation of the velocity field, due to the passage of the air bubbles, is undetectable at radial locations greater than about $\text{1}\text{1}\text{2}$ bubble radii from the centreline.The variation of bubble velocity with axial distance was aise measured by LDA for liquid height to enclosure diámeter ratios of 0.98 and 2.78. The maximum bubble velocities were about 0.1–0.2 higher than the Strokes law terminal velocity. The increase is due to the convection of the bubble column by the liquid flow. The maximum bubble velocity is established within approximately three bubble diameters of the air inlet.The motion of the liquid has been calculated by the numerical solution of the steady form of the equations of motion, with the inner boundary of the area of integration lying 1.3 bubble radii from the centerline. The boundary conditions at this surface are assumed to be steady and are taken from measurements of the time-averaged velocity components. The assumption of steady flow at this boundary is supported by experimental observation and results in calculations which are generally in close agreement with the measurements. Discrepancies are confined to the immediate vicinity of the bubble column near to the top and bottom of the enclosure. These are ascribed to a combination of small asymmetries in the experiment and inadequate numerical resolution in these regions.     

Rheology of suspensions of spherical particles in a newtonian and a non-newtonian fluid

Properties of suspensions of spherical glass beads (25–38 μm dia.) in a Newtonian fluid and a non-Newtonian (NBS Fluid 40) fluid were measured at volume fractions, φ, of 0%, 10%, 20% and 30%. Measurements were made using a modified and computerized Weissenberg Rheogoniometer. Properties measured included steady shear viscosity, η($\text{γ}\text{.}$), first normal stress difference, N₁($\text{γ}\text{.}$), linear viscoelastic properties, η′(ω) and G′(ω), shear stress relaxation, σ^? ($\text{γ}\text{.}$, t), and growth, σ⁺($\text{γ}\text{.}$, t) and normal stress relaxation, N₁^?($\text{γ}\text{.}$, t).For a the Newtonian fluid, increasing φ causes both η and η′ to increase, with η′ showing a slight frequency dependence. Both N₁ and G′ are zero and stress relaxation and growth occur essentially instantaneously. For the NBS fluid, both η and η′ increse with φ at all $\text{γ}\text{.}$ and ω, respectively, the increase being greater as $\text{γ}\text{.}$ and ω approach zero. N₁ and G′ are less affected by the presence of the particles than η and η′ with the effect on G′ being more pronounced than on N₁. For fixed $\text{γ}\text{.}$, stress relaxation and growth exhibit greater non-linear effects as φ is increased. A model for predicting a priori the linear viscoelastic properties for suspensions was found to yeild reasonable estimates up to φ = 20%.     

Non-linear boundary and eigenvalue problems for the emden-fowler equations by group methods

Two pragmatic boundary value and eigenvalue problems of the Emden-Fowler equation $\text{(t}{}^{\mathrm{\alpha }}\text{u′)′ + λt}{}^{\mathrm{\beta }}\text{?(u) = 0,?(u) = u}{}^{\mathrm{\gamma }}$ and e^u are studied using the simple one parameter group properties. In all cases boundary value problems are converted into initial value problems using the property of the invariance group. With $\text{?(u) = u}{}^{\mathrm{\gamma }}$ an eigenvalue problem is detailed and calculations presented.     

Correlation of turbulent flow rate-pressure drop data for non-newtonian solutions and slurries in pipes

Isothermal and non-isothermal flow rate-pressure drop data in turbulent flow through smooth pipes have been obtained for non-Newtonian fluids, including aqueous solutions of polymers and aqueous suspensions of titanium dioxide. It has been found that the friction factor, f, is a function of a new form of Reynolds number, Re_B, based on the parameters A, x and w of Bowen's correlation, viz.

$\text{τ}{}_{\mathrm{w}}\text{D}{}^{\mathrm{x}}\text{=A}\text{u}{}^{\mathrm{w}}$

where τ_w is the wall shear strees, ?u the mean velocity, D the pipe diameter; A, x and w are experimentally derived parameters which characterise the fluid.     

Local properties of vertical mercury-argon two-phase flow in a circular tube under transverse magnetic field

Experimental data are presented in this paper on the profiles of local void fraction, bubble impaction rate, bubble velocity and its spectrum, and also bubble length and its spectrum, of mercury-argon two-phase slug flow flowing upwards in a vertical circular tube in the presence of a transverse magnetic field. Decrease in void fraction and increase in bubble velocity are significant when the magnetic flux density is larger than $\text{0.3～0.4}$$\text{T}\text{(Ha ? 100)}$. This effect is discussed by analyzing the bubble size distribution. Recovery of local void fraction profile in the downstream of an obstacle and diffusion of void injected from only one nozzle in the presence of magnetic field are also discussed.     

A critical re-appraisal of the jet-thrust technique for normal stresses,with particular reference to axial velocity and stress rearrangement at the exit plane

It is shown that the thrust, T, exerted by a jet on the tube from which it flows, and the corresponding die-swell ratio, D, are closely related and dependent on the axial velocity and stress profiles at the exit plane. Velocity-profile data, calculated by Tanner using a finite element method, have been used to demonstrate that for a Newtonian liquid the reduction in measured thrust from the expected value arises from a re-arranged, non-parabolic axial velocity profile and the related re-arranged non-zero axial stress profile at the exit plane. The axial stress re-arrangement is the major effect.Using the correction-curve thus derived to determine the normal stresses, ν₁ + $\text{1}\text{2}$ν₂ aqueous and non-aqueous polymer solutions gives values that are higher than the “correct” results by a significant, substantial amount. The difference is not due to neglect of the second normal stress difference, ν₂, nor to the neglect of the wall pressure at the exit plane, which is shown experimentally to be very small. It is suggested that the difference, which is a function only of the shear stress (or rate of shear) at the wall, may arise from a difference in the stress profile associated with the velocity re-arrangement at the exit between Newtonian liquids and elasticoviscous liquids for which the extensional viscosity may be high.     

Prediction of mono-disperse gas–liquid turbulent flow in a vertical pipe

A two-fluid model in the Eulerian–Eulerian framework has been implemented for the prediction of gas volume fraction, mean phasic velocities, and the liquid phase turbulence properties for gas–liquid upward flow in a vertical pipe. The governing two-fluid transport equations are discretized using the finite volume method and a low Reynolds number $k-\varepsilon$ model is used to predict the turbulence field for the continuous liquid phase. In the present analysis, a fully developed one-dimensional flow is considered where the gas volume fraction profile is predicted using the radial force balance for the bubble phase. The current study investigates: (1) the turbulence modulation terms which represent the effect of bubbles on the liquid phase turbulence in the k–ε transport equations; (2) the role of the bubble induced turbulent viscosity compared to turbulence generated by shear; and (3) the effect of bubble size on the radial forces which results in either a center-peak or a wall-peak in the gas volume fraction profiles. The results obtained from the current simulation are generally in good agreement with the experimental data, and somewhat improved over the predictions of some previous numerical studies.     

Espaces d'écoulements dits « universels »Spaces of universal flows

An isochoric motion can be performed both in perfect fluid, in Newtonian fluid, in Maxwell fluid (slow motions) and in Rivlin–Ericksen fluid of second grade whatever be viscosities and viscometric coefficients, iff the motion is universal. Every universal motion with steady vorticity is a generalised Belrami flow, and fulfils the Stokes equation. If the velocity $\text{u}$ of an universal motion complies with $\text{rot}\text{[(?}{}^{\mathrm{<mtext>t</mtext>}}\text{(}\text{Δ}\text{u}\text{))}\text{u}\text{]=}\text{0}$, the motion stands for feasible motion in every second order fluid. Brothers of the potential flows, all the sets of universal motions make up bundles of linear or cono??d spaces with various dimensions, finite or infinite, issued from the rest $\text{u}\text{≡}\text{0}$. The structures appear by scanning parallel to the potential flows. To cite this article: M. Bouthier, C. R. Mecanique 331 (2003).     

Some remarks on “Drag coefficients of a slowly moving sphere in Non-Newtonian fluids”

Viscoelastic solutions were ejected vertically downwards into air and various Newtonian fluids. The measured swell increased significantly when ejected into a liquid rather than air. The observed increase is considered a result of both bouyancy and drag forces on the solution. The following dimensions expression relating the ratio of the swell diameter in liquid and air D_L/D_A to the elastic shear compliance of the ejected solution J_e was experimentally observed.(D_L/D_A)⁶-1=30(Δ?/?_s)^{?$\text{1}\text{2}$}([g²η²_N?_s]$\text{1}\text{3}$J_e)^{$\text{3}\text{5}$}, where Δ? is the density difference between the extruded and Newtonian fluid, ?_s is the solution density, g is the gravitational constant, and η_N is the Newtonian fluid viscosity. Thus with this expression a simple extrudate swell technique exists to estimate the elastic shear compliance of a viscoelastic solution.     

Dynamics of a system exhibiting the global bifurcation of a limit cycle at infinity

This paper concerns the dynamics of a class of non-linear oscillators of the form:

$\text{x″ + x ? ?x′(1?ax}{}^2\text{?bx′}{}^2\text{) = 0}$

. The non-linear term contains two parameters a and b which may be varied to give the Rayleigh and Van der Pol differential equations as special cases.The existence and approximation of limit cycles in this system are investigated using the Poincare-Bendixson theorem and the Lindstedt perturbation method. Analysis of the system at infinity is used to study the global bifurcation through which the limit cycle is created from four saddle-saddle connections between equilibrium points at infinity. Center manifold theory is used to determine the stability of the equilibrium points at infinity. Numerical integration is used to verify the analytical results.It is shown that an arbitrarily small perturbation to the damping term of the Rayleigh equation results in points close to the stable limit cycle escaping to infinity.     

A class of complete integrals of plane eikonal equations

A class of complete integrals of the plane eikonal equation

$\text{|}\text{grad}\text{u|}{}^2\text{= f(x,y)}$

for harmonic u(x,y) is determined by using complex variables. The case in which z = 0 is a singular point of the analytic function whose real part is log f is also treated. Illustrative examples are given.     

The role of shear stresses in brittle fracture

Using the three-dimensional model for brittle fracture developed earlier by S.A.F. Murrell and P.J. Digby (1970,1972) shear stress concentrations are calculated for brittle bodies and the relative roles of tensile and shear stresses in the fracture process are considered. It is found that the maximum shear stress and the maximum tensile stress occur at different places on a crack, and that there is a wide range of stress states for which they do not occur on the same crack. Furthermore, if the theoretical cleavage strength is σ_max and the theoretical shear strength is τ_max, then cleavage precedes inelastic shear and brittle fracture is possible, for suitable stress systems, when $\text{σ}{}_{\max }\text{<}\text{2τ}{}_{\mathrm{<mtext>max</mtext>}}\text{(1 ? ν)}$, where ν is the Poisson's ratio of the solid matrix. This appears to be in accordance with empirical observations.     

Investigation of droplet deposition from a turbulent gas stream

Turbulent deposition of particles from two-phase flow onto the smooth wall of a tube has been studied theoretically and experimentally. A model is proposed for the deposition motion of large particles based on turbulent diffusion in the core followed by a free flight towards the wall. The theory shows that within the Stokes regime, the dimensionless deposition velocity $\text{k-}{}_{\mathrm{<mtext>d</mtext>}}\text{/}\text{u}$* depends on Re and $\text{τ}{}^{\mathrm{+}}$ only, where u* is the friction velocity, Re is the tube Reynolds number and $\text{τ}{}^{\mathrm{+}}$ is the dimensionless particle relaxation time. Deposition data are obtained for air-water droplet flow through a 12.7-mm i.d. acrylic tubing at $\text{Re}\text{= 52,500}$ and 94,600. The proposed theory satisfactorily describes the existing deposition data as well as present measurements, covering a wide range of Re and $\text{τ}{}^{\mathrm{+}}$.     

Shock wave-induced bubble motion

There is a lack of data on transient interphase shear forces on bubbles. The present work is an experimental and analytical investigation of unsteady bubble motions induced by shock waves. Bubble trajectories were monitored in a shock tube using high-speed Schlieren photographs. Pressure histories at two locations weere used to determine the average wave speed in the bubbly mixture. An equation of motion was developed and solved for predicting the bubble trajectories and velocities. Reasonable agreement is obtained between the experimental results and the analytical predictions. It was found that for experimental range of Reynolds numbers (200 – 12,000), a constant drag coefficient $\text{c}{}_{\mathrm{D}}\text{? 2.6}$ can be utilized to estimate the drag force.     

Complexification phenomenon in an example of sensitive singular perturbationPhénomène de complexification dans un exemple de perturbation singulière sensitive

We consider singular perturbation problems depending on a parameter ε?0 such that for ε>0 the solution u^ε belongs to a Sobolev space on a domain $\text{Ω}$, but the limit u⁰ is not a distribution on $\text{Ω}$. A very simple model problem, solvable by Fourier transform allows us to study the complexification process of u^ε as $\text{ε↘0}$. The limit holds in the topology of a space of analytical functionals. To cite this article: C.A. De Souza, É. Sanchez-Palencia, C. R. Mecanique 332 (2004).     

Turbulence structure of air-water bubbly flow—III. transport properties

An experimental study of heat and bubble transport in turbulent air-water bubbly flow was carried out by means of tracer techniques. Helium tracer gas concentration data and temperature distributions were used to extract bubble and heat diffusivity information. The results indicated that the turbulent velocity components of the liquid phase play a predominant role in the turbulent transport process. A systematic increase of diffusivity of heat, ?_H, with quality and water velocity was observed. An empirical correlation for the diffusivity ratio ?_H,TP/?_H,SP is presented. The Péclet number, u′_covbar|dφ, for bubble dispersion can be approximated by 2.0, independent of the flow variables. The bubble-to-heat diffusivity ratio, φ/?_H, approaches unity with increasing quality and water velocity. Momentum transport is also discussed, based on a mixing length theory.