SpringerOpen Newsletter

Receive periodic news and updates relating to SpringerOpen.

Open Access Research

Electrogravitational stability of oscillating streaming fluid cylinder ambient with a transverse varying electric field

Alfaisal A Hasan

Author Affiliations

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

Boundary Value Problems 2011, 2011:31  doi:10.1186/1687-2770-2011-31


The electronic version of this article is the complete one and can be found online at: http://www.boundaryvalueproblems.com/content/2011/1/31


Received:29 May 2011
Accepted:11 October 2011
Published:11 October 2011

© 2011 Hasan; licensee Springer.

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

Abstract

The electrogravitational instability of a dielectric oscillating streaming fluid cylinder surrounded by tenuous medium of negligible motion pervaded by transverse varying electric field has been investigated for all the perturbation modes. The model is governed by Mathieu second-order integro-differential equation. Some limiting cases are recovering from the present general one. The self-gravitating force is destabilizing only in the axisymmetric perturbation for long wavelengths, while, the axial electric field interior, the fluid has strong destabilizing effect for all short and long wavelengths. The transverse field is strongly stabilizing. In the case of non-axisymmetric perturbation, the self-gravitating force is stabilizing for short and long waves, while the electric field has stabilizing effect on short waves.

Keywords:
electrogravitational stability; oscillating; streaming

1. Introduction

The stability of self-gravitating fluid cylinder has been studied, for the first time, by Chandrasekhar and Fermi [1]. Later on, Chandrasekhar [2] made several extensions as the fluid cylinder is acted by different forces. Radwan [3,4] studied the stability of an ideal hollow jet. Radwan [4] considered that the fluids are penetrated by constant and uniform electric fields. 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 [5,6]. Radwan and Hasan [5] studied the gravitational stability of a fluid cylinder under transverse time-dependent electric field for axisymmetric perturbations. Hasan [7,8] 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 non-axisymmetric perturbation modes. Hasan [7,8] 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 the modes of perturbations.

There are many applications of electrohydrodynamic and magnetohydrodynamic stability in several fields of science such as

1. 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 because of the motion of the liquid iron.

2. 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, etc. Many astrophysical systems are not in local thermal equilibrium, and therefore require an additional kinematic treatment to describe all the phenomena within the system.

3. Engineering applications: there are many forms in engineering sciences including oil and gas extraction process if it surrounded by electric field or magnetic field, gas and steam turbines, MHD power generation systems and magneto-flow meters, etc.

In this article, we aim to investigate 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 the axisymmetric and non-axisymmetric perturbation modes.

2. Mathematical formulation

Consider a self-gravitating fluid cylinder surrounded by a self-gravitating medium of negligible motion. The cylinder of (radius R0) dielectric constant ε(i) while the surrounding medium is being with dielectric constant ε(e). Fluid is assumed to be incompressible, inviscid, self-gravitating, and pervaded by applied longitudinal electric field.

E ¯ 0 ( i ) = ( 0 , 0 , E 0 ) (1)

The surrounding tenuous medium (being of negligible motion), self-gravitating, and penetrated by transverse varying electric field

E ¯ 0 ( e ) = 0 , β E 0 R 0 r - 1 , 0 (2)

where E0 is the intensity of the electric field in the fluid while β is some parameters satisfy certain conditions. The components of E ¯ 0 ( i ) and E ¯ 0 ( e ) are considered along the utilizing cylindrical coordinates (r, φ, z) system with z-axis coinciding with the axis of the fluid cylinder. The fluid of the cylinder streams with a periodic velocity

u ¯ 0 = 0 , 0 , U cos ω t (3)

where ω is constant and U is the speed at time t = 0.

The components of electric fields E ¯ 0 ( i ) and E ¯ 0 ( e ) are being along (r,φ,z) with the z-axis coinciding with the axis of the fluid cylinder (as shown in Figure 1).

thumbnailFigure 1. Sketch for gravitational dielectric fluid cylinder.

The basic equations for investigating the problem under consideration are being the combination of the ordinary hydrodynamic equations, Maxwell equations concerning the electromagnetic theory, and Newtonian self-gravitating equations concerning the self-gravitating matter (see [2,7-10]).

