Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “electromagnetic field calculations”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS2008

Radiation from Relativistic Jets in Turbulent Magnetic Fields

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we have investigated long-term particle acceleration associated with an relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations have been performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. The acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to the afterglow emission. We have calculated the time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique.

Nishikawa, K.-I.↗
NASA NTRS2010

Rapid Variability Generated at Relativistic Shocks Simulated by Particle-in- Cell Code

Plasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electronpositron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electrons transverse deflection behind the shock. We calculate the radiation from deflected electrons in the turbulent magnetic fields. Radiation from electrons near the trailing shock will be variable due to fluctuations of density and electromagnetic fields. The properties of this radiation may be important for rapid variability in relativistic jets such as AGN jets and blazars.

Nishikawal, Ken-Ichi↗
NASA NTRS1986

Electromagnetic instabilities and gyrophase-bunched particles

The linear theory of electromagnetic instabilities at propagation parallel or antiparallel to a uniform magnetic field B(0) is used to calculate the relative phase angle between the fluctuating velocity vector of the jth component and the fluctuating magnetic field. A criterion that determines when this phase relationship leads to observable gyrophase bunching during the linear growth phase of an instability is also derived. The theory is applied to bunching of ion components by ion beam instabilities. Among the electromagnetic cool ion beam instabilities, it is found that the right-hand resonant mode is most likely to cause observable gyrophase bunching of an ion beam during linear growth.

Gary, S. P.↗
NASA NTRS2011

Coherent Backscattering by Polydisperse Discrete Random Media: Exact T-Matrix Results

The numerically exact superposition T-matrix method is used to compute, for the first time to our knowledge, electromagnetic scattering by finite spherical volumes composed of polydisperse mixtures of spherical particles with different size parameters or different refractive indices. The backscattering patterns calculated in the far-field zone of the polydisperse multiparticle volumes reveal unequivocally the classical manifestations of the effect of weak localization of electromagnetic waves in discrete random media, thereby corroborating the universal interference nature of coherent backscattering. The polarization opposition effect is shown to be the least robust manifestation of weak localization fading away with increasing particle size parameter.

Mishchenko, Michael I.↗
NASA NTRS1997

Magnetic Damping of g-Jitter Driven Flows: 3-D Calculations

A 3-D numerical model is developed to represent the oscillating natural convection induced in a cylindrical cavity filled with Ga-doped germanium with and without the presence of an external magnetic field. The model is developed based on the penalty-finite element solution of the equations describing the transport of momentum, heat and solutal element as well as the electromagnetic field distribution in the melt pool. Automatic time step control is applied to help speed up the calculations. Numerical simulations are conducted to study the convection and magnetic damping effects as a function of frequency, directions and amplitudes of g-jitter and also the direction and magnitudes of the applied magnetic fields. The results show that the g-jitter driven flow is time dependent and exhibits a complex recirculating convection pattern in three dimensions and that an applied magnetic field can be employed to suppress this deleterious convective flow and both magnitude and orientation of the applied field are important in magnetic damping of the g-jitter induced convective flows.

Shang, D. Y.↗
DOE OSTI2021

Measurement of e + e - Momentum and Angular Distributions from Linearly Polarized Photon Collisions

The Breit-Wheeler process which produces matter and antimatter from photon collisions is experimentally investigated through the observation of 6085 exclusive electron-positron pairs in ultraperipheral Au + Au collisions at $\sqrt{s_{NN}}$ = 200 GeV. The measurements reveal a large fourth-order angular modulation of cos 4Δ$\phi$ = (16.8 ± 2.5)% and smooth invariant mass distribution absent of vector mesons ($\phi$, ω, and ρ) at the experimental limit of ≤ 0.2% of the observed yields. The differential cross section as a function of e + e - pair transverse momentum P ⟂ peaks at low value with $\sqrt\langle$$P2\atop{⟂}$$\rangle$ = 38.1 ± 0.9 MeV and displays a significant centrality dependence. These features are consistent with QED calculations for the collision of linearly polarized photons quantized from the extremely strong electromagnetic fields generated by the highly charged Au nuclei at ultrarelativistic speed. The experimental results have implications for vacuum birefringence and for mapping the magnetic field which is important for emergent QCD phenomena.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2001

A Comparison of Three-Dimensional Simulations of Traveling-Wave Tube Cold-Test Characteristics Using CST MICROWAVE STUDIO and MAFIA

