Search NASA⌕ Search

SEARCH · Search NASA

Results for “computational electromagnetics”

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 127 records · Page 7

Canonical Statistical Model for Maximum Expected Immission of Wire Conductor in an Aperture Enclosure

Prediction of the maximum expected electromagnetic pick-up of conductors inside a realistic shielding enclosure is an important canonical problem for system-level EMC design of space craft, launch vehicles, aircraft and automobiles. This paper introduces a simple statistical power balance model for prediction of the maximum expected current in a wire conductor inside an aperture enclosure. It calculates both the statistical mean and variance of the immission from the physical design parameters of the problem. Familiar probability density functions can then be used to predict the maximum expected immission for deign purposes. The statistical power balance model requires minimal EMC design information and solves orders of magnitude faster than existing numerical models, making it ultimately viable for scaled-up, full system-level modeling. Both experimental test results and full wave simulation results are used to validate the foundational model.

Computational electromagnetics↗

An MPMD approach coupling electromagnetic continuum mechanics approximations in ALEGRA

In this work, two complementary approximations for describing aspects of continuum electromagnetics in moving media are discussed: electroquasistatic and magnetoquasistatic. Each has been implemented in the finite element shock code ALEGRA for modeling dynamic electromechanical phenomena on typical engineering time scales, with fully integrated circuit coupling. The approximations can be obtained by consistent asymptotic balancing of Maxwell’s equations relative to timescales associated with magnetic diffusion, charge relaxation, and electromagnetic wave propagation. In ALEGRA, the electroquasistatic approximation is used for ferroelectric (FE) modeling, while the magnetoquasistatic approximation is used for magnetohydrodynamic (MHD) modeling. In this paper we introduce for the first time a detailed derivation of a useful quasi-steady “low-R m ” variant of the MHD approximation applicable for cases, such as with detonators, where the thermodynamic pressure arising from Joule heating dominates over magnetic forces. An additional purpose of this paper is to present a coupling mode using Multiple Program-Multiple Data (MPMD) message passing communication that allows the user to run 3D FE problems together with 2D and/or 3D MHD problems with the respective simulation domains coupled through a common circuit equation. The MPMD coupling capability is used here to model the dynamic coupling of a notional ferroelectric generator with an RP-87 exploding bridgewire detonator. The simulated bridgewire heats up and bursts under current generated by simulated depoling of the ferroelectric generator, as a demonstration of the MPMD capability.

42 ENGINEERING↗

Using the Dimensionality Reduction (DR) Approach to Treat Singular and Near-Singular Source and Test Integrals on Triangles for Moment Methods

Recently, the authors presented preliminary results of an initial study of a dimensional reduction scheme for treating (near-)singular integrals (D. R. Wilton et al., “Dimensionality Reduction Approach for Treating Singular and Near-Singular Multidimensional Integrals of Electromagnetics,” 2022 URSI USNC National Radio Science Meeting, Boulder, CO, Jan. 2022). The approach reported there refined and extended several ideas appearing in a recent paper (R. R. Chang, Z. Wang, and Q. Xie, "Fast Convergent Quadrature Method for Evaluating the RWG- and SWG-Related Convolutional Integrals," IEEE Trans. Antennas Propagat., 69, Dec. 2021). That paper reported a curated collection of important, previously published results all leading to very efficient and accurate line integral methods for handling the most common (near-)singular integrals associated with triangles and tetrahedrons with RWG and SWG bases, respectively. Our contributions provided a simpler exposition as well as a more unified, systematic, and robust overall framework for deriving and applying the approach. This presentation further elucidates our extensions to the approach and broadens our initial study to examine the convergence behaviors of the various potential forms over a wider range of triangular source element shape and observation point parameters for various smoothing transformations; we also examine convergence, not just at the subtriangle level, but also overcomplete triangles. In addition, we investigate the application of the recently reported “vertex function” concept to develop faster converging test integrals (Rivero et al., “Acceleration of the Surface Test Integral Using Vertex Functions,” 2021 IEEE Int’l Symp. Ant. Propagat. and USNC-URSI Rad. Sci. Meeting,” Singapore, 4-10 Dec. 2021). Combining these separate schemes has the potential for yielding a near-optimal approach for accurately evaluating singular and near-singular integrals for moment methods.

Singular Integrals↗

A Generic and Multifunctional Electromagnetic Transient Model for Grid-Following Inverters