For the problem under consideration, these equations are given as follows.

ρ u ¯ t + u ¯ u ¯ ( i ) = - P ( i ) + ρ V ( i ) + 1 2 ε i E ¯ i E ¯ i (4)

u ¯ ( i ) = 0 (5)

ε E ¯ i , e = 0 (6)

ε i , e E ¯ i , e = 0 (7)

2 V ( i ) = - 4 π ρ G (8)

2 V ( e ) = 0 (9)

where ρ, u ¯ , and P are the fluid density, velocity vector, and kinetic pressure, respectively, and E ¯ i and V ¯ i are the electric field intensity and self-gravitating potential of the fluid while E ¯ e and V ¯ e are these of tenuous medium surrounding the fluid cylinder, and G is the gravitational constant.

Since the motion of the fluid is irrotational, incompressible motion, the fundamental equations may be written as

2 ϕ ( i ) = 0 (10)

2 ψ ( i ) = 0 (11)

2 ψ ( e ) = 0 (12)

where ϕ and ψ are the potential of the velocity of the fluid and electrical potential.

3. Equilibrium state

In this case, the basic equations are given in the form

2 V 0 ( i ) = - 4 π ρ G (13)

2 V 0 ( e ) = 0 (14)

2 ϕ 0 ( i ) = 0 (15)

where the subscript 0 here and henceforth indicates unperturbed quantities.

Equations 12-14 are solved and moreover the solutions are matched across the fluid cylinder interface at r = R0. The non-singular solution in the unperturbed state is, finally, given as

V 0 ( i ) = - π G ρ r 2 (16)

V 0 ( e ) = - π G ρ R 0 2 1 + 2 ln r R 0 (17)

4. Linearization

For a small wave disturbance across the boundary interface of the fluid, the surface deflection at time t is assumed to be of the form as

r = R 0 + η ̃ , (18)

with

η ̃ = η t exp i k z + m φ (19)

Consequently, any physical quantity Q(r,φ,z;t) may be expressed as

Q r , φ , z , t = Q 0 r + η ̃ z , φ , t (20)

where η(t) is the amplitude of the perturbation at an instant time t, k, any real number, is the longitudinal wave number along z-direction while m, an integer, is the azimuthal wave number.

The non-singular solutions of the linearized perturbation equations give ϕ,V, and ψ as follows:

ϕ 1 ( i ) = A 1 t I m k r exp i k z + m φ , (21)

V 1 ( i ) = B 1 t I m k r exp i k z + m φ , (22)

V 1 ( e ) = B 2 t K m k r exp i k z + m φ (23)

ψ 1 ( i ) = C 1 t I m k r exp i k z + m φ (24)

ψ 1 ( e ) = C 2 t K m k r exp i k z + m φ , (25)

where A1(t), B1(t), B2(t), C1(t), and C2(t) are arbitrary functions of integrations to be determined, while Im(kr) and Km(kr) are the modified Bessel functions of the first and second kind of order m.

5. Boundary conditions

The non-singular solutions of the linearized perturbation equation given by the systems (21)-(25) and the solutions (16)-(17) of the unperturbed systems (12)-(14) must satisfy certain boundary conditions. Under the present circumstances, these appropriate boundary conditions could be applied as follows.

(i) Kinematic conditions

The normal component of the velocity vector must be compatible with the velocity of the boundary perturbed surface of the fluid at the level r = R0. This condition, yield

( t + U cos ω t z ) η ˜ = ϕ 1 ( i ) r (26)

By the use of Equations 18, 19, and 21 for the condition (26), after straight forward calculations, we get

A 1 t = 1 k I m ( x ) t + i k U η (27)

where x = k R0 is, dimensionless, the longitudinal wave number.

(ii) Self-gravitating conditions

The gravitational potential V = V0 + εV1 + ⋯ and its derivative must be continuous across the perturbed boundary fluid surface at r = R0. These conditions are given as

V 1 + η ̃ V 0 r ( i ) = V 1 + η ̃ V 0 r ( e ) , (28)

