Finite-difference Time-domain FDTD Research Articles

Tunable and anisotropic perfect absorber using graphene-black phosphorus nanoblock.

Two-dimensional (2D) materials, which have attracted attention due to intriguing optical properties, form a promising building block in optical and photonic devices. This paper numerically investigates a tunable and anisotropic perfect absorber in a graphene-black phosphorus (BP) nanoblock array structure. The suggested structure exhibits polarization-dependent anisotropic absorption in the mid-infrared, with maximum absorption of 99.73% for x-polarization and 53.47% for y-polarization, as determined by finite-difference time-domain FDTD analysis. Moreover, geometrical parameters and graphene and BP doping amounts are possibly employed to tailor the absorption spectra of the structures. Hence, our results have the potential in the design of polarization-selective and tunable high-performance devices in the mid-infrared, such as polarizers, modulators, and photodetectors.

Grounding systems under lightning surges with soil ionization for high voltage substations by using two layer capacitors (TLC) model

The grounding system behavior in the time domain for the high voltage substations is the key-parameter to ensure the system reliability during transient conditions. Many approaches are used to investigate the behavior of the grounding systems subjected to the lightning surges. These approaches aim to simplify the grounding systems by equivalent lumped circuit-parameters (RLC) or by using the electromagnetic field theory. The soil ionization phenomenon can effectively reduce the grounding electrode resistance because of changing the electrode equivalent radius located in the soils. Yet, most of these approaches have their relative assumptions with merits and demerits. The current work introduces a proposed model based on a developed grounding unit cell using the Two Layer Capacitors (TLC). The proposed cell is implemented as a unit cell to simulate the performance of the grounding system under lightning strikes with considering the soil ionization effects. Different developed models of the grounding-cells are used to simulate the behavior of the grounding system including the soil ionization process. Finite Difference Time Domain FDTD scheme is used to calculate the electromagnetic fields inside the different proposed unit cells during lightning surges. The application of the grounding unit cell is extended to be implemented for a typical high voltage substation of 500/220 kV. The grounding unit cells are presented as T neighboring sections to evaluate the overlapping between the different elementary electrodes and conductors of the complex grounding systems over the whole area of the substation. The circuit-model of the grounding unit cell is developed using EMTP-RV/MATLAB. Finally, the obtained results are compared with the similar results found in the literature with good consistency.

A non-conforming multi-element DGTD method for the simulation of human exposure to electromagnetic waves

The great majority of numerical calculations of electric energy deposition in human tissues exposed to microwaves are performed using the finite-difference time-domain FDTD method and voxel-based geometrical models of the tissues. The straightforward implementation of the method and its computational efficiency are among the main reasons for FDTD being currently the leading method for numerical assessment of human exposure to electromagnetic waves. However, the rather difficult departure from the commonly used Cartesian grid and cell size limitations regarding the discretization of very detailed structures of human tissues are often recognized as the main weaknesses of the FDTD method in this application context. In particular, interfaces between tissues where sharp gradients of the electromagnetic field may occur are hardly modeled rigorously in these studies this is generally referred as the staircasing effect. We present here an alternative numerical dosimetry methodology that is based on a DGTD method designed to work with non-conforming cubic-tetrahedral meshes for more flexibility in the discretization process of complex propagation scenes. For example, in the context of head tissues exposure to mobile phone radiation, the computation domain is divided in two regions: an unstructured tetrahedral mesh-based heterogeneous model of head tissues and an orthogonal cartesian mesh of the vacuum part of the propagation domain. We present numerical results comparing this approach with strategies based on fully Cartesian and fully tetrahedral meshes of the whole computational domain. Copyright © 2013 John Wiley & Sons, Ltd.

Numerical assessment of the combination of subgridding and the perfectly matched layer grid termination in the finite difference time domain method

