# Stability of Compressible Hollow Jet Pervaded by a Transverse Varying Magnetic Field

Authors: Samia S. Elazab; Samy A. Rahman; Alfaisal A. Hasan; Nehad A. Zidan
Journal: ISRN Mathematical Physics
 (2012)

---

## Abstract

The magnetohydrodynamic stability of an ordinary compressible hollow cylinder pervaded by a transverse varying magnetic field, under the influence of capillary, inertia, and Lorentz force, has been developed. The problem is modelized. The basic equations formulated, solved, and, upon applying appropriate boundary conditions, the singular solutions are excluded. The eigenvalue relation has been derived and discussed. The capillary force has destabilizing influence only for long wavelengths in the axisymmetric perturbation but it is stabilizing in the rest and also so in the nonaxisymmetric perturbations. The compressibility increases the stable domains and simultaneously decreases those of instability. The electromagnetic force has different effects due to the axial uniform field and varying transverse one. The axial field is stabilizing for all wavelengths in all kinds of perturbations. The transverse field is stabilizing or not according to restrictions. Here, the high compressibility increases rapidly the magnetodynamic stable domains and leads to shrinking those of instability.


---

## Body

## 1. Introduction

The stability of a fluid cylinder subjected to different external forces has been published and reported by Rayleigh [1] and Chandrasekhar [2]. The instability of the mirror case of a gas cylinder surrounded by an ideal fluid endowed with surface tension is envised by Chandrasekhar [2] for axisymmetric perturbation. See also Drazin and Reid [3] and Cheng [4]. Kendall [5] made very neat experiments to check the breaking-up of this model due to its instability as the inertia force of the exterior liquid is predominant over that of the gas cylinder. Moreover, he did write about its applications in several domains of science and attracted the attention of researchers for elaborating the stability of this model. One has to infer here that the analytical results given by Cheng [4], in (2.4) and (2.5), are incorrect in the third term. In fact, the term (1-s2-k2R02) must be in the numerator as it is clear from (2.3). See also (5.1)–(6.1) in the present work and Drazin and Reid’s result [3,  page 16], and also Chandrasekhar’s dispersion relation [2, page 538 and page 540], ((147) and (155) there). For other works extending some of the previous results, we may refer to Radwan [6–9]. In all these foregoing studies, the fluids are considered to be incompressible and nonviscous. See also Chen and Lin [10].In reality, in the model of a gas jet injected into a liquid, we have to consider the fluid is compressible. In [11, 12], Radwan et al. tried to investigate the instability of that model pervaded by constant magnetic field.The stability of different cylindrical models under the action of selfgravitating force in addition to other forces has been elaborated by Radwan and Hasan [13, 14]. Hasan [15] has discussed the stability of oscillating streaming fluid cylinder subjected to combined effect of the capillary, selfgravitating, and electrodynamic forces for all axisymmetric and nonaxisymmetric perturbation modes.There are many applications of magnetohydrodynamic stability in several fields of science such as the following.(i)
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 due to the motion of the liquid iron.(ii)
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, and so forth.  Many astrophysical systems are not in local thermal equilibrium and, therefore, require an additional kinematic treatment to describe all phenomena within the system.(iii)
Engineering applications: there are many forms in engineering sciences include oil and gas extraction process if it is surrounded by electric field or magnetic field, gas and steam turbines, MHD power generation systems, magneto-flow meters, and so forth.Here, we present a complete analysis of the stability of compressible hollow cylinder pervaded by a transverse varying magnetic field for all axisymmetric and nonaxisymmetric modes of perturbation where the fluid velocity is not solenoidal any more in the present analysis.


## 2. Formulation of the Problem

