Abstract
The magnetohydrodynamic stability criterion of selfgravitating streaming fluid cylinder under the combined effect of selfgravitating, magnetic, and capillary forces has been derived. The results are discussed analytically and some data are verified numerically for different parameters of the problem. The magnetic and capillary forces are stabilizing, but the streaming is destabilizing while the selfgravitating is stabilizing or destabilizing according to restrictions. The stable and unstable domains are identified and, moreover, the influences of the magnetic and capillary forces on the selfgravitating instability of the model have been examined. Including the magnetic force together with selfgravitating force improves the instability of the model. However, the selfgravitating instability will never be suppressed whatever the effects of the MHD force stabilizing effects are.
Keywords:
selfgravitating; magnetohydrodynamic; capillary; streamingIntroduction
The stability of a fluid cylinder under the action of the capillary or/and other forces has received the attention of several researchers (Rayleigh [1], Yuen [2], Nayfeh and Hassan [3] and Kakutani et al.[4]. The effect of the electromagnetic Lorentz force on the capillary instability has been examined in several texts by the Nobel prize winner (1986) Chandrasekhar [5]. This has been done only for small axisymmetric perturbation and with a constant magnetic field. Radwan et al.[610] extended such interesting works by studying the magnetohydrodynamic stability of a liquid jet embedded into a tenuous medium for all axisymmetric and nonaxisymmetric modes of perturbation. The stability of different cylindrical models under the action of selfgravitating force in addition to other forces has been elaborated by Radwan and Hasan [9] and [10]. They [9] studied the gravitational stability of a fluid cylinder under transverse timedependent electric field for axisymmetric perturbations. Hasan [11] discussed the stability of oscillating streaming fluid cylinder subject to the combined effect of the capillary, selfgravitating, and electrodynamic forces for all axisymmetric and nonaxisymmetric perturbation modes. He [12] 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 modes of perturbations. In [13] Hasan et al. investigated the hydromagntic stability of a selfgravitational oscillating streaming fluid jet pervaded by azimuthal varying magnetic field for all axisymmetric and nonaxisymmetric modes of perturbation. He [14] discussed 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 modes of perturbations. He [15] studied the magnetodynamic stability of a fluid jet pervaded by a transverse varying magnetic field while its surrounding tenuous medium is penetrated by uniform magnetic field.
The present work is devoted to studying the magnetogravitodynamic stability of a streaming fluid cylinder and examining the influence of capillary and magnetic forces on the selfgravitating instability of the present models. This may be carried out, for all axisymmetric and nonaxisymmetric modes of perturbation, analytically and the results will be verified numerically.
1 Formulation of the problem
We consider a uniform cylinder of an incompressible inviscid fluid of radius surrounded by a tenuous medium of negligible motion. In the initial unperturbed state, the model is assumed to be streaming uniformly with velocities
and pervaded internally and externally by the magnetic fields
Here W and U are (constants) the speed of the fluid, is the intensity of the magnetic field in the fluid, and α is some parameter. The components of , , are considered along the cylindrical coordinates with the zaxis coinciding with the axis of the cylinder as shown in Figure 1. The fluid matter of the cylinder is acted upon by the combined effects of the selfgravitating, inertial, capillary, and magnetic forces. The surrounding tenuous medium of the fluid cylinder is acted upon by the selfgravitating and magnetic forces only.
Figure 1. Sketch for a gravitational MHD fluid cylinder.
The required basic equations for such kind of study may be obtained by combining the ordinary hydrodynamic equations and those of Maxwell’s concerning the electromagnetic field theory together with Newtonian gravitational field equations.
For the problem at hand, under the present circumstances, these equations are the following.
For the fluid, we have
The curvature pressure due to the capillary force is
with
where
is the boundary surface equation at time t, while is a unit outward vector normal to the surface, T is surface tension, and is pressure due to curvature.
For the surrounding tenuous medium, the basic equations are
Here ρ, , and P are the fluid mass density, velocity vector, and kinetic pressure, respectively; , are the magnetic field intensities and V, are selfgravitating potentials, respectively, inside and outside the fluid cylinder, μ is the magnetic field permeability coefficient and G is the gravitational constant.
2 Unperturbed state
The unperturbed state is studied and the fundamental quantities of such state could be obtained. Equation (1) together with equation (3) gives
from which, taking into account equation (5), we obtain .
By integrating this equation, we get
where C is a constant of integration to be determined.
The surface pressure due to the capillary force (cf. Chandrasekhar [5]) is given by
The selfgravitating potentials and of the unperturbed state satisfy
The nonsingular solutions of equations (17) and (18) in the cylindrical coordinates with cylindrical symmetries and are given by
where , , and are constants of integration to be determined. By applying the conditions that the selfgravitational potential V and its derivative must be continuous across the unperturbed boundary surface at and choosing since the potential inside the cylinder is zero, we get
Therefore,
Moreover, by applying the condition that the total pressure must be balanced across the boundary surface at , the distribution of the fluid pressure in the unperturbed state is given by
It is worth noting that in the absence of surface tension at the boundary surface
in order that
3 Perturbation analysis
We consider small departures from an unperturbed rightcylindrical shape of an incompressible fluid. Therefore a normal mode can be expressed uniquely in terms of the deformed surface. Hence we may assume that the deformed interface is described by
with
Here is the elevation of the surface wave measured from the unperturbed position, k (real number) is the longitudinal wave number, m (integer) is the transverse wave number. The amplitude of the perturbation is given by
where (=ε at ) is the initial amplitude and σ is the temporal amplification. If σ (, ) is imaginary, then is the oscillation frequency of the propagating wave in the fluid.
As the initial streaming state is perturbed, every physical quantity may be expanded as
Here Q stands for P, u, V, , H, , and while indicates the unperturbed quantity and is a small increment of Q due to disturbances.
In view of the expansion (31), the basic equations of motion (3)(13) in the perturbation state give
where equations (33) and (34) have been used to obtain equation (35). Based on the linear perturbation technique, the linearized quantity may be expressed as
By means of the expansion (41), equations (36) and (40) give the secondorder ordinary differential equation
where stands for and . The solution of equation (42) is given in terms of the ordinary Bessel functions of imaginary argument. For the problem under consideration, apart from the singular solution, the solutions of equations (36) and (40) are finally given by
Here and are the modified Bessel functions of the first and second kind of order m, while A and B are constants of integration to be determined.
Using the spacetime dependence (41) for equation (32), we get
with
Also, equation (35) yields
By combining equations (45) and (47), we get
where
is the Alfven wave frequency defined in terms of .
By taking the divergence of both sides of equation (48) and using equation (33), we obtain
Using the space dependence (41) for equation (50) and following similar steps for the resulting differential equation as has already been done for equations (36) and (40), the solution of equation (50) could be obtained. Therefore, the nonsingular solution for is given by
where is a constant of integration to be determined.
The pressure surface in the perturbed state due to the capillary force is determined from equation (37) along with (29) in the form
where x () is the dimensionless longitudinal wavenumber.
Now, equation (34) means that the magnetic field intensity in the perturbed state may be derived from a scalar function, say, such that
By combining equations (38) and (53), we get
Similarly, as it has been done for equation (50), equation (54) is solved and its finite solution is given by
where is a constant of integration to be determined upon applying boundary conditions.
4 Boundary conditions
The solution of the basic equations (3)(13) in the unperturbed state given by (23)(25) together with (1), (2) and (6) and in the perturbed state given by (43)(55) must satisfy appropriate boundary conditions. These boundary conditions must be applied across the perturbed interface (28) at the unperturbed boundary surface .
Under the present circumstances, these boundary conditions may be stated as follows.
(i) Selfgravitating conditions.
The gravitational potential and its derivative must be continuous across the perturbed fluid interface (28) at the unperturbed boundary . These conditions at read
By substituting from equations (23), (24), (29), (43) and (44) into the conditions (56) and (57), we get
from which we obtain
where x () is the dimensionless longitudinal wave number.
(ii) Kinematic condition.
The normal component of the velocity vector u must be compatible with the velocity of the particles of the boundary surface (28) at the unperturbed surface . This condition reads
Using equations (29), (48) and (51) for the condition (62), we obtain
(iii) Magnetodynamic condition.
The jump of the normal component of the magnetic field vanishes across the fluid perturbed interface at . This means that
from which we obtain
Therefore, upon using equations (47), (48), (51), (53), and (55) for (65), we get
5 Dispersion relation
Here we apply a compatibility condition known as the compatibility dynamical condition.
The normal component of the velocity vector u must be compatible with the velocity of the particles of the boundary surface (24) at the unperturbed surface .
Mathematically, this condition could be given as
This may be rewritten, on using equation (46), in the form
By substituting from equations (2), (25), (29), (43), (51)(55), (63), (64), and (66) into the condition (68), the following dispersion relation is obtained:
6 Limiting cases
The relation (69) is the desired stability criterion of a streaming fluid cylinder under the combined effects of the capillary, inertia, selfgravitating, and magnetic forces. It is a linear combination of the dispersion relations of a streaming fluid cylinder under the influence of the selfgravitating force only, fluid cylinder under the effects of the capillary force only and the one under the electromagnetic force only.
It contains the natural quantity as well as together with , each as a unit of time. In reality the latter quantities are very interesting and have very important task as we intend to rewrite the relation (69) in a dimensionless form because σ has a unit of (time)^{−1}. This situation is exactly the same as the following cases of Chandrasekhar [5] which were performed for axisymmetric () perturbation of nonstreaming fluid cylinder:
The relation (69) relates the temporal amplification σ with the longitudinal wave number x; the modified Bessel functions and of the first and second kind of order m and with their derivatives, the magnetic field parameter α, the selfgravitating constant G, the basic magnetic field intensity , the fluid density ρ, the radius of the cylinder and with the coefficient μ of the magnetic permeability.
Since the stability criterion (69) is a general relation, we may obtain several published works as limiting cases from it.
Some approximations (, , , , and ) are required for equation (69) to yield
which is the same dispersion relation as that derived by Chandrasekhar and Fermi [16]. In fact, the authors [16] used a totally different method compared to the one used here. They used the method of representing solenoidal vectors in terms of poloidal and toroidal quantities.
If we suppose that (, , , , and ), the relation (69) yields
This relation coincides with that derived regarding the capillary instability of a full liquid jet in a vacuum by Rayleigh [1].
If we suppose that (, , , and ), the relation (69) reduces to
from which we obtain
where use has been made of the Wronskian
for . The relation (73) was established by Chandrasekhar [5] for axisymmetric disturbances.
7 Stability discussions
7.1 Capillary instability
In the absence of the magnetic field, we assume that the streaming fluid is acted upon only by the capillary force. In such a case, the dispersion relation of this model is given from the relation (69) in the form
By using the fact, for each nonzero real value of x and , that
the analytical and numerical discussions of the relation (76) reveal the following results.
In the computer for different values of M and different cases of and .
In the most important sausage mode .
The dimensionless dispersion relation is
where
The numerical data associated with correspond to the unstable states, while those associated with correspond to the stable domains. It has been found that there are many features of interest in this numerical analysis as we see in the following.
(i) For , see Figure 2.
Figure 2. Stable and unstable domains for,.
Corresponding to . It has been found that the unstable domains are
while the neighboring stable domains are
where the equalities correspond to the marginal stability states.
(ii) For , see Figure 3.
Figure 3. Stable and unstable domains for,.
Corresponding to . It has been found that the unstable domains are
while the neighboring stable domains are
where the equalities correspond to the marginal stability states.
(iii) For , see Figure 4.
Figure 4. Stable domains for,.
Corresponding to . It has been found that stable domains are .
(iv) For , see Figure 5.
Figure 5. Stable domains for,.
Corresponding to . It has been found that stable domains are .
We conclude that the streaming full fluid cylinder has stable and unstable domain for M less than 3.5 and stable domain only for M greater than this value whatever the values of velocities are. Increasing the value of M, the unstable domain is decreasing. The effect of changing velocities cases on the capillarity effect is such small that it may be considered as no effect.
7.2 Selfgravitating instability
Consider only the selfgravitating force effect, and then the dispersion relation of the model is given from equation (69) as follows:
Consider the inequalities (77) and (78) and, for each nonzero real value of x, that
the analytical and numerical discussion of the relation (79) reveal the following.
For , , it has been found that the model is gravitationally unstable in the domain ( for mode) while it is stable in the domains ( for mode) and ( for modes).
For , , it has been found that the axial flow has a strong destabilizing influence. That effect does not rely on the kind of perturbation and it is so for all short and long wavelengths. Therefore, the streaming has the effect of increasing the axisymmetric stable domain and the nonaxisymmetric domains .
We conclude that the streaming selfgravitating fluid cylinder is unstable not only for the axisymmetric mode , but also for nonaxisymmetric modes .
7.3 Magnetogravitodynamic stability
This is the case in which the streaming fluid cylinder is acted upon by the combined effects of the selfgravitating and magnetic forces. It is difficult to determine exactly in analytical ways the (un) stable domains in such a general case. However, we could determine them via the numerical discussions. Also, by means of such discussion, we may find out the effects of the magnetic field on the selfgravitating force. This could be carried out by calculating the dimensionless dispersion relation
in the computer for different values of
in the most important sausage mode .
The numerical data associated with correspond to the unstable states, while those associated with correspond to the stable domains. It has been found that there are many features of interest in this numerical analysis as we see in the following.
(i) For , see Figure 6.
Figure 6. Stable and unstable domains for,.
Corresponding to . It has been found that the unstable domains are
while the neighboring stable domains are
where the equalities correspond to the marginal stability states.
(ii) For , see Figure 7.
Figure 7. Stable and unstable domains for,.
Corresponding to . It has been found that the unstable domains are
while the neighboring stable domains are
where the equalities correspond to the marginal stability states.
(iii) For , see Figure 8.
Figure 8. Stable domains for,.
Corresponding to . It has been found that stable domains are .
(iv) For , see Figure 9.
Figure 9. Stable domains for,.
Corresponding to . It has been found that stable domains are .
We conclude that the streaming full fluid cylinder has stable and unstable domain for γ less than 0.8 and stable domain only. Increasing the value of magnetic field, the unstable domains are decreasing. The effect of changing velocities cases on magnetic effect is such small that it may be considered asno effect. If we compare these results with those of chapter two (only velocity in z direction), we observe that the existance of another velocity W in φ direction decreases the unstable domain.
7.4 Magnetogravitodynamic capillary stability
This is the general case in which the streaming fluid cylinder is acted upon by the combined effects of the selfgravitating, capillary, and magnetic forces. The dispersion relation is given in its general form by equation (69). It is difficult to determine exactly in analytical ways the (un) stable domains in such a general case. However, we could determine them via the numerical discussions. Also, by means of such discussion, we may find out the effects of capillary with a constant magnetic field on the selfgravitating force. This could be carried out by calculating the dimensionless dispersion relation
in the computer for different values of
in the most important sausage mode .
The numerical data associated with correspond to the unstable states, while those associated with correspond to the stable domains. It has been found that there are many features of interest in this numerical analysis as we see in the following.
(i) For , see Figure 10.
Figure 10. Stable and unstable domains for,.
Corresponding to and . It has been found that the unstable domains are
while the neighboring stable domains are
where the equalities correspond to the marginal stability states.
(ii) For , see Figure 11.
Figure 11. Stable and unstable domains for,.
Corresponding to and . It has been found that the unstable domains are
while the neighboring stable domains are
where the equalities correspond to the marginal stability states.
(iii) For , see Figure 12.
Figure 12. Stable domains for,.
Corresponding to and . It has been found that stable domains are
(iv) For , see Figure 13.
Figure 13. Stable domains for,.
Corresponding to and . It has been found that stable domains are
We conclude that the streaming full fluid cylinder has stable and unstable domains for M less than 3.5 and stable domain only for M greater than this value whatever the values of velocities are. The effect of changing velocities cases on capillarity effect is such small that it may be considered as no effect. Increasing M with constant magnetic field increases the unstable domain.
8 Conclusion
From the foregoing numerical results, we may deduce the following:
(1) The velocity has a strong destabilizing influence on the selfgravitating instability of the model.
(2) The capillary force has a strong stabilizing influence on the selfgravitating instability of the model.
(3) The capillary and selfgravitating modified a lot the instability of the model for all short and long wavelengths.
(4) The velocity has a strong destabilizing influence on the selfgravitating instability of the model.
(5) The magnetic force has a strong stabilizing influence on the selfgravitating instability of the model.
(6) The selfgravitating instability character has disappeared and has been dispersed, and the model has become completely stable.
(7) The velocities in two directions have a strong destabilizing influence on the selfgravitating instability of the model.
(8) The magnetic force has a strong stabilizing influence on the selfgravitating capillary instability of the model.
Competing interests
The authors declare that they have no competing interests.
Authors’ contributions
All authors read and approved the final manuscript.
Acknowledgements
We are grateful to the editor of the journal and the reviewers for their suggestion and comments of this paper.
References

Yuen, MC: Nonlinear capillary instability of a liquid jet. J. Fluid Mech.. 33, 151 (1968). Publisher Full Text

Nayfeh, A, Hassan, SD: The method of multiple scales and nonlinear dispersive waves. J. Fluid Mech.. 48, 463 (1971). Publisher Full Text

Kakutani, T, Inoue, I, Kan, T: Nonlinear capillary waves on the surface of liquid column. J. Phys. Soc. Jpn.. 37, 529 (1974). Publisher Full Text

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

Radwan, AE, Aly, FA: Selfgravitating instability of two semiinfinite streaming superposed fluids endowed with surface tension. Nuovo Cimento. 113, 601 (1998)

Radwan, AE, Ali, RM: Magnetohydrodynamic instability of a dissipative compressible rotating selfgravitating fluid medium. Nuovo Cimento. 114B, 1361 (1999)

Radwan, AE: Periodic time dependent electrogravitational instability of a fluid cylinder. Phys. Scr.. 67, 510 (2007)

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

Radwan, AE, Hasan, AA: Magnetohydrodynamic stability of selfgravitational fluid cylinder. Appl. Math. Model.. 33, 2121 (2009). Publisher Full Text

Hasan, AA: Electrogravitational stability of oscillating streaming fluid cylinder. Physica B. 406(2), 234 (2011). Publisher Full Text

Hasan, AA: Capillary electrodynamic stability of selfgravitational fluid cylinder with varying electric field. J. Appl. Mech.. 79(2), 1 (2011)

Hasan, AA, Mekheimer, KS, Azwaz, SL: Hydromagnetic stability of selfgravitational oscillating streaming fluid jet pervaded by azimuthal varying magnetic field. Int. J. Math. Arch.. 2(4), 488 (2011)

Hasan, AA: Electrogravitational stability of oscillating streaming dielectric compound jets ambient with a transverse varying electric field. Bound. Value Probl.. 2011, Article ID 31 (2011)

Hasan, AA: Hydromagnetic instability of streaming jet pervaded internally by varying transverse magnetic field. Math. Probl. Eng.. 2012, Article ID 325423 (2012)

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