Table of Contents Author Guidelines Submit a Manuscript
Mathematical Problems in Engineering
VolumeΒ 2012Β (2012), Article IDΒ 325423, 12 pages
Research Article

Magnetohydrodynamic Stability of Streaming Jet Pervaded Internally by Varying Transverse Magnetic Field

Basic and Applied Sciences Department, College of Engineering and Technology, Arab Academy for Science & Technology and Maritime Transport (AASTMT), P.O. Box 2033 Elhorria, Cairo, Egypt

Received 9 October 2011; Accepted 12 December 2011

Academic Editor: AlexeiΒ Mailybaev

Copyright Β© 2012 Alfaisal A. Hasan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


The Magnetohydrodynamic stability of a streaming cylindrical model penetrated by varying transverse magnetic field has been discussed. The problem is formulated, the basic equations are solved, upon appropriate boundary conditions the eigenvalue relation is derived and discussed analytically, and the results are verified numerically. The capillary force is destabilizing in a small axisymmetric domain 0<π‘₯<1 and stabilizing otherwise. The streaming has a strong destabilizing effect in all kinds of perturbation. The toroidal varying magnetic field interior the fluid has no direct effect at all on the stability of the fluid column. The axial exterior field has strong stabilizing effect on the model. The effect of all acting forces altogether could be identified via the numerical analysis of the stability theory of the present model.

1. Introduction

The classical stability analysis of a full fluid jet has been extensively studied (cf. Rayleigh [1], Chandrasekhar [2], Robert [3], Chenng [4], Kendall [5], and Drazin and Reid) [6] and documented by Radwan [7]. The latter author studied the hydromagnetic stability of a fluid jet pervaded by uniform magnetic field for axisymmetric perturbation. Radwan et al. [8] developed the magnetohydrodynamic stability of that model for all axisymmetric and nonaxisymmetric modes subject to electromagnetic forces generated due to constant magnetic field. The stability of different models under the action of self-gravitating force in addition to other forces has been elaborated by Radwan and Hasan [9, 10]. Hasan [11] has discussed the stability of oscillating streaming fluid cylinder subject to combined effect of the capillary, self-gravitating, and electrodynamic forces for all axisymmetric and nonaxisymmetric perturbation modes. The purpose of the present work is to investigate the magnetodynamic stability of a fluid jet pervaded by transverse varying magnetic field while its surrounding tenuous medium is penetrated by uniform magnetic field.

2. Basic State

We consider an incompressible, inviscid, and conducting fluid column of radius π‘…π‘œ in the basic state.

A transverse varying magnetic field is assumed to be pervaded interior the fluid, namely,π»π‘œ=𝐻0,π‘œπ‘Ÿπ‘…π‘œξ‚Ά,0.(2.1) The surrounding tenuous medium of the fluid jet is penetrated by the axial homogenous magnetic field:π»π‘œπ‘‘π‘›=ξ€·0,0,π›Όπ»π‘œξ€Έ,(2.2) where 𝛼 is parameter while π»π‘œ is the intensity of the transverse magnetic field across the fluid surface at π‘Ÿ=π‘…π‘œ. The fluid streams in the initial state with velocity:𝑒=ξ€·0,0,π‘ˆπ‘œξ€Έ,(2.3) where π‘ˆπ‘œ is (uniform) the amplitude of the velocity 𝑒. The components of π»π‘œandπ»π‘œπ‘‘π‘› are considered along the utilizing cylindrical coordinates (π‘Ÿ,πœ‘,𝑧) with the 𝑍-axis coinciding with the axis of the cylinder (see Figure 1).

Figure 1: Sketch for MHD fluid jet.