Consider a gas cylinder of radiusR0 surrounded by a non-viscous, compressible and perfectly conducting liquid. The liquid is pervaded by the uniform magnetic field H̲0=(0,0,H0) while the gas cylinder is penetrated by the transverse varying magnetic field H̲0g=(0,βH0r/R0,0), where β is a parameter satisfying certain restrictions and H0 is the intensity of the magnetic field in the liquid region, see Figure 1. We shall use the cylindrical coordinates (r,φ,z) with the z-axis coinciding with the axis of the gas cylinder; the components of H̲0 and H̲0g are considered along the coordinates (r,φ,z). The model is acting upon the electromagnetic, capillary, inertia and pressure gradient forces such that the liquid inertia force is paramount over that of the gas cylinder. The fundamental equations for such a study are the combination of the ordinary fluid dynamics, Maxwell electromagnetic equation and those of the perfect gas. Under the present circumstances these equations are given as follows.Figure 1
Sketch for MHD hollow jet.Equations of motion:(2.1)∂ur∂t+(u̲⋅∇)ur-μρ(H̲⋅∇)Hr=-∂Π∂r,(2.2)∂uφ∂t+(u̲⋅∇)uφ-μρ(H̲⋅∇)Hφ=-1r∂Π∂φ,(2.3)∂uz∂t+(u̲⋅∇)uz-μρ(H̲⋅∇)Hz=-∂Π∂z,
where ρΠ is the total magnetohydrodynamic pressure which is the sum of kinetic and magnetic pressures, given by(2.4)ρΠ=P+(μ2)(H̲⋅H̲).
Equation of conservation of flux:(2.5)∇⋅H̲=0.
Equation of conservation of mass of a compressible fluid:(2.6)∂ρ∂t+(u̲⋅∇)ρ+ρ(∇⋅u̲)=0.
Equation of state gives the relation between the pressure P and density ρ:(2.7)P=ρRcT,Rc is a general constant of gas.Equation of conservation of energy:(2.8)ρCv(∂∂t+u̲⋅∇)T+P(∇⋅u̲)=0.
Evaluation equations of a magnetic field (derived from the electromagnetic Maxwell equation) in a perfectly conducting fluid:(2.9)∂Hr∂t+(u̲⋅∇)Hr+(∇⋅u̲)Hr=(H̲⋅∇)ur,(2.10)∂Hφ∂t+(u̲⋅∇)Hφ+(∇⋅u̲)Hφ=(H̲⋅∇)uφ,∂Hz∂t+(u̲⋅∇)Hz+(∇⋅u̲)Hz=(H̲⋅∇)uz.
In the gas cylinder region, surrounded by a liquid, there is no current flow:(2.11)∇⋅H̲g=0,(2.12)∇ΛH̲g=0.
Along the gas liquid interface, the curvature pressure Ps due to the capillary force is given by(2.13)Ps=S(1r1+1r2),
with(2.14)r1-1+r2-1=∇⋅N̲,
where S is the surface tension coeffient, while r1 and r2 are the principle radii of curvature of the gas fluid and N̲ the unit outward drawn normal vector to the perturbed interface, F(r,φ,z,t)=0, given by(2.15)N̲=∇F(r,φ,z,t)|∇F(r,φz,t)|.
Here, (ur,uφ,uz) and (Hr,Hφ,Hz) are, respectively, the components of the velocity vector u̲ of the fluid and the magnetic field intensity H̲, while ρ,P,Cv, and T are the fluid mass density, kinetic pressure, specific heat at a constant volume, and temperature of the fluid, and H̲g is the intensity of the magnetic field in the gas cylinder.The unperturbed state is studied by simplifying the fundamental equation in view ofu̲0=0, ∂/∂φ=0 and ∂/∂z=0 and integrating the resulting differential equations. Finally, by means of the continuity of the normal component of the stress tensor across the gas-fluid interface at r=R0, the kinetic pressure of the fluid in the unperturbed state is obtained, namely,(2.16)P0=P0g+(μ2)(β2-1)H02-SR0.
Equation (2.16) gives the kinetic pressure P0 of the fluid and it is a simple linear combination of the contribution of the different forces effects.(a)
-S/R0 is the contribution due to the capillary force.(b)
(μH02/2)β2 is the contribution due to the magnetodynamic force acting in the gas region.(c)
(μH02/2)(-1) is the contribution due to the magnetodynamic force acting in the fluid region.(d)
P0g is the gas kinetic pressure and it must be suitably strong otherwise the model collapses and the gas spreads into fluid region surrounding the gas cylinder.One has to refer here that the net magnetodynamic force acting on the model has no contribution in the balance of the total pressure (cf. (2.16)) in the unperturbed state in the following cases:(1)
asH0=0,(2)
asH0≠0, β=1.To maintain and keep the model without collapsing in the unperturbed state, it must beP0>0, this means that the contribution of the capillary force S/R0 must be less than the gas pressure P0g, otherwise, the model will collapse and the gas may spread through the fluid region.In the general case, in order thatP0≥0, the gas kinetic pressure P0g in the initial state must satisfy the restriction:(2.17)P0g≥SR0+μH022(1-β2).
Otherwise, the model collapses and will be a homogeneous fluid medium.


