Analysis of the Motion and Deformation of a Bubble Rising in a Viscoelastic Fluid

2006 ◽  
Author(s):  
Shriram Pillapakkam ◽  
Pushpendra Singh ◽  
Denis L. Blackmore ◽  
Nadine Aubry
Keyword(s):  
Velocity Field ◽  
Newtonian Fluid ◽  
Viscoelastic Fluid ◽  
Bubble Radius ◽  
Polymer Concentration ◽  
Deborah Number ◽  
Two Phase ◽  
Bubble Volume ◽  
Recirculation Zones ◽  
Bubble Rising

A finite element code based on the level set method is developed for performing two and three dimensional direct numerical simulations (DNS) of viscoelastic two-phase flow problems. The Oldroyd-B constitutive equation is used to model the viscoelastic liquid and both transient and steady state shapes of bubbles in viscoelastic buoyancy driven flows are studied. The influence of the governing dimensionless parameters, namely the Capillary number (Ca), the Deborah Number (De) and the polymer concentration parameter c, on the deformation of the bubble is also analyzed. Our simulations demonstrate that the rise velocity oscillates before reaching a steady value. The shape of the bubble, the magnitude of velocity overshoot and the amount of damping depend mainly on the parameter c and the bubble radius. Simulations also show that there is a critical bubble volume at which there is a sharp increase in the bubble terminal velocity as the increasing bubble volume increases, similar to the behavior observed in experiments. The structure of the wake of a bubble rising in a Newtonian fluid is strikingly different from that of a bubble rising in a viscoelastic fluid. In addition to the two recirculation zones at the equator of the bubble rising in a Newtonian fluid, two more recirculation zones exist in the wake of a bubble rising in viscoelastic fluids which influence the shape of a rising bubble. Interestingly, the direction of motion of the fluid a short distance below the trailing edge of a bubble rising in a viscoelastic fluid is in the opposite direction to the direction of the motion of the bubble, thus creating a “negative wake”. In this paper, the velocity field in the wake of the bubble, the effect of the parameters on the velocity field and their influence on the shape of the bubble are also investigated.

Transient and steady state of a rising bubble in a viscoelastic fluid

2007 ◽  
Vol 589 ◽  
pp. 215-252 ◽  
Author(s):  
SHRIRAM B. PILLAPAKKAM ◽  
PUSHPENDRA SINGH ◽  
DENIS BLACKMORE ◽  
NADINE AUBRY
Keyword(s):  
Steady State ◽  
Vortex Ring ◽  
Newtonian Fluid ◽  
Fluid Velocity ◽  
Polymer Concentration ◽  
Critical Volume ◽  
Wake Structure ◽  
Bubble Volume ◽  
Rising Bubble ◽  
Transient And Steady State

A finite element code based on the level-set method is used to perform direct numerical simulations (DNS) of the transient and steady-state motion of bubbles rising in a viscoelastic liquid modelled by the Oldroyd-B constitutive equation. The role of the governing dimensionless parameters, the capillary number (Ca), the Deborah number (De) and the polymer concentration parameter c, in both the rising speed and the deformation of the bubbles is studied. Simulations show that there exists a critical bubble volume at which there is a sharp increase in the terminal velocity with increasing bubble volume, similar to the behaviour observed in experiments, and that the shape of both the bubble and its wake structure changes fundamentally at that critical volume value. The bubbles with volumes smaller than the critical volume are prolate shaped while those with volumes larger than the critical volume have cusp-like trailing ends. In the latter situation, we show that there is a net force in the upward direction because the surface tension no longer integrates to zero. In addition, the structure of the wake of a bubble with a volume smaller than the critical volume is similar to that of a bubble rising in a Newtonian fluid, whereas the wake structure of a bubble with a volume larger than the critical value is strikingly different. Specifically, in addition to the vortex ring located at the equator of the bubble similar to the one present for a Newtonian fluid, a vortex ring is also present in the wake of a larger bubble, with a circulation of opposite sign, thus corresponding to the formation of a negative wake. This not only coincides with the appearance of a cusp-like trailing end of the rising bubble but also propels the bubble, the direction of the fluid velocity behind the bubble being in the opposite direction to that of the bubble. These DNS results are in agreement with experiments.

Negative Wake and Velocity of a Bubble Rising in a Viscoelastic Fluid

2003 ◽  
Author(s):  
Shriram B. Pillapakkam ◽  
Pushpendra Singh
Keyword(s):  
Viscoelastic Fluid ◽  
Fluid Velocity ◽  
Polymer Concentration ◽  
Three Dimensional ◽  
Small Range ◽  
Rise Velocity ◽  
Continuous Change ◽  
Dimensional Finite Element ◽  
Bubble Rising ◽  
Three Dimensional Finite Element

