Variability of coil inductance measurements inside an interleaving structure.

The present article is a continuation of our previous works on the anomalies in the decay rates and the capacitance measurements inside a modified Faraday cage. Here we present the anomalous variations in the measurements of inductance when a coil is placed inside an interleaving structure of metal and organic material. The fluctuation in the inductance values was found to be in the range - 0.007 to 0.018 mH. Additionally, it was observed that the temperature coefficient which was 0.0062 mH/°C initially jumped to two distinct levels, 0.0095 mH/°C and 0.0145 mH/°C, respectively. A multiple factor analysis of our results revealed a very strong correlation (R2 = 95.2%) between the inductance and the combination of the temperature and the relative humidity both measured inside the cage, next to the inductor.

Introduction

It is known that many physical and chemical experiments are frequently protected from external electromagnetic disturbances by means of the well-known Faraday cages[1]. A variant of these cages is an interleaving structure, built with alternating layers of metal and organic material[2], but unlike the conventional Faraday cages this kind of structures seem to cause modification of measurements taken inside it. So, in one of the many experiments carried out with such device, strong variations in the measurements of decay rate of Ra-226 were claimed to be observed[2]. Such a strange result was at variance with the standard exponential decay. To our knowledge, no scientific literature about the repetition of such experiment was found. In this respect, we decided to perform several experiments on decay rates measurements of several nuclides (Ra-226, Tl-204 and Sr90/I90) using a Geiger–Müller (GM) counter along with a similar enclosure. Besides, we also obtained the spectrum of Cs-137 by means of a multichannel analyser. The question on whether the GM counter tube was stable enough or not was addressed in several ways[3]. As discussed in Ref.[4], the GM counter we used is stable and robust enough for such task. Besides, the GM counter tube environmental test provided by the manufacturer was considered (see the Appendix in Supplementary Material). Anomalous results when measuring decay rates of radionuclides were observed[3]. For example, the Ra-226 decay rates measurements increased up to 1.2%, and up to 0.8% for other nuclides. Regarding the Cs-137 spectrum measurement, the result was a displacement of the spectrum to lower energies and an increment of the photopeak of about 4%.

One hypothesis about the reason for such a variability was that the structure caused perturbations in the electronics of the detectors. Taking into account that one parameter that determines the duration of the pulses generated by the GM counter and the MCA is the time constant of the measuring circuit[5], another experiment was carried out by using a RC low-pass filter. The measurements of the time constant yielded an increase of about 5%.

The next step was the capacitance measurement of an ultra-stable capacitor in the inside and outside of our interleaving structure. A change in the capacitance of 0.8–1% was observed. The problem became even more puzzling when Scholkmann found that our data could be correlated with cosmic weather indices like the geomagnetic activity (GMA) and with the cosmic ray induced neutron counts (CRN)[6]. To address this problem, we first looked for the circumstances under which these anomalous behaviours took place[7,8]. Subsequently in a work by Pommè and Pelzcar[9] it was concluded that correlation with ambient humidity should also be considered. This was an indication for the improvement of our work and consequently, in a subsequent experiment[10] we used outdoor values of the relative humidity provided by the Spanish State Agency AEMET[11]. A correlation between capacitance Cp and GMA (R2 = 36%) was found.

The present work is a continuation of the previous work but, with the following variations: first, measurements of inductance are carried out for a coil placed inside and outside the interleaving structure (as the coil is a very common electronic components). Second, we looked for correlations, not only with the GMA and the CRN, but also with the environmental parameters like the temperature and the relative humidity; however, this time these variables were measured inside the cage beside the device under test (DUT). The correlated variables were the inductance (L, dependent variable), and two other variables namely, the indoors ambient temperature (Tii) and the indoors relative humidity (RHii%). This time there was no correlation, neither with the Dcx nor with the CRN, but this happened coincidentally with the low values of the GMA and the CRN, as described in[7].

The outline of this work is as follows. In the next section, the “Materials and methods” used in this work are explained. Subsequently, a section with the “Results and discussion” is presented. Finally, some “Conclusions” are drawn.

Materials and methods

