Analytical study of nonlinear vibrations of a charged viscous liquid drop

The generatrix of a nonlinearly vibrating charged drop of a viscous incompressible conducting liquid is found by directly expanding the equilibrium spherical shape of the drop in the amplitude of initial multimode deformation up to second-order terms. A fact previously unknown in the theory of nonlinear interaction is discovered: the energy of an initially excited vibration mode of a low-viscosity liquid drop is gradually (within several vibrations periods) transferred to the mode excited by only nonlinear interaction. Irrespectively of the form of the initial deformation, an unstable viscous drop bearing a charge slightly exceeding the critical Rayleigh value takes the shape of a prolate spheroid because of viscous damping of all the modes (except for the fundamental one) for a characteristic time depending on the damping rates of the initially excited modes and the further evolution of the drop is governed by the fundamental mode. In a high-viscosity drop, the rate of rise of the unstable fundamental mode amplitude does not increase continuously with time, contrary to the predictions of nonlinear analysis in terms of the ideal liquid model: it first decreases to a value slightly differing from zero (which depends on the extent of supercriticality of the charge and viscosity of the liquid), remains small for a while (the unstable mode amplitude remains virtually time-independent), and then starts growing.     

Capillary oscillations and stability of a charged drop rotating about the axis of symmetry

The stability of a charged conductive liquid drop rotating about the axis of symmetry against the pressure of the self-charge electric field and inertial force pressure is investigated in an approximation linear in oscillation amplitude and square of the spheroidal drop deformation eccentricity. It is found that the axisymmetric modes of the rotating drop are stable. Only nonaxisymmetric modes with azimuthal numbers maximal for a given mode may be unstable. The Coriolis force plays a stabilizing role.     

Instability of a spherical drop in the field of an electrical dipole

Analytical calculations show that, as a field in which an initially spherical charged conducting incompressible drop is placed becomes more and more nonuniform, coupling between the drop’s oscillation modes grows and the threshold of stability against the electrical field pressure declines. When an electrostatic parameter characterizing the electric field pressure exceeds a value that is critical for a certain mode to be unstable, the amplitude of this mode exponentially grows in an aperiodic manner and the amplitudes of modes coupled with this mode build up in an oscillatory manner, each mode having its own instability growth rate. In all cases, there exists a threshold value of the dimensionless electric parameter above which all oscillation modes are unstable.     

Nonlinear oscillations of a drop moving at a constant velocity in an insulating medium subjected to an electrostatic field

Nonlinear asymptotic calculations of the second order of smallness in the amplitude of the initial deformation of an ideally conducting liquid drop show that the laminar flow of an ideal conducting incompressible dielectric liquid flowing about the drop in an external electrostatic field parallel to the flow causes oscillation mode’s interaction in the first and second orders of smallness. Both linear and nonlinear interactions between the oscillation modes of the drop excite modes that are absent in the spectrum of modes governing the initial deformation of the drop’s equilibrium shape. In the second order of smallness, the mode interaction decreases the electrostatic field strength, as well as the velocity and density of the environment, that are critical for development of instability of the drop against the polarization charge.     

On the influence of viscosity on the characteristic time of the instability of a charged drop

A nonlinear integral equation describing the evolution of spheroidal deformation of a drop that is unstable with respect to its intrinsic charge is derived and solved for arbitrary values of viscosity. It was shown that, due to an essentially nonlinear character of the phenomenon, the characteristic time of instability develop-ment equals the time of tenfold increase in the amplitude of an initial, physically infinitesimal spheroidal deformation of an unstable drop. The dependence of the instability characteristic time on the drop viscosity is described by an increasing linear function.     

Linear (in oscillation dimensionless amplitude) interaction between the modes of a nonspherical charged drop in an external electrostatic field

