# Axisymmetric Magneto Dynamic (MHD) Stability of a Compressible Fluid Cylinder

^{*}

**Corresponding Author:**Barakat HM, Department of Mathematics, Faculty of Science and arts, Al Jouf University, Al Jouf, The Kingdom of Saudi Arabia, Tel: 966 14 624 7493, Email: [email protected]

*
Received Date: Sep 11, 2017 /
Accepted Date: Oct 13, 2017 /
Published Date: Oct 27, 2017 *

### Abstract

The axisymmetric magneto dynamic (MHD) stability of a compressible fluid cylinder under the action of inertia, and electromagnetic forces is developed. A general eigenvalue relation is derived studied analytically and the results are confirmed numerically. In absence of the effect of the electromagnetic forces interior and exterior the fluid, so the model is only subjected to the capillary force. It is found that the model is unstable in the region 0<x<1. While it is stable in the region 1<x<∞. This means that model is just unstable in small domains of axisymmetric perturbation but it stable in all other domains. For very high intensity of magnetic field the model is completely stable for all values of wave lengths. The compressibility has a stabilizing tendency.

**Keywords:**
Hydro magnetic stability; Compressible fluid cylinder; Capillary force; Magneto hydrodynamic

#### Introduction

The stability of a liquid column has been studied; the authors derived the dispersion relation. The examined the type of perturbation on the boundary of the capillary instability of liquid jet [1-4]. Such works have been extended [5-14]. The present work is different from those studied before since the velocity is not solenoidal any more (i.e., ∇*u* 6=0) and also that the density of the fluid is not uniform. In this chapter we study the stability of a full fluid cylinder of radius R0 endowed with surface tension and pervaded by axial magnetic field for all modes of perturbation. The fluid is assumed to be compressible, inviscid and perfectly conducting. The fluid (of density ρ) is pervaded by the magnetic field

(1)

The surrounding tenuous medium around the cylinder is assumed to be pervaded by the magnetic field

(2)

where α is parameter that must satisfy some physical conditions. We shall use the cylinder coordinates (*r*, *?*, *z*) system with the z-axis coinciding with the axis of the cylinder. The fluid is acting upon the capillary, gradient pressure and electromagnetic forces. The capillary force which is due to the curvature pressure is acting along the fluidtenuous interface.

#### Formulation of the Problem

The basic equations which required for studying the problem under consideration are the combination of the ordinary hydrodynamic equations together with those of Maxwell concerning the electromagnetic theory. In addition we have to add the equation concerning the curvature pressure.

Under the present circumstances, the MHD basic equations in the fluid are given be

(3)

(4)

(5)

(6)

(7)

In the vacuum region, the basic equations are

(8)

(9)

Along the fluid vacuum interface, we have

(10)

Here and *P* are the fluid velocity vector and kinetic pressure, μ and are the magnetic field permeability coefficient and intensity of magnetic field, respectively, is the magnetic field intensity in the vacuum region, P_{s} is the curvature pressure, γ and *C* are constants, S is the surface tension coefficient and is a unit outward vector normal to the cylinder interface and points like the cylinder coordinate, does.

The unit normal vector is given by

(11)

Where *F* (*γ* ,*?*, *z*,*t* ) = 0 (12)

is the equation of the boundary surface.

#### Unperturbed State

The Unperturbed state is studied by considering the basic equations system in eqns.(3) and (10) by taking into account eqns. (1) and (2), the hydromagnetics equation (3)yields

(13)

from which, we get

(14)

In order to identify this constant of integration we have to apply the balance of the pressure across the boundary surface at r=R_{0}. Taking into account that eqn. (10) in the initial state gives

(15)

Consequently, the unperturbed pressure distribution is finally given by

(16)

where the first term in the right side eqn. (16) represents the capillary force contribution while is due to the electromagnetic forces influence internal and external the fluid cylinder.

In the unperturbed state the pressure P_{0} must be positive and so, using eqn. (16), we find that α satisfies the restriction

(17)

where the equality is pertaining to the limiting case of zero fluid pressure.

#### Perturbation Analysis

We consider a small disturbance to the basic state. A perturbed quantity Q(r,0,z;t) may be expressed as

Q=Q_{0+}Q_{1+}.......,|Q_{1}| << Q_{0} (18)

