HUNGARIAN JOURNAL OF INDUSTRY AND CHEMISTRY Vol. 48(1) pp. 61–65 (2020) hjic.mk.uni-pannon.hu DOI: 10.33927/hjic-2020-09 NONLINEAR MAGNETIC PROPERTIES OF MAGNETIC FLUIDS FOR AUTOMOTIVE APPLICATIONS BARNABÁS HORVÁTH*1 AND ISTVÁN SZALAI2 1Institute of Physics and Mechatronics, University of Pannonia, Egyetem utca 10, Veszprém H-8200, HUNGARY 2Institute of Mechatronics Engineering and Research, University of Pannonia, Gasparich Márk utca 18/A, Zalaegerszeg H-8900, HUNGARY The external magnetic field required to activate a magnetic fluid in an industrial application is sufficiently large that mag- netization is no longer a linear function of the external field strength, i.e. magnetic fluids exhibit nonlinear characteristics. The aim of our research was to develop a measuring system which is capable of determining the nonlinear AC sus- ceptibility of magnetic fluids at discrete frequencies and in the presence of a high-intensity driving magnetic field. The measurement of susceptibility is based on the determination of the change in frequency of a low-intensity field, which is generated by an LC oscillator. The application of sinusoidal excitation to the material results in a variation in the suscep- tibility that modulates the frequency of the measured low-intensity field and in the appearance of higher-order harmonics of the driving field. The higher-order components of the nonlinear AC susceptibility are extracted from the measured response by Fourier analysis. By applying the measuring system, the nonlinear susceptibility of water-based ferrofluids (Ferrotec’s EMG 700) and its dependence on the magnetic field strength were investigated. Keywords: AC susceptometry, nonlinear susceptibility, ferrofluid 1. Introduction Regarding automotive applications, one of the most im- portant characteristics of magnetic fluids is their behav- ior in magnetic fields. The majority of these applications are for seals [1], vibration damping [2] and torque trans- mission [3]. In most cases, the magnitude of an external magnetic field required to activate a fluid is on the scale where the M magnetization is no longer a linear func- tion of the H field strength. This means that magnetic flu- ids show nonlinear characteristics, and the χ = ∂M/∂H susceptibility depends on the magnetic field strength. In the case of magnetic fluids, the nonlinearity is the result of two effects: normal saturation (alignment of the mag- netic dipole moments) and formation of structures (parti- cle chains). The latter influences the susceptibility as the magnetic field shifts the equilibrium between the struc- tures with different dipole moments (single particles and particle chains of different sizes). In a weak time-varying He(t) = He0 sin(ωt) mag- netic field (where t denotes the time and ω = 2πf rep- resents the angular frequency), the magnetization also changes periodically but lags behind the magnetic field because reorientation of the magnetic dipole moments is not instantaneous. In this case, the dynamic magnetic sus- ceptibility can be defined as a complex quantity χ∗(ω) = *Correspondence: bhorvath@almos.uni-pannon.hu χ′(ω) − iχ′′(ω). The real part, χ′, is related to the re- versible magnetization process and it is in-phase with the alternating field. Within the low-frequency limit, (f → 0) χ′ approaches the initial gradient of the steady-state mag- netization (initial DC susceptibility). If f → ∞, the re- orientation of the magnetic dipole moments cannot fol- low the alternating field and χ′ approaches zero. The imaginary, out-of-phase component, χ′′, is proportional to power losses due to energy absorption from the field and peaks at a characteristic frequency, fc. However, if the amplitude of the alternating field is sufficiently large, higher-order harmonics appear in the magnetization due to the nonlinear characteristics. In this case, the real part of the alternating current (AC) mag- netic susceptibility is: χ′ = χ0 + χ2ω cos(2ωt) + χ4ω cos(4ωt) + . . . , (1) where χ0 denotes the base component, and χ2ω and χ4ω represent the amplitude of the second- and fourth-order harmonics, respectively [4]. A similar equation holds for the imaginary component, χ′′. In this work, only the har- monics of the real part are considered, because the non- linear susceptibility at frequencies much higher than the characteristic frequency of the magnetic fluid is inves- tigated, where the imaginary part is very close to zero. By taking magnetic fluids into consideration, the ampli- tudes of sixth- and higher-order harmonics are so small https://doi.org/10.33927/hjic-2020-09 mailto:bhorvath@almos.uni-pannon.hu 62 HORVÁTH AND SZALAI PC LabVIEW MD analyzer oscillator C DAQ GPIB power f(t) L amplifier fm = 957 kHz Hm0 = 0.022 kA/m measuring field Helmholtz coil fe = 1 - 200 Hz He0 = 0.2 - 3.0 kA/m driving field He(t) Hm(t) FFT χ2 χ4 Figure 1: A block diagram of the nonlinear magnetic measuring system. that they can be neglected. In alternating magnetic fields, only the even-order harmonics appear, because the mag- netization of magnetic fluids exhibits inversion symme- try with respect to a change in the direction of H. Should a symmetry-breaking DC bias field be superimposed on the AC magnetic field, both the odd- and even-order har- monics appear. As the nonlinearity is caused partly by structure formation, the components of the nonlinear sus- ceptibility and their ratios are sensitive to changes in the structure of magnetic fluids. Numerous studies have discussed the theoretical de- scription of the nonlinear susceptibility of magnetic fluids [4–7], but experimental results are rather scarce. Never- theless, nonlinear magnetic susceptibility measurements are quite useful for the study of magnetic fluids and pro- vide additional information besides linear AC suscep- tometry. If the linear response to small-amplitude fields is measured, the relaxation processes of the magnetic dipoles can be studied. However, by applying the non- linear method, it is possible to obtain information about the structural evolution of an activated magnetic fluid. The aim of our work was to develop a measuring sys- tem which is capable of determining the nonlinear AC susceptibility of magnetic fluids at discrete frequencies in the presence of external magnetic fields. The dependence of the 2nd- and 4th-order nonlinear AC susceptibility of ferrofluids on the magnetic field strength was investigated by the developed susceptometer. 2. Experimental 2.1 Nonlinear measurement method and setup The measurement method of susceptibility is based on our nonlinear dielectric system [8], which is adapted to magnetic measurements. With this method, the real part of complex AC susceptibility is determined from the change in frequency of a low-intensity measuring field (Hm). An air core solenoid (L) was filled with the sam- ple of fluid, and the inductance determined the frequency of the sinusoidal measuring field generated by an LC os- cillator (Fig. 1). Therefore, any change in the susceptibil- ity of the sample caused the resonance frequency fm to shift. The susceptibility of the sample was measured at this base frequency. Inside the solenoid, the amplitude of the measuring field is so small (Hm0 = 0.022 kA/m) that within this region, the magnetization curve of the fluid can be regarded as linear, thus no significant structural change can be expected. The oscillator which generated the measuring field was a Colpitts-type parallel LC circuit. The active feed- back element was a double triode vacuum tube (ECC88). The capacitor bank C was composed of high-quality sil- ver mica capacitors with a low temperature coefficient. The measuring coil L was connected in parallel with the capacitor block, and by changing the coil, the base fre- quency of the oscillator could be varied. With different inductors, the measuring frequency was set to discrete values, namely 153 kHz, 590 kHz and 957 kHz. The di- mensions of the measuring coils were identical: 25 mm in length with an inner diameter of 7.1 mm. They were composed of enameled copper wires of different sizes and wound on coil formers made of plexiglass. The reso- nance frequency of the LC oscillator was measured by a Hewlett Packard 53310A Modulation Domain (MD) An- alyzer. The relationship between the resonance frequency and susceptibility of the sample was determined by cali- bration using different materials of known susceptibility. A high-intensity driving field (He) was generated by a pair of Helmholtz coils which were placed around the measuring coil. The axis of the Helmholtz coil and, there- fore, the direction of the field, were parallel to the axis of the solenoid. The upper limit on the amplitude of the driving magnetic field with the current setup was He0 = 6.4 kA/m, which was two orders of magnitude greater than the amplitude of the measuring magnetic field. The maximum of He0 depended on the frequency of the field: as the frequency increased, the maximum of He0 decreased. The uniformity of the magnetic field strength in the volume of the sample was better than 1 %. The Helmholtz coil was driven by a high-current function generator, which consisted of a Labworks PA-138 linear Hungarian Journal of Industry and Chemistry NONLINEAR MAGNETIC PROPERTIES OF MAGNETIC FLUIDS 63 power amplifier and a signal source. The input signal of the amplifier was provided by a multifunction data ac- quisition (DAQ) card (National Instruments PCI-6052E). Any arbitrary waveform could be generated by the signal generator, and the waveform offset by a DC value. Under the influence of the high-intensity driving field, the susceptibility of the sample changes, therefore, the frequency of the measuring field was modulated. If the driving field is sinusoidal, then the time-domain sus- ceptibility response will contain the higher-order har- monics of the field. The 2nd- and 4th-order components of the nonlinear AC susceptibility were extracted from the measured response by Fourier analysis. The Fast Fourier Transform (FFT) algorithm was implemented us- ing custom-developed LabVIEW software, which pro- vided control and data acquisition functions of the mea- suring system. 2.2 Supplemental measurements During the nonlinear measurements, the change in sus- ceptibility relative to the zero-field susceptibility was measured, so the real part of the AC susceptibility of the fluid at the base frequency had to be determined by an- other method. For this purpose, a spectrum was obtained within the frequency range of 200 Hz - 1 MHz. The ini- tial DC susceptibility of the ferrofluid was determined from the DC magnetization curve. The AC susceptibility in zero field was measured by an inductive method where the inductance of a solenoid filled with the sample was measured by an impedance analyzer (Agilent 4284A). By determining the impedances of the empty air core solenoid and when it is filled with the sample, χ′ was cal- culated. The solenoid used for these measurements was the same as the measuring coil of the nonlinear setup. 2.3 Material By applying the nonlinear magnetic measuring system, the nonlinear properties of Ferrotec’s EMG 700 ferrofluid were investigated. This material is water-based and con- tains magnetite particles with a nominal diameter of ∼ 10 nm. The volume concentration of the particles is 5.8 %(v/v). The fluid is stabilized by an anionic surfactant. For the measurements, glass-tube sample holders with an inner diameter of 3.1 mm were filled with the ferrofluid. The length of the tube extended beyond the length of the measuring coil at both ends. The volume of the sample was 0.49 cm3. 3. Results and Discussion According to the DC magnetization curve, the initial DC susceptibility of the EMG 700 ferrofluid was χDC = 12.57. This is the limiting value of the real part of the complex susceptibility if f approaches zero. The AC sus- ceptibility spectrum of the ferrofluid is shown in Fig. 2. Relaxation of χ′ was observed within the lower frequency 10-1 1 10 102 103 104 105 106 0 3 6 9 12 15 c = 0.54 cs = 12.55 f (Hz) c' He = 0 fc = 39 Hz Figure 2: Relaxation of the real part of AC susceptibility in the absence of a driving magnetic field (hollow sym- bols) by applying the Cole-Cole equation (solid line) to- gether with the estimated parameters of the relaxation. region. The experimental frequency-dependent suscepti- bility data was fitted by the Cole-Cole equation (solid line in Fig. 2). This relaxation model is suitable for describing the response of dipoles to an alternating field should the relaxation not be ideal (e.g. if instead of a single relax- ation time, a distribution of relaxation times exists). The fitting process yielded the static (f → 0) χs = 12.55 ± 0.07 and infinite frequency (f → ∞) χ∞ = 0.54±0.02, susceptibilities and the so-called central char- acteristic time of the relaxation (τ0). The reciprocal of the characteristic time yielded the characteristic frequency fc = 1/(2πτ0), which is 39 ± 5 Hz for the ferrofluid EMG 700. The spectrum shows that relaxation occurs at much lower frequencies than the base frequency of the nonlinear susceptibility measurements (fc = 39 Hz vs. fm = 957 kHz). At the high-frequency end (∼1 MHz) of the spectrum, χ′ decreased to ∼ 0.8 and changed slightly as the frequency increased in this region. At such a large distance from fc, the imaginary part of the complex sus- ceptibility was close to zero, thus it is justified to consider only the real part (see Eq. 1). To investigate the nonlinear AC susceptibility, sinu- soidal excitation was applied. The frequency of the driv- ing magnetic field was changed from 1 Hz to 200 Hz with an amplitude as high as He0 = 3 kA/m (depend- ing on fe). The nonlinear AC susceptibility was mea- sured at 957 kHz in all cases. At He0 = 0, the real part of the AC susceptibility of the ferrofluid at the measure- ment frequency of fm = 957 kHz was χ′ = 0.78. Fig. 3 shows a typical susceptibility response of the ferrofluid EMG 700 at an excitation frequency of fe = 1 Hz and using different magnetic field strengths. It can be seen that the response contained higher-order harmonics. The magnetic field always caused a decrease in the suscepti- bility (because of the saturation) regardless of its direc- tion. This inversion symmetry was expressed by the fact that the susceptibility response during the first half period 48(1) pp. 61–65 (2020) 64 HORVÁTH AND SZALAI 0.0 0.5 1.0 1.5 2.0 0.63 0.66 0.69 0.72 0.75 0.78 0.81 c'(t) He0 = 0.7 kA/m He0 = 1.1 kA/m He0 = 1.4 kA/m t (s) c' He(t) -1 0 1 2 3 4 5 6 7 H e (k A /m ) Figure 3: The susceptibility response of the ferrofluid EMG 700 at various driving field amplitudes (fe = 1 Hz). of He(t) was the same as during the second half even if the direction of the field was reversed. Therefore, the maxima of the susceptibility responses were always 0.78, which corresponded to the point when He(t) intersected the value zero. The dependence of the 2nd- and 4th-order harmonics (χ2ω and χ4ω in Eq. 1) on the magnetic field strength ex- tracted from the measured response is shown in Figs. 4a and 4b. At a given driving field strength, χ2ω was one or- der of magnitude larger than χ4ω. The amplitude of both components increased as the magnetic field strength rose (at a given driving field frequency). This occurred be- 0.00 0.02 0.04 0.06 0.08 0.10 0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.000 0.005 0.010 0.015 c 2 w fe = 1 Hz fe = 40 Hz fe = 100 Hz fe = 200 Hz a b c 4 w He0 (kA/m) Figure 4: The dependence of the components χ2ω (a) and χ4ω (b) of the nonlinear susceptibility of the ferrofluid EMG 700 on the amplitude of the driving field at different driving field frequencies. cause at higher magnetic field strengths, the nonlinearity of the magnetization curve increased as the magnetiza- tion approached the saturation level. 4. Conclusions A nonlinear magnetic measuring system was developed to determine the nonlinear AC susceptibility of magnetic fluids. The system was applied to measure the compo- nents of the nonlinear susceptibility of the water-based ferrofluid EMG 700 up to the 4th order. It has been shown that by using the aforementioned system, it is possible to determine the dependence of the nonlinear susceptibility on the magnetic field strength. If the excitation frequency is varied, the relaxation of the nonlinear susceptibility can be also studied. Our future aim is to investigate this re- laxation and other types of ferrofluids (based on organic carrier and polydisperse fluids), moreover, to study the influence of the magnetite concentration on the nonlinear magnetic properties. Acknowledgments This project was supported by the European Union and co-financed by the European Social Fund. EFOP-3.6.2- 16-2017-00002. REFERENCES [1] Ravaud, R.; Lemarquand, G.: Mechanical Proper- ties of a Ferrofluid Seal: Three-Dimensional An- alytical Study based on the Coulombian Model, Prog. Electromagn. Res. B, 2009, 13, 385–407, DOI: 10.2528/PIERB09020601 [2] Sassi, S.; Cherif, K.; Mezghani, L.; Thomas, M.; Kotrane, A.: An innovative magnetorheological damper for automotive suspension: from design to experimental characterization, Smart Mater. Struct., 2005, 14(4), 811–822, DOI: 10.1088/0964- 1726/14/4/041 [3] Rizzo, R.: An innovative multi-gap clutch based on magneto-rheological fluids and electrodynamic ef- fects: magnetic design and experimental character- ization, Smart Mater. Struct., 2017, 26(1), 015007, DOI: 10.1088/0964-1726/26/1/015007 [4] Wang, G.; Huang, J.P.: Nonlinear magnetic sus- ceptibility of ferrofluids, Chem. Phys. Lett., 2006, 421(4), 544–548, DOI: 10.1016/j.cplett.2006.02.010 [5] Huke, B.; Lücke, M.: Magnetic properties of col- loidal suspensions of interacting magnetic parti- cles, Rep. Prog. Phys., 2004, 67(10), 1731–1768, DOI: 10.1088/0034-4885/67/10/R01 [6] Szalai, I.; Nagy, S.; Dietrich, S.: Linear and non- linear magnetic properties of ferrofluids, Phys. Rev. E, 2015, 92(4), 042314, DOI: 10.1103/Phys- RevE.92.042314 Hungarian Journal of Industry and Chemistry https://doi.org/10.2528/PIERB09020601 https://doi.org/10.1088/0964-1726/14/4/041 https://doi.org/10.1088/0964-1726/14/4/041 https://doi.org/10.1088/0964-1726/26/1/015007 https://doi.org/10.1016/j.cplett.2006.02.010 https://doi.org/10.1088/0034-4885/67/10/R01 https://doi.org/10.1103/PhysRevE.92.042314 https://doi.org/10.1103/PhysRevE.92.042314 NONLINEAR MAGNETIC PROPERTIES OF MAGNETIC FLUIDS 65 [7] Ivanov, A.O.; Kuznetsova, O.B.: Magnetic prop- erties of dense ferrofluids, J. Magn. Magn. Mater., 2002, 252, 135–137, DOI: 10.1016/S0304- 8853(02)00639-X [8] Horváth, B.; Szalai, I.: Structure of electrorhe- ological fluids: A dielectric study of chain for- mation, Phys. Rev. E, 2012, 86(6), 061403, DOI: 10.1103/PhysRevE.86.061403 48(1) pp. 61–65 (2020) https://doi.org/10.1016/S0304-8853(02)00639-X https://doi.org/10.1016/S0304-8853(02)00639-X https://doi.org/10.1103/PhysRevE.86.061403 https://doi.org/10.1103/PhysRevE.86.061403 Introduction Experimental Nonlinear measurement method and setup Supplemental measurements Material Results and Discussion Conclusions