Radio waves and circuits
Diffraction and scattering of electromagnetic waves, radio wave propagation in ionized media, radio antennas, and signal processing in communications systems
SEARCH · Search NASA
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.
Diffraction and scattering of electromagnetic waves, radio wave propagation in ionized media, radio antennas, and signal processing in communications systems
Plasma effects on electromagnetic wave propagation during Apollo spacecraft reentry from lunar missions
Electromagnetic wave propagation in partially ionized paramagnetic gas in static magnetic field
Equipment and instrumentation for acoustic simulation of electromagnetic wave propagation and radar systems design studies
The small-angle scattering approximation of the scalar radiative transfer equation (RTE) is examined for the case where the extinction and scattering coefficients have a component that is a deterministic function of position along the propagation path and a component that is a random function of position transverse to the propagation direction. It is found that the resulting stochastic RTE can be reduced to a system of two stochastic integrodifferential equations for the average and fluctuating components of the radiant intensity. Two transfer equations are obtained describing the average radiant intensity and the spatial correlation function of the intensity fluctuations. The average intensity equation is then solved and applied to a simple propagation scenario. An approximate solution is also derived for the equation giving the correlation function. The developed equations can be applied to problems involving short wavelength electromagnetic wave propagation through media possessing the variable characteristics of turbulence and turbidity, such as plasmas, the atmosphere, and the ocean.
A proposal has been made to exploit the special wavelength-dispersive characteristics of devices of the type described in One-Dimensional Photonic Crystal Superprisms (NPO-30232) NASA Tech Briefs, Vol. 29, No. 4 (April 2005), page 10a. A photonic crystal is an optical component that has a periodic structure comprising two dielectric materials with high dielectric contrast (e.g., a semiconductor and air), with geometrical feature sizes comparable to or smaller than light wavelengths of interest. Experimental superprisms have been realized as photonic crystals having three-dimensional (3D) structures comprising regions of amorphous Si alternating with regions of SiO2, fabricated in a complex process that included sputtering. A photonic crystal of the type to be exploited according to the present proposal is said to be one-dimensional (1D) because its contrasting dielectric materials would be stacked in parallel planar layers; in other words, there would be spatial periodicity in one dimension only. The processes of designing and fabricating 1D photonic crystal superprisms would be simpler and, hence, would cost less than do those for 3D photonic crystal superprisms. As in 3D structures, 1D photonic crystals may be used in applications such as wavelength-division multiplexing. In the extended-range configuration, it is also suitable for spectrometry applications. As an engineered structure or artificially engineered material, a photonic crystal can exhibit optical properties not commonly found in natural substances. Prior research had revealed several classes of photonic crystal structures for which the propagation of electromagnetic radiation is forbidden in certain frequency ranges, denoted photonic bandgaps. It had also been found that in narrow frequency bands just outside the photonic bandgaps, the angular wavelength dispersion of electromagnetic waves propagating in photonic crystal superprisms is much stronger than is the angular wavelength dispersion obtained by use of conventional prisms and diffraction gratings and is highly nonlinear.
A two-dimensional particle-in-cell code for the simulation of low-frequency electromagnetic processes in laboratory plasmas has been developed. The code uses the Darwin method omitting the electromagnetic wave propagation. The Darwin method separates the electric field into solenoidal and irrotational parts. The solenoidal electric field is calculated with a new algorithm based on the equation for the electric field vorticity. The system of linear equations in the new algorithm is readily solved using a standard iterative method. The irrotational electric field is the electrostatic field calculated with the direct implicit algorithm. The code is verified by reproducing the two-stream instability, electron electromagnetic waves, and shear Alfvén waves. The code is applied to simulate an inductively coupled plasma with the driving current flowing around the plasma region. In this simulation, a ring of dense plasma forms at the initial stage but then the density becomes maximal in the center and decays monotonically toward the walls. The skin effect is in the transitional mode between local and non-local, and the electron velocity distribution function is non-Maxwellian.
There is an emerging technology of photonic (or optoelectronic) integrated circuits (PICs or OEICs). In PICs, optical and electronic components are grown together on the same chip. rib build such devices and subsystems, one needs to model the entire chip. Accurate computer modeling of electromagnetic wave propagation in semiconductors is necessary for the successful development of PICs. More specifically, these computer codes would enable the modeling of such devices, including their subsystems, such as semiconductor lasers and semiconductor amplifiers in which there is femtosecond pulse propagation. Here, the computer simulations are made by solving the full vector, nonlinear, Maxwell's equations, coupled with the semiconductor Bloch equations, without any approximations. The carrier is retained in the description of the optical pulse, (i.e. the envelope approximation is not made in the Maxwell's equations), and the rotating wave approximation is not made in the Bloch equations. These coupled equations are solved to simulate the propagation of femtosecond optical pulses in semiconductor materials. The simulations describe the dynamics of the optical pulses, as well as the interband and intraband.
For a macroscopically plane-parallel discrete random medium, the boundary conditions for the specific coherency dyadic at a rough interface are derived. The derivation is based on a modification of the Twersky approximation for a scattering system consisting of a group of particles and the rough surface, and reduces to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave propagating in a discrete random medium with non-scattering boundaries. In a matrix-form setting, the boundary conditions for the specific coherency dyadic imply the boundary conditions for specific intensity column vectors which in turn, yield the expressions for the reflection and transmission matrices. The derived expressions are shown to be identical to those obtained by applying a phenomenological approach based on a facet model to the solution of the scattering problem for a rough surface illuminated by a plane electromagnetic wave.
Ribbon waveguides made of alumina or of semiconductors (Si, InP, or GaAs) have been proposed as low-loss transmission lines for coupling electronic components and circuits that operate at frequencies from 30 to 1,000 GHz. In addition to low losses (and a concomitant ability to withstand power levels higher than would otherwise be possible), the proposed ribbon waveguides would offer the advantage of compatibility with the materials and structures now commonly incorporated into integrated circuits. Heretofore, low-loss transmission lines for this frequency range have been unknown, making it necessary to resort to designs that, variously, place circuits and components to be coupled in proximity of each other and/or provide for coupling via free space through bulky and often lossy optical elements. Even chip-to-chip interconnections have been problematic in this frequency range. Metal wave-guiding structures (e.g., microstriplines and traditional waveguides) are not suitable for this frequency range because the skin depths of electromagnetic waves in this frequency range are so small as to give rise to high losses. Conventional rod-type dielectric waveguide structures are also not suitable for this frequency range because dielectric materials, including ones that exhibit ultralow losses at lower frequencies, exhibit significant losses in this frequency range. Unlike microstripline structures or metallic waveguides, the proposed ribbon waveguides would be free of metal and would therefore not be subject to skin-depth losses. Moreover, although they would be made of materials that are moderately lossy in the frequency range of interest, the proposed ribbon waveguides would cause the propagating electromagnetic waves to configure themselves in a manner that minimizes losses.
A study to determine the feasibility of employing beamed electromagnetic energy for vehicle propulsion within and outside the Earth's atmosphere was co-funded by NASA and the Defense Advanced Research Projects Agency that began in June 2010 and culminated in a Summary Presentation in April 2011. A detailed report entitled "Beamed-Energy Propulsion (BEP) Study" appeared in February 2012 as NASA/TM-2012-217014. Of the very many nuances of this subject that were addressed in this report, the effects of transferring the required high energy-density electromagnetic fields through the atmosphere were discussed. However, due to the limitations of the length of the report, only a summary of the results of the detailed analyses were able to be included. It is the intent of the present work to make available the complete analytical modeling work that was done for the BEP project with regard to electromagnetic wave propagation issues. In particular, the present technical memorandum contains two documents that were prepared in 2011. The first one, entitled "Effects of Beaming Energy Through the Atmosphere" contains an overview of the analysis of the nonlinear problem inherent with the transfer of large amounts of energy through the atmosphere that gives rise to thermally-induced changes in the refractive index; application is then made to specific beamed propulsion scenarios. A brief portion of this report appeared as Appendix G of the 2012 Technical Memorandum. The second report, entitled "An Analytical Assessment of the Thermal Blooming Effects on the Propagation of Optical and Millimeter- Wave Focused Beam Waves For Power Beaming Applications" was written in October 2010 (not previously published), provides a more detailed treatment of the propagation problem and its effect on the overall characteristics of the beam such as its deflection as well as its radius. Comparisons are then made for power beaming using the disparate electromagnetic wavelengths of 1.06 microns and 2.0 millimeters..
We study the numerical propagation of waves through future null infinity in a conformally compactified spacetime. We introduce an artificial cosmological constant, which allows us some control over the causal structure near null infinity. We exploit this freedom to ensure that all light cones are tilted outward in a region near null infinity, which allows us to impose excision-style boundary conditions in our finite difference code. In this preliminary study we consider electromagnetic waves propagating in a static, conformally compactified spacetime.
This study reveals how dependent scattering and microstructure significantly affect electromagnetic wave propagation through aerogel monoliths, contributing to their transparency. Light scattering by particle ensembles is considered “dependent” when the scattering properties rely not only on particle size and optical constants but also on their spatial distribution, typically occurring when the average interparticle distance is small in comparison with the wavelength of incident radiation. Addressing dependent scattering requires solving Maxwell’s equations for complex heterogeneous structures, which is computationally demanding and usually limited to sample thicknesses on the same scale as the wavelength. This study combines computer-generated ambigel microstructures of fractal aggregates of polydisperse nanoparticles and the radiative transfer with reciprocal transaction method to predict the transmittance of thick ambigel slabs. Transmittance measurements of ambiently dried aerogel monoliths (ambigels) with porosities from about 50% to 90% closely matched the predicted values for their digital twins. However, ignoring dependent scattering or particle aggregation led to inaccurate predictions. This study validated the computational framework, and its findings offer insights for designing photonic metamaterials and analyzing their interactions with electromagnetic waves.
In this study we introduce EMPIRE-PIC, a finite element method particle-in-cell (FEM-PIC) application developed at Sandia National Laboratories. The code has been developed in C++ using the Trilinos library and the Kokkos Performance Portability Framework to enable running on multiple modern compute architectures while only requiring maintenance of a single codebase. EMPIRE-PIC is capable of solving both electrostatic and electromagnetic problems in two- and three-dimensions to second-order accuracy in space and time. In this paper we validate the code against three benchmark problems — a simple electron orbit, an electrostatic Langmuir wave, and a transverse electromagnetic wave propagating through a plasma. We demonstrate the performance of EMPIRE-PIC on four different architectures: Intel Haswell CPUs, Intel's Xeon Phi Knights Landing, ARM Thunder-X2 CPUs, and NVIDIA Tesla V100 GPUs attached to IBM POWER9 processors. This analysis demonstrates scalability of the code up to more than two thousand GPUs, and greater than one hundred thousand CPUs.
The numerical simulation of many aerodynamic non-periodic flows of practical interest involves discretized computational domains that often must be artificially truncated. Appropriate boundary conditions are required at these truncated domain boundaries, and ideally, these boundary conditions should be perfectly "absorbing" or "nonreflecting" so that they do not contaminate the flow field in the interior of the domain. The proper specification of these boundaries is critical to the stability, accuracy, convergence, and quality of the numerical solution, and has been the topic of considerable research. The need for accurate boundary specification has been underscored in recent years with efforts to apply higher-fidelity methods (DNS, LES) in conjunction with high-order low-dissipation numerical schemes to realistic flow configurations. One of the most popular choices for specifying these boundaries is the characteristics-based boundary condition where the linearized flow field at the boundaries are decomposed into characteristic waves using either one-dimensional Riemann or other multi-dimensional Riemann approximations. The values of incoming characteristics are then suitably modified. The incoming characteristics are specified at the in flow boundaries, and at the out flow boundaries the variation of the incoming characteristic is zeroed out to ensure no reflection. This, however, makes the problem ill-posed requiring the use of an ad-hoc parameter to allow small reflections that make the solution stable. Generally speaking, such boundary conditions work reasonably well when the characteristic flow direction is normal to the boundary, but reflects spurious energy otherwise. An alternative to the characteristic-based boundary condition is to add additional "buffer" regions to the main computational domain near the artificial boundaries, and solve a different set of equations in the buffer region in order to minimize acoustic reflections. One approach that has been used involves modeling the pressure fluctuations as acoustic waves propagating in the far-field relative to a single noise-source inside the buffer region. This approach treats vorticity-induced pressure fluctuations the same as acoustic waves. Another popular approach, often referred to as the "sponge layer," attempts to dampen the flow perturbations by introducing artificial dissipation in the buffer region. Although the artificial dissipation removes all perturbations inside the sponge layer, incoming waves are still reflected from the interface boundary between the computational domain and the sponge layer. The effect of these refkections can be somewhat mitigated by appropriately selecting the artificial dissipation strength and the extent of the sponge layer. One of the most promising variants on the buffer region approach is the Perfectly Matched Layer (PML) technique. The PML technique mitigates spurious reflections from boundaries and interfaces by dampening the perturbation modes inside the buffer region such that their eigenfunctions remain unchanged. The technique was first developed by Berenger for application to problems involving electromagnetic wave propagation. It was later extended to the linearized Euler, Euler and Navier-Stokes equations by Hu and his coauthors. The PML technique ensures the no-reflection property for all waves, irrespective of incidence angle, wavelength, and propagation direction. Although the technique requires the solution of a set of auxiliary equations, the computational overhead is easily justified since it allows smaller domain sizes and can provide better accuracy, stability, and convergence of the numerical solution. In this paper, the PML technique is developed in the context of a high-order spectral-element Discontinuous Galerkin (DG) method. The technique is compared to other approaches to treating the in flow and out flow boundary, such as those based on using characteristic boundary conditions and sponge layers. The superiority of the current PML technique over other approaches is demonstrated for a range of test cases, viz., acoustic pulse propagation, convective vortex, shear layer flow, and low-pressure turbine cascade flow. The paper is structured as follows. We first derive the PML equations from the non{linear Euler equations. A short description of the higher-order DG method used is then described. Preliminary results for the four test cases considered are then presented and discussed. Details regarding current work that will be included in the final paper are also provided.
Electromagnetic waves propagating through vacuum can polarize virtual electron–positron pairs; this polarization, in turn, nonlinearly modifies their propagation. A semi-classical nonlinear wave equation describing the propagation is derived from the Euler–Heisenberg Lagrangian density, which captures vacuum polarization effects up to the one-loop level. In this article, we present a reduced-actionintegral approach that enables rapid modeling of nonlinear phenomena arising from the Euler– Heisenberg Lagrangian. Application of the variational principle to the reduced action provides equations of motion for familiar light-pulse parameters, such as spot size, phase, polarization, and phase-front curvature, without requiring full-field simulations. Three examples demonstrate the utility of the approach: phase modulation, birefringence, and frequency mixing.
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.
The solar absorptance of a thermal receiver surface significantly affects the photothermal conversion efficiency of concentrated solar power (CSP) plants. The development of low-cost Gen3 CSP systems calls for increasing the solar absorptance of the thermal receivers at operating temperatures above 750 °C. This study presents an innovative approach to the fractal, multiscale texturing of absorber surfaces to significantly enhance solar energy absorption for Gen3 CSP applications. The absorber surface is described in terms of its fractal parameters that are uniquely determined from surface profile measurements. The interaction of solar radiation with the fractal surface is numerically simulated by solving the governing Maxwell’s equations for electromagnetic wave propagation to investigate the effect of texturing on the absorptance of the surfaces. It is shown that absorption of the solar spectrum increases with an increase in the fractal dimension and the multiscale asperity height of the surface texturing. Here, the theoretical model is demonstrated to be in close agreement with experimental measurements of spectral absorptance of electrodeposited copper (Cu), copper mono-oxide (CuO), and copper-manganese oxide (CuMnO) surfaces that are textured to produce a range of fractal parameters by tailoring the deposition parameters. Fractal surface texturing is shown to reduce reflectance by over one order of magnitude, yielding an absorptance of greater than 0.985 for CuMnO. For the first time, the study presents a significant effective means of increasing solar absorptance and a fundamental theoretical description of the underlying physics.