The Terminal Responses of the Two-Wire Line in Multiaperture Cavities Based on Electromagnetic Topology and Method of Moments

A simulation technique based on electromagnetic topology EMT theory is proposed for analyzing electromagnetic interference EMI coupling through apertures onto the two-transmission line enclosed within metallic structures. The electromagnetic interactions between apertures and the external-internal interactions were treated through the topological decomposition and the multistep iterative method. Then, the load responses of the two-wire transmission line are resolved by the the Baum-Liu-Tesche BLT equation. The simulation results both without and with the electromagnetic interaction are presented for the frequency range from 100MHz to 3GHz. These numerical results obtained by two methods imply that the electromagnetic interaction cannot be simply ignored, especially for the frequency range up to 1GHz.


Introduction
While considering analysis and design of EMC, it is important to protect electronic circuits/components from electromagnetic effects due to external illumination.For this purpose, electronic circuits/components are often shielded inside metallic cavities.However, these metallic surfaces are not perfect.On a metallic cavity, there may be apertures that become sources of the EMI problem.In most cases the cavity and apertures are usually rectangular, which has led to a number of attempts to solve the coupling problem of the rectangular cavity with a series of rectangular apertures when illuminated by a plane wave see 1-3 .So there would be interference between the fields entering through the apertures and the circuitry inside the cavities.A natural problem is a comprehensive analysis of the coupling mechanism for such structures under external EMI illumination.In this paper, electromagnetic interaction problems can be simulated through codes based on electromagnetic topology.An important element of topological analysis is the determination of a mechanism which represents the external-internal coupling through a small aperture and the subsequent propagation process.By combining methodologies suggested earlier see 4 , an equivalent source can be created to relate the electromagnetic coupling at the exterior surface and the transfer function through the free space by assuming an imaginary transmission line as a source of the aperture radiation see 5, 6 , but this method only adapts to the case where the dimension of the aperture is electrically small when compared to the wavelength.Moreover most papers in the literature neglect the electromagnetic interaction among apertures see 3, 7 , which is not reasonable.This paper provides a novel numerical technique for the electromagnetic coupling.More precisely, we first use the EMT theory and the multistep iteration to deal with problems of coupling between apertures and the external-internal interaction.Then, we employ dyadic Green's functions and the method of moments to determine the electromagnetic coupling fields inside the cavity.Finally, we apply the BLT equation to resolve the load response of the two-wire transmission line.The simulation result shows the effect of the external coupling fields on the two-wire line current depends on both the cavity with multiaperture and the electromagnetic interaction.

The Topological Decomposition of Complicated Electrical Systems
Electromagnetic interaction problems on very large and complex system, such as an aircraft, can be simulated through codes based on the EMT.The most important aspect of the EMT is the assumption that volumes can be decomposed into subvolumes that can be interacted with each other through apertures see 8 .To analyze the interaction processes for the electromagnetic coupling by the EMT, we need to establish the topological structure model and the topological diagram for the system configuration.First, we consider a rectangular cavity with R energy penetration paths multiaperture illuminated by a harmonic plane wave, see Figure 1. Figure 2 shows the topological structure model associated to Figure 1.V i denote these subvolumes, where "i" indicates the hierarchical order of the volume.S i,j is a surface separating volumes V i from V j .Figure 3 shows the topological diagram.To analyze the field coupling phenomenon in the cavity, we introduce a "field coupling junction" J 3 , which takes into account the internal EM field coupling to the two-wire line.Moreover, by Aperture interaction Field excitation of aperture

Load junction
Incident EM field Field excitation of line

Figure 2:
The system topological structure model.
Tube reciprocity, this junction also provides effects of the EM field radiation by the two-wire line see 9 .W i 0 is the outgoing wave and W s is the network cable coupling source.Signals on an entire transmission line network is expressed through the BLT Baum-Liu-Tesche equation which is the multiconductor transmission line MTL network composed of the outgoing wave supervector W 0 and the source wave supervector W s see 10 .

