Search NASA⌕ Search

SEARCH · Search NASA

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

Convergence study of wakefield simulations with GdfidL and ECHO3D

The interaction of charged particle beams with vacuum chamber components gives rise to electromagnetic wakefields, whose frequency-domain representation is known as beam coupling impedance. Geometric impedance arising from discontinuities and transitions in the vacuum chamber is the focus of this study. Minimizing this impedance is essential to mitigate adverse collective effects in modern storage rings operating with high-intensity particle beams. Accurate and reliable impedance simulations is a key factor of the vacuum chamber design. This paper presents the results of a convergence study of two widely used electromagnetic solvers, GdfidL and ECHO3D, applied to key vacuum-chamber components of the National Synchrotron Light Source II (NSLS-II) storage ring. Detailed comparisons are performed for several geometries, including flange absorbers, RF bellows, button-type beam position monitors, and an in-vacuum undulator (IVU). The results show notable differences in convergence and computational efficiency between the two codes. While GdfidL provides highly resolved results and serves as a common benchmark tool, ECHO3D yields consistent results with coarser meshes, significantly reducing simulation time and memory demands. Simulations with a full-geometry IVU model demonstrate that simplified taper-transition models can miss important impedance contributions. In conclusion, these findings provide practical guidelines for efficient and accurate impedance modeling to optimize design of vacuum chamber components for accelerators.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Black-Starting Microgrids Using Inverter-Based Distributed Energy Resources

Recent trends in power systems show increasing deployment of inverter-based distributed energy resources (DERs) at the distribution level. Efforts have also been made to use inverter-based DERs for system restoration of critical sections of distribution feeders. In the past few years, there has been active consideration of using local grid-forming (GFM) DERs, which can form their own voltage and frequency, to black-start microgrids instead of waiting for the bulk system for supply restoration. However, using current-limited inverter-based DERs for black-start operations with high inrush currents can pose some technical challenges. This paper presents detailed electromagnetic transient studies on a digital real-time simulator to assess the operational feasibility of using local inverter-based GFM DERs for microgrid black-start operations as a bottom-up approach instead of a top-down approach where the substation is used to restore critical sections of the feeder. The paper presents key modeling details for distribution system components for studying GFM DER-based microgrid black-start operations. The paper presents simulation results comparing different black-start approaches using GFM DERs for two test systems: i) a simple test microgrid and ii) a real-world utility microgrid currently being designed and deployed by Holy Cross Energy in Colorado.

black start↗
DOE OSTI2025

Contrasting Time-Frequency Representations for Unknown Waveform Detection

In real-world applications like spectrum management and interference detection, dealing with unseen electromagnetic waveforms is critical. Although some methods attempt to simulate open set data using generator models, they face challenges in generating synthetic samples for open set while simultaneously selecting an optimal discriminator for accurate classification. This results in difficulties capturing distinctive features across classes, especially in dynamic scenarios where new classes emerge. To detect unseen waveforms, we propose combining time and frequency domain features with cosine similarity loss to enhance feature distinctiveness and enabling more accurate predictions. This approach efficiently captures more comprehensive information than single-domain representations or approaches without cosine loss. Additionally, our model avoids generic feature vectors by extracting class-specific features during training, resulting in improved class representation. The experiment results show that this combined feature approach with cosine loss outperforms single-domain models and improves accuracy by 10\% over models without cosine loss.

99 - GENERAL AND MISCELLANEOUS↗
DOE OSTI2025

Effect of energetic ions on edge-localized modes in tokamak plasmas

The most efficient and promising operational regime for the International Thermonuclear Experimental Reactor tokamak is the high-confinement mode. In this regime, however, periodic relaxations of the plasma edge can occur. These edge-localized modes pose a threat to the integrity of the fusion device. Here we reveal the strong impact of energetic ions on the spatio-temporal structure of edge-localized modes in tokamaks using nonlinear hybrid kinetic–magnetohydrodynamic simulations. A resonant interaction between the fast ions at the plasma edge and the electromagnetic perturbations from the edge-localized mode leads to an energy and momentum exchange. Energetic ions modify, for example, the amplitude, frequency spectrum and crash timing of edge-localized modes. The simulations reproduce some observations that feature abrupt and large edge-localized mode crashes. The results indicate that, in the International Thermonuclear Experimental Reactor, a strong interaction between the fusion-born alpha particles and ions from neutral beam injection, a main heating and fast particle source, is expected with predicted edge-localized mode perturbations. This work advances the understanding of the physics underlying edge-localized mode crashes in the presence of energetic particles and highlights the importance of including energetic ion kinetic effects in the optimization of edge-localized mode control techniques and regimes that are free of such modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Impact of mechanical tolerances on partial turn skipping in helical flux compression generator