The stability of a heavily charged drop in a weak uniform electrostatic field (in which the equilibrium shape of the drop can be represented by a prolate spheroid) is calculated in the fourth order of smallness in the eccentricity of the spheroidal drop and in the first order of smallness in the drop oscillation dimensionless amplitude. It is found that as the order of approximation in eccentricity grows, so does the number of modes interacting with the initially excited mode. In the given order of smallness, the preferred (initially excited) mode is shown to interact with the nearest eight modes. The drop becomes unstable if such is the second mode.     

Nonlinear secondary Raman resonances of oscillations of the fundamental mode of a charged drop moving in an environment

Second-order calculations show that, when a gas flows about a charged drop, the fundamental mode of the multimode initial deformation of its equilibrium shape builds up through nonlinear secondary Raman resonant interaction with higher modes if this mode is present in the mode spectrum specifying the initial deformation. This circumstance accounts for large-amplitude spheroidal oscillations of drops in natural liquid-drop systems and provides an insight into corona initiation in the vicinity of drops in thunderstorm clouds and into lightning initiation.     

Instability of a spherical drop in a nonuniform electric field

Analytical calculation in the first order of smallness shows that the equilibrium shape of a drop in the field of a point charge is axisymmetric about the plane passing through the center of mass of the drop normally to the axis connecting the center of mass with the point charge. Whether the equilibrium shape of the drop is stable or not depends on the value of the field parameter, which, in turn, depends on the point charge and the distance to it. There is an asymptotic value of the critical parameter above which all modes become unstable. In the field of the point charge, the mode coupling grows; that is, a mode excited at the zero time generates oscillations of the six nearest modes with amplitudes proportional to that of the initially excited mode. If the initially excited mode loses stability, all the modes coupled with it also become unstable. The surface instability of the drop also develops when the initially excited mode is stable but at least one of the modes coupled with it is unstable.     

Effect of the initial deformation of a charged drop on nonlinear corrections to critical conditions for instability

A nonlinear (proportional to the vibration amplitude squared) decrease in the critical (in terms of instability) charge of a vibrating drop is found to be limited, as follows from third-order asymptotic calculations. This effect occurs when the spectrum of modes specifying the initial deformation of the drop contains, along with the fundamental mode, higher modes. The influence of the environment density on nonlinear corrections to the critical conditions for instability is analyzed.     

Nonlinear resonance interaction between the oscillation modes of a spherical liquid layer on the surface of a melting hailstone

Evolutionary equations are derived and solved that describe the time dependence of the oscillation mode amplitudes on the surface of a charged conducting liquid layer resting on a solid core. It is assumed that the layer experiences a multimode initial deformation. The equations are solved asymptotically in the second order of smallness in the small dimensionless amplitude of capillary oscillations on the surface of the layer. Mechanisms behind internal nonlinear resonance interaction between the modes of the liquid layer oscillations and behind energy transfer between the modes both in degenerate and in secondary combination resonances are investigated. It is found that in the degenerate resonance interaction between oscillation modes, the energy may be transferred not only from lower to higher modes but also vice versa if the higher mode is excited at the zero time. This conclusion is valid not only for a liquid layer on the surface of a solid core but also for a drop.     

On nonlinear resonant four-mode interaction between capillary vibrations of a charged drop

Energy transfer from higher modes of capillary vibrations of an incompressible liquid charged drop to the lowest fundamental mode under four-mode resonance is studied. The resonance appears when the problem of nonlinear axisymmetric capillary vibration of a drop is solved in the third-order approximation in amplitude of the multimode initial deformation of the equilibrium shape of the drop. Although the resonant interaction mentioned above builds up the fundamental mode even in the first order of smallness, its amplitude turns out to be comparable to a quadratic (in small parameter) correction arising from nonresonant nonlinear interaction, since the associated numerical coefficients are small.     

Nonlinear oscillations and stability of a charged drop moving relative to a dielectric medium

