Lift force acting on single bubbles in linear shear flows

Lift coefficients, C_L, of single bubbles in linear shear flows are measured to investigate effects of the bubble shape, the liquid velocity gradient and the fluid property on C_L. The range of the Morton number, M, tested is from logM = − 6.6 to − 3.2. The shapes of bubbles are spherical and ellipsoidal. A correlation of bubble aspect ratio for single bubbles in infinite stagnant liquids proposed in our previous study can give good evaluations for bubbles in the linear shear flows. The C_L of spherical bubbles at low bubble Reynolds numbers, Re, depend on the dimensionless shear rate Sr and Re and decrease with increasing Re. These characteristics agree with the Legendre-Magnaudet correlation. The use of a single dimensionless group such as Re, the Eötvös number, the Weber number and the Capillary number cannot correlate C_L of non-spherical bubbles. The trend of the critical Re for the reversal of the sign of C_L is the same as that for the onset of oscillation of bubble motion, which supports the mechanism proposed by Adoua et al., at least within the range of −6.6 ≤ logM ≤ −3.2. An experimental database of C_L is provided for validation of available C_L models and CFD.     

A computational study of the effect of initial bubble conditions on the motion of a gas bubble rising in viscous liquids

In order to examine the influence of initial bubble conditions on bubble rise motion, two-dimensional direct numerical simulations of the motion of a gas bubble rising in viscous liquids were carried out by a coupled level set/volume-of-fluid (CLSVOF) method. For dimensionless groups predicting a “spherical-cap bubble shape” (high Eötvös and low Morton numbers), we have found computationally that solutions depend on initial bubble conditions. Specifically, for spherical-cap bubble areas, we could obtain computational results of toroidal bubbles or spherical-cap bubbles depending on initial bubble conditions. On the other hand, we showed for low Eo and high M conditions that initial bubble conditions did not affect the final state of bubble rise motion.     

Terminal velocity of a Taylor drop in a vertical pipe

A scaling analysis based on the field equations for two phases and the jump conditions at the interface is carried out to deduce a balance of forces acting on a Taylor drop rising through stagnant liquid in a vertical pipe. The force balance is utilized to deduce a functional form of an empirical correlation of terminal velocity of the Taylor drop. Undetermined coefficients in the correlation are evaluated by making use of available correlations for two limiting cases, i.e. extremely high and low Reynolds number Taylor bubbles in large pipes. Terminal velocity data obtained by interface tracking simulations are also used to determine the coefficients. The proposed correlation expresses the Froude number Fr as a function of the drop Reynolds number Re_D, the Eötvös number Eo_D and the viscosity ratio μ^*. Comparisons between the correlation, simulations and experimental data confirm that the proposed correlation is applicable to Taylor drops under various conditions, i.e., 0.002 < Re_D < 4960, 4.8 < Eo_D < 228, 0 ? μ^* ? 70, 1 < N < 14700, −12 < log M < 4, and d/D < 1.6, where N is the inverse viscosity number, M the Morton number, d the sphere-volume equivalent drop diameter and D the pipe diameter.     

Computational approach for a pair of bubble coalescence process

The coalescence of bubbles has great value in mineral recovery and oil industry. In this paper, two co-axial bubbles rising in a cylinder is modelled to study the coalescence of bubbles for four computational experimental test cases. The Reynolds’ (Re) number is chosen in between 8.50 and 10, Bond number, Bo ∼4.25-50, Morton number, M 0.0125-14.7. The viscosity ratio (μ_r) and density ratio (ρ_r) of liquid to bubble are kept constant (100 and 850 respectively). It was found that the Bo number has significant effect on the coalescence process for constant Re, μ_r and ρ_r. The bubble-bubble distance over time was validated against published experimental data. The results show that VOF approach can be used to model these phenomena accurately. The surface tension was changed to alter the Bo and density of the fluids to alter the Re and M, keeping the μ_r and ρ_r the same. It was found that for lower Bo, the bubble coalesce is slower and the pocket at the lower part of the leading bubble is less concave (towards downward) which is supported by the experimental data.     