The Multistep Iterative Method
Based on the topology decomposition, we deal with problems of the electromagnetic interaction among apertures and the external-internal interaction using the multistep iterative method see 11 .We assume R apertures on the same wall located on the plane where z 0. R is the number of apertures.
Step 1. Find the zeroth-order approximation for R-aperture magnetic currents M 0 and the corresponding coupling electromagnetic fields E 0 , H 0 in a cavity due to M 0 .
Theoretically, a plane wave excitation represents the simplest electromagnetic source and therefore is particularly suitable to test numerical techniques.In this paper, without lose of practical applications, we will use a harmonic plane wave with a wide frequency spectrum to simplify the problem under consideration.Consider a rectangular cavity with R rectangular apertures illuminated by a harmonic plane wave shown in Figure 1.This field is described by angles of incidence ψ and φ, as well as a polarization angles α, which defines the E-field vector direction with respect to the vertical plane of the incidence.It is given as where r x x y y z z refers to the location of the field in the free space, the free space wave number, and Z 0 is the free space intrinsic impedance.For a unity amplitude electric field, the coefficients H 0x , H 0y , H 0z are given by We introduce equivalent magnetic currents M 0 on R apertures as where M 0r is an equivalent magnetic current on the rth aperture.The equivalent magnetic current components now are expanded as where in which x a r , y b r is the coordinates of the lower left-hand corner of the rth aperture, p, q are integers used for multiaperture modes, T p t and P q t are triangular and pulse functions defined by for p 1, 2, . . ., P r − 1, and for q 1, 2, . . ., Q r .Substituting 3.7 into expression 3.6 , we derive equivalent magnetic currents as To solve the unknown currents M 0 , we use the continuity of tangential magnetic fields across R apertures.Then, we have where H a refers to the corresponding exterior scattered magnetic field due to M 0 , and H b refers to the corresponding interior coupling magnetic field due to M 0 .The exterior scattered field can be expressed as the radiation caused by the equivalent magnetic current M 0 , say, where S r denotes the surface of the rth aperture and Γ0 r; r is the dyadic Green's function in the free space, which reads

3.15
in which r and r represent the locations of both the field and source points on the aperture, respectively.
The interior field can be formulated as the radiation due to M 0 on R apertures.Using the available dyadic Green's functions for the cavity, we can express the interior field as 16 where the dyadic Green's function is defined as where δ mn denotes the Kronecker delta, say, δ mn 1 for m or n 0, and zero otherwise.Also, the wave functions within representation 3.18 are given by

3.20
with m, n being nonnegative integers excluding m n 0.
Mathematical Problems in Engineering 7 Substituting 3.18 into expressions 3.16 and 3.17 yields all components of the electromagnetic field inside the cavity as

3.22
M 0rxpq and M 0rypq are the unknown constants of the pq mode on the rth aperture, w r , l r are the width and length of the rth aperture, respectively.Δx r , Δy r refer to the width and length of the surface element used to discretize the rth aperture.Using Galerkin's method see 12 , the integral equation to be solved for M 0rxpq , where W is a weighting function, H a M 0 and H b M 0 are given by 3.14 and 3.17 , respectively.To discretize 3.23 , the corresponding weighing functions are given by W r p q Φ r p q x, y x Ψ r p q x, y y.

3.24
where Φ r p q x, y and Ψ r p q x, y are given by 3.8 and 3.9 , respectively.Substitution of 3.14 and 3.17 into 3.23 yields a matrix equation:

3.25
The expressions of all matrix entries in 3.25 are as follows.
Self-admittance matrix Y b is represented by abk mn tan k mn c I rypqmn I r yp q mn ,

3.26
where ε m is Neumann's numbers, and I rxpqmn , I rypqmn are given by 3.22 .

3.29
where sinc x sin x /x, H 0x and H 0y are given by 3.3 and 3.4 , respectively.The solution of this matrix equation 3.25 yields the coefficients M 0rxpq , M 0rypq of R-aperture magnetic currents M 0 .A similar result can be found in 3, 11 .
Step 2. Find R-aperture magnetic current increment ΔM To solve the unknown magnetic current increment ΔM 1 , we require the tangential magnetic field be continuous in the following sense: where 31 M 0 and W are given by 3.17 and 3.24 , respectively.Substitution of 3.17 , 3.24 , and 3.31 into 3.30 yields a matrix equation: The solution of this matrix equation 3.32 will yield the coefficients ΔM 1rxpq , ΔM 1rypq of the magnetic current increment ΔM 1 .
Step 3. Find R-aperture magnetic current increment ΔM 2 due to ΔE and a matrix equation is obtained as

3.34
So the solution of this matrix equation 3.34 will yield the coefficients ΔM 2rxpq , ΔM 2rypq of the magnetic current increment ΔM 2 .
Step 4. Find R-aperture magnetic current increment ΔM 3 due to ΔE