The shielding surrounding the DUT was formed by an external modified Faraday cage (MFC), but this time two additional shielding structures were used one inside another (as shown Fig. 1a). The dimensions of the outer cage are 50 × 50 × 50 cm3. Each side of this cage contains five sheets of rough alumina and five sheets of self-adhesive elastomeric material[12].

Figure. 1

Inductance measurements inside the interleaving structure: (a) external cage contains the intermediate one; (b) 0.25 × 0.25 × 0.25 cm3 cage which contains an even smaller cage, (c) LCR measuring instrument.

The intermediate cage dimensions are 25 × 25 × 25 cm3 with layers of galvanized iron of 1 mm thickness and 2 mm wide cork sheeting. Finally, the innermost layer is a carton box of 9 × 9 × 13 cm3, in which we placed the inductance. This small box was uniformly covered by a metallic mesh for electromagnetic isolation.

As before, the inductance was measured by the same instrument i.e., a 4263B HP LCR meter[13] and the coil was a 10 mH ± 10% niquel-zinc ferrite core inductor (SBCP-80HY103HB)[14].

The excitation signal of the LCR was set to 1 kHz and 1 V as before. Under these measurement conditions, the corresponding accuracy of the equipment was 0.08%[13]. However, we performed an analysis of variance (ANOVA) with the data to determine which results were statistically different. Figure 1 shows several images of the setup, i.e., two of the cages used and the measuring instrument.

Regarding the testing of the isolation effectiveness of the structure, a test with a mobile phone yielded that the power of the signal was reduced by 40% in the biggest box. This reduction was equivalent to 48 dBm, which is an acceptable result. As stated by some authors, to consider a Faraday cage as suitable and adequate, this value should ideally be close to 50 dBm[15].

Regarding the correlations with the space weather and considering the study of Pommé and Pelzcar[9], a temperature meter and a humidity meter were placed next to the inductance. The equipment chosen was the DHT22 device[16], which was compatible with Arduino Uno Microcontroller and had a sensitivity of 0.1 °C.

Data acquisition

Since only one LCR meter was available (Fig. 1c), inductance (L, in mH) and temperature (T, in °C) values were taken outside and inside the box at different times. First, the data were taken inside the box from Oct 27 to Nov 3, 2020. Later, data were taken outside the box from Nov 3 to Nov 11, 2020. Inductance and temperature were measured every second and every minute respectively. In this way, we could build the time series of the means corresponding to 1 min and to 1 h.

Geomagnetic activity (GMA) and cosmic ray activity (CRA) data were obtained from the same sources as mentioned in previous work[7]. In previous papers we used for GMA the classic storm index Dcx; its equivalent is the Dst index, but here we replaced it by the SYM-H geomagnetic index as it has a 1-min time resolution compared to the 1-h time resolution of Dst[17,18]. The SYM-H index is essentially the same as the Dst index with a different time resolution.

Results and discussion

We measured the inductance when the system was placed indoor, but the DUT was outside the cage as shown in Fig. 2a. As expected, the inductance follows the same trend as the temperature variation (R2 = 94%).

Figure. 2

Inductance measurements vs temperature. (a) Plot of the inductance outside the box (from 2020 November 3 to 11) along with temperature beside the coil. The inductance follows the temperature variations. Time span of about two days for better visualize the synchronism of inductance and temperature (b) plot of the inductance inside the box (from 27 Oct to Nov 3); as expected, inductance deviates from temperature variations.

However, the inductance measurements did not follow the same variation with the temperature inside the box when compared with the measurements taken outside the box. (Fig. 2b)

Table 1 shows the values of the relative maxima and minima from the graphs of both the inductance and the temperature outside the cage (measurements are taken next to the inductor) (Fig. 2a). The following quantities are included in Table 1.

Table 1

Values of relative maxima and minima of temperature (T) and inductance (L) of the curves in Fig. 2a. All the peaks with the same temperature show the same inductance values. Results for the valleys are slightly different.