V 1 r + η ̃ 2 V 0 r 2 ( i ) = V 1 r + η ̃ 2 V 0 r 2 ( e ) . (29)

By utilizing Equations 18, 19, 22, and 23 for the conditions (28) and (29), we get

B 1 t = 4 π G k ρ x K m x η (30)

B 2 t = 4 π G k ρ x I m x η (31)

(iii) Electrodynamic condition

The normal component of the electric displacement current and the electric potential ψ perturbed boundary surface at the initial position r = R0. These conditions could be written in the form

ψ 1 + η ̃ ψ 0 r ( i ) = ψ 1 + η ̃ ψ 0 r ( e ) (32)

N ¯ . ( ε ( i ) E ¯ ( i ) - ε ( e ) E ¯ ( e ) ) = 0 (33)

E ¯ = E ¯ 0 + η E ¯ 0 r + E ¯ 1 (34)

While N ¯ s is, the outward unit vector normal to the interface (18) at r = R0, given by

N ¯ s = F ( r , φ , z ; t ) / | F ( r , φ , z ; t ) | (35)

F r , φ , z ; t = r - R o - η ̃ (36)

So that

N ¯ 0 = 1 , 0 , 0 , N ¯ 1 = 0 , - i m R 0 , - i k η ̃ (37)

Upon applying these conditions, we get

C 1 t = - i E 0 ε ( i ) η ξ 1 1 + m β - m β R 0 (38)

C 2 t = - i E 0 ε ( i ) η ξ 1 I m x K m x 1 + m β - m β R 0 (39)

where the quantity ξ1 is given in Appendix 1.

(iv) The dynamical stress condition

The normal component of the total stress across the surface of the coaxial fluid cylinder must be continuous at the initial position at r = R0. This condition is given as follows

ρ ϕ 1 ( i ) t + U 0 ϕ 1 ( i ) z - V 1 ( i ) - η ̃ V 0 ( i ) r + E 0 ε i ψ 1 ( i ) z - ε e 1 R 0 r β R 0 r ψ 1 e - η ̃ r E 0 E 0 e = 0 (40)

By substituting for ϕ 1 i , V 1 i , V o i , ψ 1 i , ψ 1 e and η ̃ , after some algebraic calculations, we finally obtain

d 2 η d t 2 + 2 i k U o cos ω t d η d t + G β 1 1 - i k ω U o sin ω t - k 2 U o 2 cos 2 ω t + E o 2 β 1 2 η = 0 (41)

where the quantity β11 and β12 is given in Appendix I.

In order to eliminate the first derivative term, we may use the substitution

η t = η * t e - i k U 0 ω sin ω t (42)

Equation 41 can be expressed as follows

d 2 η * d t 2 + G β 1 1 + E 0 2 β 1 2 η * = 0 (43)

Equation 43 is an integro-differential equation governing the surface displacement η*(t). By means of this relation, we may identify the (in-) stability states and also the self-gravitating and electrodynamic forces influences on the stability of the present model. However in order to do so, it is found more convenient to express this relation in the simple form

d 2 d γ 2 + b - h 2 cos 2 γ η * t = 0 , γ = ω t (44)

where

b = G β 1 1 2 (45)

h 2 = - E 0 2 β 1 2 ω 2 (46)

Equation 44 has the canonical form

d 2 d γ 2 + a - 2 q cos 2 γ η * t = 0 (47)

where

q = h 2 4 , a = b - h 2 2 (48)

Equation 47 is Mathieu differential equation. The properties of the Mathieu functions are explained and investigated by Melaclan [11]. The solutions of Equation 47, under appropriate restrictions, could be stable and vice versa. The conditions required for periodicity of Mathieu functions are mainly dependent on the correlation between the parameters a and q. However, it is well known, see [11], that (a, q)-plane is divided essentially into two stable and unstable domains separated by the characteristic curves of Mathieu functions. Thence, we can state generally that a solution of Mathieu integro-differential equation is unstable if the point (a, q) say, in the (a, q)-plane lies internal and unstable domain, otherwise it is stable.