Maintaining tight mechanical tolerances of the components used in helical flux compression generators is crucial for optimizing device performance. Moving beyond current literature, which provides various approximations for acceptable tolerances, a 3D simulation of the armature stator interaction during generator operation was developed, revealing that existing equations frequently misestimate the tolerance limits. Without restricting the generality of the presented approach, the focus was primarily on the armature’s concentricity and roundness, comparing the acceptable tolerance limits between generators of different sizes. An analytical approach for select geometries was developed, confirming the 3D results. Quantitative results reveal that a steeper Gurney angle and increased winding pitch allow for less restrictive tolerances for eccentricity and ellipticity. Even small deviations from the ideal armature positioning and shape will result in fluctuations of the output current time-derivative, observed experimentally and in simulation. Larger deviations then cause partial turn skipping, associated with loss of magnetic flux. The simulation elucidates the impact of tolerance deviations on generator performance by offering precise tolerance ranges for various generator sizes, thereby facilitating informed design decisions.

42 ENGINEERING↗
DOE OSTI2025

Three and Two Phase Rotating Field Inductive Couplers for Wireless Power Transfer with One Phase per Layer Windings

Multiphase inductive wireless charging coils have been proposed recently to improve coupler surface power density, reduce component stress and size, and provide near-constant power delivery to charge mobile electric systems. Several aspects for the fundamental characterization of multiphase coils are explored up to six phases including approximate mutual inductance with size and turn variation, induced voltage, and output power estimation. The relative component stress and size of passive components for resonant operation are compared between the multiphase variants. A combination of an experimentally validated 3D electromagnetic finite element analysis (FEA) and power electronic co-simulations are used to validate the estimated quantities approximated with a mixture of analytical equations and an artificial neural network model for mutual inductance. A novel three-phase transmitter, two-phase receiver coil pair is also proposed for electric vehicle charging to reduce the number of connections and compensation complexity on the vehicle-side with improved power output compared to a two-phase configuration.

Lewis, Donovin D. [University of Kentucky]↗
DOE OSTI2026

Design of a high power 8 MeV electron beam injector for the LERF

We present the beamline design for a high power 8 MeV electron injector to be installed at the Low Energy Recirculator Facility (LERF). The injector is designed to provide beam power of up to 80 kW for testing positron source target candidates. We discuss the results of beam dynamics simulations conducted using the General Particle Tracer (GPT) software, where electromagnetic beamline elements were optimized using multi-objective genetic algorithm tools. Additionally, we present preliminary studies on secondary particle generation from both rotating solid and free surface liquid metal jet targets described elsewhere in this workshop.

