Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FDTD methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

An Implicit Characteristic Based Method for Electromagnetics

An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

Beggs, John H.↗
NASA NTRS2005

Inference of Ice Cloud Properties from High-spectral Resolution Infrared Observations

The theoretical basis is explored for inferring the microphysical properties of ice crystal from high-spectral resolution infrared observations. A radiative transfer model is employed to simulate spectral radiances to address relevant issues. The extinction and absorption efficiencies of individual ice crystals, assumed as hexagonal columns for large particles and droxtals for small particles, are computed from a combination of the finite- difference time-domain (FDTD) technique and a composite method. The corresponding phase functions are computed from a combination of FDTD and an improved geometric optics method (IGOM). Bulk scattering properties are derived by averaging the single- scattering properties of individual particles for 30 particle size distributions developed from in situ measurements and for additional four analytical Gamma size distributions for small particles. The non-sphericity of ice crystals is shown to have a significant impact on the radiative signatures in the infrared (IR) spectrum; the spherical particle approximation for inferring ice cloud properties may result in an overest&ation of the optical thickness and an inaccurate retrieval of effective particle size. Furthermore, we show that the error associated with the use of the Henyey-Greenstein phase function can be as larger as 1 K in terms of brightness temperature for larger particle effective size at some strong scattering wavenumbers. For small particles, the difference between the two phase functions is much less, with brightness temperatures generally differing by less than 0.4 K. The simulations undertaken in this study show that the slope of the IR brightness temperature spectrum between 790-960/cm is sensitive to the effective particle size. Furthermore, a strong sensitivity of IR brightness temperature to cloud optical thickness is noted within the l050-1250/cm region. Based on this spectral feature, a technique is presented for the simultaneous retrieval of the visible optical thickness and effective particle size from high spectral resolution infrared data under ice cloudy con&tion. The error analysis shows that the uncertainty of the retrieved optical thickness and effective particle size has a small range of variation. The error for retrieving particle size in conjunction with an uncertainty of 5 K in cloud'temperature, or a surface temperature uncertainty of 2.5 K, is less than 15%. The corresponding e m r in the uncertainty of optical thickness is within 5-2096, depending on the value of cloud optical thickness. The applicability of the technique is demonstrated using the aircraft-based High- resolution Interferometer Sounder (HIS) data from the Subsonic Aircraft: Contrail and Cloud Effects Special Study (SUCCESS) in 1996 and the First ISCCP Regional Experiment - Arctic Clouds Experiment (FIRE-ACE) in 1998.

Huang, Hung-Lung↗
NASA NTRS1993

Hybrid techniques for complex aerospace electromagnetics problems

Important aerospace electromagnetics problems include the evaluation of antenna performance on aircraft and the prediction and control of the aircraft's electromagnetic signature. Due to the ever increasing complexity and expense of aircraft design, aerospace engineers have become increasingly dependent on computer solutions. Traditionally, computational electromagnetics (CEM) has relied primarily on four disparate techniques: the method of moments (MoM), the finite-difference time-domain (FDTD) technique, the finite element method (FEM), and high frequency asymptotic techniques (HFAT) such as ray tracing. Each of these techniques has distinct advantages and disadvantages, and no single technique is capable of accurately solving all problems of interest on computers that are available now or will be available in the foreseeable future. As a result, new approaches that overcome the deficiencies of traditional techniques are beginning to attract a great deal of interest in the CEM community. Among these new approaches are hybrid methods which combine two or more of these techniques into a coherent model. During the ASEE Summer Faculty Fellowship Program a hybrid FEM/MoM computer code was developed and applied to a geometry containing features found on many modern aircraft.

Aberle, Jim↗
NASA NTRS1991

Finite difference time domain calculations of antenna mutual coupling

The Finite Difference Time Domain (FDTD) technique was applied to a wide variety of electromagnetic analysis problems, including shielding and scattering. However, the method has not been exclusively applied to antennas. Here, calculations of self and mutual admittances between wire antennas are made using FDTD and compared with results obtained during the method of moments. The agreement is quite good, indicating the possibilities for FDTD application to antenna impedance and coupling.