6. Discussions and limiting cases

The appropriate solutions of Equation 47 are given in terms of what called ordinary Mathieu functions which, indeed, are periodic in time t with period π and 2π.

Corresponding to extremely small values of q, the first region of instability is bounded by the curves

a = ± q + 1 (49)

The conditions for oscillation lead to the problem of the boundary regions of Mathieu functions where Melaclan [11] gives the condition of stability as

| Δ ( 0 ) sin 2 ( π a 2 ) | 1 2 1 (50)

where Δ(0) is the Hill's determinant.

An approximation criterion for the stability near the neighborhood of the first stable domains of the Mathieu stability domains given by Morse and Feshbach [12] which is valid only for small values of h2 or q, i.e., the frequency ω of the electric field is very large.

This criterion, under the present circumstances, states that the model is ordinary stable if the restriction

h 4 - 1 6 1 - b h 2 + 3 2 b 1 - b 0 (51)

is satisfied where the equality is corresponding to the marginal stability state. The inequality (51) is a quadratic relation in h2 and could be written as

h 2 - α 1 h 2 - α 2 0 (52)

where α1 and α2 are, the two roots of the equality of the relation (51), being

α 1 = 8 1 - b - Δ (53)

α 2 = 8 1 - b + Δ (54)

with

Δ 2 = 3 2 1 - b 2 - 3 b (55)

The electrogravitational stability and instability investigations analysis should be carried out in the following two cases

(i). 0 < b < 2/3

In this case Δ2 is positive and therefore the two roots α1 and α2 of the equality (51) are real. Now, we will show that both α1 and α2 are positive. If α1 α + ve then α1 must be negative and this means that

8 1 - b b (56)

or alternatively

6 4 1 - b 2 3 2 1 - b 2 - 3 b

From which we get

2 b 3 b (57)

and this is contradiction, so α1 must be positive and consequently α2 ≥ 0 as well (noting that α2 > α1). This means that both the quantities (h2 -α1) and (h2 -α2) are negative and that in turn show that the inequality (51) is identically satisfied.

(ii). 2/3 < b < 1

In this case, in which b < 1 and simultaneously 3b > 2, it is found that Δ2 is negative, i.e., Δ is imaginary; therefore, the two roots α1 and α2 are complex. We may prove that the inequality (51) is satisfied as follows.

Let h2 - c and α1,2 = c1 - ic2 where c, c1, and c2 are real, so

h 2 - α 1 h 2 - α 2 = - c - c 1 + i c 2 - c - c 1 - i c 2 = c 2 + 2 c c 2 + c 1 2 + c 2 2 = c + c 1 2 + c 2 2 = + ve (58)

which is positive definite.

By an appeal to the cases (i) and (ii), we deduce that the model is stable under the restrictions

0 < b < 1 (59)

This means that the model is stable if there exists a critical value ω0 of the electric field frequency ω such that ω > ω0 where ω0 is given by

π G ρ i x I 0 x I 0 x I 0 x K 0 x - 1 2 > 0 (60)

One has to mention here that if ω = 0, β = 0, and E0 = 0 and we suppose that

γ t = const exp σ t (61)

The second-order integro-differential equation of Mathieu equation (41) yields

σ 2 = 4 π G ρ i x I 0 x I 0 x I 0 x K 0 x - 1 2 (62)

where σ is the temporal amplification and note by the way that 4 π G ρ i - 1 2 has a unit of time. The relation (62) is identical to the gravitational dispersion relation derived for the first time by Chandrasekhar and Fermi [1]. In fact, they [1] have used a totally different technique rather than that used here. They have used the method of representing the solenoidal vectors in terms of poloidal and toroidal vector fields for axisymmetric perturbation.

To determine the effect of ω, it is found more convenient to investigate the eigenvalue relation (62) since the right side of it is the same the middle side of (60).

Taking into account the recurrence relation of the modified Bessel's functions and their derivatives, we see, for x α 0, that

x I 0 x I 0 x > 0 (63)

and

I 0 x K 0 x > 1 2 ,  or  I 0 x K 0 x < 1 2 (64)