A three-dimensional finite element based numerical method is used to simulate the rise of a bubble in a viscoelastic fluid modeled by the Oldroyd-B model. The rise velocity is studied as a function of the bubble volume on a log-log plot. The dependence of rise velocity on the bubble shape and the viscoelastic properties of the ambient fluid are also investigated. In simulations, rather than a jump in the rise velocity at a critical volume as observed in experiments, we find that there is a steep, but continuous, change in the rise velocity over a very small range of bubble volumes. Interestingly, this steep increase in the rise velocity is exaggerated when a parameter, which is a measure of the polymer concentration, is increased, while keeping the zeroshear viscosity fixed. The wake is ‘negative’ in the sense that the direction of fluid velocity behind the bubble for this parameter range is the opposite of that for a Newtonian fluid.

Simulation of mechanically stirred two-phase liquid-gas flow

1986 ◽  
Vol 51 (5) ◽  
pp. 1001-1015 ◽  
Author(s):  
Ivan Fořt ◽  
Vladimír Rogalewicz ◽  
Miroslav Richter
Keyword(s):  
Spatial Distribution ◽  
Velocity Field ◽  
Gas Flow ◽  
Liquid Velocity ◽  
Model Parameters ◽  
Two Phase ◽  
Mean Velocity ◽  
Rushton Turbine ◽  
Low Viscosity ◽  
Volumetric Gas Flow Rate

The study describes simulation of the motion of bubbles in gas, dispersed by a mechanical impeller in a turbulent low-viscosity liquid flow. The model employs the Monte Carlo method and it is based both on the knowledge of the mean velocity field of mixed liquid (mean motion) and of the spatial distribution of turbulence intensity ( fluctuating motion) in the investigated system - a cylindrical tank with radial baffles at the wall and with a standard (Rushton) turbine impeller in the vessel axis. Motion of the liquid is then superimposed with that of the bubbles in a still environment (ascending motion). The computation of the simulation includes determination of the spatial distribution of the gas holds-up (volumetric concentrations) in the agitated charge as well as of the total gas hold-up system depending on the impeller size and its frequency of revolutions, on the volumetric gas flow rate and the physical properties of gas and liquid. As model parameters, both liquid velocity field and normal gas bubbles distribution characteristics are considered, assuming that the bubbles in the system do not coalesce.

An implicit scheme for simulation of free surface non-Newtonian fluid flows on dynamically adapted grids

2021 ◽  
Vol 36 (3) ◽  
pp. 165-176
Author(s):  
Kirill Nikitin ◽  
Yuri Vassilevski ◽  
Ruslan Yanbarisov
Keyword(s):  
Free Surface ◽  
Numerical Model ◽  
Newtonian Fluid ◽  
Viscoelastic Fluid ◽  
Fluid Flows ◽  
Numerical Results ◽  
Implicit Scheme ◽  
New Approach

Abstract This work presents a new approach to modelling of free surface non-Newtonian (viscoplastic or viscoelastic) fluid flows on dynamically adapted octree grids. The numerical model is based on the implicit formulation and the staggered location of governing variables. We verify our model by comparing simulations with experimental and numerical results known from the literature.

scholarly journals Bubble Dynamics in a Two-Phase Bubbly Mixture

2010 ◽  
Author(s):  
Arvind Jayaprakash ◽  
Sowmitra Singh ◽  
Georges Chahine
Keyword(s):  
Void Fraction ◽  
High Speed ◽  
Bubble Dynamics ◽  
Bubble Radius ◽  
Discharge Voltage ◽  
Large Bubble ◽  
Water Mixture ◽  
Two Phase ◽  
Mixture Density ◽  
Bubbly Medium

The dynamics of a primary relatively large bubble in a water mixture including very fine bubbles is investigated experimentally and the results are provided to several parallel on-going analytical and numerical approaches. The main/primary bubble is produced by an underwater spark discharge from two concentric electrodes placed in the bubbly medium, which is generated using electrolysis. A grid of thin perpendicular wires is used to generate bubble distributions of varying intensities. The size of the main bubble is controlled by the discharge voltage, the capacitors size, and the pressure imposed in the container. The size and concentration of the fine bubbles can be controlled by the electrolysis voltage, the length, diameter, and type of the wires, and also by the pressure imposed in the container. This enables parametric study of the factors controlling the dynamics of the primary bubble and development of relationships between the bubble characteristic quantities such as maximum bubble radius and bubble period and the characteristics of the surrounding two-phase medium: micro bubble sizes and void fraction. The dynamics of the main bubble and the mixture is observed using high speed video photography. The void fraction/density of the bubbly mixture in the fluid domain is measured as a function of time and space using image analysis of the high speed movies. The interaction between the primary bubble and the bubbly medium is analyzed using both field pressure measurements and high-speed videography. Parameters such as the primary bubble energy and the bubble mixture density (void fraction) are varied, and their effects studied. The experimental data is then compared to simple compressible equations employed for spherical bubbles including a modified Gilmore Equation. Suggestions for improvement of the modeling are then presented.

