Analytical Approximations for the Maximum-to-mean Ratio of the E-field in a Reverberation Chamber: A Review

Qian Xu1, Rui Jia2, Lifei Geng2, Hao Guo2, and Yongjiu Zhao1

1College of Electronic and Information Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing, 211106, China
emxu@foxmail.com, yjzhao@nuaa.edu.cn

2State Key Laboratory of Complex Electromagnetic Environment Effects on Electronic and Information System, Luo Yang, China
jiarui315@163.com

Submitted On: December 17, 2021; Accepted On: August 18, 2022

Abstract

The expected value of the maximum value of the rectangular E-field is important in radiated susceptibility testing in a reverberation chamber. In this paper, different forms of equations for the maximum-to-mean ratio of the rectangular E-field are reviewed. Important derivations are summarized and detailed. It is interesting to note that some series which could be difficult to deal with from mathematics could be solved efficiently from physical point of view. The relationship between the independent sample number N and the parameters in generalized extreme value distribution is also given.

Index Terms: maximum-to-mean ratio, reverberation chamber, statistical electromagnetics.

1 INTRODUCTION

Reverberation Chambers (RCs) have been widely used in electromagnetic compatibility (EMC) testing [1]. Compared with an anechoic chamber, no radio absorber is used on the walls, thus a rich multipath environment can be created in an RC. By rotating the mechanical stirrers in an RC, statistically uniform and isotropic electromagnetic fields can be created. A typical measurement setup is illustrated in Fig. 1. For the radiated susceptibility testing in an RC, the maximum E-field is of interest [1], [2] as the device under test (DUT) could be malfunctioned by the maximum E-field.

In an RC, statistical electromagnetics are used to characterize the field properties. The probability density function (PDF) and the cumulative distribution function (CDF) of the maximum rectangular E-field have been given in [2]–[5]. Approximate analytical equations have also been given in [4], [5] for large independent sample number N. Different forms of equations exist, and the similarities and equivalencies among them have not been summarized.

images

Figure 1: A schematic plot of radiated susceptibility in an RC.

In this paper, we review the analytical equations for the expected value and the standard deviation of the maximum E-fields in literatures. Numerical results are calculated and compared with different forms of analytical equations.

2 THEORY

The mean rectangular E-field strength in an RC is easy to estimate. However, the maximum rectangular E-field strength is of interest in many scenarios for EMC testing. The maximum-to-mean ratio of the rectangular E-field in an RC has been well studied in [2]–[5]. The PDF and the CDF of the maximum rectangular E-field in an RC are [2]–[5]:

p(x)=Nxσ2[1-exp(-x22σ2)]N-1exp(-x22σ2), (1)

and:

F(x)=[1-exp(-x22σ2)]N, (2)

where N is the independent sample number, x represents the maximum value of the rectangular E-field x=|Ex|N, and N represents the maximum value from N samples. The PDF of the normalized |Ex|N can be obtained by setting σ=2/π, as the mean value of |Ex| is σπ/2. For different N, the PDF plots are illustrated in Fig. 2. It can be observed that when N increases, the expected value of the maximum E-field increases. The confidence interval of maximum E-field can be calculated from CDF in (2) [1], [5].

images

Figure 2: PDF plots of the normalized maximum E-field (normalized to the mean value |Ex|) for different independent sample number N.

The expected value of the maximum-to-mean ratio of the rectangular E-field be expressed as [2], [5]:

α(N) =|Ex|N|Ex|
=0Nx2[1-exp(-x22)]N-1exp(-x22)dxπ/2. (3)

However, numerical integration is required to evaluate (3) using a computer, an approximate equation could be necessary to have a quick estimation for a given N. By applying the binomial theorem:

α(N)
= 0Nx2k=0N-1(N-1k)[-exp(-x22)]kexp(-x22)dxπ/2
= k=0N-1N(N-1k)(-1)k0x2[exp(-x22)]k+1dxπ/2
= k=0N-1N(N-1k)(-1)kπ(k+1)2k+2π/2
= k=0N-1N!(-1)k(N-k-1)!(k+1)!k+1
= k=0N-1s(k). (4)

Although (3) is converted to the summation of a finite series, (4) is still not easy to evaluate. Symbolic calculation is required to maintain high precision, as the magnitude of s(k) resonant drastically. A plot of |s(k)| is illustrated in Fig. 3, as can be seen, when N is large, the magnitude of the series varies from 100 to 1030. If there is no special treatment for s(k), numerical calculation is easy to diverge for large N and wrong results could be obtained.

images

Figure 3: Plots of s(k) for different N, symbolic calculation is used to keep precision.