This article presents a generic and multi-functional electromagnetic transient (EMT) dynamic model of grid following (GFL) inverter-based resource (IBR) using the PSCAD software platform. The features of the developed model includes the flexibility in selecting various types and combinations of DC sources covering PV modules, battery modules, ideal DC source module, as well as flexibility in selecting either switched or averaged model of inverter. This model also covers exhaustive lists of controller logic covering open-loop/closed-loop PQ dispatch control, DC voltage and AC terminal voltage control along with the conventional current control designed in dq-domain, ate- domain and positive-negative sequence domain. Moreover, this model is equipped with the flexibility in selecting various types of current limiting schemes that includes saturation-based as well as latching-based current limiter, anti-windup protection. Moreover, the EMT model is agnostic to the MVA rating and is suitable for interfacing transmission systems by being complaint with the IEEE Std. 2800. The generality in the power circuits and the multi-functional options in operation and control of the developed EMT model makes it suitable for both academia and industry to study various power system aspects not limited to but such as fault behavior of GFL IBR and impacts on protection system, transient stability of a system interfaced with large number of GFL IBRs etc.

14 SOLAR ENERGY↗

Numerical check on the accuracy of spherical-wave expansions.

A simple and sensitive test for evaluating the numerical accuracy of a spherical-wave expansion is presented. The basic idea is to compare two different expansions of the same electromagnetic field, computed with respect to different coordinate origins. The use of the technique is illustrated by an example based on the experimental pattern of a horn antenna.

Ludwig, A. C.↗

Rigid format alter packets for the analysis of electromagnetic field problems

The computer program NASTRAN is used to solve electromagnetic field problems. The diffusion equation and the boundary conditions valid for problems of these kinds together with a replacing potential energy function are given. The extent to which an analogy with finite element displacement and temperature approaches holds is indicated. The outputting of complex quantities is made possible after adjustment of standard rigid format 1 input data blocks to module SDR2. The applications made involve the study of the proximity effect in a system of three parallel conductors and the analysis of the magnetic field in the vicinity of the points of contact in circuit breakers.

Spreeuw, E.↗

A combined finite element-boundary element formulation for solution of two-dimensional problems via CGFFT

A method for the computation of electromagnetic scattering from arbitrary two-dimensional bodies is presented. The method combines the finite element and boundary element methods leading to a system for solution via the conjugate gradient Fast Fourier Transform (FFT) algorithm. Two forms of boundaries aimed at reducing the storage requirement of the boundary integral are investigated. It is shown that the boundary integral becomes convolutional when a circular enclosure is chosen, resulting in reduced storage requirement when the system is solved via the conjugate gradient FFT method. The same holds for the ogival enclosure, except that some of the boundary integrals are not convolutional and must be carefully treated to maintain O(N) memory requirement. Results for several circular and ogival structures are presented and shown to be in excellent agreement with those obtained by traditional methods.

Collins, Jeffery D.↗

Determination of the ocean circulation using Geosat altimetry

A spherical harmonic model of the sea surface topography complete to degree and order 10 and a model of the earth's geopotential field complete to degree and order 50 have been obtained in a simultaneous solution using Geosat altimeter data and tracking data from 14 different satellites. The sea surface topography model compares well with oceanographic models computed using hydrographic data and ship drift data. Currently, errors in the estimated gravity field model limit the determination of the spherical harmonic coefficients of the general ocean circulation to degrees 10 and lower, corresponding to a minimum wavelength of 4000 km. Error analysis indicates that the correlation between the geoid and the sea surface topography model is less than 0.2, indicating good separation of the geoid and the sea surface topography at wavelengths of 4000 km or longer. Estimates of the scale factor for the significant wave height (H1/3), which is used to compute the electromagnetic bias correction and the bias for the Geosat altimeter, are obtained. The estimate of the H1/3 correction is 3.6 + or - 1.5 percent, and the height bias estimate is zero.

Nerem, R. S.↗

Alpha particle heating in hot diamagnetic cavities

Observational data from ISEE 1 are analyzed and a one-dimensional electromagnetic hybrid computer simulation is conducted for the heating of solar wind alpha particles in hot diamagnetic cavities (HDCs). In the simulation, which envisions alpha heating by ion-ion instabilities, low beam densities excite the proton/proton right-hand resonant instability that then pitch-angle scatters the beam without significantly heating the alphas. At greater beam densities, the proton/proton nonresonant instability undergoes saturation through a trapping of all three ion components. These results support the Thomsen et al. (1988) hypothesis that the nonresonant instability is the primary source of ion heating in hot diamagnetic cavities.

Galvez, Miguel↗

Free-floating experiments in orbiter and space station

Manned spacecraft will never be free of acceleration perturbations due to necessary motions of crewmembers. However, in most cases, the perturbations can be minimized or even isolated from sensitive experiments when required. The principal disturbances and methods of isolation are described. As the crew moves about the spacecraft, no external forces are applied and the average acceleration remains near zero. External force will be exerted when jets are fired for spacecraft attitude control or translation, and these must be transferred to the microgravity experiments. Transitory acceleration and vibration are produced by crew and equipment. The most significant on Orbiter have been wall push-offs used for body translation within the spacecraft, the exercise conducted on a treadmill within the middeck area, closing doors of stowage compartments and the vibration of some machinery, such as a sample centrifuge. For short periods of perhaps a few minutes, crew and equipment motion can be largely inhibited, but not for much longer intervals. It may be necessary to mount especially sensitive experiments on an isolation table, which can greatly reduce the acceleration transferred to the experiment. Either weak mechanical springs or, even better, a computer-controlled electromagnetic suspension can be used to effectively insulate the table from vibration at frequencies above about 0.01 Hz. Another option is to free-float an experiment table within the spacecraft. All of the above-mentioned perturbations can be eliminated for some minutes if the entire experiment package is allowed to free-float within the spacecraft until a wall contact is made. Package acceleration levels should be maintainable below 10(exp -8) g sub o in this way.