## 3. Perturbation Analysis

For small departures from the unperturbed state, due to the perturbation along the gas-fluid interface, every variable quantityQ(r,φ,z,t) may be expressed as its unperturbed part plus a fluctuation part (see [11, 12], viz.)(3.1)Q(r,φ,z,t)=Q0(r)+δ(t)Q1(r,φ,z).Here,Q stands for each of u̲,H̲,ρ,P,T, and the perturbed cross-section radius distance of the gas-fluid interface. The amplitude δ of the perturbation is given by(3.2)δ=δ0exp⁡⁡(σt),
where δ0(=δ at t=0) is the initial amplitude and σ is the temporal amplification at any instant of time t.Consider a sinusoidal wave along the gas-fluid interface, for a single Fourier term, the perturbed cylindrical radial distance is described by(3.3)r=R0+R1,
with(3.4)R1=δ0exp⁡⁡[σt+i(kz+mφ)]
being the elevation of the surface wave measured from the unperturbed position, where k and m are the axial and transverse wave numbers. From the view point of  the foregoing expansions (3.1)–(3.4), the relevant perturbation equations of the present case could be deduced from the fundamental equations (2.1)–(2.15) as follows.For the fluid surrounding the gas jet,(3.5)σu̲1-(μρ0)(H̲0⋅∇)H̲1=-∇Π1,(3.6)ρΠ1=P1+(μ2)(2H̲0⋅H̲1),(3.7)∇⋅H̲1=0,(3.8)σρ1+ρ0(∇⋅u̲1)=0,(3.9)ρ0a2(∇⋅u̲1)=-σP1,(3.10)P1P0=ρ1ρ0+T1T0,(3.11)σH̲1=∇Λ(u̲1ΛH̲0),
where a is the sound speed in the gas defined by a=(γP0/ρ0)1/2,  and γ=Cp/Cv is the ratio of specific heats of the fluids.For the gas region surrounded by the fluid,(3.12)∇⋅H̲1g=0,(3.13)∇ΛH̲1g=0(thereisnocurrent).
Along the gas-liquid perturbed interface,(3.14)P1s=SR02(R1+R02∂2R1∂z2+∂2R1∂φ2).
By the aid of the series expansion (3.1) and the time-space dependence (3.4), based on the linear perturbation technique, every physical quantity Q(r,φ,z,t) could be expressed as(3.15)Q1(r,φ,z,t)=Q1*(r)exp⁡⁡[σt+i(kz+mφ)].
This means that any perturbed quantity could be expressed as an amplitude function of  r times the space-time dependence exp[σt+i(kz+mφ)].Upon utilizing the expansion, (3.15), (3.5), and (3.11) yield(3.16)(σ+ΩA2σ)u̲1=-∇Π1+iΩA2P1ka2ρ0e̲z,(3.17)σH̲1=ikH0u̲1+σH0ρ0a2P1e̲z,
where e̲z is a unit vector in the z-direction and ΩA is the Alfven wave frequency defined in terms of H0 as(3.18)ΩA=(μH02k2ρ0)1/2.
Combining the z-components of (3.16) and (3.17) which are(3.19)σ[1+(ΩAσ)2]u1z=-ikΠ1+iΩA2ka2ρ0P1,σH1z=H0(iku1z+σP1a2ρ0),
we obtain(3.20)H1z=(ik)2H0(σ2+ΩA2)[-Π1+μH02a2ρ0P1]+H0a2ρ0P1.
By substituting (3.6) into (3.20), we finally obtain(3.21)H1z=(σ2+a2k2)H0σ2a2ρ0P1.
Again, inserting (3.21) into (3.6) yields(3.22)Π1=(ξρ0)P1.
With(3.23)ξ=1+μH02(σ2+a2k2)σ2a2ρ0.
Consequently, the components u1r, u1φ, and u1z of u̲1 are given, from the vector equation (3.16), by(3.24)(σ2+ΩA2)u1r=-σ∂Π1∂r,(σ2+ΩA2)u1φ=-imσrΠ1,(σ2+ΩA2)u1z=-ikσ(-a2ξ+μH02)a2ξΠ1.
By inserting (3.24) into the nonsolenoid equation (3.9) which is(3.25)1r∂(ru1r)∂r+1r∂u1φ∂φ+∂u1z∂z=-σP1a2ρ0,
and by taking into account the space (φ,z)-dependence (3.15), we get(3.26)1rddr(rdΠ1dr)-(m2r2+η2)Π1=0,
with(3.27)η2=k2+σ2a2ξ.
Equation (3.26) is an ordinary second-order differential equation; its solution is given in terms of the ordinary Bessel functions of order m with imaginary argument. For the problem under consideration, apart from the singular solution, the finite solution of (3.26) is given by(3.28)Π1(r)=AKm(ηr).
Therefore, Π1(r,φ,z,t) is being(3.29)Π1(r,φ,z,t)=AKm(ηr)exp⁡⁡[σt+i(kz+mφ)],
where A is a constant of integration to be determined, while Km(ηr) is modified Bessel function of second kind of order m.It is worthwhile to mention here that by means of (3.29), the components of  u̲1(=(u1r,u1φ,u1z)) and H̲1(=(H1r,H1φ,H1z)) and also P1 could be identified from (3.16), (3.17), and (3.22) explicitly.Equation (3.13) means that the magnetic field H̲1g in the perturbation state can be derived by means of a scalar function that, by using (3.12), satisfies Laplace’s equation. The latter, by means of the expansion (3.15), transforms to an ordinary second-order differential equation whose solution is given in terms of cylindrical functions. For the problem under consideration, the nonsingular solution gives H̲1g in the form:(3.30)H̲1g=B∇[Im(kr)exp⁡⁡[σt+i(kz+mφ)]],
where B is an arbitrary constant of integration to be identified, while Im(kr) is the modified Bessel function of the first kind of order m.Based on the space dependence (3.15), the curvature pressure due to the capillary force along the gas-fluid interface is given, from (3.14), in the form:(3.31)P1s=SR02(1-m2-k2R02)R1.


