Numerical Simulation of a Dipole Antenna Coupling to a Thin Wire in the Near Field

1. Numerical Simulation of a Dipole Antenna Coupling to a Thin Wire in the Near Field Yaping Zhang, John Paul, C Christopoulos George Green Institute for Electromagnetics Research University of Nottingham, Nottingham NG7 2RD, UK www.nottingham.ac.uk/ggiemr/

2. Objectives of Working Group 1 within Cost 286: Near field electromagnetic coupling of a radiating antenna to a thin wire within a restricted space • Phase 1: Thin wire above a ground plane in an open area; • Phase 2: Enclosure problem, resonance; • Phase 3: Thin wire in an enclosure,cumulative effect of more sources.

3. Two configurations are simulated using the Transmission Line Modelling (TLM) method. The dipole antenna is parallel or perpendicular to the thin wire, as shown in Fig.1 (a)-(b). (a) (b) 100 50 half wave dipole antenna 100 50 0.3m half wave dipole antenna V V aligned h1=0.1m h2=0.3 m h1=0.1m h2=0.3 m 3 m 1m 1m 3 m y 1.5 m 1.5 m z 6 m 6 m x 0 Fig.1 Configurations of the simulation models: (a) a dipole antenna parallel to a thin wire above a ground plane; (b) a dipole antenna perpendicular to a thin wire above a ground plane. 100 V V Phase 1: Thin wire above a ground plane in an open area Configurations of the Simulation Models

4. The problem is considered taking into account the final configuration that needs to be tackled namely the case of coupling inside an imperfect cavity (vehicle with windows, non-metallic walls) and in the presence of people (very complex EM properties). Under these circumstances the advantages of MoM models in describing wires are no longer decisive and differential time-domain models such as TLM and FDTD appear to be advantageous! However, special techniques are required to describe thin wires (electrically small dimensions) inside cabinets (electrically large). We describe here efficient thin-wire formulations in connection with TLM.

5. Maxwell’s equations for electromagnetic wave propagation in free-space with conduction currents are described by Eq.(1): Solution of Eq.(1) is based on the 3-D TLM cell shown in Fig.2. The fields and current densities in Eq.(1) are normalised using the equivalences: wherel is the space step and 0isthe intrinsic impedance of free-space. The space and time derivative operators are normalised using: (2) Fig.2 Three dimensional TLM cell TLM Formulations

6. The maximum time-step in the three-dimensional TLM algorithm is t=l/(2c), where c is the speed of light in free-space. If the simulation is iterated at the maximum time-step, in terms of the normalisations Eqs.(2)-(3), Maxwell’s Equations Eq.(1) in free-space could be reduced to, This set of equations is then mapped into the cell in Fig.2. For free-space modelling with current densities, the calculation of the total fields at the centre of the cell can be expressed as, where the incident pulses are denoted by the superscript i. The normalised electric conduction currents {iex,iey,iez} are used to describe currents in wires. TLM Formulations (continued)

7. The incident voltages at the next time-step k+1 are founf by swapping with the reflected voltages at time step k. For example, the swap between ports 0 and 1 of the nodes located at cell indicates Z and Z-1 is, At external boundaries of a mesh (for example port 6 at Zmin), a reflection coefficient Rb is specified as, where, Rb =-1: a perfect electric conductor, Rb=1: a perfect magnetic conductor; Rb=0: a matched boundary, respectively. For description of frequency-dependent boundaries, Rbmay be specified as a frequency-dependent function. For a free-space simulations with current densities, equations (5)-(8) would apply to all nodes and ports in the mesh for the software implementation. TLM Formulations (continued) The reflected pulses are obtained from Eq. (6), where superscript r indicating reflection.

8. To aid the development of the model of field-wire coupling involving an x-directed straight wire, consider the shunt node associated with Vxderided from the3-D node. Fig.3shows this node and Fig. 4 shows the equivalent circuit. The voltage Vxis obtained as, Fig.4 Equivalent circuit of a shunt node describing Vx . Fig.3 Shunt node for evaluation of Vx 3-D TLM Implementation

9. For straight wires, the telegrapher’s equations describing an x-directed wire may be expressed as Eqs(10)-(11), where Vwis the wire potential, Iwx is the wire current and Lw /l, Rw /l and Cw /l are the wire inductance, resistance and capacitance per unit length respectively. The field-to-wire coupling is in the last terms of Eq.(10). A 3-D cell containing an x-directed wire is shown in Fig.5, which indicates the voltage pulses {V0, V1, V2 ,V3} interacting with the x-directed wire. The wire inductance Lw /l=0kw andwire capacitance Cw /l =0/kw ,where kw is a frequency-independent dimensionless geometrical factor. Fig. 5 3-D TLM cell containing an x-directed wire Formulations for Field-Wire Coupling

