Finite-difference time-domain analysis of gyrotropic media. I.Magnetized plasma

When subjected to a constant magnetic field, both plasmas and ferrites exhibit anisotropic constitutive parameters. For electronic plasmas this anisotropy must be described by using a permittivity tensor in place of the usual scalar permittivity. Each member of this tensor is also very frequency dependent. A finite-difference time-domain formulation which incorporates both anisotropy and frequency dispersion, enabling the wideband transient analysis of magnetoactive plasma, is described. Results are shown for the reflection and transmission through a magnetized plasma layer, with the direction of propagation parallel to the direction of the biasing field. A comparison to frequency-domain analytic results is included     

Finite-difference time-domain implementation of surface impedanceboundary conditions

Surface impedance boundary conditions can be utilized to avoid using small cells, made necessary by shorter wavelengths in conducting media throughout the solution volume. The standard approach is to approximate the surface impedance over a very small bandwidth by its value at the center frequency, and then use that result in the boundary condition. In this paper, two implementations of the surface impedance boundary condition are presented. One implementation is a constant surface impedance boundary condition and the other is a dispersive surface impedance boundary condition that is applicable over a very large frequency bandwidth and over a large range of conductivities. Frequency domain results are presented in one dimension for two conductivity values and are compared with exact results. Scattering width results from an infinite square cylinder are presented as a two dimensional demonstration     

Finite-difference time-domain simulation of ground penetratingradar on dispersive, inhomogeneous, and conductive soils

A three-dimensional (3D) time-domain numerical scheme for simulation of ground penetrating radar (GPR) on dispersive and inhomogeneous soils with conductive loss is described. The finite-difference time-domain (FDTD) method is used to discretize the partial differential equations for time stepping of the electromagnetic fields. The soil dispersion is modeled by multiterm Lorentz and/or Debye models and incorporated into the FDTD scheme by using the piecewise-linear recursive convolution (PLRC) technique. The dispersive soil parameters are obtained by fitting the model to reported experimental data. The perfectly matched layer (PML) is extended to match dispersive media and used as an absorbing boundary condition to simulate an open space. Examples are given to verify the numerical solution and demonstrate its applications. The 3D PML-PLRC-FDTD formulation facilitates the parallelization of the code. A version of the code is written for a 32-processor system, and an almost linear speedup is observed     

Finite-difference time-domain analysis of a complete transverseelectromagnetic cell loaded with liquid biological media in culturedishes

Transverse electromagnetic (TEM) cells can be used for exposing biological culture specimens to electromagnetic fields and observing possible anomalous effects. The uniformity of field exposure is critical to quantifying the biological response versus the electromagnetic dose. Standing waves and other electromagnetic field nonuniformities can cause nonuniform exposure. This paper reports the results of high-resolution three-dimensional finite-difference time-domain (FDTD) simulations of a complete TEM cell designed for operation at 837 MHz. Several different cases were studied in which the number of culture dishes, the depth of the culture liquid, and the orientation of the culture dishes were varied. Further, the effect of the culture-dish glass bottom thickness and the meniscus of the liquid medium were examined. The FDTD results show that there is a significant nonuniform field and specific absorption rate (SAR) distribution within the culture medium for each case examined. Hence, biological dose-response experiments using the TEM cell should account for the possibility of strong localized SAR peaking in the culture media to provide useful data in setting exposure standards for wireless communications     

电离层时域散射与逆散射的数值

采用具有四阶精度的时域有限差分法计算了电离层对时域脉冲的散射,然后应用小波变换求出了反射信号的时频分布.这一时频分布就相当于电离层的频高图.最后利用POLAN程序重建了等离子体剖面,重建剖面非常接近原始剖面.     

Finite-difference time-domain analysis of HF antennas on helicopterairframes

In this paper, the finite-difference time-domain (FDTD) method with the Berenger perfectly matched layer (PML) absorbing boundary condition (ABC) is used to model the radiation characteristics of high frequency (HF) antennas operating in the 2-30 MHz range on a full-scale helicopter. The computed input impedance of both antennas is compared with actual measurements from an operational full-scale helicopter and also with measurements on a scale model NASA generic advanced attack helicopter (GAAH). To study the coupling effects of the helicopter fuselage on the antenna systems, the S-parameters are computed and compared with measurements on the NASA GAAH scale model. Finally, computed gain patterns are compared with actual in-flight measurements of the antenna systems on an operational full-scale helicopter     

Finite-difference time-domain analysis of flip-chip interconnectswith staggered bumps

