Magnetogravitodynamic stability of streaming fluid cylinder under the effect of capillary force

The magnetohydrodynamic stability criterion of self-gravitating streaming fluid cylinder under the combined effect of self-gravitating, 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 self-gravitating 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 self-gravitating instability of the model have been examined. Including the magnetic force together with self-gravitating force improves the instability of the model. However, the self-gravitating instability will never be suppressed whatever the effects of the MHD force stabilizing effects are.

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.[6-10] extended such interesting works by studying the magnetohydrodynamic stability of a liquid jet embedded into a tenuous medium for all axisymmetric and non-axisymmetric modes of perturbation. 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 [9] and [10]. They [9] studied the gravitational stability of a fluid cylinder under transverse time-dependent electric field for axisymmetric perturbations. Hasan [11] discussed the stability of oscillating streaming fluid cylinder subject to the combined effect of the capillary, self-gravitating, and electrodynamic forces for all axisymmetric and non-axisymmetric perturbation modes. He [12] 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 modes of perturbations. In [13] Hasan et al. investigated the hydromagntic stability of a self-gravitational oscillating streaming fluid jet pervaded by azimuthal varying magnetic field for all axisymmetric and non-axisymmetric modes of perturbation. He [14] discussed 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 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 self-gravitating instability of the present models. This may be carried out, for all axisymmetric and non-axisymmetric modes of perturbation, analytically and the results will be verified numerically.

We consider a uniform cylinder of an incompressible inviscid fluid of radius R 0 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, H 0 is the intensity of the magnetic field in the fluid, and α is some parameter. The components of u 0 , H 0 , H 0 e x are considered along the cylindrical coordinates ( r , φ , z ) with the z-axis 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 self-gravitating, inertial, capillary, and magnetic forces. The surrounding tenuous medium of the fluid cylinder is acted upon by the self-gravitating and magnetic forces only.

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 N ̲ s is a unit outward vector normal to the surface, T is surface tension, and P s is pressure due to curvature.

For the surrounding tenuous medium, the basic equations are

Here ρ, u 0 , and P are the fluid mass density, velocity vector, and kinetic pressure, respectively; H 0 , H 0 e x are the magnetic field intensities and V, V e x are self-gravitating potentials, respectively, inside and outside the fluid cylinder, μ is the magnetic field permeability coefficient and G is the gravitational constant.

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 ∇ ( ρ V 0 − P 0 − ( μ 8 π ) H 0 2 ) = 0 .

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 self-gravitating potentials V 0 and V 0 e x of the unperturbed state satisfy

The non-singular solutions of equations (17) and (18) in the cylindrical coordinates ( r , φ , z ) with cylindrical symmetries ( ∂ ∂ φ ) = 0 and ( ∂ ∂ z ) = 0 are given by

where C 1 , C 2 , and C 3 are constants of integration to be determined. By applying the conditions that the self-gravitational potential V and its derivative must be continuous across the unperturbed boundary surface at r = R 0 and choosing C 1 = 0 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 r = R 0 , 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

We consider small departures from an unperturbed right-cylindrical 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 R 1 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 ε ( t ) of the perturbation is given by

where ε 0 (=ε at t = 0 ) is the initial amplitude and σ is the temporal amplification. If σ ( = i ω , i = − 1 ) is imaginary, then ω / 2 π is the oscillation frequency of the propagating wave in the fluid.

As the initial streaming state is perturbed, every physical quantity Q ( r , φ , z ; t ) may be expanded as

Here Q stands for P, u, V, V e x , H, H e x , and N s while Q 0 indicates the unperturbed quantity and Q 1 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 Q 1 ( r , φ , z ; t ) may be expressed as

By means of the expansion (41), equations (36) and (40) give the second-order ordinary differential equation

where ϕ 1 ( r ) stands for V 1 ( r ) and V 1 e x ( r ) . 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 I m ( k r ) and K m ( k r ) 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 space-time 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 H 0 .

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 non-singular solution for Π 1 ( r , φ , z ; t ) is given by

where C 4 is a constant of integration to be determined.

The pressure surface P 1 s in the perturbed state due to the capillary force is determined from equation (37) along with (29) in the form

where x ( = k R 0 ) is the dimensionless longitudinal wavenumber.

Now, equation (34) means that the magnetic field intensity H 1 e x in the perturbed state may be derived from a scalar function, ψ 1 e x 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 C 5 is a constant of integration to be determined upon applying 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 r = R 0 .

Under the present circumstances, these boundary conditions may be stated as follows.

(i) Self-gravitating conditions.

The gravitational potential and its derivative must be continuous across the perturbed fluid interface (28) at the unperturbed boundary r = R 0 . These conditions at r = R 0 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 ( = k R 0 ) 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 r = R 0 . 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 r = R 0 . This means that

from which we obtain

Therefore, upon using equations (47), (48), (51), (53), and (55) for (65), we get

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 r = R 0 .

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:

The relation (69) is the desired stability criterion of a streaming fluid cylinder under the combined effects of the capillary, inertia, self-gravitating, and magnetic forces. It is a linear combination of the dispersion relations of a streaming fluid cylinder under the influence of the self-gravitating 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 ( T / ρ R 0 3 ) − 1 2 as well as ( μ H 0 2 / 4 π ρ R 0 2 ) − 1 2 together with ( 4 π G ρ ) − 1 2 , 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)⁻¹. This situation is exactly the same as the following cases of Chandrasekhar [5] which were performed for axisymmetric ( m = 0 ) perturbation of nonstreaming fluid cylinder:

The relation (69) relates the temporal amplification σ with the longitudinal wave number x; the modified Bessel functions I m ( x ) and K m ( x ) of the first and second kind of order m and with their derivatives, the magnetic field parameter α, the self-gravitating constant G, the basic magnetic field intensity H 0 , the fluid density ρ, the radius R 0 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 ( α = 0 , H 0 = 0 , U = 0 , W = 0 , T = 0 and m = 0 ) 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 ( α = 0 , H 0 = 0 , U = 0 , W = 0 , G = 0 and m = 0 ), 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 ( α = 1 , U = 0 , W = 0 , T = 0 and m = 0 ), the relation (69) reduces to

from which we obtain

where use has been made of the Wronskian

for m = 0 . The relation (73) was established by Chandrasekhar [5] for axisymmetric disturbances.

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 non-zero real value of x and m ≥ 0 , 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 U ∗ and W ∗ .

In the most important sausage mode m = 0 .

The dimensionless dispersion relation is

where

The numerical data associated with σ / ( 4 π G ρ ) 1 2 correspond to the unstable states, while those associated with ω / ( 4 π G ρ ) 1 2 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 M = 0.5 , 1.0 , 1.5 , 2.5 , 3.0 , 3.5 , see Figure 2.

Corresponding to U ∗ = W ∗ = 0.2 . 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 M = 0.5 , 1.0 , 1.5 , 2.5 , 3.0 , 3.5 , see Figure 3.

Corresponding to U ∗ = W ∗ = 0.5 . 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 M = 3.6 , 4.0 , 4.5 , 5.0 , 5.5 , 6.0 , see Figure 4.

Corresponding to U ∗ = W ∗ = 0.2 . It has been found that stable domains are 0 ≤ x < ∞ .

(iv) For M = 3.6 , 4.0 , 4.5 , 5.0 , 5.5 , 6.0 , see Figure 5.

Corresponding to U ∗ = W ∗ = 0.5 . It has been found that stable domains are 0 ≤ x < ∞ .

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 Self-gravitating instability

Consider only the self-gravitating 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 non-zero real value of x, that

the analytical and numerical discussion of the relation (79) reveal the following.

For U = 0 , W = 0 , it has been found that the model is gravitationally unstable in the domain ( 0 < x < 1.0667 for m = 0 mode) while it is stable in the domains ( 1.0667 ≤ x ≤ ∞ for m = 0 mode) and ( 0 ≤ x ≤ ∞ for m ≥ 1 modes).

For U ≠ 0 , W ≠ 0 , 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 1.0667 ≤ x ≤ ∞ and the non-axisymmetric domains 0 < x < ∞ .

We conclude that the streaming self-gravitating fluid cylinder is unstable not only for the axisymmetric mode m = 0 , but also for non-axisymmetric modes m ≥ 1 .

7.3 Magnetogravitodynamic stability

This is the case in which the streaming fluid cylinder is acted upon by the combined effects of the self-gravitating 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 self-gravitating 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 m = 0 .

The numerical data associated with σ / ( 4 π G ρ ) 1 2 correspond to the unstable states, while those associated with ω / ( 4 π G ρ ) 1 2 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 γ = 0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.7 , see Figure 6.

Corresponding to U ∗ = W ∗ = 0.2 . 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 γ = 0.1 , 0.2 , 0.3 , 0.4 , 0.5 , 0.7 , see Figure 7.

Corresponding to U ∗ = W ∗ = 0.5 . 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 γ = 0.8 , 0.9 , 1.0 , 1.2 , 1.3 , 1.4 , see Figure 8.

Corresponding to U ∗ = W ∗ = 0.2 . It has been found that stable domains are 0 ≤ x < ∞ .

(iv) For γ = 0.8 , 0.9 , 1.0 , 1.2 , 1.3 , 1.4 , see Figure 9.

Corresponding to U ∗ = W ∗ = 0.5 . It has been found that stable domains are 0 ≤ x < ∞ .

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 self-gravitating, 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 self-gravitating 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 m = 0 .

The numerical data associated with σ / ( 4 π G ρ ) 1 2 correspond to the unstable states, while those associated with ω / ( 4 π G ρ ) 1 2 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 M = 0.5 , 1.0 , 1.5 , 2.5 , 3.0 , 3.5 , see Figure 10.

Corresponding to γ = 0.5 and U ∗ = W ∗ = 0.2 . 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 M = 0.5 , 1.0 , 1.5 , 2.5 , 3.0 , 3.5 , see Figure 11.

Corresponding to γ = 0.5 and U ∗ = W ∗ = 0.5 . 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 M = 0.5 , 1.0 , 1.5 , 2.5 , 3.0 , 3.5 , see Figure 12.

Corresponding to γ = 0.5 and U ∗ = W ∗ = 0.2 . It has been found that stable domains are

(iv) For M = 0.5 , 1.0 , 1.5 , 2.5 , 3.0 , 3.5 , see Figure 13.

Corresponding to γ = 0.5 and U ∗ = W ∗ = 0.5 . 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.

From the foregoing numerical results, we may deduce the following:

(1) The velocity has a strong destabilizing influence on the self-gravitating instability of the model.

(2) The capillary force has a strong stabilizing influence on the self-gravitating instability of the model.

(3) The capillary and self-gravitating modified a lot the instability of the model for all short and long wavelengths.

(4) The velocity has a strong destabilizing influence on the self-gravitating instability of the model.

(5) The magnetic force has a strong stabilizing influence on the self-gravitating instability of the model.

(6) The self-gravitating 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 self-gravitating instability of the model.

(8) The magnetic force has a strong stabilizing influence on the self-gravitating capillary instability of the model.

The authors declare that they have no competing interests.

All authors read and approved the final manuscript.

We are grateful to the editor of the journal and the reviewers for their suggestion and comments of this paper.

1. Rayleigh, JM: The Theory of Sound, Dover, New York (1945)

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

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

4. 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

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

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

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

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

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

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

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

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

13. 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)

14. 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)

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

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