The basic magnetodynamic equations appropriate for studying such a problem may be formulated as follows interior the fluid cylinder:πœŒπ‘‘π‘’ξ€·π‘‘π‘‘=βˆ’βˆ‡π‘ƒ+πœ‡βˆ‡βˆ§π»ξ€Έβˆ§π»,(2.4)βˆ‡β‹…π‘’=0,(2.5)πœ•π»ξ€·π‘’πœ•π‘‘=βˆ‡βˆ§βˆ§π»ξ€Έ,(2.6)βˆ‡β‹…π»π‘ƒ=0,(2.7)𝑠1=π‘‡π‘Ÿ1+1π‘Ÿ2ξ‚Άξ‚΅=π‘‡βˆ‡β‹…βˆ§π‘ξ‚Ά,(2.8) where π‘Ÿ1andπ‘Ÿ2 are the principle radii of curvature and 𝑇 is the surface tension coefficient. 𝑁𝑠 is, the unit outward vector normal to the surface, given by βˆ§π‘π‘ =βˆ‡π‘“(π‘Ÿ,πœ‘,𝑧;𝑑)||||βˆ‡π‘“(π‘Ÿ,πœ‘,𝑧;𝑑).(2.9) The space βˆ‡ and the time 𝑑/𝑑𝑑 operators are given byξ‚΅πœ•βˆ‡β‰‘,1πœ•π‘Ÿπ‘Ÿπœ•,πœ•πœ•πœ‘ξ‚Ά,π‘‘πœ•π‘§=πœ•π‘‘π‘‘+ξ€·π‘’πœ•π‘‘ξ€Έβ‹…βˆ‡.(2.10) Here 𝜌,𝑒,and𝑃 are the fluid mass density, velocity vector and kinetic pressure, 𝐻and𝐻𝑑𝑛 are the magnetic field intensity of fluid and of tenuous medium surrounding the fluid jet. πœ‡ is the magnetic field permeability coefficient.

In the surrounding tenuous medium outside the fluid jet, the basic equations are βˆ‡β‹…π»π‘‘π‘›=0,(2.11)βˆ‡βˆ§π»π‘‘π‘›=0,thereisnocurrent.(2.12)

3. Unperturbed State

The unperturbed state is studied and the kinetic fluid pressure π‘ƒπ‘œ is identified.

From now on, the quantities with index π‘œ refer to equilibrium quantities.

In such state, (2.4) givesβˆ‡ξ‚΅π‘ƒπ‘œ+πœ‡2𝐻2π‘œπ‘Ÿ2𝑅2π‘œξ‚Ά=0.(3.1) From which, we getπ‘ƒπ‘œ+πœ‡2𝐻2π‘œπ‘Ÿ2𝑅2π‘œ=π‘˜(const.).(3.2) By applying the balance of the pressure across the fluid tenuous interface at π‘Ÿ=π‘…π‘œ, we finally obtainπ‘ƒπ‘œ=π‘‡π‘…π‘œ+πœ‡π»2π‘œ2𝛼2ξ€Έβˆ’1,(3.3) where 𝛼 satisfies the restriction 𝛼β‰₯1 in order that the pressure π‘ƒπ‘œ must be positive as we neglect the surface tension contribution in the unperturbed state.

4. Perturbation Analysis

For small departures from the initial state, every physical quantity could be expressed as𝑄(π‘Ÿ,πœ‘,𝑧;𝑑)=π‘„π‘œ(π‘Ÿ)+πœ€(𝑑)𝑄1(π‘Ÿ,πœ‘,𝑧),(4.1) where 𝑄 stands for 𝑒,𝑃,𝐻,and𝐻𝑑𝑛 while, πœ€(𝑑) is the amplitude of the perturbation, (cf. [2]) beingπœ€(𝑑)=πœ€π‘œexp(πœŽπ‘‘).(4.2) Here 𝜎 is the growth rate at time 𝑑 and πœ€π‘œ(=πœ€at𝑑=0) that is the initial amplitude.

By inserting the expansion (4.1) into (2.4)–(2.11), the linearized perturbation equations are given by πœ•π‘’1+ξ€·π‘’πœ•π‘‘π‘œξ€Έπ‘’β‹…βˆ‡1βˆ’πœ‡πœŒπ»ξ€·ξ€·1ξ€Έπ»β‹…βˆ‡π‘œ+ξ€·π»π‘œξ€Έπ»β‹…βˆ‡1=βˆ’βˆ‡1,(4.3)βˆ‡β‹…π‘’1=0,(4.4)πœ•π»1=ξ€·π»πœ•π‘‘π‘œξ€Έπ‘’β‹…βˆ‡1βˆ’ξ€·π‘’1ξ€Έπ»β‹…βˆ‡π‘œβˆ’ξ€·π‘’π‘œξ€Έπ»β‹…βˆ‡1,(4.5)βˆ‡β‹…π»1=0,(4.6)βˆ‡β‹…π»1𝑑𝑛=0,(4.7)βˆ‡βˆ§π»1𝑑𝑛𝑃=0,(4.8)1𝑠=βˆ’π‘‡π‘…2π‘œξ‚΅π‘…1+πœ•2𝑅1πœ•πœ‘2+𝑅2π‘œπœ•2𝑅1πœ•π‘§2ξ‚Ά,(4.9) where 𝜌∏1(=π‘ƒπ‘œ+(πœ‡/2)(𝐻⋅𝐻)1) (see [3]) is the total magnetohydrodynamics pressure which is the sum of kinetic and magnetic pressures.

