Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ELECTROMAGNETIC FIELD”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Application of mesh refinement to relativistic magnetic reconnection

During relativistic magnetic reconnection, antiparallel magnetic fields undergo a rapid change in topology, releasing a large amount of energy in the form of non-thermal particle acceleration. This work explores the application of mesh refinement to 2D reconnection simulations to efficiently model the inherent disparity in length-scales. We have systematically investigated the effects of mesh refinement and determined necessary modifications to the algorithm required to mitigate non-physical artifacts at the coarse–fine interface. We have used the ultrahigh-order pseudo-spectral analytical time-domain Maxwell solver to analyze how its use can mitigate the numerical dispersion that occurs with the finite-difference time-domain (or “Yee”) method. Absorbing layers are introduced at the coarse–fine interface to eliminate spurious effects that occur with mesh refinement. We also study how damping the electromagnetic fields and current density in the absorbing layer can help prevent the non-physical accumulation of charge and current density at the coarse–fine interface. Using a mesh refinement ratio of 8 for two-dimensional magnetic reconnection simulations, we obtained good agreement with the high-resolution baseline simulation, using only 36% of the macroparticles and 71% of the node-hours needed for the baseline. The methods presented here are especially applicable to 3D systems where higher memory savings are expected than in 2D, enabling comprehensive, computationally efficient 3D reconnection studies in the future.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Numerical simulation for the design of induction heating based radio frequency reactor for ethylene production

Ethylene is a vital petrochemical compound produced in vast amounts yearly by manufacturers that have enough scale to overcome the inherent thermodynamic inefficiencies of the process. In order to address the inefficiencies that prevent smaller scale or intermittent production ethylene, investigation of new production methods are required. In this work, we investigate the use of a radio frequency (RF) based reactor system that generates heat internally as opposed to applying heat externally via steam or direct combustion of fossil fuels. In order to guide the design of an electromagnetic based reactor system, we have created a macroscale model capable of capturing heat generation at the susceptors from a produced electromagnetic field and subsequent chemical reactions. Several susceptor and induction coil designs were investigated to better understand the performance of the system.

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

Imaging Photonic Resonances within an All‐Dielectric Metasurface via Photoelectron Emission Microscopy

Dielectric metasurfaces, through volume‐type photonic resonances, enable precise control of light‐matter interactions for applications including imaging, holography, and sensing. The application space of dielectric metasurfaces has extended from infrared to visible wavelengths by incorporating high refractive index materials, such as titanium dioxide (TiO 2 ). Understanding the fundamental and fabrication limits for these applications requires metrology with nanoscale resolution, sensitivity to electromagnetic fields within the meta‐atom volume, and far‐field excitation. In this work, photoelectron emission microscopy (PEEM) is used to image field distributions of photonic resonances in a TiO 2 metasurface excited with far‐field, visible‐wavelength illumination. The local volumetric field variations within the meta‐atoms are analyzed as a function of illumination angle and polarization by comparing photoelectron images to finite‐difference time‐domain simulations. This study determines the inelastic mean free path of very low‐energy (<1 eV) photoelectrons to be 35 ± 10 nm, which is comparable to the meta‐atom height thereby highlighting PEEM sensitivity to resonances within the volume. Additionally, the simulations reveal high sensitivity of PEEM images to an in‐plane component of the illumination k ‐vector. These results demonstrate that photoelectron imaging with subwavelength resolution offers unique advantages for examining light‐matter interactions in volume‐type (as opposed to surface) photonic modes within dielectric nanophotonic structures.

electron inelastic mean free path↗
DOE OSTI2024

Electrobaryonic axion: hair of neutron stars

Axion-like particles are predicted in many physics scenarios beyond the Standard Model (SM). Their interactions with SM particles may arise from the triangle anomaly of the associated global symmetry, along with other SM global and gauge symmetries, including anomalies with the global baryon number and electromagnetic gauge symmetries. We initiate the phenomenological study of the corresponding “electrobaryonic axion”—a particle that couples with both the baryon chemical potential and the electromagnetic field. Neutron stars, particularly magnetars, possessing high baryon density and strong magnetic fields, can naturally develop a thin axion hair around their surface. In this study, we calculate this phenomenon, considering the effects of neutron star rotation and general relativity. For axion particles lighter than the neutron star rotation frequency, the anomalous interaction can also induce the emission of axion particles from the neutron star. In the light axion regime, this emission can significantly contribute to the neutron star cooling rate.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part II: Enforcing the Lorenz Gauge Condition

