Oldroyd-B model

The Oldroyd-B model is a constitutive model used to describe the flow of viscoelastic fluids. This model can be regarded as an extension of the Upper Convected Maxwell model and is equivalent to a fluid filled with elastic bead and spring dumbbells. The model is named after its creator James G. Oldroyd.[1]

The model can be written as:

${\displaystyle \mathbf {T} +\lambda _{1}{\stackrel {\nabla }{\mathbf {T} }}=2\eta _{0}(\mathbf {D} +\lambda _{2}{\stackrel {\nabla }{\mathbf {D} }})}$

where:

• ${\displaystyle \mathbf {T} }$ is the stress tensor;
• ${\displaystyle \lambda _{1}}$ is the relaxation time;
• ${\displaystyle \lambda _{2}}$ is the retardation time = ${\displaystyle {\frac {\eta _{s}}{\eta _{0}}}\lambda _{1}}$;
• ${\displaystyle {\stackrel {\nabla }{\mathbf {T} }}}$ is the Upper convected time derivative of stress tensor:

${\displaystyle {\stackrel {\nabla }{\mathbf {T} }}={\frac {\partial }{\partial t}}\mathbf {T} +\mathbf {v} \cdot \nabla \mathbf {T} -((\nabla \mathbf {v} )^{T}\cdot \mathbf {T} +\mathbf {T} \cdot (\nabla \mathbf {v} ))}$;

• ${\displaystyle \mathbf {v} }$ is the fluid velocity;
• ${\displaystyle \eta _{0}}$ is the total viscosity composed of solvent and polymer components, ${\displaystyle \eta _{0}=\eta _{s}+\eta _{p}}$;
• ${\displaystyle \mathbf {D} }$ is the deformation rate tensor or rate of strain tensor, ${\displaystyle \mathbf {D} ={\frac {1}{2}}\left[{\boldsymbol {\nabla }}\mathbf {v} +({\boldsymbol {\nabla }}\mathbf {v} )^{T}\right]}$.

The model can also be written split into polymeric (viscoelastic) part separately from the solvent part:[2]

${\displaystyle \mathbf {T} =2\eta _{s}\mathbf {D} +\mathbf {\tau } }$.

where

${\displaystyle \mathbf {\tau } +\lambda _{1}{\stackrel {\nabla }{\mathbf {\tau } }}=2\eta _{p}\mathbf {D} }$

Whilst the model gives good approximations of viscoelastic fluids in shear flow, it has an unphysical singularity in extensional flow, where the dumbbells are infinitely stretched. This is, however, specific to idealised flow; in the case of a cross-slot geometry the extensional flow is not ideal, so the stress, although singular, remains integrable, i.e. the stress is infinite in a correspondingly infinitely small region.[3]

If the solvent viscosity is zero, the Oldroyd-B becomes the Upper Convected Maxwell model.

References

1. Oldroyd, James (Feb 1950). "On the Formulation of Rheological Equations of State". Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences. 200 (1063): 523–541. Bibcode:1950RSPSA.200..523O. doi:10.1098/rspa.1950.0035.
2. Owens, R. G.,Phillips, T. N. (2002). Computational Rheology. Imperial College Press. ISBN 978-1-86094-186-3.CS1 maint: multiple names: authors list (link)
3. Poole, Rob (Oct 2007). "Purely elastic flow asymmetries". Physical Review Letters. 99 (16): 164503. Bibcode:2007PhRvL..99p4503P. doi:10.1103/PhysRevLett.99.164503. hdl:10400.6/634.