Based on the expansion (4.1), the perturbed fluid-tenuous interface may be expressed asπ‘Ÿ=π‘…π‘œ+𝑅1(πœ‘,𝑧,𝑑),(4.10) with𝑅1=πœ€π‘œexp(𝑖(π‘˜π‘§+π‘šπœ‘)+πœŽπ‘‘).(4.11) Here 𝑅1 is the elevation of the surface wave measured from the unperturbed position, where π‘˜ is the longitudinal wavenumber and π‘š is the azimuthal wavenumber. By an appeal to the perturbation technique, we may write𝑄1(π‘Ÿ,πœ‘,𝑧)=π‘„βˆ—1(π‘Ÿ)exp(𝑖(π‘˜π‘§+π‘šπœ‘)).(4.12) Consequently, the linearized system of (4.3)–(4.8) is solved and the fluctuation parts are determined.

In the radial direction, we obtainξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έπ‘’1π‘Ÿ=πœ‡π»π‘œπœŒπ‘…π‘œξ€·π‘–π‘šπ»1π‘Ÿξ€Έβˆ’πœ•βˆ1,ξ€·πœ•π‘ŸπœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έπ»1π‘Ÿ=π‘–π‘šπ»π‘œπ‘…π‘œξ€·π‘’1π‘Ÿξ€Έ.(4.13) In the transverse direction, we haveξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έπ‘’1πœ‘=βˆβˆ’π‘–π‘š1π‘Ÿ+πœ‡πœŒπ‘…π‘œξ€·π‘–π‘šπ»1πœ‘+𝐻1π‘Ÿξ€Έπ»π‘œ,ξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έπ»1πœ‘=π‘–π‘šπ»π‘œπ‘…π‘œπ‘’1πœ‘βˆ’π»π‘œπ‘…π‘œπ‘’1π‘Ÿ.(4.14) In the longitudinal 𝑧-direction, we getξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έπ‘’1𝑧=βˆ’π‘–π‘˜1+π‘–πœ‡π»π‘œπ‘šπœŒπ‘…π‘œπ»1𝑧,ξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έπ»1𝑧=π‘–π‘šπ»π‘œπ‘…π‘œπ‘’1𝑧.(4.15) From (4.13)–(4.15), we obtain𝑒1π‘Ÿ=βˆ’ξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έξ‚€ξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έ2+πœ‡π»2π‘œπ‘š2/πœŒπ‘…2π‘œξ‚πœ•βˆ1,π‘’πœ•π‘Ÿ1πœ‘=ξ€·βˆ’π‘–π‘šπœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έξ‚€ξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έ2+πœ‡π»2π‘œπ‘š2/πœŒπ‘…2π‘œξ‚βˆ1π‘Ÿ,𝑒1π‘Ÿ=ξ€·βˆ’π‘–π‘˜πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έξ‚€ξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έ2+πœ‡π»2π‘œπ‘š2/πœŒπ‘…2π‘œξ‚βˆ1.(4.16) By taking the divergence of (4.3) and using (4.4), we obtainβˆ‡21=0.(4.17) Apart from the singular solution, the finite solution (as π‘Ÿβ†’0) of this equation is given by 1=π΄πΌπ‘š(π‘˜π‘Ÿ)exp(πœŽπ‘‘+𝑖(π‘˜π‘§+π‘šπœ‘)).(4.18) Now (4.8) means that 𝐻1𝑑𝑛 could be derived from a scalar function, say πœ™1, such that 𝐻1𝑑𝑛=βˆ’βˆ‡πœ™1.(4.19) Combining (4.19) with (4.7), we getβˆ‡2πœ™1=0.(4.20) The nonsingular solution of (4.20) is identified and given by𝐻1𝑑𝑛=βˆ‡π΅πΎπ‘šξ€»(π‘˜π‘Ÿ)exp(πœŽπ‘‘+𝑖(π‘˜π‘§+π‘šπœ‘)).(4.21) Here πΌπ‘š(π‘˜π‘Ÿ)andπΎπ‘š(π‘˜π‘Ÿ) are the modified Bessel functions, respectively, of the first and second kind of order π‘š while 𝐴 and 𝐡 are constants of integration.