The finite-difference time-domain FDTD method is one of the well-established numerical techniques to model electromagnetic problems. In order to be able to efficiently include the finer geometrical details, subgridding is used to avoid the use of a fine mesh over the whole region of interest. Our work builds on and extends the systematic framework provided by Chilton in his PhD H-refinement, P-refinement, and T-refinement strategies for the FDTD method developed via finite element principles, Ohio State University, 2008.

Performance Optimization of Massively Parallel FDTD Computations

The advancements of General Purpose Graphic Processing Units GPGPUs, have paved the way for computationally intensive scientific calculations to be done on an off the shelf massively parallel graphic processors GPUs, rather than the use of expensive solutions such as High Performance Clusters HPCs or Supercomputers. In this paper the NVIDIA’s Compute Unified Device Architecture CUDA on an NVIDIA GeForce processor will be used to solve a Finite Difference Time Domain FDTD computation. An FDTD computation has unfeasible running time using a single processor, and attempts to use multi-core CPUs have been made, but with the high overhead of network traffic in HPCs or synchronizations among cores. Multiple attempts have been made to utilize GPU for solving FDTD computations; but in this paper the focus will be on optimizing the efficiency of the algorithm by maximizing the throughput through affective use of the fast on-chip shared memory, and avoid using the slow off-chip global memory.

Electromagnetic simulation of a synchronous impedance transformer via the hybrid FDTD-MoL method

The paper aims the electromagnetic simulation of a synchronous impedance transformer using a combination of the finite difference time domain FDTD and the method of lines MoL. The study demonstrates that the hybrid technique called FDTD-MoL can be used for the characterization of passive microstrip circuits. The capability to combine different methods of electromagnetic simulation provides a powerful tool for interpreting the behavior of microwave circuits. The hybrid FDTD-MoL method can characterize microstrip planar structures and is comparable with previous case studies of electromagnetic simulation.

Study on the stability of the higher-order FDTD (2, 4) including lumped inductors

The paper presents a study of the numerical stability of the higher-order finite-difference time domain FDTD (2, 4) technique for calculating the lumped-inductor-loaded region. Six different cases including explicit, implicit and semi-implicit schemes based on the differential and integral equations including the lumped inductor's current-voltage characteristics are analysed in detail. Five stability condition formulas and an unconditional instability conclusion are obtained using the Routh–Hurwitz criterion. The theoretical results are validated by the numerical results gained by means of actual FDTD (2, 4) simulations.

Domain decomposition FDTD algorithm for the analysis of a new type of E-plane sectorial horn with aperture field distribution optimization

In this paper, a universal and efficient approach of domain decomposition finite-difference time-domain FDTD (DD-FDTD) is presented for the analysis of a new type of horn antenna-E-plane sectorial horn with field amplitude taper and phase correction in the aperture. The power fed into the horn is redistributed to achieve the optimal field amplitude distribution in the aperture, and meanwhile the field phase is corrected by metal lens. Compared with conventional E-plane sectorial horns, the new horn antenna takes the advantages of low sidelobe level, short physical length and wide flare angle etc. Moreover, the most important property of this horn is the weak coupling with each other when it is used as the element of a phased array antenna. The field analysis of such a horn antenna is an extremely complicated three-dimensional EM boundary value problem. The domain decomposition FDTD method is presented in this paper to break through the drawback. The whole horn is decomposed into several subdomains and the meshes are created in local coordinates. In the iteration procedure of FDTD, the data are exchanged between adjacent subdomains with overlapped meshes. The aperture field distribution, voltage standing-wave ratio and pattern calculated by the DD-FDTD method are in good agreement with the experimental data.

A time-domain near- to far-field transformation for FDTD in two dimensions

This paper proposes an algorithm based on the equi alence principle to calculate, straightforwardly in the time domain, the transient far field response of a two-dimensional structure. The algorithm implementation in ol es the numerical e aluation of a time integral along the nearto far-field con ersion contour of the equi alent currents at the ( ) Huygens surface, obtained with the finite-difference time-domain FDTD method. 2000 John Wiley & Sons, Inc. Microwave Opt Technol Lett 27: 427 432, 2000.

