Skip to main content

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

Abstract

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.

Introduction

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

1 Formulation of the problem

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

u ̲ 0 =(0,W,U)
(1)

and pervaded internally and externally by the magnetic fields

H ̲ 0 =(0,0, H 0 ), H ̲ 0 e x =(0,0,α H 0 ).
(2)

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.

Figure 1
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

(3)
(4)
(5)
(6)
(7)

The curvature pressure due to the capillary force is

P s =T( N ̲ s )
(8)

with

N ̲ s =F/|F|,
(9)

where

F(r,φ,z)=0
(10)

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

(11)
(12)
(13)

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.

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

ρ V 0 P 0 + ( μ 4 π ) ( H ̲ 0 ) H ̲ 0 ( μ 8 π ) ( H ̲ 0 H ̲ 0 )=0,
(14)

from which, taking into account equation (5), we obtain (ρ V 0 P 0 ( μ 8 π ) H 0 2 )=0.

By integrating this equation, we get

P 0 =ρ V 0 ( μ 8 π ) H 0 2 +C,
(15)

where C is a constant of integration to be determined.

The surface pressure due to the capillary force (cf. Chandrasekhar [5]) is given by

P 0 s =T/ R 0 .
(16)

The self-gravitating potentials V 0 and V 0 e x of the unperturbed state satisfy

(17)
(18)

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

(19)
(20)

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

(21)
(22)

Therefore,

(23)
(24)

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

P 0 = ( T R 0 ) +πG ρ 2 ( R 0 2 r 2 ) + ( μ 8 π ) H 0 2 ( α 2 1 ) .
(25)

It is worth noting that in the absence of surface tension at the boundary surface

α0,
(26)

in order that

P0.
(27)

3 Perturbation analysis

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

r= R 0 +ε(t) R 1 +
(28)

with

R 1 =exp ( i ( k z + m φ ) ) .
(29)

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

ε(t)= ε 0 exp(σt),
(30)

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

Q(r,φ,z,t)= Q 0 (r)+ Q 1 (r,φ,z,t).
(31)

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

(32)
(33)
(34)
(35)
(36)
(37)
(38)
(39)
(40)

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

Q 1 (r,φ,z;t)= ε 0 q 1 (r)exp ( σ t + i ( k z + m φ ) ) .
(41)

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

( 1 r ) ( d d r ) ( r ( d ϕ 1 ( r ) d r ) ) ( ( m 2 r 2 ) + k 2 ) ϕ 1 (r)=0,
(42)

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

(43)
(44)

Here I m (kr) and K m (kr) 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

(σ+imW+ikU) u 1 ̲ (iμk/4πρ) H 0 H 1 ̲ = Π 1
(45)

with

Π 1 = ( P 1 ρ ) V 1 +(μ/4πρ)( H ̲ 0 H ̲ 1 ).
(46)

Also, equation (35) yields

H ̲ 1 = ( i k H 0 σ + i m W + i k U ) u 1 ̲ .
(47)

By combining equations (45) and (47), we get

u 1 ̲ = ( ( σ + i m W + i k U ) ( ( σ + i m W + i k U ) 2 + Ω A 2 ) ) Π 1 ,
(48)

where

Ω A = ( μ k 2 H 0 2 4 π ρ ) 1 / 2
(49)

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

2 Π 1 =0.
(50)

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

Π 1 = C 4 ε 0 I m (kr)exp ( σ t + i ( k z + m φ ) ) ,
(51)

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

P 1 s = ( T R 0 2 ) ( 1 m 2 x 2 ) exp ( σ t + i ( k z + m φ ) ) ,
(52)

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

H 1 e x = ψ 1 e x .
(53)

By combining equations (38) and (53), we get

2 ψ 1 e x =0.
(54)

Similarly, as it has been done for equation (50), equation (54) is solved and its finite solution is given by

