Abstract
We present a Navier-Stokes/Oldroyd-B immersed boundary algorithm that captures the interaction of a flexible structure with a viscoelastic fluid. In particular, we study the effects of bulk viscoelasticity on freely decaying shape oscillations of an Oldroyd-B fluid droplet suspended in an Oldroyd-B matrix. Our numerical data indicate that if the fluid viscosity is low, viscoelasticity plays a modulating role in the drop shape relaxation; specifically, it increases the oscillation frequency and decreases the decay rate when the fluid relaxation time is above a critical value. In the high-viscosity limit, i.e., when a Newtonian droplet is expected to return to a spherical shape with an aperiodic decay, an increase in the relaxation time eventually results in the reappearance of the oscillations. Both these results are in line with the prediction of small deformation theory for viscoelastic droplet oscillations. The algorithm was also validated by direct comparison with linear asymptotics.