with,

Q_{1}(r,0,z,t)=Q(r)δ(t)exp(i(kz)), (19)

where Q(r,0,z,t) stands for and while the subscripts (0)and (1) are pertaining to the unperturbed and perturbed quantities respectively and δ(0) is the initial amplitude of perturbation.

*δ(t)=δ _{0}exp(σt)*, (20)

and hence,

Q_{1}(r,z;t)=Q(r)δ_{0} exp(σt+ikz), (21)

where σ is the growth rate. If is imaginary then is the wave oscillation frequency. The perturbed radial distance of the fluid cylinder may be expressed as

r=R_{0}+R_{1}, R_{1}=δ_{0} exp(*ikz*+*σt*) (22)

Here *R*_{1} is the elevation of the surface wave measured from the unperturbed position, k is the longitudinal wave number and σ is the growth rate. By the use of the expansion (22), the perturbed equations are given by

(23)

(24)

(25)

(26)

(27)

while those equations in the region surrounding the fluid cylinder are

(28)

(29)

and along the fluid-vacuum interface we have

(30)

In view of the expansion in eqns. (21) and (22), the relevant perturbation eqns. (23-30) may be simplified and solved. Under the present circumstances the non-singular solution is given by

(31)

(32)

(33)

(34)

(35)

(36)

(37)

By substituting from eqn. (32) into eqn. (31), we get

(38)

with

(39)

By using eqn. (38) taking into account the dependence in eqn. (22), the solution of eqn. (38) is given in terms of Bessel functions in the form

(40)

where I_{0}(ηr) is the modified Bessel function of the first kind of order zero and A is an arbitrary constant to be determined .

The perturbed magnetic field in the region surrounding the fluid cylinder is obtained by solving the relevant perturbation eqns. (28) and (29).

Eqn. (28) means that the perturbed magnetic field in the region surrounding the field can be derived from scalar function ψ_{1}, say

(41)

By combining eqns. (28) and (41), we get

(42)

By using expansion in eqn. (21) and taking into account the space dependence in eqn. (22), the solution in eqn. (42) is given in terms of Bessel functions. Under the present circumstances, we have

(43)

and consequently

(44)

where K_{0}(kr) is the modified Bessel function of the second kind of order zero and B is an arbitrary constant. Finally inserting eqn. (22) into eqn. (30), yields

(45)

#### Boundary Conditions

The solution gives by eqns. (31,32,33, 34 and 45) must satisfy some boundary conditions appropriate to the problem under consideration. these conditions may be written as follows.

1. The normal component of the velocity vector u must be compatible with the velocity of the fluid particles across the boundary surface at r=R_{0}.

(46)

By use of the eqns. (22) and (32) for the conditions in eqn. (46), we get

(47)

where y=ηR_{0} is the dimensionless longitudinal wave number. then, we get

(48)

from which

(49)

with

(50)

2. The normal component of the magnetic field H must be continuous across the perturbed boundary (22) at r=R_{0}. This condition may be written as

(51)

The unit vector n normal to the cylinder interface and the magnetic fields H and H^{vac} are given by

(52)

then

n=(1,0,0)+(0,0,−ik)R_{1} (53)

(54)

By substituting from eqns. (33) and (44) into eqn. (52), we get

(55)

Therefore, the perturbed magnetic field external the fluid cylinder is given by

(56)

3. The normal component of the stresses due to the kinetic pressure of the fluid and the magnetic pressure of electromagnetic forces acting inside and outside column must be discontinuous by the curvature pressure. This may be written as

(57)

Therefore, the condition in eqn. (57) may be rewritten in the form

(58)

Substituting from eqns. (35, 44 and 45) into eqn. (58), after lengthy calculations, we finally arrive to

(59)

This dispersion relation is valid for all axisymmetric m=0 and nonaxisymmetric (m ≥ 1) modes of perturbation.

#### General Discussions

Eqn. (59) is the capillary dispersion relation of a fluid cylinder acting upon Lorentz force with a homogeneous magnetic field. It is a simple linear combination of the dispersion relation of a fluid cylinder subject to the capillary force only and that one subject upon electromagnetic force with uniform magnetic fields.