ψ 1 e x = C 5 ε 0 K m (kr)exp ( σ t + i ( k z + m φ ) ) ,
(55)

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

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

  1. (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

(56)
(57)

By substituting from equations (23), (24), (29), (43) and (44) into the conditions (56) and (57), we get

(58)
(59)

from which we obtain

(60)
(61)

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

  1. (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

u ̲ 1 r = R 1 t +U R 1 z +W R 1 φ .
(62)

Using equations (29), (48) and (51) for the condition (62), we obtain

C 4 = ( ( σ + i m W + i k U ) 2 + Ω A 2 ) ( R 0 / x I m ( x ) ) .
(63)
  1. (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

N s ̲ H ̲ N s ̲ H ̲ e x =0at r= R 0 ,
(64)

from which we obtain

H ̲ 1 r H ̲ 1 r e x =ik R 1 H 0 (1α).
(65)

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

C 5 = i α H 0 K m ( x ) .
(66)

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

Mathematically, this condition could be given as

P 1 + R 1 P 0 r + μ 4 π ( H ̲ 0 H ̲ 1 ) μ 4 π ( H ̲ 0 H ̲ 1 ) e x = P 1 s .
(67)

This may be rewritten, on using equation (46), in the form

ρ( Π 1 + V 1 )= P 1 s + μ 4 π ( H ̲ 0 H ̲ 1 ) e x R 1 P 0 r .
(68)

By substituting from equations (2), (25), (29), (43), (51)-(55), (63), (64), and (66) into the condition (68), the following dispersion relation is obtained:

( σ + i m W + i k U ) 2 = 4 π G ρ x I m ( x ) I m ( x ) [ I m ( x ) K m ( x ) 1 2 ] + T ρ R 0 3 ( 1 m 2 x 2 ) x I m ( x ) I m ( x ) + μ H 0 2 ( 4 π ρ R 0 2 ) [ x 2 + α 2 x 2 I m ( x ) K m ( x ) I m ( x ) K m ( x ) ] .
(69)

6 Limiting cases

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)−1. 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:

(T=0,G0, H 0 0),(T0,G=0, H 0 0)and(T0,G0, H 0 =0).

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

σ 2 =4πGρ x I 1 ( x ) I 0 ( x ) [ I 0 ( x ) K 0 ( x ) 1 2 ] , I 0 (x)= I 1 (x),
(70)

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

σ 2 = T ρ R 0 3 x I 1 ( x ) I 0 ( x ) ( 1 x 2 ) .
(71)

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

σ 2 =4πGρ x I 1 ( x ) I 0 ( x ) [ I 0 ( x ) K 0 ( x ) 1 2 ] + μ H 0 2 σ R 0 2 [ x 2 + x 2 I 1 ( x ) K 0 ( x ) I 0 ( x ) ( K 1 ( x ) ) ] ,
(72)

from which we obtain

σ 2 =4πGρ x I 1 ( x ) I 0 ( x ) [ I 0 ( x ) K 0 ( x ) 1 2 ] + μ H 0 2 ρ R 0 2 x I 0 ( x ) K 1 ( x ) ,
(73)

where use has been made of the Wronskian

I m (x) K m (x) I m (x) K m (x)= x 1
(74)

for m=0. 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

( σ + i m W + i k U ) 2 = T ρ R 0 3 x I 1 ( x ) I 0 ( x ) ( 1 m 2 x 2 ) .
(75)

By using the fact, for each non-zero real value of x and m0, that

(76)
(77)

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

( σ + i m W + i k U ) 2 4 π G ρ = x I m ( x ) I m ( x ) [ I m ( x ) K m ( x ) 1 2 ] +M [ ( 1 m 2 x 2 ) x I m ( x ) I m ( x ) ] ,
(78)

where

M= [ T ( 4 π G ρ 2 R 0 3 ) ] , U = [ i k U ( 4 π G ρ ) 1 / 2 ] , W = [ i m W ( 4 π G ρ ) 1 / 2 ] .

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.

  1. (i)

    For M=0.5,1.0,1.5,2.5,3.0,3.5, see Figure 2.

Figure 2
figure 2

Stable and unstable domains for γ=0 , U = W =0.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.

  1. (ii)

    For M=0.5,1.0,1.5,2.5,3.0,3.5, see Figure 3.

Figure 3
figure 3

Stable and unstable domains for γ=0 , U = W =0.5 .

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.

  1. (iii)

    For M=3.6,4.0,4.5,5.0,5.5,6.0, see Figure 4.

Figure 4
figure 4

Stable domains for γ=0 , U = W =0.2 .

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

  1. (iv)

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

Figure 5
figure 5

Stable domains for γ=0 , U = W =0.5 .

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

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:

( σ + i m W + i k U ) 2 =4πGρ x I m ( x ) I m ( x ) [ I m ( x ) K m ( x ) 1 2 ] .
(79)

Consider the inequalities (77) and (78) and, for each non-zero real value of x, that

K m (x)>0
(80)

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.0667x for m=0 mode) and (0x for m1 modes).

For U0, W0, 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.0667x 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 m1.

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

( σ + i m W + i k U ) 2 4 π G ρ = x I m ( x ) I m ( x ) [ I m ( x ) K m ( x ) 1 2 ] +γ [ x 2 + α 2 x 2 I m ( x ) K m ( x ) I m ( x ) K m ( x ) ]
(81)

in the computer for different values of

γ ( = ( μ / 16 π 2 G ) ( H 0 / ρ R 0 ) 2 ) and U = [ i k U ( 4 π G ρ ) 1 / 2 ]

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.

  1. (i)

    For γ=0.1,0.2,0.3,0.4,0.5,0.7, see Figure 6.

Figure 6
figure 6

Stable and unstable domains for M=0 , U = W =0.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.

  1. (ii)

    For γ=0.1,0.2,0.3,0.4,0.5,0.7, see Figure 7.

Figure 7
figure 7

Stable and unstable domains for M=0 , U = W =0.5 .

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.

  1. (iii)

    For γ=0.8,0.9,1.0,1.2,1.3,1.4, see Figure 8.

Figure 8
figure 8

Stable domains for M=0 , U = W =0.2 .

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

  1. (iv)

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

Figure 9
figure 9

Stable domains for M=0 , U = W =0.5 .

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

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

( σ + i m W + i k U ) 2 4 π G ρ = x I m ( x ) I m ( x ) [ I m ( x ) K m ( x ) 1 2 ] + M [ ( 1 m 2 x 2 ) x I m ( x ) I m ( x ) ] + γ [ x 2 + α 2 x 2 I m ( x ) K m ( x ) I m ( x ) K m ( x ) ]
(82)

in the computer for different values of

M= [ T ( 4 π G ρ 2 R 0 3 ) ] ,γ ( = ( μ / 16 π 2 G ) ( H 0 / ρ R 0 ) 2 ) and U = [ i k U ( 4 π G ρ ) 1 / 2 ]

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.

  1. (i)

    For M=0.5,1.0,1.5,2.5,3.0,3.5, see Figure 10.

Figure 10
figure 10

Stable and unstable domains for γ=0.5 , U = W =0.2 .

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.

  1. (ii)

    For M=0.5,1.0,1.5,2.5,3.0,3.5, see Figure 11.

Figure 11
figure 11

Stable and unstable domains for γ=0.5 , U = W =0.5 .

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.

  1. (iii)

    For M=0.5,1.0,1.5,2.5,3.0,3.5, see Figure 12.

Figure 12
figure 12

Stable domains for γ=0.5 , U = W =0.2 .

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

0x<.
  1. (iv)

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

Figure 13
figure 13

Stable domains for γ=0.5 , U = W =0.5 .

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

0x<.

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. (1)

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

  2. (2)

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

  3. (3)

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

  4. (4)

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

  5. (5)

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

  6. (6)

    The self-gravitating instability character has disappeared and has been dispersed, and the model has become completely stable.

  7. (7)

    The velocities in two directions have a strong destabilizing influence on the self-gravitating instability of the model.

  8. (8)

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

References

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

    MATH  Google Scholar 

  2. Yuen MC: Non-linear capillary instability of a liquid jet. J. Fluid Mech. 1968, 33: 151. 10.1017/S0022112068002429

    Article  MATH  Google Scholar 

  3. Nayfeh A, Hassan SD: The method of multiple scales and non-linear dispersive waves. J. Fluid Mech. 1971, 48: 463. 10.1017/S0022112071001708

    Article  MATH  MathSciNet  Google Scholar 

  4. Kakutani T, Inoue I, Kan T: Nonlinear capillary waves on the surface of liquid column. J. Phys. Soc. Jpn. 1974, 37: 529. 10.1143/JPSJ.37.529

    Article  Google Scholar 

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

    Google Scholar 

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

    Google Scholar 

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

    Google Scholar 

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

    Article  Google Scholar 

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

    MathSciNet  Google Scholar 

  10. Radwan AE, Hasan AA: Magnetohydrodynamic stability of self-gravitational fluid cylinder. Appl. Math. Model. 2009, 33: 2121. 10.1016/j.apm.2008.05.014

    Article  MATH  MathSciNet  Google Scholar 

  11. Hasan AA: Electrogravitational stability of oscillating streaming fluid cylinder. Physica B 2011, 406(2):234. 10.1016/j.physb.2010.10.050

    Article  Google Scholar 

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

    Google Scholar 

  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. 2011, 2(4):488.

    Google Scholar 

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

    Google Scholar 

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

    Google Scholar 

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

    Article  MathSciNet  Google Scholar 

Download references

Acknowledgements

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

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Alfaisal A Hasan.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 2.0 International License (https://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Hasan, A.A., Abdelkhalek, R.A. Magnetogravitodynamic stability of streaming fluid cylinder under the effect of capillary force. Bound Value Probl 2013, 48 (2013). https://doi.org/10.1186/1687-2770-2013-48

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1687-2770-2013-48

Keywords