# On the Designing, Efficiency, and Stirring Force of an Induction Coil for the Processing of Prototype Al Based Nanocomposites

Authors: Muhammad Mansoor; Muhammad Shahid
Journal: Journal of Metallurgy
 (2014)

---

## Abstract

Induction melting could be a potential fabrication method for aluminum matrix nanocomposites, owing to its characteristics stirring action and rapid heating. In the present work, an induction coil for small scale composite fabrication was designed through geometric, thermal, and electromagnetic analyses. Subsequently, the designed coil was simulated for its thermal efficiency and stirring force using finite element method magnetics software. The designed coil yielded more than 60% of the total energy supplied into thermal efficiency with a stirring force <3 mN, which could be significant for efficient melting and stirring to fabricate the aluminum matrix nanocomposites.


---

## Body

## 1. Introduction

All over the word, aluminum based nanocomposites are gaining interest of the researchers exponentially, owing to the facts of their higher mass to strength ratio, ease of fabrication, lower processing temperatures, wide compositional range, and lower costs. However, due to the density differences and nonwetting behavior of molten aluminum, dispersion of strengthening phases or nanoparticles in aluminum matrix is one of the prime issues. Since last decade, many researchers tried to improve the second phase dispersion by powder precursor casting [1], stir casting [2], compo-casting [3], ultrasonic assisted casting [4], and squeeze casting [5]. Although the efforts were successful with respect to the dispersion of strengthening phases certain limitations were always associated with them, for example, contamination from the stirrer in case of stir and compo-castings, limited melt size in case of ultrasonic assisted casting, and reduced melt penetration in larger section for squeeze casting. Moreover, all of the techniques require the specialized type of equipment and capabilities to handle and process molten aluminum, consequently, increasing capital and processing costs.Alternatively, induction melting could be a solution to the problem featuring compositional consistency, liberty of melt quantity, ease of processing, lower costs, and innate stirring. The consequential stirring action of induction melting eliminates any need of external stirrer. In the present work, design parameters of an induction coil were calculated for prototype fabrication of aluminum based nanocomposites. Furthermore, heating efficiency and stirring force induced by the designed coil were evaluated using “finite element method magnetic” version 4.2 (FEMM).


## 2. Background

An induction melting unit comprises a workpiece, an induction coil, load-matching station, and power supply. Besides other parameters, optimum heating could be achieved by designing a compatible and a specific heating coil for a particular application. The heating process involves many factors like electromagnetism, heat transfer, and metallurgical phenomena. Physical properties of the workpiece (usually a metallic material) change with temperature; hence variation in induced electromagnetic and heat transfer will occur [6].The basic mechanism of induction heating is quite simple, where alternating current is produced in a coil (i.e., induction coil) by applying alternating voltage. The alternating current produces magnetic field (of the same frequency) in the vicinity of the coil. The strength of the magnetic field produced depends upon various parameters of the coil like current, geometry, and distance. A workpiece, placed in the coil, experiences eddy currents due to the varying magnetic fields around it. The coil currents and induced eddy currents have the same frequencies but opposite in direction. The induced eddy currents generate separate magnetic fields in opposite direction to the magnetic fields of the coil. Therefore, the induction coil produces a net magnetic field of source and induced magnetism. The induced eddy currents in the workpiece generate heat in a thin surface layer, called skin depth, by ohmic or Joule effect [6, 7]. Besides heating, the induced eddy currents (j) generate Lorentz force (F) in the presence of magnetic flux (B) according to
(1)
F
=
j
×
B
.
In molten form, the conductor behaves as an incompressible fluid and due to the Lorentz force it moves, resulting in continuous stirring of the melt [8, 9]. A schematic of the stirring action of an induction furnace is shown in Figure 1.Figure 1
Schematic diagram of an induction furnace showing stirring action in the molten pool.Electrical resistivity (ρ) and relative magnetic permeability (μ
r) are very important properties of the materials affecting all major parameters of an induction system, for example, skin depth, stirring, heat distribution, coil impedance, and coil efficiency. Induction heating is handy, rapid, localized, and economical heating method, frequently used in industry [10].