This paper presents finite-difference time-domain (FDTD) analysis of flip-chip interconnects. Transitions between coplanar waveguides on the chip and the mother board are investigated over a broad band of frequency by means of Fourier transform of the time-domain results. Objectives of the analysis include the evaluation of bump reflection and insertion loss as well as the reconfiguration of the transition to improve package performance. Novel designs have been developed and presented to reduce the effects of package discontinuities and asymmetry. Staggering the bumps has been found to reduce reflection and insertion loss over a broad band of frequency. A reduction In bump reflection of up to 8 dB per transition can be achieved by staggering the ground and signal connects. The degradation in package performance due to structure asymmetry is also studied. The present designs have also been found to reduce the effects of flip-chip asymmetry on insertion and reflection losses     

A numerical simulation of scattering from one-dimensionalinhomogeneous dielectric random surfaces

An efficient numerical solution for the scattering problem of inhomogeneous dielectric rough surfaces is presented. The inhomogeneous dielectric random surface represents a bare soil surface and is considered to be comprised of a large number of randomly positioned dielectric humps of different sizes, shapes, and dielectric constants above an impedance surface. Clods with nonuniform moisture content and rocks are modeled by inhomogeneous dielectric humps and the underlying smooth wet soil surface is modeled by an impedance surface. In this technique, an efficient numerical solution for the constituent dielectric humps over an impedance surface is obtained using Green's function derived by the exact image theory in conjunction with the method of moments. The scattered field from a sample of the rough surface is obtained by summing the scattered fields from all the individual humps of the surface coherently ignoring the effect of multiple scattering between the humps. The statistical behavior of the scattering coefficient σ° is obtained from the calculation of scattered fields of many different realizations of the surface. Numerical results are presented for several different roughnesses and dielectric constants of the random surfaces. The numerical technique is verified by comparing the numerical solution with the solution based on the small perturbation method and the physical optics model for homogeneous rough surfaces. This technique can be used to study the behavior of scattering coefficient and phase difference statistics of rough soil surfaces for which no analytical solution exists     

Finite-difference, time-domain analysis of lossy transmission lines

An active and efficient method of including frequency-dependent conductor losses into the time-domain solution of the multiconductor transmission line equations is presented. It is shown that the usual A+B√s representation of these frequency-dependent losses is not valid for some practical geometries. The reason for this the representation of the internal inductance the at lower frequencies. A computationally efficient method for improving this representation in the finite-difference time-domain (FDTD) solution method is given and is verified using the conventional time-domain to frequency-domain (TDFD) solution technique     

Finite-difference time-domain analysis of the Celestron-8 telescope

Electromagnetic interference analyses of large complex systems demand large computational resources and give limited information on general types of systems. A finite-difference time-domain (FDTD) code was used to determine the response of a “generic” optical system to microwave radiation. A plane wave with a Gaussian pulse excitation was used along with “point sensors” within the system model to determine time and frequency response. In the low-frequency region, ramped sinusoidal excitation from a point within the sensor was used to determine angles of high sensitivity and field distributions within the sensor. From these field distributions, resonance modes were identified that are similar to those found in a simple cylindrical cavity     

NUMERICAL SIMULATION OF TIME-DOMAIN EM SCATTERING BY SHALLOW SUBSURFACE OBJECTS

The time-domain ElectroMagnetic(EM) scattering by buried objects in dispersive media is calculated with FD-TD method. The FD-TD formula in Debye dispersive media (both the complex permeability and the complex permittivity are described by Debye equations) are deduced, and the absorbing boundary condition is given. The validity of FD-TD method in lossy media is verified through comparing the FD-TD's results and the other ones. The propagation of transient pulses in dispersive media is studied in detail. The scattering pulses and the wiggle traces for typical buried objects are given.     

Finite-difference time-domain analysis of laser diodes integratedwith tapered beam-expanders

Waveguide properties of laser diodes integrated with horizontally tapered beam-expanders are analyzed by the finite-difference time-domain method. The taper length dependence of the radiation loss and fiber-coupling efficiency are clarified. Lower loss and higher fiber-coupling efficiency are achieved in the exponentially tapered beam-expander compared with a linearly tapered one having an equivalent length     

Finite-difference time-domain solution of Maxwell's equations forthe dispersive ionosphere

The application of the finite-difference time-domain (FDTD) technique to problems in ionospheric radio wave propagation is complicated by the dispersive nature of the ionospheric plasma. In the time domain, the electric displacement is the convolution of the dielectric tensor with the electric field, and thus requires information from the entire signal history. It is shown that this difficulty can be avoided by returning to the dynamical equations from which the dielectric tensor is derived. By integrating these differential equations simultaneously with the Maxwell equations, temporal dispersion is fully incorporated. An FDTD approach utilizing the vector wave equation is also presented. The accuracy of the method is shown by comparison for a special case for which an analytic solution is available. The method is demonstrated with examples of pulse propagation in one and two dimensions. The computational limitations of present-generation computers are discussed. The application of this approach to the study of wave propagation in randomly structured ionization is addressed     

Finite-difference time-domain (FDTD) analysis of magnetic diffusion