It seems we cannot go further from (4). Instead, we can approach it physically. Suppose an E-field probe (the Rx antenna) is placed inside a well-stirred RC, it can be found that the maximum received power (max(PrRC)=PrRCN) can be expressed as [2]:

max(PrRC)max(|Ex|)ληRx8π5, (5)

where |Ex| represents the tangential component of the E-field in an RC (|Ex|2=4|Ex|2/π, max(|Ex|)=|Ex|N), λ is the wavelength, and ηRx is the efficiency (including mismatch loss) of the Rx antenna respectively. It is interesting to compare it with received power in an anechoic chamber (AC):

PrAC=EACληRx8π52D/3, (6)

where EAC is the magnitude of the incident E-field and D is the directivity of the Rx antenna in the direction of the incident wave. Not surprisingly, D plays an important role in an AC. it can be found that when D=3/2, (5) and (6) give the same results. This means that if the E-field probe is calibrated in an AC, only when the directivity of the E-field probe is 3/2 (electrically small dipole antenna), the measured E-field in an RC is statistically accurate. When the frequency is high, the E-field probe is no longer an electrically small antenna; the measured mean E-field in an RC using the free space antenna factor is no longer accurate and ηRx is actually the key parameter. This effect is well known in standards related to RC measurements [1]. At high frequencies, an antenna is normally used to monitor the E-fields inside an RC instead of using an E-field probe.

From statistical analysis in [2]:

max(PrRC)PrRCk=1N1k. (7)

Note that [2], [7]:

PrRC =12E02ηλ24πηRx=32|Ex|2ηλ24πηRx
=6|Ex|2πηλ24πηRx. (8)

From (5), (7) and (8) we have:

max(|Ex|)8π5ληRx6|Ex|2πηλ24πηRxk=1N1k. (9)

Thus for large N:

α(N)max(|Ex|)|Ex|4πk=1N1k. (10)

The problem is now converted to the partial sum of the harmonic series. The following equation can be used [8]:

k=1N1k=γ+ψ0(N+1), (11)

where γ0.5772 is the Euler-Mascheroni constant, ψ0(x) is the digamma function:

ψ0(N+1) lnN+12N-n=1B2n2nN2n
=lnN+12N-112N2+1120N4-1252N6+,

and B2n is the Bernoulli number x(ex-1)=n=0Bnxnn!. Thus we have:

α(N)4π(γ+lnN+12N-112N2+1120N4-1252N6+). (12)

Now we have obtained an analytical approximation for (3). In [4], the CDF value of 0.5 is used to calculate the expected value:

α(N)4πln1[1-0.51N]. (13)

In [1] and [5], a similar expression is given with α(N)4π[0.5772+ln(N+1)-12(N+1)]with slightly different parameters. Another approximation is given in [8] for the partial sum of harmonic series

k=1N1k(2N+1)tan-1(12N+1)
+12ln(2N2+2N+12)
+56N6+168N5+140N4-42N2-14N-12520(2N2+2N+1)7
-6N2+6N+1180(2N2+2N+1)3+γ-1. (14)

images

Figure 4: (a) Plots of (4), (12)-(14) for different N, (b) error plot of (12)-(14), units are in linear.

A comparison of α(N) is given in Fig. 4 (a) using (12)-(14), (12) is truncated with 3 terms, while the exact value is calculated using symbolic calculation in (4). The difference between (12)-(14) and the exact values (the error plots) are illustrated in Fig. 4 (b), respectively. As can be seen, (12) and (14) give smaller errors than (13). However, they all give good approximations when N is large.

For the approximate value of the relative standard deviation of |Ex|N, we can start from the definition

stdrel(|Ex|N) =0x2p(x)dx-π2α(N)2α(N)π2
= k=0N-1N(N-1K)(-1)k0x3[exp(-x22)]k+1dx-α(N)2α(N)π2
= k=0N-1N!(-1)k2(N-k-1)!(k+1)!(k+1)-π2α(N)2α(N)π2
= 2k=1N1k-π2α(N)2α(N)π2
= 4πk=1N1kα(N)2-1. (15)

It seems we cannot go further from (15). To obtain the approximate value of (15), a different approach can be used [5]. We start from the relative standard deviation of |Ex|2N. From [2], we have:

stdrel(|Ex|2N) =std(|Ex|2N)|Ex|2N
= k=1N1k2k=1N1k, (16)

and we know that the uncertainty of power is 2 times the uncertainty of E-field (from P|Ex|2we have dP/P2d|Ex|/|Ex|) [5], thus the relative standard deviation of |Ex|N can be approximated as:

stdrel(|Ex|N) 12stdrel(|Ex|2N)12π26-1Nγ+ψ0(N+1)
=12π26-1Nγ+lnN+12N-112N2+, (17)
π261γ+lnN. (18)

In [5], a similar expression is given with 12π26-N+1N(N+2)0.5772+ln(N+1)-12(N+1).

The plot of (15),(17) and (18) are illustrated in Fig. 5 (a), (15) is used as the reference values (the exact values), and the difference between (17), (18) and (15) are presented in Fig. 5 (b), respectively. It is interesting to note that the convergence speed of (18) is different from the mean value (by using the central limit theorem) which is 1N [5].

images

Figure 5: (a) Plots of the approximated relative standard deviation in dB units (20log10(1+stdrel)), (b) error plot of (17), (18), units are in dB.

images

Figure 6: Comparison plots of exact PDFs and GEV PDFs for different N.

We can also link the parameters to the generalized extreme value distribution (GEV) [10]–[18] when N. By applying [10]:

s=6πstd=6πstdrel(|Ex|N)α(N), (19)
m=α(N)-γs. (20)

The approximated GEV CDF is:

GGEV(x)=exp[-exp(-x-ms)], (21)

Table 1: A summary of maximum-to-mean ratio expressions

the PDF is:

pGEV(x)=1sexp[m-x-sexp(-x-ms)s] (22)

When N, (19) and (20) can be approximated as:

s=6πstd1πγ+lnN+12N, (23)
m4π(γ+lnN+12N)-γπγ+lnN+12N. (24)

The plot for the exact PDFs (1) and GEV PDFs (22) of the normalized maximum rectangular E-field are given in Fig. 6. As expected, when N, two PDFs are close to each other.

For the PDF of the maximum-to-mean ratio of the received power |Ex|2N, we have the same pGEV(x) when N. The only difference is: the standard deviation and the expected value (α(N)) in (19) and (20) are replaced by the values of k=1Nk-2 and k=1Nk-1, respectively.

3 CONCLUSIONS

We have reviewed different analytical equations for the expected value and the standard deviation of the normalized maximum rectangular E-field in a well-stirred RC. Useful derivations are revisited and similar results are obtained. Numerical results show that: for small independent sample number N, we can calculate the results using the finite series directly; when N>5, very good approximations can be obtained for both the expected value and the standard deviation. Table 1 summarizes the expressions in different forms. The approximated analytical expressions are also linked to the parameters of the GEV parameters when N.

ACKNOWLEDGMENT

This work was supported in part by the Fund of Prospective Layout of Scientific Research for NUAA (Nanjing University of Aeronautics and Astronautics).

REFERENCES

[1] IEC 61000-4-21, Electromagnetic compatibility (EMC) – Part 4-21: Testing and measurement techniques – Reverberation chamber test methods, IEC Standard, Ed 2.0, 2011-01.

[2] J. Ladbury, G. Koepke, and D. Camell, Evaluation of the NASA Langley Research Center mode-stirred chamber facility, National Institute of Standards and Technology, Gaithersburg, MD, USA, Tech. Note 1508, Jan. 1999.

[3] T. H. Lehman and G. J. Freyer, “Characterization of the maximum test level in a reverberation chamber,” IEEE International Symposium on Electromagnetic Compatibility, pp. 44-47, 1997.

[4] L. R. Arnuat and P. D. West, Electric Field Probe Measurements in the NPL Untuned Stadium Reverberation Chamber, NPL CETM13, Sep. 1999.

[5] L. R. Arnaut, Measurement Uncertainty in Reverberation Chambers – I. Sample Statistics, NPL Report TQE 2, Dec. 2008.

[6] Q. Xu and Y. Huang, Anechoic and Reverberation Chambers: Theory, Design and Measurements, Wiley-IEEE, 2019.

[7] D. A. Hill, Electromagnetic Fields in Cavities: Deterministic and Statistical Theories, Hoboken, NJ, USA, Wiley, 2009.

[8] E. W. Weisstein, Harmonic Series, From MathWorld – A Wolfram Web Resource. https://mathworld.wolfram.com/HarmonicSeries.Html

[9] T. N. Vu, A New Approximation Formula for Computing the N-th Harmonic Number (Update), from Series Math Study Resource. http://www.seriesmathstudy.com/sms/ApproximateFormulaHarmonicNum, 11/11/2011.

[10] G. Gradoni and L. R. Arnaut, “Generalized extreme-value distributions of power near a boundary inside electromagnetic reverberation chambers,” IEEE Transactions on Electromagnetic Compatibility, vol. 52, no. 3, pp. 506-515, Aug. 2010.

[11] G. Orjubin, “Maximum field inside a reverberation chamber modeled by the generalized extreme value distribution,” IEEE Transactions on Electromagnetic Compatibility, vol. 49, no. 1, pp. 104-113, Feb. 2007.

[12] C. E. Hager and J. D. Rison, “On applying the generalized extreme value distribution to undermoded cavities within a reverberation chamber,” IEEE International Symposium on Electromagnetic Compatibility & Signal/Power Integrity (EMCSI), pp. 681-686, 2017.

[13] N. Nourshamsi, J. C. West, and C. F. Bunting, “Effects of aperture dimension on maximum field level inside nested chambers by applying the generalized extreme value distribution,” IEEE Symposium on Electromagnetic Compatibility, Signal Integrity and Power Integrity (EMC, SI & PI), pp. 628-633, 2018.

[14] N. Nourshamsi, J. C. West, and C. F. Bunting, “Investigation of electromagnetic complex cavities by applying the generalized extreme value distribution,” IEEE Symposium on Electromagnetic Compatibility, Signal Integrity and Power Integrity (EMC, SI & PI), pp. 622-627,2018.

[15] N. Nourshamsi, J. C. West, C. E. Hager, and C. F. Bunting, “Generalized extreme value distributions of fields in nested electromagnetic cavities,” IEEE Transactions on Electromagnetic Compatibility, vol. 61, no. 4, pp. 1337-1344, Aug.2019.

[16] P. Hu, Z. Zhou, X. Zhou, J. Li, J. Ji, M. Sheng, and P. Li, “Generalized extreme value distribution based framework for shielding effectiveness evaluation of undermoded enclosures,” International Symposium on Electromagnetic Compatibility-EMC EUROPE, pp. 1-6,2020.

[17] P. Hu, X. Zhou, and Z. Zhou, “On the modelling of maximum field distribution within reverberation chamber using the generalized extreme value theory,” IEEE MTT-S International Conference on Numerical Electromagnetic and Multiphysics Modeling and Optimization (NEMO), pp. 1-4, 2020.

[18] P. Hu, Z. Zhou, X. Zhou, and M. Sheng, “Maximum field strength within reverberation chamber: A comparison study,” 6th Global Electromagnetic Compatibility Conference (GEMCCON), pp. 1-4, 2020.

BIOGRAPHIES

images

Qian Xub (Member, IEEE) received the B.Eng. and M.Eng. degrees from the Department of Electronics and Information, Northwestern Polytechnical University, Xi’an, China, in 2007 and 2010, and received the PhD degree in electrical engineering from the University of Liverpool, U.K, in 2016. He is currently an Associate Professor at the College of Electronic and Information Engineering, Nanjing University of Aeronautics and Astronautics, China.

He was as a RF engineer in Nanjing, China in 2011, an Application Engineer at CST Company, Shanghai, China in 2012. His work at University of Liverpool was sponsored by Rainford EMC Systems Ltd. (now part of Microwave Vision Group) and Centre for Global Eco-Innovation. He has designed many chambers for the industry and has authored the book Anechoic and Reverberation Chambers: Theory, Design, and Measurements(Wiley-IEEE, 2019). His research interests include statistical electromagnetics, reverberation chamber, EMC and over-the-air testing.

images

Rui Jia received the B.Eng. degrees from the Department of Electronics and Information, Zhengzhou University, Zhengzhou, China, in 2008, and received the M.Eng. and Ph.D. degrees in electronic engineering from Mechanical Engineering College, Shijiazhuang, China, in 2010 and 2014, respectively. His research interests include statistical electromagnetics, reverberation chamber, and EMC.

images

Lifei Geng received the B.Eng., M.Eng. and Ph.D. degrees from the Department of Electronics and Information, Mechanical Engineering College, Shijiazhuang, China, in 2006, 2008 and 2013, respectively. His research interests include electromagnetic theory, reverberation chamber, and EMC.

images

Hao Guo received the B.Sc. and M.Sc. degrees from the Department of Physics, Dalian University of Technology, Dalian, China, in 2004, and 2009, respectively. His research interests include plasma, reverberation chamber, and EMC.

images

Yongjiu Zhao received the M.Eng. and Ph.D. degrees in electronic engineering from Xidian University, Xi’an, China, in 1990 and 1998, respectively. Since March 1990, he has been with the Department of Mechano-Electronic Engineering, Xidian University where he was a professor in 2004. From December 1999 to August 2000, he was a Research Associate with the Department of Electronic Engineering, The Chinese University of Hong Kong. His research interests include antenna design, microwave filter design and electromagnetic theory.

Abstract

1 INTRODUCTION

images

2 THEORY

images

images

images

images

images

3 CONCLUSIONS

ACKNOWLEDGMENT

REFERENCES

BIOGRAPHIES