_{1}

^{*}

Two-dimensional thermomagnetic convection of magnetic fluid possessing internal spin and relaxation of magnetization with high thermal sensitivity is numerically analyzed under a non-uniform magnetic and thermal field, using a spectral difference scheme, magnetic fluid is assumed to be placed between concentric cylinders, an azimuthal magnetic field being produced by an electric current though the inner cylinder, which is an adiabatic wall, whereas a half of the outer cylinder is kept at a high constant temperature and the rest half is at a lower constant temperature. As a result, in case of ▽T × ▽|H| ≠ 0, thermomagnetic convection is found to be produced, the convection pattern depends on electric current distribution, depending on which multiple or single circulation is produced, thus, it shows that thermomagnetic convection pattern is controlled by changing relative direction of temperature gradient to that of a magnetic field.

A familiar example of thermal convection is buoyancy-induced convection, where a driving body force resides in a force of gravity, another example is established in a heated magnetic fluid subject to a magnetic body force, which depends on a thermodynamic constitutive equation. The enhancement of convective heat trans- fer in magnetic fluids could be applied in cooling current-carrying conductors, and transmissionlines, especially in gravity-free space. Characteristics of thermomagnetic convection are investigated by many researchers, e.g. Finlayson (1970) [

In this paper, two-dimensional thermomagnetic convection of magnetic fluid possessing internal spin and relaxation of magnetization with high thermal sensitivity is numerically analyzed under a non-uniform magnetic and thermal field, using a spectral finite difference schemes, and the thermomagnetic convection pattern of magnetic fluid within an annulus is numerical investigated with spatial distributions of magnetic field variations.

Analyzed is Two-dimensional thermomagnetic convection followed by a steady- state natural convection of magnetic fluid as shown in _{0} with a direct current I_{0} is placed in the center of an annular space filled with a magnetic fluid of magnetic permeability μ_{0}. It is assumed that magnetic fluid is magnetized by a magnetic field H due to an external magnetic field from a wire current and an induction magnetic field of the magnetic fluid, thosemagnetic fluid is assumed to be at rest up to time t = 0 with a uniform temperature T_{l}, and at t > 0 a half part of the outside circumference is assumed to be fixed at a high temperature (T_{h}), and the rest half part of the outside circumference at a low temperature(T_{l}).

The following assumption applies:

1) The radius of wire is much smaller than that of the inside circle ( R 0 ≪ R 1 ).

2) The effect of magnetic field is limited within the distance 50R_{2}.

3) The magnetic permeability of magnetic fluid is equal to that of vacuum.

4) The magnetic fluid is incompressible fluid.

The equation of continuity is

∇ ⋅ ν = 0 , (1)

The momentum equation neglecting gravity force is

ρ D ν D t = − ∇ p + η ∇ 2 v + ∇ ⋅ ( B H ) + I t s ∇ × ( Ω − ω ) , (2)

where t is time, v: velocity, ρ: mass density of magnetic fluid, p: pressure, η: viscosity of magnetic fluid in absence of a magnetic field, I: average inertia moments of particles per unit volume, t_{s}: the relaxation time of internal spin rotation, H: magnetic field, B: magnetic induction, Ω: internal spin rate, ω: effective rate of rotation of a fluid element ≡ ( 1 / 2 ) ∇ × ν .

The internal angular momentum equation is

I D Ω D t = γ I ∇ 2 Ω − I t s ( Ω − ω ) + μ 0 M × H , (3)