10. The velocity of propagation of charge variations along the wire is the speed of light, i.e. . The characteristic impedance of the wire is . The wire current is normalised using , where iwx has the dimension of volts. The time and space derivative operators are normalised using . Using the notations given above, Eqs.(10)-(11) may be rewritten as equations (12), where the normalised wire resistance rw=Rw/Zw . Equations (12) may be converted to the travelling wave format equations (13) by using the equivalences Eq.(9), The TLM equivalent circuit for the discrete-time solution of the equations (10)-(11) can therefore be written as, Formulations for Field-Wire Coupling (continued)

11. Fig. 7 Normalised equivalent circuit for field-wire coupling for an x-directed wire. Fig. 6 Equivalent circuit for field-wire coupling for an x-directed wire. Equivalent Circuit for Field-Wire Coupling Equivalent circuit for field-wire coupling for an x-directed wire is shown in Fig.6, and the normalised equivalent circuit for field-wire coupling for an x-directed wire is shown in Fig.7.

12. Simulation Results of a Dipole Antenna Coupling to a Thin Wire The configurations of the simulation models are shown in Fig.1 (a)-(b). The currents on the thin wire at points 25 cm, 50 cm and 75 cm from the junction of the wire, and the feed current at the gap of the dipole antenna excited by a sinusoid of 900 MHz or a Gaussian pulse with a halfwidth of 0.8 nm are calculated in both the parallel and perpendicular cases respectively by TLM method. The simulation results of the current transfer ratio (Iw/I_feed) for both configurations, excited by a sinusoid or a Gaussian function, are listed in Table 1. Table1 : Simulation Results of a Dipole Antenna Coupling to a Thin Wire Table 1 shows the current transfer ratios obtained from the steady-state and pulse excitations are in good agreement. The electromagnetic coupling is greater when a dipole is parallel to a thin wire.

13. (a) (b) (c) (d) Sinusoid Dipole Coupling to a Thin Wire (Parallel Configuration) The currents on the thin wire at points 25 cm, 50 cm and 75 cm to the wire junction and the feed current at the gap of the dipole antenna, excited by sinusoid source in a parallel configuration, are calculated and the results are shown in Figs.8(a)-(d). The current transfer ratios at points 25 cm, 50 cm and 75 cm may be obtained from these graphs and listed in Table 1. Fig.8 Currents excited by a sinusoid dipole parallel to a thin wire: (a) feed current at the gap of the sinusoid dipole; (b) current on the thin wire at a point 25 cm from the junction ; (c) current on the thin wire at a point 50 cm from the junction; (d) current on the thin wire at a point 75 cm from the junction

14. (a) (b) (c) (d) Sinusoid Dipole Coupling to a Thin Wire (Perpendicular Configuration) The currents on the thin wire at points 25 cm, 50 cm and 75 cm to the wire junction and the feed current at the gap of the dipole antenna, excited by sinusoid source in a perpendicular configuration, are calculated and the results are shown in Figs.9(a)-(d). The current transfer ratios at points 25 cm, 50 cm and 75 cm may be obtained from these graphs and listed in Table 1. Fig.9 Currents excited by a sinusoid dipole perpendicular to a thin wire: (a) feed current at the gap of the sinusoid dipole; (b) current on the thin wire at a point 25 cm from the junction ; (c) current on the thin wire at a point 50 cm from the junction; (d) current on the thin wire at a point 75 cm from the junction.

15. Hz Field Distribution in the Dipole-Thin Wire Plane

16. (a) (b) (c) (d) (e) (f) Fast Fourier Transform Analyses Results The current transfer ratios at points 25cm, 50cm and 75cm on the thin wire excited by a Gaussian function are obtained for both the parallel and perpendicular configurations respectively. Results obtained over the frequency range 600-1200 MHz are shown in Figs.10 (a)-(f). Fig.10 FFT analyses of the thin wire coupling to a dipole antenna excited by a Gaussian function in both parallel and perpendicular cases: Parallel case: (a) Current transfer ratio at 25cm; (b) Current transfer ratio at 50cm; (c) Current transfer ratio at 75cm; Perpendicular case: (a) Current transfer ratio at 25cm; (b) Current transfer ratio at 50cm; (c) Current transfer ratio at 75cm.

17. The thin-wire formulation technique employed here to simulate this problem is based on a quasi-static solution for the field around the wire. This is a symmetrical solution and therefore can only be used to place the wire centrally in the a computational cell. In cases where an offset wire description is required a more complete solution may be used which takes account of more modes (not just the quasi-static term included by Holland and Simpson) in the wire solution. The total field at the sampling points of the computational cell is decomposed into its modal components, these are then reflected by the correct modal impedance (see next slide) and then re-combined to obtain the reflected total field for transmission to the rest of the numerical solution.

18. “An accurate thin-wire model for 3D TLM simulation”, P Sewell, Y K Choong, C Christopoulos, TEMC Trans on EMC, 45(2), 2003, pp 207-217

20. Conclusions The electromagnetic coupling of a half wave dipole antenna to a thin wire is simulated using the TLM method. Two configurations of the dipole positions are considered, in which the antenna is positioned either parallel to, or perpendicular to, the thin wire. Results obtained by using the steady-state and pulse excitations agree closely. Results are obtained by combining the versatility of the TLM with the efficiency of powerful thin-wire formulations.