Search NASA⌕ Search

SEARCH · Search NASA

Results for “finite-difference time-domain method”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Analysis of microstrip dipoles and slots transversely coupled to a microstrip line using the FDTD method

Printed dipole elements and their complement, linear slots, are elementary radiators that have found use in low-profile antenna arrays. Low-profile antenna arrays, in addition to their small size and low weight characteristics, offer the potential advantage of low-cost, high-volume production with easy integration with active integrated circuit components. The design of such arrays requires that the radiation and impedance characteristics of the radiating elements be known. The FDTD (Finite-Difference Time-Domain) method is a general, straight-forward implementation of Maxwell's equations and offers a relatively simple way of analyzing both printed dipole and slot elements. Investigated in this work is the application of the FDTD method to the analysis of printed dipole and slot elements transversely coupled to an infinite transmission line in a multilayered configuration. Such dipole and slot elements may be used in dipole and slot series-fed-type linear arrays, where element offsets and interelement line lengths are used to obtain the desired amplitude distribution and beam direction, respectively. The design of such arrays is achieved using transmission line theory with equivalent circuit models for the radiating elements. In an equivalent circuit model, the dipole represents a shunt impedance to the transmission line, where the impedance is a function of dipole offset, length, and width. Similarly, the slot represents a series impedance to the transmission line. The FDTD method is applied to single dipole and slot elements transversely coupled to an infinite microstrip line using a fixed rectangular grid with Mur's second order absorbing boundary conditions. Frequency-dependent circuit and scattering parameters are obtained by saving desired time-domain quantities and using the Fourier transform. A Gaussian pulse excitation is applied to the microstrip transmission line, where the resulting reflected signal due to the presence of the radiating element is used to determine the equivalent element impedance.

Tulintseff, A. N.↗
NASA NTRS1987

The finite-difference time-domain (FD-TD) method for electromagnetic scattering and interaction problems

The formulation and recent applications of the finite-difference time-domain (FD-TD) method for the numerical modeling of electromagnetic scattering and interaction problems are considered. It is shown that improvements in FD-TD modeling concepts and software implementation often make it a preferable choice for structures which cannot be easily treated by conventional integral equations and asymptotic approaches. Recent FD-TD modeling validations in research areas including coupling to wires and wire bundles in free space and cavities, scattering from surfaces in relativistic motion, inverse scattering, and radiation condition theory, are reviewed. Finally, the advantages and disadvantages of FD-TD, and guidelines concerning when FD-TD should and should not be used in high-frequency electromagnetic modeling problems, are summarized.

Taflove, A.↗
NASA NTRS1992

Finite difference time domain modeling of spiral antennas

The objectives outlined in the original proposal for this project were to create a well-documented computer analysis model based on the finite-difference, time-domain (FDTD) method that would be capable of computing antenna impedance, far-zone radiation patterns, and radar cross-section (RCS). The ability to model a variety of penetrable materials in addition to conductors is also desired. The spiral antennas under study by this project meet these requirements since they are constructed of slots cut into conducting surfaces which are backed by dielectric materials.

Penney, Christopher W.↗
NASA NTRS1992

Frequency-dependent FDTD methods using Z transforms

While the frequency-dependent finite-difference time-domain, or (FD)2TD, method can correctly calculate EM propagation through media whose dielectric properties are frequency-dependent, more elaborate applications lead to greater (FD)2TD complexity. Z-transform theory is presently used to develop the mathematical bases of the (FD)2TD method, simultaneously obtaining a clearer formulation and allowing researchers to draw on the existing literature of systems analysis and signal-processing.

Sullivan, Dennis M.↗
NASA NTRS1989

Review of FD-TD numerical modeling of electromagnetic wave scattering and radar cross section

Applications of the finite-difference time-domain (FD-TD) method for numerical modeling of electromagnetic wave interactions with structures are reviewed, concentrating on scattering and radar cross section (RCS). A number of two- and three-dimensional examples of FD-TD modeling of scattering and penetration are provided. The objects modeled range in nature from simple geometric shapes to extremely complex aerospace and biological systems. Rigorous analytical or experimental validatons are provided for the canonical shapes, and it is shown that FD-TD predictive data for near fields and RCS are in excellent agreement with the benchmark data. It is concluded that with continuing advances in FD-TD modeling theory for target features relevant to the RCS problems and in vector and concurrent supercomputer technology, it is likely that FD-TD numerical modeling will occupy an important place in RCS technology in the 1990s and beyond.

Taflove, Allen↗
NASA NTRS1993

