Complete modes of star-shaped oscillating drops

TL;DR

This paper provides a comprehensive analysis of star-shaped oscillations in water drops, combining surface and azimuthal modes to improve frequency prediction accuracy.

Contribution

It introduces a complete model coupling surface and azimuthal modes, advancing understanding of star-shaped drop oscillations beyond previous quasi-2D approaches.

Findings

Coupling of surface and azimuthal modes explains star-shaped oscillations.

New dispersion relation improves frequency prediction accuracy.

Surface motion patterns are driven by parametric instability.

Abstract

The star-shaped oscillation of water drops has been observed in various physical situations with different sources of vertical excitation. Previous studies apply a quasi-2D model to analyze the resonance frequency and only consider the azimuthal oscillation mode. In this paper, we find that the upper surface of water drops also develop different motion patterns due to parametric instability and it is the coupling of surface motion and azimuthal oscillation that leads to star-shaped oscillations. We will introduce the analysis of the surface mode, combining with the azimuthal mode to give a complete description of the motion of water drops. We propose a new dispersion relation based on the complete description, which provides a significant increase of accuracy in predicting the oscillating frequencies.

Figures

Equations

$${\omega }_n$$

$$\frac{\partial \varphi }{\partial t}+\frac{1}{2}(\nabla \varphi {)}^2+\frac{p}{\rho }+U=0$$

$$p=\sigma ({\gamma }_1+{\gamma }_2)$$

$${\left(\frac{\partial \varphi }{\partial t}\right)}_{z=H+h}+\frac{\sigma }{\rho }{\nabla }^{2}h+\left(g-Fcos(\Omega t)\right)h=0$$

$${\left(\frac{\partial \varphi }{\partial t}\right)}_{r=R+a}-\frac{\sigma }{\rho R^2}\left(a+\frac{{\partial }^{2}a}{\partial {\theta }^2}\right)=0$$

$${\left(\frac{\partial \varphi }{\partial z}\right)}_{z=0}=0$$

$${\left(\frac{\partial \varphi }{\partial z}\right)}_{z=H+h}$$

$${\left(\frac{\partial \varphi }{\partial r}\right)}_{r=R+a}$$

$$\frac{{\partial }^2\varphi }{\partial r^2}+\frac{1}{r}\frac{\partial \varphi }{\partial r}+\frac{1}{r^2}\frac{{\partial }^2\varphi }{\partial r^2}+\frac{{\partial }^2\varphi }{\partial z^2}=0$$

$$\varphi$$

$$h=\sum _{k,n}{J}_n(kr)sin(n\theta )B_{k,n}(t)$$

$$a=\sum _{k,n}cosh(kz)sin(n\theta )C_{k,n}(t)$$

$$\frac{d{C}_{k,n}(t)}{dt}$$

$$\frac{d{B}_{k,n}(t)}{dt}$$

$$\begin{array}{rl}\frac{d^2{B}_{k,n}}{d{t}^2}+ktanh& (kH)\left(g+\frac{\sigma k^2}{\rho }-Fcos\Omega t\right)B_{k,n}(t)=0\end{array}$$

$$p$$

$$q$$

$$\frac{d^2{B}_{k,n}}{d{t}^2}+\left(p+qcos\Omega t\right)B_{k,n}(t)=0$$

$$\frac{{\Omega }^2}{4}-\frac{q}{2}\text{textless}{\omega }^2\text{textless}\frac{{\Omega }^2}{4}+\frac{q}{2}$$

$$\frac{d^2{C}_{k,n}}{d{t}^2}+\frac{(n^2-1)k\sigma J_n^{\prime }(kR)}{\rho R^2{J}_n(kR)}{C}_{k,n}(t)=0$$

$$\begin{array}{r}\frac{(n^2-1)k\sigma J_n^{\prime }(kR)}{\rho R^2{J}_n(kR)}=ktanh(kH)\left(g+\frac{\sigma k^2}{\rho }\right)\end{array}$$

