Next Article in Journal
A Clean Process for Obtaining High-Quality Cellulose Acetate from Cigarette Butts
Previous Article in Journal
The Impact of Hemp Shives Impregnated with Selected Plant Oils on Mechanical, Thermal, and Insulating Properties of Polyurethane Composite Foams
Previous Article in Special Issue
Weak Antilocalization Tailor-Made by System Topography in Large Scale Bismuth Antidot Arrays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Determination of Q Quality Factor and Resonance Frequency in Sensors Based on the Magnetoelastic Resonance Through the Fitting to Analytical Expressions

1
Departamento de Electricidad y Electrónica, Universidad del País Vasco/Euskal Herriko Unibertsitatea (UPV/EHU), Barrio Sarriena s/n, 48940 Leioa, Spain
2
BCMaterials, Basque Center for Materials, Applications and Nanostructures, UPV/EHU Science Park, 48940 Leioa, Spain
3
Departamento de Matemática Aplicada y Estadística e Investigación Operativa, Universidad del País Vasco UPV/EHU, P.O. Box 644, 48080 Bilbao, Spain
*
Author to whom correspondence should be addressed.
Materials 2020, 13(21), 4708; https://doi.org/10.3390/ma13214708
Submission received: 17 September 2020 / Revised: 15 October 2020 / Accepted: 19 October 2020 / Published: 22 October 2020
(This article belongs to the Special Issue Advances in Magnetic Materials and Magneto-Elastic Sensors)

Abstract

:
The resonance quality factor Q is a key parameter that describes the performance of magnetoelastic sensors. Its value can be easily quantified from the width and the peak position of the resonance curve but, when the resonance signals are small, for instance when a lot of damping is present (low quality factor), this and other simple methods to determine this parameter are highly inaccurate. In these cases, numerical fittings of the resonance curves allow to accurately obtain the value of the quality factor. We present a study of the use of different expressions to numerically fit the resonance curves of a magnetoelastic sensor that is designed to monitor the precipitation reaction of calcium oxalate. The study compares the performance of both fittings and the equivalence of the parameters obtained in each of them. Through these numerical fittings, the evolution of the different parameters that define the resonance curve of these sensors is studied, and their accuracy in determining the quality factor is compared.

1. Introduction

Magnetoelastic resonance sensors are usually made of amorphous, ribbon-shaped ferromagnetic alloys [1]. Magnetostriction in these materials provides a strong coupling of their magnetic and mechanical properties, so that when an external alternating magnetic field is applied to a magnetoelastic ribbon, elastic waves are induced in it. Conversely, the magnetic state of these materials is highly dependent on the application of external mechanical forces and loads, and therefore these mechanical variations generate changes in the magnetic flux that can be remotely detected by a detection coil placed on the ribbon. Matching the physical dimensions and elastic properties of the material used in the sensor, the phenomenon of magnetoelastic resonance occurs at specific frequencies that can be expressed as:
f n = n 2 L E ρ
where L is the length of the ribbon, E is the Young’s modulus of the material, ρ its density, and n = 1 corresponds to the fundamental resonant frequency. This resonant behavior is highly sensitive to different external parameters, which has generated a great interest in the use of these materials as a main part of different sensing systems [2]. Magnetoelastic sensors have been used to measure several environmental parameters such as pressure, humidity or temperature [3,4], liquid viscosity and density [5,6], chemical agents and pH [7,8], and biological agents [9,10,11]. Specifically, and with respect to its use as a mass sensor, when the magnetoelastic ribbon is loaded with a uniform and small amount of mass ( δ m ), a shift in its resonance frequency occurs, that can be evaluated from Equation (1) to be, to first order approximation:
δ f = 1 2 δ m m 0 f 0
where f 0 is the resonance frequency of the bare material and m 0 the initial mass of the sensor.
The performance of this type of mass sensor is mainly determined by two parameters: the sensitivity of the sensor ( S ) and the quality factor ( Q ) [12]. The sensitivity indicates the minimum detectable frequency change when the mass of the sensor increases:
S = δ f δ m
The quality factor Q is a dimensionless parameter that describes how sharp the resonance curve is, and decreases when the resonator is damped. Damping influences the shape of the magnetoelastic resonance curve and the resonance frequency itself, considered as the frequency of the maximum of the curve. A high Q value indicates a lower rate of energy loss, a sharp magnetoelastic resonance peak (Figure 1) and therefore, good resolution when determining the resonance frequency and frequency shifts. For developing high performance sensors, the highest S and Q are desirable.
The quality factor can be estimated from the resonance curve as the experimentally measured resonance frequency f m a x (at which the amplitude is maximum, A m a x ), divided by the width of the signal Δ f (measured as the full width at half maximum, corresponding to amplitudes A m a x / 2 ) [13]:
Q = f m a x Δ f
But the accuracy of this estimation method is very low (it can lead to errors up to 20 % in the determination of Q ) [14], so finding alternative and robust methods to determine this parameter with greater accuracy is of great importance for the study of the performance of these sensors, especially in cases where the curves have considerable noise or damping, as it is the case when the magnetoelastic sensor operates immersed in a fluid. Some studies have been made on the different methods of determining this quality factor Q [13,15], and Lopes et al. [16] reported an extensive study of the determination of the quality factor of a magnetoelastic sensor using different strategies. In this last work, a numerical fitting of the experimental data to an expression of the magnetic susceptibility of the sample around the resonance is carried out and used as a suitable technique to calculate the value of Q .
In the present work, the numerical fittings of the resonance curves obtained with a magnetoelastic sensor have been carried out using two different expressions: on the one hand, the expression used by Lopes et al. [16], that appears as Equation (6) below in this work and, on the other hand, a phenomenological approach that describes the frequency response of these magnetoelastic sensors and that has already been used to fit these resonance curves [17,18], appearing as Equation (9) within this work. The performance of each of the fitting expressions to determine the quality factor Q has been studied, and their comparison and equivalence carried out.
The resonance curves used in the fit analysis we will show in the following, correspond to those obtained in the monitoring of the precipitation reaction of calcium oxalate salt crystals (inorganic salts resulting from various metabolic activities in humans) using a magnetoelastic sensor. This monitoring process was first described by Bouropoulos et al. [19], and the measurements obtained in this reaction sensing and used in the present work are explained in detail in a previous publication [20]. Through the numerical fitting of these curves, we have studied the evolution of the different parameters that define the resonance curve as the precipitation reaction progress, which provides a better understanding of the effect that the changes of mass in the sensor and damping have on the resonance signal obtained.