In a previous paper Christlieb et al. (A particle-in-cell method for plasmas with a generalized momentum formulation, part I: Model formulation, 2024), we developed a new particle-in-cell (PIC) method for the relativistic Vlasov–Maxwell system in which the electromagnetic fields and the equations of motion for the particles were cast in terms of scalar and vector potentials through a Hamiltonian formulation. This new method evolved the potentials under the Lorenz gauge using integral equation methods. New methods to construct spatial derivatives of the potentials that converge at the same rates as the fields were also presented. The new particle method was compared against standard explicit discretizations, including the well-known FDTD-PIC method, for a range of applications involving sheaths and particle beams. Here, this paper extends this new class of methods by focusing on the enforcement the Lorenz gauge condition in both exact and approximate forms using co-located meshes. A time-consistency property of the proposed field solver for the vector potential form of Maxwell’s equations is established, which is shown to preserve the equivalence between the semi-discrete Lorenz gauge condition and the analogous semi-discrete continuity equation. Using this property, we present three methods to enforce a semi-discrete gauge condition. The first method introduces an update for the continuity equation that is consistent with the discretization of the Lorenz gauge condition. Both the finite difference and spectral implementations satisfy this discrete gauge condition to machine precision. The second approach we propose enforces a semi-discrete continuity equation using the boundary integral solution to the field equations. The potential benefit of this approach is that it eliminates spatial derivatives that appear on the particle data, namely the current density, which is often calculated by linear combinations of low-order spline basis functions. This method is ideally suited to boundary integral equation methods that invert multi-dimensional operators without dimensional splitting techniques and will be the subject of future work. The third approach introduces a gauge correcting method that makes direct use of the gauge condition to modify the scalar potential and uses local maps for both the charge and current densities. This results in a gauge error, as the maps do not enforce the continuity equation. The vector potential coming from the current density is taken to be exact, and using the Lorenz gauge, we compute a correction to the scalar potential that makes the two potentials satisfy the gauge condition. This method also enforces the gauge condition to machine precision. We demonstrate two of the proposed methods in the context of periodic domains. Problems defined on bounded domains, including those with complex geometric features remain an ongoing effort. However, this work shows that it is possible to design computationally efficient methods that can effectively enforce the Lorenz gauge condition in a non-staggered PIC formulation.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Nyström type exponential integrators for strongly magnetized charged particle dynamics

Solving for charged particle motion in electromagnetic fields (i.e. the particle pushing problem) is a computationally intensive component of particle-in-cell (PIC) methods for plasma physics simulations. This task is especially challenging when the plasma is strongly magnetized due numerical stiffness arising from the wide range of time scales between highly oscillatory gyromotion and long term macroscopic behavior. A promising approach to solve these problems is by a class of methods known as exponential integrators that can solve linear problems exactly and are A-stable. This work extends the standard exponential integration framework to derive Nyström-type exponential integrators that integrates the Newtonian equations of motion as a second-order differential equation directly. In particular, we derive second-order and third-order Nyström-type exponential integrators for strongly magnetized particle pushing problems. Numerical experiments show that the Nyström-type exponential integrators exhibit significant improvement in computation speed over the standard exponential integrators.

general physics↗
DOE OSTI2025

An explicit, energy-conserving particle-in-cell scheme

We present an explicit temporal discretization of particle-in-cell schemes for the non-relativistic Vlasov equation that results in exact energy conservation when combined with an appropriate spatial discretization. The scheme is inspired by a simple, second-order explicit scheme that conserves energy exactly in the Eulerian context. We show that direct translation to particle-in-cell does not result in strict conservation, but derive a simple correction based on an analytically solvable optimization problem that recovers conservation. While this optimization problem is not guaranteed to have a real solution for every particle, we provide a correction that makes imaginary values extremely rare and still admits $\mathcal{O}$(10 –12 ) fractional errors in energy for practical simulation parameters. We present the scheme in both electrostatic – where we use the Ampère formulation – and electromagnetic contexts. With an electromagnetic field solve, the field update is most naturally linearly implicit, but the more computationally intensive particle update remains fully explicit. Here, we also show how the scheme can be extended to use the fully explicit leapfrog and pseudospectral analytic time-domain (PSATD) field solvers. The scheme is tested on standard kinetic plasma problems, confirming its conservation properties.