A finite-difference time-domain electromagnetic field solver for vertical-cavity lasers

ABSTRACT: We present a full- ¤ ector Maxwell equation sol ¤ er for theanalysis of ¤ ertical-ca ¤ ity surface-emitting lasers VCSELs with circular()symmetry using the finite-difference time-domain FDTD method. We()apply this model to the air-post index-guided 1.55 m m VCSEL structure,and determine the resonant wa ¤ elengths and the quality factor forselected lowest order modes. The results of this study enable quantitati ¤ eand qualitati ¤ e descriptions of the trans erse-mode competition andemission spectra of VCSELs. Q 1998 John Wiley & Sons, Inc.Microwave Opt Technol Lett 18: 385]387, 1998. Key words: electromagnetic field sol ¤ er; finite-difference time-domainmethods; ¤ ertical-ca ity lasers Vertical-cavity surface-emitting lasers VCSELs demonstrate.many properties desirable for optical communications andinterconnects, such as low-beam divergence and single-modeoperation which result in superior coupling efficiency intosingle-mode fibers 1 , the possibility of wafer-level testing,wxand ease for two-dimensional array fabrication 2 . In manywxapplications, and for long-haul optical communications, inparticular, single-transverse-mode operation is desired. Forthis reason, a number of theoretical analyses have beenrecently reported in an attempt to understand the transverse-mode competition and the emission spectra of these lasers.The scalar beam-propagation 4 and the effective index 5wx wxmethods have been successfully used, but since they neglectthe vector nature of the electromagnetic field, these methodsare unable to model the full transverse-mode spectrum ofthese lasers. The transverse modes of vertical-cavity lasersare rather complicated, while their spatial shape and polar-ization properties depend on the type of mode definitionscheme. Index-guided, gain-guided, and oxide-apertured lasersshow distinctly different transverse modes and their depen-dence on injection current 6, 7 . Because these cavities arewxsmall on the order of a few wavelengths in all directions , the.application of beam-propagation methods and effective-indexmethods is unable to provide a sufficiently accurate descrip-tion of the optical field in the cavity. Moreover, the beam-propagation method does not account for multiple reflectionsin quarter-wave mirrors. Particularly difficult structures tomodel are the index-guided air-post devices and the oxide-apertured devices in which the large lateral refractive indexvariations severely add to the complexity of the problem 8 .wxFor this reason, it is necessary to include the vector nature ofthe electromagnetic field to accurately analyze vertical-cavitylasers. In this paper, we present an analysis of the passivecavity of a 1.55 mm air-post index-guided VCSEL using afull-vector finite-difference time-domain FDTD method..The analyzed VCSEL structure is a 1550 nm double-fusedcavity that consists of an InGaAsP cavity sandwiched betweentwo AlGaAs AlGaAs.rGaAs quarter-wavelength distributedBragg reflectors DBR . A detailed description and the physi-.cal and optical structure parameters of this device are givenin 9 , while the simplified structure is shown in Figure 1 a .wx .The transverse modes of this resonator are determined byindex guiding in the top mirror, while the cavity loss occursdue to diffraction of the wave in the unguided sections of theactive layer and the bottom mirror. In a high-quality res-onator, such as the VCSEL, one may assume that the modesare not significantly perturbed by free-space propagation, andthe approximate cavity spectrum may be obtained by model-ing the VCSEL as a perfect cylindrical dielectric resonatorwx10 . However, this model does not provide the quality factorof each individual mode, namely, the discrimination betweenthe modes. For this purpose, one has to resort to a moreinvolved analysis: we use a full-vector finite-difference time-domain FDTD applicable to cylindrical high-quality open.resonators with dielectric andror metal boundaries. The pro-gram is used to determine the resonance wavelengths, thequality factor, and the mode spatial distribution for selectlowest order modes in the VCSEL cavity.The principle of operation of the FDTD program is basedon determining the passive-resonator temporal response tonarrowband noise delivered by an antenna source within thecavity. The resonator temporal response, detected at several