Next Article in Journal
Fundamental Research on the Invention of Breast Support Underwear with a Focus on Women’s Comfort
Next Article in Special Issue
Automatic Classification of Rotor Faults in Soft-Started Induction Motors, Based on Persistence Spectrum and Convolutional Neural Network Applied to Stray-Flux Signals
Previous Article in Journal
Interrogation Method with Temperature Compensation Using Ultra-Short Fiber Bragg Gratings in Silica and Polymer Optical Fibers as Edge Filters
Previous Article in Special Issue
Experimental Investigation and CFD Analysis of Pressure Drop in an ORC Boiler for a WHRS Implementation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Compensation of Thermal Gradients Effects on a Quartz Crystal Microbalance

by
Marianna Magni
1,2,
Diego Scaccabarozzi
1 and
Bortolino Saggin
1,*
1
Department of Mechanical Engineering Politecnico di Milano, Polo Territoriale di Lecco, Via G. Previati 1/c, 23900 Lecco, Italy
2
Rebel Dynamics, Via Carlo Porta 38, Cesana Brianza, 23861 Lecco, Italy
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(1), 24; https://doi.org/10.3390/s23010024
Submission received: 8 November 2022 / Revised: 15 December 2022 / Accepted: 16 December 2022 / Published: 20 December 2022
(This article belongs to the Special Issue Feature Papers in Physical Sensors 2022)

Abstract

:
Quartz Crystal Microbalances (QCM) are widely used instruments thanks to their stability, low mass, and low cost. Nevertheless, the sensitivity to temperature is their main drawback and is often a driver for their design. Though the crystal average temperature is mostly considered as the only disturbance, temperature affects the QCM measurements also through the in-plane temperature gradients, an effect identified in the past but mostly neglected. Recently, it has been shown that this effect can prevail over that of the average temperature in implementations where the heat for thermal control is released directly on the crystal through deposited film heaters. In this study, the effect of temperature gradients for this kind of crystal is analyzed, the sensitivity of frequency to the average temperature gradient on the electrode border is determined, and a correction is proposed and verified. A numerical thermal model of the QCM has been created to determine the temperature gradients on the electrode borders. The frequency versus temperature-gradient function has been experimentally determined in different thermal conditions. The correction function has been eventually applied to a QCM implementing a crystal of the same manufacturing lot as the one used for the characterization. The residual errors after the implementation of the correction of both average temperature and temperature gradients were always lower than 5% of the initial temperature disturbance. Moreover, using the correlation between the heater power dissipation and the generated temperature gradients, it has been shown that an effective correction strategy can be based on the measurement of the power delivered to the crystal without the determination of the temperature gradient.

1. Introduction

Quartz Crystal Microbalances are continuously expanding their application field because of the high and tunable sensitivity to mass deposition, stability in time, low cost, low power consumption, ease of reading, and capability of operating in very harsh environments, including high vacuum, cryogenic temperatures, and zero gravity [1,2]. Although QCMs have been used in many ways for a long time, the interest in this kind of sensor is still high because many different measurements can be done, either designing specific coatings selectively trapping specific chemical species [3,4,5,6] or using them to identify condensation or sublimation temperatures as in dew point sensors or for Thermo Gravimetric Analysis [7,8]. In most QCM applications, the crystal temperature is a relevant parameter, and, in many, it must be actively controlled to achieve the specific measurement goals. Temperature, nevertheless, is also the major disturbance in QCM measurements. Minimization of the temperature effects has been sought since the beginning of quartz crystals oscillators usage. The identification of the optimal crystal cut angle has been addressed [9,10], and configurations allowing for better temperature compensation are continuously investigated [11,12]. The average crystal temperature affects the crystal resonance; also, the temperature gradients within the crystal have similar effects [13,14,15]. Compensation of temperature effects based on the dual crystal configuration exploits the beating frequency of a reference crystal and the active one. The system is conceived to remove the effect of the average crystal temperature, assuming that both crystals have the same temperature [11,12], but the effect of temperature gradients is neglected. The uniform temperature condition is hardly achieved when the crystal is actively thermally controlled [16], and the gradient effect may become dominant when the heat is delivered directly on the crystal through film heaters deposited on it [17]. This work, starting from the findings reported in [17], derives the correlation between temperature gradient and frequency change, providing a relationship that can be used to compensate for the effect. An experimental campaign was performed on a QCM, measuring the resonance frequency in different thermal conditions. For each test condition, the temperature gradient was evaluated through a validated thermal model. The effect of the average crystal temperature was preliminarily removed, leaving the effect of the temperature gradient alone.