$${\omega }_{m,n}$$

$$\underset{k\to 0}{\lim }\frac{k{J}_n^{\prime }(kR)}{J_n(kR)}=\frac{n}{R}$$

$$S_{n,m}=\frac{{\omega }_{n,m}^2}{{\omega }_n^2}=\frac{R}{n}\frac{k_m{J}_n^{\prime }(k_{m}R)}{J_n(k_{m}R)}$$

Full text

††preprint: APS/123-QED

I Introduction

Liquid drops driven by vertical excitation can develop star-shaped oscillations under certain conditions, including drops on a vertically vibrating hydrophobic plate PRL-Triplon-Modes-of-Puddles , drops on a pulsating air cushion PRL-Gas-Film-Levitated-Liquids-Shape-Fluctuations-of-Viscous-Drops , acoustically levitating drops PRE-Parametrically-excited-sectorial-oscillation-of-liquid-drops-floating-in-ultrasound , metal drops placed on an oscillating magnetic field JFM-Free-surface-horizontal-waves-generated-by-low-frequency-alternating-magnetic-fields . These oscillation systems are considered to undergo parametric resonance because the water drops are observed to oscillate in half of the frequency of the excitation source experimentally JPSJ-Self-Induced-Vibration-of-a-Water-Drop-Placed-on-an-Oscillating-Plate ; EPJST-Star-drops-formed-by-periodic-excitation-and-on-an-air-cushion-A-short-review ; PRE-Parametrically-excited-sectorial-oscillation-of-liquid-drops-floating-in-ultrasound . Even in the Leidenfrost effect Annual-Review-of-Fluid-Mechanics-Leidenfrost-Dynamics ; JPSJ-Vibration-of-a-Flattened-Drop.-I.-Observation ; JPSJ-Vibration-of-a-Flattened-Drop.-II.-Normal-Mode-Analysis ; JASA-Vibrations-of-Evaporating-Liquid-Drops ; Physica-A-Nitrogen-stars-morphogenesis-of-a-liquid-drop that the water drops are levitated by water vapor on a very hot plate, recent studies have found the periodically vibrating pressure of the vapor layer under the water drops can serve as parametric excitation PRF-Star-shaped-oscillations-of-Leidenfrost-drops ; Physics-of-Fluids-The-many-faces-of-a-Leidenfrost-drop . Despite these experimental evidences, the mechanism that induces the parametric instability is still unclear.

In the past, the oscillation modes of the flattened water drops are described by Rayleigh equation, proposed for inviscid cylindrical drops. The eigen frequencies for small deformations are given by Rayleigh-On-the-Capillary-Phenomena-of-Jets ; JPSJ-Vibration-of-a-Flattened-Drop.-II.-Normal-Mode-Analysis :

[TABLE]

where $n$ and $\omega_{n}$ represent azimuthal mode and corresponding eigen frequency respectively, and $R$ is the balanced radius of the water drop, $\sigma$ the surface tension constant, $\rho$ the density of water. In the equation, the motion of the water drops is simplified as two-dimensional. Though the dispersion relation 1 gives right trend between the oscillation mode number $n$ and the frequency, more detailed studies by Bouwhuis et al. PRE-Oscillating-and-star-shaped-drops-levitated-by-an-airflow found that the measured resonance frequencies of star-shaped drops levitating on a steady ascending airflow are lower than that given by 1, overestimating of the stiffness coefficient of the oscillating drop.

The two-dimensional model apparently oversimplifies the oscillation modes and loses accuracy when upper surface motion is appreciable. Shen et al. have already found that the upper surface develops petal-shaped motion patterns PRE-Parametrically-excited-sectorial-oscillation-of-liquid-drops-floating-in-ultrasound ; JFM-Free-surface-horizontal-waves-generated-by-low-frequency-alternating-magnetic-fields ; Observation-of-the-shape-of-a-water-drop-on-an-oscillating-Teflon-plate . In our experiments, we also find that, with the increase of mode number $n$, the upper surface of the water drop becomes unstable and forms various patterns.