5. Eigenvalue Relation

The solutions which are given by (3.3) and (4.13)–(4.21) of the basic equations (2.4)–(2.12) must satisfy certain boundary conditions. Under the present circumstances, these boundary conditions are the following.

(i) The normal component of the velocity vector 𝑒 of the fluid must be compatible with the velocity of the perturbed fluid across the fluid-tenuous interface at π‘Ÿ=π‘…π‘œ. This condition givesβˆ’ξ€·π΄=𝜎+π‘–π‘˜π‘ˆπ‘œξ€Έ2+πœ‡π»2π‘œπ‘š2/πœŒπ‘…2π‘œπ‘˜πΌβ€²π‘š(π‘₯),(5.1) where π‘₯(=π‘˜π‘…π‘œ) is the nondimension longitudinal wave number.

(ii) The normal component of the magnetic field must be continuous across the fluid-tenuous interface at π‘Ÿ=π‘…π‘œ. This condition yields𝐡=π‘–π›Όπ»π‘œπΎξ…žπ‘š(π‘₯).(5.2)

(iii) The total magnetohydrodynamic pressure ∏(𝜌1) interior the fluid must be discontinuous with magnetic pressure due to the electromagnetic force in the tenuous medium surrounding the fluid jet across the interface (4.10) at π‘Ÿ=π‘…π‘œ by the contribution of the capillary force. This condition reads𝑃1+𝑅1πœ•π‘ƒπ‘œξ€·π»πœ•π‘Ÿ+πœ‡π‘œβ‹…π»1ξ€Έ+𝑅1πœ•ξ‚€πœ‡πœ•π‘Ÿ2ξ€·π»π‘œβ‹…π»π‘œξ€Έξ‚ξ€·π»βˆ’πœ‡π‘œβ‹…π»1ξ€Έπ‘‘π‘›βˆ’π‘…1πœ•ξ‚€πœ‡πœ•π‘Ÿ2ξ€·π»π‘œβ‹…π»π‘œξ€Έπ‘‘π‘›ξ‚=𝑃1𝑠.(5.3)

Upon substitution from (3.3)–(5.2) into (5.3), following eigenvalue relation is derived ξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έ2=πœ‡π»2π‘œπœŒπ‘…2π‘œξ‚Έβˆ’π‘š2+𝛼2π‘₯2πΎπ‘š(π‘₯)πΌξ…žπ‘š(π‘₯)πΎξ…žπ‘š(π‘₯)πΌπ‘šξ‚Ή+𝑇(π‘₯)πœŒπ‘…3π‘œξ‚Άξ‚΅π‘₯πΌξ…žπ‘š(π‘₯)πΌπ‘šξ‚Άξ€·(π‘₯)1βˆ’π‘š2βˆ’π‘₯2ξ€Έ.(5.4)

6. Discussions

Equation (5.4) is the desired capillary dispersion relation of a streaming fluid column pervaded by varying transverse magnetic field (cf. (2.1)) and surrounded by uniform magnetic field. This relation is valid for all axisymmetric and nonaxisymmetric perturbation modes. It relates the growth rate 𝜎 with the wavenumbers π‘˜andπ‘š, the problem parameters πœ‡,π»π‘œ,𝜌,π‘…π‘œ and the speed π‘ˆπ‘œ of the streaming fluid and with the modified Bessel functions πΌπ‘šandπΎπ‘š and their derivatives.