2. Temperature Gradient Determination from the Thermal Model

The thermal model of the QCM was developed to compute the temperature field on the crystal with the required accuracy and spatial resolution. The direct measurement of temperature through thermal mapping, in fact, was not accurate enough to derive the temperature gradient [17] because of the background disturbance, the low emissivity of the deposited electrodes, and the coarse spatial resolution of the infrared camera. As a consequence, the thermal model was implemented and preliminarily validated by correlation with the experimental tests described in [17]. The validated model has afterward been used to analyze crystals under different thermal conditions. The model geometry is shown in Figure 1.
A FE model based on the commercial software PTC CREO® 9 was developed, representing the crystal and the mechanical mounting holding it. A relevant parameter for the modeling is the contact thermal resistance between the crystal and the supports. The contact was modeled as bonded surfaces for the contact area determined with the Hertzian contact model. The mounting thermal resistance, nevertheless, was one of the parameters for which experimental validation was mandatory. The materials used in the modeling were quartz for the crystal, AISI316 stainless steel for the supports, and a ceramic material, Macor®, for the support’s spacers. All materials were modeled as linear and isotropic. This is a simplification for the quartz crystal [18] that was adopted because, for the purposes of the analysis, the properties in the crystal plane are dominant. Table 1 reports the material properties used in the model; the AISI316 was in the CREO database, while the MACOR parameters were derived from the manufacturer datasheet and the quartz crystal ones from [14].

2.1. Analysis Conditions

To validate the model, the same conditions of the experimental tests (heating power and environmental temperature) were simulated. The reference analysis considered a power dissipation of the device of 0.445 W localized on the heater area at the crystal bottom surface. The total heat load was distributed considering the power density due to the film resistance shape. The power is non-uniformly distributed throughout the heater path, shaped with different film cross-sections to improve the temperature uniformity of the crystal. There are different power densities in three regions, as shown in Figure 2, according to the resistor shape.

2.2. Boundary Conditions

The thermal analysis is considered as boundary conditions for the convective and radiative exchange factors with the environment and the temperature of the spacers. The convection coefficient was computed considering the natural convection on a horizontal plate for the crystal [19]:
h c o n v = C 1 ( Δ T M ) 1 / 4
where ΔT is the difference between surface and environment temperatures, M is the area-to-perimeter ratio; C1 is a constant that for turbulent flow is 1.71 W m−1.75 K−1.25. The above equation was used for the crystal top surface, while the bottom corresponding value for downward-facing plates was exploited. The radiative exchanges make the model non-linear. However, both radiation and convection constraints were applied to the same surfaces with a linearization process. Starting from the equation:
W = σ ε A ( T w 4 T 4 )
where Tw is the temperature of the surface and T the air one. After linearization it becomes:
W = h r A ( T w T )
where hr is the linearized radiation coefficient, given by:
h r = σ ε ( T w 2 + T 2 ) ( T w + T )
Eventually, the global heat exchange coefficient was obtained by adding the convection and the radiation coefficients:
h = h c o n v + h r
The last constraint imposed on the model was the temperatures of the supports’ spacer derived from the thermal images. In Figure 3, the values for the three supports are reported.

2.3. Results

Figure 4 shows both the temperature and the temperature gradient maps on the crystal determined with the thermal model.
The plot in Figure 5 shows the temperature profile along the electrode edge. It can be noticed that the lower temperature value is reached on support 1. This is consistent with the observation that this area, the electrode pad, is the farthest from the heater. Considering the plot of Figure 6, the maximum gradient is between the supports 1–3 and 1–2, where there is the largest heat density. The minimum gradient on the electrode border is where there is no heat dissipation because the heater is interrupted because of the support 1 electrode pad. Furthermore, the thermal field is almost symmetric with respect to an axis passing through the center of the disk surface and contact 1, i.e., rotated 60° from the x-axis, as shown in Figure 7. It can be noticed that the hottest areas are on the heater but not in the regions at the highest power density, where the cooling of support 1 is prevailing. Concerning the crystal, the lowest temperature is about 68 °C, reached on the regions of contact with the supports, cooled by conduction.