Outside
TM1,o	TM2,o	TM3,o	TM4,o	TM5,o	TM6,o	TM7,o	TM8,o	TM9,o
30.9	30.9	30.9	30.9	30.9	30.9	30.9	30.9	30.9
LM1,o	LM2,o	LM3,o	LM4,o	LM5,o	LM6,o	LM7,o	LM8,o	LM9,o
11.342	11.342	11.342	11.342	11.342	11.342	11.342	11.342	11.342
Tm1,o	Tm2,o	Tm3,o	Tm4,o	Tm5,o	Tm6,o	Tm7,o	Tm8,o	Tm9,o
29.3	29.4	29.3	29.4	29.4	29.3	29.3	29.4	29.4
Lm1,o	Lm2,o	Lm3,o	Lm4,o	Lm5,o	Lm6,o	Lm7,o	Lm8,o	Lm9,o
11.326	11.326	11.326	11.326	11.326	11.326	11.325	11.326	11.326

TMj,i (LMj,i), the maximum (M) temperature (T) (inductance (L)) value of peak j, (j = 1, 2,…, 9), inside (i) the box. It is given in Celsius degrees (millihenries).

Tmj,i (Lm): the minimum (m) temperature (inductance L) value of valley j inside (i) the box.

In Table 1, the same inductance values (LMj,o = 11.342 mH) correspond to the same temperature values (TMj,o = 30.9 °C). Similar observation applies for valley values (third and fourth row in Table 1). However, when the coil was placed inside the cage, we observed some unexpected results, which are presented in Table 2, where the values of relative maxima of temperature (T) and inductance (L) correspond to the graphs in Fig. 2b. Specifically, the following peaks have the same value: TM3,i = TM5,i = TM7,i = TM8,i = 30.9 °C, whereas the corresponding inductance values are no longer equal to LMj,o = 11.342 mH (which was the value outside), and are also different from each other (2nd row in Table 2). The differences with LMj,o = 11.342 mH are shown in the third row. Moreover, the inductance LMj,o = 11.342 mH (which remained constant all the time when the coil was placed outside the box), fluctuated between 11.334 and 11.346 mH when placed inside. Similar fluctuations can be observed in the inductance corresponding to the temperatures Tm2,i = Tm4,i = 29.2 °C and to Tm3,i = Tm6,i = Tm7,i = 29.3 °C.

Table 2

Values of relative maxima and minima of the curves in Fig. 2b.

Inside
TM1,i	TM2,i	TM3,i	TM4,i	TM5,i	TM6,i	TM7,i	TM8,i
30.6	30.5	30.9	31.0	30.9	30.7	30.9	30.9
LM1,i	LM2,i	LM3,i	LM4,i	LM5,i	LM6,i	LM7,i	LM8,i
11.342	11.345	11.346	11.341	11.337	11.334	11.335	11.334
		+ 0.004		− 0.005		− 0.007	− 0.008
Tm1,i	Tm2,i	Tm3,i	Tm4,i	Tm5,i	Tm6,i	Tm7,i	Tm8,i
28.6	29.2	29.3	29.2	29.1	29.3	29.3	29.4
Lm1,i	Lm2,i	Lm3,i	Lm4,i	Lm5,i	Lm6,i	Lm7,i	Lm8,i
11.326	11.331	11.327	11.322	11.318	11.319	11.320	11.318

The observed behavior was very different in each scenario. In the first one (outside), it was clearly observed that the only variable that affected the inductance variation was the temperature. Thus, the temperature and the inductance increased and decreased at the same time. However, in the second case (inside), in addition to this phenomenon, a different behavior was apparent. This can be better understood by taking the data of Fig. 2b and rearranging them from minimum to maximum values of the temperature and performing an ANOVA analysis of the rearranged data[19], which is shown in Fig. 3. This shows how the inductance increases with the temperature both outside and inside the MFC cage. The uncertainty bars correspond to the least significant difference (LSD), which were calculated by the statistical program Statgraphics[20]. This could be considered as an estimation of the inductance temperature coefficient (TCL).

Figure. 3

Inductance temperature coefficient. (a) Outside the box. (b) Inside the box. The increase from 27.5 to 28.5 °C took place during the days 4th and 5th of November 2020.