2. Materials and Methods

2.1. Magnetoelastic Material of the Sensor

The material used as sensor platform is a magnetoelastic ribbon of composition Fe 73 Cr 5 Si 10 B 12 provided by Vacuumschmelze GmbH & Co. KG, Hanau Germany. This material has an excellent corrosion resistance behavior due to the fact that it contains a small amount of chromium ( 5 %   at . ), which forms a superficial passivation layer in the material. Magnetic, magnetoelastic and corrosion resistance properties of this material were studied by Sagasti et al. [21], and are summarized in Table 1. The strip was laser cut with dimensions 20   mm   × 2   mm × 25   μ m (length to width ratio R = 10 ) .

2.2. Precipitation Reaction Measurements

The measurements used in this work to make the numerical fittings correspond to the monitoring of the precipitation reaction of calcium oxalate crystals ( CaC 2 O 4 ) :
CaCl 2 ( a q ) + H 2 C 2 O 4 ( a q )   CaC 2 O 4 ( s ) + 2 HCl ( a q )
The formation of the salt crystals during the precipitation process was tracked by detecting changes in the sensor resonance frequency f m a x , as fully described by Sisniega et al. [20]. The magnetoelastic sensor described in the previous subsection was placed in a vial with a mixture of equal parts of the reagents ( CaCl 2 and H 2 C 2 O 4 ) at the same concentration (for three different concentrations: 30, 50, and 100 mM) and the changes in the resonance frequency of the sensor (determined by Equation (2)) were monitored along time, as the crystals precipitate and settle onto the sensor (increasing its mass). The experimental setup used to register the resonance curves consists of three coaxial solenoids: one to apply a constant bias field; a second one to produce the alternating field to magnetostrictively excite the sample; and the third one, consisting of a compensated pick-up coil, to detect the induced magnetization oscillations on the sensor. A spectrum analyzer (HP 3589A, Hewlett-Packard, Palo Alto, CA, USA) working in swept mode is used to produce the excitation and to receive the signal induced in the pick-up coil. After recording the resonance curves, the frequency of the maximum f m a x and its amplitude are measured using the built-in analysis procedures of the analyzer and transmitted to a control computer. The sweep speed was set looking for a compromise between good curve resolution and the necessary speed to be able to record the whole resonance-antiresonance curve fast enough to follow the precipitation process. With this configuration, the experimental uncertainty in the frequency determination is 100 Hz, and the measurement noise can be up to 1 kHz when the resonance curves are wide and with low amplitude.
Figure 2 presents an example of the results obtained in the precipitation monitoring measurements. The curves in Figure 2a show the decrease of the sensor resonance frequency ( f m a x ) as the precipitation takes place, for the three solutions with different concentrations. Figure 2b shows, for one particular solution (30 mM), the evolution of the complete resonance curve as the reaction progresses (curves taken at different reaction times). It can be seen that, in addition to the decrease in the resonance frequency, there is also a reduction of the quality factor Q that characterizes the resonance (or increase in damping), and a decrease in the amplitude of the sensor signal. The resonance curves displayed in Figure 2b were used in the present work to carry out the numerical fittings, together with analogous curves corresponding to the solutions of concentration 50 and 100 mM. Successive experiments in the same conditions yield similar results. Here we present the analysis in a single run for each concentration. Nine resonance curves corresponding to different reaction times were fitted for each concentration.

2.3. Numerial Fitting of the Sensor Resonance Curves