Similarly
and a matrix equation is obtained as: The two-wire line So the solution of this matrix equation 3.36 will yield the coefficients ΔM 3rxpq , ΔM 3rypq of the magnetic current increment ΔM 3 , and so on.Finally, this iterative process approximates the magnetic currents on R apertures and the total coupling electromagnetic field inside the cavity as

3.39
We must point out that the four-step algorithm admits high accuracy for the approximation to the interaction among apertures and the external-internal interaction.

The Application of the BLT Equation
The EMT can simulate the response of inner circuits of the electrically large complex system.The BLT equation is the key equation of the EMT to express the signals on the entire transmission line network 13 .Hence we use the BLT equation to solve the terminal responses of the two-wire line in a multiaperture cavity illuminated by a plane wave.We consider the case where a lossless two-wire line is illuminated by a coupling EM field.In order to simplify the notation, let us place the two-wire transmission line in the plane where y y 0 , which parallels to the x-axis as shown in Figure 4. We put the reference wire at the plane where z z 0 , and the other wire at z z 0 d so that there is distance d between two wires.We assume that d is greater than the wire radius r a d r a .ρ i Z i − Z c / Z i Z c i 1, 2 are the reflection coefficients at each node of the line.L is the length of the line.Z 1 and Z 2 are the load impedances at the ends x x 0 and x x 0 L, respectively, and Z c 120 ln d/r a is the characteristic impedance of the line.The transmission line theory shows that the wave propagation constant γ jk 0 .

Mathematical Problems in Engineering
In this paper, we only consider the transmission line mode current.The BLT equation of the load currents and the total voltages at the loads can be expressed in the matrix form as

4.1
We choose the Taylor formulation see 14, 15 , since it consists of both voltage and current sources that are appropriate for localized excitations through apertures and the source vector given by

Simulation
We present some numerical examples to demonstrate the formulation given in the previous section.We consider the coupling in a cavity with two rectangular apertures on the plane z 0. There is a two-wire transmission line along the x direction in this cavity, see Figure 1.The two-wire transmission line can be excited by the interior coupling electromagnetic field.By employing the BLT equation, Figure 5 shows the results of the terminal response at Junction 4 with those parameters.This example indicates that the effect of the external coupling fields on the two-wire line current depends on both the multiaperture cavity and the electromagnetic interaction among apertures.More precisely, the difference between the results with the electromagnetic interaction and those without the electromagnetic interaction is actually not small at all.Therefore, this electromagnetic interaction cannot be simply neglected for the frequency range up to 1 GHz.Moreover, using the BLT equation, we can also conclude that the currents and voltages at arbitrary point on the two-wire line can be derived from the interior coupling EM fields.

Conclusion
In this paper, the field coupling phenomenon between external fields and a two-wire line in the multiaperture cavity is studied.For the issue of the field penetration through apertures, the EMT theory and the multistep iteration are used to deal with these problems.This method can also be used for electromagnetic interaction problems on more complex system.Then, using the Modal Green's Function and the method of moments, the electromagnetic total coupling fields are determined inside a multiaperture cavity.Finally, the load response of the two-wire transmission line could be resolved by the BLT equation.The results with the electromagnetic interaction are largely different from those without the electromagnetic interaction for the low frequency.

Figure 1 :
Figure 1: The rectangular cavity with R rectangular apertures.

Figure 4 :
Figure 4: The isolated two-wire line excited by the interior coupling field.

2 inSubstituting 3 .
which voltage and current sources V s x and I s x are given by 38 into 4.4 and 3.39 into 4.3 then into 4.2 can determine the source vector.Thereafter, using the BLT equation, we obtain the induced voltage and current of a two-wire line at each load in a multiaperture cavity under excitation.

Figure 5 :
Figure 5: Comparison of induced currents of a two-wire line at load under excitation.
and I rxpqmn 8sin 2 k xm Δx r /2 sin k yn Δy r /2 rpq dk xm dk yn , 3.27 where F rpq k xm , k yn and G rpq k xm , k yn are the Fourier transforms of Ψ rpq x, y and Φ rpq x, y , respectively.By replacing k xm , k yn by −k xm , −k yn , the expressions F * rpq and G * rpq can be obtained from expressions of F rpq and G rpq , respectively.Admittance matrix Y a b is represented by