based on the values of x.

Now, returning to the relation (62), we deduce that the determining of the sign σ2/(4πGρi) is identified if the sign of the quantity

Q o x = I 0 x K 0 x - 1 2 (65)

is identified.

Here, it is found that the quantity Q0 (x) may be positive or negative depending on x α 0 values. Numerical investigations and analysis of the relation (62) reveal that σ2 is positive for small values of x while it is negative in all other values of x. In more details, it is unstable in the domain 0 < x < 1.0667 while it is stable in the domains 1.0667 ≤ x < ∞ where the equality is corresponding to the marginal stability state.

From the foregoing discussions, investigations, and analysis, we conclude (on using (65) for (62)) that the quantity

L 2 = x I 0 x I 0 x I 0 x K 0 x - 1 2 , L = σ 4 π G ρ 1 2 (66)

has the following properties

L 2 0 in the ranges 1.0667 x < L 2 > 0 in the range 0 < x < 1.0667 } (67)

Now, returning to the relation (60) concerning the frequency ω0 of the periodic electric field

ω 2 4 π G ρ > x I 0 x I 0 x 1 2 - I 0 x K 0 x > 0 . (68)

Therefore, we deduce that the electrodynamic force (with a periodic time electric field) has stabilizing influence and could predominate and overcoming the self-gravitating destabilizing influence of the dielectric fluid cylinder dispersed in a dielectric medium of negligible motion.

However, the self-gravitating destabilizing influence could not be suppressed whatever is the greatest value of the magnitude and frequency of the periodic electric field because the gravitational destabilizing influence will persist.

7. Numerical discussions

If we assume that ω = 0 and consider the condition (61), then the second-order integro-differential equation of Mathieu equation (47) yields

σ 2 4 π G ρ = x I 0 x I 0 x I 0 x K 0 x - 1 2 - M x I 0 x I 0 x x I 0 x K 0 x I 0 x K 0 x - ε I 0 x K 0 x - ε e β 2 = 0 (69)

where

M = E 0 E s 2 , E s 2 = 4 π G ( ρ ) 2 R 0 2 ε i (70)

and

ε = ( ε ( e ) / ε ( i ) ) (71)

To verify and confirm the foregoing analytical results, the relation (69) has been inserted in the computer and computed. This has been done for several values of β as β < 1, β = 1, and β > 1 in the wide domain 0 ≤ x ≤ 0.5. The numerical data of instability corresponding σ / ( 4 π G ρ i ) 1 2 and those of stability corresponding to ζ / ( 4 π G ρ i ) 1 2 are collected and tabulated and presented graphically (see Figures 2, 3, 4, 5, and 6). There are many features and properties in this numerical presentation as we see in the following:

thumbnailFigure 2. Electrogravitational stable and unstable domains for β = 0.5.

thumbnailFigure 3. Electrogravitational stable and unstable domains for β = 1.0.

thumbnailFigure 4. Electrogravitational stable and unstable domains for β = 1.5.

thumbnailFigure 5. Electrogravitational stable domains for β = 2.5.

thumbnailFigure 6. Electrogravitational stable domains for β = 3.0.

(i) For β = 0.5 corresponding to M = 0.1, 0.3, 0.5, 0.7, 1.0, and 1.5 it is found that the electrogravitational unstable domains are 0 < x < 1.1175, 0 < x <1.19759, 0 < x < 1.27235, 0 < x 1.29599, 0 < x < 1.362741, and 0 < x < 1.3978, the neighboring stable domains are 1.1175 ≤ x < ∞, 1.19759 ≤ x < ∞, 1.27235 ≤ x < ∞, 1.29599 ≤ x < ∞, 1.362741 ≤ x < ∞, and 1.3978 ≤ x < ∞, where the equalities correspond to the marginal stability states (see Figure 2).

