Propagation of Current Waves along Randomly Located Multiconductor Transmission Lines inside a Rectangular Resonator

Propagation of current waves along randomly located multiconductor transmission lines inside a rectangular resonator is investigated. Two different line configurations are considered: parallel wires and perpendicular lines. The analytical solutions for both types of lines are obtained, using the method of symmetrical lines inside the resonator. Computer realization of this method allows it to obtain a fast solution. Therefore, it becomes possible to investigate statistical properties of induced currents. It could be shown that the mutual coupling of the currents in the wires is strong when the frequency is near resonances in the system “cavity and wires” and relatively small far from these resonances. Stochastization leads to a broadening of the resonance curves.


Introduction
Transmission lines of different kinds are widely used in different branches of modern technique and industry, especially in power engineering and electronics.They can serve as a channel for the propagation of useful signals and energy, as well as a receiver of different kind of intentional and unintentional electromagnetic interferences.Thus, the study of signal propagation along transmission lines and the external electromagnetic field coupling with transmission lines and wiring systems is an important task of electromagnetic compatibility.This problem is complex itself, especially for the high frequency case, when the transversal dimension of the transmission line is comparable with the wavelength, and also for the treatment of multiconductor lines.Moreover, this task is additionally complicated by the fact that transmission lines are often placed inside resonator-like objects: car bodies, aircraft fuselages, computer cases, shielded rooms, etc. Due to the resonator behavior there is a strong coupling of the penetrated field with the wires inside and vice versa; this can greatly affect the signal propagation on conductors.In addition, often the exact positions of the lines and the corresponding electrical parameters per-unit length are unknown and have randomly distributed values.This leads to the necessity to study stochastic transmission lines in the resonator, including the more complex multiconductor variant.
In literature, there are a number of papers which deal with stochastic cavities, i.e., with cavities whose parameters are only known stochastically.This activity is mainly connected with the practical necessity to study Mode Stirred Chambers (MSC) and to clarify the correspondence between measurements in a MSC and in free space.For those studies papers of S. Hemmadi and T.M. Antonsen and their groups [1,2] are well known, which are based on a deep investigation of the mathematical properties of stochastic cavities [3], when they investigate the response of components in stochastic cavities.The authors of [4] consider stochastically located small antennas to create stochastic fields in a MSC.Other papers (see, e.g., [5][6][7]) are devoted to the investigation of coupling of a random electromagnetic field to deterministic transmission lines inside stochastic cavities (MSC).
In contrast, we are interested in stochastic lines inside the resonator with deterministic parameters.This view is linked to the need for a statistical analysis of the vulnerability of electronic devices to intentional and nonintentional interferences, mutual coupling interferences, etc.We only know one paper which deals with a single-wire line of stochastic geometry inside a resonator [8].In [8], the transmission line is considered inside a metallic rectangular box with a rectangular aperture.The system is excited by an incident plane wave on the aperture.The transmission line has a fixed length and positions of connectors but different layouts (or pats) for each statistical realization.The finitedifference time-domain method (FDTD) and fast Fourier transform were used to obtain the voltages at the terminals and their statistical parameters: average value and dispersion.It was shown that maxima of these averaged voltages occurred next to the cavity resonances.It is physically obvious that the dispersion of these values is defined by the stochastic characteristic of the line geometry.However, the direct numerical methods are very time consuming.Therefore in [8] only 12 statistical realizations were chosen.It seems to be clear that for a deep investigation of this model some analytical or analytic-numerical method has to be developed.
Recently, we proposed an exact solution method for a loaded line in a resonator, which can be excited by an arbitrary electromagnetic field [9,10].From the computational point of view, the solution is reduced to the calculation of a double modal sum with 200-300 terms for each sum direction and for each considered frequency.This method was applied to a multiloaded line or to two parallel lines [11] excited by a lumped source(s) and by distributed sources, such as additional internal radiating antennas or by EM fields which penetrated into the resonator through apertures [12].
Note that our method is applicable to a line with geometric symmetry of the resonator: the line connects two opposite rectangular walls and it is parallel to the four other walls (see Figure 1).Nevertheless, we have applied this method for the investigation of stochastic lines inside the resonator.
In [13,14], one single-wire line in the resonator was considered.The stochastization was created by multiple loading with random positions along the line or by random values of the reactive loads.With an appropriate choice of parameters, this model can serve to describe a line with stochastic geometry.It was shown that for a stochastic line in the resonator, unlike in free space, the current wave on the "average" passes through the transmission line.This is explained by the multiple rereflection of the signal from walls of the resonator.Also, it was shown that the PDF of the absolute value of the transmission coefficient for current wave for the case of strong stochastization (when the average value of transmission coefficient equals zero) can be described by the binominal distribution, as the consequence of the principle of maximum entropy.
In the present paper we extend the method [9, 10] to multiconductor lines in the resonator.Stochastization is created by random positions of the lines.The rest of the paper is organized as follows: In Section 2 we derive equations for the EM coupling of parallel multiconductor loaded lines in a rectangular resonator.Then the corresponding solution will be applied to the lines with random transverse coordinates and for the calculation of the corresponding lower statistical moments.In Section 3 we consider two perpendicular wires in the resonator.Again, the general equations for the currents induced in both wires for arbitrary excitations are obtained.Unlike the previous case, the symmetry of the problem is not conserved.The modal matrix equations which describe the induced currents do not split.For a solution operations with infinite matrices are required.However, for a relatively small number of accounted modes (about 300-500) one can obtain the solution quite fast.The obtained solutions again are applied for the investigation of statistical moments of the induced currents for random positions of the lines.Section 4 concludes the paper.

