Unsteady motion of a spherical bubble in a complex fluid: Mathematical modelling and simulation

The nonlinear response of an oscillatory bubble in a complex fluid is studied. The bubble is immersed in a Newtonian liquid, which may have a dilute volume fraction of anisotropic additives such as fibers or few ppm of macromolecules. The constitutive equation for the fluid is based on a Maxwell mod...

Full description

Saved in:
Bibliographic Details
Published in:Applied mathematical modelling 2013-11, Vol.37 (20-21), p.8972-8984
Main Authors: Albernaz, D.L., Cunha, F.R.
Format: Article
Language:English
Subjects:
Additives
Asymptotic properties
Bubble
Bubbles
Collapse
Complex fluid
Computational fluid dynamics
Computer simulation
Fluid flow
Fluids
Mathematical models
Modelling
Orientation
Viscoelastic
Citations: Items that this one cites
Items that cite this one
Online Access:Get full text
Tags: Add Tag
No Tags, Be the first to tag this record!
Description
Summary:The nonlinear response of an oscillatory bubble in a complex fluid is studied. The bubble is immersed in a Newtonian liquid, which may have a dilute volume fraction of anisotropic additives such as fibers or few ppm of macromolecules. The constitutive equation for the fluid is based on a Maxwell model with an extensional viscosity for the viscous contribution. The model is considered new in the study of bubble dynamics in complex fluids. The numerical computation solves a system of three first order ordinary differential equations, including the one associated with the solution of the convolution integral, using a fifth order Runge–Kutta scheme with appropriated time steps. Asymptotic solutions of governing equation are developed for small values of the pressure forcing amplitude and for small values of the elastic parameter. A study of the bubble collapse radius is also presented. We compare the results predicted by our model with other model in the literature and a good agreement is observed. The calculated asymptotic solutions are also used to test the results of the numerical simulations. In addition, the orientation of the additives is considered. The angular probability density function is assumed to be a normal distribution. The results show that the model based on the fully aligned additives with the radial direction overestimates the tendency of the additives to stabilize the bubble motion, since the effect of extensional viscosity occurs due to the particle resistance to the movement throughout its longitudinal direction.
ISSN:0307-904X
DOI:10.1016/j.apm.2013.03.065