A novel methodology for time-domain characterization of a full anechoic chamber for antennas measurements and exposure evaluation

1. International Journal of Electrical and Computer Engineering (IJECE) Vol. 11, No. 4, August 2021, pp. 3285∼3292 ISSN: 2088-8708, DOI: 10.11591/ijece.v11i4.pp3285-3292 r 3285 A novel methodology for time-domain characterization of a full anechoic chamber for antennas measurements and exposure evaluation Chakib Taybi, Mohammed Anisse Moutaouekkil, Bachir Elmagroud, Abdelhak Ziyyat Electronic and Systems Laboratory, Faculty of Sciences, Mohammed First University, Oujda, Morocco Article Info Article history: Received Aug 11, 2020 Revised Dec 15, 2020 Accepted Jan 19, 2021 Keywords: Anechoic chamber Antennas measurement Autoregressive filter Electromagnetics measurement Exposure evaluation Finite technique method Time domain simulation ABSTRACT In this paper we present a novel methodology for time-domain characterization of a full anechoic chamber using the finite integral method. This approach is considered fast, accurate and not intensive for computer resources. The validation of this approach is carried out on CST-microwave studio for a full anechoic chamber intended for antennas measurement applications and electromagnetic exposure evaluation for cellular network. Low, medium and high gain sources are used in this study. The simulations are realized on a personal computer of medium performances (i7 CPU and 16 GB of RAM). The stability and the convergence of our approach are obtained thanks to local mesh and auto-regressive linear filtering techniques. The minimization of the simulation time is based on use of the Huygens sources in the place of the antennas. The maximum error of the chamber as well as the wave depolarization into the chamber are at one with the previous work and the catalogs of the principles chambers manufacturers for the proposed tests in this paper. The Full simulations time is about 15 hours in average. This is an open access article under the CC BY-SA license. Corresponding Author: Chakib Taybi Electronic and Systems Laboratory Faculty of Sciences Mohammed First University P.O. Box: 717, Oujda, Morocco Email: [email protected] 1. INTRODUCTION Today, anechoic chamber becomes a perfect solution for antenna measurements [1] and for the as- sessment of exposition to electromagnetic fields (EMF) [2]. They allow to reproduce OATS conditions [3] in a limited, closed and protected environment. Different forms and types of anechoic chamber are presented in previous work [4, 5]. However, the rectangular form is the most one used for antennas measurement and EMF evaluation. Typically, an anechoic chamber essentially consists of a shielding structure [6], protects the volume test, called quiet zone generally, from the noise sources and a matrix of electromagnetic absorber [7] minimizes the reflections on the walls of the chamber. Special antennas type, called probes [8], are used to sample the EM field radiated by the antenna under test. The characterization of the anechoic chambers and its constituent dates back to the 1975 [4]. In this era, the Ray-Tracing [9, 10] method was used to model objects whose dimensions are higher than the Journal homepage: http://ijece.iaescore.com

2. 3286 r ISSN: 2088-8708 wavelength. Years later, other methods benefit, from the computing and storage power of personal computer, and from the development of numerical codes, which are implemented for the characterization of chambers whose dimensions are in the lower of wavelength. In this context, the finite elements method (FEM) presented in [11, 12] was used for modeling a full anechoic chamber for antenna measurements. The simulations were launched on a super calculation station (dual Xeon CPU and 256 GB of RAM) and they required 196 GB of memory and a calculation time of 15 hours for a frequency of 1 GHz. For higher frequencies the proposed method requires more memory and more of the computing time. The finite differences in time domain FDTD were also used to analyze the behavior of a full anechoic and compact chambers [13, 14]. The obtained results show that this method is more suitable for dielectric objects such as electromagnetic absorbers. For metallic structures such as walls, probes, antennas and scanners system, it requires a very fine mesh. Both a larger memory space and a longer computation time are needed for this purpose. An alternative solution based on the transmission lines matrix (TLM) was presented by Gwenael [15] in the order to characterize a full and a semi-anechoic chamber. The obtained results show the need for a special type of perfect matched layers (PML) [15] to reduce the IT resource requirements and also the simulation time. Therefore, the simulations were lunched out on a personal computer and lasted about 150 hours in average. In addition, the work presented in [16, 17] demonstrates the need to combining TLM with an other method to ensure accuracy at lower frequencies. Also, vision-based technology are used for modeling many physical systems [18]. These methods proved their utility for wireless channel characterization [19], vehicle problem and RFID signal measurements [20]. However, for electromagnetic problems, these methods reverent incapable to describe Maxwell equations in a finite cell and the boundary conditions for a large structures such as anechoic chamber. In order to optimize IT resource requirements and to properly control the simulation time, we have presented in [21] a novel approach for the characterization of a full anechoic chamber. Based on the finite integrals method (FIT) and on linear regression filtering technique. This novel method results show that it’s accurate, fast and not intensive in computing resources. The validation was performed for a chamber for antenna measurements and with a medium gain source. In this paper, we are expanding the validation domain of our novel approach. The chamber behavior was evaluated for a low, medium and high gain sources. This three sources types, represent the possible cases encountered in an anechoic chamber for antenna measurements. The chamber error and the wave depolarization in this later was evaluated as function of the three proposed sources. The calculation time and the memory size used are also compared in this paper. 2. METHOD AND MATERIAL 2.1. Method The characterization of the structure of the full anechoic chamber and its elements was carried out with the finite integrals technique (FIT) [22, 23] available under the CST microwave studio suite. This technique allows Maxwell equations resolution in the integral form described by (1-4). I ∂A E.dS = − Z A ∂B ∂t .dA (1) I ∂A H.dS = − Z A (J + ∂D ∂t ).dA (2) I ∂V B.dA = 0 (3) I ∂V D.dA = Z V ρ.dV (4) A cubic elementary mesh as shown in Figure 1 discretization is applied to the structure of the chamber or to one of its elements. Equation (1) can be written on a mesh in the form of (5). êl + êk + êj + êi = − d dt b̂n (5) Int J Elec & Comp Eng, Vol. 11, No. 4, August 2021 : 3285 – 3292

3. Int J Elec & Comp Eng ISSN: 2088-8708 r 3287 The scalars ê and b̂n represent respectively the E field and the B flux on the edge and the face of a mesh cell. Equation (5) can also be written as (6). By similar demonstration, the (2), (3) and (4) can be written in the form of the (7), (8) and (9). C.ê = − d dt .b̂n (6) C̄.ĥ = ĵ + d dt .b̂ (7) S.b̂ = 0 (8) S̄. ˆ d = ρ (9) Where C and S are topological matrix (-1,0,1) and represent the edge-to-face incidence. Time domain resolution (or transit solver in CST) of the previous (6), (7), (8), (9) allows the field e and the flux b to be written in an iterative form Figure 2 given by (10), (11). hi+1 = hi − ∆t.M−1 µ .C.ei+ 1 2 (10) ei+ 3 2 = ei+ 1 2 + ∆t.M−1 [C̄.hi+1 − ji+1 ] (11) Where M and Mµ represent the discrete matrix of permittivity and permeability. Figure 1. Cubic elementary mesh discretization Figure 2. Iterative time domain resolution of Maxwell equation For a large and highly resonant structure (case of an anechoic chamber), the FIT becomes unstable and more intensive both in computing resources and in computing time due to the very fine mesh requirements. Calculation time itself is a function of the number of computation cells, of the process of excitation of the probe and the AUT. To overcome this problem, reduce and minimize both the IT resources and also the computation time, the EM field of the source is calculated on another simulation file and with an adequate mesh. These EM fields are calculated separately and are exported into the main simulation file as a Huygens sources. This technique reduce significantly the simulation as will discussed in the results part. However, the instability and divergence of FIT remains one of the main problems dealt with in this method. This is mainly due to the high resonance of the structure and to the FIT conditioning problems. A novel methodology for time-domain characterization of a full anechoic chamber for... (Chakib Taybi)

4. 3288 r ISSN: 2088-8708 To solve this problem, we have used auto-regressive (order 3) filter in parallel [24, 25]. This filter types are capable of calculating spectra (12) for time signals which have not yet reached their equilibrium state. These time signals are essentially those of the EM field monitors. x(n) = − p X k=1 akx(n − k) (12) Where p is the filter order. 2.2. Simulation materials In order to validate our approach, we studied in this paper the behavior of a full rectangular anechoic chamber for near and/or far field antenna measurement applications. The dimensions of this chamber are 500×300×300 (L(cm)-W(cm)-H(cm)) and it consists as shows the Figure 3 of the following elements: − Shielding structure: Modeled by a metal cage of conductivity 5.8.108 S.m−1 and thickness of 0.5 cm. This cage has five apertures associated to their waveguides for the passage of cables. The shielding efficiency (SE) of this cage was studied in detail in our previous article [21]. − Electromagnetic absorbers: Two types of absorbers were used in this study. The first one is flat absorbers (AT), used in the corners of the chamber. The second one is pyramidal shape (APM21, APM30 and APM 45), used in high incidence areas. The dimensions and the reflectivity of the absorbents used are reported in [21]. − Antennas sources: Three types of antennas are tested in this paper. Half-wave dipole, open ended rectangular wave guide (OERWG) and finally a pyramidal horn. These types of sources are the most used in an anechoic chamber for antennas measurement applications. The dimensions and the radiation pattern of the antenna sources at a frequency of 1 GHz are reported in literature and our previous works [5, 21]. (a) (b) Figure 3. Illustration of the full chamber form and measurements planes, (a) dimensions and the components of the chamber, (b) Measurement planes inside the chamber Because our method is fast and doesn’t require a lot of computer resources, the simulations are launched on a personal computer of 3.4 GHz of CPU and 16 GB of RAM. A recent version of CST microwave studio is installed on this computer. The performances of the anechoic chamber were evaluated using two principles parameter. The first one is the chamber error (13), defined by the ratio of the field radiated by the source inside the chamber and the field radiated in free space (without chamber). The second one is the axial ratio, represents the wave depolarization (14) inside the chamber. This last one is defined by the ratio between the vertical to the horizontal component of the field radiated inside the chamber. These two parameters (error and axial ratio) are generally given in dB value (0 dB is assigned to a perfect chamber). Practically, a value lower than 4 dB is generally supported by antennas measurement systems. Err (dB) =

9. s P |Fieldchm|2 P |Fieldfree|2

14. (13) Int J Elec Comp Eng, Vol. 11, No. 4, August 2021 : 3285 – 3292

15. Int J Elec Comp Eng ISSN: 2088-8708 r 3289 AR (dB) =

20. sP |FieldH|2 P |FieldV |2

25. (14) The components of the field inside the chamber and in the free space are calculated for a frequency of 1 GHz and on two orthogonal planes (Plane E and Plane H) localized at center of the chamber as shown the Figure 3. 3. RESULTS AND DISCUSSION In this section, we discuss the results validating our novel approach. As mentioned in the previous section, the first step in our approach consists in calculating the Huygens sources necessary to characterize the EM field in the chamber. The Table 1 summarizes for each source the directivity D in dBi, the number of cells mesh, the size of the memory and the simulation times. Table 1. Resume of Huygens calculation for each source Source D Mesh RAM Time Dipole 2.5 2560 ≤ 1GB ≤ 1min OERWG 6.2 15450 1.3 GB 5 min Horn 15 265400 3.2 GB 15 min To validate our approach and to analyze the behavior of the chamber, we present on the Figures 4-7, for each of the sources, the error (13) and the axial ratio (14), in the full length of the chamber and also in its quiet zone. The chamber error and the axial ratio are analyzed on both planes and also as function of the source polarization. For each of the sources studied in this paper, the obtained results show that the error is maximum in the structure of the absorbers, in zones with great reflections and in areas near to the radiating source. This error, which varies as function of the source and the previous regions, is about of 50 dB for both horizontal and vertical polarizations. In the quiet zone of the chamber, the maximum error is 5 dB, 1 dB and 3 dB for dipole, OERWG and horn source respectively for the two polarizations. In the same way, the axial ratio, which represents the depolarization of the wave in the chamber is maximum in the absorbers structures and near to the source and it is minimum in the quiet zone in the chamber. For the dipole source, it is by means of 3 dB, 1.2 dB for the OERWG and 2 dB for the horn, for horizontal as for vertical polarization. Analysis of the results presented in the figures below shows that the dimensions of the chamber proposed in this study are compatible with the sources with moderate directivity. For other types of sources, a calibration process is necessary to correct the error and the depolarization of the wave in the chamber. On the other hand, the chamber error and the axial ratio calculated with our method are comparable with the results published in the bibliography [11, 12, 26] and in the catalogs of the chamber manufacturers [27, 28]. Unlike other methods, which require either good hardware or enormous computation time, our simulations are performed on a desktop computer with acceptable performances. The Table 2 summarizes for each source the total number of cells, the memory required and finally the total calculation time. Table 2. Resume simulations details of each sources Source Mesh RAM Time Dipole 58000000 12.5 GB 12.5 h OERWG 61000000 14.8 GB 14.4 h Horn 64000000 15.4 GB 5.5 h The results presented in the Table 2 show the capacity and the solidity of our new approach for the characterization of a full anechoic chamber on a machine with acceptable performances, in a reasonable time and with an error and an axial ratio tolerable by the most of the antenna measurement systems. The slight difference of the number of mesh cells, memory and of the computation time from source to another is mainly due to the spherical surface used to calculate the Huygens sources. A novel methodology for time-domain characterization of a full anechoic chamber for... (Chakib Taybi)

26. 3290 r ISSN: 2088-8708 Err(PlanE,Chm,H) Err(PlanH,Chm,H) Err(PlanE,Chm,V) Err(PlanH,Chm,V) Err(PlanE,Zone,H) Err(PlanH,Zone,H) Err(PlanE,Zone,V) Err(PlanH,Zone,V) Figure 4. Error in the full length of chamber and in the quiet zone at 1GHz for a dipole source Err(PlanE,Chm,H) Err(PlanH,Chm,H) Err(PlanE,Chm,V) Err(PlanH,Chm,V) Err(PlanE,Zone,H) Err(PlanH,Zone,H) Err(PlanE,Zone,V) Err(PlanH,Zone,V) Figure 5. Error in the full length of chamber and in the quiet zone at 1 GHz for a OERWG source Err(PlanE,Chm,H) Err(PlanH,Chm,H) Err(PlanE,Chm,V) Err(PlanH,Chm,V) Err(PlanE,Zone,H) Err(PlanH,Zone,H) Err(PlanE,Zone,V) Err(PlanH,Zone,V) Figure 6. Error in the full length of chamber and in the quiet zone at 1 GHz for a dipole source Int J Elec Comp Eng, Vol. 11, No. 4, August 2021 : 3285 – 3292

27. Int J Elec Comp Eng ISSN: 2088-8708 r 3291 AR(PlanE,Chm,Dipole) AR(PlanE,Zone,Dipole) AR(PlanH,Chm,Dipole) AR(PlanH,Zone,Dipole) AR(PlanE,Chm,OERWG) AR(PlanE,Zone,OERWG) AR(PlanH,Chm,OERWG) AR(PlanH,Zone,OERWG) AR(PlanE,Chm,Horn) AR(PlanE,Zone,Horn) AR(PlanH,Chm,Horn) AR(PlanH,Zone,Horn) Figure 7. Axial ratio in the full length of chamber and in the quiet zone at 1 GHz for a dipole, OERWG and Horn source 4. CONCLUSION In this paper, we have presented a novel methodology for the modeling and characterization of a full anechoic chamber for antenna measurements. This approach is based on the time domain resolution of Maxwell’s equations with the method of finite integrals technique. The reduction in computing time and IT resources was achieved by replacing the traditional process of excitation of the antennas by direct excitation of the structure of the chamber by the equivalent Huygens sources of these antennas. This novel technique allowed us to reduce the total simulation time and also the need for IT resources by 90%. The stability and convergence of our approach were ensured by an autoregressive filter. The results obtained show that the chamber error calculated by the ratio between the field in free space and the field in the chamber as well as the axial ratio given by the ratio of vertical components to horizontal components are comparable with those published in the literature and catalogs of anechoic chamber manufacturers. In order to prove the validity of our approach, the behavior of the chamber will be studied, in the low frequencies where the absorbers are of large size and in the high frequency where the structure requires a very fine mesh. In addition, the chamber performances will be investigated as function of the geometric form of the absorbers. Others anechoic chambers forms, as semi-anechoic and compact chambers are under test to proven our novel approach. REFERENCES [1] Taybi, C., Moutaouekkil, M. A., Rodrigues, K. K., Elmagroud, B., and Ziyyat, A, “Probes correction for planar near field antennas measurements,” in 2015 IEEE 15th Mediterranean Microwave Symposium (MMS), jan 2015, pp. 1-4. A novel methodology for time-domain characterization of a full anechoic chamber for... (Chakib Taybi)

28. 3292 r ISSN: 2088-8708 [2] Kwate, R. K., Derkaoui, A., Elmagroud, B., Taybi, C., Ziyyat, A., and Picard, D, “Non-invasive estimation of whole body averaged SAR by the personal RF dosimeter measurements,” in 2019 7th Mediterranean Congress of Telecom- munications (CMT), Fès, Morocco, 2019, pp. 1-4. [3] Meng, D., Li, X., and Gao, X, “Near-field corrections to antenna factor of dipole antennas calibrated in very narrow OATS,” in International Conference on Microwave and Millimeter Wave Technology, ICMMT 2010, Chengdu, China, 2010, pp. 1869-1872. [4] L. H. Hemming, “Electromagnetic Anechoic Chambers: A Fundamental Design and Specification Guide,” Wiley Interscience, 2002. [5] C. A. Balanis, “Antenna Theory: Analysis and Design,” 3rd ed. Wiley, 2012. [6] S. Celozzi, R. Araneo, and G. Lovat, “Electromagnetic Shielding. Wiley,” 2008. [7] V. Rodriguez, “Validation of a Method for Predicting Anechoic Chamber Performance: A Technique That Uses Polynomial Approximations for RF Absorber Reflectivity,” IEEE Antennas and Propagation Magazine, vol. 60, no. 4, pp. 31–40, aug 2018. [8] Taybi, C., Derkaoui, A., Khalid, S., Moutaouekkil, M. A., Elmagroud, B., and Ziyyat, A, “Geometrical and Ana- lytical Methods for Probes Correction applied to Near Field Antennas Measurement Systems,” in 7th Mediterranean Congress of Telecommunications 2019, CMT 2019, Fès, Morocco, 2019, pp. 1-5. [9] M. Gillette, “RF anechoic chamber design using ray tracing,” in Antennas and Propagation Society International Symposium, Stanford, CA, USA, jun 1977, pp. 246–249. [10] M. K. Mansour and J. Jarem, “Anechoic chamber design using ray tracing and theory of images,” in Conference Proceedings-IEEE SOUTHEASTCON, vol. 2, New Orleans, LA, USA, 1990, pp. 689–695. [11] D. Campbell, G. Gampala, M. Vogel, and C. Reddy, “Simulating antenna measurements in an anechoic chamber,” Microwave Journal, vol. 57, no. 7, pp. 90–100, 2014. [12] D. Campbell and Gampal, “Modeling and analysis of anechoic chamber using CEM tools,” Applied Computational Electromagnetics Society Journal, vol. 28, no. 9, pp. 755–762, 2013. [13] Kawabata, M., Ishida, Y., Shimada, K., and Kuwabara, N, “FDTD Method for Site Attenuation Analysis of Compact Anechoic Chamber Using Large- Cell Concept,” IEEE Transactions on Fundamentals and Materials, vol. 125, no. 10, pp. 784–790, 2005. [14] Gradoni, G., Moglie, F., Pastore, A. P., and Primiani, V. M, “Numerical and experimental analysis of the field to enclosure coupling in reverberation chamber and comparison with anechoic chamber,” IEEE Transactions on Elec- tromagnetic Compatibility, vol. 48, pp. 203–211, 2006. [15] G. Dun, “Modélisation et optimisation de chambres anécho ıques pour applications CEM,” Ph.D. disser- tation, Bretagne Occidentale, 2007. [16] F. Bellamine, “Analysis of semi-anechoic chambers using the frequency dependent transmission line modeling,” in European Space Agency, vol. 626 SP, 2006. [17] B. F, “Simulation of anechoic chambers using a combined TLM and homogenisation method in the frequency range 30-200 MHz,” in First European Conference on Antennas and Propagation, nov 2008, pp. 1–6. [18] R. Maleewan and S. Thitirat, “Deep Learning Based Gesture Classification for Hand Physical Therapy Interactive Program,” Lecture Notes in Computer Science, vol. 12198, pp. 349–358, 2020. [19] X. Tianqi, “A Computer Vision Based Beamforming Scheme for Millimeter Wave Communication in LOS Scenar- ios,” in 7th International Conference on Computer Science and Network Technology, Oct 2019, pp. 1–6. [20] Xie, L., Wang, C., Bu, Y., Sun, J., Cai, Q., Wu, J., and Lu, S, “An RFID-Based Approach for Recognition of Multiple Tagged Objects in Augmented Reality Systems,” IEEE Transactions on Mobile Computing, vol. 15, pp. 1188–1202, 2020. [21] Taybi, C., Ziyyat, A., Moutaouekkil, M. A., and Elmagroud, B, “Novel Approach for Time-Domain Modeling a Full Anechoic Chamber for NearIFar-Field Antennas Measurement,” in International Symposium on Advanced Electrical and Communication Technologies, ISAECT 2018, Rabat, Morocco, 2019, pp. 1-6. [22] M. Clemens and T. Weiland, “Discrete electromagnetism with the finite integration technique,” Progress in Electro- magnetics Research, vol. 32, pp. 65–87, 2001. [23] I. Munteanu and T. Weiland, “RF Microwave Simulation with the Finite Integration Technique-From Component to System Design,” in Scientific Computing in Electrical Engineering. Springer, Berlin, Heidelberg, 2007, pp. 247–260. [24] T. Wittig, R. Schuhmann, and T. Weiland, “Model order reduction for large systems in computational electromagnetics,” Linear Algebra and Its Applications, vol. 415, no. 2-3, pp. 499–530, 2006. [25] M. Mocker, S. Hipp, F. Spinnler, H. Tazi, and T. Eibert, “Comparison of electromagnetic solvers for antennas mounted on vehicles,” Advances in Radio Science, vol. 13, pp. 49–55, 2015. [26] A. Smirnov, “Procedural aspects of the certification of anechoic chambers,” Measurement Techniques, vol. 53, no. 5, pp. 556–562, 2010. [27] TDK, “Fully Compliant Anechoic Chambe,” Tech. Rep., 2016. [28] ETS-Lindgren, “ETS-Lindgren Anechoic Chambers,” Tech. Rep., 2016. Int J Elec Comp Eng, Vol. 11, No. 4, August 2021 : 3285 – 3292

A novel methodology for time-domain characterization of a full anechoic chamber for antennas measurements and exposure evaluation

More Related Content

What's hot (13)

Similar to A novel methodology for time-domain characterization of a full anechoic chamber for antennas measurements and exposure evaluation (20)

More from IJECEIAES (20)

Recently uploaded (20)

A novel methodology for time-domain characterization of a full anechoic chamber for antennas measurements and exposure evaluation