## 4. Boundary Conditions

The solution of the relevant perturbation equations (3.5)–(3.14) represented by (3.15)–(3.31) and the solution of  the unperturbed system of equations represented by (2.16) must satisfy appropriate boundary conditions. Under the present circumstances for the problem at hand, these boundary conditions are given as follows.(i)
The normal component of the magnetic field must be continuous across the gas-fluid interface (3.3) at r=R0. This condition may be formulated as follows:
(4.1)N̲o⋅H̲1+N̲1⋅H̲o=N̲o⋅H̲1g+N̲1⋅H̲og,
where N̲(=N̲0+εN̲1) is a unit outward vector normal to the gas-fluid interface given by (2.15), with
(4.2)F(r,φ,z,t)=r-R0-R1=0
being the equation of the perturbed interface. Then,
(4.3)N̲o=(1,0,0),N̲1=(0,-imRo,-ik)R1.By substituting from (4.3) and the required expressions concerning H̲0,H̲1,H̲0g, and H̲1g into (4.1), we obtain
(4.4)B=imβH0xIm′(x),
where x(=kR0) is the dimensionless longitudinal wave number.(ii)
The normal componentur of the velocity u̲ must be compatible with the velocity of the deformed gas-fluid interface (3.3) at r=R0. This condition reads
(4.5)N̲0⋅u̲1+N̲1⋅u̲0=∂r∂t.By substitutingu̲0,u̲1,N̲0,N̲1,and r into (4.5), this yields
(4.6)A=-(σ2+ΩA2)ηKm′(ηR0).As we seen, complete nonsingular solutions for the variables of the problem have been obtained. Here, for the aim of stability theory, one has to make one more step to identify the stability criterion upon applying some compatibility condition.(iii)
This compatibility condition states that the normal component of the total stress tensor concerning the kinetic and magnetic pressures must be discontinuous by the curvature pressure due to the capillary force, across the gas-fluid interface (3.3) at r=R0.Mathematically, this reads(4.7)P1+R1∂P0∂r+(μ2)[(2H̲0⋅H̲1)+R1∂∂r(H̲0⋅H̲0)]=P1s+(μ2)[(2H̲0g⋅H̲1g)+R1∂∂r(H̲0g⋅H̲0g)].Upon substituting into the condition (4.7) about the different variables, the following eigenvalue relation is obtained:
(4.8)σ2=-SρR03(1-m2-x2)yKm′(y)Km(y)+μH02ρoRo2{-x2+[-β2+(mβ)2Im(x)xIm′(x)]yKm′(y)Km(y)},
where y(=ηR0) is the dimensionless longitudinal wave number due to compressibility.