## 3. Coil Design Analyses

“No Toil, No Coil” is a common saying of the induction coil designers. Whatever the type of induction heating (with core or coreless) is under consideration, coil is a crucial part, always demanding utmost efforts and time. A major reason for the criticality of coil design is the need of a particular type of coil for a particular application to have optimum results. In the present work, it is required to melt aluminum in a small crucible with maximum stirring action. Therefore, specific design parameters (i.e., geometrical, thermal, and electromagnetic) are analyzed to optimize the coil design.The design is based upon the conditions given in Table1. These conditions are set according to the experimental requirements and limitations posed by available induction generator.Table 1
Preconditions used for the designing of the induction coil.
Material to melt
Aluminum
Capacity of the furnace
0.1 Kg
Crucible material
Sintered alumina
Crucible shape
Cylindrical
Wall thickness of the crucible
2 mm
Gap between coil and the crucible
5 mm
Material of coil pipe
ETP* copper
External diameter of copper pipe
5 mm
Internal diameter of copper pipe
4 mm
Frequency of power supply
10 KHz
Current of power supply
100 A
*Electrolytic tough pitch copper.


### 3.1. Geometrical Analysis

As described in Table1, aluminum (100 g) is required to melt in a cylindrical alumina crucible. The dimensions of the crucible [11] and volume of the melt are calculated using (2) and (3), respectively, as follows:
(2)
H
D
=
1.6
-
2.0
(3)
V
=
π
d
2
H
4
,
where H and D are height and diameter of molten metal in the crucible, respectively.The internal diameter (D
°) of the induction coil is a function of crucible diameter (D
c), crucible wall thickness (B
w), and insulation thickness (B
ins) and is determined using [11, 12]
(4)
D
°
=
D
c
+
B
w
+
B
ins
.
Height of the inductor coil (H
°) depends upon the height of the workpiece or the melt (H), which is usually kept 1.1 to 1.2 times of the H [12]:
(5)
H
°
=
1.1
~
1.2
H
.


### 3.2. Thermal Analysis

The heat energy required (Q
r) to process the aluminum includes energy to melt aluminum (Q
a), energy to superheat (Q
s), energy to melt slag (Q
g), energy for endothermic reaction (Q
e), and energy liberated by any exothermic reaction (Q
x) [13]:
(6)
Q
r
=
Q
a
+
Q
s
+
Q
g
+
Q
e
+
Q
x
,
(7)
Q
a
=
m
C
T
°
-
T
+
E
,
where m is mass of charge in Kg, C is specific heat capacity of aluminum and is equal to 1100 J/Kg·K, T
° is melting temperature (933 K for aluminum), and T is room temperature (298 K). E is the heat energy required for any phase transformations, which would be zero for pure aluminum.Similarly,(8)
Q
s
=
m
C
°
T
s
,
where C
° is heat capacity of aluminum in molten form (992 J/kg·K) and T
s is the superheating temperature (323 K).Also,(9)
Q
g
=
m
g
C
g
,
where m
g is the amount of slag, which is considered 8% of the total melt, and C
g is slag heat energy (18 KJ/kg).


### 3.3. Electromagnetic Analysis

The heat energy produced (Q
p) to melt the charge in the form of eddy currents is a function of frequency of the power supply (f), height of the workpiece (H), maximum magnetic flux density (B), diameter of the workpiece (D), and resistivity (ρ) of the work piece (which is 2.83
×
10
-
8
Ωm for aluminum). The function is shown in [11, 14]
(10)
Q
p
=
π
3
f
2
H
B
2
D
4
8
ρ
or
(11)
B
=
8
ρ
Q
p
π
3
f
2
D
4
H
and also
(12)
Q
p
=
Q
r
t
,
where t is time in seconds to attain maximum flux.The current density in the inductor is(13)
I
a
=
I
A
.
To produce required amount of magnetic flux density (B), number of turns of the inductor coil could be estimated by [11–15]
(14)
B
=
μ
μ
°
N
I
L
or
(15)
N
=
B
L
μ
°
I
,
where μ and μ
° are relative permeability and permeability of workpiece material and free space, respectively. Relative permeability of nonmagnetic materials is equal to 1, and permeability of free space is 4
π
×
10
-
7 Tm−1. L is length of the inductor, which is H in present work. The current passing through the inductor is I, in amperes.