where γ is a dissipation coefficient of inertia spin moment, M: magnetization of the magnetic fluid. The magnetization relaxation equation, Shliomis (1974) [

D M D t = Ω × M − 1 t b ( M − M 0 ) , (4)

where t_{b} is a relaxation time of the particle rotation by Brownian rotation motion, and M_{0} is the equilibrium magnetization of the magnetic fluid, and magnitude is given by

| M 0 | ( | H | , T ) = ϕ M s ( T ) L ( π μ 0 M s ( T ) d 3 | H | 6 k T ) , (5a)

where L ( x ) = coth x − 1 / x is a Langevin function, T: temperature, k: Boltzmann constant, d: diameter of the particle, ϕ : volume fraction of ferromagnetic particles, M_{s}: magnetic moment of a particle, if T ≤ 0.8T_{c} (T_{c}: temperature of Curie point),

M s ( T ) ≈ M s ( 0 ) ( 1 − k T n M B H M ) , (5b)

where M_{B}: Bohrmagneton, n: number of spin perunitatom of ferromagntic metal, H_{M}: molecular field.

The energy equation can be derived by Cowley and Rosensweig [

ρ C m D T D t − μ 0 ρ T Κ ( D | H | D t ) = k 1 ∇ 2 T + Φ , (6)

where C_{m} is the heat capacity at constant volume with a magnetic field, K: pyromagnetic coefficient of magnetic fluid, k_{1}: the thermal conductivity, Φ: a viscous dissipation term usually negligible. The Maxwell’s equation for a non-conducting fluid with no displacement current becomes

∇ × H = 0, ∇ ⋅ B = 0 , (7)

B = μ 0 ( H + M ) . (8)

Cartesian coordinate (x, y) is fixed in the center of an annular pipe, the y-axis being vertically upward, Relationship between cylindrical and Cartesian coordinates is

x + i y = R e i θ , (9)

where, a spectral finite difference method is applied for solving Equation (1) - Equation (8), and primary variables are expressed in a Fourier series to the circumferential θ-direction.

Introduced are a magnetic induction stream functions, ψ i m :

B i x = ∂ ψ i m ∂ y , B i y = − ∂ ψ i m ∂ x , (10)

ψ i m ≡ ψ i m ¯ − μ 0 I 0 2 π ln R , (11)

where i = 1, 2, 3 correspond to the following three regions:

・ R_{0} ≤ R ≤ R_{1}:

∇ 2 ψ 1 m ¯ = 0 , (12a)

・ R_{1} ≤ R ≤ R_{2}:

1 μ 0 ∇ 2 ψ 2 m ¯ + ( M θ R + ∂ M θ ∂ R − 1 R ∂ M R ∂ θ ) = 0 , (12b)

・ R_{2} ≤ R ≤ 50R_{2}:

∇ 2 ψ 3 m ¯ = 0. (12c)

where M_{R} and M_{θ} are R- and θ-component of magnetization of the magnetic fluid, M.

The boundary conditions used are as follows:

・ At R = R_{1}, no-slip, no-spin, and thermal insulation:

[ v ] R = R 1 = 0, [ Ω ] R = R 1 = 0, ( ∂ T ∂ R ) R = R 1 = 0 , (13)

・ At R=R_{2}, no-slip, no-spin, and Dirichlet temperature condition:

[ v ] R = R 2 = 0 , [ Ω ] R = R 2 = 0 , (14)

T ( θ , R 2 ) = T h for θ ∈ ( π 2 , 3 2 π ) , (15a)

T ( θ , R 2 ) = T l for θ ∈ [ 0 , π 2 ] ∪ [ 3 2 π , 2 π ] . (15b)

The normal component of the magnetic induction and the tangential component of the magnetic field are continuous across the interface between dissimilar regions,

at R = R 1 : [ B 1 ] n = [ B 2 ] n , [ H 1 ] t = [ H 2 ] t , (16a)

at R = R 2 : [ B 2 ] n = [ B 3 ] n , [ H 2 ] t = [ H 3 ] t , (16b)

where the subscript n and t denote normal and tangential component. Then the boundary condition of Equations 12(a)-(c) are given by

φ 1 m ¯ = 0 , ∂ φ 1 m ¯ ∂ R = 0 ( R = 0 ) , (17a)

∂ φ 1 m ¯ ∂ R = ∂ φ 2 m ¯ ∂ R + μ 0 M θ ( R = R 1 ) , (17b)

∂ φ 3 m ¯ ∂ R = ∂ φ 2 m ¯ ∂ R + μ 0 M θ ( R = R 2 ) , (17c)

φ 3 m ¯ = 0 ( R = 5 0 R 2 ) . (17d)

Assuming that the magnetic particle of magnetic fluid is iron oxide, and that the carrier liquid is hydrocarbon, the following constants and parameters of magnetic fluid are used

C m = 2.09 × 1 0 3 J ⋅ K − 1 ⋅ kg − 1 , d = 1.0 × 10 − 8 m , H M = 1.7 × 10 9 A ⋅ m − 1 , I = 1.0 × 10 − 4 kg ⋅ m − 1 , k = 1.38 × 10 − 23 J ⋅ K − 1 , k 1 = 0.22 J ⋅ K − 1 ⋅ m − 1 ⋅ s − 1 , K = 220 A ⋅ m − 1 ⋅ K − 1 , M B = 1.17 × 10 − 29 A ⋅ m − 1 , M s ( 0 ) = 3.18 × 10 4 A ⋅ m − 1 , n = 2.2 , ϕ = 0. 041 , η = 0.2 Pa ⋅ s , t b = 7. 6 × 10 − 7 s , t s = 4.17 × 10 − 11 s , μ 0 = 4 π × 1 0 − 7 H ⋅ m − 1 , ρ = 1. 22 × 10 3 kg ⋅ m − 3 .

Since the relaxation time t_{s} (~10^{−11} s) of internal spin rotation and the dissipation coefficient γ ( ≈ t b / 2 ~ 10 − 7 s − 1 ) of inertia spin moment is very small, Equation (3) can be reduced as

I t s ( Ω − ω ) = μ 0 M × H , (18)

and substituting Equation (18) into Equation (2) leads to an approximate momentum equation:

ρ D v D t = − ∇ p + η ∇ 2 v + μ 0 ( M ⋅ ∇ ) H , (19)

and then the curl of the magnetic body force is given by

∇ × μ 0 ( M ⋅ ∇ ) H = μ 0 ∂ χ ∂ T [ ∇ T × ( H ⋅ ∇ ) H ] , (20)

where χ (= χ(T, |H|)) is the magnetic susceptibility of magnetic fluid, and the magnetization of the magnetic fluid can be written as M = χH. As for the Equation (20), the producing condition of thermomagnetic convection is derived in the gravity-free space as ∇ T × ∇ | H | ≠ 0 .

In the case of I_{0} =10 A, T_{h} − T_{l} = 20 K and R_{1}/R_{2} = 0.2, the isotherms, streamlines and dimensionless magnetization field strengths with the position of electric current wire varies are displayed as shown in Figures 2-4, where all quantities are made dimensionless base on a length L_{0} (=R_{2}), a temperature T_{0} (=T_{h} − T_{l}),

and reference velocity V_{0} ( = H 0 μ 0 / ρ ), H_{0} ( = I 0 2 π L 0 ) is representative mag-

netic field. It shows that the convection pattern depends on electric current distribution, depending on which multiple or single circulation is produced. Thus, it is found that thermomagnetic convection is controlled by changing relative direction of temperature gradient to that of a magnetic field.

1) For producing thermomagnetic convection of a magnetic fluid, a spatially non-uniform of temperature and a spatially non-uniform of external magnetic field are required so that ∇ T × ∇ | H | ≠ 0 .

2) In the case that the electric current wire is fixed in center of annular pipe, two symmetrical circulating flows are produced.