The relation in eqn. (59) is valid for all axisymmetric and nonaxisymmetric modes of perturbation. By means of that relation the hydrodynamic, magnetodynamic and magnetohydrodynamic (MHD) characteristics of the present model may be determined. The eigenvalue relation in eqn. (59) relates the growth rate σ, or rather the oscillation frequency ω as σ=iω is imaginary with transverse wave number.

Some reported works may be obtained from relation in eqn. (59) under appropriate choices.

If we impose H_{0}=0, eqn. (59) reduce to

(60)

As the fluid is incompressible as a → ∞ and in such case y → x, therefore eqn. (59) degenerates to

(61)

If we impose H_{0}=0 eqn. (61) reduce to that derived by Rayleigh [15]

(62)

#### Hydrodynamic Instability

If the model under consideration is acting upon the capillary force only and all other forces are neglected, the stability criterion for such a case is given by the relation in eqn.(60). The discussions of this relation showed that the model is capillary stable in the nonaxisymmetric modes for all short and long wave lengths. In the axisymmetric mode the incompressible model is stable as long as x ≥ 1 and only unstable if 0<x<1 where x=1 corresponding to the marginally stability state.

It is deduced that the fluid cylinder is capillary unstable only in the axisymmetric mode as long as the perturbed wave length is longer than the circumference of the fluid jet.

#### Magnetodynamic Stability

In the absence of the capillary force and the fluid cylinder is only subjected to the electromagnetic forces interior and exterior the fluid, the dispersion relation in eqn. (59) yields

(63)

The influence of the longitudinal magnetic field interior the fluid in eqn. (1) is represented by the term −x^{2} following the fundamental quantity in eqn. (63).

It has a strong stabilizing influence for all short and long wave lengths.

Now we prove that the exterior magnetic field in eqn. (2) is also stabilizing.

By the use of the recurrence relations of the modified Bessel functions

(64)

and the fact that

I_{m}(y) > 0 and K_{m}(y) > 0 (65)

For every non-zero value of x and y, we may prove that

(66)

For all x≠ 0, y≠0 and m ≥ 0.

From the point of view of the relation in eqn. (64) and the inequalities in eqn. (66) it can be shown that the magnetic field acting exterior the fluid cylinder is stabilizing. That character is valid not only in the non-axisymmetric modes but also in the axisymmetric mode for all short and long wave lengths.

We conclude that the electromagnetic forces acting interior and exterior the fluid cylinder are stabilizing. This may be due to the fact that the applicable magnetic fields are axial and uniform.

#### Magneto Hydrodynamic Stability

In order to discuss this general case in which the fluid cylinder is acting upon the capillary and electromagnetic forces we have to discuss the relation in eqn. (59) in its general form we have the following

(i) The capillary force is stabilizing if 0 ≤ x<∞, 1 ≤ y<∞ where the equality corresponds to the neutral stability, while it is destabilizing as long as 0 ≤ x<1, 0 ≤ y<1 (ii) The lorentz force is stabilizing for all values of x 6=0, y 6=0 i.e., for all short and long wave lengths.

We deduce that the model under consideration is MHD stabilizing for small wave length.

#### Numerical Discussions

In order to verify the results which obtained analytical concerning the acting different forces effect on the present model. It is found very important to discuss the dispersion relation in eqn. (59) numerically for axisymmetric perturbation mode. In order to do that we have to rewrite this dispersion relation in non-dimensional form, so we may insert this relation in the computer for making the numerical computation.

Based on the input data, if the values of are positive then we have unstable regions. As values are negative we put *ó* = *iw* so are the values which concerning the stable domains. In the transition from the negative values to positive values of σ2 we have to path with the values of σ=0 which means marginal stability states. The points at which the transition from stability regions to those of instability called the critical points.

The dispersion relation in eqn. (59) has been formulated in the dimensionless form upon using the quantity which has a unit of times as and that the quantity say Hs has a unit of the intensity of magnetic field. Consequently, equation in eqn. (58) takes the form

(67)

The numerical calculations have been performed by inserting the dispersion relation in eqn. (67) in the computer and computed for the important mode of perturbation. Taking into account that and that . The calculations have been carried out for all short and long wave lengths as 0.1 ≤ x = 0,0.2,0.4,0.6,0.8,1.0and 2.0 and 2.0.