### 3.4. Efficiency of the Designed Induction Heating System

After defining the design parameters it is necessary to evaluate the efficiency of the designed coil prior to its practical implementation. Walter and Ceglia [16] proposed a mathematical model (16) for the evaluation of heating efficiency of induction coil:
(16)
η
°
=
R
s
R
s
+
R
p
×
100
,
where η
° is percentage heat efficiency of the coil, R
p is the primary coil resistance, and R
s is load (or secondary coil) resistance. The losses within the induction coil are associated with the resistance of the coil (R
p) and, therefore, kept to minimum values. However, R
s defines the extent of induced energy in the workpiece and should be maximized.


### 3.5. Simulation

To find the values ofR
p and R
s, simulation based upon calculated parameter of the induction coil was carried out using FEMM. Geometry of the heating system was defined in FEMM according to Figure 2, where external diameter of the workpiece (aluminum) was taken as the diameter of the crucible neglecting the actual material of the crucible (alumina), which was inactive for magnetic inductions, whereas copper was taken as coil material.(a) Actual crucible used for melting, (b) transversal cutting of the crucible used for simulation, and (c) FEMM drawing.
(a)
(b)
(c)The simulation was required for axisymmetric condition; therefore only a traversal cutting of the crucible was drawn (Figure2(c)). According to the program’s algorithm, FEMM first determined the magnetic contributions of 2D slice and then evaluated remaining by full revolution imparting the solid.All the regions were defined in the preprocessor of the FEMM with their respective materials to realize the working problem; that is, workpiece was defined as “Aluminum 1100” and conductors of the coil were defined as “Copper” having electrical conductivity values of 34.45 and 58 MS/m, respectively. The copper coil was defined in one circuit named “coil” and environment around the workpiece and inductor was defined as “air.” The inside area of the copper tube was not used in the simulation hence defined as “No mesh” (Figure3).(a) Definition of all the working areas and respective materials and (b) preprocessor parameters used to define the problem.
(a)
(b)


## 3.1. Geometrical Analysis

As described in Table1, aluminum (100 g) is required to melt in a cylindrical alumina crucible. The dimensions of the crucible [11] and volume of the melt are calculated using (2) and (3), respectively, as follows:
(2)
H
D
=
1.6
-
2.0
(3)
V
=
π
d
2
H
4
,
where H and D are height and diameter of molten metal in the crucible, respectively.The internal diameter (D
°) of the induction coil is a function of crucible diameter (D
c), crucible wall thickness (B
w), and insulation thickness (B
ins) and is determined using [11, 12]
(4)
D
°
=
D
c
+
B
w
+
B
ins
.
Height of the inductor coil (H
°) depends upon the height of the workpiece or the melt (H), which is usually kept 1.1 to 1.2 times of the H [12]:
(5)
H
°
=
1.1
~
1.2
H
.


## 3.2. Thermal Analysis

The heat energy required (Q
r) to process the aluminum includes energy to melt aluminum (Q
a), energy to superheat (Q
s), energy to melt slag (Q
g), energy for endothermic reaction (Q
e), and energy liberated by any exothermic reaction (Q
x) [13]:
(6)
Q
r
=
Q
a
+
Q
s
+
Q
g
+
Q
e
+
Q
x
,
(7)
Q
a
=
m
C
T
°
-
T
+
E
,
where m is mass of charge in Kg, C is specific heat capacity of aluminum and is equal to 1100 J/Kg·K, T
° is melting temperature (933 K for aluminum), and T is room temperature (298 K). E is the heat energy required for any phase transformations, which would be zero for pure aluminum.Similarly,(8)
Q
s
=
m
C
°
T
s
,
where C
° is heat capacity of aluminum in molten form (992 J/kg·K) and T
s is the superheating temperature (323 K).Also,(9)
Q
g
=
m
g
C
g
,
where m
g is the amount of slag, which is considered 8% of the total melt, and C
g is slag heat energy (18 KJ/kg).