The quality factor Q of the resonance can be quantified from the width and the peak position of the resonance curve. However, when the resonance signals are heavily damped this and other simple methods are unable to determine this parameter or are highly inaccurate. We propose to use a frequency response fitting method to determine the quality factor and other parameters that characterize the resonance curves. In 1978, Savage et al. [22] derived an expression for the susceptibility of the magnetoelastic resonator around the magnetoelastic resonance given by
χ ( ω ) = χ 0 [ 1 8 k 2 π 2 n 1 n 2 × 1 1 ω n 2 ω 2 + j Q 1 ω n ω ]
where k is the magnetoelastic coupling coefficient ( k = π 2 8 ( 1 ( f r f a ) 2 ) , f r and f a being the resonance and anti-resonance frequencies, respectively), ω n = 2 π f n is the resonance frequency of the nth harmonic of the excited fundamental mode ( n = 1 ,   f 1 = f r ), Q 1 is a phenomenological damping coefficient, χ 0 is the magnetic susceptibility measured at a frequency far below the resonance and j = 1 . Due to the experimental procedure used to measure the resonance curves [20], the voltage induced in the pick-up coil (which is displayed in Figure 1 and Figure 2b) is proportional to the magnetic susceptibility. The parameters to be fitted in Equation (6) are: ω r , ω a , χ 0 , and Q . The numerical fitting of the experimental data to the modulus or magnitude of the Equation (6) for its first resonant mode ( n = 1 ) was performed using Mathematica© software. The parameters ω r , ω a , and χ 0 , were allowed to vary around their experimental corresponding values (extracted from the experimental data), until the value of Q that minimizes the norm between the fit and the experimental data was found. The norm (or residual) is defined as:
= 1 N i = 1 , N ( χ e x p , i χ f i t , i χ m a x ) 2
where χ m a x = max ( χ m a x , e x p , χ m a x , f i t ) , N is the number of experimental points ( N = 401 , in our measurements), and 0 1 , meaning values close to 0 that the fitting is good.
Additionally, a phenomenological expression using the formalism of linear systems was developed to describe the frequency response of a magnetoelastic material [17]. The magnetoelastic resonance is characterized by maximum amplitude at the resonance frequency f r , and a null minimum at the anti-resonance frequency f a . The analytical expression of the transference function ( G ) that describes a system with such frequency response is composed of a couple of complex conjugated poles in the denominator (describing the resonance) and a couple of complex conjugated zeros in the numerator (describing the anti-resonance):
G ( s ) = ω r 2 ω a 2 s 2 + 2 δ a ω a s + ω a 2 s 2 + 2 δ r ω r s + ω r 2
where ω r = 2 π f r is the resonance frequency, ω a = 2 π f a the anti-resonance frequency, and δ r and δ a are damping parameters. The transfer function is expressed as a Laplace operator with s = j ω , where j = 1 and ω is the frequency.
The amplitude of the frequency response of the system, when submitted to harmonic excitation, based on the frequency response described on the transference function is:
V ( ω ) = A ω r 2 ω a 2 | ω 2 2 j δ a ω a ω ω a 2 ω 2 2 j δ r ω r ω ω r 2 | + a ω + b
where A accounts for the amplitude and a and b for the background signal of the experiment (assumed to be linear). The parameters to be fitted in the expression are: ω r , ω a , δ r , δ a , A , a and b . Note that in this expression, the quality factor Q does not appear explicitly. Instead, the damping parameters δ r and δ a carry the relevant information. The numerical fitting of Equation (9) to the experimental data was carried out using MATLAB® software and a nonlinear least-squares fitting. The parameters ω r , ω a , and A to be fitted are initially taken as their values obtained from the experimental curve. All parameters are fitted until those that best fit the experimental curve are found, in this case those that correspond to a local minimum of a function that is a sum of squared residuals (being the residual the difference between the experimental value of the dependent variable and the value predicted by the fitting model).

3. Results and Discussion

3.1. Numerical Fitting of the Resonance Curves and Resonance Frequency