3) If the position of electric current wire is different from the center, the center of the circulating flows moves to the side of lower or higher temperature depending on whether the location of electric current wire is to the lower or higher temperature side. Especially in the case of the position of electric current wire deviates greatly from the median line of temperature, only one circulating flow zone is produced.

Liu, J. (2017) Thermomagnetic Convection of Magnetic Fluid in an Annular Space under a Non- Uniform Magnetic and Thermal Field. Ap- plied Mathematics, 8, 655-662. https://doi.org/10.4236/am.2017.85052

B: magnetic induction

C_{m}: heat capacity at constant volume with a magnetic field

d: diameter of the particle

H: magnetic field

H_{0}: external magnetic field

H_{M}: molecular field

I: average inertia moments of particles per unit volume

I_{0:} a direct current

k: Boltzmann constant

k_{1}: thermal conductivity

K: pyromagnetic coefficient of magnetic fluid

M: magnetization of the magnetic fluid

M_{0}: equilibrium magnetization of the magnetic fluid

M_{B}: Bohr magneton

M_{s}: magnetic moment of a particle

n: number of spin per unit atom of ferromagntic metal

p: pressure

R_{0}: radius of conductive wire

R, θ: cylindrical coordinate system

t: time

T: temperature

T_{c}: temperature of Curie point

t_{b}: relaxation time of the particle rotation by Brownian rotation motion

t_{s}: relaxation time of internal spin rotation

v: velocity

x, y: Cartesian coordinate system

γ: dissipation coefficient of inertia spin moment

η: viscosity of magnetic fluid in absence of a magnetic field

μ_{0}: permeability of the magnetic fluid

ρ: mass density of magnetic fluid

ψ : stream function

ψ i m : a magnetic induction stream functions

ϕ : volume fraction of ferromagnetic particles

Φ: viscous dissipation term

χ: magnetic susceptibility of magnetic fluid

ω: effective rate of rotation of a fluid element

Ω: internal spin rate