Exact Analytical Solution of EFIE for Short-Circuited Lines.
We consider N parallel wires inside a rectangular resonator with sides a, b, and h.The wires keep the symmetry of the resonator: they connect two opposite walls of the resonator and they are parallel to all other walls and to the z axis (see Figure 1).Further, assuming that all wires are thin-the radius of the each wire r 0 is essentially smaller than all linear dimensions of the problem including wavelength -we can apply so called thin-wire approximation.In this approximation only axial components of the current are taken into account.Moreover, it is assumed that all current densities are concentrated along the wire axis; however, the zero-boundary conditions for the perfectly conducting wires have to be satisfied for the total tangential electric field on the boundary of the each wire   .
The scattered field in the considered geometry can be calculated using zz components of the resonator's Green's function "electrical current->electrical field",    .Thus, we can establish a system of integrodifferential equations that determine the induced current in the wires: The zz component Green function for the electric field, in turn, can be evaluated from the zz component of Green's function for the vector potential    : where We have introduced the following notations: and In the denominator of ( 5) the summand  sign()with small  >0 shifts the pole  2 =  2 V in the complex frequency plane in the upper half plane of the complex variable k fl /.This is necessary, e.g., during the calculation of the system response in time domain by the inverse Fourier transform, when one has to perform the integration in the complex plane by applying the residue theorem.In the paper [15], e.g., this small quantity  is assigned a physical interpretation, e.g., small radiation losses through an aperture or small losses in the walls.In other words, this quantity is connected with quality factor of the resonator Q by  =  2 / [16].Since further in this paper we will only consider resonators without losses in the frequency domain, the factor  will be omitted.

