The electrogravitational instability of a dielectric oscillating streaming fluid cylinder surrounded by tenuous medium of negligible motion pervaded by transverse varying electric field has been investigated for all the perturbation modes. The model is governed by Mathieu second-order integro-differential equation. Some limiting cases are recovering from the present general one. The self-gravitating force is destabilizing only in the axisymmetric perturbation for long wavelengths, while, the axial electric field interior, the fluid has strong destabilizing effect for all short and long wavelengths. The transverse field is strongly stabilizing. In the case of non-axisymmetric perturbation, the self-gravitating force is stabilizing for short and long waves, while the electric field has stabilizing effect on short waves.
Keywords:electrogravitational stability; oscillating; streaming
The stability of self-gravitating fluid cylinder has been studied, for the first time, by Chandrasekhar and Fermi . Later on, Chandrasekhar  made several extensions as the fluid cylinder is acted by different forces. Radwan [3,4] studied the stability of an ideal hollow jet. Radwan  considered that the fluids are penetrated by constant and uniform electric fields. The stability of different cylindrical models under the action of self-gravitating force in addition to other forces has been elaborated by Radwan and Hasan [5,6]. Radwan and Hasan  studied the gravitational stability of a fluid cylinder under transverse time-dependent electric field for axisymmetric perturbations. Hasan [7,8] has discussed the stability of oscillating streaming fluid cylinder subject to combined effect of the capillary, self-gravitating, and electrodynamic forces for all axisymmetric and non-axisymmetric perturbation modes. Hasan [7,8] studied the instability of a full fluid cylinder surrounded by self-gravitating tenuous medium pervaded by transverse varying electric field under the combined effect of the capillary, self-gravitating, and electric forces for all the modes of perturbations.
There are many applications of electrohydrodynamic and magnetohydrodynamic stability in several fields of science such as
1. Geophysics: the fluid of the core of the Earth and other theorized to be a huge MHD dynamo that generates the Earth's magnetic field because of the motion of the liquid iron.
2. Astrophysics: MHD applies quite well to astrophysics since 99% of baryonic matter content of the universe is made of plasma, including stars, the interplanetary medium, nebulae and jets, stability of spiral arm of galaxy, etc. Many astrophysical systems are not in local thermal equilibrium, and therefore require an additional kinematic treatment to describe all the phenomena within the system.
3. Engineering applications: there are many forms in engineering sciences including oil and gas extraction process if it surrounded by electric field or magnetic field, gas and steam turbines, MHD power generation systems and magneto-flow meters, etc.
In this article, we aim to investigate the stability of oscillating streaming self-gravitating dielectric incompressible fluid cylinder surrounded by tenuous medium of negligible motion pervaded by transverse varying electric field for all the axisymmetric and non-axisymmetric perturbation modes.
2. Mathematical formulation
Consider a self-gravitating fluid cylinder surrounded by a self-gravitating medium of negligible motion. The cylinder of (radius R0) dielectric constant ε(i) while the surrounding medium is being with dielectric constant ε(e). Fluid is assumed to be incompressible, inviscid, self-gravitating, and pervaded by applied longitudinal electric field.
The surrounding tenuous medium (being of negligible motion), self-gravitating, and penetrated by transverse varying electric field
where E0 is the intensity of the electric field in the fluid while β is some parameters satisfy certain conditions. The components of and are considered along the utilizing cylindrical coordinates (r, φ, z) system with z-axis coinciding with the axis of the fluid cylinder. The fluid of the cylinder streams with a periodic velocity
where ω is constant and U is the speed at time t = 0.
The components of electric fields and are being along (r,φ,z) with the z-axis coinciding with the axis of the fluid cylinder (as shown in Figure 1).
Figure 1. Sketch for gravitational dielectric fluid cylinder.
The basic equations for investigating the problem under consideration are being the combination of the ordinary hydrodynamic equations, Maxwell equations concerning the electromagnetic theory, and Newtonian self-gravitating equations concerning the self-gravitating matter (see [2,7-10]).
For the problem under consideration, these equations are given as follows.
where ρ, , and P are the fluid density, velocity vector, and kinetic pressure, respectively, and and are the electric field intensity and self-gravitating potential of the fluid while and are these of tenuous medium surrounding the fluid cylinder, and G is the gravitational constant.
Since the motion of the fluid is irrotational, incompressible motion, the fundamental equations may be written as
where ϕ and ψ are the potential of the velocity of the fluid and electrical potential.
3. Equilibrium state
In this case, the basic equations are given in the form
where the subscript 0 here and henceforth indicates unperturbed quantities.
Equations 12-14 are solved and moreover the solutions are matched across the fluid cylinder interface at r = R0. The non-singular solution in the unperturbed state is, finally, given as
For a small wave disturbance across the boundary interface of the fluid, the surface deflection at time t is assumed to be of the form as
Consequently, any physical quantity Q(r,φ,z;t) may be expressed as
where η(t) is the amplitude of the perturbation at an instant time t, k, any real number, is the longitudinal wave number along z-direction while m, an integer, is the azimuthal wave number.
The non-singular solutions of the linearized perturbation equations give ϕ,V, and ψ as follows:
where A1(t), B1(t), B2(t), C1(t), and C2(t) are arbitrary functions of integrations to be determined, while Im(kr) and Km(kr) are the modified Bessel functions of the first and second kind of order m.
5. Boundary conditions
The non-singular solutions of the linearized perturbation equation given by the systems (21)-(25) and the solutions (16)-(17) of the unperturbed systems (12)-(14) must satisfy certain boundary conditions. Under the present circumstances, these appropriate boundary conditions could be applied as follows.
(i) Kinematic conditions
The normal component of the velocity vector must be compatible with the velocity of the boundary perturbed surface of the fluid at the level r = R0. This condition, yield
By the use of Equations 18, 19, and 21 for the condition (26), after straight forward calculations, we get
where x = k R0 is, dimensionless, the longitudinal wave number.
(ii) Self-gravitating conditions
The gravitational potential V = V0 + εV1 + ⋯ and its derivative must be continuous across the perturbed boundary fluid surface at r = R0. These conditions are given as
By utilizing Equations 18, 19, 22, and 23 for the conditions (28) and (29), we get
(iii) Electrodynamic condition
The normal component of the electric displacement current and the electric potential ψ perturbed boundary surface at the initial position r = R0. These conditions could be written in the form
While is, the outward unit vector normal to the interface (18) at r = R0, given by
Upon applying these conditions, we get
where the quantity ξ1 is given in Appendix 1.
(iv) The dynamical stress condition
The normal component of the total stress across the surface of the coaxial fluid cylinder must be continuous at the initial position at r = R0. This condition is given as follows
By substituting for and , after some algebraic calculations, we finally obtain
where the quantity β11 and β12 is given in Appendix I.
In order to eliminate the first derivative term, we may use the substitution
Equation 41 can be expressed as follows
Equation 43 is an integro-differential equation governing the surface displacement η*(t). By means of this relation, we may identify the (in-) stability states and also the self-gravitating and electrodynamic forces influences on the stability of the present model. However in order to do so, it is found more convenient to express this relation in the simple form
Equation 44 has the canonical form
Equation 47 is Mathieu differential equation. The properties of the Mathieu functions are explained and investigated by Melaclan . The solutions of Equation 47, under appropriate restrictions, could be stable and vice versa. The conditions required for periodicity of Mathieu functions are mainly dependent on the correlation between the parameters a and q. However, it is well known, see , that (a, q)-plane is divided essentially into two stable and unstable domains separated by the characteristic curves of Mathieu functions. Thence, we can state generally that a solution of Mathieu integro-differential equation is unstable if the point (a, q) say, in the (a, q)-plane lies internal and unstable domain, otherwise it is stable.
6. Discussions and limiting cases
The appropriate solutions of Equation 47 are given in terms of what called ordinary Mathieu functions which, indeed, are periodic in time t with period π and 2π.
Corresponding to extremely small values of q, the first region of instability is bounded by the curves
The conditions for oscillation lead to the problem of the boundary regions of Mathieu functions where Melaclan  gives the condition of stability as
where Δ(0) is the Hill's determinant.
An approximation criterion for the stability near the neighborhood of the first stable domains of the Mathieu stability domains given by Morse and Feshbach  which is valid only for small values of h2 or q, i.e., the frequency ω of the electric field is very large.
This criterion, under the present circumstances, states that the model is ordinary stable if the restriction
is satisfied where the equality is corresponding to the marginal stability state. The inequality (51) is a quadratic relation in h2 and could be written as
where α1 and α2 are, the two roots of the equality of the relation (51), being
The electrogravitational stability and instability investigations analysis should be carried out in the following two cases
(i). 0 < b < 2/3
In this case Δ2 is positive and therefore the two roots α1 and α2 of the equality (51) are real. Now, we will show that both α1 and α2 are positive. If α1 α + ve then α1 must be negative and this means that
From which we get
and this is contradiction, so α1 must be positive and consequently α2 ≥ 0 as well (noting that α2 > α1). This means that both the quantities (h2 -α1) and (h2 -α2) are negative and that in turn show that the inequality (51) is identically satisfied.
(ii). 2/3 < b < 1
In this case, in which b < 1 and simultaneously 3b > 2, it is found that Δ2 is negative, i.e., Δ is imaginary; therefore, the two roots α1 and α2 are complex. We may prove that the inequality (51) is satisfied as follows.
Let h2 - c and α1,2 = c1 - ic2 where c, c1, and c2 are real, so
which is positive definite.
By an appeal to the cases (i) and (ii), we deduce that the model is stable under the restrictions
This means that the model is stable if there exists a critical value ω0 of the electric field frequency ω such that ω > ω0 where ω0 is given by
One has to mention here that if ω = 0, β = 0, and E0 = 0 and we suppose that
The second-order integro-differential equation of Mathieu equation (41) yields
where σ is the temporal amplification and note by the way that has a unit of time. The relation (62) is identical to the gravitational dispersion relation derived for the first time by Chandrasekhar and Fermi . In fact, they  have used a totally different technique rather than that used here. They have used the method of representing the solenoidal vectors in terms of poloidal and toroidal vector fields for axisymmetric perturbation.
To determine the effect of ω, it is found more convenient to investigate the eigenvalue relation (62) since the right side of it is the same the middle side of (60).
Taking into account the recurrence relation of the modified Bessel's functions and their derivatives, we see, for x α 0, that
based on the values of x.
Now, returning to the relation (62), we deduce that the determining of the sign σ2/(4πGρi) is identified if the sign of the quantity
Here, it is found that the quantity Q0 (x) may be positive or negative depending on x α 0 values. Numerical investigations and analysis of the relation (62) reveal that σ2 is positive for small values of x while it is negative in all other values of x. In more details, it is unstable in the domain 0 < x < 1.0667 while it is stable in the domains 1.0667 ≤ x < ∞ where the equality is corresponding to the marginal stability state.
From the foregoing discussions, investigations, and analysis, we conclude (on using (65) for (62)) that the quantity
has the following properties
Now, returning to the relation (60) concerning the frequency ω0 of the periodic electric field
Therefore, we deduce that the electrodynamic force (with a periodic time electric field) has stabilizing influence and could predominate and overcoming the self-gravitating destabilizing influence of the dielectric fluid cylinder dispersed in a dielectric medium of negligible motion.
However, the self-gravitating destabilizing influence could not be suppressed whatever is the greatest value of the magnitude and frequency of the periodic electric field because the gravitational destabilizing influence will persist.
7. Numerical discussions
If we assume that ω = 0 and consider the condition (61), then the second-order integro-differential equation of Mathieu equation (47) yields
To verify and confirm the foregoing analytical results, the relation (69) has been inserted in the computer and computed. This has been done for several values of β as β < 1, β = 1, and β > 1 in the wide domain 0 ≤ x ≤ 0.5. The numerical data of instability corresponding and those of stability corresponding to are collected and tabulated and presented graphically (see Figures 2, 3, 4, 5, and 6). There are many features and properties in this numerical presentation as we see in the following:
Figure 2. Electrogravitational stable and unstable domains for β = 0.5.
Figure 3. Electrogravitational stable and unstable domains for β = 1.0.
Figure 4. Electrogravitational stable and unstable domains for β = 1.5.
Figure 5. Electrogravitational stable domains for β = 2.5.
Figure 6. Electrogravitational stable domains for β = 3.0.
(i) For β = 0.5 corresponding to M = 0.1, 0.3, 0.5, 0.7, 1.0, and 1.5 it is found that the electrogravitational unstable domains are 0 < x < 1.1175, 0 < x <1.19759, 0 < x < 1.27235, 0 < x 1.29599, 0 < x < 1.362741, and 0 < x < 1.3978, the neighboring stable domains are 1.1175 ≤ x < ∞, 1.19759 ≤ x < ∞, 1.27235 ≤ x < ∞, 1.29599 ≤ x < ∞, 1.362741 ≤ x < ∞, and 1.3978 ≤ x < ∞, where the equalities correspond to the marginal stability states (see Figure 2).
(ii) For β = 1.0 corresponding to M = 0.1, 0.3, 0.5, 0.7, 1.0, and 1.5 it is found that the electrogravitational unstable domains are 0 < x < 1.22669, 0 < x < 1.5266, 0 < x < 1.750969, 0 < x < 1.90513, 0 < x < 2.05422, and 0 < x < 2.19341, the neighboring stable domains are 1.22669 ≤ x < ∞, 1.5266 ≤ x < ∞, 1.750969 ≤ x < ∞, 1.90513 ≤ x < ∞, 2.05422 ≤ x < ∞, and 2.19341 ≤ x < ∞, where the equalities correspond to the marginal stability states (see Figure 3).
(iii) For β = 1.5 corresponding to M = 0.1, 0.3, 0.5, 0.7, 1.0, and 1.5 it is found that the electrogravitational unstable domains are 0 < x < 1.35924, 0 < x < 1.9735, 0 < x < 2.3982, 0 < x < 2.6563, 0 < x < 2.8835, and 0 < x < 3.0798, the neighboring stable domains are 1.35924 ≤ x < ∞, 1.9735 ≤ x < ∞, 2.3982 ≤ x < ∞, 2.6563 ≤ x < ∞, 2.8835 ≤ x < ∞, and 3.0798 ≤ x < ∞, where the equalities correspond to the marginal stability states (see Figure 4).
(iv) For β = 2.5, corresponding to M = 0.1, 0.3, 0.5, 0.7, 1.0, and 1.5 it is found that the electrogravitational fluid cylinder is completely stable not only for short wavelengths, but also for very long wavelengths and the gravitational unstable domains are completely suppressed (see Figure 5).
(v) For β = 3.0, corresponding to M = 0.1, 0.3, 0.5, 0.7, 1.0 and 1.5 it is found that the electrogravitational fluid cylinder is completely stable not only for short wavelengths, but also for very long wavelengths and the gravitational unstable domains are completely suppressed (see Figure 6).
From the presented numerical results, we may deduce the following. For the same value of M, it is found that the unstable domains are increasing with increasing of β values. This means that the influence of electric field has a destabilizing effect for all short and long wavelengths.
If β > 2.0, then the model is completely stable not only for short wave lengths, but also for long wave lengths.
The authors declare that they have no competing interests.
We are grateful to the Editor of the Journal and the Reviewers for their suggestions and comments on this article.
Radwan, AE: MHD stability of gravitational compressible fluid cylinder. J Phys Soc Jpn. 58, 1225–1227 (1989). Publisher Full Text
Radwan, AE, Hasan, AA: Magnetohydrodynamic stability of selfgravitational fluid cylinder. Appl Math Model. 33(4), 2121–2131 (2009). Publisher Full Text
Hasan, AA: Electrogravitational stability of oscillating streaming fluid cylinder. Physica B. 406(2), 234–240 (2011). Publisher Full Text
Lin, SH: Oxygen diffusion in a spherical cell with nonlinear oxygen uptake kinetics. Theor Biol. 60, 449–457 (1976). Publisher Full Text
Mestel, AJ: The electrohydrodynamic cone-jet at high Reynolds number. J Aerosol Sci. 25(6), 1037–1047 (1994). Publisher Full Text