Problems with very slow waveforms or very long diffusion times may be difficult to treat using finite-difference time-domain techniques because of the Courant stability condition. Problems of this class, however, often prove to have a response which does not depend on the speed of light. The examples presented here show cases where internal fields do not change if c is reduced by as much as five orders of magnitude. This permits Δt to be proportionally increased. For simplicity much of this paper is restricted to one dimension, although generalization to three dimensions is also presented. The author considers an aluminum enclosure. Initially, the transient field will induce eddy currents on the enclosure which exactly cancel the external field and exclude it from the enclosure interior. This scheme has been, in fact, proposed to shield large systems which contain magnetic memories     

Finite-difference time-domain modeling of curved surfaces [EMscattering]

The finite-difference-time-domain (FDTD) method is generalized to include the accurate modeling of curved surfaces. This generalization, the contour path CP), method, accurately models the illumination of bodies with curved surfaces, yet retains the ability to model corners and edges. CP modeling of two-dimensional electromagnetic wave scattering from objects of various shapes and compositions is presented     

Multidomain pseudospectral time-domain simulations of scattering by objects buried in lossy media

A multidomain pseudospectral time-domain (PSTD) method with a newly developed well-posed PML is introduced as an accurate and flexible tool for the modeling of electromagnetic scattering by 2-D objects buried in an inhomogeneous lossy medium. Compared with the previous single-domain Fourier PSTD method, this approach allows for an accurate treatment of curved geometries with subdomains, curvilinear mapping, and high-order Chebyshev polynomials. The effectiveness of the algorithm is confirmed by an excellent agreement between the numerical results and analytical solutions for perfectly conducting as well as permeable dielectric cylinders. The algorithm has been applied to model various ground-penetrating radar (GPR) applications involving curved objects in a lossy half space with an undulating surface. This multidomain PSTD algorithm is potentially a very useful tool for simulating antennas near complex objects and inhomogeneous media.     

基于散射中心模型的目标时域回波快速仿真

为快速获取超宽带电磁脉冲激励下雷达目标在时域上的电磁响应,提出一种基于属性散射中心正向建模的方法用于目标时域回波仿真。从目标几何模型出发,利用射线追踪、分集技术对空间中所有射线进行标记与归类,分离并定量表征目标的强散射源。基于属性散射中心模型,正向确定模型参数,构建出目标属性散射中心模型,在选定的辐射源激励下进行仿真计算,快速获取目标时域回波信号。以典型目标简化坦克为例,选取不同形式的超宽带电磁脉冲信号作为辐射源,基于正向建模方法构建简化坦克的散射中心模型,快速获取给定电磁脉冲激励下的雷达回波信号,并与利用高频仿真方法得到的一维距离像对比,结果具有较好的一致性,从而验证了利用散射中心模型快速进行不同辐射源激励下回波仿真的有效性。     

Discontinuous galerkin time-domain method for GPR simulation in dispersive media

This paper presents a newly developed high-order discontinuous Galerkin time-domain (DGTD) method for solving Maxwell's equations in linear dispersive media with UPML boundary treatment. A unified formulation is derived for linear dispersive media of Debye type and the artificial material in the UPML regions with the help of auxiliary differential equations. The DGTD employs finite-element-type meshes, and uses piecewise high-order polynomials for spatial discretization and Runge-Kutta method for time integrations. Arbitrary high-order accuracy can be obtained for scattering of various objects in dispersive media. After validating the numerical convergence of the DGTD method together with the second-order Yee's scheme, we apply this new method to the ground-penetrating radar for the detection of buried objects in a lossy half space.     

基于FDTD的罗兰-C信号地波传播特性的时域分析

罗兰-C信号为载波调制的高斯脉冲,而传统方法对其传播特性预测均是基于100 kHz单频信号结果。文章采用FDTD方法计算了实际罗兰信号的时域特性,并和100 kHz连续单频信号结果进行了比较。结果显示:对于地形起伏不大的传播路径,两种信号结果吻合得很好,而对于地形起伏较大路径,脉冲信号结果与单频信号结果存在较大差异。该结果对地波传播特性工程测量与分析具有一定指导意义。     

Finite-difference time-domain calculation of the spontaneousemission coupling factor in optical microcavities

We present a general method for the β factor calculation in optical microcavities. The analysis is based on the classical model for atomic transitions in a semiconductor active medium. The finite-difference time-domain method is used to evolve the electromagnetic fields of the system and calculate the total radiated energy, as well as the energy radiated into the mode of interest. We analyze the microdisk laser and compare our result with the previous theoretical and experimental analyses. We also calculate the β factor of the microcavity based on a two-dimensional (2-D) photonic crystal in an optically thin dielectric slab. From the β calculations, we are able to estimate the coupling to radiation modes in both the microdisk and the 2-D photonic crystal cavity, thereby showing the effectiveness of the photonic crystal in suppressing in-plane radiation modes