Thus, one can write the function 𝐺
where we introduced the notation for the double sum: First we consider short-circuited wires (see Figure 1, where all loads are zero), where the induced currents satisfy the zero Neumann boundary conditions:   ()/| =0 = ()/| =ℎ = 0 with n=1,2..N.This allows decomposing each current into the Fourier transform: or Then we decompose the exciting electrical field into a Fourier series for each wire in (3).(One can make this decomposition due to the zero boundary conditions for the tangential components of electrical field on the walls) where the Fourier coefficients (10a), (10b) can be obtained by integration: Then, substituting (9a), (9b), (10a), (10b), and ( 7)-( 8) into (3) and taking into account the orthogonality of the cosine functions, we can obtain the linear system for the Fourier images of the unknown currents: Introducing the matrix of two-dimensional scattering S(,   ), we can rewrite (13) in the matrix form: where [S(,   )] , 1 = (,   ,   +  0 ,   ,   1 ,   1 ).
Then, in the definition of the function (,   ,  →  ,  →  1 ) (8) one summation can be done analytically (see [17]) and, after some manipulations, we can write expressions for the diagonal and nondiagonal elements of the matrix S(,   ), as a one-dimensional sum [9]: Two wires in resonator One wire in resonator One wire, TL approximation The solutions for the modal amplitude and the z-dependence of the vector-current can be obtained using ( 13), (14), and (15): With the aid of (3), (17a) and (17b) can be rewritten as follows: where we introduced Green's function matrix for the multiconductor transmission line Equations ( 18)- (19) together with ( 15)-( 16) yield the solution for the current in the short-circuited multiconductor transmission line for an arbitrary exciting field.One can show that if the frequency is far from cavity resonances and the wires are close to one of the cavity walls, these equations yield the classical TL approximation solution for the short-circuited multiconductor line obtained by Fourier transformation. In The transverse coordinates of the second wire are x 1 =20 cm, y 1 =50 cm.The second wire is also shortcircuited at left and right terminal.
(For the numerical simulation we used the values Z left 2 = Z right 2 =0.2 Ω.)The radius of both wires is r 0 =1 mm.In Figure 2 the current at the right terminal of the first line is presented, calculated by an analytical method with and without accounting for the second wire, also in the transmission line approximation [9].In Figures 3 and 4 the currents at the right terminals of the first and second wire are shown, respectively, together with a comparison of the results obtained by the code CST MICROWAVE STUDIO [18].Since the model includes a resonator, the frequency domain solver has been used based on the finite element method (FEM).From the analysis of the calculations one can observe the following: The analytical method for the resonator and the direct numerical calculation yield a good agreement for the considered system.Far from the cavity resonances the current in the second, passive wire is smaller in comparison to the current in the first one, but near the resonances they can have the same order of magnitude.In turn, the influence of the second wire on the first one is quite small for a frequency interval far from cavity resonances.In this case the current in the first wire can be described by the usual TL approximation.At the same time, near the resonances the influence of the second wire is essential.The response of the first wire alone and in the presence of the second one is strongly different.

Exact Analytical Solution for MTL: Accounting for Lumped
Excitations and Loads.Consider now the excitation of the transmission line by lumped sources which have the same coordinate z 1 , but different amplitudes: where Then, from (18) the current is given by Now it is possible to take into account also lumped impedance at the left and right terminals of the line.The lumped loads at the left terminal at the point z 1 =Δ (Δ →0) can be considered as controlled lumped voltage sources with unknown amplitudes: or where The same trick can be done for the lumped loads at the right terminal z 2 =h-Δ (Δ →0): where ].Then, the column vector for total exciting field for the loaded multiconductor line with lumped sources at the left terminal is which yields matrix equation for the current: After that calculating column vector of currents at the points z 1 and z 2 , we obtain linear system (26) for the values I( 1 )and I( 2 ),which, in turn, defines current in arbitrary point by ( 25): Below we presented a numerical example for the response of the loaded lines (see Figures 5 and 6).We consider a system with the same parameters resonator and lines, except the nonzero load at the right terminal of the first line.It is loaded by the matched load for the single line for TEM waves in resonator: As in a previous example, the comparison has shown a good agreement between developed theory and numerical simulation.Far from the resonances current in the first wire is approximately constant, i.e., has TL structure (see [9,10]).Also, as in previous example, the mutual influence of the wires is strong near resonances.Note that near the resonances current in the second wire can be comparable and even more than the current in the first wire.

Randomly Located Multiconductor-Transmission Line.
In this section to realize stochastization of multiconductor TL, we consider random positions of the parallel wires.We consider a similar resonator with dimensions a=1.5 m, b=1.2 m, and h=0.9 m.As in the previous case, the first line has a left lumped source; both lines have zero loads at both terminals.The radiuses of the wires are 1 mm.The transversal positions of the both wires in the plane (x,y) are defined by  01 =  ⋅ (0.05 + 0.9 ⋅ ) ,  01 =  ⋅ (0.05 + 0.9 ⋅ )  02 =  ⋅ (0.05 + 0.9 ⋅ ) ,  02 =  ⋅ (0.05 + 0.9 ⋅ ) where r is a stochastic value uniformly distributed in the interval 0 ≤  ≤ 1.
The simulation was carried out in the frequency band 150-550 MHz with 500 frequency points.For each frequency point we made 100 statistical events.All calculation required about 4 hours on the computer HP ENVY 17 Notebook PC with processor Intel(R) Core(TM) i7-5500U@ 2.40 GHz.Code was realized on the Fortran programming language.
The results of simulations for the right terminal of the first (active) and the second (receiving) wires are represented at Figures 7-10 both for the average absolute value of the induced current for unit source (response function) and for the variance of these values.For comparison, also we present response functions for the deterministic lines with parameters x 1 =9 cm, y 1 = 39 cm, x 2 =20 cm, and y 2 =50 cm.From the analysis of the curves one can conclude the following.For the active wire the stochastization changes the positions of the response curves near cavity resonances (system cavity + two wires), but not near the transmission line resonances (see Figure 7).For the receiving wire the influence of stochastization is more strong, but, in principle, all features of previous case are kept (see Figure 9).The statistical variance of both response functions, for the active and receiving wires, has the same order of magnitude: see Figures 8 and 10.Note that from the relative small number of statistical events the frequency dependence of the average values and variances are not smooth.