FD-TD modeling of 2-D dielectric waveguides for propagation and scattering of femtosecond optical solitons

Experimentalists have produced all-optical switches capable of 100-fs responses. To adequately model such switches, nonlinear effects in optical materials (both instantaneous and dispersive) must be included. In principle, the behavior of electromagnetic fields in nonlinear dielectrics can be determined by solving Maxwell's equations subject to the assumption that the electric polarization has a nonlinear relation to the electric field. However, until our previous work, the resulting nonlinear Maxwell's equations have not been solved directly. Rather, approximations have been made that result in a class of generalized nonlinear Schrodinger equations (GNLSE) that solve only for the envelope of the optical pulses. In this paper, we present first-time calculations from the vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional systems of dielectric waveguides exhibiting the Kerr and Raman quantum effects. We use the finite-difference time-domain (FD-TD) method in an extension of our 1-D work. There, in a fundamental innovation, we treated the linear and nonlinear convolutions for the electric polarization as new dependent variables. By differentiating these convolutions in the time domain, we derived an equivalent system of coupled, nonlinear second-order ODE's. These equations together with Maxwell's equations form the system that is solved to determine the electromagnetic fields in inhomogeneous nonlinear dispersive media. Backstorage in time is limited to only that needed by the time-integration algorithm for the ODE's, rather than that needed to store the time-history of the kernel functions of the convolutions (1000-10,000 time steps). Thus, a 2-D nonlinear optics model from Maxwell's equations is now feasible.

Joseph, Rose↗
NASA NTRS1985

Unsteady transonic flows - Introduction, current trends, applications

The computational treatment of unsteady transonic flows is discussed, reviewing the historical development and current techniques. The fundamental physical principles are outlined; the governing equations are introduced; three-dimensional linearized and two-dimensional linear-perturbation theories in frequency domain are described in detail; and consideration is given to frequency-domain FEMs and time-domain finite-difference and integral-equation methods. Extensive graphs and diagrams are included.

Yates, E. C., Jr.↗
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 NTRS1992

Finite-difference algorithms for the time-domain Maxwell's equations - A numerical approach to RCS analysis

The applications of two CFD-based finite-difference methods to computational electromagnetics are investigated. In the first method, the time-domain Maxwell's equations are solved using the explicit Lax-Wendroff scheme and in the second method, the second-order wave equations satisfying the Maxwell's equations are solved using the implicit Crank-Nicolson scheme. The governing equations are transformed to a generalized curvilinear coordinate system and solved on a body-conforming mesh using the scattered-field formulation. The induced surface current and the bistatic radar cross section are computed and the results are validated for several two-dimensional test cases involving perfectly-conducting scatterers submerged in transverse-magnetic plane waves.

Vinh, Hoang↗
NASA NTRS1993

Shape optimization for aerodynamic efficiency and low observability

Field methods based on the finite-difference approximations of the time-domain Maxwell's equations and the potential-flow equation have been developed to solve the multidisciplinary problem of airfoil shaping for aerodynamic efficiency and low radar cross section (RCS). A parametric study and an optimization study employing the two analysis methods are presented to illustrate their combined capabilities. The parametric study shows that for frontal radar illumination, the RCS of an airfoil is independent of the chordwise location of maximum thickness but depends strongly on the maximum thickness, leading-edge radius, and leadingedge shape. In addition, this study shows that the RCS of an airfoil can be reduced without significant effects on its transonic aerodynamic efficiency by reducing the leading-edge radius and/or modifying the shape of the leading edge. The optimization study involves the minimization of wave drag for a non-lifting, symmetrical airfoil with constraints on the airfoil maximum thickness and monostatic RCS. This optimization study shows that the two analysis methods can be used effectively to design aerodynamically efficient airfoils with certain desired RCS characteristics.

Vinh, Hoang↗
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 NTRS1993

Issac, Jason Cherian ses in transonic flow