Luebbers, Raymond J.↗
NASA NTRS1991

Finite difference time domain calculations of antenna mutual coupling

The Finite Difference Time Domain (FDTD) technique has been applied to a wide variety of electromagnetic analysis problems, including shielding and scattering. However, the method has not been extensively applied to antennas. In this short paper calculations of self and mutual admittances between wire antennas are made using FDTD and compared with results obtained using the Method of Moments. The agreement is quite good, indicating the possibilities for FDTD application to antenna impedance and coupling.

Luebbers, Raymond J.↗
NASA NTRS1994

Theoretical and Experimental Study of Microstrip-to-Slot Line Uniplanar Transition

Recent advances in MMCI technology make it possible to construct transitions from CPW-to-microstrip with via hole, microstrip-to-slot line and microshield line-to-CPW all of which have potential applications in the feed network of antennas. In this study we investigate the characteristics of the microstrip-to-slot line uniplanar transition using the finite element methods (FEM) and finite difference time domain (FDTD) techniques, and compared the theoretical results with the measurements. In both cases, the results agree with the measurements within a few percent.

Yook, Jong-Gwan↗
NASA NTRS2005

Characterization of Meta-Materials Using Computational Electromagnetic Methods

An efficient and powerful computational method is presented to synthesize a meta-material to specified electromagnetic properties. Using the periodicity of meta-materials, the Finite Element Methodology (FEM) is developed to estimate the reflection and transmission through the meta-material structure for a normal plane wave incidence. For efficient computations of the reflection and transmission over a wide band frequency range through a meta-material a Finite Difference Time Domain (FDTD) approach is also developed. Using the Nicholson-Ross method and the Genetic Algorithms, a robust procedure to extract electromagnetic properties of meta-material from the knowledge of its reflection and transmission coefficients is described. Few numerical examples are also presented to validate the present approach.

Deshpande, Manohar↗
NASA NTRS1993

Numerical modeling for an electric-field hyperthermia applicator

Hyperthermia, in conjunction with radiation and chemotherapy for treatment of cancers, is an area of current concern. Experiments have shown that hyperthermia can increase the potency of many chemotherapy drugs and the effectiveness of radiation for treating cancer. A combination of whole body or regional hyperthermia with chemotherapy or radiation should improve treatment results. Conventional methods for inducing whole body hyperthermia, such as exposing a patient in a radiant cabinet or under a hot water blanket, conduct heat very slowly from the skin to the body core. Thus a more efficient system, such as the three-plate electric-field hyperthermia applicator (EHA), is developed. This three-plate EHA has one top plate over and two lower plates beneath the patient. It is driven at 27.12 MHz with 500 Watts through a matching circuit. Using this applicator, a 50 kg pig was successfully heated to 42 C within 45 minutes. However, phantom and animal studies have indicated non-uniform heating near the side of the body. In addition, changes in the size and distance between the electrode plates can affect the heating (or electromagnetic field) pattern. Therefore, numerical models using the method of moments (MOM) or the finite difference time domain (FDTD) technique are developed to optimize the heating pattern of this EHA before it is used for human trials. The accuracy of the numerical modeling has been achieved by the good agreement between the MOM and FDTD results for the three-plate EHA without a biological body. The versatile FDTD technique is then applied to optimize the EHA design with a human body. Both the numerical and measured data in phantom blocks will be presented. The results of this study will be used to design an optimized system for whole body or regional hyperthermia.

Wu, Te-Kao↗
NASA NTRS1991

Conical cut radar cross section calculations for a thin, perfectly conducting plate

Radar Cross Section (RCS) calculations for flat, perfectly conducting plates are readily available through the use of conventional frequency domain techniques such as the Method of Moments. However, if time domain scattering or wideband frequency domain results are desired, then the Finite Difference Time Domain (FDTD) technique is a suitable choice. In this paper, we present the application of the Finite Difference Time Domain (FDTD) technique to the problem of electromagnetic scattering and RCS calculations from a thin, perfectly conducting plate for a conical cut in the scattering angle phi. RCS calculations versus angle phi will be presented and discussed.