Energy conservation↗
DOE OSTI2025

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗
DOE OSTI2025

No evidence of a reduction in the dipole polarisability at shell closures

Photoabsorption reactions have a long and productive history as probes of the response of nuclei to electromagnetic fields. In recent years, there has been a suggestion that the dipole polarisability may be reduced at shell closures, potentially due to the reduced level density. In the current paper, more recent photoabsorption strengths are reanalysed including nuclei on closed proton shells to search for a similar reduction. The previously reported reduction is not reproduced. We conclude that the previously observed reduction may be due to incomplete evaluations of the photoabsorption due to missing cross section contributions and experimental artefacts. More experimental data are necessary to draw a firm conclusion about whether shell effects persist into high excitation energies, especially data which can obviate any systematic effects and look at relative photoabsorption strengths between nuclei.

Giant dipole response↗
DOE OSTI2024

Plasmons in the Kagome metal CsV 3 Sb 5

Plasmon polaritons, or plasmons, are coupled oscillations of electrons and electromagnetic fields that can confine the latter into deeply subwavelength scales, enabling novel polaritonic devices. While plasmons have been extensively studied in normal metals or semimetals, they remain largely unexplored in correlated materials. In this paper, we report infrared (IR) nano-imaging of thin flakes of CsV3Sb5, a prototypical layered Kagome metal. We observe propagating plasmon waves in real-space with wavelengths tunable by the flake thickness. From their frequency-momentum dispersion, we infer the out-of-plane dielectric function that is generally difficult to obtain in conventional far-field optics, and elucidate signatures of electronic correlations when compared to density functional theory (DFT). We propose correlation effects might have switched the real part of ϵ c from negative to positive values over a wide range of middle-IR frequencies, transforming the surface plasmons into hyperbolic bulk plasmons, and have dramatically suppressed their dissipation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Towards higher frequencies in a compact prebunched waveguide THz-FEL

Free-electron-lasers fill a critical gap in the space of THz-sources as they can reach high average and peak powers with spectral tunability. Using a waveguide in a THz FEL significantly increases the coupling between the relativistic electrons and electromagnetic field enabling large amounts of radiation to be generated in a single passage of electrons through the undulator. In addition to transversely confining the radiation, the dispersive properties of the waveguide critically affect the velocity and slippage of the radiation pulse which determine the central frequency and bandwidth of the generated radiation. In this paper, we characterize the spectral properties of a compact waveguide THz FEL including simultaneous lasing at two different frequencies and demonstrating tuning of the radiation wavelength in the high frequency branch by varying the beam energy and ensuring that the electrons injected into the undulator are prebunched on the scale of the resonant radiation wavelength.

Experimental particle physics↗
DOE OSTI2026

Purcell enhancement of directional edge photocurrent in a van der Waals self-cavity

Cavities provide a means to manipulate the optical and electronic responses of quantum materials by selectively enhancing light-matter interaction at specific frequencies and momenta. While cavities typically involve external structures, exfoliated flakes of van der Waals (vdW) materials can form intrinsic self-cavities due to their small finite dimensions, confining electromagnetic fields into plasmonic cavity modes, characterized by standing-wave current distributions. While cavity-enhanced phenomena are well-studied at optical frequencies, the impact of self-cavities on nonlinear electronic responses—such as directional photocurrent—remains largely unexplored, particularly in the terahertz regime, critical for emerging ultrafast optoelectronic technologies. Here, we report a self-cavity-induced Purcell enhancement of directional photocurrents in the vdW semimetal WTe 2 . Using ultrafast optoelectronic circuitry, we measured coherent near-field THz emission resulting from nonlinear photocurrents excited at the sample edges. We observed enhanced emission at finite frequencies, tunable via excitation fluence and sample geometry, which we attribute to plasmonic interference effects controlled by the cavity boundaries. We developed an analytical theory that captures the cavity resonance conditions and spectral response across multiple devices. Our findings establish WTe 2 as a bias-free, geometry-tunable THz emitter and demonstrate the potential of self-cavity engineering for controlling nonlinear, nonequilibrium dynamics in quantum materials.