Numerical fittings were performed to the resonance curves measured at different times in the precipitation reaction, for solution with different concentrations (30, 50 and 100 mM). Figure 3 shows an example of the numerical fittings made with Equation (6), while Figure 4 shows the same fitting but using Equation (9). It can be seen that both models fit the experimental data reasonably well, resulting in a low value of the corresponding residuals.
However, when the deposited mass is high (at high reaction times or higher concentration of reagents), the fitting to Equation (9) behaves slightly better (with lower error between the experimental data and the values obtained through the model) than the one using the Equation (6).
If we analyze the evolution of the parameters obtained by the fitting procedure, we can see that the two damping parameters in Equation (9) ( δ r and δ a ) increase progressively as the mass deposited in the sensor increases (Figure 5). As the reaction advances, the crystals are formed in the solution of the vial and are deposited on the surface of the sensor, increasing its mass, and this increases the damping suffered by this sensor. In the same manner, as the concentration of the reagents increases (30, 50, 100 mM), the mass of calcium oxalate crystals is greater and, therefore, so is the damping. The numerical fittings show this direct relationship between the deposited mass and the damping in the resonance curves obtained. Similarly, it can be seen that the Q parameter obtained by fitting the data to Equation (6) decreases as the precipitation reaction progresses and the mass is deposited on the sensor (Figure 6).
Figure 6 shows the inverse relation between the quality factor Q and the damping parameters ( δ r and δ a ). All these parameters show a clear and marked tendency as the precipitation reaction occurs.
The evolution of resonance parameters, preferentially the resonance frequency, can be used to monitor the reaction by the changes caused by the mass deposition on the sensor signal. Numerical fittings improve the information obtained through the changes observed in the measured resonance frequency, since the result comes from the whole resonance curve, not only from the point of maximum amplitude and the width at half maximum. In addition, the fitting procedure is useful to reduce the noise, especially when the signal is low or highly damped, as evidenced in Figure 7a. This plot displays the frequency of the maxima (that is, the resonance frequency f m a x ) of all the resonance curves measured during the precipitation process of the solution with concentration 50 mM. The data obtained directly from the measurements present a considerable ripple with an amplitude close to 1 kHz, which severely limits the resolution of the sensor. It is mainly caused by the limited resolution with which frequency differences are discriminated by the measuring system. The evolution of the frequency of the maxima obtained from the fitting of the resonance curves displays a smooth monotonous tendency, which facilitates the calibration and use of the sensor.
One important consequence of the fitting procedure is that the values of the f r parameter obtained in the fittings do not correspond to f m a x , the position of the maxima of the resonance curves. Figure 7b compares the fitted values from both fitting expressions with the experimental resonance frequency ( f m a x , taken as the maximum of the resonance curve).
This tells us that what we usually take as the resonance frequency, that is the frequency of the maximum of the resonance curve, is not really the physical resonance frequency of the sensor, f r which can be obtained from the fitting to the analytical expressions. The frequency of the maximum of the measured curves, f m a x , is systematically lower than the resonance frequency, f r , since it is affected by the damping of the curve. The same applies to the frequency of anti-resonance (minimum of the curve), which, when affected by damping due to the mass deposition, increases with respect to the value, f a , which is obtained from the fits. Figure 8a illustrates this discrepancy, showing the frequency response of the system, for ideal resonance and anti-resonance separately, and for the full curve (sum of resonance and anti-resonance) affected by damping. In Figure 8b, the difference between the frequency of the maximum of the resonance f m a x and the resonance frequency f r is shown as a function of the damping, and a linear relation between them is observed.
The curves in Figure 8a were obtained with the magnitude of the frequency response corresponding to the transfer functions representing systems of: a pure resonance ( G 1 ), a pure anti-resonance ( G 2 ), and the combination of both ( G 3 ), which represents the observed resonance-antiresonance behaviour.
G 1 ( s ) = ω r 2 s 2 + 2 δ r ω r s + ω r 2
G 2 ( s ) = s 2 + 2 δ a ω a s + ω a 2 ω a 2
G 3 ( s ) = G 1 G 2 = ω r 2 ω a 2 s 2 + 2 δ a ω a s + ω a 2 s 2 + 2 δ r ω r s + ω r 2
The frequency values obtained for the resonance (maximum value) and anti-resonance (minimum value) of the respective curves in Figure 8a, coincide with those obtained in the fittings ( f r   F i t and f a   F i t in Figure 8a). In the combined curve, however, the maximum and minimum are shifted with respect to these values, resulting in the values of the resonance and anti-resonance frequencies ( f r   D a m p e d and f a   D a m p e d in Figure 8a), which match the experimental data.
The mass increase on the resonator produces both a decrease of the maximum of the resonance frequency ( f m a x ) and a decrease of the quality factor. Figure 9 illustrates the correlation between both parameters. However, the observed reduction of f m a x combines the effect of mass loading (Equation (2)) and the effect described in Figure 8a. The analysis of the resonance curves through the fitting procedure allows the determination of f r , which is, in principle, the one that is considered in Equation (2).

3.2. Comparison of the Results Obtained with Both Numerical Fittings

As it is to be expected, Equations (6) and (9) are analytically equivalent, when only the first resonant mode (n = 1) is considered in Equation (6). By rearranging the Equations (6) and (9) and after simple transformations, the following equivalences between the equation’s parameters are found:
δ r = 1 2 Q
δ a = ω a ω r 1 2 Q
A = χ 0
Figure 10 represents the difference between the fitted parameter δ r obtained with both numerical fittings, when it is related by means of the above equivalence expressions, for the data corresponding to a precipitation experiment. As can be seen, the difference between the damping parameter increases as the precipitation reaction advances. This is due to the fact that, as the mass of calcium oxalate is deposited in the resonator, the damping increases and the amplitude of the signal becomes smaller. This decrease of the signal causes that the background signal of the measurement gains in importance with respect to the resonance signal. It should be noted that Equation (9) explicitly takes this background into account, unlike Equation (6). For this reason, as the reaction progresses, the effect of the background signal is more noticeable in the fitted parameters. In order to elucidate the role that this background has in the values of the parameters obtained and to be able to compare both expressions under equal conditions (when neither of the two expressions takes this background into account), a fitting to a similar expression to Equation (9) where the background contribution is removed was carried out. The comparison of the damping parameters δ r obtained with this expression and with Equation (6) is shown in Figure 11. It can be seen that, indeed, in this case both models provide parameters of almost equal value.
Bearing this in mind, we can establish that the background plays an important role to make a proper fitting of the resonance curves of these sensors, especially when working under conditions where the damping is considerable and the signal amplitude is low.

3.3. Determination of the Quality Factor