Luebbers, Raymond J.↗
NASA NTRS1991

Conical cut radar cross section calculations for a thin, perfectly conducting plate

Radar cross section (RCS) calculations for flat, perfectly conducting plates are readily available through the use of conventional frequency domain techniques such as the method of moments. However, if time domain scattering or wideband frequency domain results are desired, then the finite difference time domain (FDTD) technique is a suitable choice. We present the application of the FDTD technique to the problem of electromagnetic scattering and RCS calculations from a thin, perfectly conducting plate for a conical cut in the scattering angle phi. RCS calculations versus angle phi are presented and discussed.

Luebbers, Raymond J.↗
NASA NTRS2012

GPR Imaging for Deeply Buried Objects: A Comparative Study Based on FDTD Models and Field Experiments

Conventional use of Ground Penetrating Radar (GPR) is hampered by variations in background environmental conditions, such as water content in soil, resulting in poor repeatability of results over long periods of time when the radar pulse characteristics are kept the same. Target objects types might include voids, tunnels, unexploded ordinance, etc. The long-term objective of this work is to develop methods that would extend the use of GPR under various environmental and soil conditions provided an optimal set of radar parameters (such as frequency, bandwidth, and sensor configuration) are adaptively employed based on the ground conditions. Towards that objective, developing Finite Difference Time Domain (FDTD) GPR models, verified by experimental results, would allow us to develop analytical and experimental techniques to control radar parameters to obtain consistent GPR images with changing ground conditions. Reported here is an attempt at developing 20 and 3D FDTD models of buried targets verified by two different radar systems capable of operating over different soil conditions. Experimental radar data employed were from a custom designed high-frequency (200 MHz) multi-static sensor platform capable of producing 3-D images, and longer wavelength (25 MHz) COTS radar (Pulse EKKO 100) capable of producing 2-D images. Our results indicate different types of radar can produce consistent images.

Tilley, roger↗
NASA NTRS1993

A physical optics/equivalent currents model for the RCS of trihedral corner reflectors

The scattering in the interior regions of both square and triangular trihedral corner reflectors is examined. The theoretical model presented combines geometrical and physical optics (GO and PO), used to account for reflection terms, with equivalent edge currents (EEC), used to account for first-order diffractions from the edges. First-order, second-order, and third-order reflection terms are included. Calculating the first-order reflection terms involves integrating over the entire surface of the illuminated plate. Calculating the second- and third-order reflection terms, however, is much more difficult because the illuminated area is an arbitrary polygon whose shape is dependent upon the incident angles. The method for determining the area of integration is detailed. Extensive comparisons between the high-frequency model, Finite-Difference Time-Domain (FDTD) and experimental data are used for validation of the radar cross section (RCS) of both square and triangular trihedral reflectors.

Balanis, Constantine A.↗
NASA NTRS1995

Development of Fast Algorithms Using Recursion, Nesting and Iterations for Computational Electromagnetics

In the first phase of our work, we have concentrated on laying the foundation to develop fast algorithms, including the use of recursive structure like the recursive aggregate interaction matrix algorithm (RAIMA), the nested equivalence principle algorithm (NEPAL), the ray-propagation fast multipole algorithm (RPFMA), and the multi-level fast multipole algorithm (MLFMA). We have also investigated the use of curvilinear patches to build a basic method of moments code where these acceleration techniques can be used later. In the second phase, which is mainly reported on here, we have concentrated on implementing three-dimensional NEPAL on a massively parallel machine, the Connection Machine CM-5, and have been able to obtain some 3D scattering results. In order to understand the parallelization of codes on the Connection Machine, we have also studied the parallelization of 3D finite-difference time-domain (FDTD) code with PML material absorbing boundary condition (ABC). We found that simple algorithms like the FDTD with material ABC can be parallelized very well allowing us to solve within a minute a problem of over a million nodes. In addition, we have studied the use of the fast multipole method and the ray-propagation fast multipole algorithm to expedite matrix-vector multiplication in a conjugate-gradient solution to integral equations of scattering. We find that these methods are faster than LU decomposition for one incident angle, but are slower than LU decomposition when many incident angles are needed as in the monostatic RCS calculations.