condensed-matter physics↗
DOE OSTI2025

Symmetry-protected electronic metastability in an optically driven cuprate ladder

Optically excited quantum materials exhibit non-equilibrium states with remarkable emergent properties, but these phenomena are usually transient, decaying on picosecond timescales and limiting practical applications. Advancing the design and control of non-equilibrium phases requires the development of targeted strategies to achieve long-lived, metastable phases. Here, in this study, we report the discovery of symmetry-protected electronic metastability in the model cuprate ladder Sr 14 Cu 24 O 41 . Using femtosecond resonant X-ray scattering and spectroscopy, we show that this metastability is driven by a transfer of holes from chain-like charge reservoirs into the ladders. This ultrafast charge redistribution arises from the optical dressing and activation of a hopping pathway that is forbidden by symmetry at equilibrium. Relaxation back to the ground state is, hence, suppressed after the pump coherence dissipates. Our findings highlight how dressing materials with electromagnetic fields can dynamically activate terms in the electronic Hamiltonian, and provide a rational design strategy for non-equilibrium phases of matter.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Emerging supersolidity in photonic-crystal polariton condensates

A supersolid is a counter-intuitive phase of matter in which its constituent particles are arranged into a crystalline structure, yet they are free to flow without friction. This requires the particles to share a global macroscopic phase while being able to reduce their total energy by spontaneous, spatial self-organization. The existence of the supersolid phase of matter was speculated more than 50 years ago1, 2, 3–4. However, only recently has there been convincing experimental evidence, mainly using ultracold atomic Bose–Einstein condensates (BECs) coupled to electromagnetic fields. There, various guises of the supersolid were created using atoms coupled to high-finesse cavities5,6, with large magnetic dipole moments7, 8, 9, 10, 11, 12–13, and spin–orbit-coupled, two-component systems showing stripe phases14, 15–16. Here we provide experimental evidence of a new implementation of the supersolid phase in a driven-dissipative, non-equilibrium context based on exciton–polaritons condensed in a topologically non-trivial, bound state in the continuum (BiC) with exceptionally low losses, realized in a photonic-crystal waveguide. We measure the density modulation of the polaritonic state indicating the breaking of translational symmetry with a precision of several parts in a thousand. Direct access to the phase of the wavefunction allows us to also measure the local coherence of the supersolid. We demonstrate the potential of our synthetic photonic material to host phonon dynamics and a multimode excitation spectrum.

Trypogeorgos, Dimitrios↗
DOE OSTI2025

X-ray free-electron lasing in a flying-focus undulator

Laser-driven free-electron lasers (LDFELs) replace magnetostatic undulators with the electromagnetic fields of a laser pulse. Because the undulator period is half the wavelength of the laser pulse, LDFELs can amplify X rays using lower electron energies and over shorter interaction lengths than a traditional free-electron laser. In LDFELs driven by conventional laser pulses, the undulator uniformity required for high gain necessitates large laser-pulse energies. Here, we show that a flying-focus pulse provides the undulator uniformity required to reach high gain with a substantially lower energy than a conventional pulse. The flying-focus pulse features an intensity peak that travels in the opposite direction of its phase fronts. This enables an LDFEL configuration where an electron beam collides head-on with the phase fronts and experiences a near-constant undulator strength as it co-propagates with the intensity peak. Three-dimensional simulations of this configuration demonstrate the generation of megawatts of coherent X-ray radiation with 20 × less energy than a conventional laser pulse.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A cavity-enhanced spectroscopist's lens on molecular polaritons

Polariton chemistry has been hailed as a potential new route to direct molecular processes with electromagnetic fields. To make further strides, the community needs to clarify which unusual polaritonic phenomena are true hallmarks of cavity quantum electrodynamics and which can be rationalized with classical optical physics. Here, in this work, we provide a tutorial perspective on the formation, spectroscopy, and behavior of molecular polaritons using classical optics. Where possible, we draw connections to cavity-enhanced spectroscopy and recast open questions in terms that may be more familiar to the broader community of physical chemists.

Cavity quantum electrodynamics↗
DOE OSTI2024

Towards modelling AR Sco: generalized particle dynamics and strong radiation-reaction regimes

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.

79 ASTRONOMY AND ASTROPHYSICS↗