Motion of singles bubbles moving under a slightly inclined surface through stationary liquids

The influence of the liquid properties on the dynamical bubble shape and on the bubble motion has been investigated for bubbles moving under a downward facing inclined surface. The Morton number Mo varied from 2.59 × 10⁻¹¹ to 2.52 × 10⁺⁰¹. The Bond number Bo covered the range from 10 to 150 and the surface inclination angle θ was varied from 2° to 6°. To cover the wide range of Mo, several liquids such as glycerine, propanediol, water and isopropanol were used. The results have shown that the relation Fr = Fr(Bo, Mo, θ) is not adequate to describe the bubble motion, where Fr is the terminal Froude number. The choice of the terminal Reynolds number Re as the dependent parameter, allowed the clarification of the role of the Morton number on the bubble motion. At a given Bond number, the bubble Reynolds number decreases monotonously with the Morton number. Furthermore, an empirical correlation Re = Re(Bo, Mo, θ) is given that can be readily used in the mathematical modelling of bubble laden flows under solids.     

Ellipsoidal model of the rise of a Taylor bubble in a round tube

The rise velocity of long gas bubbles (Taylor bubbles) in round tubes is modeled by an ovary ellipsoidal cap bubble rising in an irrotational flow of a viscous liquid. The analysis leads to an expression for the rise velocity which depends on the aspect ratio of the model ellipsoid and the Reynolds and Eötvös numbers. The aspect ratio of the best ellipsoid is selected to give the same rise velocity as the Taylor bubble at given values of the Eötvös and Reynolds numbers. The analysis leads to a prediction of the shape of the ovary ellipsoid which rises with same velocity as the Taylor bubble.     

Bubble interaction in low-viscosity liquids

An experimental study investigated how freely rising ellipsoidal bubbles approach each other, make contact and coalesce or breakup. Pulsed planar swarms of 10–20 bubbles with Eötvös numbers from 6.0 to 27.5 were released simultaneously in aqueous solutions of 0–48 wt% sugar with Morton numbers from 3.2 × 10⁻¹¹ to 3.7 × 10⁻⁶. Bubble interaction was recorded by a video camera following the rising bubbles. Essentially, all coalescence and breakup events occurred after, not during, wake-induced collisions by a complex process related to the bubble vortex shedding cycle. This same process was also found in multi-bubble clusters and may account for excess turbulent kinetic energy generation in bubbly flow.     

Effect of incident shock wave strength on the decay of Richtmyer–Meshkov instability-introduced perturbations in the refracted shock wave