Previously, it was shown that MAFIA (solutions of Maxwell's equations by the Finite Integration Algorithm), a three-dimensional simulation code, can be used to produce accurate cold-test characteristics including frequency-phase dispersion, interaction impedance, and attenuation for traveling-wave tube (TWT) slow-wave structures. In an effort to improve user-friendliness and simulation time, a model was developed to compute the cold-test parameters using the electromagnetic field simulation software package CST MICROWAVE STUDIO (MWS). Cold-test parameters were calculated for several slow-wave circuits including a ferruled coupled-cavity, a folded waveguide, and a novel finned-ladder circuit using both MWS and MAFIA. Comparisons indicate that MWS provides more accurate cold-test data with significantly reduced simulation times. Both MAFIA and MWS are based on the finite integration (FI) method; however, MWS has several advantages over MAFIA. First, it has a Windows based interface for PC operation, making it very user-friendly, whereas MAFIA is UNIX based. MWS uses a new Perfect Boundary Approximation (PBA), which increases the accuracy of the simulations by avoiding stair step approximations associated with MAFIA's representation of structures. Finally, MWS includes a Visual Basic for Applications (VBA) compatible macro language that enables the simulation process to be automated and allows for the optimization of user-defined goal functions, such as interaction impedance.

Chevalier, C. T.↗
NASA NTRS1991

A theoretical analysis of the electromagnetic environment of the AS330 super Puma helicopter external and internal coupling

Numerical techniques such as Finite Difference Time Domain (FDTD) computer programs, which were first developed to analyze the external electromagnetic environment of an aircraft during a wave illumination, a lightning event, or any kind of current injection, are now very powerful investigative tools. The program called GORFF-VE, was extended to compute the inner electromagnetic fields that are generated by the penetration of the outer fields through large apertures made in the all metallic body. Then, the internal fields can drive the electrical response of a cable network. The coupling between the inside and the outside of the helicopter is implemented using Huygen's principle. Moreover, the spectacular increase of computer resources, as calculations speed and memory capacity, allows the modellization structures as complex as these of helicopters with accuracy. This numerical model was exploited, first, to analyze the electromagnetic environment of an in-flight helicopter for several injection configurations, and second, to design a coaxial return path to simulate the lightning aircraft interaction with a strong current injection. The E field and current mappings are the result of these calculations.

Flourens, F.↗
NASA NTRS1972

The electromagnetic structure of interplanetary space

A method to calculate the three-dimensional structure of the interplanetary magnetic field is presented. The integrations are based on magnetograph recordings of longitudinal magnetic fields in the solar photosphere. The program by Altschuler and Newkirk is used to calculate the radial component of the magnetic field on the source surface, situated at r = 2.6r. This determines the inner boundary conditions for the integration outwards of the interplanetary field equations by means of the method that is described. Computer-drawn plots of the interplanetary field lines out to the earth's orbit are presented for the periods around the total solar eclipses of November 12, 1966, and March 7, 1970. During the former period the interplanetary field exhibited a clean, dipole-type structure, while during the latter period the field was more complicated and had four sectors in the equatorial plane. A movie was presented, showing how the interplanetary field structure rotates as seen by a magnetometer in the sun's equatorial plane.

Stenflo, J. O.↗
DOE OSTI2024

Nature of polariton transport in a Fabry-Perot cavity

Fabry-Perot microcavities can strongly enhance interactions between light and molecules, leading to the formation of hybrid light-matter states known as polaritons. Polaritons possess much smaller effective masses and much larger group velocities when the molecules are resonant with cavity modes that have finite (nonzero) in-plane wave vectors, giving rise to the possibilities of long-range and ultrafast ballistic transport. In this paper, we present the results of numerical simulations of the ultrafast ballistic transport phenomenon in real space and time during and after initialization with a short, spatially localized pulse. Here, we address this problem with two approaches: the standard transfer-matrix method (TMM) for planar structures and a second simulation based on the numerical solution of the Maxwell-Bloch equations that can be used for general configurations (complex metasurfaces, for instance) and excitation modes. The agreement between the TMM and the full numerical calculation with the Maxwell-Bloch equations when applied to multilayer planar configurations provides proof of the validity of both approaches for the present analyses. Overall, we find that the transport of the molecular excitons inside the cavity synchronize with the evolution of the enhanced electromagnetic field inside the cavity. Moreover, the synchronized transport rate is in good agreement with the group velocities predicted from a calculated dispersion relation across a wide range of frequencies. Finally, we relate the group velocity to the Hopfield coefficient obtained from quantum modeling and suggest that the dependence of light-matter coupling on the in-plane wave vector can be an important but overlooked factor for understanding the transport behavior of polaritons. These simulations provide an intuitive tool for understanding the collective motion of light and excitons and helps us to better understand how experimental observations of polaritons should be interpreted.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Gauge-symmetrization method for energy-momentum tensors in high-order electromagnetic field theories

For electromagnetic field theories, canonical energy-momentum conservation laws can be derived from the underpinning spacetime translation symmetry according to the Noether procedure. However, the canonical energy-momentum tensors (EMTs) are neither symmetric nor gauge-symmetric (gauge invariant). The Belinfante-Rosenfeld (BR) method is a well-known procedure to symmetrize the EMTs, which also renders them gauge symmetric for first-order field theories. High-order electromagnetic field theories appear in the study of gyrokinetic systems for magnetized plasmas and the Podolsky system for the radiation reaction of classical charged particles. For these high-order field theories, gauge-symmetric EMTs are not necessarily symmetric and vice versa. In the present study, we develop a new gauge-symmetrization method for EMTs in high-order electromagnetic field theories. The Noether procedure is carried out using the Faraday tensor $F_{μν}$, instead of the 4-potential $A_{μ}$, to derive a canonical EMT $T^{μν}_{N}$. We show that the gauge-dependent part of $T^{μν}_{N}$ can be removed using the displacement-potential tensor $F^{σμν}$ ≡ $D^{σμ}A^{ν}/4π$, where $D^{σμ}$ is the antisymmetric electric displacement tensor. This method gauge-symmetrizes the EMT without necessarily making it symmetric, which is adequate for applications not involving general relativity. For first-order electromagnetic field theories, such as the standard Maxwell system, $F^{σμν}$ reduces to the familiar BR superpotential $S^{σμν}$, and the method developed can be used as a simpler procedure to calculate $S^{σμν}$ without employing the angular momentum tensor in 4D spacetime. When the electromagnetic system is coupled to classical charged particles, the gauge-symmetrization method for EMTs is shown to be effective as well.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Forward doubly-virtual Compton scattering off an unpolarized deuteron in pionless effective field theory

In this work, we calculate the forward unpolarised doubly-virtual Compton scattering (VVCS) off the deuteron in the framework of pionless effective field theory, up to next-to-next-to-next-to-leading order (N3LO) for the longitudinal and next-to-leading order (NLO) for the transverse amplitude. The charge elastic form factor of the deuteron, obtained from the residue of the longitudinal VVCS amplitude, is used to extract the value of the single unknown two-nucleon one-photon contact coupling that enters the longitudinal amplitude at N3LO. We also study the lowest spin-independent generalised polarisabilities of the deuteron. The calculated unpolarised VVCS amplitude provides a high-precision model-independent input for a future calculation of the two-photon-exchange correction to the Lamb shift of muonic deuterium.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1995

Computing Radiation From Multiple-Circular-Aperture Antennas

Mutual Coupling Program for Circular Waveguide-fed Aperture Array (CWG) computer program developed to calculate electromagnetic interactions among elements of antenna arrays with circular apertures with specified distributions of electromagnetic fields in apertures. Distributions assumed superpositions of electromagnetic modes existing in circular waveguide. Various external media included to provide flexibility of use; for example, flexibility to determine effects of dielectric covers upon impedances of aperture-type antennas. Written in VAX FORTRAN for DEC VAX-series computers running VMS (LAR-15236) and IBM PC-series and compatible computers running MS-DOS (LAR-15226).

Bailey, M. C.↗
NASA NTRS1992

Electromagnetic scattering calculations on the Intel Touchstone Delta

During the first year's operation of the Intel Touchstone Delta system, software which solves the electric field integral equations for fields scattered from arbitrarily shaped objects has been transferred to the Delta. To fully realize the Delta's resources, an out-of-core dense matrix solution algorithm that utilizes some or all of the 90 Gbyte of concurrent file system (CFS) has been used. The largest calculation completed to date computes the fields scattered from a perfectly conducting sphere modeled by 48,672 unknown functions, resulting in a complex valued dense matrix needing 37.9 Gbyte of storage. The out-of-core LU matrix factorization algorithm was executed in 8.25 h at a rate of 10.35 Gflops. Total time to complete the calculation was 19.7 h-the additional time was used to compute the 48,672 x 48,672 matrix entries, solve the system for a given excitation, and compute observable quantities. The calculation was performed in 64-b precision.

Cwik, Tom↗
NASA NTRS2003

Large- and Small-Scale Ring Current Electrodynamic Coupling

In this talk we will address the two primary issues of ring current (RC) electrodynamic coupling: 1. RC self-consistent magnetosphere-ionosphere coupling that includes calculation of the magnetospheric electric field (large scale electrodynamic coupling); and 2. RC self-consistent coupling with electromagnetic ion cyclotron (EMIC) waves (small scale electrodynamic coupling). Our study will be based on two RC models that we have recently developed in our group. The first model by Khazanov et al. [2002] couples the system of two kinetic equations: one equation which describes the RC ion dynamics and another equation which describes the energy density evolution of EMIC waves. The second model by Khazanov et al. [2003] deals with large scale electrodynamic coupling processes and provides a self-consistent simulation of RC ions and the magnetospheric electric field. There is presently no model that addresses both of these issues simultaneously in a self-consistent calculation. However, the need exists for such a model, because these two processes directly influence each other, with the mesoscale coupling changing the drift paths of the thermal and energetic particle populations in the inner magnetosphere, thereby changing the wave interactions, and the microscale coupling altering the pitch angle distributions and ionospheric conductivities (through increased precipitation), thus changing the field-aligned currents and electric potential structure. The initial thrust of the work will be the development of a combined kinetic model of micro- and meso-scale RC electrodynamic coupling processes and to examine their interactions with each other on a global scale.

Khazanov, G. V.↗
DOE OSTI2022

Role of electromagnetic gauge-field fluctuations in the selection between chiral and nematic superconductivity

Motivated by the observation of nematic superconductivity in several systems, we revisit the problem of the leading pairing instability of two-component unconventional superconductors on the triangular lattice—such as (p x , p y )-wave and (d x 2 –y 2 , d xy )-wave superconductors. Such a system has two possible superconducting states: the chiral state (e.g., p + ip or d + id), which breaks time-reversal symmetry, and the nematic state (e.g., p + p or d + d), which breaks the threefold rotational symmetry of the lattice. Weak-coupling calculations generally favor the chiral over the nematic superconducting state, raising the question of what mechanism can stabilize the latter. Here, we show that the electromagnetic field fluctuations can play a crucial role in selecting between these two states. Specifically, we derive and analyze the effective free energy for the two-component superconducting order parameter after integrating out the gauge-field fluctuations, which is formally justified if the spatial order parameter fluctuations can be neglected. A nonanalytic cubic term arises, as in the case of a conventional s -wave superconductor. However, unlike the latter, the cubic term depends on the relative phase and on the relative amplitudes between the two order parameter components, in such a way that it generally favors the nematic state. This result is a direct consequence of the fact that the stiffness of the superconducting order parameter is not isotropic. Competition with the quartic term, which favors the chiral state, leads to a renormalized phase diagram in which the nematic state displaces the chiral state over a wide region in the parameter space. Here, we analyze the stability of the fluctuation-induced nematic phase, generalize our results to tetragonal lattices, and discuss their applicability to candidate nematic superconductors, including twisted bilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Electromagnetic signatures of a chiral quantum spin liquid

Quantum spin liquids (QSL) have emerged as a captivating subject within interacting spin systems that exhibit no magnetic ordering even at the lowest temperature accessible experimentally. However, definitive experimental evidence remains elusive. In light of the recent surge in theoretical and experimental interest in the half-filled Hubbard model on a triangular lattice, which offers the potential for stabilizing a chiral QSL, we investigate the electromagnetic signatures of this phase to facilitate experimental detection. Utilizing a combination of parton mean-field theory and unbiased density-matrix renormalization group calculations, we systematically examine the electrical charge and orbital electrical current associated with a spinon excitation in the chiral QSL. Additionally, we calculate the longitudinal and transverse optical conductivities below the Mott gap. Furthermore, employing quantum field theory analysis, we unravel the connection between spinon excitations and emergent as well as physical gauge fields. Our results demonstrate that the chiral QSL phase exhibits a distinct electromagnetic response, even within a Mott insulator regime. This finding holds great potential for enabling the experimental detection of this long-sought-after phase.

36 MATERIALS SCIENCE↗
NASA NTRS2022

Preferential Acceleration of Heavy Ions in a Spontaneously Fragmenting Flare Current Sheet

We study the ion acceleration in a mesoscale, spontaneously fragmenting flare current sheet (SFCS) characterized by the presence of a plasmoid cascade. The main subject of our investigation is to determine whether and how plasmoid cascades at intermediate scales in a fragmented current sheet of a solar flare can impact the (preferential) acceleration of specific ions. The time evolution of the SFCS is obtained from high-resolution 2.5D MHD simulations. The ion trajectories (in the background fields resulting from the MHD model), energies, and pitch angles are calculated using a relativistic test-particle code based on the half-acceleration–rotation–half-acceleration method. For light ions, the main acceleration effects of electromagnetic fields within the SFCS are analyzed using the guiding center approximation. We identify regions with the most-efficient ion acceleration within the SFCS, the accelerator efficiency, and spectra of the accelerated ions. The influence of the charge-to-mass ratio on ion behavior is also studied and resulting ion abundances are compared with observational data. The main ion acceleration takes place in the regions with a strong polarization term, which is part of the first-order Fermi acceleration. Because the term is mass dependent, heavier ions undergo preferential acceleration. The ion energy spectra, abundance-enhancement factors, and differential fluxes, obtained from the model, exhibit power-law profiles, in agreement with observed solar energetic particle events. Nonetheless, the obtained slopes for the abundance-enhancement factor do not exactly match the observed data. The computed slopes and profiles are not sensitive to changes in the initial plasma temperature.

David Kramoliš↗