## 3.3. Electromagnetic Analysis

The heat energy produced (Q
p) to melt the charge in the form of eddy currents is a function of frequency of the power supply (f), height of the workpiece (H), maximum magnetic flux density (B), diameter of the workpiece (D), and resistivity (ρ) of the work piece (which is 2.83
×
10
-
8
Ωm for aluminum). The function is shown in [11, 14]
(10)
Q
p
=
π
3
f
2
H
B
2
D
4
8
ρ
or
(11)
B
=
8
ρ
Q
p
π
3
f
2
D
4
H
and also
(12)
Q
p
=
Q
r
t
,
where t is time in seconds to attain maximum flux.The current density in the inductor is(13)
I
a
=
I
A
.
To produce required amount of magnetic flux density (B), number of turns of the inductor coil could be estimated by [11–15]
(14)
B
=
μ
μ
°
N
I
L
or
(15)
N
=
B
L
μ
°
I
,
where μ and μ
° are relative permeability and permeability of workpiece material and free space, respectively. Relative permeability of nonmagnetic materials is equal to 1, and permeability of free space is 4
π
×
10
-
7 Tm−1. L is length of the inductor, which is H in present work. The current passing through the inductor is I, in amperes.


## 3.4. Efficiency of the Designed Induction Heating System

After defining the design parameters it is necessary to evaluate the efficiency of the designed coil prior to its practical implementation. Walter and Ceglia [16] proposed a mathematical model (16) for the evaluation of heating efficiency of induction coil:
(16)
η
°
=
R
s
R
s
+
R
p
×
100
,
where η
° is percentage heat efficiency of the coil, R
p is the primary coil resistance, and R
s is load (or secondary coil) resistance. The losses within the induction coil are associated with the resistance of the coil (R
p) and, therefore, kept to minimum values. However, R
s defines the extent of induced energy in the workpiece and should be maximized.


## 3.5. Simulation

To find the values ofR
p and R
s, simulation based upon calculated parameter of the induction coil was carried out using FEMM. Geometry of the heating system was defined in FEMM according to Figure 2, where external diameter of the workpiece (aluminum) was taken as the diameter of the crucible neglecting the actual material of the crucible (alumina), which was inactive for magnetic inductions, whereas copper was taken as coil material.(a) Actual crucible used for melting, (b) transversal cutting of the crucible used for simulation, and (c) FEMM drawing.
(a)
(b)
(c)The simulation was required for axisymmetric condition; therefore only a traversal cutting of the crucible was drawn (Figure2(c)). According to the program’s algorithm, FEMM first determined the magnetic contributions of 2D slice and then evaluated remaining by full revolution imparting the solid.All the regions were defined in the preprocessor of the FEMM with their respective materials to realize the working problem; that is, workpiece was defined as “Aluminum 1100” and conductors of the coil were defined as “Copper” having electrical conductivity values of 34.45 and 58 MS/m, respectively. The copper coil was defined in one circuit named “coil” and environment around the workpiece and inductor was defined as “air.” The inside area of the copper tube was not used in the simulation hence defined as “No mesh” (Figure3).(a) Definition of all the working areas and respective materials and (b) preprocessor parameters used to define the problem.
(a)
(b)


## 4. Results and Discussion