Finally, a comparison of the quality factor obtained by the different procedures was carried out. Figure 12 shows the results of the quality factors obtained, on the one hand by the method of the half maximum bandwidth (Equation (4)), and on the other hand by the numerical fitting, calculating Q from Equation (13) in the case of the fitting to Equation (9). As can be seen, the value of the quality factor Q estimated from Equation (4) is not reliable for these measurements, where the curves are highly damped. In fact, this method is no longer applicable when the bandwidth, Δ f , cannot be calculated since the resonance signal becomes small and starts to be buried by the background signal, which occurs for times greater than t = 100 s in the case of reagents at a concentration of 100 mM, and for times above t = 300 s for concentrations of 50 mM (in continuous line of pale color in Figure 12). Table 2 shows a comparison of the Q values obtained through Equation (4) and through the fitting to Equations (6) and (9), for each concentration of reagents and different reaction times, and relative error between Q obtained through Equation (4) and Equation (9).
In addition, on behalf of a quick measurement, the frequency range in which the resonance curve is measured can be reduced. The fitting procedure can still provide, under these conditions, a valid value of the damping parameter, which the method of the half maximum bandwidth cannot give. As a proof of this, Figure 13 shows the fitting of the resonance curve in a reduced frequency range. The Q value obtained in this way, Q = 17.4, is practically the same as the one obtained with the fitting of the complete curve, in a wider frequency range, Q = 17.5.
Therefore, the quality factors obtained by means of numerical fittings make up for the lack of measurement frequency range and address the difficulty of estimating this parameter when the signals have low amplitude and are highly damped, providing a more accurate determination of Q (as already demonstrated by Lopes et al. [16]). This can be seen in Figure 12 and in data of Table 2, where Q values obtained by the different methods are more similar when the resonance curves are sharper, with less damping (at the beginning of the reaction), and differ as the reaction progresses and the signal quality decreases. In these cases, numerical fitting methods provide a better way to determine the quality factor Q . The result of both expressions (Equations (6) and (9)) used for the fitting is similar, and the slight difference between them arises from considering (or not) the background part (parameters a and b ) of the experimental measurements, as shown in Figure 11. Taking this into account, the fitting to Equation (9) is more precise and performs slightly better than the fitting to Equation (6).

4. Conclusions

In the present work, we have carried out the numerical fittings of the resonance curves of a magnetoelastic sensor designed to monitor the precipitation reaction of calcium oxalate. Through these fittings, different parameters describing the resonance curves, such as damping parameters, or the quality factor, have been obtained. It has been shown that these parameters are affected by the increase of mass of the sensor when the precipitation reaction advances, and therefore they give us information of how this mass deposition affects the sensor, increasing the damping suffered by it and decreasing its quality factor. These parameters have proven to have potential for monitoring the reaction, as they evolve as the reaction progresses.
The fitting procedure shows that the frequency at which the resonance curve is maximum ( f m a x ), does not coincide with the resonance frequency of the system, f r , as included in the analytical expressions used for the fits. The difference is produced because the damping of the resonance displaces the maximum of the resonance towards lower frequencies.
Clear analytical relationships have been established between the different parameters of the two fitting expressions used. Both expressions have proven to be suitable for making numerical fittings to the resonance curves and provide a very low residual in relation to the experimental data. However, there is a difference between them, and it is mainly due to the adjustment of the background signal, which when taken into account, provides a more precise fitting of the experimental data.
Resonance curve fits have proven to be a great alternative for determining the quality factor Q , and have shown good performance, especially when the resonance curves are highly damped or have low amplitude, where other determination methods have proven to be inaccurate.

Author Contributions

B.S., J.G. and A.G.-A. conceived and designed the work and needed measurements; B.S. performed the experiments; B.S. and V.M. performed all numerical fitting work; B.S. wrote the manuscript. All authors analyzed the data and discussed the results and implications, and commented on the manuscript at all stages. All authors have read and agreed to the published version of the manuscript.

Funding

B.S. thanks the Basque Government for funding under the FPI grant PRE_2019_1_0200. The authors would like to thank the financial support from the Basque Government under µ4indust project (KK-2019/00101, Elkartek program) and University Basque Research Groups Funding (IT1245-19).

Acknowledgments