2.4. Model Validation

The thermal field predicted by the thermal analysis was compared with the one measured with a thermal mapper [17]. Five points of the electrode were identified on the model and compared with those corresponding to the thermal images. The same comparison was made for the heater region considering four points. Figure 8 and Figure 9 show the temperature maps, whereas Table 2 summarizes the prediction errors.
The temperatures measured during the test and those predicted by the thermal model are quite close: the maximum error on the electrode temperature is 1 °C, while the one on the heater is 0.34 °C. Considering that the uncertainty in the electrode temperature measurement was 0.85 °C and on the heater 0.42 °C [17], the contribution of the model uncertainty is of the order of 0.5 °C, a value considered acceptable for this study. The above comparison, therefore, validated the numerical model.

3. Temperature Gradient Sensitivity

To correlate the temperature gradient with the frequency variation, different thermal loads and boundary conditions used during the experimental activity were simulated, assessing the temperature gradients as a result. The different gradients were obtained by feeding the heater with different powers and keeping constant the other conditions. Generating different temperature gradients while maintaining the same average temperature would have been quite complex, requiring the contemporary application of heating and cooling on different crystal regions. Thus, the resulting set of conditions was characterized by different temperature gradients and different average temperatures. Assuming that the effect of the temperature gradients and that of the average temperature can be superimposed, the dependency of the frequency on the temperature gradient was determined by preliminarily removing the effect of the average temperature. The correction for the average temperature was based on the crystal characterization reported in [17], where the behavior of the crystal in a uniform temperature environment was assessed and discussed.
Low-temperature tests, below 50 °C, were eventually discarded because of the interfering effect of contaminants condensation on the crystal.
The thermal gradients were evaluated for different power dissipations of the heaters and used to assess the effects on the crystal frequency are shown in Figure 10.
The average thermal gradient was correlated with the frequency change after the application of the correction for the average temperature. The results are reported in Figure 11.
The correlation coefficient for the thermal gradient versus frequency fitting with a 3rd-order polynomial was 0.964 with a root mean square error of 55 Hz. Eventually, the effectiveness of the correction procedure was verified with a new experimental campaign on a standard QCM, including a crystal of the same manufacturing lot as the one used for characterization. By using the model of the thermal gradient effects, i.e., the polynomial regression of the gradient values of Figure 11, along with that of the uniform temperature, the frequency variation of the crystal with more general temperature fields can be predicted. To validate the model, the resonance frequency was measured while supplying different powers to the crystal resistor. Thermal analyses were performed to compute the corresponding gradient, and the corresponding corrective factor was evaluated. The average temperature on the crystal border was measured with the RTD deposited on the surface of the crystal and compared with the thermal model predictions, verifying that the differences were within the already identified range, i.e., within ±0.5 °C.
The average temperature was then used to compute the corresponding frequency change under uniform temperature. Eventually, by adding the frequency variation due to the thermal gradient with the previous one, the predicted frequency variation was obtained. The predicted frequency, along with the measured one, and the effect of the uniform temperature are reported in Figure 12.
It can be noticed that at low temperatures (between 40 °C and 50 °C), the correction significantly overestimates the measured frequency change. This has been attributed to the condensation of contaminants that was also observed in other low-temperature tests [17], so these points were not considered in the assessment of the model error. With the crystal temperature increasing, the frequency change becomes more relevant, but it is still quite accurately predicted by the model, allowing it to correct for the effect of the gradient with residual errors smaller than 5%.
The above procedure relies on the thermal model to determine the temperature gradient on the crystal. Nevertheless, considering that the heat flux is proportional to the temperature gradient, a correlation between the average gradient on the electrode border and the power dissipated by the heater deposited right there is expected. The mean thermal gradients, computed while performing five heating cycles, were plotted against the heating power in Figure 13. The correlation is evident from the linear plot. Moreover, the least square fitting leads to a linearity standard deviation of 0.027 [°C/mm]. The latter figure would be the standard uncertainty of the temperature gradient when it was determined from the power dissipation. The propagation of this contribution on the correction function leads to a negligible contribution in comparison with the observed 5% error. Therefore, the correction procedure based directly on the measurement of the heating power performs with the same accuracy as the one shown above, based on the evaluation of the temperature gradient but without requiring the implementation of the crystal thermal analysis.