(ii) For β = 1.0 corresponding to M = 0.1, 0.3, 0.5, 0.7, 1.0, and 1.5 it is found that the electrogravitational unstable domains are 0 < x < 1.22669, 0 < x < 1.5266, 0 < x < 1.750969, 0 < x < 1.90513, 0 < x < 2.05422, and 0 < x < 2.19341, the neighboring stable domains are 1.22669 ≤ x < ∞, 1.5266 ≤ x < ∞, 1.750969 ≤ x < ∞, 1.90513 ≤ x < ∞, 2.05422 ≤ x < ∞, and 2.19341 ≤ x < ∞, where the equalities correspond to the marginal stability states (see Figure 3).

(iii) For β = 1.5 corresponding to M = 0.1, 0.3, 0.5, 0.7, 1.0, and 1.5 it is found that the electrogravitational unstable domains are 0 < x < 1.35924, 0 < x < 1.9735, 0 < x < 2.3982, 0 < x < 2.6563, 0 < x < 2.8835, and 0 < x < 3.0798, the neighboring stable domains are 1.35924 ≤ x < ∞, 1.9735 ≤ x < ∞, 2.3982 ≤ x < ∞, 2.6563 ≤ x < ∞, 2.8835 ≤ x < ∞, and 3.0798 ≤ x < ∞, where the equalities correspond to the marginal stability states (see Figure 4).

(iv) For β = 2.5, corresponding to M = 0.1, 0.3, 0.5, 0.7, 1.0, and 1.5 it is found that the electrogravitational fluid cylinder is completely stable not only for short wavelengths, but also for very long wavelengths and the gravitational unstable domains are completely suppressed (see Figure 5).

(v) For β = 3.0, corresponding to M = 0.1, 0.3, 0.5, 0.7, 1.0 and 1.5 it is found that the electrogravitational fluid cylinder is completely stable not only for short wavelengths, but also for very long wavelengths and the gravitational unstable domains are completely suppressed (see Figure 6).

8. Conclusion

From the presented numerical results, we may deduce the following. For the same value of M, it is found that the unstable domains are increasing with increasing of β values. This means that the influence of electric field has a destabilizing effect for all short and long wavelengths.

If β > 2.0, then the model is completely stable not only for short wave lengths, but also for long wave lengths.

Competing interests

The authors declare that they have no competing interests.

Appendix I

ξ 1 = ε ( i ) I m ( x ) K m ( x ) - ε ( e ) I m ( x ) K m ( x ) β 1 1 = 2 π ρ R o k I m ( x ) - 4 π ρ x I m ( x ) K m ( x ) I m x β 1 2 = k 2 ε ( i ) 2 I m ( x ) ξ 1 ρ 1 + m β - m β R o + 2 β 2 k I m ( x ) R o ρ I m ( x ) + i ε ( i ) ε ( e ) I m ( x ) k ξ 1 ρ m β R o 1 + m β - m β R o

Acknowledgements

We are grateful to the Editor of the Journal and the Reviewers for their suggestions and comments on this article.

References

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

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

  3. Radwan, AE: MHD stability of gravitational compressible fluid cylinder. J Phys Soc Jpn. 58, 1225–1227 (1989). Publisher Full Text OpenURL

  4. Radwan, AE: Electrogravitational instability of an annular fluid jet coaxial with a very dense fluid cylinder under radial varying fields. J Plasma phys. 44, 455–465 (1991)

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

  6. Radwan, AE, Hasan, AA: Magnetohydrodynamic stability of selfgravitational fluid cylinder. Appl Math Model. 33(4), 2121–2131 (2009). Publisher Full Text OpenURL

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

  8. Hasan, AA: Electrogravitational stability of oscillating streaming fluid cylinder. J Appl Mech ASME (2011, in press)

  9. Lin, SH: Oxygen diffusion in a spherical cell with nonlinear oxygen uptake kinetics. Theor Biol. 60, 449–457 (1976). Publisher Full Text OpenURL

  10. Mestel, AJ: The electrohydrodynamic cone-jet at high Reynolds number. J Aerosol Sci. 25(6), 1037–1047 (1994). Publisher Full Text OpenURL

  11. Melaclan, N: Theory and Application of Mathieu Function. Oxford University Press, New York (1964)

  12. Morse, PM, Feshbach, H: Methods of Theoretical Physics Part I. McGraw-Hill, New York (1953)