For axisymmetric perturbation (π‘š=0), the general dispersion (5.4) reduces to ξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έ2=ξ‚΅βˆ’πœ‡π»2π‘œπœŒπ‘…2π‘œξ‚Άξ€·π›Ό2π‘₯2𝐾0(π‘₯)𝐼1(π‘₯)𝐾1(π‘₯)𝐼0+𝑇(π‘₯)πœŒπ‘…3π‘œξ‚Άξ‚΅π‘₯𝐼1(π‘₯)𝐼0ξ‚Άξ€·(π‘₯)1βˆ’π‘₯2ξ€Έ.(6.1)

For the lowest non-symmetric perturbation mode (π‘š=1), the dispersion relation (5.4) yields ξ€·πœŽ+π‘–π‘˜π‘ˆπ‘œξ€Έ2=πœ‡π»2π‘œπœŒπ‘…2π‘œξƒ¬βˆ’1+𝛼2π‘₯2𝐾1(π‘₯)𝐼′1(π‘₯)πΎξ…ž1(π‘₯)𝐼1(ξƒ­βˆ’ξ‚΅π‘‡π‘₯)πœŒπ‘…3π‘œξ‚Άξƒ©π‘₯3πΌξ…ž1(π‘₯)𝐼1ξƒͺ(π‘₯).(6.2)

For the higher nonaxisymmetric perturbation π‘šβ‰₯2, the eigenvalue relation of such case may be easily obtained from (5.4) as π‘šβ‰₯2.

By the use of the recurrence relations (cf. [12]) 2πΌξ…žπ‘š(π‘₯)=βˆ’πΌπ‘š+1(π‘₯)+πΌπ‘šβˆ’1(π‘₯),2πΎξ…žπ‘š(π‘₯)=βˆ’πΎπ‘šβˆ’1(π‘₯)βˆ’πΎπ‘š+1,(π‘₯)(6.3) and that πΌπ‘š(π‘₯) is monotonically increasing and positive definite for every nonzero value of π‘₯, we have πΌπ‘š(π‘₯)>0.(6.4) The function πΎπ‘š(π‘₯) is monotonically decreasing but never negative for since π‘₯β‰ 0 valuesπΎπ‘š(π‘₯)>0.(6.5) Upon using (6.4) and (6.5) for the relations (6.3), we getπΌξ…žπ‘š(π‘₯)>0,πΎξ…žπ‘š(π‘₯)<0,(6.6) for all axisymmetric and nonaxisymmetric perturbation modes π‘šβ‰₯0. By an appeal to the recurrence relations (6.3) and the inequalities (6.6) for the stability criterion (5.4), we deduce the following.(1)The toroidal varying magnetic field interior the fluid has no influence at all on the stability of the fluid column.(2)The uniform exterior magnetic field is stabilizing for all short and long wavelengths in all kinds of perturbation π‘šβ‰₯0.(3)The streaming is strongly destabilizing.

Therefore, the magnetohydrodynamic streaming fluid column is not completely stable in the axisymmetric mode π‘š=0, but there will be exist some unstable domains.

The stability discussions for the case of the lowest nonaxisymmetric mode (π‘š=1), may be carried out by utilizing the relation (6.2). It is found that both the interior toroidal and exterior axial magnetic fields are stabilizing. Note that the streaming has a strong destabilizing effect.

Therefore, the present model is (MHD) unstable in the lowest nonaxisymmetric mode (π‘š=1) of perturbation.

In a similar way, one may show that the fluid column is MHD unstable for any higher nonaxisymmetric mode π‘šβ‰₯2 of perturbation.

Physically, the stabilizing effect of the exterior magnetic field is expected because it has been assumed that the pervading magnetic field is uniform.

Moreover, the stabilizing effect of the toroidal magnetic field in the fluid region is due to the influence of Lorentz force that comes out from the interaction of the magnetic induction and the electric current produced due to the pervading magnetic field. Indeed such electromagnetic force causes the following stresses [7].

The magnetic pressure (πœ‡/2)(𝐻⋅𝐻) per unit area acting in all directions of the fluid (resistivity is neglected here) and an equal magnetic tension (πœ‡/2)(𝐻⋅𝐻) per unit area acting along the magnetic lines of force. Due to these stresses the lines of force are able to endow the fluid with a sort of rigidity. The magnetic fields exert strong influence not only to the axisymmetric mode (π‘š=0) that causes only the bending of the magnetic lines of force but also to nonaxisymmetric modes that lead to twisting of the lines of force.