Exact Analytical Solution of EFIE for Short-Circuited
Perpendicular Lines.In the present section, we consider another configuration of multiconductor lines with symmetrical geometry in rectangular resonator, which can be interesting for the practical applications on the one hand and allows the analytical consideration on the other hand.Namely, we consider two perpendicular lines in resonator (see Figure 11).As in the previous case, we expand the exciting tangential electrical fields and induced currents into the Fourier series, but now the basic cosines functions are different for the wires of different directions: Then the z and x components of the scattered fields in the arbitrary point with coordinates (x,y,z) inside resonator can be written as follows: where    (  →  ,  →   ) are components of tensor Green's function "electrical current -> electrical field" for the rectangular resonator.In particular, the component    (  →  ,  →   ) is given by ( 7)-( 8) and and Then one can write zero boundary conditions for tangential components of the total electric field on the boundary of the first and second wire by a standard way (of course, it is again assumed that the thin-wire approximation is valid): and, taking into account the orthogonality of the system of cosines functions, one can obtain the system of electrical field integral equations for the two-wire system in the modal representation: Here the elements   are given by During the calculation of matrices S summation over one index can be carried out analytically (like in ( 15) and ( 16)).Note that, unlike the case of two parallel wires, considered in the previous section, the amplitudes of current modes for different wire couple with each other; i.e., the system of equations (34a) and (34b) does not split.By this way, to solve this system two methods can be proposed.
The first method, analytical, is an iteration approach.Due to the lack of space in this paper, we will only briefly describe this method.On the first step, the nondiagonal terms in (34a) and (34b) are moved in the right side.Then, as the zero iteration one can consider the noninteracting wires, i.e., assuming that the nondiagonal terms equals zero.Then, each next iteration includes these nondiagonal terms with current modes of previous iteration as sources.
The second method, more numerical, reduces the equation to the linear system for current amplitudes.Restricting the summation by some index M max (the calculations have shown that the value M max =500 describes the current with good accuracy), one can obtain a linear system of the 2M max × 2M max order: where we use the following notations for the M max × M max matrices and M max -dimension column vectors: By this way, the solution of (36) yields currents induced in the short-circuited perpendicular lines with symmetrical geometry in rectangular resonator by arbitrary (lumped or distributed) sources.The loads at the terminals of the line can be introduced by the usual way, as controlled voltage sources, as it was described in Section 2.2.The numerical results for the deterministic lines are presented in Figures 12 and 13 (1) Our method and numerical calculations have a satisfactory agreement for the considered system of two wires.
(2) Far from the cavity resonances the influence of the second wire on the first one is quite small.The current in the first wire can be described by usual TL approximation (see Figure 12).At the same time, near the resonances the influence of the second wire is essential and response of the first wire alone and in the presence of the second one is strongly different.
(3) Far from the cavity resonances (|- V | ≪ Δ, where value Δ is a shift of resonator eigenfrequency caused by the presence of transmission lines in resonator; for the case of one wire in resonator it can be evaluated as Δ the current in the second (receive wire) is smaller in comparison with the current in the first one, but near the resonances they can have the same order of magnitude.
(4) The influence of the cavity is essential in the frame of the width of the resonance Δ (in the system cavity + two wires).For the large number of resonances the distance between them becomes smaller, but the widths of resonances are the same.By this way, the cavity defines the frequency dependence of coupling for large frequencies.

Randomly Located Perpendicular Wires.
Using the developed method one can investigate (within a reasonable time) the statistical properties of induced currents and fields for the case of stochastically arranged symmetrical lines in rectangular resonators.The results of such investigations are presented in Figures 14-17.As in the previous case, the positions of both lines were arranged stochastically:   1 st transmission line ‖ z-axis, x 01 =a⋅(0.05+0.9⋅r),y 01 =b⋅(0.05+0.9⋅r) 2 nd transmission line ‖ x-axis, y 02 =b⋅(0.05+0.9⋅r),z 02 =h⋅(0.05+0.9⋅r),where 0<r <1 is uniformly distributed stochastic variable.Again, for each frequency point we made 100 statistical events.Analysis of the curves has shown the next characteristic features.As in the case of parallel wires, the stochastization changes the positions of resonances near the cavity ones, but not the positions of transmission line resonances.For the receiving wire the influence of stochastization is more strong (see Figures 16 and 17).As in deterministic case, the average response of the receiving wire is much smaller (one order of magnitude) than the response of the active wire.At the same time, the statistical variance for the active and receiving wire has the same order of magnitude (Figures 15 and 17).

Conclusion
In this work the propagation of current waves along stochastic transmission lines inside a resonator was investigated.Two wire configurations were considered: parallel straight wires and perpendicular straight wires.These configurations keep the symmetry of the resonator: the wires connect two opposite walls of the resonator and are parallel to the other four walls.The stochastic line was created by randomizing the positions of the wires.To find the current, an earlier developed method for symmetrical wires inside the resonator was applied for the multiconductor problem.This method gives the general solver for arbitrary excitations, including the considered case of excitation of one wire by a lumped voltage source.
For both geometries it was shown that the influence of the receiving line on the active line is small, when the frequency is far from cavity resonances, but it is essential near these resonances.These effects are kept, for both deterministic and stochastic lines.
The stochastization changes the positions of resonances of the system "wires in resonator" near the cavity ones, but not the positions of transmission line resonances.For the receiving wire the influence of stochastization on the resonant picture is stronger.Unlike the mean value, the statistical variance for the active and receiving wire has the same order of magnitude.Qualitatively the frequency dependency of the variance corresponds to those of paper [8].
The obtained results have several fields of application.They can be used to evaluate the statistical properties of parasitic mutual coupling for the propagation of the signal along a TL in resonators; one can calculate the response of the lines to the external excitation by the EM field penetrating into the resonator through slits and apertures and can evaluate the corresponding statistical properties.Also the results can be used to investigate the damping of the penetrated field by the scattering of the loaded lines (see [19]).

Figure 1 :
Figure 1: Two-wire symmetrically loaded line in a resonator.

Figure 2 :
Figure 2: Deterministic parallel lines in a resonator.Response at the end of the first line calculated by TL approximation and the exact solution with and without accounting for the second line.

Figure 3 :Figure 4 :
Figure 3: Deterministic parallel lines in resonator.Response at the end of the first line calculated by the exact analytical solution with CST code.

Figure 5 :
Figure 5: Deterministic parallel lines in resonator.Response at the end of the first, loaded, line calculated by the exact analytical solution and by the CST code.

Figure 6 :
Figure 6: Deterministic parallel lines in resonator.Response at the end of the second, short-circuited, line calculated by the exact analytical solution and by the CST code.

Figure 7 : 2 )Figure 8 :
Figure 7: Stochastic parallel lines in resonator.Averaged response function at the end of the first, active, line compared with the deterministic result.

Figure 9 : 2 )Figure 10 :
Figure 9: Stochastic parallel lines in resonator.Averaged response function at the end of the second, receiving, line compared with the deterministic result.

Figure 11 :
Figure 11: Perpendicular wires with symmetrical geometry in resonator.

Figure 12 :
Figure 12: Current in the first, active wire at the point z=0 with and without the second, receiving wire.

Figure 13 :
Figure 13: Comparison of the currents in the perpendicular wires at the points z=0 and x=0, correspondingly.

Figure 14 :Figure 15 :Figure 16 :Figure 17 :
Figure 14: Stochastic perpendicular lines in resonator.Averaged response function at the end of the first, active, line compared with the deterministic result.
The positions of the n th wire in the xy plane is given by  →   = (  ,   ).It is assumed that there is an