Technical and human support provided by SGIker (UPV/EHU, MICINN, GV/EJ, ESF) is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Grimes, C.A.; Roy, S.C.; Rani, S.; Cai, Q. Theory, Instrumentation and Applications of Magnetoelastic Resonance Sensors: A Review. Sensors 2011, 11, 2809–2844. [Google Scholar] [CrossRef] [PubMed]
  2. Grimes, C.A.; Mungle, C.S.; Zeng, K.; Jain, M.K.; Dreschel, W.R.; Paulose, M.; Ong, K.G. Wireless Magnetoelastic Resonance Sensors: A Critical Review. Sensors 2002, 2, 294–313. [Google Scholar] [CrossRef] [Green Version]
  3. Kouzoudis, D.; Grimes, C.A. The frequency response of magnetoelastic sensors to stress and atmospheric pressure. Smart Mater. Struct. 2000, 9, 885. [Google Scholar] [CrossRef]
  4. Jain, M.K.; Schmidt, S.; Ong, K.G.; Mungle, C.; Grimes, C.A. Magnetoacoustic remote query temperature and humidity sensors. Smart Mater. Struct. 2000, 9, 502. [Google Scholar] [CrossRef]
  5. Grimes, C.A.; Kouzoudis, D.; Mungle, C. Simultaneous measurement of liquid density and viscosity using remote query magnetoelastic sensors. Rev. Sci. Instrum. 2000, 71, 3822–3824. [Google Scholar] [CrossRef]
  6. Bravo-Imaz, I.; García-Arribas, A.; Gorritxategi, E.; Arnaiz, A.; Barandiarán, J.M. Magnetoelastic Viscosity Sensor for On-Line Status Assessment of Lubricant Oils. IEEE Trans. Magn. 2013, 49, 113–116. [Google Scholar] [CrossRef]
  7. Cai, Q.Y.; Cammers-Goodwin, A.; Grimes, C.A. A wireless, remote query magnetoelastic CO2 sensor. J. Env. Monit. 2000, 2, 556–560. [Google Scholar] [CrossRef]
  8. Cai, Q.Y.; Grimes, C.A. A remote query magnetoelastic pH sensor. Sens. Actuators B Chem. 2000, 71, 112–117. [Google Scholar] [CrossRef]
  9. Puckett, L.G.; Barrett, G.; Kouzoudis, D.; Grimes, C.; Bachas, L.G. Monitoring blood coagulation with magnetoelastic sensors. Biosens. Bioelectron. 2003, 18, 675–681. [Google Scholar] [CrossRef]
  10. Lakshmanan, R.S.; Guntupalli, R.; Hu, J.; Kim, D.J.; Petrenko, V.A.; Barbaree, J.M.; Chin, B.A. Phage immobilized magnetoelastic sensor for the detection of Salmonella typhimurium. J. Microbiol. Methods 2007, 71, 55–60. [Google Scholar] [CrossRef]
  11. Ruan, C.; Zeng, K.; Varghese, O.K.; Grimes, C.A. Magnetoelastic Immunosensors: Amplified Mass Immunosorbent Assay for Detection of Escherichia coli O157: H7. Anal. Chem. 2003, 75, 6494–6498. [Google Scholar] [CrossRef] [PubMed]
  12. Zhang, K.; Zhang, L.; Fu, L.; Li, S.; Chen, H.; Cheng, Z.Y. Magnetostrictive resonators as sensors and actuators. Sens. Actuator A Phys. 2013, 200, 2–10. [Google Scholar] [CrossRef]
  13. Petersan, P.J.; Anlage, S.M. Measurement of resonant frequency and quality factor of microwave resonators: Comparison of methods. J. Appl. Phys. 1998, 84, 3392–3402. [Google Scholar] [CrossRef] [Green Version]
  14. Kaczkowski, Z. Piezomagnetic parameters of the magnetostrictive materials. Arch. Acoust. 2014, 23, 307–329. [Google Scholar]
  15. Miller, J.M.L.; Ansari, A.; Heinz, D.B.; Chen, Y.; Flader, I.B.; Shin, D.D.; Kenny, T.W. Effective quality factor tuning mechanisms in micromechanical resonators. Appl. Phys. Rev. 2018, 5, 041307. [Google Scholar] [CrossRef]
  16. Lopes, A.C.; Sagasti, A.; Lasheras, A.; Muto, V.; Gutiérrez, J.; Kouzoudis, D.; Barandiarán, J.M. Accurate determination of the Q quality factor in magnetoelastic resonant platforms for advanced biological detection. Sensors 2018, 18, 887. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. García-Arribas, A.; Gutiérrez, J.; Kurlyandskaya, G.V.; Barandiarán, J.M.; Svalov, A.; Fernández, E.; Lasheras, A.; de Cos, D.; Bravo-Imaz, I. Sensor applications of soft magnetic materials based on magneto-impedance, magneto-elastic resonance and magneto-electricity. Sensors 2014, 14, 7602–7624. [Google Scholar] [CrossRef] [Green Version]
  18. Bravo, I.; Arnaiz, A.; García-Arribas, A. Damping of Magnetoelastic Resonance for Oil Viscosity Sensing. IEEE Trans. Magn. 2019, 55, 1–5. [Google Scholar] [CrossRef]
  19. Bouropoulos, N.; Kouzoudis, D.; Grimes, C. The real-time, in situ monitoring of calcium oxalate and brushite precipitation using magnetoelastic sensors. Sens. Actuators B Chem. 2005, 109, 227–232. [Google Scholar] [CrossRef]
  20. Sisniega, B.; Sagasti Sedano, A.; Gutiérrez, J.; García-Arribas, A. Real Time Monitoring of Calcium Oxalate Precipitation Reaction by Using Corrosion Resistant Magnetoelastic Resonance Sensors. Sensors 2020, 20, 2802. [Google Scholar] [CrossRef]
  21. Sagasti, A.; Palomares, V.; Porro, J.M.; Orúe, I.; Sánchez-Ilárduya, M.B.; Lopes, A.C.; Gutiérrez, J. Magnetic, Magnetoelastic and Corrosion Resistant Properties of (Fe–Ni)-Based Metallic Glasses for Structural Health Monitoring Applications. Materials 2020, 13, 57. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Savage, H.; Abbundi, R. Perpendicular susceptibility, magnetomechanical coupling and shear modulus in Tb.27Dy.73Fe2. IEEE Trans. Magn. 1978, 14, 545–547. [Google Scholar] [CrossRef]