7. Numerical Analysis

The numerical analysis has been carried out in order to identify and examine the magnetic field influence and surface tension and also the effect of the streaming on stability of the model. In addition to that the oscillation states and the transition points from these states to those of instability may be also determined for given values of the magnetic field intensity.

This has been elaborated by computing the nondimension dispersion relationπœŽξ”π‘‡/πœŽπ‘…3π‘œ=ξƒŽξ‚΅π»π‘œπ»π‘ π‘šξ‚Άξ‚΅2+𝛼2π‘₯2πΎπ‘š(π‘₯)πΌξ…žπ‘š(π‘₯)πΎξ…žπ‘š(π‘₯)πΌπ‘šξ‚Ά+ξ‚΅(π‘₯)π‘₯πΌξ…žπ‘š(π‘₯)πΌπ‘šξ‚Άξ€·(π‘₯)1βˆ’π‘š2βˆ’π‘₯2ξ€Έ+π‘ˆβˆ—(7.1) in the computer simulation for the most dangerous sausage mode π‘š=0 for the different values of (𝛼)and(π»π‘œ/𝐻𝑠), where 𝐻𝑠(=(𝑇/πœ‡π‘…π‘œ)) has the dimension of magnetic field while 𝜎(=π‘–πœ”) is the growth rate and πœ” is the oscillation frequency.

The numerical data associated with πœ”/(𝑇/πœŒπ‘…3π‘œ)1/2 corresponding to the stable states and those associated with 𝜎/(𝑇/πœŒπ‘…3π‘œ)1/2 corresponding to the unstable states are collected, tabulated, and presented graphically. See Figures 2, 3, 4, and 5. There are many features of interest in these numerical illustration. In the absence of the streaming effect it is found that the domains of instability are decreasing with increasing (π»π‘œ/𝐻𝑠) values even for different values of 𝛼.

Figure 2: MHD stable and unstable domains as 𝛼=1,π‘š=0,π‘ˆβˆ—=0.
Figure 3: MHD stable and unstable domains as 𝛼=1,π‘š=0,π‘ˆβˆ—=0.3.
Figure 4: MHD stable and unstable domains as 𝛼=1,π‘š=0,π‘ˆβˆ—=0.7.
Figure 5: MHD stable and unstable domains as 𝛼=1,π‘š=0,π‘ˆβˆ—=1.0.

For 𝛼=1 see Figure 2: corresponding to β„Ž=(π»π‘œ/𝐻𝑠)=0.1,0.3,0.7,1.0,1.5,and3.5, it is found that the unstable domains are 0≀π‘₯<0.9965,0≀π‘₯<0.9693,0≀π‘₯<0.8496,0≀π‘₯<0.7307,0≀π‘₯<0.5436,0≀π‘₯<0.2083, while the neighboring stable domains, respectively, are 0.9965≀π‘₯<∞,0.9693≀π‘₯<∞,0.8496≀π‘₯<∞,

0.7307≀π‘₯<∞,0.5436≀π‘₯<∞,and 0.2083≀π‘₯<∞.

For (𝛼,π‘ˆβˆ—)=(1,0.3): corresponding to β„Ž=(π»π‘œ/𝐻𝑠)=0.1,0.3,0.7,1.0,1.5,and3.5, it is found that the unstable domains are 0≀π‘₯<1.0819,0≀π‘₯<1.0561,0≀π‘₯<0.9491,0≀π‘₯<0.8464,0≀π‘₯<0.6791,0≀π‘₯<0.3678, while the neighboring stable domains, respectively, are 1.0819≀π‘₯<∞,1.0561≀π‘₯<∞,0.9491≀π‘₯<∞,

0.8464≀π‘₯<∞,0.6791≀π‘₯<∞,0.3678≀π‘₯<∞,see Figure 3.