In this paper, we will present a complete theoretical model to derive the surface normal modes of the drop and combine the azimuthal oscillation modes to give a complete description of the motion of water drops. We verify that the vertical driving force induces an instability on the upper surface due to parametric resonance. The coupling of the surface motion and azimuthal oscillation leads to the appearance of star-shaped drops oscillating at half of the driving frequency. The surface oscillation will be described in much the way as the Faraday waves. Faraday waves are patterns of standing waves observed at the free surface, when a layer of liquid in a container is driven to vibrate vertically Faraday-Waves ; Rayleigh-On-maintained-vibrations ; The-Royal-Society-The-stability-of-the-plane-free-surface-of-a-liquid-in-vertical-periodic-motion ; RSI-Patterning-of-particulate-films-using-Faraday-waves .

We will introduce a new mode number $m$, the surface mode number, and propose a new dispersion relation, in which the oscillation frequency of the water drop depends on both mode numbers $n$ and $m$. The stiffness coefficient is reduced due to the additional surface mode, so that the eigen frequencies are lower than that 1 predicts. We also conduct experiments by exciting star-shaped oscillating drops on a hydrophobic vibrating substrate. Our dispersion relation can predict the oscillation frequencies with high accuracy.

II Theory

We consider the oscillation of water drops on a hydrophobic substrate, which vibrates periodically and serves as a vertical excitation. A water drop on a stable hydrophobic substrate is approximately a flattened cylinder with thickness about twice the capillary length ($2l_{c}$), determined by a balance between the surface tension and the gravity Physics-of-Fluids-Leidenfrost-drops .

We take cylindrical coordinates $(r,\theta,z)$ with respect to the frame of reference moving with the vibrating substrate, in which the undisturbed upper surface and periphery of the drop are $z=H$ and $r=R$ respectively. For an oscillating star-shaped drop, we use $h(r,\theta,t)$ and $a(z,\theta,t)$ to describe the small deviation of the upper surface and the periphery, hence the boundaries of the drop are described using $z=H+h$, $r=R+a$.

For irrotational liquid the velocity potential $\phi(r,\theta,z,t)$  An-Introduction-to-Fluid-Dynamics satisfies Eulerian equation

[TABLE]

where $\nabla\phi$ is the fluid velocity, $p$ is the internal pressure, $U$ is the potential of external forces. For small movements of the water drop, we linearize the equation by omitting the term $(\nabla\phi)^{2}$.

The pressure at the free liquid surface Hydrodynamics is

[TABLE]

where $\sigma$ is the surface tension and $\gamma_{1}$, $\gamma_{2}$ are the principal radius of curvature (ROC) of the surface.

The ROCs of the upper surface can be obtained from membrane theory rayleigh1894theory , by which we obtain $p=\sigma\nabla^{2}h$. So the boundary condition of the upper surface is

[TABLE]

Applying 3 to the periphery, we can get the boundary condition as

[TABLE]

The bottom surface is in contact with the hydrophobic surface, thus the boundary condition is

[TABLE]

Furthermore, at the liquid surfaces, $z=H+h$ and $r=R+a$, there are kinematic boundary conditions

[TABLE]

Considering the incompressibility, the equation of continuity can be written as

[TABLE]

For small deformations, Eq. 9 and Eqs. 4, 5, 6 are variable-separable. The oscillation of a flattened cylinder drop can be decomposed on the Bessel basis. Hence the general solution of the velocity potential is

[TABLE]

where $J_{n}$ refers to Bessel function. Here $A_{k,n}(t)$ represents the time-related coefficient. Similarly, $h(r,\theta,t)$ and $a(z,\theta,t)$ can be expressed as

[TABLE]

Inserting Eqs. 10, 11 and 12 into Eqs. 7 and 8, the coefficients $A_{k,n}(t)$, $B_{k,n}(t)$ and $C_{k,n}(t)$ satisfy