Experiment and Numerical Simulation of Bubble Behaviors in Argon Gas Injection Into Lead-Bismuth Pool

2006 ◽  
Author(s):  
Yumi Yamada ◽  
Toyou Akashi ◽  
Minoru Takahashi
Keyword(s):  
Direct Contact ◽  
Two Phase Flow ◽  
Experimental Result ◽  
Feed Water ◽  
Force Model ◽  
Phase Flow ◽  
Two Dimensional ◽  
Two Phase ◽  
Argon Gas ◽  
Bubble Rising

In a lead-bismuth alloy (45%Pb-55%Bi) cooled direct contact boiling water fast reactor (PBWFR), steam can be produced by direct contact of feed water with primary Pb-Bi coolant in the upper core plenum, and Pb-Bi coolant can be circulated by buoyancy forces of steam bubbles. As a basic study to investigate the two-phase flow characteristics in the chimneys of PBWFR, a two-dimensional two-phase flow was simulated by injecting argon gas into Pb-Bi pool in a rectangular vessel (400mm in length, 1500mm in height, 50mm in width), and bubble behaviors were investigated experimentally. Bubble sizes, bubble rising velocities and void fractions were measured using void probes. Argon gas was injected through five nozzles of 4mm in diameter into Pb-Bi at two locations. The experimental conditions are the pressure of atmospheric pressure, Pb-Bi temperatures of 443K, and the flow rate of injection Ar gas is 10, 20, and 30 NL/min. The measured bubble rising velocities were distributed in the range from 1 to 3 m/s. The average velocity was about 0.6 m/s. The measured bubble chord lengths were distributed from 1mm up to 30mm. The average chord length was about 7mm. An analysis was performed by two-dimensional and two-fluid model. The experimental results were compared with the analytical results to evaluate the validity of the analytical model. Although large diameter bubbles were observed in the experiment, the drag force model for spherical bubbles performed better for simulation of the experimental result because of high surface tension force of Pb-Bi.

A comparative study of MHD fluid-particle suspension induced by metachronal wave under the effects of lubricated walls

2021 ◽  
pp. 2150204
Author(s):  
Mubbashar Nazeer ◽  
Farooq Hussain ◽  
Laiba Shabbir ◽  
Adila Saleem ◽  
M. Ijaz Khan ◽  
...  
Keyword(s):  
Heat Transfer ◽  
Magnetic Fields ◽  
Multiphase Flow ◽  
Thermal Radiation ◽  
Newtonian Fluid ◽  
Closed Form Solution ◽  
Casson Fluid ◽  
Form Solution ◽  
Phase Flow ◽  
Two Phase

In this paper, the two-phase flow of non-Newtonian fluid is investigated. The main source of the flow is metachronal waves which are caused by the back and forth motion of cilia attached to the opposite walls of the channel. Magnetohydrodynamics (MHD) of Casson fluid experience the effects of transverse magnetic fields incorporated with the slippery walls of the channel. Thermal effects are examined by taking Roseland’s approximation and application of thermal radiation into account. The heat transfer through the multiphase flow of non-Newtonian fluid is further, compared with Newtonian bi-phase flow. Since the main objective of the current study is to analyze heat transfer through an MHD multiphase flow of Casson fluid. The two-phase heated flow of non-Newtonian fluid is driven by cilia motion results in nonlinear and coupled differential equations which are transformed and subsequently, integrated subject to slip boundary conditions. A closed-form solution is eventually obtained form that effectively describes the flow dynamics of multiphase flow. A comprehensive parametric study is carried out which highlights the significant contribution of pertinent parameters of the heat transfer of Casson multiphase flow. It is inferred that lubricated walls and magnetic fields hamper the movement of multiphase flow. It is noted that a sufficient amount of additional thermal energy moves into the system, due to the Eckert number and Prandtl number. While thermal radiation acts differently by expunging the heat transfer. Moreover, Casson multiphase flow is a more suitable source of heat transfer than Newtonian multiphase flow.

scholarly journals A MATHEMATICAL MODEL OF NONSTATIONARY MOTION OF A VISCOELASTIC FLUID IN ROLLER BEARINGS