4. Conclusions

The full characterization of the temperature effects on QCMs must investigate along with the effect of the average temperature, the one related to the temperature gradients. In some crystal configurations, like those recently proposed with built-in film heaters and temperature sensors, the temperature gradients may have a dominant effect. In this study, the relationship between the average temperature gradient on the electrode border and the frequency change was experimentally determined for QCM crystals of this type. The temperature gradient was determined for different test conditions through a specifically developed and validated thermal model. The modeling of both effects of temperature, i.e., the average and the spatial gradients, was performed for temperatures of the crystal between the ambient one and 100 °C and with average gradient amplitudes from 2 to 8 °C/mm.
The correction procedure of the temperature effects based on this model allowed for achieving residual errors below 5% for crystals temperatures above 50 °C. More significant discrepancies were observed at lower temperatures of the crystal. Nevertheless, in these conditions, the condensation of contaminants on the QCM jeopardizes the testing repeatability. Testing under a high vacuum is expected to improve the repeatability also in this condition; nevertheless, it is less interesting because the effect is not very large. Although the procedure has been developed using the temperature gradient determined from the thermal model of the crystal, it was additionally shown that evaluating the temperature gradient through the correlation with the power dissipation of the heaters allows for performing the correction without significant accuracy loss. This implementation strongly reduces the complexity and the effort required for the correction because all the needed information is directly available from the standard QCM monitoring parameters. This study allows for an easy correction of the temperature gradient effect when the environmental temperature is close to 20 °C; QCM is often operating at cryogenic temperatures, and extending the characterization in these conditions is a natural extension. Considering that reducing the disturbance source would ease any correction procedure, the other evolution of this research will be in the direction of changing the crystal design to reduce the temperature disturbance. This can be achieved by optimizing the power distribution on the crystal, with a proper shaping of the deposited film heater, for the purpose of minimizing the temperature gradients.

Author Contributions

