Abstract
The paper presents a new approach to modeling of magnetic properties and the influence of temperature on hysteresis loops in La(FeCoSi)13 magnetocaloric compound. The GRUCAD hysteresis model has been used for analysis of magnetic behavior in the vicinity of the magnetic phase transition. The modeling results are in a qualitative agreement with experiment.
1 Introduction
Modern cooling based on magnetocaloric effect is an environment-friendly technique used to lower the temperature. The well-known magnetocaloric materials are pure gadolinium and its alloys, especially Gd5Ge2Si2 discovered by Pecharsky and Gschneidner Jr. in 1997 [1]. Due to the fact that Gd is lanthanide, its price is relatively high. Relatively cheap pseudo-binary La(FeCoSi)13 compounds have been the subject of considerable interest of engineering community as the magnetocaloric effect occurs in these materials near room temperature [2, 3]. This group of alloys is based on face centered cubic NaZn13 structure and it is responsible for excellent magnetic properties. The La atoms occupy 8a (in Wyckoff notation) positions and Fe atoms two non-equivalent positions (8b and 96i). The La atoms may be substituted by other lanthanides (Ce [4], Pr or Ho [5]). Moreover, the Fe atoms placed in position 96i may be substituted by i.e. Co [6], Ni [7] or Al [8]. Such modifications of chemical composition have an important influence on the value of the Curie point and change in magnetic entropy.
Most of the papers concerning the La(Fe,Si)13 alloys deliver information about their processing, microstructure and magnetocaloric properties. From the scientific point of view, it is also important to investigate various effects in the vicinity of phase transition observed in the material. An interesting technique used to study the response of magnetic material near the phase transition is modeling of hysteresis curves [9].
The present paper considers the possibility to model the hysteresis curves of La(FeCoSi)13 compound near the transition temperature with the use of a phenomenological description. It is an extension of the previous work [10], where the variations of shape of the anhysteretic (hysteresis-free) curves were considered.
The paper is structured as follows. Section 2 provides basic pieces of information on contemporary hysteresis models used in studies of magnetocaloric effect and the GRUCAD description applied for this purpose in this paper. Section 3 discusses in detail the measurement methods and results. Section 4 is devoted to modeling. Finally, Section 5 comments on the obtained results.
2 Phenomenological hysteresis models
Hysteresis modeling may be carried out with an appropriate description like the bottom-up approach by Preisach, extended by Mayergoyz [11] or the phenomenological proposal advanced by Jiles and Atherton [12]. The latter model has attracted much attention of the scientific community in the last thirty years due to the uncomplicated structure of its equations and physical meaning attributed to its parameters. The Jiles-Atherton (JA) model still remains one of the most used descriptions for modeling magneto-structural first order transitions in magneto-caloric materials as well [13]. However, some authors have pointed out a number of its drawbacks, to mention the necessity to “patch” the model equations with an artificially introduced term suppressing the irreversible magnetization term after a sudden field reversal [14], non-closure of minor hysteresis loops if the values of model parameters are kept constant [15], the ambiguities due to the use of modified Langevin function as the description of the anhysteretic curve [16, 17] (“modified” means that the argument of the function includes the mean field term).
The model developed by the GRUCAD research group [18, 19] is a valuable alternative to the original JA formulation. In the GRUCAD description the reversible and irreversible processes are decoupled. The irreversibility is implemented by offsetting the loop branches along the H-axis from the anhysteretic (truly reversible in thermodynamic sense) curve, given with Langevin function. There is no need to introduce additional conditions, which are aimed at the modification of model equations after a field reversal. This description is similar in spirit to some models consistent with the laws of irreversible thermodynamics [20, 21, 22, 23]. The equations that define the model are as follows:
Here δ = sign (dB/dt) = ± 1, whereas α, a, γ, HHS and Ms are model parameters. Their values may be obtained using any robust optimization routine. Two first equations refer to the reversible, whereas the third and the fourth one – to the irreversible magnetization processes.
3 Material and measurements
Magnetocaloric bulk samples of La(FeCoSi)13 applied for the tests are made of reactively sintered fine-grained powders manufactured in Vacuumschmelze Company (Figure 1) [24]. Curie temperature Tc and peak entropy changes ΔSM(T) can be tailored through introduction of cobalt into structure La(FexCoySi1−x−y)13 [3, 24]. The study and modeling were carried out with composition tuned-up for the Curie point of 300 K (Figure 2).
The final-form of plates with overall dimension 18 × 36 × 5 mm has been obtained in machining process [24]. Nominal parameters of examined composition are comprised in Table 1. Temperature and magnetic dependencies of magnetic entropy changes |ΔSM(T)| are depicted in Figure 2.
Parameter | Designation and value |
---|---|
Composition | La(Fe10.8Co1.1Si1.1) |
Max. entropy changes (at 111kA/m) | |ΔSM(T)| = 44.7 kJ/m3K |
Curie temperature | Tc = 300 K |
Temperature of peak entropy | Tpeak = 296 K at H = 111 kA/m |
Volumetric mass density | ρ = 7288 kg/m3 |
The experimental study of magnetic properties concerns hysteresis loops, magnetic polarization and power losses in the vicinity of Curie temperature. The presented results of measurements have been obtained with application of IEC 60404-6:2004 method. The reference measurements of power losses have been performed simultaneously with the use of Unbalanced Bridge Method (UBM) [25]. Instrumentation and primal uncertainties of the measurement methods are presented in Table 2.
Uncertainty | Designation and value |
---|---|
H(t) B-type uncertainty | uBREL(H) = 0.4 |
J(t) B-type uncertainty | uBREL(J) = 2.5 |
UUBM(t) B-type uncertainty | uBREL(J) = 2.4 |
Ps std. rel. uncertainty | uBREL(PS) = 5 |
T std. uncertainty | u(T) = 0.25 (K) |
Measurements of magnetic waveforms H(t), J(t), uB(t) and specific power losses PS have been carried out in a closed magnetic core. The core was made by setting of four homogeneous LaFe10.8Co1.1Si1.1 (total mass = 0.094 kg) plates in non-magnetic yoke. Each component of the core was equipped with a section of pickup coil W2 (4 × 150 turns) and exciting coil W1 (4 × 72 turns) (see Figure 4). Comprehensive description of the measurement setup was presented previously in Ref. [26], whereas the most significant uncertainty components are collected in Table 2.
Investigations of thermal effects on magnetic properties and power losses were evaluated through analysis of the magnetic hysteresis loops J(H) measured at the frequency of 0.1 Hz (quasi-static conditions - cf. Figure 3).
Transitional stage of magnetic coercive field HC, remanent magnetization Jr and saturation Jmax of La(FeCoSi)13 core are shown in Figure 4. The measurements confirm a strong influence of the temperature on hysteresis loops of La(FeCoSi)13 compound.
The measured hysteresis loops (Figure 3) and the magnetization curve Jmax(T) in Figure 4 clearly indicate a transition from ferromagnetic to paramagnetic regime close to 300 K.
4 Modeling
In the previous work [10], we have focused on the description of the anhysteretic curve near the transition temperature. The shape parameter in the Langevin function used in modeling was found to be directly related to temperature. This fact allowed us to obtain a unified description for anhysteretic curves for temperatures both below and above the transition point. For temperatures above 300 K the hysteresis loops practically disappeared and the slope of the anhysteretic curve changed abruptly. In the present paper we focus on modeling hysteresis curves for temperatures below 300 K.
The values of GRUCAD model parameters were determined for successive temperatures, keeping the values of fixed to those determined previously [10]. The results of modeling are presented in Figure 5.
In order to check how well the GRUCAD model is able to describe the hysteresis loops of La(FeCoSi)13 compound near the transition temperature, we have compared the measured and the modeled values of coercive field strength and remanence induction for quasi-static magnetization conditions (f = 0.1 Hz). Moreover we have integrated numerically the area of the modeled hysteresis loops in order to find out the approximate value of power loss density. The results for chosen temperatures are summarized in Table 3.
Parameter | ||||||
---|---|---|---|---|---|---|
Temp. | Hcmeas. (A/m) | Hcmodel (A/m) | Brmeas. (mT) | Brmodel (mT) | Psmeas. (mW/kg) | Psmodel (mW/kg) |
295.5 (K) | 47.3 | 39.7 | 38.7 | 74.7 | 1.26 | 1.37 |
299.5 (K) | 33.4 | 33.8 | 24.0 | 25.4 | 0.48 | 0.50 |
It can be stated, that the measured and modeled values of coercive field strength and power loss density are in qualitative agreement. The maximum relative error understood as the absolute value of the difference between the measured and the modeled values of the quantity referred to the measured value does not exceed 17.1%. However it can be noticed that the model overestimated significantly the remanence value at 295.5 K. Figure 6 depicts the hysteresis loops for temperatures considered in Table 3 in a more detailed way.
5 Conclusions
The presented modeling results clearly indicate that the GRUCAD model might be a valuable tool supporting prediction of the effects of temperature variations in La(FeCoSi)13 magnetocaloric compound near the transition temperature. The obtained results are in qualitative agreement with the experiment and for coercive field strength and power loss density there is also an excellent quantitative agreement.
Future work shall be aimed at fine-tuning of the GRUCAD description. In the present contribution it was assumed that the compound may be described with a single set of model parameters. It should be born in mind that in reality magneto-caloric materials exhibit a highly complicated internal structure due to the presence of multiple phases. Therefore a more adequate description should include distributions of model parameters to take this fact into account.
Acknowledgement
This research was partly supported by the grants No. 6370/B/T02/2011/40 and 4889/B/T02/2011/40, from the National Science Centre of Poland.
References
[1] Tishin A.M., Spichkin Y.I., Themagnetocaloric effect and its application, Institute of Physics Publishing, London, 200310.1887/0750309229Search in Google Scholar
[2] Liu J., Moore J. D., Skokov K. P., Krautz M., Löwe K., Barcza A., et al., Exploring La(Fe,Si)13-based magnetic refrigerants towards application, Scripta Mat., 2012, 67, 584-589.10.1016/j.scriptamat.2012.05.039Search in Google Scholar
[3] M’nassri R., Cheikhrouhou A., Magnetocaloric effect in LaFe10.7Co0.8Si1.5 compound near room temperature, J. Supercond. Nov. Mater., 2014, 27, 1059-1064.10.1007/s10948-013-2375-1Search in Google Scholar
[4] Fujita A., Fujieda S., Fukamichi K., Realtive cooling power of La(FexSi1−x)13 after controlling the Curie temperature by hydrogenation and partial substitution of Ce, J. Magn. Magn. Mat., 2007, 310, e1006-e1007.10.1016/j.jmmm.2006.10.962Search in Google Scholar
[5] Gębara P., Kovac J., Magnetocaloric effect of the LaFe11.2Co0.7Si1.1 modified by partial substitutiton of La by Pr or Ho, J. Mat. Des., 2017, 129, 111-115.10.1016/j.matdes.2017.05.038Search in Google Scholar
[6] Yan A., Muller K.H., Gutfleisch O., Magnetocaloric effect in the LaFe11.8−xCoxSi1.2 melt-spun ribbons, J. All. Compd., 2008, 450, 18-21.10.1016/j.jallcom.2006.10.093Search in Google Scholar
[7] Gębara P., Pawlik P., Hasiak M., Alteration of negative lattice expansion of the La(Fe,Si)13-type phase in LaFe11.14−xCo0.66NixSi1.2 alloys, J. Magn. Magn. Mat., 2017, 422, 61-65.10.1016/j.jmmm.2016.08.068Search in Google Scholar
[8] Gębara P., Pawlik P., Michalski B., Wysłocki J.J., Measurements of magnetocaloric effect in LaFe11.14Co0.66Si1.2-xAlx (x = 0.1, 0.2, 0.3) alloys, Acta Phys. Pol. A, 2015, 127, 576-578.10.12693/APhysPolA.127.576Search in Google Scholar
[9] Basso V., Kuepferling M., Sasso C. P., LoBue M., Modeling hysteresis of first-order magneto-structural phase transformations, IEEE Trans. Magn., 2008, 44, 3177-3180.10.1109/TMAG.2008.2002796Search in Google Scholar
[10] Gozdur R., Chwastek K., Najgebauer M., Lebioda M., Bernacki Ł., Wodzynski A., Scaling of anhysteretic curves for LaFeCoSi alloy near the transition point, Acta Phys. Pol. A, 2017, 131, 801-803.10.12693/APhysPolA.131.801Search in Google Scholar
[11] Mayergoyz I.D., Mathematical models of hysteresis and their applications, 2nd Ed., Elsevier Academic Press, Amsterdam, 200310.1016/B978-012480873-7/50005-0Search in Google Scholar
[12] Jiles D.C., Atherton D.L., Theory of ferromagnetic hysteresis, J. Magn. Magn. Mater., 1986, 61, 48-60.10.1016/0304-8853(86)90066-1Search in Google Scholar
[13] Melikhov Y., Hadimani R.L., Raghunathan A., Phenomenological modelling of first order transitions in magnetic systems, J. Appl. Phys., 2014, 115, 18390210.1063/1.4875678Search in Google Scholar
[14] Steentjes S., Chwastek K., Petrun M., Dolinar D., Hameyer K., Sensitivity analysis and modeling of symmetric minor hysteresis loops using the GRUCAD description, IEEE Trans. Magn., 2014, 50, 730080410.1109/TMAG.2014.2323250Search in Google Scholar
[15] Benabou A., Leite J. V., Clénet S., Simão C., Sadowski N., Minor loops modelling with a modified Jiles-Atherton model and comparison with the Preisach Model, J. Magn. Magn. Mater., 2008, 320, e1034-e1038.10.1016/j.jmmm.2008.04.092Search in Google Scholar
[16] Zirka S. E., Moroz Yu. I., Harrison R. G., Chwastek K., On physical aspects of Jiles-Atherton models, J. Appl. Phys., 2012, 112, 04391610.1063/1.4747915Search in Google Scholar
[17] Qingyou Liu, Xu Luo, Haiyan Zhu, Yiwei Han, Modified magnetomechanical model in the constant and low intensity magnetic field based on J-A theory, Chinese Phys. B, 2017, 26, 07750210.1088/1674-1056/26/7/077502Search in Google Scholar
[18] Koltermann P. I., Righi L. A., Bastos J. P. A., Carlson R., Sadowski N., Batistela N. J., A modified Jiles method for hysteresis computation including minor loops, Physica B, 2000, 275, 233-237.10.1016/S0921-4526(99)00770-XSearch in Google Scholar
[19] Righi L. A., Sadowski N., Carlson R., Bastos J. P. A., Batistela N. J., A new approach for iron losses calculation in voltage fed time stepping Finite Elements, IEEE Trans. Magn., 2001, 37, 3353-3356.10.1109/20.952612Search in Google Scholar
[20] Bergqvist A., Magnetic vector hysteresis model with dry friction-like pinning, Physica B, 1997, 233, 342-347.10.1016/S0921-4526(97)00319-0Search in Google Scholar
[21] Henrotte F., Nicolet A., Hameyer K., An energy-based vector hysteresis model for ferromagnetic materials, COMPEL, 2006, 25, 71-80.10.1108/03321640610634344Search in Google Scholar
[22] Tavakoli H., Bormann D., Ribberfjärd D., Engdahl G., Comparison of a simple and a detailed model of magnetic hysteresis with measurements on electrical steel, COMPEL, 2009, 28, 700-710.10.1108/03321640910940954Search in Google Scholar
[23] François-Lavet V., Henrotte F., Stainer L., Noels L. , Geuzaine C., An energy-based variational model of ferromagnetic hysteresis for finite element computations, J. Comp. Appl. Math., 2013, 246, 243-250.10.1016/j.cam.2012.06.007Search in Google Scholar
[24] Katter M., Zellmann V., Reppel G.W., Uestuener K., Magnetocaloric Properties of La(Fe,Co,Si)13 Bulk Material Prepared by Powder Metallurgy, IEEE Trans. Magn., 2008, 44, 3044-3047.10.1109/TMAG.2008.2002523Search in Google Scholar
[25] Gozdur R., Majocha A., Power losses measurements of nanocrystalline and amorphous magnetic cores, Electr. Machin. - Trans. J., 2013, 100, 175-179.Search in Google Scholar
[26] Gozdur R., Lebioda M., Bernacki Ł., Power losses in LaFexCoySi1.1 intermetallics near the magnetic phase transition, Acta Phys. Pol. A., 2015, 128, 98-103.10.12693/APhysPolA.128.98Search in Google Scholar
© 2018 Roman Gozdur et al., published by De Gruyter
This work is licensed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 License.