EEE 431 Computational Methods in Electrodynamics Lecture 14

EEE 431 Computational Methods in Electrodynamics Lecture 14 By Dr. Rasime Uyguroglu Rasime. uyguroglu@emu. edu. tr 1

FINITE DIFFERENCE TIME DOMAIN METHOD (FDTD) Lossy Material 2

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l When a material has a finite conductivity , a conduction current term is added to Ampere’s Law (different than the source term). Thus: 3

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l As before assuming the x component of E and the variation only in the z direction: 4

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l l This equation can be expanded about a point : to find FDTD update equation. However, when loss is present, the undifferentiated electric field appears on the left side of the equation. 5

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l l With the assumption that the electric field point is , there is no electric field at this space-time point. This problem can be solved by using the average of the electric field to the either side of the desired point: 6

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l Thus a suitable discretization of Ampere’s Law when loss is present is: 7

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l As before this can be solved for , which the present field, in term of purely past fields. The result is: 8

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l Now, consider the Faraday’s Law: l When, no magnetic loss is assumed 9

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l Now considering that there is varying with z direction: , 10

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l Discretized form: 11

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l l Solving for Yields the following update equation: 12

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l Consider a lossy dielectric half-space which starts at node 100. The relative permittivity is 9 as before. However there is also an electric loss present such that: 13

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l A previously employed simple ABC does not work at the right edge of the grid. It can be removed. But the left side of the grid can be terminated by using the same ABC as before. 14

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) l The magnetic field update remained the same. That us observe the pulse propagation for different time steps. (i. e. multiples of 10) 15

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 16

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 17

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 18

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 19

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 20

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 21

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 22

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 23

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 24

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 25

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 26

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 27

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 28

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 29

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 30

FINITE DIFFERENCE TIME DOMAIN METHOD ( Lossy Material) 31

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 32

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Material) 33

FINITE DIFFERENCE TIME DOMAIN METHOD (Lossy Dielectric) 34

FINITE DIFFERENCE TIME DOMAIN METHOD (TFSF-Lossy Material) l l l Conclusion: The pulse decays as it propagates in the lossy region and eventually decays to a rather negligible value. The lack of an ABC at the right side of the grid is not really a concern in this particular instance. 35