Figure 1. Two magnetoelastic resonance curves from the same resonator with different Q values, caused by different amount of mass loading. The peak corresponding to the higher Q (blue) is sharper and narrower, the resonance curve with a lower Q is wider (red).
Figure 1. Two magnetoelastic resonance curves from the same resonator with different Q values, caused by different amount of mass loading. The peak corresponding to the higher Q (blue) is sharper and narrower, the resonance curve with a lower Q is wider (red).
Materials 13 04708 g001
Figure 2. (a) Temporal evolution of the magnetoelastic resonance frequency of the sensor during the reaction of precipitation for different concentrations (30 mM, 50 mM and 100 mM) and a control curve (sensor in a vial with distilled water). The value of the resonance frequency, f m a x , is obtained as the frequency at the maximum amplitude of measured resonance curves (Figure 2b); (b) Measured magnetoelastic resonance curve signals of the sensor at different times during the precipitation process for the concentration 30 mM. Data taken from Ref. [20].
Figure 2. (a) Temporal evolution of the magnetoelastic resonance frequency of the sensor during the reaction of precipitation for different concentrations (30 mM, 50 mM and 100 mM) and a control curve (sensor in a vial with distilled water). The value of the resonance frequency, f m a x , is obtained as the frequency at the maximum amplitude of measured resonance curves (Figure 2b); (b) Measured magnetoelastic resonance curve signals of the sensor at different times during the precipitation process for the concentration 30 mM. Data taken from Ref. [20].
Materials 13 04708 g002
Figure 3. Fit to Equation (6). (a) Measured magnetoelastic resonance curves of the sensor at different times during the precipitation process for reactants of concentration 30 mM, and corresponding numerical fittings (in dashed lines) to Equation (6); (b) Measured magnetoelastic resonance curve, fitting to Equation (6) (in dashed line), and corresponding error (in red) of the sensor at time t = 125 s for reactants of concentration 30 mM.
Figure 3. Fit to Equation (6). (a) Measured magnetoelastic resonance curves of the sensor at different times during the precipitation process for reactants of concentration 30 mM, and corresponding numerical fittings (in dashed lines) to Equation (6); (b) Measured magnetoelastic resonance curve, fitting to Equation (6) (in dashed line), and corresponding error (in red) of the sensor at time t = 125 s for reactants of concentration 30 mM.
Materials 13 04708 g003
Figure 4. Fit to Equation (9). (a) Measured magnetoelastic resonance curves of the sensor at different times during the precipitation process for reactants of concentration 30 mM, and corresponding numerical fittings (in dashed lines) to Equation (9); (b) Measured magnetoelastic resonance curve, fitting to Equation (9) (in dashed line), and corresponding error (in red) of the sensor at time t = 125 s for reactants of concentration 30 mM.
Figure 4. Fit to Equation (9). (a) Measured magnetoelastic resonance curves of the sensor at different times during the precipitation process for reactants of concentration 30 mM, and corresponding numerical fittings (in dashed lines) to Equation (9); (b) Measured magnetoelastic resonance curve, fitting to Equation (9) (in dashed line), and corresponding error (in red) of the sensor at time t = 125 s for reactants of concentration 30 mM.
Materials 13 04708 g004
Figure 5. Evolution of the damping parameters obtained through the numerical fitting to Equation (9) of the resonance curves during the precipitation time for different concentration of reagents (30, 50, and 100 mM). (a) Damping parameter δ r ; (b) Damping parameter δ a .
Figure 5. Evolution of the damping parameters obtained through the numerical fitting to Equation (9) of the resonance curves during the precipitation time for different concentration of reagents (30, 50, and 100 mM). (a) Damping parameter δ r ; (b) Damping parameter δ a .
Materials 13 04708 g005
Figure 6. Evolution of the quality factor Q of the resonator during the precipitation process for different reagents concentrations (30, 50, and 100 mM) obtained through the numerical fittings to Equation (6) of the resonance curves.
Figure 6. Evolution of the quality factor Q of the resonator during the precipitation process for different reagents concentrations (30, 50, and 100 mM) obtained through the numerical fittings to Equation (6) of the resonance curves.
Materials 13 04708 g006
Figure 7. (a) Evolution of the frequency corresponding to the maximum of the resonance curves resulting from the fitting to Equation (9) of the resonance curves measured during the precipitation, when the concentration of the reactants is 50 mM, compared to experimental data ( f m a x ); (b) Evolution of the f r obtained in the numerical fitting to Equations (6) and (9) of the measured resonance curves, compared to experimental data ( f m a x ).
Figure 7. (a) Evolution of the frequency corresponding to the maximum of the resonance curves resulting from the fitting to Equation (9) of the resonance curves measured during the precipitation, when the concentration of the reactants is 50 mM, compared to experimental data ( f m a x ); (b) Evolution of the f r obtained in the numerical fitting to Equations (6) and (9) of the measured resonance curves, compared to experimental data ( f m a x ).
Materials 13 04708 g007
Figure 8. (a) Frequency response of the system for resonance and anti-resonance behavior, and the total response affected by the damping. Corresponding resonance and anti-resonance frequencies are indicated; (b) Evolution of the difference between the experimental resonance frequency f m a x and the resonance frequency obtained through the fitting to Equation (9) as the damping parameter δ r increases, for the case of concentration 50 mM. The shift in the resonance frequency becomes greater as the damping increases.
Figure 8. (a) Frequency response of the system for resonance and anti-resonance behavior, and the total response affected by the damping. Corresponding resonance and anti-resonance frequencies are indicated; (b) Evolution of the difference between the experimental resonance frequency f m a x and the resonance frequency obtained through the fitting to Equation (9) as the damping parameter δ r increases, for the case of concentration 50 mM. The shift in the resonance frequency becomes greater as the damping increases.
Materials 13 04708 g008
Figure 9. Relation between the resonance frequency f m a x obtained experimentally and the quality factor Q obtained through the fitting to Equation (9) at different times during the precipitation reaction. The resonance curves corresponding to times t = 6 s (high Q ) and t = 258 s (low Q ) are shown.
Figure 9. Relation between the resonance frequency f m a x obtained experimentally and the quality factor Q obtained through the fitting to Equation (9) at different times during the precipitation reaction. The resonance curves corresponding to times t = 6 s (high Q ) and t = 258 s (low Q ) are shown.
Materials 13 04708 g009
Figure 10. Comparison of the damping parameter δ r obtained with both numerical fittings of the resonance curves (for reactants of concentration 50 mM) using the equivalences found between δ r and Q , Equation (13).
Figure 10. Comparison of the damping parameter δ r obtained with both numerical fittings of the resonance curves (for reactants of concentration 50 mM) using the equivalences found between δ r and Q , Equation (13).
Materials 13 04708 g010
Figure 11. Comparison of the damping parameter δ r obtained by the fitting to Equation (9) without taking into account the background parameters ( a and b ), and the damping parameter δ r obtained through the fitting to Equation (6).
Figure 11. Comparison of the damping parameter δ r obtained by the fitting to Equation (9) without taking into account the background parameters ( a and b ), and the damping parameter δ r obtained through the fitting to Equation (6).
Materials 13 04708 g011
Figure 12. Comparison of the quality factor Q obtained by the fitting to Equations (6) and (9), and the one obtained through Equation (4).
Figure 12. Comparison of the quality factor Q obtained by the fitting to Equations (6) and (9), and the one obtained through Equation (4).
Materials 13 04708 g012
Figure 13. Magnetoelastic resonance curve measured at time t = 203 s during the precipitation reaction with 50 mM reactants (in dashed red line), and corresponding fitting to Equation (9) (in dashed black line) for a reduced range of frequencies (in blue). The value of the quality factor obtained with the reduced fitting is Q = 17.4.
Figure 13. Magnetoelastic resonance curve measured at time t = 203 s during the precipitation reaction with 50 mM reactants (in dashed red line), and corresponding fitting to Equation (9) (in dashed black line) for a reduced range of frequencies (in blue). The value of the quality factor obtained with the reduced fitting is Q = 17.4.
Materials 13 04708 g013
Table 1. Magnetic, magnetoelastic, and corrosion behaviour parameters of the Fe 73 Cr 5 Si 10 B 12 sample [21]. M s is the spontaneous magnetization, λ s the saturation magnetostriction, Δ E the change in the Young’s modulus with the applied magnetic field (or Δ E effect), k the magnetoelastic coupling coefficient and E c o r r the corrosion potential.
Table 1. Magnetic, magnetoelastic, and corrosion behaviour parameters of the Fe 73 Cr 5 Si 10 B 12 sample [21]. M s is the spontaneous magnetization, λ s the saturation magnetostriction, Δ E the change in the Young’s modulus with the applied magnetic field (or Δ E effect), k the magnetoelastic coupling coefficient and E c o r r the corrosion potential.
Composition μ 0 M s   ( T ) λ s (ppm) Δ E   ( % ) k E c o r r
(mV)
Corrosion Rate
(μm/Year)
Fe73Cr5Si10B121.1214170.41470.035
Table 2. Comparison of the Q values obtained through Equation (4) and the fitting to Equations (6) and (9) for each concentration of reactants and different times of reaction; and relative error between Q obtained through Equation (4) and Equation (9).
Table 2. Comparison of the Q values obtained through Equation (4) and the fitting to Equations (6) and (9) for each concentration of reactants and different times of reaction; and relative error between Q obtained through Equation (4) and Equation (9).
Concentration (mM)Time (s)Q (Equation (4))Q (Equation (6))Q (Equation (9))Relative Error * (%)
3010.247.850.452.48.7
4345.847.950.28.7
7040.543.545.611.2
12531.535.237.014.8
31716.822.524.130.3
42612.619.320.939.7
506.345.647.750.08.8
2838.942.144.011.6
7224.629.831.321.4
12115.521.823.433.7
312-12.915.0-
477-12.114.2-
100638.941.744.212
3914.621.122.936.2
993.611.213.473.1
154-8.911.3-
368-4.97.7-
478-4.16.9-
* Relative error between Q factors obtained through Equation (4) and through the fitting to Equation (9) calculated as: 100 × ( Q E q .   ( 9 ) Q E q . ( 4 ) ) / Q E q . ( 9 ) . The cells marked with a dash correspond to the cases where Equation (4) is not applicable.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sisniega, B.; Gutiérrez, J.; Muto, V.; García-Arribas, A. Improved Determination of Q Quality Factor and Resonance Frequency in Sensors Based on the Magnetoelastic Resonance Through the Fitting to Analytical Expressions. Materials 2020, 13, 4708. https://doi.org/10.3390/ma13214708

AMA Style

Sisniega B, Gutiérrez J, Muto V, García-Arribas A. Improved Determination of Q Quality Factor and Resonance Frequency in Sensors Based on the Magnetoelastic Resonance Through the Fitting to Analytical Expressions. Materials. 2020; 13(21):4708. https://doi.org/10.3390/ma13214708

Chicago/Turabian Style

Sisniega, Beatriz, Jon Gutiérrez, Virginia Muto, and Alfredo García-Arribas. 2020. "Improved Determination of Q Quality Factor and Resonance Frequency in Sensors Based on the Magnetoelastic Resonance Through the Fitting to Analytical Expressions" Materials 13, no. 21: 4708. https://doi.org/10.3390/ma13214708

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop