Articles | Volume 18
10 Dec 2020
 | 10 Dec 2020

Characteristic Mode Analysis of surface current distributions on metallic structures exposed to HIRF- and DCI-excitations

Jan Ückerseifer and Frank Gronwald

This paper treats Characteristic Mode Analyses of three-dimensional test objects in the context of EMC. Based on computed Characteristic Modes and mode-specific physical quantities, series expansions for HIRF- and DCI-induced surface currents are deduced. The contribution of single Characteristic Modes to surface currents at different test frequencies is analyzed. HIRF- and DCI-excitations are compared with regard to their surface current distributions in their resonance region determined by Characteristic Mode Analysis.

1 Introduction

Radiated immunity tests in anechoic chambers using HIRF (High Intensity Radiated Field) are a common EMC test procedure in the automotive sector described in DIN EN 61000-4-3 (2011). Related test environments include TEM waveguides specified by DIN EN 61000-4-20 (2010) and reverberation chambers, see DIN EN 61000-4-21 (2011), all of which are based on radiation coupling to a device under test (DUT).

In the recent past, a conducted alternative called DCI (Direct Current Injection), also refered to as Direct Drive in ED-107A (2015) found application in immunity testing. With the help of adapters galvanically coupled to a DUT, currents are directly injected onto its surface, see Leat (2007) for an exemplary test scenario. In this manner, DCI seeks to approximate the surface currents on DUTs generated during HIRF tests. As mathematical background serves the surface equivalence principle proving that identical charge or current distributions on an object's surface cause equal fields within its interior volume, see Harrington (2001). In general, EMC testing is most critical when a DUT is excited at its resonance frequencies. Characteristic Mode Analysis (CMA) is a tool that allows to determine the resonances of a system, characterized by so called Characteristic Modes. Since both HIRF and DCI test setups differ from each other geometrically, both sets of Characteristic Modes cannot coincide with each other, thus making DCI only an approximation to HIRF.

Other than most publications in which CMA has been applied to antenna design, see e.g. Vogel et al. (2015), EMC testing is addressed in this paper. Herefore, explicit mathematical expressions for surface currents of HIRF and DCI setups expanded in Characteristic Modes are deduced. In doing so, it is quantitatively investigated in how far Characteristic Modes determine surface current distributions in EMC tests at resonance and beyond.

The following sections are introduced by a summary on the theoretical background of CMA. It is subsequently applied on various test objects in HIRF- and DCI-configuration, starting off with a cylinder as generic test object. Finally, as more realistic test object, a car geometry is investigated.

2 Characteristic Mode Analysis

The mathematical foundations of Characteristic Mode Analysis are covered in Harrington and Mautz (1971). The basic problem to solve in CMA is the generalized eigenvalue problem

(1) X J n = λ n R J n

with Jn,nN signifying real-valued eigenfunctions called Characteristic Modes to determine, as are the eigenvalues λn[-,]. Matrices R and X are real and imaginary parts of the system's impedance matrix Z¯=R+jX.All Characteristic Modes Jn are orthonormal to one another with reference to the resistance matrix R,

(2) J m , R J n = δ m n = 1 for m = n 0 for m n ,

and can thus serve as basis for a Fourier series expansion of the total surface current J¯(r) generated in HIRF and DCI tests according to

(3) J ¯ ( r ) = n = 1 α ¯ n J n ( r ) n = 1 N α ¯ n J n ( r ) ,

where the α¯n denote Fourier coefficients and r=[x,y,z] the position on a DUT's surface. Because in practical computations only a finite number N of modes can be considered, the Fourier series can only be approximated.

Figure 1HIRF test setup.


Figure 2Characteristic Mode |J1H(r)| for f=fres,1.


As quantity to describe contributions of each mode to the total surface current dependent on the frequency f, the modal significance

(4) s n ( f ) = 1 1 + j λ n ( f ) = 1 1 + λ n 2 ( f )

is introduced with sn[0,1]. Especially important for EMC applications is the fact that CMA yields resonance frequencies fres,n of a system, which are indicated by values sn(fres,n)=1. In contrast to the system-inherent modal significance, the modal excitation coefficients

(5) V ¯ n i = J n , E ¯ i = S J n E ¯ i d S

indicate to which extent each mode is excited in dependency on the incident electrical field E¯i scattered at the DUT's surface S. It is thus no direct result of CMA, since the latter solely exploits information of the system itself.

Modal excitation coefficients V¯ni and Fourier coefficients α¯n correlate via

(6) α ¯ n = V ¯ n i 1 + j λ n

such that Eq. (3) can be rewritten as

(7) J ¯ ( r ) = n = 1 α ¯ n J n ( r ) n = 1 N V ¯ n i J n 1 + j λ n .

Note that CMA requires all investigated DUTs to be perfect electric conductors (PEC) causing the sum of incident and scattered tangential electrical field components to vanish on their surface.

Figure 3Characteristic Mode |J2H(r)| for f=fres,1.


Figure 4Characteristic Mode |J3H(r)| for f=fres,1.


Figure 5Characteristic Mode |J4H(r)| for f=fres,1.


3 Cylinder model

The first test setup in HIRF configuration consists of a PEC cylinder measuring l=5 m and r=1 m, situated in a height h=1 m above an infinitely extended ground plane, see Fig. 1. Hence, the conditions in a semi-anechoic chamber are tried to be approximated, which is the most common test site with regard to HIRF testing.

A linearly polarized plane wave characterized by an electrical field Ei at a given frequency f is used as excitation. Propagation vector k heads in 45-direction to the horizontal line, whereas magnetic field vector Hi is tangential to the ground plane.

Initially, a CMA of the surface current is carried out using the integral equation solver of CST Microwave Studio resulting in the Characteristic HIRF Modes JnH,n=1,,5 depicted in Figs. 2 to 6 as magnitude value. The triangular surface mesh used consists of nt=402 triangles and np=223 grid points.

Figure 6Characteristic Mode |J5H(r)| for f=fres,1.


Each mode corresponds to a certain resonance of the system such as fundamental mode J1H and second harmonic J2H, which indicate a λ∕2- and λ-resonance respectively.

The precise resonance frequencies of each mode are deduced from the modal significance sn(f) illustrated in Fig. 7. Herein, only the fundamental mode shows a typical resonance-like curve with fres,1≈20.7 MHz, whereas higher order modes either show no unique maximum or no maximum at all in the chosen frequency domain. One possible reason is a decreased quality factor due to radiation losses.

Figure 7Modal significance sn(f).


Figure 8 illustrates the magnitude of the modal excitation coefficients V¯ni(f) numerically determined using FEKO, since CST does not provide these quantities. As can be seen, only the first two modes exhibit a relevant excitation throughout the entire frequency range.

Figure 8Modal excitation coefficients |V¯ni(f)|.


After the cylinder in HIRF-configuration has been analyzed by CMA, an exemplary excitation at f=40 MHz with |Ei|=1 V m−1 according to the model in Fig. 1 is investigated. The frequency was chosen such that no single Characteristic Mode dominates over others.

Then, the HIRF current J¯HIRF(r) on the cylinder arranges like in Fig. 9. Apparently, the current density is highest in the two bottom parts to the right and left, accompanied by a minimum close to the center.

Figure 9HIRF current |J¯HIRF(r)| for f=40 MHz.


Taking the discussed results of CMA, it is possible to expand the total surface current of the preceeding HIRF excitation in terms of Characteristic Modes. Figure 10 outlines the associated Fourier series

(8) J ¯ HIRF r n = 1 N α ¯ n J n H r = 0.07 + j 0.28 J 1 H + 0.20 + j 0.24 J 2 H + 1.3 × 10 - 5 - j 2.9 × 10 - 5 J 3 H - 0.08 + j 0.10 J 4 H - 3.6 × 10 - 3 + j 1.9 × 10 - 4 J 5 H

after Eq. (3) taking into account only N=5 modes. Although a non-resonance frequency was used as excitation, only three Characteristic Modes have significant influence on the total current. Distribution as well as amplitude fit the HIRF current to an acceptable degree except for the top left cylinder area. Raising the number of modes up to N=15 could not eliminate this circumstance, however, all other maxima could be reproduced.

Figure 10Fourier series of |J¯HIRF(r)| for N=5 and f=40 MHz.


Since it is intended to approximate HIRF- by DCI-excitations, a comparison to a DCI setup's surface current is desirable. As test frequency serves the cylinder's first resonance obtained by CMA, which occurs at fres,1=20.7 MHz correspondent to Fig. 7. Like Fig. 11 points out, the HIRF current J¯HIRF(r) accumulates at the bottom with a slight deviation from the center to the right.

Figure 11HIRF current |J¯HIRF(r)| for f=fres,1.


Figure 12Fourier series of |J¯HIRF(r)| for N=5 and f=fres,1.


Figure 13DCI current |J¯DCI(r)| for f=fres,1.


The way the HIRF current is composed of single Characteristic Modes is revealed by the series expansion

(9) J ¯ HIRF r n = 1 N α ¯ n J n H r = 0.19 - j 0.47 J 1 H - 0.20 + j 0.10 J 2 H - 7.3 × 10 - 6 - j 3.1 × 10 - 6 J 3 H - 1.1 × 10 - 7 - j 7.6 × 10 - 6 J 4 H - 0.02 - j 5.7 × 10 - 4 J 5 H

showing a dominant fundamental mode, because its resonance frequency is used as excitation, but also a certain contribution of the second harmonic. Figure 12 stresses the aforementioned statement with a distribution looking similar to Characteristic Mode J1H in Fig. 2.

A DCI-configuration is specified by galvanically connecting both cylinder ends to the ground plane via two wire-shaped adapters. If the left one acts as excitation source with 50 Ω internal resistance and a power of 1 W at f=fres,1 and the right one as termination with a 1 MΩ resistance, the resulting DCI current J¯DCI(r) distributes like in Fig. 13.

Figure 14Fourier series of |J¯DCI(r)| for N=5 and f=fres,1.


Figure 15Absolute error ε(ri) of HIRF and DCI for f=fres,1.


Its decomposition into Characteristic DCI Modes JnD is outlined by the corresponding series expansion

(10) J ¯ DCI r n = 1 N α ¯ n J n D r = - 0.06 + j 0.16 J 1 D - 0.02 - j 0.11 J 2 D + 1.6 × 10 - 6 - j 2.3 × 10 - 5 J 3 D - 8.9 × 10 - 6 + j 1.9 × 10 - 4 J 4 D - 1.6 × 10 - 5 - j 6.4 × 10 - 4 J 5 D

depicted in Fig. 14 for the cylinder without its adapters. The attained coincidence in space and magnitude appears different. This is explainable with the fact that the DCI setup is not excited at its own resonance frequency, but with the HIRF setup's one. Consequently, the DCI fundamental mode is of decreased influence for the total DCI current.

A quantitative comparison of HIRF and DCI is possible, when both currents are subtracted from each other at corresponding triangles located at ri,i=1,,nt of the simulation mesh yielding their absolute error

(11) ε ( r i ) = J ¯ HIRF ( r i ) - J ¯ DCI ( r i ) ,

whose graphical representation is given in Fig. 15. Apart from a small residual current at the right cylinder face its approximation to HIRF seems satisfying.

Figure 16Characteristic Mode |J1H(r)| for f=fres,1.


Figure 17Characteristic Mode |J2H(r)| for f=fres,1.


Figure 18Characteristic Mode |J3H(r)| for f=fres,1.


Figure 19Characteristic Mode |J4H(r)| for f=fres,1.


4 Car model

As more complex three-dimensional object a PEC car body is investigated. Its length is approximately 4 m with a maximum height of 1 m. The HIRF setup is geometrically the same as presented in Fig. 1, only with the DUT having been exchanged.

Given these preconditions, a CMA is processed for the first five surface current modes JnH, whose magnitude is illustrated in Figs. 16 to 20. Similar to the cylinder, basic resonance patterns such as half-, full- and higher-wave resonances can be observed by analyzing maxima and minima of the surface current. The triangular surface mesh used for all following simulations consists of nt=12 412 triangles and np=6208 grid points.

Figure 20Characteristic Mode |J5H(r)| for f=fres,1.


Surprisingly, the modal significance sn calculated is comparable to those of the cylinder in Fig. 7, like Fig. 21 conveys. As an example, the first two harmonics basically have the same qualitative shape, but exhibit a frequency shift to slightly higher frequencies. This is due to the fact that the car is geometrically shorter implying higher resonance frequencies. Consequently, it can be stated that both test objects have similar electromagnetic properties regarding their frequency behavior, although being geometrically distinct. The fundamental mode's resonance frequency is fres,127.32 MHz, whereas higher harmonics do not expose a dedicated resonance frequency, as already noted for the cylinder.

Figure 21Modal significance sn(f).


Likewise, the modal excitation coefficients V¯ni of the first couple of modes, see Fig. 22, appear to be comparable to their cylinder equivalents in Fig. 8, apart from the frequency shift already discussed for the modal significance. Again, the first two modes are excited most and, additionally, J4H also has a broadbanded influence on the total surface current. Nevertheless, the second mode is excited less at most frequency points compared to the cylinder implying a smaller Fourier coefficient α¯2.

Figure 22Modal excitation coefficients |V¯ni(f)|.


Once again, an attempt to expand the total surface current in Characteristic Modes is made. Herefore, the car is excited identically to the cylinder, i.e. with the same electrical field strength and frequency, the latter being selected intentionally different from a resonance. The resulting surface current, see Fig. 23, exposes a maximum in the rear wheel region, which should be reproduced because of its potential EMC relevance.

Figure 23HIRF current |J¯HIRF(r)| for f=40 MHz.


With the help of the modal excitation coefficients preceedingly discussed the desired Fourier series is set up to

(12) J ¯ HIRF r n = 1 N α ¯ n J n H r = - 0.22 + j 0.10 J 1 H - 0.05 - j 3.1 × 10 - 3 J 2 H + 4.1 × 10 - 6 + j 1.3 × 10 - 6 J 3 H + 2.4 × 10 - 3 + j 1.5 × 10 - 3 J 4 H - 1.8 × 10 - 5 - j 2.7 × 10 - 5 J 5 H

using just five summands. It reveals larger errors in contrast to the cylinder expansion as depicted in Fig. 24. Beside the maximum in the rear wheel region, a further unintended maximum is generated at the front wheels. Apart from this region the current is quantitatively reproduced well by the series expansion. A possible explanation concerning the existing deviations is the larger complexity of the car's geometry such as the edges at which field excesses arise, which are more difficult to approximate.

Figure 24Fourier series of |J¯HIRF(r)| for N=5 and f=40 MHz.


An equivalence consideration between HIRF and DCI is initiated by a field excitation with the car's first resonance frequency fres,1, its result being illustrated in Fig. 25. The surface current J¯HIRF(r) is dominated by Characteristic Mode J1H because of a lack in significance of higher harmonics at this frequency.

Figure 25HIRF current |J¯HIRF(r)| for f=fres,1.


This is also supported by the series expansion in Characteristic Modes resulting in

(13) J ¯ HIRF r n = 1 N α ¯ n J n H r = - 0.31 - j 0.29 J 1 H + 0.04 + j 0.02 J 2 H - 5.7 × 10 - 7 - j 1.0 × 10 - 6 J 3 H + 3.3 × 10 - 3 - j 5.6 × 10 - 3 J 4 H - 2.1 × 10 - 6 - j 1.6 × 10 - 6 J 5 H

from which the mentioned mathematical dominance of the first mode can been also seen in Fig. 26 when being compared to Characteristic Mode J1H in Fig. 16.

Figure 26Fourier series of |J¯HIRF(r)| for N=5 and f=fres,1.


As DCI setup the same excitation scenario as with the cylinder is selected meaning a 50 Ω-source to the right and a high-impedance termination on the left hand side.

By visual comparison to the DCI current J¯DCI(r) in Fig. 27 an overtesting of the conducted setup can be noticed, e.g. at both wheels, but also in the area where the side doors are located.

Figure 27DCI current |J¯DCI(r)| for f=fres,1.


An expansion of the DCI-current in its affiliated Characteristic Modes gives