The calculated geometrical parameters of the required coil for the melting of 100 g aluminum are given in Table2. A total clearance of 7 mm was kept to accommodate alumina crucible having 2 mm wall thickness and working clearance of 5 mm. No additional insulation layer was used between the crucible and the coil. Height of the coil was determined using (5) with upper allowable limit to induce maximum magnetic flux in the workpiece; however, not more than 20% of the height of the molten surface remained.Table 2
Results of the design analyses.
Geometrical analysis
Height of molten metal (H)
5.5 cm
Diameter of molten metal (D)
3.0 cm
Volume of the molten melt (V)
38.9 cm³
Internal diameter of the induction coil (D
°)
4.5 cm
Height of the inductor coil (H
°)
6.0 cm
Thermal analysis
Energy to melt aluminum (Q
a)
69850 J
Energy to superheat molten aluminum (Q
s)
32042 J
Energy to melt slag (Q
g)
144 J
Total heat energy required (Q
r)
102036 J
Electromagnetic analysis
Maximum magnetic flux density (B)
0.00409 T
Current density through the coil (I
a)
5.67 A⋅mm−2
Number of turns of the coil (N)
4.47 turnsTotal thermal energy required by the system is ~105 KJ, as calculated by (6). The system under consideration entails melting of pure metal and incorporation of strengthening particles, which will remain unmelted and/or inert within the melt; therefore, for the calculation, endothermic (Q
e) and exothermic (Q
x) energy changes were considered equal (i.e., Q
e
=
Q
x theoretically).A very critical aspect of an induction coil is “number of turns,” which was determined after calculating induced magnetic flux and current densities. The calculated turns were rounded-off to five and a pitch of 8 mm was maintained during fabrication the coil. Results of the theoretical analyses are given in Table2.As operational and geometrical parameters were defined (Tables1 and 2), the simulation of the heating system was generated using FEMM for the distribution of magnetic field intensities. Magnetic flux distribution in the workpiece and around the inductor is shown in Figure 4(a). The electrical parameters of the heating circuit were evaluated using a postprocessor facility provided by FEMM. The results are shown in Figure 4(b). The parameter “Voltage/Current” directly yields the value of (R
p
+
R
s) of  (16). To obtain R
p value, an additional simulation was carried out where no load or workpiece was used (Figures 4(c) and 4(d)), which was done by eliminating workpiece while keeping all other parameters [16]. Therefore,
(17)
R
p
+
R
s
=
0.01274
Ω
or
≅
12.7
m
Ω
,
R
p
=
0.004768
Ω
or
≅
4.77
m
Ω
.
Thus, according to (16),
(18)
η
°
=
R
s
R
p
+
R
s
×
100
=
62.6
%
.
Similarly, Lorentz force exerted by the induced magnetic flux in the workpiece was determined using postprocessor option of FEMM. Lorentz force is the cross product of magnetic flux density and induced current, that is, (1). To determine the Lorentz force induced in the workpiece, contour points of the workpiece were selected to define the area. Subsequently, block integral (a provision of FEMM postprocessor) was carried out by selecting “Lorentz force” option, which calculated a force of 2.45 mN generated in the workpiece under the given preprocessor conditions. When the workpiece (aluminum) will become molten, the generated force may stir the melt smoothly and avoid settling of the reinforcing nanoparticles.(a) Magnetic flux distribution in the designed coil system. (b) postprocessor results of the circuit “coil.” (c) and (d) same as “(a) and (b)” but workpiece is eliminated.
(a)
(b)
(c)
(d)As pure aluminum is excellent conductor of electricity and heat, small amount of eddy currents will be generated during induction heating. However, the generated heat will transfer rapidly to the areas away from the proximity of the induction coil. Therefore, it would be beneficial to use relatively larger pieces of aluminum than smaller ones. Additionally, large pieces having lower surface areas will be less exposed to oxidation, resulting in ease of melting along with low entrapped oxides.


## 5. Conclusions

An induction coil was designed and analyzed for its geometrical, thermal, and electromagnetic parameters. The heating efficiency as well as induced stirring force of the coil was determined using FEMM. At the induction generator parameters (i.e., frequency of 10 KHz and current of 100 A) a heating efficiency >60% and stirring force <3 mN were achievable.



---

*Source: 637031-2014-12-31.xml*