Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Tuning a magnetic field to generate spinning ferrofluid droplets with controllable speed via nonlinear periodic interfacial waves

Zongxin Yu and Ivan C. Christov
Phys. Rev. E 103, 013103 – Published 12 January 2021
PDFHTMLExport Citation

Abstract

Two-dimensional free surface flows in Hele-Shaw configurations are a fertile ground for exploring nonlinear physics. Since Saffman and Taylor's work on linear instability of fluid-fluid interfaces, significant effort has been expended to determining the physics and forcing that set the linear growth rate. However, linear stability does not always imply nonlinear stability. We demonstrate how the combination of a radial and an azimuthal external magnetic field can manipulate the interfacial shape of a linearly unstable ferrofluid droplet in a Hele-Shaw configuration. We show that weakly nonlinear theory can be used to tune the initial unstable growth. Then, nonlinearity arrests the instability and leads to a permanent deformed droplet shape. Specifically, we show that the deformed droplet can be set into motion with a predictable rotation speed, demonstrating nonlinear traveling waves on the fluid-fluid interface. The most linearly unstable wave number and the combined strength of the applied external magnetic fields determine the traveling wave shape, which can be asymmetric.

  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Figure
  • Received 14 September 2020
  • Revised 17 December 2020
  • Accepted 21 December 2020

DOI:https://doi.org/10.1103/PhysRevE.103.013103

©2021 American Physical Society

Physics Subject Headings (PhySH)

Fluid DynamicsNonlinear DynamicsPolymers & Soft Matter

Authors & Affiliations

Zongxin Yu* and Ivan C. Christov

  • School of Mechanical Engineering, Purdue University, West Lafayette, Indiana 47907, USA

Article Text (Subscription Required)

Click to Expand

Supplemental Material (Subscription Required)

Click to Expand

References (Subscription Required)

Click to Expand
Issue

Vol. 103, Iss. 1 — January 2021

Reuse & Permissions
Access Options
CHORUS

Article Available via CHORUS

Download Accepted Manuscript
Author publication services for translation and copyediting assistance advertisement

Authorization Required


×

Images

  • Figure 1
    Figure 1

    Schematic illustration of a Hele-Shaw cell confining a ferrofluid droplet, initially circular with radius R. The azimuthal magnetic field Ha is produced by a long wire conveying an electric current I. The radial magnetic field Hr is produced by a pair of anti-Helmholtz coils with equal currents IAH in opposite directions. The combined external magnetic field H deforms the droplet, and its interface is given by h(θ,t). In comparison, the fluid exterior to the droplet (e.g., air) is assumed to have negligible viscosity and velocity.

    Reuse & Permissions
  • Figure 2
    Figure 2

    (a) The evolution of the first five harmonic modes from fully nonlinear simulation (solid) and weakly nonlinear approximation (dashed), for NBa=1.0, NBr=37, R=1, and χ=1 [same parameters for (b), (c) and (d)]. (b) The fully nonlinear evolution of the interface from a small perturbation of the flat base state into a permanent traveling wave (rotating droplet). (c) Comparison between the final shape from fully nonlinear simulation (solid) and weakly nonlinear approximation (dashed). (d) Stability diagram based on the first two harmonic modes of the final shape (marked with ) shown in (b); (resp. ×) denotes the stable (resp. unstable) initial conditions, solid (resp. dashed) curves track the stable (resp. unstable) evolution trajectories. The unstable region is shaded, and the “f” superscript represents the final harmonic mode amplitude.

    Reuse & Permissions
  • Figure 3
    Figure 3

    (a) Comparison of the propagation velocity predicted by linear theory vfL (dashed), weakly nonlinear theory vfW (empty symbols), and fully nonlinear simulation vfF (filled symbols). The circles represent results for case (i) R=1 fixed and NBa[0,103,102,101,1,3,5], the triangles represent results for case (ii) NBa=1 fixed with NBr varying according to R[0.8,0.9,1.1,1.2], and the squares represent case (iii) km[5,6,7,8,9], R=1 and NBa,NBr determined so that vfL=85.16. (b) The skewness Sk of the fully nonlinear profile. (c) The permanent wave shape (only one wavelength shown).

    Reuse & Permissions
  • Figure 4
    Figure 4

    (a) Time evolution of the wave profile asymmetry for different combinations of NBa[0,0.1,1,3,5,6,7,8,9] and NBr varying so that km=7 for R=1. Solid curves represent stable cases yielding a propagating profile; dashed curves represent unstable cases in which the profile distorts and grows without bound. (b) Permanent wave profiles that emerge and propagating in a stable manner. (c) Stable (top, with NBa=1,NBr=37) and unstable (bottom, NBa=8,NBr=41) evolution of the profile. The instants of time [at which the shapes in (c) are shown] are marked with white dots in (a), superimposed on the asymmetry profiles.

    Reuse & Permissions
  • Figure 5
    Figure 5

    The unstable evolution of the first five harmonic modes (k=4,8,12,16,20) from fully nonlinear simulation (solid) and their stable evolution from weakly nonlinear approximation (dashed) for (a) nonrotating (km=4, NBa=0) and (b) rotating (km=4, NBa=1) shapes.

    Reuse & Permissions
  • Figure 6
    Figure 6

    Flow chart of the vortex-sheet algorithm using the Crank-Nicolson method for time advancement and fixed-point iteration for resolving the implicit nonlinear terms.

    Reuse & Permissions
  • Figure 7
    Figure 7

    Grid convergence study for fundamental modes km=7 (black) and km=9 (blue) with NBa=1 and N=256 (dotted), N=512 (dot-dashed), N=1024 (dashed), and N=2048 (solid). (a) Spectral energy of harmonic modes (k, 2k, 3k, ). (b) The root-mean-square error taking the N=2048 solution ẑ as “exact.” (c) Grid convergence of the evolution of the skewness Sk(t). (d) Grid convergence of the evolution of the asymmetry As(t).

    Reuse & Permissions
×

Sign up to receive regular email alerts from Physical Review E

Log In

Cancel
×

Search


Article Lookup

Paste a citation or DOI

Enter a citation
×