[TABLE]

Then substituting Eqs. 10, 11 and 12 to upper surface motion equation Eq. 4, we obtain

[TABLE]

It’s a parametric resonance equation, describing the response of the drop surface to the vertical excitation $Fcos\Omega t$. Here $k$ can be regarded as a wave-vector of the upper surface. If we denote $p$ and $q$ as follow

[TABLE]

Then Eq. 15 can be transformed to the standard form of Mathieu’s equation

[TABLE]

The natural oscillating frequency $\omega$ of the system can be obtained by setting $\Omega=0$. The principle subharmonic resonance occurs when Stability-chart-for-the-delayed-Mathieu-equation

[TABLE]

In practical cases, $q$ is much smaller than $p$ so thagt the resonance condition is simplified as $\omega\approx\Omega/2$. When the driving frequency is twice of the system natural frequency, the free surface becomes unstable and the surface oscillation mode increases in amplitude until it is finally restricted by damping or non-linear effects.

On the lateral surface, substituting Eqs. 10 and 12 into motion equation 5, we obtain

[TABLE]

Equation 20 defines a harmonic oscillating frequency of the drop periphery.

From 13 and 14 we know that, the induced surface motion due to parametric instability drives azimuthal motion in the same frequency $\omega$ and an azimuthal mode whose eigen frequency is in the vicinity of $\omega$ will be excited. Apparently, the two frequencies defined in Eqs. 15 and 20 should be equal, thus we obtain the eigen equation of $k$

[TABLE]

The eigen value $k$ determines oscillation patterns of the surface. We use $F(k)$ and $G(k)$ to represent the left and right side of the equation 21, which are plotted together in figure 1. We find that $F(k)$ increases monotonically with $k$ while $G(k)$ changes quasi-periodically. $F(k)$ and $G(k)$ have multiple intersections, and the $m$ $th$ nonzero intersection determines the $m$ $th$ eigen value $k_{m}$, which can be solved numerically. $m$ also represents $m$ $th$ eigen mode of the upper surface. The pattern of the water drop is presented in figure 2, in which the height variation $h(r,\theta)$ of the upper surface is denoted with different colors for a $n=6$ drop with $m$ increasing from $1$ to $4$. Along the radial direction, the number of wave nodes is proportional to surface mode number $m$.

From Eqs. 20 and 21, $n$ and $m$ determine the resonance frequency of the water drop, so we have the dispersion relation

[TABLE]

When the surface motion is ignored, i.e. $k=0$, we have the following equation

[TABLE]

Then equation 22 naturally reduces to Rayleigh equation 1. We define a factor $S_{n,m}$

[TABLE]

It reflects the influence of the surface mode on the resonance frequency.

In figure 1, we can find that $m=1$ and higher modes $m$ have different oscillation frequencies. In practical cases, each oscillation frequency $\omega$ should correspond to one mode $m$. Usually, for the upper surface of the drop, higher-order modes are suppressed due to the dissipation of the system and an appreciable disturbance at the surface occurs only when a low-order mode is excited The-Royal-Society-The-stability-of-the-plane-free-surface-of-a-liquid-in-vertical-periodic-motion . We numerically calculate $k_{m}$ and $S_{m,n}$ for $m=1$ and with $n$ increasing from $2$ to $10$. As shown in figure 3, for drops with the same radius $R=10.8mm$, the larger $n$, the larger the wave vector $k_{m}$, and the smaller the factor $S_{n,m}$ (all of them are smaller than $1$), which corresponds to a more significant softening to resonance frequency.

For an oscillation system, the oscillating frequency is determined by the stiffness coefficient. The existence of surface modes would decrease the surface tension potential, consequently decreases the stiffness coefficient for the oscillation. Ignoring surface modes and reducing the motion of the water drop to two-dimensional overestimates the potential energy as well as the resonance frequency. We will compare the frequency to experimental results in the following section. And the correction also explains the discrepancy of oscillation frequency found in Ref. PRE-Oscillating-and-star-shaped-drops-levitated-by-an-airflow .