The effect of incident shock wave strength on the decay of interface introduced perturbations in the refracted shock wave was studied by performing 20 different simulations with varying incident shock wave Mach numbers (M ~ 1.1? 3.5). The analysis showed that the amplitude decay can be represented as a power law model shown in Eq.7, where A is the average amplitude of perturbations (cm), B is the base constant (cm^?(E?1), S is the distance travelled by the refracted shockwave (cm), and E is the power constant. The proposed model fits the data well for low incident Mach numbers, while at higher mach numbers the presence of large and irregular late time oscillations of the perturbation amplitude makes it hard for the power law to fit as effectively. When the coefficients from the power law decay model are plotted versus Mach number, a distinct transition region can be seen. This region is likely to result from the transition of the post-shock heavy gas velocity from subsonic to supersonic range in the lab frame. This region separates the data into a high and low Mach number region. Correlations for the power law coefficients to the incident shock Mach number are reported for the high and low Mach number regions. It is shown that perturbations in the refracted shock wave persist even at late times for high incident Mach numbers.     

The effect of the VOF–CSF static contact angle boundary condition on the dynamics of sliding and bouncing ellipsoidal bubbles

The static contact angle is the only empiricism introduced in a Volume of Fluid–Continuum Surface Force (VOF–CSF) model of bubbly flow. Although it has previously been shown to have a relatively limited effect on the accuracy of velocity and shape predictions in the case of large gas bubbles sliding under inclined walls (e.g. Cook and Behnia, 2001), it may have a more determining influence on the numerical prediction of the dynamics of smaller ellipsoidal bubbles which were shown by Tsao and Koch (1997) to bounce repeatedly when sliding under inclined walls at certain wall inclinations. The present paper reports on the influence of surface tension and the static contact angle on the dynamics of an ellipsoidal air bubble of equivalent diameter D_e = 3.4 mm. The bubble Eötvös and Morton numbers are Eo = 1.56 and Mo = 2 × 10⁻¹¹ respectively. The computational results are achieved with a Piecewise Linear Construction (PLIC) of the interface and are reviewed with reference to experimental measurements of bubble velocity and interface shape oscillations recorded using a high speed digital camera. Tests are performed at plate inclination angles θ ∈ {10°, 20°, 30°, 45°} to the horizontal and computational models consider three static contact angles θ_c ∈ {10°, 20°, 30°}. The static contact angle has been found to have a significant effect on the bubble dynamics but to varying degree depending on the plate inclination. It is shown to promote lift off and bouncing when the plate inclination angle reaches 30° in a way that does not necessarily reflect experimental observations.     

Hyperbolicity and change of type in the flow of viscoelastic fluids through channels

We consider steady flow of an upper convected Maxwell fluid through a channel with wavy walls. The vorticity of this flow will change type when the velocity in the center of the channel is larger than a critical value defined by the propagation of shear waves. There is then a central region of the channel in which the vorticity equation is hyperbolic and a low speed region near the walls where the vorticity equation is elliptic. We linearize the problem for small amplitude waviness and the linearized problem is solved in detail. The characterstic nets depend on the viscoelastic “Mach” number which is the ratio (M = U/c) of the unperturbed maximum velocity U to the speed of shear waves c into the fluid at rest and the elasticity number E. There is a supercritical (hyperbolic) region around the center of the channel when M > 1. When M ? 1, the thickness of this hyperbolic region is small when E is large, and large when E is small. Regions of positive and negative vorticity are swept out along forward facing characteristics in the hyperbolic region. There is rapid damping of vorticity in the hyperbolic region away from the boundary when M ? 1 and the Weissenberg number

. (The Weissenberg number is proportional to the relaxation time of the fluid.)The rate of damping of vorticity decreases as W is increased. Flows with high M appear to be more “elastic” when W is large in the sense that the damping is suppressed as the relaxation time of the fluid is increased.     

The non-uniform transformation strain problem for an anisotropic ellipsoidal inclusion

The problem of an anisotropic ellipsoidal inclusion which undergoes a stress-free transformation strain (in the sense of J.D. Eshelby) is considered, and the following theorem is proved: If an ellipsoidal region in an infinite anisotropic linear elastic medium undergoes, in the absence of its surroundings, a stress-free transformation strain which is a polynomial of degree M in the position coordinates x_t, then the final stress and strain state in the transformed inclusion, when constrained by its surroundings, is also a polynomial of degree M in x_t.     

Terminal velocities of clean and fully-contaminated drops in vertical pipes

Terminal velocities and shapes of drops rising through vertical pipes in clean and fully-contaminated systems are measured by using a high-speed video camera and an image processing method. Silicon oils and glycerol water solutions are used for the dispersed and continuous phases, respectively. Triton X-100 is used for surfactant. Clean and contaminated drops take either spherical, spheroidal or deformed spheroidal shapes when the diameter ratio λ is less than a critical value, λ_C, whereas they take bullet shapes for λ > λ_C (Taylor drops). The applicability of available drag and Froude number correlations is examined through comparisons with the measured data. Effects of surfactant on the shape and terminal velocity of a Taylor drop are also discussed based on the experimental data and interface tracking simulations. The conclusions obtained are as follows: (1) drag and Froude number correlations proposed so far give reasonable estimations of the terminal velocities of clean drops at any λ, (2) the terminal velocities of contaminated drops are well evaluated by making the viscosity ratio μ^* infinity in the drag correlation for clean drops in the viscous force dominant regime, (3) the effects of surfactant on the shape and terminal velocity of a Taylor drop become significant as the Eötvös number, Eo_D, decreases and μ^* increases, and (4) the reduction in surface tension due to the addition of surfactant would be the cause of the increase in the terminal velocity and elongation of a contaminated Taylor drop.     

A note on magnetoconvection in a vertical enclosure

The combined effect of viscous and ohmic dissipations on magnetoconvection in a vertical enclosure heated at the vertical side walls in the presence of applied electric field parallel to gravity and magnetic field normal to gravity is investigated. The coupled non-linear equations governing the motion are solved both analytically valid for small buoyancy parameter N and numerically valid for large N. Solutions for large N reveal a marked change in velocity profile, mass flow rate, skin friction and rate of heat transfer. These results are presented for various Hartmann number M, electric field loading parameter E and buoyancy parameter N. It is shown in the case of open circuit (i.e. E ≠ 0) that the effect of magnetic field is to increase both the velocity and temperature in contrast with the short circuit case (i.e. E = 0). The results for the case when the walls are maintained at the same temperatures (i.e. T₁ = T₂) are obtained as a particular case.     

Experimental study of single bubble motion in a liquid metal column exposed to a DC magnetic field

The motion of single Argon bubbles rising in the eutectic alloy GaInSn under the influence of a DC longitudinal magnetic field (parallel to the direction of bubble motion) was examined. The magnetic field strength was varied up to 0.3 T corresponding to a magnetic interaction parameter N (which measures the ratio of electromagnetic forces to inertial forces) slightly greater than 1. The liquid metal was at rest in a cylindrical container. Bubble and liquid velocities were measured using ultrasound Doppler velocimetry (UDV). The measured bubble terminal velocity showed oscillations indicating a zigzag movement of ellipsoidal bubbles. For small bubbles (d_e ⩽ 4.6 mm) an increase of the drag coefficient with increasing magnetic interaction parameter N was observed, whereas for larger bubbles (d_e ⩾ 5.4 mm) the application of the magnetic field reduces the drag coefficient. The measurements revealed a distinct electromagnetic damping of the bubble induced liquid velocity leading to more rectilinear bubble trajectories when the magnetic field is applied. Moreover, significant modifications of the bubble wake structure were observed. Raising of the magnetic field strength caused an enlargement of the eddies in the wake. The Strouhal number decreases with increasing magnetic interaction parameter N.     

Effect of incomplete thermal accommodation on intensive subsonic gas condensation

The direct Monte Carlo simulation method is used for investigating the effect of the thermal accommodation coefficient α _E on the relation for the Knudsen layer in the presence of intensive subsonic condensation. It is shown that the deviation of α _E from unity may significantly affect the flow parameters, in particular, M_lim, the Mach number value limiting for subsonic condensation. It is shown that a decrease in α _E leads to an increase in M_lim (M_lim < 1) if the relative flow temperature (ratio of the outer Knudsen layer boundary to the surface temperature) T < 1 and to a decrease in M_lim if T > 1. It is shown that for mirror reflection of molecules from the surface this effect may intensify.     

Image decomposition method for the analysis of a mixed dislocation in a general multilayer

The elastic solutions for a mixed dislocation in a general multilayer with N dissimilar anisotropic layers are obtained via a generalized image decomposition method. The original problem is decomposed into N homogeneous subproblems with strategically placed continuously distributed image (virtual) dislocations which satisfy the consistency conditions for degenerate N − M (M < N) layer problems. The image dislocations are used to satisfy the interface or free surface conditions, and represent the unknowns of the problem. The resulting singular Cauchy integral equations are transformed into non-singular Fredholm integral equations of the second kind using certain H- and I-integral transforms. The Fredholm integral equations are then solved via the classical Nyström method. The general decomposition and the elimination of all singular integrals yield an exact formulation of the problem; the approximation arises only in the Nyström method. The dislocation mixity and the number of layers dissimilar in thickness and elastic anisotropy can be handled without difficulty, constrained only by the number of linear algebraic equations in the Nyström method for large N. For the numerical study, image forces on a dislocation in two- and three-layer systems are calculated. The accuracy of the results is verified by checking the boundary conditions and by comparison with previous results. The dependence of the image force on the dislocation position and mixity, and on the layer thicknesses and elastic anisotropies, is also illustrated via numerical investigations.     

Comment on “Correlation of entrainment for annular flow in horizontal pipes”, by Pan,L., Hanratty,T.J., Int. J. Multiphase flow, 28(3), (2002), pp. 385–408

The above referenced paper, published in International Journal of Multiphase Flow (Pan and Hanratty, 2002), proposed an entrainment fraction correlation for annular flow in horizontal pipes. The entrainment fraction in annular flow is defined as the ratio of the mass flow rate of the liquid droplets in the gas to the total mass flow of liquid, F_E = W_LE/W_L. The proposed correlation was verified with experimental data for liquids with viscosities close to that of water. The proposed entrainment fraction correlation includes another correlation for the critical film flow rate, W_F,cr to estimate a maximum entrainment fraction F_E,max. It is shown that the critical film flow rate correlation can result in negative maximum entrainment fraction values, for low liquid flow rates.     

Efficient methods for calculating the number of states,levels and lines in atomic configurations

Exact or statistical methods for determining the distribution of the M_J values (projection of total angular momentum J) in an electron configuration are presented. This distribution, noted P(M_J), is used to calculate the allowed values of J and the number of electric-dipolar (E1) lines between two configurations. First, the difficulty to account for the Pauli exclusion principle for equivalent electrons is stressed. Showing the limit of the usual exact approach, a very efficient recursive technique is proposed for determining exactly the distribution P(M_J). Second, the statistical approach of Bauche and Bauche-Arnoult [J. Phys. B Atom. Mol. Opt. Phys. 20 (1987) 1659] is extended in order to account for configurations with a high-? spectator. In this case, identical consecutive values may exist in the center of P(M_J), which can neither be modeled by a Gaussian nor by a Gram–Charlier type function. It is shown that the Generalized Gaussian function, with the exponent constrained by the kurtosis (reduced fourth-order centered moment) of P(M_J), is more suited in these situations. A new analytical formula for the evaluation of the number of E1 lines with a larger range of applicability is then proposed.     

Characteristic lengths and maximum entropy estimation from probe signals in the ellipsoidal bubble regime

The bubble size, surface and volume distributions in two and three phase flows are essential to determine energy and mass transfer processes. The traditional approaches commonly use a conditional probability density function of chord-lengths to calculate the bubble size distribution, when the bubble size, shape and velocity are known. However, the approach used in this paper obtains the above distributions from statistical relations, requiring only the moments inferred from the measurements given by a sampling probe. Using image analysis of bubbles injected in a water tank, and placing an ideal probe on the image, a sample of bubble diameter, shape factor and velocity angle are obtained. The samples of the bubble chord-length are synthetically generated from these variables. Thus, we propose a semi-parametric approach based on the maximum entropy (MaxEnt) distribution estimation subjected to a number of moment constraints avoiding the use of the complex backward transformation. Therefore, the method allows us to obtain the distributions in close form. The probability density functions of the most important length scales (D, D₂₀, D₃₀, D₃₂), obtained applying the semi-parametric approach proposed here in the ellipsoidal bubble regime, are compared with experimental measurements.