2019 ◽  
Vol 81 (4) ◽  
pp. 501-512
Author(s):  
I.A. Zhurba Eremeeva ◽  
D. Scerrato ◽  
C. Cardillo ◽  
A. Tran
Keyword(s):  
Mathematical Model ◽  
Viscoelastic Fluid ◽  
Viscoelastic Properties ◽  
Fluid Model ◽  
Maxwell Fluid ◽  
Initial Boundary ◽  
Deborah Number ◽  
Nonstationary Motion ◽  
Friction Forces ◽  
Roller Bearings

Nowadays, the emergence of new lubricants requires an enhancement of the rheological models and methods used for solution of corresponding initial boundary-value problems. In particular, models that take into account viscoelastic properties are of great interest. In the present paper we consider the mathematical model of nonstationary motion of a viscoelastic fluid in roller bearings. We used the Maxwell fluid model for the modeling of fluid properties. The viscoelastic properties are exhibited by many lubricants that use polymer additives. In addition, viscoelastic properties can be essential at high fluid speeds. Also, viscoelastic properties can be significant in the case of thin gaps. Maxwell's model is one of the most common models of viscoelastic materials. It combines the relative simplicity of constitutive equations with the ability to describe a stress relaxation. In addition, viscoelastic fluids also allow us to describe some effects that are missing in the case of viscous fluid. An example it is worth to mention the Weissenberg effect and a number of others. In particular, such effects can be used to increase the efficiency of the film carrier in the sliding bearings. Here we introduced characteristic assumptions on the form of the flow, allowing to significantly simplify the solution of the problem. We consider so-called self-similar solutions, which allows us to get a solution in an analytical form. As a result these assumptions, the formulae for pressure and friction forces are derived. Their dependency on time and Deborah number is analyzed. The limiting values of the flow characteristics were obtained. The latter can be used for steady state of the flow regime. Differences from the case of Newtonian fluid are discussed. It is shown that viscoelastic properties are most evident at the initial stage of flow, when the effects of non-stationarity are most important.

Direct Numerical Simulation of Drops in Three Dimensional Viscoelastic Simple Shear Flows

2001 ◽  
Author(s):  
Shriram B. Pillapakkam ◽  
Pushpendra Singh
Keyword(s):  
Numerical Simulation ◽  
Shear Flow ◽  
Direct Numerical Simulation ◽  
Viscoelastic Fluid ◽  
Simple Shear ◽  
Polymer Concentration ◽  
Three Dimensional ◽  
Simple Shear Flow ◽  
Splitting Technique

Abstract A three dimensional finite element scheme for Direct Numerical Simulation (DNS) of viscoelastic two phase flows is implemented. The scheme uses the Level Set Method to track the interface and the Marchuk-Yanenko operator splitting technique to decouple the difficulties associated with the governing equations. Using this numerical scheme, the shape of Newtonian drops in a simple shear flow of viscoelastic fluid and vice versa are analyzed as a function of Capillary number, Deborah number and polymer concentration. The viscoelastic fluid is modeled via the Oldroyd-B model. The role of viscoelastic stresses in deformation of a drop subjected to simple shear flow and its effect on the steady state shape is analyzed. Our results compare favorably with existing experimental data and also help in understanding the role of viscoelastic stresses in drop deformation.

scholarly journals Insight of thermally stratified Jeffrey fluid flow inside porous medium subject to chemical species and melting heat transfer

2019 ◽  
Vol 11 (9) ◽  
pp. 168781401987618
Author(s):  
Mubashar Javed ◽  
Muhammad Farooq ◽  
Aisha Anjum ◽  
Shakeel Ahmad
Keyword(s):  
Heat Transfer ◽  
Velocity Field ◽  
Thermal Stratification ◽  
Heterogeneous Reaction ◽  
Concentration Field ◽  
Chemical Species ◽  
Internal Heat ◽  
Nonlinear Problems ◽  
Deborah Number ◽  
Jeffrey Fluid

This article concentrates on two-dimensional magnetohydrodynamic stagnation flow of Jeffrey liquid on a nonlinearly stretching sheet which possesses variable thickness. Simultaneous impact of melting as well as thermal stratification is specifically investigated in this study due to their tremendous involvement in plenty of natural and industrial processes. Internal heat generation and presence of chemical species are considered to ponder at heat transfer properties. Series solution has been obtained by solving the developed nonlinear problems. Physical behavior of various controlling parameters such as velocity, thermal, and concentration fields are investigated. It has been found that temperature field decays due to higher intensity of thermal stratification parameter, but thickness of thermal boundary layer boosts up. Larger Deborah number results in incremented velocity field. For uplifted wall thickness parameter, velocity field depreciates. Concentration field declines for enhanced parameters of homogeneous as well as heterogeneous reaction. Moreover, velocity is decreasing function of porosity parameter.

Sign in / Sign up

Export Citation Format

Share Document