SEARCH · Engineering Papers
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.
We demonstrate a method for generation of beams of light with large angular momenta. The method utilizes whispering gallery mode resonators that transform a plane electromagnetic wave into high order Bessel beams. Interference pattern among the beams as well as shadow pictures induced by the beams are observed and studied.
When galaxies merge, the supermassive black holes in their centers may form binaries and emit low-frequency gravitational radiation in the process. In this paper, we consider the galaxy 3C 66B, which was used as the target of the first multimessenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational-wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C 66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C 66B to less than (1.65±0.02)×109Me using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data over “blind” pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.
A consequence of an exoatmospheric nuclear burst is an electromagnetic pulse (EMP) radiated from it. In a region far enough away from the burst, where nonlinear effects can be ignored, the EMP can be represented by a large-amplitude narrow-time-width plane-wave pulse. If the ionosphere intervenes the origin and destination of the EMP, frequency dispersion can cause significant changes in the original pulse upon reception. A method of computing these dispersive effects of transient wave propagation is summarized. The method described is different from the standard transform techniques and provides physical insight into the transient wave process. The method, although exact, can be used in approximating the early-time transient response of an ionospheric region by a simple integration with only explicit knowledge of the electron density, electron collision frequency, and electron gyrofrequency required. As an illustration of the method, it is applied to a simple example and contrasted with the corresponding transform solution.
#he modal element method has been employed to determine the scattered field from a plane acoustic wave impinging on a two dimensional body. In the modal element method, the scattering body is represented by finite elements, which are coupled to an eigenfunction expansion representing the acoustic pressure in the infinite computational domain surrounding the body. The present paper extends the previous work by developing the algorithm necessary to calculate the acoustics scattering cross section by the modal element method. The scattering cross section is the acoustical equivalent to the Radar Cross Section (RCS) in electromagnetic theory. Since the scattering cross section is evaluated at infinite distance from the body, an asymptotic approximation is used in conjunction with the standard modal element method. For validation, the scattering cross section of the rigid circular cylinder is computed for the frequency range 0.1 is less than or equal to ka is less than or equal to 100. Results show excellent agreement with the analytic solution.
A semiautomatic method of fitting transition curves to X-ray film optical density measurements of electromagnetic particle cascades is described. Several hundred singly and multiple interacting cosmic ray events from the JACEE 8 balloon flights were analyzed using this procedure. In addition to greatly increased speed compared to the previous manual method, the semiautomatic method offers increased accuracy through maximum likelihood fitting.
This manual is a combination of a user manual, theory manual, and programmer manual. The reader is assumed to have some previous exposure to the finite element method. This manual is written with the idea that the CSTEM (Coupled Structural Thermal Electromagnetic-Computer Code) user needs to have a basic understanding of what the code is actually doing in order to properly use the code. For that reason, the underlying theory and methods used in the code are described to a basic level of detail. The manual gives an overview of the CSTEM code: how the code came into existence, a basic description of what the code does, and the order in which it happens (a flowchart). Appendices provide a listing and very brief description of every file used by the CSTEM code, including the type of file it is, what routine regularly accesses the file, and what routine opens the file, as well as special features included in CSTEM.
Rapid inspection of aircraft structures for flaws is of vital importance to the commercial and defense aircraft industry. In particular, inspecting thin aluminum structures for flaws is the focus of a large scale R&D effort in the nondestructive evaluation (NDE) community. Traditional eddy current methods used today are effective, but require long inspection times. New electromagnetic techniques which monitor the normal component of the magnetic field above a sample due to a sheet of current as the excitation, seem to be promising. This paper is an attempt to understand and analyze the magnetic field distribution due to a current sheet above an aluminum test sample. A simple theoretical model, coupled with a two dimensional finite element model (FEM) and experimental data will be presented in the next few sections. A current sheet above a conducting sample generates eddy currents in the material, while a sensor above the current sheet or in between the two plates monitors the normal component of the magnetic field. A rivet or a surface flaw near a rivet in an aircraft aluminum skin will disturb the magnetic field, which is imaged by the sensor. Initial results showed a strong dependence of the flaw induced normal magnetic field strength on the thickness and conductivity of the current-sheet that could not be accounted for by skin depth attenuation alone. It was believed that the eddy current imaging method explained the dependence of the thickness and conductivity of the flaw induced normal magnetic field. Further investigation, suggested the complexity associated with the mutual inductance of the system needed to be studied. The next section gives an analytical model to better understand the phenomenon.
Described are the theoretical development and computer implementation of reliable and efficient methods for the analysis of coupled mechanical problems that involve the interaction of mechanical, thermal, phase-change and electromag subproblems. The focus application has been the modeling of superconductivity and associated quantum-state phase change phenomena. In support of this objective the work has addressed the following issues: (1) development of variational principles for finite elements, (2) finite element modeling of the electromagnetic problem, (3) coupling of thermel and mechanical effects, and (4) computer implementation and solution of the superconductivity transition problem. The main accomplishments have been: (1) the development of the theory of parametrized and gauged variational principles, (2) the application of those principled to the construction of electromagnetic, thermal and mechanical finite elements, and (3) the coupling of electromagnetic finite elements with thermal and superconducting effects, and (4) the first detailed finite element simulations of bulk superconductors, in particular the Meissner effect and the nature of the normal conducting boundary layer. The theoretical development is described in two volumes. Volume 1 describes mostly formulation specific problems. Volume 2 describes generalization of those formulations.
We use the superposition T-matrix method to compare the far-field scattering matrices generated by spheroidal and spherical volumes of discrete random medium having the same volume and populated by identical spherical particles. Our results fully confirm the robustness of the previously identified coherent and diffuse scattering regimes and associated optical phenomena exhibited by spherical particulate volumes and support their explanation in terms of the interference phenomenon coupled with the order-of-scattering expansion of the far-field Foldy equations. We also show that increasing non-sphericity of particulate volumes causes discernible (albeit less pronounced) optical effects in forward and backscattering directions and explain them in terms of the same interference/multiple-scattering phenomenon.
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.
The paper studies the inversion of the radiative transfer equation describing the interaction of electromagnetic radiation with atmospheric aerosols. The interaction can be considered as the propagation in the aerosol medium of two light beams: the direct beam in the line-of-sight attenuated by absorption and scattering, and the diffuse beam arising from scattering into the viewing direction, which propagates more or less in random fashion. The latter beam has single scattering and multiple scattering contributions. In the former case and for single scattering, the problem is reducible to first-kind Fredholm equations, while for multiple scattering it is necessary to invert partial integrodifferential equations. A nonlinear minimization search method, applicable to the solution of both types of problems has been developed, and is applied here to the problem of monitoring aerosol pollution, namely the complex refractive index and size distribution of aerosol particles.
Traditionally microstrip or printed reflectarrays are designed using the transmit mode technique. In this method, the size of each printed element is chosen so as to provide the required value of the reflection phase such that a collimated beam results along a given direction. The reflection phase of each printed element is approximated using an infinite array model. The infinite array model is an excellent engineering approximation for a large microstrip array since the size or orientation of elements exhibits a slow spatial variation. In this model, the reflection phase from a given printed element is approximated by that of an infinite array of elements of the same size and orientation when illuminated by a local plane wave. Thus the reflection phase is a function of the size (or orientation) of the element, the elevation and azimuth angles of incidence of a local plane wave, and polarization. Typically, one computes the reflection phase of the infinite array as a function of several parameters such as size/orientation, elevation and azimuth angles of incidence, and in some cases for vertical and horizontal polarization. The design requires the selection of the size/orientation of the printed element to realize the required phase by interpolating or curve fitting all the computed data. This is a substantially complicated problem, especially in applications requiring a computationally intensive commercial code to determine the reflection phase. In dual polarization applications requiring rectangular patches, one needs to determine the reflection phase as a function of five parameters (dimensions of the rectangular patch, elevation and azimuth angles of incidence, and polarization). This is an extremely complex problem. The new method employs the reciprocity principle and reaction concept, two well-known concepts in electromagnetics to derive the receive mode analysis and design techniques. In the "receive mode design" technique, the reflection phase is computed for a plane wave incident on the reflectarray from the direction of the beam peak. In antenna applications with a single collimated beam, this method is extremely simple since all printed elements see the same angles of incidence. Thus the number of parameters is reduced by two when compared to the transmit mode design. The reflection phase computation as a function of five parameters in the rectangular patch array discussed previously is reduced to a computational problem with three parameters in the receive mode. Furthermore, if the beam peak is in the broadside direction, the receive mode design is polarization independent and the reflection phase computation is a function of two parameters only. For a square patch array, it is a function of the size, one parameter only, thus making it extremely simple.
The T-matrix method is one of the most versatile, efficient, and accurate theoretical techniques widely used for numerically exact computer calculations of electromagnetic scattering by single and composite particles, discrete random media, and particles imbedded in complex environments. This paper presents the fifth update to the comprehensive database of peer-reviewed T-matrix publications initiated by us in 2004 and includes relevant publications that have appeared since 2012. It also lists several earlier publications not incorporated in the original database, including Peter Waterman's reports from the 1960s illustrating the history of the T-matrix approach and demonstrating that John Fikioris and Peter Waterman were the true pioneers of the multi-sphere method otherwise known as the generalized Lorenz - Mie theory.
The design, development, application, and capabilities of a variable frequency microwave radiometer are described. This radiometer demonstrated the versatility, accuracy, and stability required to provide contributions to the geophysical understanding of ocean and ice processes. A closed-loop feedback method was used, whereby noise pulses were added to the received electromagnetic radiation to achieve a null balance in a Dicke switched radiometer. Stability was achieved through the use of a constant temperature enclosure around the low loss microwave front end. The Dicke reference temperature was maintained to an absolute accuracy of 0.1 K using a closed-loop proportional temperature controller. A microprocessor based digital controller operates the radiometer and records the data on computer compatible tapes. This radiometer exhibits an absolute accuracy of better than 0.5 K when the sensitivity is 0.1 K. The sensitivity varies between 0.0125 K and 1.25 K depending upon the bandwidth and integration time selected by the digital controller. Remote sensing experiments were conducted from an aircraft platform and the first radiometeric mapping of an ocean polar front; exploratory experiments to measure the thickness of lake ice; first discrimination between first year and multiyear ice below 10 GHz; and the first known measurements of frequency sensitive characteristics of sea ice.
The applicability of a parallel architecture to the solution of large electromagnetic scattering problems is demonstrated. Two techniques, finite difference and the method of moments, are used to provide insight into the comparative speedups which can be attained for several algorithms. The flexibility of the hypercube architecture for different analysis algorithms is illustrated.
Electro-Impulse De-Icing (EIDI) is a method of ice removal by sharp blows delivered by a transient electromagnetic field. Detailed results are given for studies of the electrodynamic phenomena. Structural dynamic tests and computations are described. Also reported are ten sets of tests at NASA's Icing Research Tunnel and flight tests by NASA and Cessna Aircraft Company. Fabrication of system components are described and illustrated. Fatigue and electromagnetic interference tests are reported. Here, the necessary information for the design of an EIDI system for aircraft is provided.
This paper presents the major characteristics of a simulation package capable of performing a complete narrow and wideband analysis of the mobile satellite communication channel in urban environments for any given orbital configuration. The wavelength-to-average urban geometrical dimension ratio has required the use of the Geometrical Theory of Diffraction (GTD). For the RF frequency range, the model has been designed to be (1 up to 60 GHz) extended to include effects of non-perfect conductivity and surface roughness. Taking advantage of the inherent capabilities of such a high frequency method, we are able to provide a complete description of the electromagnetic field at the mobile terminal. Using the information made available at the ray-tracer and GTD solver outputs, the Land Mobile Satellite (LMS) urban model can also give a detailed description of the communication channel in terms of power delay profiles, Doppler spectra, channel scattering functions, and so forth. Statistical data, e.g. cumulative distribution functions, level crossing rates or distributions of fades are also provided. The user can access the simulation tool through a Design-CAD user-friendly interface by means of which she can effectively design her own urban layout and run consequently all the envisaged routines. The software is optimized in its execution time so that numerous runs can be achieved in a considerably short time.
The recent discovery of High Temperature Superconductors (HTS) with superconducting transition temperature, T(sub c), above the boiling point of liquid nitrogen has opened the door for using these materials in new and practical applications. These materials have zero resistance to electric current, have the capability of carrying large currents and as such have the potential to be used in high magnetic field applications. One of the space applications that can use superconductors is electromagnetic launch of payloads to low-earth-orbit. An electromagnetic gun-type launcher can be used in small payload systems that are launched at very high velocity, while sled-type magnetically levitated launcher can be used to launch larger payloads at smaller velocities. Both types of launchers are being studied by NASA and the aerospace industry. The use of superconductors will be essential in any of these types of launchers in order to produce the large magnetic fields required to obtain large thrust forces. Low Temperature Superconductor (LTS) technology is mature enough and can be easily integrated in such systems. As for the HTS, many leading companies are currently producing HTS coils and magnets that potentially can be mass-produced for these launchers. It seems that designing and building a small-scale electromagnetic launcher is the next logical step toward seriously considering this method for launching payloads into low-earth-orbit. A second potential application is the use of HTS to build sensitive portable devices for the use in Non Destructive Evaluation (NDE). Superconducting Quantum Interference Devices (SQUID's) are the most sensitive instruments for measuring changes in magnetic flux. By using HTS in SQUID's, one will be able to design a portable unit that uses liquid nitrogen or a cryocooler pump to explore the use of gradiometers or magnetometers to detect deep cracks or corrosion in structures. A third use is the replacement of Infra-Red (IR) sensor leads on Earth Orbit Systems (EOS) with HTS leads. IR detectors on these EOS missions are cooled to a 4.2K to improve their signal to noise ratio. They are connected to data acquisitions systems using manganin wires (low thermal conductors) to reduce the heat load on the cryogen. Replacing these wires with HTS leads will increase the lifetime of these missions by about 50 percent. This is a promising application that is ready for actual implementation on such systems. The analysis also show that an the number of IR detectors increase in larger EOS systems, substantial increase in the lifetime of each mission will be realized by using HTS leads instead of the manganin ones.