Flutter analysis of a two degree of freedom airfoil in compressible flow is performed using a state-space representation of the unsteady aerodynamic behavior. Indicial response functions are used to represent the normal force and moment response of the airfoil. The structural equations of motion of the airfoil with bending and torsional degrees of freedom are coupled to the unsteady air loads and the aeroelastic system so modelled is solved as an eigenvalue problem to determine the stability. The aeroelastic equations are also directly integrated with respect to time and the time-domain results compared with the results from the eigenanalysis. A good agreement is obtained. The derivatives of the flutter speed obtained from the eigenanalysis are calculated with respect to the mass and stiffness parameters by both analytical and finite-difference methods for various transonic Mach numbers. The experience gained from the two degree of freedom model is applied to study the sensitivity of the flutter response of a wing with respect to various shape parameters. The parameters being considered are as follows: (1) aspect ratio; (2) surface area of the wing; (3) taper ratio; and (4) sweep. The wing deflections are represented by Chebyshev polynomials. The compressible aerodynamic state-space model used for the airfoil section is extended to represent the unsteady aerodynamic forces on a generally laminated tapered skewed wing. The aeroelastic equations are solved as an eigenvalue problem to determine the flutter speed of the wing. The derivatives of the flutter speed with respect to the shape parameters are calculated by both analytical and finite difference methods.

Issac, Jason Cherion↗
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 NTRS1988

Theory and application of radiation boundary operators

A succinct unified review is provided of the theory of radiation boundary operators. With the recent introduction of the on-surface radiation condition (OSRC) method and the continued growth of finite-difference and finite-element techniques for modeling electromagnetic wave scattering problems, the understanding and use of radiation boundary operators has become increasingly important. Results are presented to illustrate the application of radiation boundary operators in both these areas. Recent OSRC results include analysis of the scattering behavior of both electrically small and large cylinders, a reactively loaded acoustic sphere, and a simple reentrant duct. Radiation boundary operator results include the demonstration of the effectiveness of higher-order operators in truncating finite-difference time-domain grids.

Moore, Thomas G.↗
NASA NTRS1991

Advanced electromagnetic methods for aerospace vehicles

During this period the research program addressed the following three topics: (1) composite materials; (2) precipitation static (P-Static); and (3) antenna technology. On the topic of Composite Materials our main efforts were directed toward making measurements on several new samples of composite materials made available to ASU by Stanford Research Institute (SRI) through the efforts of Mr. Frank Casler of AVRADA. These samples can be classified into three distinct materials with each material having its own distinct electrical properties. In addition, attempts were made to make predictions of the effects on antenna patterns by composite materials. This will take a greater emphasis in the next reporting period. In Precipitation Static (P-Static), the main effort was devoted toward developing a Voltage Finite-Difference Time-Domain computer code to account for the voltage variation on a conducting body as the primary source of corona discharge, instead of the electric field. Due to complexities stemming from the interactions between the potentials, the fields, and current sources, the decision was made to begin with a simple two-dimensional problem without the corona discharge and check our programs in a series of simple models, culminating in the full corona discharge problem. This report deals with the first stage of such development. During this reporting period, the main effort in Antenna Technology was toward the design, fabrication, and testing of a cavity-backed slot antenna using ferrite material. Using the ferrite material available to us during this period, the resonances of this antenna were around 5 and 8 GHz. Attempts will be made to model such an antenna and to lower its resonance down into the VHF and UHF bands.

Balanis, Constantine A.↗
NASA NTRS1984

Application of the Green's function method for 2- and 3-dimensional steady transonic flows

A Time-Domain Green's function method for the nonlinear time-dependent three-dimensional aerodynamic potential equation is presented. The Green's theorem is being used to transform the partial differential equation into an integro-differential-delay equation. Finite-element and finite-difference methods are employed for the spatial and time discretizations to approximate the integral equation by a system of differential-delay equations. Solution may be obtained by solving for this nonlinear simultaneous system of equations in time. This paper discusses the application of the method to the Transonic Small Disturbance Equation and numerical results for lifting and nonlifting airfoils and wings in steady flows are presented.

Tseng, K.↗
NASA NTRS1985

A first-order time-domain Green's function approach to supersonic unsteady flow

A time-domain Green's Function Method for unsteady supersonic potential flow around complex aircraft configurations is presented. The focus is on the supersonic range wherein the linear potential flow assumption is valid. The Green's function method is employed in order to convert the potential-flow differential equation into an integral one. This integral equation is then discretized, in space through standard finite-element technique, and in time through finite-difference, to yield a linear algebraic system of equations relating the unknown potential to its prescribed co-normalwash (boundary condition) on the surface of the aircraft. The arbitrary complex aircraft configuration is discretized into hyperboloidal (twisted quadrilateral) panels. The potential and co-normalwash are assumed to vary linearly within each panel. Consistent with the spatial linear (first-order) finite-element approximations, the potential and co-normalwash are assumed to vary linearly in time. The long range goal of our research is to develop a comprehensive theory for unsteady supersonic potential aerodynamics which is capable of yielding accurate results even in the low supersonic (i.e., high transonic) range.

Freedman, M. I.↗