Methodology, B.S.; Validation, D.S.; Investigation, M.M. and D.S.; Data curation, M.M.; Writing—original draft, M.M.; Writing—review & editing, D.S. and B.S.; Supervision, B.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wallace, D.A.; Wallace, S.A. Realistic performance specifications for flight quartz crystal microbalance instruments for contamination measurement on spacecraft. In Proceedings of the AIAA 23rd Thermophysics, Plasmadynamics and Lasers Conference, San Antonio, CA, USA, 27–29 June 1988. [Google Scholar]
  2. Wallace, D.A. Performance of a flight qualified, thermoelectrically temperature controlled QCM sensor with power supply, thermal controller and signal processor. In Proceedings of the NASA Conference Publication, Hampton, VI, USA, 5–7 November 1980; pp. 257–268. [Google Scholar]
  3. Yan, W.; Zhang, D.; Liu, X.; Chen, X.; Yang, C.; Kang, Z. Guar Gum/Ethyl Cellulose-Polyvinyl Pyrrolidone Composite-Based Quartz Crystal Microbalance Humidity Sensor for Human Respiration Monitoring. ACS Appl. Mater. Interfaces 2022, 14, 31343–31353. [Google Scholar] [CrossRef] [PubMed]
  4. Yamaguchi, T.; Matsuoka, T. Measurement of complex shear viscosity up to 3 GHz using an electrodeless AT-cut quartz transducer. Jpn. J. Appl. Phys. 2022, 61, SG1021. [Google Scholar] [CrossRef]
  5. Kumar, A.; Prajesh, R. The potential of acoustic wave devices for gas sensing applications. Sens. Actuators A Phys. 2022, 339, 113498. [Google Scholar] [CrossRef]
  6. Zhang, D.; Kang, Z.; Liu, X.; Guo, J.; Yang, Y. Highly sensitive ammonia sensor based on PSS doped ZIF-8-derived porous carbon/polyaniline hybrid film coated on quartz crystal microbalance. Sens. Actuators B Chem. 2022, 357, 131419. [Google Scholar] [CrossRef]
  7. Johnson, W.L.; Mansfield, E. Thermogravimetric analysis with a heated quartz crystal microbalance. In Proceedings of the 2012 IEEE International Frequency Control Symposium Proceedings (IFCS 2012), Baltimore, MD, USA, 21–24 May 2012; pp. 609–613. [Google Scholar]
  8. Dirri, F.; Palomba, E.; Longobardo, A.; Biondi, D.; Boccaccini, A.; Galiano, A.; Zampetti, E.; Saggin, B.; Scaccabarozzi, D.; Martin-Torres, J. VISTA Instrument: A PCM-Based sensor for organics and volatiles characterization by using thermogravimetric technique. In Proceedings of the 5th IEEE International Workshop on Metrology for AeroSpace, MetroAeroSpace 2018-Proceedings, Rome, Italy, 20–22 June 2018; pp. 150–154. [Google Scholar]
  9. Mason, W.P. Low Temperature Coefficient Quartz Crystals. Bell Labs Tech. J. 1940, 19, 74–93. [Google Scholar] [CrossRef]
  10. Mason, W.P.; Sykes, R.A. Low-Frequency Quartz-Crystal Cuts Having Low Temperature Coefficients. Proc. IRE 1944, 32, 208–215. [Google Scholar] [CrossRef]
  11. Tsuchiya, Y.; Kukita, H.; Shiobara, T.; Yukumatsu, K.; Miyazaki, E. Temperature Controllable QCM Sensor with Accurate Temperature Measurement for Outgas and Contamination Assessment. In Proceedings of the IEEE Sensors, Montreal, QC, Canada, 27–30 October 2019. [Google Scholar]
  12. Scaccabarozzi, D.; Saggin, B.; Tarabini, M.; Palomba, E.; Longobardo, A.; Zampetti, E. Thermo-mechanical design and testing of a microbalance for space applications. Adv. Space Res. 2014, 54, 2386–2397. [Google Scholar] [CrossRef]
  13. Valentin, J.P.; Theobald, G.; Gagnepain, J.J. Frequency Shifts Arising from in-Plane Temperature Gradient Distribution in Quartz Resonators. In Proceedings of the Annual IEEE International Frequency Control Symposium, Philadelphia, PA, USA, 29 May–1 June 1984; pp. 157–163. [Google Scholar]
  14. Valentin, J.P.; Théobald, G.; Gagnepain, J.J. Temperature-induced frequency shifts in quartz resonators. J. Appl. Phys. 1985, 58, 1388–1396. [Google Scholar] [CrossRef]
  15. Sinha, B.K. On the Interaction of Temperature Gradients with Surface Wave Propagation in Anisotropic Solids. IEEE Trans. Sonics Ultrason. 1985, 32, 583–591. [Google Scholar] [CrossRef]
  16. Li, N.; Meng, X.; Nie, J.; Lin, L. A QCM Dew Point Sensor With Active Temperature Control Using Thermally Conductive Electrodes. IEEE Sens. J. 2018, 18, 5715–5722. [Google Scholar] [CrossRef]
  17. Magni, M.; Scaccabarozzi, D.; Palomba, E.; Zampetti, E.; Saggin, B. Characterization of Thermal Gradient Effects on a Quartz Crystal Microbalance. Sensors 2022, 22, 7256. [Google Scholar] [CrossRef] [PubMed]
  18. Gibert, B.; Mainprice, D. Effect of crystal preferred orientations on the thermal diffusivity of quartz polycrystalline aggregates at high temperature. Tectonophysics 2009, 465, 150–163. [Google Scholar] [CrossRef]
  19. Rohsenow, W.R.; Hartnett, J.P.; Cho, Y.I. Handbook of Heat Transfer; McGraw Hill: New York, NY, USA, 1998. [Google Scholar]