For (𝛼,π‘ˆβˆ—)=(1,0.7): corresponding to β„Ž=(π»π‘œ/𝐻𝑠)=0.1,0.3,0.7,1.0,1.5,and3.5, it is found that the unstable domains are 0≀π‘₯<1.2985,0≀π‘₯<1.2769,0≀π‘₯<1.1735,0≀π‘₯<1.0867,0≀π‘₯<0.9277,0≀π‘₯<0.5519, while the neighboring stable domains, respectively, are 1.2985≀π‘₯<∞,1.2769≀π‘₯<∞,1.1735≀π‘₯<∞,

1.0867≀π‘₯<∞,0.9277≀π‘₯<∞,0.5519≀π‘₯<∞,see Figure 4.

For (𝛼,π‘ˆβˆ—)=(1,1.0): corresponding to β„Ž=(π»π‘œ/𝐻𝑠)=0.1,0.3,0.7,1.0,1.5,and3.5, it is found that the unstable domains are 0≀π‘₯<1.4661,0≀π‘₯<1.4453,0≀π‘₯<1.3455,0≀π‘₯<1.2500,0≀π‘₯<1.0882,0≀π‘₯<0.6765, while the neighboring stable domains, respectively, are 1.4661≀π‘₯<∞,1.4453≀π‘₯<∞,1.3455≀π‘₯<∞,1.2500≀π‘₯<∞,1.0882≀π‘₯<∞,0.6765≀π‘₯<∞,see Figure 5.

8. Conclusion

From the presented numerical results, we may deduce the following. For the same value of π”βˆ— it is found that the unstable domains are decreasing with increasing of 𝐑 values. This means that the influence of magnetic field has a stabilizing effect for all short and long wavelengths. For the same value of 𝐑, it is found that the unstable domains are increasing with increasing of π”βˆ— values. This means that the streaming is strongly destabilizing.


  1. B. J. W. S. Rayleigh, The Theory of Sound, Dover, New York, NY, USA, 2nd edition, 1945.
  2. S. Chandrasekhar, Hydrodynamic and Hydromagnetic Stability, Dover, New York, NY, USA, 1981.
  3. P. H. Robert, An Introduction to MHD, Longman, London, UK, 1967.
  4. L. Chenng, β€œInstability of a gas jet in liquid,” Physics of Fluids, vol. 28, article 2614, 3 pages, 1985. View at Publisher Β· View at Google Scholar
  5. J. M. Kendall, β€œExperiments on annular liquid jet instability and on the formation of liquid shells,” Physics of Fluids, vol. 29, no. 7, pp. 2086–2094, 1986. View at Publisher Β· View at Google Scholar
  6. P. G. Drazin and W. H. Reid, Hydrodynamic Stability, Cambridge Mathematical Library, Cambridge University Press, Cambridge, Mass, USA, 2nd edition, 2004. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  7. A. E. Radwan, β€œEffect of magnetic fields on the capillary instability of an annular liquid jet,” Journal of Magnetism and Magnetic Materials, vol. 72, no. 2, pp. 219–232, 1988. View at Publisher Β· View at Google Scholar
  8. A. E. Radwan, M. A. Elogail, and N. E. Elazab, β€œLarge hydromagnetic axisymmetric instability of a streaming gas cylinder surrounded by bounded fluid with non uniform field,” Kyungpook Mathematical Journal, vol. 47, no. 4, pp. 455–471, 2007. View at Google Scholar
  9. A. E. Radwan and A. A. Hasan, β€œAxisymmetric electrogravitational stability of fluid cylinder ambient with transverse varying oscillating field,” IAENG International Journal of Applied Mathematics, vol. 38, no. 3, pp. 113–120, 2008. View at Google Scholar
  10. A. E. Radwan and A. A. Hasan, β€œMagneto hydrodynamic stability of self-gravitational fluid cylinder,” Applied Mathematical Modelling, vol. 33, no. 4, pp. 2121–2131, 2009. View at Publisher Β· View at Google Scholar Β· View at Zentralblatt MATH
  11. A. A. Hasan, β€œElectrogravitational stability of oscillating streaming fluid cylinder,” Physica B, vol. 406, no. 2, pp. 234–240, 2011. View at Publisher Β· View at Google Scholar
  12. M. Abramowitz and I. Stegun, Handbook of Mathematical Functions, Dover, New York, NY, USA, 1970.