Abstract
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 secondorder integrodifferential equation. Some limiting cases are recovering from the present general one. The selfgravitating 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 nonaxisymmetric perturbation, the selfgravitating force is stabilizing for short and long waves, while the electric field has stabilizing effect on short waves.
Keywords:
electrogravitational stability; oscillating; streaming1. Introduction
The stability of selfgravitating fluid cylinder has been studied, for the first time, by Chandrasekhar and Fermi [1]. Later on, Chandrasekhar [2] made several extensions as the fluid cylinder is acted by different forces. Radwan [3,4] studied the stability of an ideal hollow jet. Radwan [4] considered that the fluids are penetrated by constant and uniform electric fields. The stability of different cylindrical models under the action of selfgravitating force in addition to other forces has been elaborated by Radwan and Hasan [5,6]. Radwan and Hasan [5] studied the gravitational stability of a fluid cylinder under transverse timedependent electric field for axisymmetric perturbations. Hasan [7,8] has discussed the stability of oscillating streaming fluid cylinder subject to combined effect of the capillary, selfgravitating, and electrodynamic forces for all axisymmetric and nonaxisymmetric perturbation modes. Hasan [7,8] studied the instability of a full fluid cylinder surrounded by selfgravitating tenuous medium pervaded by transverse varying electric field under the combined effect of the capillary, selfgravitating, 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 magnetoflow meters, etc.
In this article, we aim to investigate the stability of oscillating streaming selfgravitating dielectric incompressible fluid cylinder surrounded by tenuous medium of negligible motion pervaded by transverse varying electric field for all the axisymmetric and nonaxisymmetric perturbation modes.
2. Mathematical formulation
Consider a selfgravitating fluid cylinder surrounded by a selfgravitating medium of negligible motion. The cylinder of (radius R_{0}) dielectric constant ε^{(i) }while the surrounding medium is being with dielectric constant ε^{(e)}. Fluid is assumed to be incompressible, inviscid, selfgravitating, and pervaded by applied longitudinal electric field.
The surrounding tenuous medium (being of negligible motion), selfgravitating, and penetrated by transverse varying electric field
where E_{0 }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 zaxis 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 zaxis 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 selfgravitating equations concerning the selfgravitating matter (see [2,710]).
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 selfgravitating 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 1214 are solved and moreover the solutions are matched across the fluid cylinder interface at r = R_{0}. The nonsingular solution in the unperturbed state is, finally, given as
4. Linearization
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
with
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 zdirection while m, an integer, is the azimuthal wave number.
The nonsingular solutions of the linearized perturbation equations give ϕ,V, and ψ as follows:
where A_{1}(t), B_{1}(t), B_{2}(t), C_{1}(t), and C_{2}(t) are arbitrary functions of integrations to be determined, while I_{m}(kr) and K_{m}(kr) are the modified Bessel functions of the first and second kind of order m.
5. Boundary conditions
The nonsingular 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 = R_{0}. 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 R_{0 }is, dimensionless, the longitudinal wave number.
(ii) Selfgravitating conditions
The gravitational potential V = V_{0 }+ εV_{1 }+ ⋯ and its derivative must be continuous across the perturbed boundary fluid surface at r = R_{0}. 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 = R_{0}. These conditions could be written in the form
While is, the outward unit vector normal to the interface (18) at r = R_{0}, given by
So that
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 = R_{0}. 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 integrodifferential equation governing the surface displacement η*(t). By means of this relation, we may identify the (in) stability states and also the selfgravitating 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
where
Equation 44 has the canonical form
where
Equation 47 is Mathieu differential equation. The properties of the Mathieu functions are explained and investigated by Melaclan [11]. 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 [11], 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 integrodifferential 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 [11] 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 [12] which is valid only for small values of h^{2 }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 h^{2 }and could be written as
where α_{1 and }α_{2 }are, the two roots of the equality of the relation (51), being
with
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
or alternatively
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 (h^{2 }α_{1}) and (h^{2 }α_{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 h^{2 } c and α_{1,2 }= c_{1 } ic_{2 }where c, c_{1}, and c_{2 }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 E_{0 }= 0 and we suppose that
The secondorder integrodifferential 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 [1]. In fact, they [1] 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
and
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
is identified.
Here, it is found that the quantity Q_{0 }(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 selfgravitating destabilizing influence of the dielectric fluid cylinder dispersed in a dielectric medium of negligible motion.
However, the selfgravitating 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 secondorder integrodifferential equation of Mathieu equation (47) yields
where
and
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).
8. Conclusion
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.
Competing interests
The authors declare that they have no competing interests.
Appendix I
Acknowledgements
We are grateful to the Editor of the Journal and the Reviewers for their suggestions and comments on this article.
References

Chandrasekhar, S, Fermi, E: Problems of gravitational stability in the presence of a magnetic field. Astrophys J. 118, 116–141 (1953)

Chandrasekhar, S: Hydrodynamic and Hydromagnetic Stability. Dover, New York (1981)

Radwan, AE: MHD stability of gravitational compressible fluid cylinder. J Phys Soc Jpn. 58, 1225–1227 (1989). Publisher Full Text

Radwan, AE: Electrogravitational instability of an annular fluid jet coaxial with a very dense fluid cylinder under radial varying fields. J Plasma phys. 44, 455–465 (1991)

Radwan, AE, Hasan, AA: Axisymmetric electrogravitational stability of fluid cylinder ambient with transverse varying oscillating field, (IAENG). Int J Appl Math. 38(3), 113–120 (2008)

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

Hasan, AA: Electrogravitational stability of oscillating streaming fluid cylinder. J Appl Mech ASME (2011, in press)

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 conejet at high Reynolds number. J Aerosol Sci. 25(6), 1037–1047 (1994). Publisher Full Text

Melaclan, N: Theory and Application of Mathieu Function. Oxford University Press, New York (1964)

Morse, PM, Feshbach, H: Methods of Theoretical Physics Part I. McGrawHill, New York (1953)