Figure 1. Model of crystal and mounting system with detail of the mounting contacts.
Figure 1. Model of crystal and mounting system with detail of the mounting contacts.
Sensors 23 00024 g001
Figure 2. Different heater areas for power distribution.
Figure 2. Different heater areas for power distribution.
Sensors 23 00024 g002
Figure 3. Temperatures of the supports’ spacers.
Figure 3. Temperatures of the supports’ spacers.
Sensors 23 00024 g003
Figure 4. (a) Electrode temperature mapping; (b) Temperature gradient.
Figure 4. (a) Electrode temperature mapping; (b) Temperature gradient.
Sensors 23 00024 g004
Figure 5. Temperature profile on the electrode edge.
Figure 5. Temperature profile on the electrode edge.
Sensors 23 00024 g005
Figure 6. Temperature gradient magnitude on the electrode edge.
Figure 6. Temperature gradient magnitude on the electrode edge.
Sensors 23 00024 g006
Figure 7. Thermal field of the complete model.
Figure 7. Thermal field of the complete model.
Sensors 23 00024 g007
Figure 8. Temperature map on the electrode: comparison between the measured thermal image (left) and the model predictions (right).
Figure 8. Temperature map on the electrode: comparison between the measured thermal image (left) and the model predictions (right).
Sensors 23 00024 g008
Figure 9. Temperature map on the heater: comparison between the measured thermal image (left) and the model predictions (right).
Figure 9. Temperature map on the heater: comparison between the measured thermal image (left) and the model predictions (right).
Sensors 23 00024 g009
Figure 10. The average thermal gradient (TG_ mean) on the electrode edge.
Figure 10. The average thermal gradient (TG_ mean) on the electrode edge.
Sensors 23 00024 g010
Figure 11. Frequency variation vs. thermal gradient.
Figure 11. Frequency variation vs. thermal gradient.
Sensors 23 00024 g011
Figure 12. Comparison between the relative frequency changes predicted by the model and the measured ones.
Figure 12. Comparison between the relative frequency changes predicted by the model and the measured ones.
Sensors 23 00024 g012
Figure 13. Temperature gradient vs. power dissipation.
Figure 13. Temperature gradient vs. power dissipation.
Sensors 23 00024 g013
Table 1. Material properties used in the model.
Table 1. Material properties used in the model.
ComponentsDensityYoung ModulusPoisson CoefficientCTEThermal Conductivity
[kg/m3][MPa] [1/°C][W/(m °C)]
QUARTZ264976,5000.177.10 × 10−66.2
AISI3168000193,0000.291.60 × 10−514
MACOR252066,9000.299.00 × 10−61.5
Table 2. Model errors: the difference between model-predicted temperatures and measured ones on the selected points on the electrode and heater.
Table 2. Model errors: the difference between model-predicted temperatures and measured ones on the selected points on the electrode and heater.
Electrode PointsError [°C]Heater PointsError [°C]
111−0.07
2−0.520.33
30.230.17
4040.34
50.4
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Magni, M.; Scaccabarozzi, D.; Saggin, B. Compensation of Thermal Gradients Effects on a Quartz Crystal Microbalance. Sensors 2023, 23, 24. https://doi.org/10.3390/s23010024

AMA Style

Magni M, Scaccabarozzi D, Saggin B. Compensation of Thermal Gradients Effects on a Quartz Crystal Microbalance. Sensors. 2023; 23(1):24. https://doi.org/10.3390/s23010024

Chicago/Turabian Style

Magni, Marianna, Diego Scaccabarozzi, and Bortolino Saggin. 2023. "Compensation of Thermal Gradients Effects on a Quartz Crystal Microbalance" Sensors 23, no. 1: 24. https://doi.org/10.3390/s23010024

APA Style

Magni, M., Scaccabarozzi, D., & Saggin, B. (2023). Compensation of Thermal Gradients Effects on a Quartz Crystal Microbalance. Sensors, 23(1), 24. https://doi.org/10.3390/s23010024

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