Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electromagnetic Simulation”

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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Systematic Construction of Time-Dependent Hamiltonians for Microwave-Driven Josephson Circuits

Time-dependent electromagnetic drives are fundamental for controlling complex quantum systems, including superconducting Josephson circuits. In these devices, accurate time-dependent Hamiltonian models are imperative for predicting their dynamics and designing high-fidelity quantum operations. Existing numerical methods, such as black-box quantization (BBQ) and energy-participation ratio (EPR), excel at modeling the static Hamiltonians of Josephson circuits. However, these techniques do not fully capture the behavior of driven circuits stimulated by external microwave drives, nor do they include a generalized approach to account for the inevitable noise and dissipation that enter through microwave ports. Here, we introduce numerical techniques that leverage classical microwave simulations, efficiently executable in finite-element solvers, to obtain the time-dependent Hamiltonian of microwave-driven superconducting circuits with arbitrary geometries under charge, flux, or mixed electromagnetic modulation. Importantly, our techniques do not rely on a lumped-element description of the superconducting circuit, in contrast to previous approaches to tackling this problem. We demonstrate the versatility of our approach by characterizing the driven properties of realistic circuit devices in complex electromagnetic environments, including coherent dynamics due to charge and flux modulation, as well as drive-induced relaxation and dephasing. Our techniques offer a powerful toolbox for optimizing circuit designs and advancing practical applications in superconducting quantum computing.

