Fdtd algorithm
The algorithm can be expressed in a fully explicit form, i.e., equation 6 is written asĪ detailed description of the entire FDTD updating equations for electric and magnetic field components is given in Allen Taflove’s standard FDTD book. The x component of the electric field within equation 3 is formulated for FDTD by the application of central differences as follows: O and O indicate the spatial and temporal second order error terms. The indices i, j, k indicate integer values specifying a certain node in the spatial FDTD grid and u n any temporal dependent function at the time nΔt, i.e., the electric or magnetic field. Achieving 2nd order accuracy using a 2nd order finite differences approximation for the differential equations 2 and 3 leads toįor time and space as shown for the x direction. The original FDTD formulation proposed by Yee was intended for homogeneous, isotropic and lossless media on a uniform grid by application of a cartesian coordinate system. The computation of the magnetic components in FDTD is performed at t = (n + l + 1/2)Δt, (l = 0, 1, 2. ) with respect to a globally defined timestep Δt. Thereby the electric components are calculated at incrementing time t = (n+l)Δt, (l = 0, 1, 2. The discretization in time for FDTD is performed in a leap frog manner by application of a temporally shifted updating for E- and H-field components as shown in Fig. Each three E- and H-field components are assigned to a node i, j, k within the three dimensional (3-D) FDTD grid. The fields are located in a way in which each E component is surrounded by four H components and vice versa, which leads to a spatially coupled system of field circulations corresponding to the law of Faraday and Ampere. 2 depicts the position of the electric (yellow) and magnetic (green) field components for each primary and secondary Yee cell, respectively, allocated within the staggered cartesian grid. The FDTD algorithm proposed by Yee is based on a description of the temporarily coupled system as described in equations 2 and 3 on the basis of a finite central difference approximation. The exponentially increasing availability of computational power has made FDTD the most popular numerical method for a broad range of applications and has resulted in the release of several public and commercial FDTD-based simulation platforms.īy assuming a spatial environment without any electric or magnetic sources, the relationsĮnable a definition of the time dependent Maxwell’s curl equations in differential form as follows: 1, approximately 4500 scientific papers related to FDTD have been published, mainly from the late 90s to the present (1999-2001: counts not complete). Within the last two decades, FDTD gained rapidly increasing interest, mainly in electromagnetics (EM), for the simulation of complex and largely inhomogeneous structures due to its straightforward and explicit approach. The Finite-Difference Time-Domain (FDTD) method is based on a spatial and temporal discretization of Maxwell’s equations, commonly within a rectilinear cartesian grid originally proposed by Yee in 1966. Taflove, ed., Advances in Computational Electrodynamics: The Finite-Difference Time-Domain Method, Boston: Artech House, 1998.The Finite-Difference Time-Domain (FDTD) Technique Introduction to FDTD Potter, "FDTD Discrete Planewave (FDTD-DPW) Formulation for a Perfectly Matched Source in TFSF Simulations," IEEE Transactions on Antennas and Propagation, 58(8), 2010 pp. The program lets you plot the electromagnetic field components, and and also some intermediary fields. Since we have considered -periodicity, the parameter cannot be set to zero. Where the parameters and are multiples of and, respectively. The propagation angle is defined as the angle between the axis and the wave vector of the source and is given by the relation At, we set electric and magnetic currents in order to cancel the incoming plane wave. The periodicity of the simulation area is considered in the direction such that a periodic source is set in the plane defined by the equation. The parameters and denote the number of mesh points per unit on the and axes, respectively. The simulation area is defined by the and parameters. The resulting equations are updated using similar leapfrog algorithms, such as the finite-difference time-domain (FDTD) algorithm. This problem is avoided by introducing a phase shift and by splitting the obtained plane wave into two components. For an oblique incident plane wave, the periodic boundary condition needs some special attention, due to the fact that an advanced time field is needed. The study of periodic structures can often be reduced to the analysis of a single basic cell.