As expected, when the coil was placed outside the cage, TCL was positive and practically linear, with a slope of about 0.0062 mH/°C (Fig. 3a). However, when the coil was placed inside the cage (Fig. 3b), the obtained data lacked the characteristic linearity of the first, which agrees with the greater variability of Fig. 3b. This time the inductance temperature coefficient shows high variation. The first ascending part of the graph (corresponding to the low temperatures) has an approximate slope of 0.0145 mH/°C. This is followed by a negative TCL of − 0.0302 mH/°C. Finally, the slope of the second ascending part is approximately 0.0095 mH/°C.

By extrapolation of the temperature values in Fig. 3a to 27.5 °C (which is the lowest temperature in Fig. 3b) a straight line can be obtained to be used as a reference. Then, the maximum difference between the inductance values in Fig. 3b and corresponding to the same temperature of the straight line in Fig. 3b can be approximately 0.025 mH (at 28.5 °C).

Besides, it is worth noting that in Fig. 3b at 28 °C the inductance is 11.34 mH and RHii ~ 35% whereas at 31 °C, L = 11.34 mH but RHii ~ 29.6%.

Correlations

In addition to the above results, correlation between the inductance L and the other independent variables were calculated and analysed. The independent variables were:

Temperature outdoor (To); temperature indoor inside the box (Tii).

Relative humidity outdoor (RHo); relative humidity indoor, measured inside the box (RHii).

Geomagnetic activity GMA (represented by Dcx) and,

Cosmic ray activity (CRA) (represented by cosmic neutrons counts, N (cps)).

The values of the inductance L along with Dcx and N are presented in Fig. 4a,b on an hourly basis. This kind of variability shown by both Dcx and N was observed in[7] to be coincident with low or no-correlation between these variables and the magnitudes under study. Specifically, Dcx remains between − 25 and + 15 nT (low level Dcx), and N remains lower than 70 cps.

Figure 4

(a) Inductance versus Dcx, the latter ranging between − 25 and + 15 nT; (b) Inductance versus N, where N < 70 cps. Plots in (a) and (b) are obtained on an hourly basis (means were calculated every hour). The variability of Dst and N resembles the circumstances when there was lack of correlation with the DUT.

On the other hand, the contribution of each variable (and/or the combination of variables) to the total variability of the inductance can be determined by simple or multiple regression analysis[21,22] applied to the time series. The time series of L consists of the mean values taken every hour (or every minute). The first row in Table 3 shows a multiple regression coefficient of 95.4% between L and the variables considered so far when the setup is placed inside the cage, namely: Tii, RHii, Dcx (equivalent to SYM_H) and N. A regression coefficient of 95.2% between L and Tii and RHii is obtained. Finally, the correlation coefficient is small when one considers only Dcx and N.

Table 3

Correlation between the different variables studied.

	L vs Tii, RHii, Dcx, N	R2 (%)
1	L = L(Tii, RHii, Dcx, N)	95.4
2	L = L(Tii, RHii)	95.2
3	L = L(Dcx, N)	22.8

Conclusions

The premise of this study is that the measurements of inductance of a coil can vary linearly with temperature and other external variables involved. This occurs mainly due to the increase in the resistance of the coil itself as the temperature increases. The resistance of copper wire per unit length (main material used in the construction of these electronic components) increases considerably with temperature, which causes a higher voltage to be needed to obtain the same current passing through it. Under this premise, it would have expected that the inductance value of the coil would have varied uniformly with the temperature. However, this happened only when it was measured outside the MFC. Within it, the performance was anomalous, as in the case of the capacitor and of the nuclides mentioned in[1]. The result was the observation of a variability of 0.025 mH, which is much greater than instrument accuracy (0.008 mH), and much greater than the LSD intervals calculated by the statistical program. Moreover, it was seen that the same inductance values could be measured at two different temperatures. The inductance variations were analyzed for correlations with the independent variables Tii, RHii, Dcx and N, thus obtaining a correlation coefficient of R2 = 95.4%. In this specific experiment, under the current circumstances, we found a weak correlation between L and the space weather variables Dcx and N. Therefore, almost the same correlation coefficient was obtained after erasing the superfluous independent variables (Dcx, N). The low RH values would not explain a dependence of L on RH.

Supplementary Information.