Incorporating non-linear lumped elements in FDTD: the equivalent source method

Abstract

A new approach is presented for modelling three‐dimensional lumped elements in the finite‐difference time‐domain (FDTD) solution of Maxwell's equations. The finite‐difference equations for the lumped element's circuit behaviour are derived to produce discrete relationships between the device's terminal currents and voltages. These difference equations are then implemented in the FDTD grid with equivalent voltage and current sources based on static field approximations. The method can be used for a wide range of applications which include lumped elements that may be passive, active, linear, non‐linear, single, and multiple port devices. The underlying advantages of modelling a lumped element in this manner as compared to the more traditional extended FDTD technique is that the lumped element is easily extended across multiple FDTD cells and the need for a separate transcendental equation solver is eliminated. The method is demonstrated by a specific example based on the Curtice–Cubic non‐linear model for a MESFET. The s‐parameters from an FDTD simulation are compared to the manufacturer's measured data and to results from simulations using Hewlett‐Packard's Microwave Design System (MDS). Finally, the method is used to model the MESFET in a three port oscillator. Copyright © 1999 John Wiley & Sons, Ltd.

Implementing voltage controlled current source in electromagnetic full-wave simulation using the FDTD method

“…To the best of the authors' knowledge, the implementation of dependent sources using FDTD has not been adequately addressed before. In addition, in most of the previous work the implementation of nonlinear devices such as transistors has been handled using FDTD-SPICE models or by importing the S-parameters from another technique to the FDTD simulation [2]- [7]. In this work the implementation of the VCCS in FDTD will be used to simulate a MOSFET with its equivalent model without the use of external tools, the entire simulation can be done using the FDTD.…”

“…4. The values of the equivalent circuit parameters according to [7] are as follows: C gd = 0.06pF, C gs = 0.96pF, C ds = 0.26pF, g m =10 mS and R ds = 197 Ω. The FDTD grid cell size is ∆x = ∆y = ∆z = 0.25 mm, the circuit is simulated with 5×1×2 in x, y and z-directions respectively, 1000 time steps are used, CPML is used as the absorbing boundaries of the computational domain with 5 cells air buffer in every direction.…”

Computation of Complex Eigenmodes for Resonators Filled With Gyrotropic Materials

“…In general, lumped elements can be incorporated into the FIT by formally adding a supplementary current j L due to the lumped element in Ampère's circuital law (3.16) [62,63]:…”

Combined Electromagnetic and Drift Diffusion Models for Microwave Semiconductor Device

“…In order to derive the iteration procedure, at each time step we relate Gummel's method to direct UMFPACK method which is well known to converge quadratically [7]. The Equations (10), (11) and (12) are solved in a decoupled manner. Poisson equation is solved at all grid points, followed by electron continuity equation, and then by hole continuity one.…”

“…These equations are first-order linear the field vectors at any point in the space at any time. The current density distribution is obtained from the solution of the driftdiffusion model [11,12].…”