(14) J ¯ DCI r n = 1 N α ¯ n J n D r = - 0.06 + j 0.07 J 1 D + 0.02 - j 0.03 J 2 D - 7.6 × 10 - 7 + j 1.9 × 10 - 6 J 3 D + 3.5 × 10 - 3 + j 7.8 × 10 - 3 J 4 D + 2.6 × 10 - 7 + j 6.7 × 10 - 7 J 5 D

with Fig. 28 as graphical representation. Just like in the cylinder's case, the first two Characteristic Modes contribute most to the total current. However, the current maxima at the wheels and the door are not equally reproduced.

Figure 28Fourier series of |J¯DCI(r)| for N=5 and f=fres,1.


Computing the absolute error ε(ri) between HIRF- and DCI-excitation according to Eq. (11) followed by plotting the corresponding simulation mesh in Fig. 29 allows a numerical comparison. Unlike in the cylinder's case, a reasonable coincidence cannot be achieved. There still exists a significant surface current especially at the front wheel section and the door area.

Figure 29Absolute error ϵ(ri) of HIRF and DCI for f=fres,1.


Table 1 additionally lists the absolute values |α¯n| of all Fourier coefficients in Eqs. (8) to (10) and (12) to (14) allowing the contribution of single Characteristic Modes to be judged more clearly. Their order of magnitude directly indicates the degree of contribution to the total surface current despite a rather coarse truncation, which is of minor importance.

Table 1Absolute values of the Fourier coefficients α¯n.

Download Print Version | Download XLSX

5 Conclusions

Characteristic Mode Analysis was introduced and applied to three-dimensional test objects. Based on computed Characteristic Modes for HIRF- and DCI-excitations, series expansions for surface currents were set up. Its coefficients gave a quantitative insight into the contribution of single Characteristic Modes to surface current distributions at different test frequencies. It could be shown for geometrically simple objects that HIRF excitations can be characterized by only a very limited number of Characteristic Modes both at resonance and beyond. Furthermore, CMA revealed geometrically different DUTs to be of comparable electrical characteristics due to similar modal significances, provided their overall dimensions coincide. Additionally, CMA turned out to be useful in finding resonance frequencies of complex systems, whose excitation is of relevance for EMC tests. Comparable HIRF and DCI excitations were given and a quantitative comparison of both current distributions was made based on absolute errors.

Data availability

The data presented in this article are available from the authors upon request.

Author contributions

The computer simulations and data processing including all graphical representations were carried out by JÜ. Editorial hints and mathematical suggestions were given by FG.

Competing interests

The authors declare that they have no conflict of interest.

Special issue statement

This article is part of the special issue “Kleinheubacher Berichte 2019”. It is a result of the Kleinheubacher Berichte 2019, Miltenberg, Germany, 23–25 September 2019.


The authors thank Martin Aidam and Markus Rothenhäusler for helpful discussions.

Financial support

This research was supported by the Federal Ministry of Transport and Digital Infrastructure of Germany BMVI, in the funding program Automated and connected driving, LINKTEST (grant no. 16AVF2138A).

Review statement

This paper was edited by Lars Ole Fichte and reviewed by two anonymous referees.


61000-4-20: Emission and immunity testing in transverse electromagnetic (TEM) waveguides, DIN, 2010. a

61000-4-21: Testing and measurement techniques – Reverberation chamber test methods, DIN, 2011. a

61000-4-3: Testing and measurement techniques – Radiated, radio-frequency, electromagnetic field immunity test, DIN, 2011. a

ED-107A: Guide to Certification of Aircraft in a High Intensity Radiated Field (HIRF) Environment, European Organization for Civil Aviation Electronics, 2015. a

Harrington, R. and Mautz, J.: Theory of characteristic modes for conducting bodies, IEEE T. Antenn. Propag., 19, 622–628,, 1971. a

Harrington, R. F.: Time-Harmonic Electromagnetic Fields, Wiley-IEEE Press, 496 pp.,, 2001.  a

Leat, C.: The Safety of Aircraft Exposed to Electromagnetic Fields: HIRF Testing of Aircraft using Direct Current Injection, Australian Government Department of Defence, 2007. a

Vogel, M., Gampala, G., Ludick, D., Jakobus, U., and Reddy, C. J.: Characteristic Mode Analysis: Putting Physics back into Simulation, IEEE Antenn. Propag. M., 57, 307–317,, 2015. a