## 5. Limiting Cases

Equation (4.8) is the required capillary instability eigenvalue relation for a compressible gas cylinder surrounded by a liquid and pervaded by a transverse varying magnetic field. It relates the temporal amplification σ or alternatively the oscillation frequency ω (i.e., if σ(=iω) is imaginary) with the dimensionless (ordinary) longitudinal x and y (compressible) wave numbers, azimuthal wave number m, the fundamental quantity (S/ρ0R03)-1/2, as well as (μH02/ρ0R02)-1/2 as a unit of time, the modified Bessel functions Im and Km of the first and second kind of order m and their derivatives of different arguments, the parameter β of the magnetic field  pervading the interior of the gas cylinder, and with the parameters S,ρ0,R0,μ, and H0 of the problem.Several reported works may be obtained as limiting cases from the general result (4.8).For incompressible(a→∞), nonconducting (μ=0,H0=0) fluid and axisymmetric perturbation m=0, the relation (4.8) yields(5.1)σ2=-SρR03(xK0′(x)K0(x))(1-x2),
where (a→∞), we find (y→x). Upon using the recurrence relation,(5.2)K0′(x)=-K1(x),
the relation (5.1) becomes(5.3)σ2=SρR03(xK1(x)K0(x))(1-x2).
The relation (5.3) is indicated by Chandrasekhar [2, page 540], for the first time as a dispersion relation for the mirror case of a full fluid cylinder surrounded by vacuum. For the discussion of this relation, we may refer to Radwan and Elazab [6] as we neglect the contribution of the viscosity there.If we suppose that(a→∞) so (y→x), H0=0 for m≠0, the relation (4.8) gives(5.4)σ2=-SρR03(xKm′(x)Km(x))(1-m2-x2).
This relation coincides with the relation given by Drazin and Reid [3].The capillary eigenvalue relation of a compressible hollow cylinder is given from (4.8) as H0=0, in the form:(5.5)σ2=-SρR03(1-m2-x2)yKm′(y)Km(y).
The magnetodynamic eigenvalue relation of a compressible hollow cylinder is given from the general relation (4.8) as =0, in the form(5.6)σ2=μH02ρoRo2{-x2+[-β2+(mβ)2Im(x)xIm′(x)]yKm′(y)Km(y)}.


## 6. Stability Discussions