For each values of , the calculations performed for several values of (α,a)=(1,5), (2,5), (3,5), (4,5), (1,30), (1,20),(1,10), (2,20), (3,10) and (4,10).

In each case the data are collected, tabulated in tables and presented in figures, (**Figures 1-4**). In each figures and table it is printed clearly the stable and unstable domains corresponding to a given value of and several values of (α,a).

For =0.0, we obtain the same unstable domain which is due to the capillary force. Whatever are the values of (α,a) we found the unstable is 0<xc=1. Not that as H_{0}=0 then the values of (α,a) has no effect at all (**Figure 1**).

For =0.4 corresponding to (α,a)=(1,30), (2,20), (3,10) and (4,5), it is found that the unstable domain are given by 0.1<x<0.718, 0.1<xc=0.718, 0.654, 0.481, 0.356 (**Figure 2**).

For =1.0 , corresponding to (α,a)=(1,30), (2,20), (3,10) and (4,5), it is found that there is no any unstable domain. The model is completely stable (**Figure 3**).

For =2.0, corresponding to (α,a)=(1,30), (2,20), (3,10) and (4,5), it is found that there is no any unstable domain. The model is completely stable (**Figure 4**).

#### Conclusion

From the forgoing discussions, we may conclude the following results.

1. In absence of the effect of the electromagnetic forces interior and exterior the fluid, so the model is only subjected to the capillary force. It is found that the model is unstable in the region 0<x<1. While it is stable in the region 1<x<∞. This means that model is just unstable in small domains of axisymmetric perturbation but it stable in all other domains.

2. For very high intensity of magnetic field the model is completely stable for all values of wave lengths.

3. The compressibility has a stabilizing tendency.

#### References

- Abramowitz M, Stgun IA (1970) Handbook of mathematical functions: with Formulas, Graphs, and Mathematical Tables. Dover Publisher.
- Donnelly RJ, Glaberson W (1966) Experiments on the capillary instability of a liquid jet. Proceedings of the Royal Society 290: 547.
- Drazin PG (1970) Kelvin-Helmholtz Instability of finite Amplitude. Journal Fluid Mechanics 42: 321-335.
- Drazin GP, Reid WH (1980) Hydrodynamic Stability. Cambridge University Press, p: 605.
- Drazin PG (2002) Introduction to Hydrodynamic Stability. Cambridge University Press.
- Korovin VM (2001) The capillary instability of a thin cylinder of viscous Ferro fluid in a longitudinal magnetic field. Journal of Applied Mathematics and Mechanics 65: 243-251.
- Radwan AE, Callebaut DK (1986) Magnetohydrodynamic stability of a hollow jet. Proceedings of the Physical Society 21: 32-49.
- Radwan AE, Elazab SS (1988) Magnetohydrodynamic stability of an infinite self-gravitating. Journal of the Physical Society of Japan 57: 461-476.
- Radwan AE (1989) Non-axisymmeric magnetohydrodynamic stability of a streaming annular jet containing very dense fluid cylinder and surrounded by a tenuous medium pervaded by a varying azimuthal magnetic field. Physical Society of Japan 58: 1225-1239.
- Radwan AE, Hoda AH (2004) Magnetohydrodynamic instability of KelvinHelmholtz with annulus cylinderical interface pervaded by helical varying magnetic fields. Applied Mathematics and computation 147: 191-210.
- Radwan AE (2005) Instability of a compressibile annular magnetized fluid. Applied Mathematics and Computation 160: 283-301.
- Radwan AE (2007) Electrohydrodynamic stability of a dielectric fluid cylinder. Physica Scripta 75: 510-526.
- Radwan AE, Hasan AA (2008) Magneto hydrodynamic stability of self-gravitational fluid cylinder. Applied Mathematical Modeling 38: 113-124.
- Radwan AE, Hasan AA (2009) Axisymmetric electro-gravitational stability of fluid cylinder ambient with transverse varying oscillating field. IAENG International Journal of Applied Mathematics 33: 2121-2132.
- Rayleigh JM (1945) The Theory of sound. Dover Publication.

Citation: Barakat HM (2017) Axisymmetric Magneto Dynamic (MHD) Stability of a Compressible Fluid Cylinder. J Appl Computat Math 6: 372. DOI: 10.4172/2168-9679.1000372

Copyright: © 2017 Barakat HM. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.