Lizarraga-Rubio, Victor [Thomas Jefferson National↗
DOE OSTI2026

Design of a high power 8 MeV electron beam injector for the LERF

Here, we present the beamline design for a high power 8 MeV electron injector to be installed at the Low Energy Recirculator Facility (LERF). The injector is designed to provide beam power of up to 80 kW for testing positron source target candidates. We discuss the results of beam dynamics simulations conducted using the General Particle Tracer (GPT) software, where electromagnetic beamline elements were optimized using multi-objective genetic algorithm tools. Additionally, we present preliminary studies on secondary particle generation from both rotating solid and free surface liquid metal jet targets described elsewhere in this workshop.

Lizarraga-Rubio, Victor Manuel [Thomas Jefferson N↗
DOE OSTI2026

Purpose, Design, and Analysis of the NSTX-U Centerstack Casing Lateral Support Shims

The National Spherical Torus Experiment Upgrade (NSTX-U) centerstack casing surrounds the central magnet core of NSTXU. It provides the inner vacuum boundary of the plasma chamber and supports the radially inboard plasma facing components lining the plasma chamber. The casing is supported like a cantilever, fixed at the lower flange by a connection to the TF Bundle through a G-10 ring. It is flexibly supported at the top by bellows, and insulated shims intended to react the largest part of the inventory of the disruption loads at the upper support, thereby protecting the upper and lower connections from excessive loads. A large part of the qualification of the casing and tiles requires addressing and quantifying disruption loads. These have been estimated based on extensions of the behavior of NSTX during its run period roughly 15 years ago. Structural qualification of the Upgrade used these estimates of the loading as input to calculations to demonstrate adequacy for a “most probable” disruption. The estimates are included in a requirements document to guide conservative design of all components loaded by disruptions. However, the changes made for the upgrade made the extrapolation of disruption behavior uncertain. Disruptions include axisymmetric loading from axisymmetric currents induced in the vessel shells. In addition, more complex currents called Halo currents are a result of plasma poloidal currents impinging on the casing. The halo currents add another layer of uncertain behavior and net lateral loading to the casing. As a result, NSTX-U is equipped with a number of instruments intended to assess the nature and severity of disruptions. The casing shims are required for lateral support of the casing, but these have also been converted to load cells to provide another measure of the magnitude of symmetric and asymmetric loading. Some of the logic in the development of the requirements document are presented. Electromagnetic analysis of axisymmetric and halo loading is described. Global simulations of the casing transient dynamic response, and local analysis and testing of the shim/load cell are included. The shim load cells and other instrumentation will be an important guide in planning NSTX-U experiments beyond its next run period expected in late 2026 or early 2027.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Dielectric Resonator Design for Low Power and Low Temperature Microwave Plasma

Waveguide-based microwave plasmas generally operate at high temperatures (2000 - 6000K)[1], making it difficult to directly interface solid materials with the plasma without significant thermal damage. Dielectric microwave resonators (DMRs), long studied for wave-based manipulation of electromagnetic radiation for telecom and optics, can focus radiation to extremely small mode volumes, creating intense localized fields with low-power input.[2] This phenomenon can be used for applications ranging from efficient plasma electronics to near-ambient plasma-materials interactions. Such DMR-based plasmas have been demonstrated a handful of times in the literature, but the majority of research towards this utilize the lowest frequency resonance mode.[3], [4], [5] By carefully controlling the geometry of cylindrical resonators, a variety of electromagnetic modes can be excited. In this work, COMSOL Multiphysics simulations are used to study the electric field enhancement and absorption properties of CaTiO3 DMRs as a function of geometry and excitation frequency. Whereas previous studies have utilized the HEM111 resonance frequency to drive low power plasma excitation, we find that higher order resonance frequencies are more effective at field enhancement and result in less power loss within the dielectric material, hence less wasted heating. The effectiveness of these modes is also geometry dependent and can be computationally optimized for plasma generation. Complementing these computational efforts, we demonstrate a new closed-system reactor design built in a WR-650 waveguide and experimentally demonstrate the formation of atmospheric argon microwave plasma using < 30 W input power on DMR dimers. We observe a shifting resonance frequency as the DMRs heat in response to microwave excitation and develop a Python-based lock-in mechanism to effectively track the DMR resonance over time, leading to stable plasma operation. We use infrared thermal imaging to monitor the temperature of the DMR dimers and surrounding quartz chamber, demonstrating thermal temperatures < 60 degreesC. Finally, we utilize optical emission spectroscopy (OES) to probe the plasma properties as a function of the resonance mode.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Collaborative Research: Vlasov-Maxwell Simulations to Resolve Electron Heating and Dissipation, in Quasi-Perpendicular Shocks (Final Report)

Collisionless shocks are a grand challenge problem in plasma physics and have been the subject of study for more than six decades. Shocks are ubiquitous phenomena in the universe and are responsible for transforming high energy and momentum flows into thermal energy and energetic particles. Understanding how shocks operate is of primary importance to understand the sun-earth coupling, protecting manned missions and spacecrafts from high energy particles, inertial confinement fusion, and radiation observed from astrophysical plasmas, such as supernova remnants and astrophysical jets. A major unanswered question on this frontier is, how does a collisionless plasma transform flow energy and momentum into electron thermal energy? Many potential mechanisms have been proposed to perform the conversion between flow and thermal energy in a collisionless plasma, but the answer has been elusive. To solve the mystery, non-relativistic perpendicular and quasi-perpendicular shocks are studied via a coordinated approach employing cutting-edge numerical simulations with the continuum kinetic code Gkeyll and spacecraft observations supplied by NASA's Magnetospheric Multiscale (MMS) mission. In both sets of data, we employ novel techniques to identify how and where energy is exchanged between electromagnetic fields and plasma particles. By combining cutting-edge spacecraft data and numerical simulations with these diagnostics, we have been able to identify some of the primary mechanisms responsible electron heating and dissipation in heliospheric shocks. In addition to electron dynamics, we have complete access to proton simulation and in situ data, enabling us to characterize proton heating and acceleration in quasi-perpendicular shocks. Many astrophysical quasi-perpendicular shocks have similar parameters to those in this study, including supernova remnants and shocks in galaxy clusters. Additionally, the simulation code, Gkeyll, and associated analysis scripts are openly available for community use. This award has supported a collaborative effort between Princeton University, the University of Arizona - Tucson, and the University of Maryland - College Park. The award funded three early career scientists, including the funding and training of one female postdoctoral researcher, who was involved in all aspects of the MMS data analysis and some portions of the simulation analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical modeling of electromagnetic field spatiotemporal evolution to evaluate the effects on calcium carbonate crystallization

Calcium carbonate (CaCO 3 ) scaling is a significant impediment to water systems. Electromagnetic field (EMF) treatment is a promising approach to control scaling owing to its simplicity and low or no energy requirements. However, the underlying mechanisms by which EMF impacts CaCO 3 crystallization remain unclear due to the challenges in measuring the EMFs in feed solutions and the lack of a fundamental understanding of the applied EMFs and the observed physicochemical phenomena. To fill this knowledge gap, a high-fidelity COMSOL model was first developed to simulate EMFs in bulk solutions for three alternating current-induced EMF devices with different configurations and properties. These were then integrated with experimental data to unveil the underlying mechanism by which applied EMFs alter the physicochemical processes. The study revealed that even low-strength EMFs (e.g., electric fields <0.15 V/m and magnetic fields <0.03 mT) promoted CaCO 3 precipitation in bulk solutions. The electric fields created by these EMF devices resulted in higher Lorentz force compared to their induced magnetic fields. In conclusion, the methodology of this study offers the capability to predict the effectiveness of different EMF devices in facilitating crystallization processes, and these mechanistic insights lay the foundation for the smart design of EMF devices for diverse water treatment applications.

COMSOL simulation↗
DOE OSTI2025

Particle‐In‐Cell Simulations of Starfish Prime

The Starfish Prime high altitude nuclear test created a transient diamagnetic cavity in the Earth's magnetic field above Johnston Island and launched an electromagnetic pulse (EMP) that was detected around the globe. For this work, we use the ion‐kinetic particle‐in‐cell code Topanga to simulate diamagnetic cavity evolution and the E3 EMP signal for Starfish Prime out to over a minute of physical time. The simulation domain has a longitudinal and latitudinal extent of 60° X 60° and a vertical extent of 2,000 km from the surface of the Earth. We compare our simulated results to magnetometer measurements taken in space and on the ground, finding good agreement in both cases. The diamagnetic cavity in the simulation forms in about a second, while the associated debris flux tube takes approximately 30 s to decay. The debris flux tube undergoes significant motion during this timeframe, rising upward. The measured E3 EMP signal on the ground consists of several components, all of which are present in our simulations. We discuss the physical origin of these components in relation to E3a (blast) and E3b (heave).

Belyaev, Mikhail A. [Lawrence Livermore National L↗
DOE OSTI2025

Origin of Pulsed Radio Emission from Magnetars

Extended periods of radio pulsations have been observed for six magnetars, displaying characteristics different from those of ordinary pulsars. In this Letter, we argue that radio emission is generated in a closed, twisted magnetic flux bundle originating near the magnetic pole and extending beyond 100 km from the magnetar. The electron–positron flow in the twisted bundle has to carry electric current and, at the same time, experiences a strong drag from the radiation field of the magnetar. This combination forces the plasma into a “radiatively locked” state with a sustained two-stream instability, generating radio emission. We demonstrate this mechanism using novel first-principles simulations that follow the plasma behavior by solving the relativistic Vlasov equation with the discontinuous Galerkin method. First, using one-dimensional simulations, we demonstrate how radiative drag induces the two-stream instability, sustaining turbulent electric fields. When extended to two dimensions, the system produces electromagnetic waves, including superluminal modes capable of escaping the magnetosphere. We measure their frequency and emitted power and incorporate the local simulation results into a global magnetospheric model. The model explains key features of the observed radio emission from magnetars: its appearance after an X-ray outburst, wide pulse profiles, luminosities ∼10 30 erg s −1 , and a broad range of frequencies extending up to ∼100 GHz.

Astronomical simulations↗
DOE OSTI2025

Extracting symplectic maps for space-charge dominated beams

Symplecticity of transfer maps is important for reliable evaluation of space-charge dominated beams in accelerators. Unfortunately, most simulation codes that include collective effects, such as space charge, do not use canonical phase-space variables and therefore are not symplectic in the presence of electromagnetic fields. In this paper, we present a numerical method to extract local linear symplectic transfer maps using particle tracking simulation code for space-charge dominated beams. We demonstrate this method for the photoinjector (113 MHz SRF gun) section of the Coherent electron Cooling (CeC) Proof of Principle (POP) experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Advanced fuel fusion, phase space engineering, and structure-preserving geometric algorithms

Non-thermal advanced fuel fusion trades the requirement of a large amount of recirculating tritium in the system for that of large recirculating power. Phase space engineering technologies utilizing externally injected electromagnetic fields can be applied to meet the challenge of maintaining non-thermal particle distributions at a reasonable cost. The physical processes of the phase space engineering are studied from a theoretical and algorithmic perspective. It is emphasized that the operational space of phase space engineering is limited by the underpinning symplectic dynamics of charged particles. The phase space incompressibility according to the Liouville theorem is just one of many constraints, and Gromov's non-squeezing theorem determines the minimum footprint of the charged particles on every conjugate phase space plane. In this sense and level of sophistication, the mathematical abstraction of phase space engineering is symplectic topology. To simulate the processes of phase space engineering, such as the Maxwell demon and electromagnetic energy extraction, and to accurately calculate the minimum footprints of charged particles, recently developed structure-preserving geometric algorithms can be used. The family of algorithms conserves exactly, on discretized spacetime, symplecticity and thus incompressibility, non-squeezability, and symplectic capacities. The algorithms apply to the dynamics of charged particles under the influence of external electromagnetic fields as well as the charged particle–electromagnetic field system governed by the Vlasov–Maxwell equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

COSMIC DAWN: Distributed Analysis of Wireless at Nextscale

Distributed Analysis of Wireless at Nextscale (DAWN) is a novel simulation framework for large-scale design-space exploration (DSE) of unmodified software-defined radio (SDR) applications interacting in a scalable, high-fidelity, virtual physics environment. The software-defined nature of the coupled software-physics simulation leverages hardware emulation to permit in-depth examination and modification of not only the electromagnetic environment, including each signal in flight, but also the precise state of system software and components. DAWN supports modular, customizable physics environments allowing realistic propagation effects so that computationally efficient empirical models, reduced order/surrogate models, or large-scale, high-fidelity, site-specific simulations can be used as a propagation medium based on scenario requirements. This paper introduces DAWN’s design and initial implementation, detailing key architectural components, including the Physics Realization Engine (PhyRE), Runtime Infrastructure for Simulation Environments (RISE), and the design space exploration (DSE) suite. It concludes with demonstrations using unmodified 4G/LTE software available from srsRAN on computing resources ranging from a small cluster to ORNL’s Frontier Exascale system.

Wise, Mike [ORNL] (ORCID:0000000266120641)↗
DOE OSTI2024

Investigation of two-dimensional radio-frequency sheath properties using a microscale fluid model

Abstract In previous work (Kohno and Myra 2023 Comput. Phys. Commun. 291 108841), we developed a numerical scheme based on a two-dimensional microscale radio-frequency (RF) sheath model with periodically curved wall boundaries. Here, we expand the capability of this scheme through modification of the boundary conditions (BCs) on the conducting walls, which allows the ion flow to turn back to the plasma at locations on the walls where the electromagnetic force on the ions is reversed from its usual direction. Numerical simulations are carried out to investigate the dependences of the surface-integrated admittances on the wall bump height, ion magnetization, ion mobility, and the magnetic field angle, and to visualize the sheath structures in several cases. One of the main results is the ion cyclotron admittance resonance observed under the condition of low ion mobility (high normalized frequency). It is shown that the amplitude of the resonance peak depends on the wall bump height and the ion velocity is reversed on the sides of the bump in an RF cycle for the resonance cases. Furthermore, the differences in the admittances between the one- and two-dimensional microscale models are assessed for the purpose of understanding non-locality of the sheath near the wall surface for the parameters considered in this study. This information will be essential for improving the sheath BC for macroscale calculations in the future.

Kohno, H. (ORCID:0000000332690010)↗