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)  , Polevikov & Fertman (1977)  , Busse & Riahi (1982)  , Stiles & Kagan (1990)  , Zebib (1996)  , Fruh (2005)  , Yamaguchi (2009)  and Siddiqa (2013)  .
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.
2.1. Physical Model
Analyzed is Two-dimensional thermomagnetic convection followed by a steady- state natural convection of magnetic fluid as shown in Figure 1, where a conductive wire of radius R0 with a direct current I0 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 Tl, and at t > 0 a half part of the outside circumference is assumed to be fixed at a high temperature (Th), and the rest half part of the outside circumference at a low temperature(Tl).
The following assumption applies:
1) The radius of wire is much smaller than that of the inside circle ( ).
2) The effect of magnetic field is limited within the distance 50R2.
3) The magnetic permeability of magnetic fluid is equal to that of vacuum.
4) The magnetic fluid is incompressible fluid.
Figure 1. The analytical region.
2.2. Governing Equations of an Incompressible Magnetic Fluid
The equation of continuity is
The momentum equation neglecting gravity force is
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, ts: the relaxation time of internal spin rotation, H: magnetic field, B: magnetic induction, Ω: internal spin rate, ω: effective rate of rotation of a fluid element .
The internal angular momentum equation is
where γ is a dissipation coefficient of inertia spin moment, M: magnetization of the magnetic fluid. The magnetization relaxation equation, Shliomis (1974)  , is
where tb is a relaxation time of the particle rotation by Brownian rotation motion, and M0 is the equilibrium magnetization of the magnetic fluid, and magnitude is given by
where is a Langevin function, T: temperature, k: Boltzmann constant, d: diameter of the particle, : volume fraction of ferromagnetic particles, Ms: magnetic moment of a particle, if T ≤ 0.8Tc (Tc: temperature of Curie point),
where MB: Bohrmagneton, n: number of spin perunitatom of ferromagntic metal, HM: molecular field.
The energy equation can be derived by Cowley and Rosensweig  as
where Cm is the heat capacity at constant volume with a magnetic field, K: pyromagnetic coefficient of magnetic fluid, k1: the thermal conductivity, Φ: a viscous dissipation term usually negligible. The Maxwell’s equation for a non-conducting fluid with no displacement current becomes
2.3. Coordinate System
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
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.
2.4. Equation of the Distribution of Magnetic Field
Introduced are a magnetic induction stream functions, :
where i = 1, 2, 3 correspond to the following three regions:
・ R0 ≤ R ≤ R1:
・ R1 ≤ R ≤ R2:
・ R2 ≤ R ≤ 50R2:
where MR and Mθ are R- and θ-component of magnetization of the magnetic fluid, M.
2.5. Boundary Conditions
The boundary conditions used are as follows:
・ At R = R1, no-slip, no-spin, and thermal insulation:
・ At R=R2, no-slip, no-spin, and Dirichlet temperature condition:
The normal component of the magnetic induction and the tangential component of the magnetic field are continuous across the interface between dissimilar regions,
where the subscript n and t denote normal and tangential component. Then the boundary condition of Equations 12(a)-(c) are given by
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
, , , , , , , , , , , , , , , .
3.1. Magnetic Body Force of Magnetic Fluid
Since the relaxation time ts (~10−11 s) of internal spin rotation and the dissipation coefficient γ ( ) of inertia spin moment is very small, Equation (3) can be reduced as
and substituting Equation (18) into Equation (2) leads to an approximate momentum equation:
and then the curl of the magnetic body force is given by
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 .
3.2. Numerical Results
In the case of I0 =10 A, Th − Tl = 20 K and R1/R2 = 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 L0 (=R2), a temperature T0 (=Th − Tl),
and reference velocity V0 ( ), H0 ( ) 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 .
Figure 2. Isotherms (at an interval of 0.1), streamlines(at an interval of 1.0 × 10−7 stream function) and dimensionless magnetization field strengths (at an interval of 1.0 × 10−3 with from 0.001 to 0.007) in the case of the electric current wire fixed at (x, y) = (0, 0), where the direction of velocity along the center is leftward.
Figure 3. Isotherms (at an interval of 0.1), streamlines (at an interval of 2.0 × 10−8 stream function) and dimensionless magnetization field strengths (at an interval of 2.0 × 10−3 with from 0.002 to 0.02) in the case of the electric current wire fixed at (x, y) = (0.15, 0), where the direction of velocity along the center is leftward.
Figure 4. Isotherms (at an interval of 0.1), streamlines (at an interval of 1.0 × 10−7 stream function) and dimensionless magnetization field strengths (at an interval of 2.0 × 10−3 with from 0.002 to 0.02) in the case of the electric current wire fixed at (x, y) = (0, 0.15), where the direction of flow is clockwise.
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.
B: magnetic induction
Cm: heat capacity at constant volume with a magnetic field
d: diameter of the particle
H: magnetic field
H0: external magnetic field
HM: molecular field
I: average inertia moments of particles per unit volume
I0: a direct current
k: Boltzmann constant
k1: thermal conductivity
K: pyromagnetic coefficient of magnetic fluid
M: magnetization of the magnetic fluid
M0: equilibrium magnetization of the magnetic fluid
MB: Bohr magneton
Ms: magnetic moment of a particle
n: number of spin per unit atom of ferromagntic metal
R0: radius of conductive wire
R, θ: cylindrical coordinate system
Tc: temperature of Curie point
tb: relaxation time of the particle rotation by Brownian rotation motion
ts: relaxation time of internal spin rotation
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
: 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