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.
A theory is presented and a calculational procedure is outlined for evaluating transition amplitudes of multiphoton excitations of vibrational-rotational levels in diatomic molecules. This theory can be utilized in studying behavior of molecules in intense electromagnetic fields.
In this work, we study the different ways of introducing light-matter interaction in first-principle tight-binding (TB) models. The standard way of describing optical properties is the velocity gauge, defined by linear coupling to the vector potential. In finite systems a transformation to represent the electromagnetic radiation by the electric field instead is possible, albeit subtleties arise in periodic systems. The resulting dipole gauge is a multi-orbital generalization of the Peierls substitution. In this work we investigate the accuracy for both pathways, with particular emphasis on gauge invariance, for TB models constructed from maximally localized Wannier functions. This approach accurately captures the light-matter interaction close to the Fermi level. Focusing on paradigmatic two-dimensional materials, we construct first-principle models and calculate the response to electromagnetic fields in linear response and for strong excitations. Benchmarks against fully converged first-principle calculations allow for ascertaining the accuracy of the TB models. We find that the dipole gauge provides a more accurate description than the velocity gauge in all cases. The main deficiency of the velocity gauge is an imperfect cancellation of paramagnetic and diamagnetic current. Formulating a corresponding sum rule however provides a way to explicitly enforce this cancellation. This procedure corrects the TB models in the velocity gauge, yielding excellent agreement with dipole gauge and thus restoring gauge invariance.
Theoretical derivations, computer analysis and test data are provided to demonstrate that the cavity model is a feasible one to analyze thin-substrate, rectangular-patch microstrip antennas. Seven separate antennas were tested. Most of the antennas were designed to resonate at L-band frequencies (1-2 GHz). One antenna was designed to resonate at an S-band (2-4 GHz) frequency of 2.025 GHz. All dielectric substrates were made of Duroid, and were of varying thicknesses and relative dielectric constant values. Theoretical derivations to calculate radiated free space electromagnetic fields and antenna input impedance were performed. MATHEMATICA 2.2 software was used to generate Smith Chart input impedance plots, normalized relative power radiation plots and to perform other numerical manipulations. Network Analyzer tests were used to verify the data from the computer programming (such as input impedance and VSWR). Finally, tests were performed in an anechoic chamber to measure receive-mode polar power patterns in the E and H planes. Agreement between computer analysis and test data is presented. The antenna with the thickest substrate (e(sub r) = 2.33,62 mils thick) showed the worst match to theoretical impedance data. This is anticipated due to the fact that the cavity model generally loses accuracy when the dielectric substrate thickness exceeds 5% of the antenna's free space wavelength. A method of reducing computer execution time for impedance calculations is also presented.
A mathematical representation for the electromagnetic force field and the fluid flow field in a coreless induction furnace is presented. The fluid flow field was represented by writing the axisymmetric turbulent Navier-Stokes equation, containing the electromagnetic body force term. The electromagnetic body force field was calculated by using a technique of mutual inductances. The kappa-epsilon model was employed for evaluating the turbulent viscosity and the resultant differential equations were solved numerically. Theoretically predicted velocity fields are in reasonably good agreement with the experimental measurements reported by Hunt and Moore; furthermore, the agreement regarding the turbulent intensities are essentially quantitative. These results indicate that the kappa-epsilon model provides a good engineering representation of the turbulent recirculating flows occurring in induction furnaces. At this stage it is not clear whether the discrepancies between measurements and the predictions, which were not very great in any case, are attributable either to the model or to the measurement techniques employed.
Proton orbits are calculated in the electromagnetic vacuum field of a magnetic point dipole rotating with its angular velocity omega perpendicular to its dipole moment mu by numerical integration of the Lorentz-Dirac equation. Trajectories are shown and discussed for various initial conditions. A critical surface is shown separating initial positions of protons which finally hit the pulsar in the polar region from those which finally recede to infinity.
This final report covers nine years of research on future tether applications and on the actual flights of the Small Expendable Deployment System (SEDS). Topics covered include: (1) a description of numerical codes used to simulate the orbital and attitude dynamics of tethered systems during station keeping and deployment maneuvers; (2) a comparison of various tethered system simulators; (3) dynamics analysis, conceptual design, potential applications and propagation of disturbances and isolation from noise of a variable gravity/microgravity laboratory tethered to the Space Station; (4) stability of a tethered space centrifuge; (5) various proposed two-dimensional tethered structures for low Earth orbit for use as planar array antennas; (6) tethered high gain antennas; (7) numerical calculation of the electromagnetic wave field on the Earth's surface on an electrodynamically tethered satellite; (8) reentry of tethered capsules; (9) deployment dynamics of SEDS-1; (10) analysis of SEDS-1 flight data; and (11) dynamics and control of SEDS-2.
A paper discusses the need to perform accurate dielectric property measurements on larger sized samples, particularly liquids at microwave frequencies. These types of measurements cannot be obtained using conventional cavity perturbation methods, particularly for liquids or powdered or granulated solids that require a surrounding container. To solve this problem, a model has been developed for the resonant frequency and quality factor of a cylindrical microwave cavity containing concentric cylindrical samples. This model can then be inverted to obtain the real and imaginary dielectric constants of the material of interest. This approach is based on using exact solutions to Maxwell s equations for the resonant properties of a cylindrical microwave cavity and also using the effective electrical conductivity of the cavity walls that is estimated from the measured empty cavity quality factor. This new approach calculates the complex resonant frequency and associated electromagnetic fields for a cylindrical microwave cavity with lossy walls that is loaded with concentric, axially aligned, lossy dielectric cylindrical samples. In this approach, the calculated complex resonant frequency, consisting of real and imaginary parts, is related to the experimentally measured quantities. Because this approach uses Maxwell's equations to determine the perturbed electromagnetic fields in the cavity with the material(s) inserted, one can calculate the expected wall losses using the fields for the loaded cavity rather than just depending on the value of the fields obtained from the empty cavity quality factor. These additional calculations provide a more accurate determination of the complex dielectric constant of the material being studied. The improved approach will be particularly important when working with larger samples or samples with larger dielectric constants that will further perturb the cavity electromagnetic fields. Also, this approach enables the ability to have a larger sample of interest, such as a liquid or powdered or granulated solid, inside a cylindrical container.
Abstract For designing high-field electromagnets, the Lorentz force on coils needs to be computed in order to design suitable support structures, and the inductance should be computed to evaluate the stored energy and dynamics. Also, the magnetic field and its variation inside the conductor is of interest for computing stress and strain, and due to superconducting quench limits. For these force, inductance, energy, and internal field calculations, the coils cannot be naively approximated as infinitesimally thin filaments due to divergences when the source and evaluation points coincide, so more computationally demanding calculations are usually required, resolving the finite cross-section of the conductors. Here, we present a new alternative method that enables the internal magnetic field vector, self-force, and self-inductance to be computed rapidly and accurately within a 1D filament model. The method is applicable to coils for which the curve center-line can have general noncircular shape, as long as the conductor width is small compared to the radius of curvature. This paper extends a previous calculation for circular-cross-section conductors (Hurwitz et al 2024 IEEE Trans. Magn. ) to consider the case of rectangular cross-section. The reduced model is derived by rigorous analysis of the singularity, regularizing the filament integrals such that they match the true high-dimensional integrals at high coil aspect ratio. The new filament model exactly recovers analytic results for a circular coil, and is shown to accurately reproduce full finite-cross-section calculations for a non-planar coil of a stellarator magnetic fusion device. Due to the efficiency of the model here, it is well suited for use inside design optimization.
Understanding microwave-material interactions will help facilitate the utilization of microwave technology in gasification and renewable energy production. In this study, cellulose was used as a model compound to simulate the organic matter in biomass, and its catalytic decomposition under a microwave (MW) field was studied to identify structural changes from the reaction. The study was conducted using a MW source coupled to a fixed-bed gas-flow reactor, mass spectroscopic, and Fourier transform infrared spectroscopy post-reaction analysis. Zeolite 13X was chosen as a microwave absorber to study the catalytic enhancement of the decomposition of cellulose. Density functional theory (DFT) was used to gain insights into the molecular transformations occurring in the presence of a static electric field, which was used to simulate the electric field component of the microwave electromagnetic radiation. Theoretical calculations demonstrated that both the positive and negative portion of the electric field interact with the permanent dipoles of the cellulose effecting the decomposition through the glycosidic bond breaking mechanism. The theoretical result was verified using infrared spectroscopic analysis of the pure cellulose during microwave heating. The theoretical calculations suggest only the positive electric field component may be active for the cellulose decomposition in the presence of Zeolite 13X. Physically mixing Zeolite 13X with the cellulose led to a significant enhancement in the decomposition rate of the glycosidic bond. Zeolite 13X enhanced the glycosidic O-C decomposition at lower MW power (lower temperatures), whereas the O-H functional group required higher MW power (higher temperature) for its decomposition. The DFT study coupled with the reaction studies revealed that the electric field polarizability is dependent upon both the direction and the orientation of the cellulose. Gas products revealed that applying 250 W of MW power led to the production of CO, H 2 , along with some CO 2 , CH 4 , and benzene at 305 °C. Finally, reaction under 500, 750, and 1000 W of power at constant temperature (305 °C) revealed that higher power led to the complete decomposition of cellulose to mostly CO, H 2 , and CO 2 .
Fields and power transfer for electromagnetic waves in hot plasma calculated using Fourier integral and Fourier series methods
Amorphous spherical titanium dioxide (TiO{sub 2}) nanoparticles, decorated gold nanoparticles, are shown to be formed under irradiation of commercial TiO{sub 2} nanopowders dispersed in aqueous solution of hydrogen tetrachloroaurate (HAuCl{sub 4}) by second-harmonic nanosecond pulses of an Nd : YAG laser. It is found that these hybrid nanostructures are formed during laser-induced remelting of the initial nanoparticles, stimulated by gold nanoclusters (providing more efficient absorption of visible laser radiation by titanium dioxide) recovered on their surface. The morphology and chemical composition of the newly formed hybrid nanomaterials are investigated in detail by electron microscopy, Raman spectroscopy, and energy-dispersive analysis. It is shown that the average size and number of gold nanoparticles recovered on the TiO{sub 2} surface can be controlled by varying the concentration ratio of the initial nanomaterial and HAuCl{sub 4}. A spectroscopic analysis of light scattering by single hybrid nanoparticles and the results of numerical calculation of the structure of electromagnetic fields near their surface indicate good prospects of these hybrid nanomaterials for various applications in modern optics, optoelectronics, and nanophotonics; for example, for designing chemo- and biosensor platforms and new-generation solar cells. (nanoparticles)
Thin metal plate surface impedance during excitation by HF electromagnetic field as function of magnetic field, calculating line shapes
Determining transient electromagnetic fields in antennas with nonlinear loads is a challenging problem. Typical methods used involve calculating frequency domain parameters at a large number of different frequencies, then applying Fourier transform methods plus nonlinear equation solution techniques. If the antenna is simple enough so that the open circuit time domain voltage can be determined independently of the effects of the nonlinear load on the antennas current, time stepping methods can be applied in a straightforward way. Here, transient fields for antennas with more general geometries are calculated directly using Finite Difference Time Domain (FDTD) methods. In each FDTD cell which contains a nonlinear load, a nonlinear equation is solved at each time step. As a test case, the transient current in a long dipole antenna with a nonlinear load excited by a pulsed plane wave is computed using this approach. The results agree well with both calculated and measured results previously published. The approach given here extends the applicability of the FDTD method to problems involving scattering from targets, including nonlinear loads and materials, and to coupling between antennas containing nonlinear loads. It may also be extended to propagation through nonlinear materials.
A formula is proposed for calculating the rate of field ionisation of an atom or ion, taking into account both tunnel ionisation and barrier suppression ionisation. Compared with the previous formula proposed in 2018, it more accurately describes the transition region between both ionisations and is determined mainly by two parameters: the ionisation potential of an atom or ion and the amplitude of an external electric field. This makes the presented formula suitable for use in numerical packages simulating the interaction of high-power laser radiation with matter by the particle-in-cell (PIC) method. (paper)
The measurement of the anomalous electron magnetic moment g – 2 through quantum transitions of a single trapped electron is the most stringent test of quantum field theory. These experiments are now so precise that they must account for the effects of the cavity containing the electron. Classical calculations of this “cavity shift” must subtract the electron’s divergent self-field, and thus require knowledge of the exact Green’s function for the cavity’s electromagnetic field. We perform the first fully quantum calculation of the cavity shift in a closed cavity, which instead involves subtracting linearly divergent cavity mode sums and integrals. Using contour integration methods, we find perfect agreement with existing classical results for both spherical and cylindrical cavities, justifying their current use. Moreover, our mode-based results can be naturally generalized to account for systematic effects, necessary to push future measurements to the next order of magnitude in precision.
The Magnetospheric Multiscale mission (MMS) consists of four identical spacecraft forming a closely separated (less than or equal to 10 km) and nearly regular tetrahedron. This configuration enables the decoupling of spatial and temporal variations and allows the calculation of the spatial gradients of plasma and electromagnetic field quantities. We make full use of the well cross-calibrated MMS magnetometers and fast plasma instruments measurements to calculate both the magnetic and plasma forces in flux transfer events (FTEs) and evaluate the relative contributions of different forces to the magnetopause momentum variation. This analysis demonstrates that some but not all FTEs, consistent with previous studies, are indeed force-free structures in which the magnetic pressure force balances the magnetic curvature force. Furthermore, we contrast these events with FTE events that have non-force-free signatures.
The RMF (Rotating Magnetic Field) code is designed to calculate the motion of a charged particle in a given electromagnetic field. It integrates Hamilton’s equations in cylindrical coordinates using an adaptive predictor-corrector double-precision variable-coefficient ordinary differential equation solver for speed and accuracy. RMF has multiple capabilities for the field. Particle motion is initialized by specifying the position and velocity vectors. Here, the six-dimensional state vector and derived quantities are saved as functions of time. A post-processing graphics code, XDRAW, is used on the stored output to plot up to 12 windows of any two quantities using different colors to denote successive time intervals. Multiple cases of RMF may be run in parallel and perform data mining on the results. Recent features are a synthetic diagnostic for simulating the observations of charge-exchange-neutral energy distributions and RF grids to explore a Fermi acceleration parallel to static magnetic fields.
Numerical simulations of electromagnetic instabilities, driven by a cool tenuous ion beam propagating along an ambient magnetic field, have been conducted in one and two spatial dimensions. The calculations employ particle ions, fluid electrons, and a predictor corrector scheme for solving the electromagnetic field in two dimensions that is described in some detail. While the principal features of the one-dimensional calculations (which reproduce previous work) are retained, the two-dimensional simulations show some reduction of the overall level of the magnetic field fluctuations. Enhancement of the heating of the beam ions at the expense of the core ions also occurs in the case where the beam density is sufficiently large that the right-hand nonresonant instability dominates. Implications of the results for modeling of the ion foreshock and quasi-parallel shocks are discussed.