The general eigenvalue relation (4.8) of the model under consideration is a quadratic relation in σ2. Therefore, we have to distinguish between the following different cases.(i)
The model will be ordinary stable, asσ2 is negative.(ii)
The model will be ordinary unstable, asσ2 is positive.(iii)
The model will be marginally stable, asσ2 is zero.In order to judge such cases and investigate the stability of the present model, we have to write down about the character and behavior of the modified Bessel functionsIm  and Km and their derivatives.Consider the recurrence relations (cf. [16]):(6.1)2Im′(x)=Im-1(x)+Im+1(x),2Km′(x)=-Km-1(x)-Km+1(x).
For each nonzero real value of x. Also Im(x) is positive definite and monotonic increasing, while Km(x) is monotonic decreasing but never negative:(6.2)Im(x)>0,Km(x)>0.
Therefore, we may see, for a nonzero real value of x, that(6.3)Im′(x)>0,Km′(x)<0.In view of (3.27), we see that y→x as a→∞, where the fluid in this case is an incompressible one. Consequently, the stability of the hollow cylinder under the action of the combined effects of the capillary and electromagnetic forces could be discussed.As the model of a hollow cylinder is acted by the capillary force while the effect of the electromagnetic force is neglected, the eigenvalue relation (4.8)  reduces to (5.5). By means of the relation, the capillary stable and unstable regions of the hollow cylinder could be identified. By an appeal to the relations (6.1)  and the inequalities (6.2) and (6.3), the sign of  σ2 depends on the sign of the quantity (1-m2-x2). Henceforth, we have the different cases:(6.4){σ*2=0,asx=1,σ*2<0as1<x<∞σ*2>0as0<x<1,form=0,σ*2<0,as0≤x<∞,form≠0,
where(6.5)σ*=σ(Sρ0R03)-1/2
is the dimensionless growth rate at the instant of time t.This means that the hollow cylinder is stable in the case1≤x<∞ for m=0 and 0≤x<∞ for m≥1, while it is unstable only as 0≤x≤1 for m=0. This results, of course, as a-1→0. However, for ordinary values of a: in such case, we have the modified Bessel functions with argument y which in turn is a function of the compressibility factor a(=(γP0/ρ0)1/2).  Now, for x≠0,y≠0 and y>x, we have(6.6)Im(y)>Im(x),Km(x)>Km(y),Im′(y)>Im′(x),Km′(x)>Km′(y).
In view of these inequalities, the discussion of the relation (5.5) reveals that the compressibility has a stabilizing effect in all (m=0 and m≠0) modes of perturbation for  all wavelengths. The stabilizing effect becomes stronger for high compressibility and, therefore, it acts to overcome the capillary destabilizing influence.As the hollow cylinder model is acted upon the electromagnetic force due to the pervading magnetic fields in the initial states and pressure gradient forces asS=0, the eigenvalue relation of such case is given, in its general form, by the relation (5.6). The effect of the axial magnetic field in the liquid region is represented by the term -x2 following μH02/ρ0R02. It contributes as a negative part in σ2, that is, σ is imaginary in this case. This means that it has a stabilizing effect on the model. The transverse varying magnetic field pervaded into the gas region is presented by the terms -β2yKm′(y)/Km(y) and (m2β2)(yIm(x)Km′(y)/xIm′(x)Km(y)) following the natural quantity μH02/ρ0R02. In the axisymmetric mode, m=0, the transverse magnetic field is purely destabilizing. In the nonaxisymmetric modes, m≥1, the transverse magnetic field is purely destabilizing in the term -β2yKm′(y)/Km(y) while it is stabilizing due to the other terms. Therefore, we conclude that the transverse magnetic field is stabilizing or destabilizing according to restrictions.Consequently, the electromagnetic forces in the gas and liquid regions  are strongly stabilizing  the model in the axisymmetric modem=0, while they have a destabilizing influence in m≥1 modes in the gas region only. Also, in such magnetodynamic case the compressibility has a stabilizing tendency not only in m=0 mode but also in the non-axisymmetric modes m≥1. Its influence is decreasing the magnetodynamic unstable domain in m≥1 and simultaneously increasing those of stability in m=0 mode.One has to mention here it is argued that when the effects of materials compressibility are considered, the growth rate value may be reduced in comparison with the incompressible case. This is due to the fact that compression absorbs some of the energy which would, otherwise, go into fluid motion and causes more instability. For this reason, it is stated that the compressibility has a stabilizing tendency.It is worth to mention here that the stabilizing effect of the electromagnetic forceμ(∇ΛH̲)ΛH̲ may be interpreted as follows.This force is interpreted as arising from the action on the fluid of Maxwell's stresses: a magnetic tensionμ(H̲·H̲)/2 per unit area along the magnetic lines of force and equal magnetic pressure acting in all the directions in the conducting fluid. Taking into account that the latter is not perpendicular to the magnetic lines of force and acting in all directions because the diffusion term is neglected in the evolution equations of the magnetic field (2.9). Due to these stresses, the lines of force are able to endow the fluid with a sort of rigidity.


## 7. Numerical Analysis

The main purpose of the numerical analysis is that one could determine exactly where are the MHD stable and unstable domains of the model under consideration influenced by the combined effects of the electromagnetic and capillary forces. In addition, we may identify the critical points which separate the stable and unstable domains for different values of the problem parameters.To perform that, we have to write down the eigenvalue relation (4.8) in dimensionless form in the most important mode m=0 of perturbation. Now, since the natural quantity (S/ρ0R03)-1/2 has a unit of (time)−1, we have σ*=σ(S/ρ0R03)-1/2 which is the dimensionless growth rate. For m=0, the relation (4.8) leads to(7.1)σ*2=(1-x2)yK1(y)K0(y)+(H0Hs)2{-x2+β2yK1(y)K0(y)},
where σ* is the nondimensional growth rate.Ifσ*2 is negative, we write σ*=iω* (with i=-1 being the imaginary factor) where ω/2π is the oscillation frequency. The notation Hs2=S/μR0 is with Hs as a unit of magnetic field, and we have used the notation I0′=I1 and K0′=-K1. The dispersion relation (7.1) has been computed for all short and long wavelengths in the new form in relation (7.2) which is(7.2)σ*2=(1-x2)yK1(y)K0(y)+M{-x2+β2yK1(y)K0(y)}
such calculations have been elaborated for different values of β and M for regular values of y. The values of  σ* corresponding the unstable domains and those of corresponding the stable domains are collected, tabulated, and presented graphically.It is found that there are many features of interest in this numerical analysis as we see in the following.(i)
Forβ=0.7 and x=1, see Figure 2. Corresponding to  M=0.5, 1, 1.5, 2, 2.5, and 3, it is found that the unstable domains are 0<y<1.6209,0<y<1.6209, 0<y<1.6209, 0<y<1.6209, 0<y<1.6209, and 0<y<1.6209 while the neighbouring stable domains are 1.6209≤y<∞,  1.6209≤y<∞,  1.6209≤y<∞, 1.6209≤y<∞, 1.6209≤y<∞,  and 1.6209≤x<∞ where the equalities correspond to the marginal stability states.(ii)
Forβ=0.8 and x=1, see Figure 3. Corresponding to  M=0.5, 1, 1.5, 2, 2.5, and 3, it is found that the unstable domains are 0<y<1.1383,0<y<1.1383, 0<y<1.1383, 0<y<1.1383, 0<y<1.1383, and 0<y<1.1383 while the neighbouring stable domains are 1.1383≤y<∞,  1.1383≤y<∞,  1.1383≤y<∞, 1.1383≤y<∞, 1.1383≤y<∞,  and 1.1383≤x<∞ where the equalities correspond to the marginal stability states.(iii)
Forβ=0.9 and x=1, see Figure 4. Corresponding to  M=0.5, 1, 1.5, 2, 2.5, and 3, it is found that the unstable domains are 0<y<0.829,0<y<0.829, 0<y<0.829, 0<y<0.829, 0<y<0.829, and 0<y<0.829 while the neighbouring stable domains are 0.829≤y<∞,  0.829≤y<∞,  0.829≤y<∞, 0.829≤y<∞, 0.829≤y<∞,  and 0.829≤x<∞ where the equalities correspond to the marginal stability states.(iv)
Forβ=1 and x=1, see Figure 5. Corresponding to  M=0.5, 1, 1.5, 2, 2.5, and 3, it is found that the unstable domains are 0<y<0.6179,0<y<0.6179, 0<y<0.6179, 0<y<0.6179, 0<y<0.6179, and 0<y<0.6179 while the neighbouring stable domains are 0.6179≤y<∞,  0.6179≤y<∞,  0.6179≤y<∞, 0.6179≤y<∞, 0.6179≤y<∞,  and 0.6179≤x<∞ where the equalities correspond to the marginal stability states.(v)
Forβ=1.1 and x=1, see Figure 6. Corresponding to  M=0.5, 1, 1.5, 2, 2.5, and 3, it is found that the unstable domains are 0<y<0.4439,0<y<0.4439, 0<y<0.4439, 0<y<0.4439, 0<y<0.4439, and 0<y<0.4439 while the neighbouring stable domains are 0.4439≤y<∞,  0.4439≤y<∞,  0.4439≤y<∞, 0.4439≤y<∞, 0.4439≤y<∞,  and 0.4439≤x<∞ where the equalities correspond to the marginal stability states.Figure 2
Stable and unstable domains forβ=0.7 and x=1.Figure 3
Stable and unstable domains forβ=0.8 and x=1.Figure 4
Stable and unstable domains forβ=0.9 and x=1.Figure 5
Stable and unstable domains forβ=1 and x=1.Figure 6
Stable and unstable domains forβ=1.1 and x=1.


## 8. Conclusion

From the numerical discussion, we deduce that the compressibility has a strong stabilizing tendency for all wavelengths. The capillary force is destabilizing for a small domain of long wavelengths in the axisymmetric mode while it is stabilizing in all the rest. The electromagnetic force interior the gas is stabilizing. The electromagnetic force in the liquid region is stabilizing also. These results are in good agreement with the analytical discussions of relation (4.8).



---

*Source: 531250-2012-03-05.xml*