Abstract
The abundance of low-grade waste heat necessitates energy harvesting devices to convert thermal energy to electrical energy. Through magnetic transduction, thermomagnetics can perform this conversion at reasonable efficiencies. Thermomagnetic materials use thermal energy to switch between magnetic and non-magnetic states and convert thermal energy into electrical energy. In this study, we numerically analyzed an active thermomagnetic device for thermal energy harvesting composed of gadolinium (Gd) and neodymium iron boron (NdFeB). A parametric study to determine the device efficiency was conducted by varying the gap distance, heat source temperature, and Gd thickness. Furthermore, the effect of the thermal conductance and applied field was also evaluated. It was found that the relative efficiency for smaller gap distances ranges from ∼15% to 28%; the largest allowable volume of Gd should be used and higher applied field leads to higher efficiencies.
Introduction
With an ever-growing need for alternative energy sources, the abundance of thermal energy has spawned an array of conversion technologies for various applications. In 2017, the USA increased its clean energy focus recognizing the waste heat generated by methods such as chemical processing, vehicle operation, or other devices amounted to 66.7 quadrillion BTU.2 According to Advanced Research Projects Agency-Energy (ARPA-E), for a temperature range of 50–700 °C, nearly 75% of the waste heat is low-grade (≤230 °C). Additionally, assuming a cold side temperature of 25 °C, 45% of the work potential is accessible using low-grade waste heat converted at a Carnot efficiency of ∼40% [1]. Therefore, waste heat represents a pervasive form of energy that can be converted into useful electrical energy.
Various thermal conversion technologies have been investigated but thermoelectric (Seebeck) devices remain the primary focus due to their long history of directly transducing temperature gradients to an electric energy [2–4]. Though numerous research efforts investigating Seebeck devices have been conducted, their performance for low-temperature (<100 °C) energy conversion remains limited due to relatively low conversion efficiency of <20% of Carnot and a power conversion efficiency of less than 5% [3]. Another approach that has received modest attention is thermomagnetic devices which convert thermal energy to electrical energy through magnetic transduction [4–22].
Thermomagnetic transduction consists of a system containing a soft magnet and a hard magnet where the soft magnet undergoes a thermally induced second-order ferromagnetic to paramagnetic phase transition. The transition temperature between the phases is known as the Curie point or Curie temperature. Energy is harvested using a passive or active thermomagnetic device [23]. For passive devices, the soft magnet is attached to a piezoelectric material that converts the mechanical oscillations into electricity through a two-step process. The temperature change causes the soft magnet to oscillate between a hot and a cold surface where the hard magnet is located at the hot side with mechanical motion producing electrical energy, i.e., piezoelectric effect. Active devices use a soft magnet as a shunt in a permanent magnet. A coil surrounds the configuration and the thermal transition switches the magnetic field, converting thermal energy directly to electric energy [23] and thus avoids the two-step process in passive devices.
Devices were initially developed and patented in the late 19th century [5,16]. However, it was not until 1948 when Brillouin and Iskenderian determined that thermally oscillating ferromagnetic materials above and below the Curie Point could produce large relative efficiencies (i.e., 55% of Carnot) [8]. The analytical predictions were substantially larger than Seebeck devices, but they required large magnetic fields and did not account for the disorder present in the system when transitioning across the second-order phase transition limiting operating efficiency. In 1959, Elliott et al. experimentally evaluated gadolinium (Gd) for use in a thermomagnetic energy harvesting system operating at room temperature [24]. However, it failed to produce significant power outputs due to the lack of sufficiently strong magnets to bias the structures with high fields.
The advent of high remanence rare earth magnet catalyzed further development. In the 1980s and early 1990s, several device iterations were developed to increase efficiency bias field variance or through regeneration [6,9,13]. More recently, Ujihara et al. developed a strain mediated multiferroic device for thermal to electrical transduction using Gd [12] relying on a two-step process that requires piezoelectric springs for transduction from thermal to electrical energy. This particular device was evaluated for power generation rather than efficiency. Later, Bulgrin et al. analytically optimized Ujihara’s device and determined that the usage of ferromagnetic materials which transition over smaller temperature intervals increases the efficiencies [10]. Furthermore, Hsu experimentally evaluated poly-crystalline Gd and analytically postulated that using nanoscale single domain Gd elements could increase relative efficiencies to 30% of Carnot [25]. Chun et al. further developed a thermo-magneto-electric generator and evaluated the device for waste heat recovery from a central processing unit and reported promise in electricity generation as well as heat dissipation (cooling) [22].
In this study, an active thermomagnetic device, which directly transduces thermal to electrical energy, was developed using a finite element analysis (FEA) model. The model was numerically evaluated based on several device parameters that affect operational efficiency for low-grade waste heat recovery [10,21]. The focus of this study centers on the effects of operational conditions assuming maximum transduction efficiency. The transduction efficiency is the conversion of thermal energy to magnetic energy and is strongly sensitive to material selection and processing [15]. The maximum transduction efficiency is used to decouple the materials’ conversion ability due to processing from its performance through varying external operational conditions. Thus, the lens is placed on the conversion of thermal to electrical energy using a highly efficient material.
Parametric studies varying the distance between the hot and cold sources, the hot source temperature and the thickness of Gd for applied magnetic field, and normal and enhanced contact conductance were used to assess their impact on efficiencies. The conclusion of this study found that the development of a device with the smallest gap distance, while utilizing the largest volume of ferromagnetic material that operates with a Curie temperature close to the hot source temperature, is optimal for increasing the thermal to magnetic efficiencies.
Thermomagnetic Device Model Development
This active energy harvesting system is based on the motion of the ferromagnetic material within a solenoid coil as illustrated in Fig. 1. The mechanical motion of the ferromagnetic material coupled with its magnetic changes during this motion produces electrical current in the coil which can be harvested. The system modeled in this study consists of a ferromagnetic material, Gd, a neodymium iron boron (NdFeB) permanent magnet, a solenoid coil, mechanical springs, and a hot/cold source.
As seen in Fig. 1(a), the ferromagnetic Gd is attached to mechanical springs connected to the cold source at the bottom of the figure. The duty cycle of the device is explained with the four stage process illustrated in Fig. 1(b) representing a cross section of Fig. 1(a). In stage I, the Gd is at room temperature and in contact with the cold surface to initiate the thermal energy harvesting cycle. At room temperature, the Gd is below the Curie temperature Tc and is thus ferromagnetic. At this temperature, the dipolar fields attract the Gd to the NdFeB permanent magnet. In stage II, the magnetic force between the Gd and NdFeB exceeds the spring forces and any gravitational forces to translate the Gd toward the NdFeB permanent magnet. In accordance with Faraday’s law, this movement generates a voltage within the surrounding coil. At stage III, the Gd is in contact with the hot side and is heated above the Tc to become paramagnetic. At stage IV, the spring and gravitational forces exerted on the Gd exceed the magnetic force (i.e., Gd is paramagnetic at this stage) and the Gd moves toward the cold source.
The above energy harvesting system was analyzed using comsol 5.1 with a FEA model using a 2D axisymmetric approach as shown in Fig. 2(a). A 360-deg revolution of the 2D axisymmetric model results in the 3D formation (Figs. 2(b) and 2(c)). Triangular mesh elements were used for every domain except the air gap above and below the Gd, as seen in Fig. 3. The air gap contains a mapped mesh with rectangular elements that allows for expansion and compression during the movement of the Gd domain. A smaller element size is assigned to the vertical boundary between the coil and the NdFeB/Gd/C domains. The smaller element size increases numerical accuracy at the interface due to large gradient terms. A time-dependent, fully coupled study was used with the pardiso solver.