Garriott, Owen K.↗

Fly-by-light technology development plan

The driving factors and developments which make a fly-by-light (FBL) viable are discussed. Documentation, analyses, and recommendations are provided on the major issues pertinent to facilitating the U.S. implementation of commercial FBL aircraft before the turn of the century. Areas of particular concern include ultra-reliable computing (hardware/software); electromagnetic environment (EME); verification and validation; optical techniques; life-cycle maintenance; and basis and procedures for certification.

Todd, J. R.↗

Perfectly Matched Layer for Linearized Euler Equations in Open and Ducted Domains

Recently, perfectly matched layer (PML) as an absorbing boundary condition has widespread applications. The idea was first introduced by Berenger for electromagnetic waves computations. In this paper, it is shown that the PML equations for the linearized Euler equations support unstable solutions when the mean flow has a component normal to the layer. To suppress such unstable solutions so as to render the PML concept useful for this class of problems, it is proposed that artificial selective damping terms be added to the discretized PML equations. It is demonstrated that with a proper choice of artificial mesh Reynolds number, the PML equations can be made stable. Numerical examples are provided to illustrate that the stabilized PML performs well as an absorbing boundary condition. In a ducted environment, the wave mode are dispersive. It will be shown that the group velocity and phase velocity of these modes can have opposite signs. This results in a confined environment, PML may not be suitable as an absorbing boundary condition.

Tam, Christopher K. W.↗

Interfacing WIPL-D with Mechanical CAD Software

of almost any popular CAD format, e.g. IGES, Parasolid, DXF, ACIS etc. The solid models are processed (simplified) and meshed in GiD(R), and then converted into WIPL-D Pro input file by simple Fortran or Matlab code. This algorithm allows the user to control the mesh of imported geometry, and to assign electric p~perties to metalic and dielectric surfaces. Implementation of the algorithm is demonstrated by examples obtained from the NASA Discovery mission, Phoenix Lander 2008. Results for radiation pattern of Phoenix Lander UHF relay antenna with effect of Martian surface, both simulated in WIPL-D Pro and measured, are shown for comparison.

method of moments (MoM)↗

Comprehensive T-Matrix Reference Database: A 2007-2009 Update

The T-matrix method is among the most versatile, efficient, and widely used theoretical techniques for the numerically exact computation of electromagnetic scattering by homogeneous and composite particles, clusters of particles, discrete random media, and particles in the vicinity of an interface separating two half-spaces with different refractive indices. This paper presents an update to the comprehensive database of T-matrix publications compiled by us previously and includes the publications that appeared since 2007. It also lists several earlier publications not included in the original database.

Mishchenko, Michael I.↗

Comprehensive T-matrix Reference Database: A 2009-2011 Update

The T-matrix method is one of the most versatile and efficient theoretical techniques widely used for the computation of electromagnetic scattering by single and composite particles, discrete random media, and particles in the vicinity of an interface separating two half-spaces with different refractive indices. This paper presents an update to the comprehensive database of peer-reviewed T-matrix publications compiled by us previously and includes the publications that appeared since 2009. It also lists several earlier publications not included in the original database.

electromagnetic scattering↗

Comprehensive Thematic T-Matrix Reference Database: A 2014-2015 Update

The T-matrix method is one of the most versatile and efficient direct computer solvers of the macroscopic Maxwell equations and is widely used for the computation of electromagnetic scattering by single and composite particles, discrete random media, and particles in the vicinity of an interface separating two half-spaces with different refractive indices. This paper is the seventh update to the comprehensive thematic database of peer-reviewed T-matrix publications initiated by us in 2004 and includes relevant publications that have appeared since 2013. It also lists a number of earlier publications overlooked previously.

electromagnetic scattering↗

Comprehensive Thematic T-Matrix Reference Database: A 2015-2017 Update

The T-matrix method pioneered by Peter C. Waterman is one of the most versatile and efficient numerically exact computer solvers of the time-harmonic macroscopic Maxwell equations. It is widely used for the computation of electromagnetic scattering by single and composite particles, discrete random media, periodic structures (including metamaterials), and particles in the vicinity of plane or rough interfaces separating media with different refractive indices. This paper is the eighth update to the comprehensive thematic database of peer-reviewed T-matrix publications initiated in 2004 and lists relevant publications that have appeared since 2015. It also references a small number of earlier publications overlooked previously.

T-matrix method↗