Nonlinear calculations to within the second order of smallness with respect to the initial deformation of a liquid drop show that a stream of an ideal incompressible dielectric liquid streamlining the charged ideally conducting drop causes interaction between modes both in the first and second orders of smallness. Both the linear and nonlinear interactions of the oscillation modes result in the excitation of modes absent in the spectrum of the initial drop deformation. The relative motion of the drop and the medium leads to broadening of the spectrum of modes excited in the second order of smallness. The presence of the flow streamlining the drop and the intermode interaction result in decreasing the critical magnitudes of the drop charge and the velocity and density of the medium determining drop instability development.     

Critical conditions for instability of a highly charged oblate spheroidal drop

The spectrum of capillary oscillations of a charged oblate spheroidal drop is calculated in neglect of the interaction between modes by means of a perturbation expansion in the small deviation of the equilibrium shape of the drop from spherical. The critical conditions for instability of its nth mode with respect to the self-charge are calculated in the form of an analytical function describing how the dimensionless Rayleigh parameter characterizing the stability of the drop depends on the value of the spheroidal deformation. Zh. Tekh. Fiz. 69, 10–14 (July 1999)     

Nonlinear oscillations of a charged drop accelerated in an electrostatic field

An analytical asymptotic solution to the problem of nonlinear oscillations of a charged drop moving with acceleration through a vacuum in a uniform electrostatic field is found. The solution is based on a quadratic approximation in two small parameters: the eccentricity of the equilibrium spheroidal shape of the drop and the amplitude of the initial deformation of the equilibrium shape. In the calculations carried out in an inertial frame of reference with the origin at the center of mass of the drop, expansions in fractional powers of the small parameter are used. Corrections to the vibration frequencies are always negative and appear even in the second order of smallness. They depend on the stationary deformation of the drop in the electric field and nonlinearly reduce the surface charge critical for development of the drops’s instability. It is found that the evolutions of the shapes of nonlinearly vibrating unlike-charged drops differ slightly owing to inertial forces.     

Stability of a charged drop near a conductor wall