III Experiment and results

Experimental setup is shown in Figure 4. A water-repellent cloth is attached horizontally to a loudspeaker cone (So-Voioe SVF149WR). For each drop the measured contact angle is more than $120^{\circ}$ so that the hysteretic behavior can be avoided. The loudspeaker is connected to a signal generator (Right SG1020P), which applies a vertical excitation to the drop. We use a injector to control the volume and place the water drop on the cloth. The top view of the water drop is recorded using a high-speed video camera (Metalab 300C-U3) at the rate of 400 frames/s. When a drop performs stationary star-shaped oscillation, the camera records a sequence of images. The oscillating frequency $f$ and mean radius $R$ of the drop are analyzed frame by frame.

We set the input sinusoidal waves at the driving frequency $f_{0}=\Omega/2\pi=60Hz,90Hz,120Hz$ respectively. The volume of the water drop is increased by $0.05ml$ each time, measured by the injector. For $f_{0}=120Hz$, star-shaped oscillating drops from $n=3$ to $n=11$ are observed, as shown in figure 5. We can see the petal-like patterns at the upper surface, illustrating the existence of surface motion patterns. For $f_{0}=90Hz$, star-shaped drops $n=3,5,6,7,8,9$ are observed; for $f_{0}=60Hz$, star-shaped drops $n=4,5,6,7$ are observed. The surface mode computed with our model is compared to the experimental photograph for a $n=3$ drop in figure 6. The radius of the drop is $2.7mm$ and the calculated wave-vector $k_{1}$ is $0.99mm^{-}1$. The drop is captured from very close distance and we enhance the image contrast for visibility. The theoretical image resembles the experimental one. We show the relation between the measured oscillating frequency $f$ and half of the excitation frequency $f_{0}/2$ in figure 7. The subharmonic parametric resonance condition 19 is satisfied.

In dispersion relation 1, the predicted value of $\omega^{2}$ varies as $n(n^{2}-1)$, while in the new dispersion relation 22 we propose, it varies as $n(n^{2}-1)$ multiplied by an additional factor $S_{n,m}$. For each azimuthal mode $n$, we calculate the eigen value of $k_{m}$ by solving equation 21 numerically and we take mode number $m=1$. For $f_{0}=120Hz$,$90Hz$,$60Hz$, $k_{1}$ are calculated and equal to $1.4mm^{-1}$,$1.0mm^{-1}$,$0.8mm^{-1}$ respectively. We find $k_{1}$ appears to be independent of the azimuthal mode $n$, only determined by the excitation frequency $f_{0}$. In figure 8 we plot the theoretical results according to 1 and 22, as well as the experimental results of $\omega^{2}R^{3}\sigma/\rho$. The yellow dots are experimental values. The blue dashed lines represent the values predicted by 1, obviously overestimating the oscillating frequency. And the orange solid lines represent those predicted by 22, which fits the experimental data much better. For some azimuthal number $n$ the observed oscillation frequency consist of a set of data points. Actually these drops have the same azimuthal mode $n$ but are slightly different in radius $R$ and surface factor $S_{n,m}$. while we only substitute in 22 the average R observed in experiment.

IV Conclusion

In this paper, we derive a complete theoretical model to include both surface modes and azimuthal modes of water drops under vertical excitation. The model we propose is applicable to drops under any form of vertical excitation. We prove that the star-shaped oscillation originates from parametric instability of the upper surface and explain the mechanism that leads to the star-shaped oscillation. We also propose a new dispersion relation based on the combination of surface modes and azimuthal modes, which gives a lower resonance frequency due to the additional surface mode. The dispersion relation explains the discrepancy of oscillation frequencies found in previous studies and is in good agreement with our experimental data. These results enhance our understanding of the dynamics of water drops.

Acknowledgements

The authors are grateful to Mr. Y. Luo for beneficial discussions.