Model of the thermomagnetic device. (a) 2D axisymmetric view, (b) 3D view, and (c) 3D top view where the colored circle is in the top cross-sectional area of Gd.
![Magnetization of a gadolinium as a function of field at various temperatures [26]](https://asmedc.silverchair-cdn.com/asmedc/content_public/journal/energyresources/142/8/10.1115_1.4046273/1/m_jert_142_8_082102_f004.png?Expires=1655032278&Signature=2FQHcTIzgMRtEGzMmJENYV-mGTZ9fAe5mwqCGlM3zhKP4MFmytnLNmHpG-Z-LdKx4rALo9G6Cvc0qI9kKvXdf8B7DmltiKNYFARnSu4APELF328~SLDinfQib3anU~wKftuEfR5zP~FZ-byyn573CUpoQVYpHm7m4QwL9zsr97-O7ZH3zLZCiQfSF-2-3SqPmIggXRNHMzq36O4yBuGrf-igv9hPMu-Ta~OIQxxbMRMA3YRL8JxI7NvaCjm2tCUO13TENSeSxYBlYYXhbG3r8c~b89aXMavNOwFCBz6uYr97vOLru~2Xu1xX5Gkx8C1z98nXp12nMde6i8goujlpzg__&Key-Pair-Id=APKAIE5G5CRDK6RD3PGA)
Magnetization of a gadolinium as a function of field at various temperatures [26]
An overview of the parametric study performed is shown in Fig. 5. The base model was designed using typical geometry (gap distance, Gd thickness) and operational (contact conductance, applied field, heat source temperature) inputs. The first parametric study was compared using a 1.32 T versus 2 T NdFeB generated applied field. The second parametric study compared the normal versus enhanced contact conductance. The influence of gap distance (zg), thickness of the Gd (thGd), and hot source temperature (TH) was evaluated for each parametric study as well as the influence on efficiency.
Model inputs
Base | Range | |
---|---|---|
Gap distance (mm) | 2 | 1–2 |
Gd thickness (mm) | 0.1 | 0.1–1.0 |
Gd initial temperature (K) | 286 | – |
Heat source (K) | 300 | 300–320 |
Heat sink (K) | 285 | – |
Spring constant (N/m) | 1 | – |
# of coil turns | 800 | – |
Applied field (T) | 1 (normal) | 2 (enhanced) |
h (W/m2 K) | 650 (normal) | 2 × 103 (enhanced) |
Ω[ohms] | 10.53 | 10.53 |
Base | Range | |
---|---|---|
Gap distance (mm) | 2 | 1–2 |
Gd thickness (mm) | 0.1 | 0.1–1.0 |
Gd initial temperature (K) | 286 | – |
Heat source (K) | 300 | 300–320 |
Heat sink (K) | 285 | – |
Spring constant (N/m) | 1 | – |
# of coil turns | 800 | – |
Applied field (T) | 1 (normal) | 2 (enhanced) |
h (W/m2 K) | 650 (normal) | 2 × 103 (enhanced) |
Ω[ohms] | 10.53 | 10.53 |
Results and Discussion
Prior to discussing the study results, it is helpful to understand representative outputs from the model. Figure 7 plots Gd position and temperature versus time, magnetic force versus time, and voltage versus time from a model simulated using the base properties in Table 1. These data were collected and used to determine the frequency, power, and relative efficiency per cycle. In Fig. 7(a), when the position is 0 mm, the lower boundary of the Gd was in contact with the heat sink. Then, the Gd moved from the heat sink to the heat source, stopping at the heat source. In Fig. 7(a), the Gd stopped at 1.75 mm, where the top boundary of the Gd was in contact with the heat source. Once the Gd reached the heat source, its temperature increased (i.e., T = 299 K) near the Curie temperature and the magnetization approached “0” allowing the Gd to return to the heat sink. The effects of the changes in magnetization can be seen closer in Fig. 7(b). The magnetic force is dependent upon the magnetization of the Gd which is dependent on the temperature of the Gd. Initially, the temperature of the Gd is below the Curie temperature which enables full magnetization. This magnetization contributes to the magnetic force (as described in Eq. (5)) which attracted mechanically the Gd toward the NdFeB heat source. The heating of the Gd increased the temperature to the Cure temperature, where the magnetization reached “0” and subsequently the magnetic force decreased to 0 mN. The spring force restored the Gd to the heat sink. As can be seen in the figures, the system exhibits a delay in the initialization of the return to the heat sink due to remanence magnetization.

(a) Gd position and temperature versus time, (b) magnetic force felt by Gd versus time, and (c) voltage generated versus time
The case results for frequency, power, and relative efficiency are displayed in the remainder of this section. For the applied field parametric study, the case results will be referred to as 1.32 T or 2 T. For the contact conductance parametric study, the case results will be referred to as NORMAL or ENHANCED. As a reminder, the parameters used in the studies are presented in Table 1.
Frequency
Figure 8 displays the frequency versus the three parameters (1) gap distance (zg), (2) hot source temperature (TH), and (3) Gd thickness (thGd) for both applied fields and both contact conductance cases. Though frequency was calculated from voltage induction, it was dependent on the thermal cycle. An inverse correlation between frequency and each parameter (zg, TH, and thGd) was expected. The increase in zg increased the distance between the heat source and the sink, and thus a longer thermal cycle, which leads to a lower frequency. In Fig. 8(a), this occurred for the 1.32 T case, and the 2 T case showed a resonant frequency at the gap distance of 1.75 mm. Generally, when operating at the resonant frequency, power generation for a given time duration increases due to an increase in the number of cycles [7]. However, the magnitude of the voltage generated is still magnetic field gradient dependent. Additionally, the emergence of the resonant frequency for the 2 T and absence for the 1.32 T case gives insight to a correlation between the gap distance, Gd volume, and applied field. This information could possibly foreshadow the scalability of the device. Further investigation is needed to determine the extent of its relevance and relationship to efficiency. As the TH increased, the heat flux into the Gd was increased and the magnetic force decreases similarly, increasing the cycle duration. The larger cycle duration is directly responsible for reducing the frequency. Figure 8(b) displays the frequency for 1.32 T and 2 T as a result of that performance. A reduction in thGd decreased the time it took the magnetic material shuttled through the coil but it also reduced the thermal cycle. The reduced thermal cycles result in a net increase in frequency. This is observed in Fig. 8(c) for both 1.32 T and 2 T. They both had the same frequencies except at the smallest thGd. This could be due to the existence of a minimal amount of magnetic material within the 0.25–0.1 mm range for each applied field.

Frequency versus gap distance, heat source temperature, and Gd thickness for (a)–(c) the applied field study and (d)–(f) the contact conductance study
Furthermore, Ujihara et al. and Bulgrin et al. have demonstrated that increasing the contact conductance will increase the frequency. In Fig. 8(d), as the Zg increased, the frequency for ENHANCED is higher (as expected) but decreased (∼18.3%) compared to (∼10.5%) for NORMAL. ENHANCED decreased almost two times faster because the higher contact conductance affects the thermomagnetic cycle by increasing the rate of thermal activation. Thus, the thermal activation driven by the heat source and heat sink is more sensitive to the distance between them. The frequencies for NORMAL and ENHANCED (Fig. 8(e)) were the highest at TH = 300 K or 305 K, respectively, and decreased with minimal inflections as the temperature increased. This behavior is evident of the need for TH to be selected as close to the Curie temperature as possible to ensure a higher frequency. Additionally, a high contact conductance will allow operation at a higher temperature but within certain limits. Determining the optimal temperature operation range as a function of contact conductance could represent the next steps in device design. ENHANCED had a higher frequency than NORMAL at each thickness but with the same correlation. ENHANCED had the highest frequency for each parameter further exemplifying the influence of contact conductance on generating higher frequency. However, thGd had the largest range for each comparison case (∼2.5 to 6 Hz for ENHANCED and ∼1 to 4 Hz for NORMAL). The magnetic force is dictated by the thGd. Though a reduction in thickness means a reduction in magnetic material and thus magnetic force, influencing the travel velocity, it also means faster conduction. This can be seen clearly in Fig. 8(f) as ENHANCED was a combination of higher conduction and contact conductance. Overall, the contact conductance is more influential than the applied field in targeting higher frequencies. Moreover, smaller zg values lead to higher frequencies.
Power
Figure 9 provides plots of power versus the parameters (gap distance (zg), hot source temperature (TH), and Gd thickness (thGd)) for the applied field and contact conductance cases. Power and gap distance are inversely correlated to an increase in zg. The applied field case demonstrated that an increase in magnetic field generated more power while the contact conductance case was less definitive which differs from the frequency response. All the 2 T cases generated higher power than the applied 1.32 T case. In Fig. 9(a), the largest difference in power occurred at the smallest zg and the smallest difference occurred at the largest zg. The applied magnetic field had more influence at smaller distances. Figure 9(b) shows maxima TH at 305 K and 310 K for 1.32 T and 2 T, respectively. This higher power is attributed to either an extension of the ferromagnetic region, a decrease in frequency or both. However, it can be observed that 2 T case increased the power output by ∼128%. The larger thGd had the highest power generation due to the larger amount of magnetic material that is critical to voltage generation.

Power versus gap distance, heat source temperature, and Gd thickness for (a)–(c) the applied field study and (d)–(f) the contact conductance study
In Figs. 9(d) and 9(f), the power output is nearly the same for NORMAL and ENHANCED. The power steadily declines until the 1.5 mm distance (Fig. 9(d)). When zg is larger than 1.5 mm, the correlation weakened. For the TH, NORMAL had a maxima at 305 K and the power generation for ENHANCED increased after 310 K. The ENHANCED case exhibits a slight frequency increase after 310 K, which provides a bump in the power generation. The power output increased by ∼14% compared to the NORMAL case. This further leads to the concept that increased frequency can increase power generation and it also demonstrates the possibility of expanding the heat source temperature at which the device operates. Similar to the applied field case, the largest thGd had the highest power generation. Overall, the conclusion remains that higher applied fields generate higher power.
Relative Efficiency
Figure 10 plots the relative efficiency versus the parameters (gap distance (zg), hot source temperature (TH), and Gd thickness (thGd)) for the two applied field and two contact conductance cases. The relative efficiency for both cases is inversely related to zg, TH, and thGd. The efficiency values for zg follow a similar same trend as the gap distance cases observed in power generation (Figure 9). In Fig. 10, the 2 T case had the highest efficiency for all parameters studied compared to the lower 1.32 T case. Figures 10(a)–10(c) show relative efficiencies of ∼28% for the smallest zg, and ∼8% for both 300 K TH as well as 0.1 mm thGd. However, for all three plots (Figs. 10(a)–10(c)), the efficiency decreases as the parameters increase. Furthermore, the effective efficiency differences between 1.32 T and 2 T for these three parameters decrease as the parameters increase. These results clearly show the benefits of electing a higher magnetic field in conjunction with other parameters to maximize the efficiency outcomes.

Relative efficiency versus gap distance, heat source temperature, and Gd thickness for (a)–(c) the applied field study and (d)–(f) the contact conductance study
Comparatively, the magnitudes and trend of NORMAL versus ENHANCED conductance (Fig. 10(d)) for zg were nearly equal. For TH (Fig. 10(e)), NORMAL and ENHANCED had the same efficiency at 310 K but NORMAL was more efficient at temperatures below 310 K and ENHANCED was more efficient at temperatures above 310 K due to power generation. For thGd, NORMAL and ENHANCED had nearly the same efficiency for each thGd except 0.1 mm. Additionally, as observed in Fig. 10(f), each case reached a point where efficiency stopped decreasing as thGd increased. The minimum in Fig. 10(f) can be regarded as the minimum relative efficiency for both conductance cases. The minimum relative efficiency can be used to help curate design specifications based on a balance between compact design, desired power output, and efficiency. Overall, gap distance had a much larger influence on relative efficiency than TH and thGd. However, a higher applied field was more important to achieve higher efficiencies.
The parametric studies revealed the influence of key design features on device operation. The zg is the more dominant feature. As the zg reduces, more power is generated due to a reduction in travel time. The TH remains stationary as the gap distance changes. Therefore, we determined that the efficiency is dominated by the power. The culmination of each study framed an optimal device that operates close to its Curie temperature and consists of a smaller gap distance with a larger volume of Gd to achieve higher efficiencies. The amount of power reported in Fig. 9 is sufficient for low power devices such as wearables. Additionally, we found that the higher applied field is more important than the enhanced contact conductance. Though the enhanced contact conductance leads to higher frequencies, a higher applied field generates more power/efficiency in the thermomagnetic device. This differs from passive devices where higher contact conductance equates to higher frequencies and thus larger power generation. This is attributed to the coupling of the magnetic force directly to voltage induction. Whereas in a passive device, the mechanical response induces voltage. Further optimization will be achieved through incorporating hysteresis, additional ferromagnetic materials, atomistic thermal transition, and other losses.
Conclusion
This study explored the design aspects of an active thermomagnetic device used to transduce thermal energy to electrical energy. The base system consisted of a NdFeB hard magnet and a Gd soft magnet as the working material that oscillated between a heat sink and a heat source. A parametric study was conducted to understand the effects of varying the gap distance (zg), heat source temperature (TH), and Gd thickness (thGd). It was found that when selecting a design, the zg should be as small as possible and include the largest volume of Gd possible. A higher applied field will achieve favorable efficiencies but should be balanced with an enhanced contact conductance. Additionally, the Curie temperature of the chosen soft magnet should be close to the operating temperature. Future work will include determining the effect of spin reorientation on device efficiency and experimental work to fully assess the promise of this technology.
Footnote
Acknowledgment
This material was based upon work completed at the University of California Los Angeles and supported by the National Science Foundation (Grant No. 1306683) and the Center for Engineering Excellence and Diversity.