The effect of conductor boundaries on the deformation and stability of a charged drop is presented. The motivation for such a study is the occurrence of a charged conductor drop near a conductor wall in experiments (Millikan-like set-up in studies on Rayleigh break-up) and applications (such as electrospraying, ink-jet printing and ion mass spectroscopy). In the present work, analytical (linear stability analysis (LSA)) and numerical methods (boundary element method (BEM)) are used to understand the instability. Two kinds of boundaries are studied: a spherical, conducting, grounded enclosure (similar to a spherical capacitor) and a planar conducting wall. The LSA of a charged drop placed at the center of a spherical cavity shows that the Rayleigh critical charge (corresponding to the most unstable l = 2 Legendre mode) is reduced as the non-dimensional distance ?d = (b - a)/a decreases, where a and b are the radii of the drop and spherical cavity, respectively. The critical charge is independent of the assumptions of constant charge or constant potential conditions. The trans-critical bifurcation diagram, constructed using BEM, shows that the prolate shapes are subcritically unstable over a much wider range of charge as [Formula: see text] decreases. The study is then extended to the stability of a charged conductor drop near a flat conductor wall. Analytical theory for this case is difficult and the stability as well as the bifurcation diagram are constructed using BEM. Moreover, the induced charges in the conductor wall lead to attraction of the drop to the wall, thereby making it difficult to conduct a systematic analysis. The drop is therefore assumed to be held at its position by an external force such as the electric field. The case when the applied field is much smaller than the field due to inherent charge on the drop ((a(3)ρg)/(3ε(0)Ψ(2)) ? 1 is considered. The wall breaks the fore-aft symmetry in the problem, and equilibrium, predominantly prolate shapes corresponding to the legendre mode, l = 2 , are observed. The deformation increases with increasing charge on the drop. The breakup of the prolate equilibrium shapes is independent of the legendre modes of the initial perturbations. The prolate perturbations are subcritically unstable. Since the equilibrium prolate shapes cannot continuously exchange instability with equilibrium oblate shapes, an imperfect transcritical bifurcation is observed. A variety of highly deformed equilibrium oblate shapes are predicted by the BEM calculations.     

Deformation and nonlinear corrections to critical conditions of charged drop instability in an electrostatic suspension

Nonlinear asymptotic analysis of a charged drop placed in electrostatic and gravitational fields reveals a correction to the oscillation frequency and, accordingly, to the critical Rayleigh parameter. The analysis uses approximations quadratic in oscillation amplitude and linear in dimensionless equilibrium deformation of the drop. The correction is found to be proportional to the product of the oscillation amplitude and deformation. It is natural to name this correction deformational. In computations of the third order of smallness in oscillation dimensionless amplitude, a correction to the frequency and Rayleigh parameter appears, which is due to a nonlinear interaction between oscillation modes. This correction is larger than the deformational one in magnitude. Deformational corrections can be eliminated by experimenting under no-gravity conditions, but corrections due to the nonlinearity of hydrodynamic equations cannot be eliminated in this way. It is these corrections that are responsible for a critical Rayleigh parameter measurement inaccuracy.     

Nonlinear oscillations of an uncharged conducting drop in an external uniform electrostatic field

The problem of nonlinear oscillations of the finite amplitude of an uncharged drop of an ideal incompressible conducting liquid in an external uniform electrostatic field is solved for the first time by analytical asymptotic methods. The problem is solved in an approximation quadratic in amplitude of the initial deformation of the equilibrium shape of the drop and in eccentricity of its equilibrium spheroidal deformation. Compared with the case of nonlinear oscillations of charged drops in the absence of the field, the curvature of the vertices of uncharged drops nonlinearly oscillating in the field is noticeably higher, whereas the number of resonant situations (in the sense of internal resonant interaction of modes) is much smaller.     

On the structure of the velocity field under the free spheroidal surface of a viscous liquid drop oscillating in an electrostatic field

An asymptotic analytical solution to an initial boundary-value problem considering (i) the time evolution of the capillary oscillation amplitude as applied to a viscous spheroidal liquid drop placed in a uniform electrostatic field and (ii) the liquid flow velocity field inside the drop is found. The problem is solved in an approximation that is linear in two small parameters: the dimensionless oscillation amplitude and the dimensionless field-induced constant deformation of the equilibrium (spherical) shape of the drop. Terms proportional to the product of the small parameters are retained. In this approximation, interaction between oscillation modes is revealed. It is shown that the intensity of the eddy component of the oscillation-related velocity field depends on the liquid viscosity and the external uniform electrostatic field strength. The intensity of the eddy component decays rapidly with distance from the free surface. The depth to which the eddy flow (which is caused by periodical flows on the free surface) penetrates into the drop is a nonmonotonic function of the polar angle and increases with dimensionless viscosity and field strength.     

Influence of charge relaxation on the capillary oscillations of a charged viscous spheroidal drop

A dispersion relation is derived for the spectrum of capillary modes of a charged spheroidal drop of a viscous liquid with allowance for charge relaxation. It is shown that the finite charge transport rate leads to lowering of the instability growth rates for various capillary modes of a spheroidal drop of a low-viscosity liquid. As the degree of deformation of the drop increases, the magnitude of the absolute change in the growth rate caused by the finite rate of charge redistribution decreases. Zh. Tekh. Fiz. 69, 28–36 (August 1999)     

Influence of the viscosity of a nonlinearly vibrating charged liquid drop on the positions of internal resonances

An asymptotic analytical expression for the generatrix of a viscous charged liquid drop is for the first time derived in the second order of smallness in the axisymmetric initial deformation of the drop. The expression is represented as an infinite series in the roots of the dispersion relation and a finite sum of the numbers of modes specifying the initial deformation. In some of the terms of the analytical expression, the denominators involve the differences between the mode frequencies. These differences may become small under certain values of the charge, causing internal nonlinear resonant mode interaction. Analytical and numerical investigations of the effect of viscosity on the vibrating frequency show that the resonant values of the self-charge of the drop tend to increase with increasing viscosity. The viscosity of the liquid does not affect the spectrum of modes excited via nonlinear mode interaction.