Chew, W. C.↗
NASA NTRS1991

Time domain scattering and radar cross section calculations for a thin, coated perfectly conducting plate

Radar cross section (RCS) calculations for flat, perfectly conducting plates are readily available through the use of conventional frequency domain techniques such as the Method of Moments (MOM). However, if the plate is covered with a dielectric material that is relatively thick in comparison with the wavelength in the material, these frequency domain techniques become increasingly difficult to apply. We present the application of the Finite Difference Time Domain (FDTD) Technique to the problem of electromagnetic scattering and RCS calculations from a thin, perfectly conducting plate that is coated with a thick layer of lossless dielectric material. Both time domain and RCS calculations are presented and disclosed.

Luebbers, Raymond J.↗
NASA NTRS1991

Time domain scattering and radar cross section calculations for a thin, coated perfectly conducting plate

Radar cross section (RCS) calculations for flat, perfectly conducting plates are readily available through the use of conventional frequency domain techniques such as the Method of Moments (MOM). However, if the plate is covered with a dielectric material that is relatively thick in comparison with the wavelength in the material, these frequency domain techniques become increasingly difficult to apply. The application is presented of the Finite Difference Time Domain (FDTD) technique to the problem of electromagnetic scattering and RCS calculations from a thin, perfectly conducting plate that is coated with a thick layer of lossless dielectric material. Both time domain and RCS calculations are presented and discussed.

Luebbers, Raymond J.↗
NASA NTRS1993

High-frequency techniques for RCS prediction of plate geometries and a physical optics/equivalent currents model for the RCS of trihedral corner reflectors

Part 1 of this report continues the investigation, initiated in previous reports, of scattering from rectangular plates coated with lossy dielectrics. The hard polarization coefficients given in the last report are incorporated into a model, which includes second- and third-order diffractions, for the coated plate. Computed results from this model are examined and compared to measured data. A breakdown of the contribution of each of the higher-order terms to the total radar cross section (RCS) is given. The effectiveness of the uniform theory of diffraction (UTD) model in accounting for the coating effect is investigated by examining a Physical Optics (PO) model which incorporates the equivalent surface impedance approximation used in the UTD model. The PO, UTD, and experimental results are compared. Part 2 of this report presents a RCS model, based on PO and the Method of Equivalent Currents (MEC), for a trihedral corner reflector. PO is used to account for the reflected fields, while MEC is used for the diffracted fields. Single, double, and triple reflections and first-order diffractions are included in the model. A detailed derivation of the E(sub theta)-polarization, monostatic RCS is included. Computed results are compared with finite-difference time-domain (FDTD) results for validation. The PO/MEC model of this report compares very well with the FDTD model, and it is a much faster model in terms of computational speed.

Balanis, Constantine A.↗
NASA NTRS2003

Diffraction Efficiency of Thin Film Holographic Beam Steering Devices

Dynamic holography has been demonstrated as a method for correcting aberrations in space deployable optics, and can also be used to achieve high-resolution beam steering in the same environment. In this paper, we consider some of the factors affecting the efficiency of these devices. Specifically, the effect on the efficiency of a highly collimated beam from the number of discrete phase steps per period is considered for a blazed thin film beam steering grating. The effect of the number of discrete phase steps per period on steering resolution is also considered. We also present some result of Finite-Difference Time-Domain (FDTD) calculations of light propagating through liquid crystal "blazed" gratings. Liquid crystal gratings are shown to spatially modulate both the phase and amplitude of the propagating light.

Titus, Charles M.↗