Y. Dimakopoulos

Steady bubble rise and deformation in Newtonian and viscoplastic fluids and conditions for bubble entrapment

We examine the buoyancy-driven rise of a bubble in a Newtonian or a viscoplastic fluid assuming axial symmetry and steady flow. Bubble pressure and rise velocity are determined, respectively, by requiring that its volume remains constant and its centre of mass remains fixed at the centre of the coordinate system. The continuous constitutive model suggested by Papanastasiou is used to describe the viscoplastic behaviour of the material. The flow equations are solved numerically using the mixed finite-element/Galerkin method. The nodal points of the computational mesh are determined by solving a set of elliptic differential equations to follow the often large deformations of the bubble surface. The accuracy of solutions is ascertained by mesh refinement and predictions are in very good agreement with previous experimental and theoretical results for Newtonian fluids. We determine the bubble shape and velocity and the shape of the yield surfaces for a wide range of material properties, expressed in terms of the Bingham Bn=τy*ρ*g*R b Bond B o = ρ*g*R* 2 b /γ* and Archimedes A r =ρ *2 g*R *3 b /μ * o 2 numbers, where ρ* is the density, μ * o the viscosity, γ* the surface tension and τ * y the yield stress of the material, g* the gravitational acceleration and R * b the radius of a spherical bubble of the same volume. If the fluid is viscoplastic, the material will not be deforming outside a finite region around the bubble and, under certain conditions, it will not be deforming either behind it or around its equatorial plane in contact with the bubble. As Bn increases, the yield surfaces at the bubble equatorial plane and away from the bubble merge and the bubble becomes entrapped. When Bo is small and the bubble cannot deform from the spherical shape the critical Bn is 0.143, i.e. it is a factor of 3/2 higher than the critical Bn for the entrapment of a solid sphere in a Bingham fluid, in direct correspondence with the 3/2 higher terminal velocity of a bubble over that of a sphere under the same buoyancy force in Stokes flow. As Bo increases allowing the bubble to squeeze through the material more easily, the critical Bingham number increases as well, but eventually it reaches an asymptotic value. Ar affects the critical Bn value much less.

Flow of two immiscible fluids in a periodically constricted tube: Transitions to stratified, segmented, churn, spray, or segregated flow

We study the flow of two immiscible, Newtonian fluids in a periodically constricted tube driven by a constant pressure gradient. Our volume-of-fluid algorithm is used to solve the governing equations. First, the code is validated by comparing its predictions to previously reported results for stratified and pulsing flow. Then, it is used to capture accurately all the significant topological changes that take place. Initially, the fluids have a core-annular arrangement, which is found to either remain the same or change to a different arrangement depending on the fluid properties, the pressure driving the flow, or the flow geometry. The flow-patterns that appear are the core-annular, segmented, churn, spray, and segregated flow. The predicted scalings near pinching of the core fluid concur with similarity predictions and earlier numerical results [I. Cohen et al., “Two fluid drop snap-off problem: Experiments and theory,” Phys. Rev. Lett. 83, 1147–1150 (1999)]. Flow-pattern maps are constructed in terms of the Reynolds and Weber numbers. Our result provides deeper insights into the mechanism of the pattern transitions and is in agreement with previous studies on core-annular flow [Ch. Kouris and J. Tsamopoulos, “Core-annular flow in a periodically constricted circular tube, I. Steady state, linear stability and energy analysis,” J. Fluid Mech. 432, 31–68 (2001) and Ch. Kouris et al., “Comparison of spectral and finite element methods applied to the study of interfacial instabilities of the core-annular flow in an undulating tube,” Int. J. Numer. Methods Fluids 39(1), 41–73 (2002)], segmented flow [E. Lac and J. D. Sherwood, “Motion of a drop along the centreline of a capillary in a pressure-driven flow,” J. Fluid Mech. 640, 27–54 (2009)], and churn flow [R. Y. Bai et al., “Lubricated pipelining—Stability of core annular-flow. 5. Experiments and comparison with theory,” J. Fluid Mech. 240, 97–132 (1992)].

Viscous effects on the oscillations of two equal and deformable bubbles under a step change in pressure

According to linear theory and assuming the liquids to be inviscid and the bubbles to remain spherical, bubbles set in oscillation attract or repel each other with a force that is proportional to the product of their amplitude of volume pulsations and inversely proportional to the square of their distance apart. This force is attractive, if the forcing frequency lies outside the range of eigenfrequencies for volume oscillation of the two bubbles. Here we study the nonlinear interaction of two deformable bubbles set in oscillation in water by a step change in the ambient pressure, by solving the Navier–Stokes equations numerically. As in typical experiments, the bubble radii are in the range 1–1000 μm. We find that the smaller bubbles (~5 μm) deform only slightly, especially when they are close to each other initially. Increasing the bubble size decreases the capillary force and increases bubble acceleration towards each other, leading to oblate or spherical cap or even globally deformed shapes. These deformations may develop primarily in the rear side of the bubbles because of a combination of their translation and harmonic or subharmonic resonance between the breathing mode and the surface harmonics. Bubble deformation is also promoted when they are further apart or when the disturbance amplitude decreases. The attractive force depends on the Ohnesorge number and the ambient pressure to capillary forces ratio, linearly on the radius of each bubble and inversely on the square of their separation. Additional damping either because of liquid compressibility or heat transfer in the bubble is also examined.

How viscoelastic is human blood plasma

Blood plasma has been considered a Newtonian fluid for decades. Recent experiments (Brust et al., Phys. Rev. Lett., 2013, 110) revealed that blood plasma has a pronounced viscoelastic behavior. This claim was based on purely elastic effects observed in the collapse of a thin plasma filament and the fast flow of plasma inside a contraction-expansion microchannel. However, due to the fact that plasma is a solution with very low viscosity, conventional rotational rheometers are not able to stretch the proteins effectively and thus, provide information about the viscoelastic properties of plasma. Using computational rheology and a molecular-based constitutive model, we predict accurately the rheological response of human blood plasma in strong extensional and constriction complex flows. The complete rheological characterization of plasma yields the first quantitative estimation of its viscoelastic properties in shear and extensional flows. We find that although plasma is characterized by a spectrum of ultra-short relaxation times (on the order of 10-3-10-5 s), its elastic nature dominates in flows that feature high shear and extensional rates, such as blood flow in microvessels. We show that plasma exhibits intense strain hardening when exposed to extensional deformations due to the stretch of the proteins in its bulk. In addition, using simple theoretical considerations we propose fibrinogen as the main candidate that attributes elasticity to plasma. These findings confirm that human blood plasma features bulk viscoelasticity and indicate that this non-Newtonian response should be seriously taken into consideration when examining whole blood flow.