The effect of applied magnetic fields on the deformation of a biocompatible hydrophobic ferrofluid drop suspended in a viscous medium is investigated numerically and compared with experimental data. A numerical formulation for the time-dependent simulation of magnetohydrodynamics of two immiscible non-conducting fluids is used with a volume-of-fluid scheme for fully deformable interfaces. Analytical formulae for ellipsoidal drops and near-spheroidal drops are reviewed and developed for code validation. At low magnetic fields, both the experimental and numerical results follow the asymptotic small deformation theory. The value of interfacial tension is deduced from an optimal fit of a numerically simulated shape with the experimentally obtained drop shape, and appears to be a constant for low applied magnetic fields. At high magnetic fields, on the other hand, experimental measurements deviate from numerical results if a constant interfacial tension is implemented. The difference can be represented as a dependence of apparent interfacial tension on the magnetic field. This idea is investigated computationally by varying the interfacial tension as a function of the applied magnetic field and by comparing the drop shapes with experimental data until a perfect match is found. This estimation method provides a consistent correlation for the variation in interfacial tension at high magnetic fields. A conclusion section provides a discussion of physical effects which may influence the microstructure and contribute to the reported observations.

a b s t r a c tTypical VOF algorithms rely on an implicit slip that scales with mesh refinement, to allow contact lines to move along no-slip boundaries. As a result, solutions of contact line phenomena vary continuously with mesh spacing; this paper presents examples of that variation. A mesh-dependent dynamic contact angle model is then presented, that is based on fundamental hydrodynamics and serves as a more appropriate boundary condition at a moving contact line. This new boundary condition eliminates the stress singularity at the contact line; the resulting problem is thus well-posed and yields solutions that converge with mesh refinement. Numerical results are presented of a solid plate withdrawing from a fluid pool, and of spontaneous droplet spread at small capillary and Reynolds numbers.Published by Elsevier Inc.

a b s t r a c tThe influence of matrix and droplet viscoelasticity on the steady deformation and orientation of a single droplet subjected to simple shear is investigated microscopically. Experimental data are obtained in the velocity-vorticity and velocity-velocity gradient plane. A constant viscosity Boger fluid is used, as well as a shear-thinning viscoelastic fluid. These materials are described by means of an Oldroyd-B, Giesekus, Ellis, or multi-mode Giesekus constitutive equation. The drop-to-matrix viscosity ratio is 1.5. The numerical simulations in 3D are performed with a volume-of-fluid algorithm and focus on capillary numbers 0.15 and 0.35. In the case of a viscoelastic matrix, viscoelastic stress fields, computed at varying Deborah numbers, show maxima slightly above the drop tip at the back and below the tip at the front. At both capillary numbers, the simulations with the Oldroyd-B constitutive equation predict the experimentally observed phenomena that matrix viscoelasticity significantly suppresses droplet deformation and promotes droplet orientation. These two effects saturate experimentally at high Deborah numbers. Experimentally, the high Deborah numbers are achieved by decreasing the droplet radius with other parameters unchanged. At the higher capillary and Deborah numbers, the use of the Giesekus model with a small amount of shear-thinning dampens the stationary state deformation slightly and increases the angle of orientation. Droplet viscoelasticity on the other hand hardly affects the steady droplet deformation and orientation, both experimentally and numerically, even at moderate to high capillary and Deborah numbers.

SUMMARYRecent work on a consistent representation of surface tension and on the accurate computation of interface curvature has extended the applicability of the volume-of-fluid, or VOF, method to surface tension-driven phenomena. We have extended these methodologies to interfaces near solid surfaces; specifically, we show how height functions can be used to enforce a contact angle boundary condition at a contact line, for the full range of contact angles. As such, this work may be viewed as following up on the methodology of