WG adaptor | Side a (mm) | Side b (mm) | Freq. Range (GHz) | Cutoff Freq. (GHz) |
R9LS | 247.7 | 123.8 | 0.70–1.15 | 0.606 |
R22LS | 109.9 | 54.6 | 1.72–2.60 | 1.365 |
R26LS | 86.4 | 43.2 | 2.20–3.30 | 1.736 |
R58LS | 40.7 | 20.2 | 4.90–7.05 | 3.686 |
Citation: Dianyuan Qi, Zongying Yu, Jing Zhao, Fangzhu Zou. A calibration of E-field probe based on tissue-equivalent liquid waveguide and evaluation of uncertainty[J]. Mathematical Biosciences and Engineering, 2019, 16(6): 7911-7920. doi: 10.3934/mbe.2019397
[1] | Daniele Bernardo Panaro, Andrea Trucchia, Vincenzo Luongo, Maria Rosaria Mattei, Luigi Frunzo . Global sensitivity analysis and uncertainty quantification for a mathematical model of dry anaerobic digestion in plug-flow reactors. Mathematical Biosciences and Engineering, 2024, 21(9): 7139-7164. doi: 10.3934/mbe.2024316 |
[2] | Ping Zhou, Yuting Zhang, Yunlei Yu, Weijia Cai, Guangquan Zhou . 3D shape measurement based on structured light field imaging. Mathematical Biosciences and Engineering, 2020, 17(1): 654-668. doi: 10.3934/mbe.2020034 |
[3] | Enrico Bertino, Régis Duvigneau, Paola Goatin . Uncertainty quantification in a macroscopic traffic flow model calibrated on GPS data. Mathematical Biosciences and Engineering, 2020, 17(2): 1511-1533. doi: 10.3934/mbe.2020078 |
[4] | Yongjie Zhao, Sida Li, Zhiping Huang . TI-ADC multi-channel mismatch estimation and calibration in ultra-high-speed optical signal acquisition system. Mathematical Biosciences and Engineering, 2021, 18(6): 9050-9075. doi: 10.3934/mbe.2021446 |
[5] | Chuangxia Huang, Shumin Tan, Le Li, Jie Cao . Interbank complex network and liquidity creation: Evidence from European banks. Mathematical Biosciences and Engineering, 2023, 20(11): 19416-19437. doi: 10.3934/mbe.2023859 |
[6] | Yan Wang, Guichen Lu, Jiang Du . Calibration and prediction for the inexact SIR model. Mathematical Biosciences and Engineering, 2022, 19(3): 2800-2818. doi: 10.3934/mbe.2022128 |
[7] | Peng-Yeng Yin, Chih-Chun Tsai, Rong-Fuh Day, Ching-Ying Tung, Bir Bhanu . Ensemble learning of model hyperparameters and spatiotemporal data for calibration of low-cost PM2.5 sensors. Mathematical Biosciences and Engineering, 2019, 16(6): 6858-6873. doi: 10.3934/mbe.2019343 |
[8] | Qinghua Liu, Jing Wang, Keying Song . Displacement identification of oil storage tank and calibration of tank capacity table model based on indicator function integral. Mathematical Biosciences and Engineering, 2024, 21(2): 1872-1883. doi: 10.3934/mbe.2024082 |
[9] | Zheng Dai, I.G. Rosen, Chuming Wang, Nancy Barnett, Susan E. Luczak . Using drinking data and pharmacokinetic modeling to calibrate transport model and blind deconvolution based data analysis software for transdermal alcohol biosensors. Mathematical Biosciences and Engineering, 2016, 13(5): 911-934. doi: 10.3934/mbe.2016023 |
[10] | Floriane Lignet, Vincent Calvez, Emmanuel Grenier, Benjamin Ribba . A structural model of the VEGF signalling pathway: Emergence of robustness and redundancy properties. Mathematical Biosciences and Engineering, 2013, 10(1): 167-184. doi: 10.3934/mbe.2013.10.167 |
SAR testing is a measure, applied to the wireless communication devices, to quantify the exposure to radiofrequency (RF) fields [1,2,3], and the peak spatial-average SAR is commonly used for safety compliance of the RF exposure. The core component of the current SAR measurement system is the Dosimetric E-field Probe that allows to be read when immersed in the tissue-equivalent liquid. In order to confirm its reliability, the E-field probe needs to be calibrated regularly, and the calibration is usually required to be performed with the E-field probe as a whole (such complexity is because of its encapsulation). IEC 62209-1 has introduced the easier way to calibrate the probe, but there is no detail given. The method in this paper managed to simplify the calibration process while maintaining precision.
More specifically, for the liquid waveguide with a known E-filed, there will be a change in its E-field distribution if the dosimetric E-field probe inserts into the waveguide. By analyzing the difference between these distributions, the conversion factor can be gained effectively.
An E-field probe generally consists of three orthogonally arranged dipoles directly loaded with a detector diode and connected with resistive lines to the readout electronics. Probe robustness and protection of the components are achieved by mounting the sensors on a core and encapsulating the entire probe.
The total electrical field can be evaluated by the root sum square (RSS) of the three E-field components. In the square-law region of the diode response, the sensor output voltage is proportional to the mean square of the corresponding field component. Also, linearization is required within the dynamic region, as the output voltage will be compressed beyond this square-law region. Manufacturing tolerances between the sensors and diodes may produce different sensitivities for each sensor [3], therefore, for the purpose of having accurate measurement data, it is necessary to have the sensitivities evaluation during calibration.
The total field shall be evaluated according to the Eq. (1).
|E|2=3∑i=1|Ei|2=3∑i=1fi(Vi)ηiγi …… | (1) |
Where fi(Vi) is the linearization function of the rectified sensor signal Vi, ηi is the absolute sensitivities [μV/(v/m)2] of the dipole sensors in free space, ϒiis the ratios of the sensitivities of the probe sensors in media to their sensitivities in free space.
The separation of the E-field probe sensitivity into two factors (ηi and ϒi) allows the probe to be calibrated in free-space and tissue-equivalent liquid, respectively. That is, the frequency response and linearity can be estimated in free space, while it is able to make the assessment of conversion factor (ConvF), isotropy and boundary compensation in the tissue-equivalent liquid. In this way, the complexity of calibration can be considerably reduced.
According to the IEC 62209-1, the liquid waveguide was used as the calibration platform in this paper (see table 1). The liquid waveguide is a self-contained system in which the cross-sectional field distributions are not dependent upon reflections.
WG adaptor | Side a (mm) | Side b (mm) | Freq. Range (GHz) | Cutoff Freq. (GHz) |
R9LS | 247.7 | 123.8 | 0.70–1.15 | 0.606 |
R22LS | 109.9 | 54.6 | 1.72–2.60 | 1.365 |
R26LS | 86.4 | 43.2 | 2.20–3.30 | 1.736 |
R58LS | 40.7 | 20.2 | 4.90–7.05 | 3.686 |
A dielectric spacer should be of low-loss material (δ < 0.02) with low water absorption characteristics (<0.1%). The dielectric spacer is placed at a distance > λ from the feeding coupler in order to provide a good impedance match (>10 dB return loss) (see Figure 1). The dimensions of the spacer block should be adjusted to the dimensions of the actual waveguide used and the edges sealed after placing it in the final location (e.g. with silicon). The medium depth must be three times larger than skin depths in order to have negligible reflections at the upper surface of the liquid.
For accurate power assessments (i.e. high-end power meter and sensors), the most critical parameter in the calibration procedure is the assessment of the net RF power dissipated into the waveguide. This implies the precise measurement of the incident power (Pf), adaptor loss (transloss) and return loss at the waveguide input port (γ). The net power inside the lossy liquid can then be calculated as the Eq. (2):
Pnet=Pfw−Pbw=Pfw⋅(1−|r|2) …… | (2) |
Figure 2 shows the recommended setup. The P1 and attenuator placed in front of the P1 sensor head measure the forward power at the location of the system check dipole connector. The signal generator is adjusted for the desired forward power at the dipole connector and the power meter P2 is read at that level. After connecting the cable to the dipole, the signal generator is readjusted for the same reading at power meter P2. If the signal generator does not allow a setting in 0.01 dB steps, the remaining dierence at P2 must be noted and considered in the normalization of the system check results. Power meter P3 is optional. If available, it will provide a constant monitoring of the backward power.
The total power absorbed by the lossy liquid can be accurately determined by measurement of the forward and reflected powers. Because of the low cutoff frequency in the lossy media, the wave is almost a TEM wave such that the decay in the liquid can be assessed by measuring along a vertical axis from the liquid surface and exponential curve fitting. As all power entering the lossy liquid is absorbed by the liquid, the SAR along the vertical z-axis can be determined as the Eq. (3):
SAR(z)=4(Pfw−Pbw)ρabδe(−2z/δ) …… | (3) |
where a, b is the cross-sectional area of the waveguide, Pfw and Pbw the forward and backward (reflected) powers inside the waveguide, δ the penetration depth inside the lossy liquid, and ρ the density of the lossy liquid conventionally assumed to be 1000 kg/m3.
Filling the liquid waveguide with tissue-equivalent liquid at the calibrated frequency and inserting the E-field probe vertically in the center of the waveguide. The probe shall be rotated around its axis from 0° to 360° by the probe positioner and measurement data (Eθ) is recorded in steps of 10°. After the rotation, the probe is moved to the max data of the angle, then along with the Z-axis, the measurement data is recorded in steps of 1mm [4].
By comparing the measurement data from the Z-scan and the target value SARcal by every step, it was getting the minimum of SAR values (eq 4):
min=z∑[SAR(z,l,f)−SARcal(z,l,f)]2 …… | (4) |
In the liquid waveguide, as the waveguide dimensions (a, b) and the net power absorbed in the liquid (Pnet) are known, Eq. (4) can be reformulated using the carrier frequency (f) and the dielectric (l) parameters of the medium. ConFi(f, l) and δ(f, l) then are expressed as Eq. (5):
m=z∑[V2iConvFi(l,f)⋅σρ⋅1000⋅cos(θ)2−4Pnet(l,f)abδ(l,f)e(−2z/δ(l,f))]2…… | (5) |
Moreover, the maximum values are below the linearized range and that lowest values are sufficiently above the noise floor, e.g., > 20 dB. As the measured values close to the waveguide layer are distorted by the boundary effect, only values that were measured at positions larger than twice the probe tip diameter from the surface of the waveguide spacer should be considered. The boundary parameters β(f, l) and α(f, l) are determined by the Eq.6:
m=z∑[SAR(z,l,f)(1+β(f,l)e(−zα(f,l))cos(π−zλ))−SARcal(z,l,f)]2…… | (6) |
A dosimetric E-field probe is chosen to perform calibration. The calibration frequency is 900MHz by the head liquid, so the liquid waveguide filled with 900 MHz head tissue-equivalent liquid, which conforms to the standard requirements. According to Table 1, the waveguide chosen in this paper is R9LS.
In conformity with the calibration procedure, after the E-field probe inserted into R9LS, the Z-scan job retracts the probe from the waveguide and acquires the voltage at each step. Since the reference filed of the waveguide is known, its relation to the measured field showing in figure 3 can be drawn. it is straightforward to find a considerable deviation between them.
According to chapter Ⅳ, the ConvF is imported to the measured field, then the deviation is shown in figure 4.
For the Z-scan, after the measurement data had been compensated by ConvF, the error of the E-field will be less than ± 0.1 dB, which meets the standard requirements. Hence, as shown in figure 5, it has been proved that this procedure is feasible and practicable.
In accordance with the procedures in IEEE 1528, IEC 62209-1 and the uncertainty under the actual circumstance, the uncertainty components of the ConvF in this calibration please see table 2.
No. | Source of the uncertainty |
1 | RF power |
2 | Repeatability |
3 | Probe mechanical positioning |
4 | Field uniformity |
5 | Probe linearity |
6 | Tissue-equivalent liquid |
1) According to the calibration certificate of the power meter, the uncertainty u1a = 3.4% (k = 2), so the standard uncertainty u1a' = 1.7% (k = 1).
2) The mismatch error between the power meter and cable is smaller than 0.05 dB, so the uncertainty introduced by mismatch error α1 = 1.2%, U-shape distribution, the coverage factor k=√2,
u1b=1.2%√2=0.9% |
3) The uncertainty of net power delivered to the waveguide, the maximum allowable drift is 0.20 dB, so the uncertainty introduced by mismatch error α3 = 4.7%, U-shape distribution, the coverage factor k=√2,
u1c=4.7%√2=3.3% |
4) Due to the u1a', u1b and u1c are uncorrelated, so uncertainty introduced by the RF power:
u1=√u21a+u21b+u21c≈3.5% |
For this source, the calibration system has been repeated 10 times for the ConvF at 900 MHz. The data please see table 3. According to table 3, the u2 = 0.3%
NO. | Value |
1 | 9.12 |
2 | 9.18 |
3 | 9.16 |
4 | 9.15 |
5 | 9.12 |
6 | 9.13 |
7 | 9.13 |
8 | 9.11 |
9 | 9.17 |
10 | 9.20 |
Average | 9.15 |
standard deviation | 0.03 |
The mechanical tolerance of the field probe positioner can introduce probe positioning uncertainty. The resulting SAR uncertainty is assessed by comparing the SAR obtained according to the specifications of the probe positioner with respect to the actual position defined by the geometric center of the probe sensors. According to the standard IEC 602209-1, The tolerance is determined as:
ErrorSAR(%)=100×dssδ/2 |
The specified repeatability of the Staubli robot arm used in DASY system is ±25 μm. The absolute accuracy for short distance movements is better than ±0.1 mm, i.e., below 3 GHz the ErrorSAR(%) is better than 0.4%, rectangular distribution, coverage factor, k=√3,
u3=0.3% |
At the geometric center of the R9LS waveguide, along with the vertical and horizontal respectively, the field data are recorded in steps of 1mm, and they can be seen in table 4. According to table 4, the average is 1631.8 mV/m, the standard deviation is 4.045, so the u4 = 0.3%
Vertical distance from the center point (mm) | Field Value (mV/m) | Horizontal distance from the center point (mm) | Field Value (mV/m) |
-4 | 1642 | -4 | 1630 |
-3 | 1636 | -3 | 1632 |
-2 | 1638 | -2 | 1633 |
-1 | 1635 | -1 | 1633 |
0 | 1634 | 0 | 1634 |
1 | 1633 | 1 | 1632 |
2 | 1632 | 2 | 1626 |
3 | 1632 | 3 | 1628 |
4 | 1632 | 4 | 1622 |
5 | 1634 | 5 | 1620 |
According to the manufacturing characteristic of the E-field probe, the maximum deviation is ± 0.2 dB, rectangular distribution, coverage factor, k=√3,
u5=1.4% |
According to the IEC 62209-1, the maximum allowable drift of conductivity is 5 %, rectangular distribution, coverage factor k=√3, so u6a = 2.9 %.
Likewise, the maximum allowable drift of permittivity is 5 %, rectangular distribution, coverage factor k=√3, so u6b = 2.9 %.
The u6a and u6b are uncorrelated, so uncertainty introduced by the Tissue-equivalent liquid:
u6=√u26a+u26b≈4.1% |
Concluding the contents of the chapters 6.1 to 6.7, the standard combined uncertainties of ConvF is:
uc=√6∑i=1u2i≈5.8% |
So the expanded uncertainties of ConvF are:
U=kuc=11.6%,k=2 |
The convF of E-filed probe has been obtained through the method proposed in this paper. In addition, the result values show no excess of the standard requirement. The calibration procedure, therefore, is achievable and practicable.
However, the calibration of the E-field probe involves numerous parameters. Further studies and researches will be necessary.
The authors declare no conflict of interest in this paper.
[1] | IEEE Std C95.1-2005, IEEE standard for safety levels with respect to human exposure to radio frequency electromagnetic fields, 3 kHz to 300 GHz. 2005. |
[2] | IEEE Std 1528, IEEE recommended practice for determining the peak spatial-average specific absorption rate (SAR) in the human head from wireless communications devices: Measurement techniques, Sept. 2013. |
[3] | IEC Std 62209-1, Measurement procedure for the assessment of specific absorption rate of human exposure to radio frequency fields from hand-held and body-mounted wireless communication devices-Part 1: Devices used next to the ear (Frequency range of 300 MHz to 6 GHz). 2016. |
[4] | IEC Std 62209-2, Human exposure to radio frequency fields from hand-held and body-mounted wireless communication devices-Human models, instrumentation, and procedures-Part 2: Procedure to determine the specific absorption rate (SAR) for wireless communication devices used in close proximity to the human body (frequency range of 30 MHz to 6 GHz). 2010. |
[5] | ICNIRP, International Commission on Non-Ionizing Radiation Protection. Guidelines for limiting exposure to time-varying electric, magnetic, and electromagnetic fields (up to 300 GHz). 1998. |
[6] | ISO/TS 14253-3, Geometrical product specifications (GPS)-Inspection by measurement of workpieces and measuring equipment-Part 3: Guidelines for achieving agreements on measurement uncertainty statements. 2002. |
[7] | K. Poković, "Advanced electromagnetic porbe for near-field evaluations," ph.D. dissertation, Swiss Federal Institute of Technology, Zurich, Switzerland, 1999. |
[8] | H. I. Bassen and G. S. Smith, Electric field probes-A Review, IEEE Transact. Anten. Propagat., 31 (1983), 710–718. |
[9] | S. Pirbhulal, H. Y. Zhang, W. Q. Wu, et al., Heartbeats based biometric random binary sequences generation to secure wireless body sensor networks, IEEE Transact. Biomed. Eng. 65 (2018), 2751–2759. |
[10] | W. Q. Wu, S. Pirbhulal, H. Y. Zhang, et al., Quantitative assessment for self-tracking of acute stress based on triangulation principle in a wearable sensor system, IEEE J. Biomed. Health Inform. 99 (2018), 1. |
[11] | S. Pirbhulal, H. Y. Zhang, S. Mukhopadhyay, et al., An efficient biometric-based algorithm using heart rate variability for securing body sensor networks, Sensors, 15 (2015), 15067–15089. |
[12] | S. Pirbhulal, H. Y. Zhang, M. E. E. Alahi, et al., A novel secure IoT-based smart home automation system using a wireless sensor network, Sensors, 17 (2017), 69. |
[13] | W. Q. Wu, H. Y. Zhang, S. Pirbhulal, et al., Assessment of biofeedback training for emotion management through wearable textile physiological monitoring system, IEEE Sensors J., 15 (2015), 7087–7095. |
[14] | H. Viltanen, J. Keshvari and R. Lappalainen, Interaction of radio frequency electromagnetic fields and passive metallic implants-A brief review, Bioelectromagnetics, 27(2006), 431–439. |
[15] | V. Hombach, K. Meier, M. Burkhardt, et al., The dependence of EM energy absorption on human head modeling at 900 MHz, IEEE Transact. Microwave Theor. Techn, 44 (1996), 1865–1873. |
[16] | A. Kyriacou, A. Christ, E. Neufeld, et al., Maximum Local Tissue Temperature Increase by Implanted Medical Devices Due to Electromagnetic Field Exposure, Bioelectromagnetics Annual Meeting, Seoul, Korea, June 2010. |
[17] | L. Aberbour, B. Derat and A. Cozza, Analysis of vector E-field sensor array for real-time SAR assessment, Proceedings of European Conference on Antenna and Propagation (EuCAP), 2013. |
[18] | K. Poković, T. Schmid and N. Kuster, Millimeter-resolution E-field probe for isotropic measurement in lossy media between 100 MHz and 20 GHz, IEEE Transact. Instrum. Meas., 49 (2000), 873–878. |
[19] | T. Iyama, K. Kiminami and T. Onishi, Applicability of three-axis electro-optic (EO) probe for specific absorption rate (SAR) measurement. IEICE Transact. Commun., 92 (2009), 1414–1417. |
[20] | O. Merckel, G. Fleury and J. C. Bolomey, Rapid SAR measurement via parametric modeling. In 5th International Congress of the European Bioelectromagnetics Association, Helsinki, Finland, September 2001. |
[21] | C. C. Davis and Q. Balzano, The international intercomparison of SAR measurements on cellular telephones, IEEE Trans. Electromagn. Compat., 51 (2009), 210–216. |
WG adaptor | Side a (mm) | Side b (mm) | Freq. Range (GHz) | Cutoff Freq. (GHz) |
R9LS | 247.7 | 123.8 | 0.70–1.15 | 0.606 |
R22LS | 109.9 | 54.6 | 1.72–2.60 | 1.365 |
R26LS | 86.4 | 43.2 | 2.20–3.30 | 1.736 |
R58LS | 40.7 | 20.2 | 4.90–7.05 | 3.686 |
No. | Source of the uncertainty |
1 | RF power |
2 | Repeatability |
3 | Probe mechanical positioning |
4 | Field uniformity |
5 | Probe linearity |
6 | Tissue-equivalent liquid |
NO. | Value |
1 | 9.12 |
2 | 9.18 |
3 | 9.16 |
4 | 9.15 |
5 | 9.12 |
6 | 9.13 |
7 | 9.13 |
8 | 9.11 |
9 | 9.17 |
10 | 9.20 |
Average | 9.15 |
standard deviation | 0.03 |
Vertical distance from the center point (mm) | Field Value (mV/m) | Horizontal distance from the center point (mm) | Field Value (mV/m) |
-4 | 1642 | -4 | 1630 |
-3 | 1636 | -3 | 1632 |
-2 | 1638 | -2 | 1633 |
-1 | 1635 | -1 | 1633 |
0 | 1634 | 0 | 1634 |
1 | 1633 | 1 | 1632 |
2 | 1632 | 2 | 1626 |
3 | 1632 | 3 | 1628 |
4 | 1632 | 4 | 1622 |
5 | 1634 | 5 | 1620 |
WG adaptor | Side a (mm) | Side b (mm) | Freq. Range (GHz) | Cutoff Freq. (GHz) |
R9LS | 247.7 | 123.8 | 0.70–1.15 | 0.606 |
R22LS | 109.9 | 54.6 | 1.72–2.60 | 1.365 |
R26LS | 86.4 | 43.2 | 2.20–3.30 | 1.736 |
R58LS | 40.7 | 20.2 | 4.90–7.05 | 3.686 |
No. | Source of the uncertainty |
1 | RF power |
2 | Repeatability |
3 | Probe mechanical positioning |
4 | Field uniformity |
5 | Probe linearity |
6 | Tissue-equivalent liquid |
NO. | Value |
1 | 9.12 |
2 | 9.18 |
3 | 9.16 |
4 | 9.15 |
5 | 9.12 |
6 | 9.13 |
7 | 9.13 |
8 | 9.11 |
9 | 9.17 |
10 | 9.20 |
Average | 9.15 |
standard deviation | 0.03 |
Vertical distance from the center point (mm) | Field Value (mV/m) | Horizontal distance from the center point (mm) | Field Value (mV/m) |
-4 | 1642 | -4 | 1630 |
-3 | 1636 | -3 | 1632 |
-2 | 1638 | -2 | 1633 |
-1 | 1635 | -1 | 1633 |
0 | 1634 | 0 | 1634 |
1 | 1633 | 1 | 1632 |
2 | 1632 | 2 | 1626 |
3 | 1632 | 3 | 1628 |
4 | 1632 | 4 | 1622 |
5 | 1634 | 5 | 1620 |