Lu, Yao [Yale U.; Yale U. (main); Fermilab] (ORCID↗
DOE OSTI2025

Bounce-averaged theory in arbitrary multi-well plasmas: solution domains and the graph structure of their connections

Bounce-averaged theories provide a framework for simulating relatively slow processes, such as collisional transport and quasilinear diffusion, by averaging these processes over the fast periodic motions of a particle on a closed orbit. This procedure dramatically increases the characteristic time scale and reduces the dimensionality of the modelled system. The natural coordinates for such calculations are the constants of motion (COM) of the fast particle motion, which by definition do not change during an orbit. However, for sufficiently complicated fields – particularly in the presence of local maxima of the electric potential and magnetic field – the COM are not sufficient to specify the particle trajectory. In such cases, multiple domains in COM space must be used to solve the problem, with boundary conditions enforced between the domains to ensure continuity and particle conservation. Previously, these domains have been imposed by hand, or by recognising local maxima in the fields, limiting the flexibility of bounce-averaged simulations. Here, we present a general set of conditions for identifying consistent domains and the boundary condition connections between the domains, allowing the application of bounce-averaged theories in arbitrarily complicated and dynamically evolving electromagnetic field geometries. We also show how the connections between the domains can be represented by a directed graph, which can help to succinctly represent the trajectory bifurcation structure.

fusion plasma↗
DOE OSTI2020

Solitary Magnetic Structures at Quasi–Parallel Collisionless Shocks: Formation

Solitary magnetic structures known as SLAMS (short large-amplitude magnetic structures) have been considered as essential elements of collisionless shocks with quasi-parallel geometries. Yet the physics underlying their formation remains an open question. In this paper, we use measurements from the magnetospheric multiscale mission combined with fully kinetic simulations to study the formation of SLAMS. We find that gyro-resonance between solar wind ions and right-hand circularly polarized electromagnetic waves results in magnetic field amplification. Gyro-trapping by the growing magnetic field builds up the plasma density that further enhances the current and field. The solitary nature of SLAMS stems from a beat-like magnetic field envelope where the maximum sets the initial location for nonlinear growth. Our results present a conceptual advance on SLAMS, and may shed new light on the open question of magnetic field amplification at astrophysical shocks.

58 GEOSCIENCES↗
DOE OSTI2026

Calibrated particle identification for Belle II

We present several efforts aimed at improving charged particle identification at the Belle II experiment. We define an ablation test to quantify and evaluate the impact of each sub-detector on the global particle identification performance. We demonstrate that the performance of the identification scheme can be improved via a simple calibration of the sub-detector likelihoods through a set of per hypothesis, per sub-detector weights. Finally, we present preliminary results on an improved definition of the likelihoods contributed by the electromagnetic calorimeter. A set of multiclass boosted decision trees is trained to exploit the shape of energy depositions of different particle species. In simulated $B\bar{B}$ events, the pion-to-electron and muon fake rates are reduced by 55% and 31% respectively at low-medium momentum.

Hohmann, Marcel [Univ. of Melbourne, Parkville, VI↗
DOE OSTI2022

Localized Excitation of Electromagnetic Ion Cyclotron Waves From Anisotropic Protons Filtered by Magnetic Dips

The excitation of electrostatic and/or electromagnetic waves in the plasma universe is often associated with anisotropic velocity distributions of charged particles. In Earth's inner magnetosphere, this anisotropy can gradually develop as particles injected from the magnetotail drift around the Earth at different speeds depending on their energy and pitch angle. Here, we show that the perpendicular-moving and bouncing ions can be separated more abruptly near the injection front. These pitch-angle filters are localized magnetic dip structures formed by the diamagnetic behavior of the injected particles, which can trap perpendicular-moving ions and allow bouncing ions to overtake. The resulting ion anisotropy facilitates the rapid generation of electromagnetic ion cyclotron (EMIC) waves, which in turn can largely reshape the Van Allen radiation belts. This scenario is examined by case and statistical observations, together with numerical simulations that reproduce most of the observational signatures, to support the causal relationship between magnetic dips, anisotropic ion distributions, and localized excitation of EMIC waves. Finally, our study highlights the important roles of magnetic dips in the inner magnetospheric dynamics, as pitch-angle filters of the injected ions and traveling hotspots of EMIC wave activities.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Celeritas Midterm SciDAC Report

Celeritas is a new Monte Carlo (MC) code that helps satisfy the increasing demand for high energy physics (HEP) detector simulation, using Graphics Processing Unit (GPU) hardware on high performance computing (HPC) systems to model Large Hadron Collider (LHC) experiments and beyond. This report details the project’s progress midway through its SciDAC funding period, highlighting the first complete implementation of standard electromagnetic (EM) physics on GPUs, initial results for performance and scalability on Leadership Computing Facilities (LCFs), and preliminary integration into the CMS and ATLAS experiments. By integrating HEP domain knowledge with expertise in MC transport, Celeritas has catalyzed a shift in the HEP community’s perception of GPU platforms as the future for HPC simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Ponderomotive electron physics captured in a single-fluid extended magnetohydrodynamics model

The ponderomotive force, arising from the interaction between electromagnetic waves and plasma, plays a critical role in laser fusion, astrophysical plasmas, and laser diagnostics. Traditionally, modeling this force requires multi-fluid or particle-in-cell simulations due to its strong coupling to electron-scale dynamics. In this work, we demonstrate that a one-fluid, two-temperature extended magnetohydrodynamics (XMHD) model—augmented with a generalized Ohm's law (GOL) including electron inertia—can accurately reproduce key ponderomotive effects. We derive the ponderomotive force within this framework using a phasor-based approach and then validate its nonlinear manifestations through direct numerical simulations in the PERSEUS code, where steepening and density modulation phenomena typically associated with kinetic-scale models are reproduced. These results and prior work establish XMHD as a robust and efficient alternative for modeling nonlinear laser–plasma dynamics, bridging the gap between ideal MHD and fully kinetic approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Effect of Density Irregularities on Radio Frequency Wave Propagation in Ionospheric Plasmas

Density irregularities play a vital role in determining how radio frequency (RF) waves travel through plasmas. In the Earth’s ionosphere, these density irregularities also impact radio communication. In this study, we conduct a detailed numerical analysis of RF wave propagation in small-scale ionospheric density irregularities using the advanced Petra-M code. We focus specifically on high-frequency (HF) waves, ranging from 3 to 30 MHz, which are essential for military, amateur radio operators, and emergency communications. By introducing density structures, such as equatorial plasma bubbles derived from fluid simulations, we demonstrate that HF waves can scatter in multiple directions when they encounter these irregularities. Additionally, we observe significant mode conversion, where incoming electromagnetic waves transform into electrostatic modes within the density gradient layer. This shows that smaller density irregularities can greatly weaken signals or cause complete signal loss for receivers, emphasizing the need for increased awareness and innovative solutions in radio communication transmission.

Kim, Eun-Hwa [Princeton Plasma Physics Laboratory ↗
DOE OSTI2020

Direct electron acceleration in multi-kilojoule, multi-picosecond laser pulses

The physics mechanism behind the acceleration of electrons to energies much higher than the laser ponderomotive potential is reported. While electrons with these energies have been observed for the last twenty years at the one percent level, it is only with the recent advent of high power lasers with extended pulse durations, large focal spots, and high intensities that significant numbers (>10%) of the so-called “superponderomotive” electrons have been observed. Evidence is provided that the dominant acceleration mechanism of superponderomotive electrons is direct acceleration by the electromagnetic (EM) wave of the laser as they co-propagate in an under-dense plasma; one-dimensional particle-in-cell simulations yield electron spectra with typically ten times the ponderomotive temperature. Direct acceleration requires tens of plasma wavelengths of under-dense plasma, as well as pulse durations of tens of plasma periods, conditions that can be found in large focal spot multi-picosecond laser pulses. It is shown that the plasma simply provides the source of accelerated electrons and a variety of methods to inject a significant population of electrons into the EM waves at relativistic energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Wavelength-tunable infrared chiral metasurfaces with phase-change materials

Optical phase-change materials exhibit tunable permittivity and switching properties during phase transition, which offers the possibility of dynamic control of optical devices. Here, a wavelength-tunable infrared chiral metasurface integrated with phase-change material GST-225 is demonstrated with the designed unit cell of parallelogram-shaped resonator. By varying the baking time at a temperature above the phase transition temperature of GST-225, the resonance wavelength of the chiral metasurface is tuned in the wavelength range of 2.33 µm to 2.58 µm, while the circular dichroism in absorption is maintained around 0.44. The chiroptical response of the designed metasurface is revealed by analyzing the electromagnetic field and displacement current distributions under left- and right-handed circularly polarized (LCP and RCP) light illumination. Moreover, the photothermal effect is simulated to investigate the large temperature difference in the chiral metasurface under LCP and RCP illumination, which allows for the possibility of circular polarization-controlled phase transition. The presented chiral metasurfaces with phase-change materials offer the potential to facilitate promising applications in the infrared regime, such as chiral thermal switching, infrared imaging, and tunable chiral photonics.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Radiation and Polarization Signatures from Magnetic Reconnection in Relativistic Jets. I. A Systematic Study

Blazars are relativistic magnetized plasma outflows from supermassive black holes that point very close to our line of sight. Their emission is nonthermal-dominated and highly variable across the entire electromagnetic spectrum. Relativistic magnetic reconnection has been proposed as the driver of particle acceleration during blazar flares. While recent particle-in-cell (PIC) simulations have self-consistently studied the evolution of magnetic reconnection and particle acceleration therein, the resulting radiation signatures have not been systematically explored. In particular, the polarization signatures, which directly reflect the characteristic strongly dynamical evolution of magnetic field during reconnection, have not been carefully investigated. In this paper, we present a systematic study of radiation and polarization signatures arising from magnetic reconnection in blazars, based on combined PIC and polarized radiation transfer simulations with various physical parameters. We identify a harder-when-brighter trend in the spectral evolution. Moreover, higher-frequency bands (ultraviolet to X-ray) tend to flare earlier than lower-frequency bands (infrared to optical) in the synchrotron spectral component. Most importantly, polarization signatures appear more variable with higher frequencies. We find that the variation in temporal polarization depends strongly on the guide field strength. Specifically, reconnection with a significant guide field component leads to a very high polarization degree that contradicts typical blazar observations, while large polarization angle rotations are unique signatures of magnetic reconnection between nearly antiparallel magnetic field lines. These rotations are at least 90° and can extend to >180°, and they may be in either direction. Furthermore, these results imply that blazars that have shown large polarization angle rotations intrinsically have more nearly antiparallel magnetic field morphology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Simulations of plasma current induced by toroidal field ramping down on tokamaks

In tokamak experiments, the BT ramping down can drive off-axis and parallel inductive current. This approach leads to a decrease in q95 and an increase in normalized beta, βN. The off-axis inductive current also broadens the current profile. Typically, the 1.5D transport code is used to simulate the time evolution of the plasma current profile, which is based on the flux-surface-averaged Faraday’s law. The ONETWO code is one of those transport codes. However, this code cannot simulate the situation of the evolving toroidal field BT. This article proposes a new time-dependent model to take into account the BT ramping down situation. A modified formula of flux-surface-averaged Faraday’s law was derived to consider the effect of the BT ramping rate on the current evolution, and it was implemented in the ONETWO code. Then, the modified ONETWO code was used to simulate the BT ramping down experiment on DIII-D. The simulation result of the plasma current evolution with BT ramping down shows a broader current profile with smaller ohmic current induced by the poloidal field, compared to that without BT ramping down.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Port Impedance Measurement and Current Injection Response Analysis for PLCs

Programmable Logic Controllers (PLCs) are used to control devices throughout the power system since they have fast control capabilities and can utilize multiple types of communication interfaces. Therefore, studying and mitigating their vulnerabilities to electromagnetic pulse is important for the reliability of PLC operations. Here, in this paper, an effective impedance measurement scheme is proposed and demonstrated for three PLCs to estimate their susceptibility to an electromagnetic pulse. The equivalent non-uniform transmission line model is established to eliminate the impact of the fixture in the de-embedding process. Then different parameters of the impedance measurement setup are explored. Based on the measured impedance, the equivalent circuit is established to calculate the response of the device when subjected to the electromagnetic pulse. The voltage and current responses of different interfaces are compared utilizing the developed pulsed Current Injection (CI) method. Finally, the impedance measurement scheme is verified through testing using three measuring instruments. And the CI simulation experiments reveal the characteristics and susceptibilities of different PLCs interfaces, indicating that some protection measures are required for the reliable operation of the PLC.

42 ENGINEERING↗
DOE OSTI2020

End-to-End Physics Event Classification with CMS Open Data: Applying Image-Based Deep Learning to Detector Data for the Direct Classification of Collision Events at the LHC

This paper describes the construction of novel end-to-end image-based classifiers that directly leverage low-level simulated detector data to discriminate signal and background processes in pp collision events at the Large Hadron Collider at CERN. To better understand what end-to-end classifiers are capable of learning from the data and to address a number of associated challenges, we distinguish the decay of the standard model Higgs boson into two photons from its leading background sources using high-fidelity simulated CMS Open Data. We demonstrate the ability of end-to-end classifiers to learn from the angular distribution of the photons recorded as electromagnetic showers, their intrinsic shapes, and the energy of their constituent hits, even when the underlying particles are not fully resolved, delivering a clear advantage in such cases over purely kinematics-based classifiers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Modeling of a chain of three plasma accelerator stages with the WarpX electromagnetic PIC code on GPUs

The fully electromagnetic particle-in-cell code WarpX is being developed by a team of the U.S. DOE Exascale Computing Project (with additional non-U.S. collaborators on part of the code) to enable the modeling of chains of tens to hundreds of plasma accelerator stages on exascale supercomputers, for future collider designs. The code is combining the latest algorithmic advances (e.g., Lorentz boosted frame and pseudo-spectral Maxwell solvers) with mesh refinement and runs on the latest computer processing unit and graphical processing unit (GPU) architectures. In this paper, we summarize the strategy that was adopted to port WarpX to GPUs, report on the weak parallel scaling of the pseudo-spectral electromagnetic solver, and then present solutions for decreasing the time spent in data exchanges from guard regions between subdomains. In Sec. IV, we demonstrate the simulations of a chain of three consecutive multi-GeV laser-driven plasma accelerator stages.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Computational prediction of dielectric breakdown strength of a transformer paper in oil with uncertainty quantification

The determination of the dielectric breakdown strengths of microstructurally heterogeneous materials has been a primarily experimental endeavor. We report the development of a microstructure-level model for computationally predicting the breakdown strength and analyzing the interactions between electromagnetic pulses (EMP) and the constituents in a composite of cellulose-based paper and mineral oil found in electrical transformers. The model allows explicit simulation of the material breakdown process by tracking the transition of dielectric constituents from non-conductive to conductive states. The focus is on the electric fields induced in the materials and the overall conditions for dielectric breakdown (defined as the onset of avalanche) caused by the electric field induced in the composite. Responses to three distinct pulse shapes, i.e., Steep Front (SF), Lightning (L), and AC with spectra spanning 60–9 × 105 Hz are considered. It is found that the breakdown strength of the material is significantly affected by microstructure heterogeneities, the spatial variations of the constituent properties, and the pulse shapes. A probabilistic characterization of the breakdown strength is computationally obtained and compared with experimental measurements. Although one particular material is analyzed, the model and approach are applicable to other heterogeneous materials as well.

breakdowns↗
DOE OSTI2020

Reconstruction and Measurement of $\mathcal{O}$(100) MeV Energy Electromagnetic Activity from $\pi^0 \rightarrow \gamma\gamma$ Decays in the MicroBooNE LArTPC

We present results on the reconstruction of electromagnetic (EM) activity from photons produced in charged current νμ interactions with final state π0s. We employ a fully-automated reconstruction chain capable of identifying EM showers of &calO;(100) MeV energy, relying on a combination of traditional reconstruction techniques together with novel machine-learning approaches. These studies demonstrate good energy resolution, and good agreement between data and simulation, relying on the reconstructed invariant π0 mass and other photon distributions for validation. The reconstruction techniques developed are applied to a selection of νμ + Ar → μ + π0 + X candidate events to demonstrate the potential for calorimetric separation of photons from electrons and reconstruction of π0 kinematics.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Quantum Electromagnetic Transients Program

This work devises a quantum electromagnetic transients program (QEMTP) which is the first attempt to tackle the computational challenges in solving EMTP through quantum computing. The main contributions lie in: (1) A quantum-enabled EMTP formulation with Dommel's model encoded in quantum states; 2) An Harrow-Hassidim-Lloyd (HHL)-based QEMTP to solve the discrete-time nodal equations; 3) An iteration-based HHL revision for mitigating temporal errors to achieve high accuracy of QEMTP with limited quantum resources. Case studies verify the correctness and efficacy of QEMTP in simulating both fast and slow dynamics.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