Search NASA⌕ Search

SEARCH · Search NASA

Results for “EM 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

PHASE: Personalized Head-based Automatic Simulation for Electromagnetic properties in 7T MRI

Accurate and individualized human head models are becoming increasingly important for electromagnetic (EM) simulations. These simulations depend on precise anatomical representations to realistically model electric and magnetic field distributions, particularly when evaluating Specific Absorption Rate (SAR) within safety guidelines. State of the art simulations use the Virtual Population due to limited public resources and the impracticality of manually annotating patient data at scale. Here, this paper introduces Personalized Head-based Automatic Simulation for EM properties (PHASE), an automated open-source toolbox that generates high-resolution, patient-specific head models for EM simulations using paired T1-weighted (T1w) magnetic resonance imaging (MRI) and computed tomography (CT) scans with 14 tissue labels. To evaluate the performance of PHASE models, we conduct semi-automated segmentation and EM simulations on 15 real human patients, serving as the gold standard reference. The PHASE model achieved comparable global SAR and localized SAR averaged over 10 grams of tissue (SAR-10g), demonstrating its potential as a promising tool for generating large-scale human model datasets in the future. The code and models of PHASE toolbox have been made publicly available: https://github.com/hrlblab/PHASE.

Deep learning↗
DOE OSTI2026

Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion (Final Scientific/Technical Report)

This report summarizes the technical accomplishments of the four-year research project “Multiphysics and Multiscale Simulation Methods for Electromagnetic Energy Assisted Fossil Fuel to Hydrogen Conversion” (Award No. DE-FE0032092), conducted at Howard University and the University of Houston (subawardee) from September 2021 to August 2025. The project successfully achieved all four major objectives: 1. 3D Structural Characterization – Developed 3D optical imaging and mechanical sectioning methods to characterize catalyst distribution and support morphology in nickel foam substrates. Successfully reconstructed 3D geometries and imported them into COMSOL Multiphysics for electromagnetic simulations. 2. EM Hotspot Simulation – Created all-frequency stable electromagnetic formulations and 3D nodal discontinuous Galerkin (NDG) methods for coupled electromagnetic-thermal-fluid problems in multiscale catalytic media. Demonstrated stable solutions from DC to microwave frequencies. 3. Multiphysics Coupling – Developed multiscale simulation methods coupling FEM electromagnetic solvers with thermal transport equations. Reactive molecular dynamics (ReaxFF MD) simulations were performed to investigate catalytic reaction mechanisms at the atomistic level. Demonstrated electromagnetic-thermal co-simulation capabilities for porous catalyst structures. 4. System Optimization – Designed and optimized EM-assisted catalytic systems using nickel foam and carbon foam structures, demonstrating significant temperature increases due to microwave heating. Observed and characterized plasma generation in carbon fiber catalysts. Investigated multiple reaction chamber geometries for improved microwave energy deposition. The project produced significant scientific contributions including 15+ peer-reviewed publications, trained multiple Ph.D. students and undergraduate researchers, and advanced the understanding of microwave-assisted hydrogen production from fossil fuels.

08 HYDROGEN↗
DOE OSTI2025

Structural Analysis of a Serpentine Superconducting Magnet for the Interaction Region of the Electron Ion Collider

The electron ion collider (EIC) is under design to collide the high energy and highly polarized hadron with electron beams with luminosities up to 10 34 cm −2 s −1 . The superconducting magnet designs at the interaction region (IR) are quite challenging due to the close proximity of the hadron and electron beams. Several serpentine types of superconducting magnets have been designed due to the space restrictions at the IR. Here, to validate the designs in terms of the mechanical strength, structural integrity, as well as the allowed deformations, a detailed structural mechanical analysis for a 12-layer superconducting quadrupole serpentine magnet (Q1BpR) in the rear side hadron beams has been carried out using the finite element method (FEM) considering the pretension stress, shrinking due to cool down, and the electromagnetic force on the conductors. Anisotropy material properties and individual roving tension have been considered. The electro-magnetic (EM) simulations were performed using COMSOL, confirmed with RAT, the structural mechanical analysis were performed using COMSOL by considering the contact elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1993

Solar wind-magnetosphere interaction as simulated by a 3D, EM particle code

The results of simulating the solar wind-magnetosphere interaction with a three dimensional, electromagnetic (EM) particle code are presented. Hitherto such global simulations were done with magnetohydrodynamic (MHD) codes while lower dimensional particle or hybrid codes served to account for microscopic processes and such transport parameters as have to be introduced ad hoc in MHD. The kinetic model combines macroscopic and microscopic tasks. It relies only on the Maxwell curl equations and the Lorentz equation for particles. The preliminary results are for an unmagnetized solar wind plasma streaming past a dipolar magnetic field. The results show the formation of a bow shock and a magnetotail, the penetration of energetic particles into cusp and radiation belt regions, and dawn to dusk asymmetries.

Buneman, O.↗
NASA NTRS1992

Solar wind-magnetosphere interaction as simulated by a 3-D EM particle code

We present here our first results of simulating the solar wind-magnetosphere interaction with a new three-dimensional electromagnetic particle code. Hitherto such global simulations were done with MHD codes while lower-dimensional particle or hybrid codes served to account for microscopic processes and such transport parameters as have to be introduced ad hoc in MHD. Our kinetic model attempts to combine the macroscopic and microscopic tasks. It relies only on the Maxwell curl equation and the Lorentz equation for particles, which are ideally suited for computers. The preliminary results shown here are for an unmagnetized solar wind plasma streaming past a dipolar magnetic field. The results show the formation of a bow shock and a magnetotail, the penetration of energetic particles into cusp and radiation belt regions, and dawn-dusk asymmetries.

Buneman, Oscar↗
NASA NTRS1994

Loop heating by D.C. electric current and electromagnetic wave emissions simulated by 3-D EM particle zone

We have shown that a current-carrying plasma loop can be heated by magnetic pinch driven by the pressure imbalance between inside and outside the loop, using a 3-dimensional electromagnetic (EM) particle code. Both electrons and ions in the loop can be heated in the direction perpendicular to the ambient magnetic field, therefore the perpendicular temperature can be increased about 10 times compared with the parallel temperature. This temperature anisotropy produced by the magnetic pinch heating can induce a plasma instability, by which high-frequency electromagnetic waves can be excited. The plasma current which is enhanced by the magnetic pinch can also excite a kinetic kink instability, which can heat ions perpendicular to the magnetic field. The heating mechanism of ions as well as the electromagnetic emission could be important for an understanding of the coronal loop heating and the electromagnetic wave emissions from active coronal regions.

Sakai, J. I.↗
NASA NTRS1994

Coalescence of two current loops with a kink instability simulated by a three-dimensional electromagnetic particle code

We have studied the dynamics of a coalescence of current loops using three-dimensional electromagnetic (EM) particle simulation code. Our focus is the investigation of such kinetic processes as energy trasnfer, heating particles, and electromagnetic emissions associated with a current loop coalescence which cannot be studied by MHD simulations. First, the two loops undergo a pinching oscillation due to a pressure imbalance between the inside and outside of the current loop. During the pinching oscillation, a kinetic kink instability is excited and electrons in the loops are heated perpendicularly to an ambient magnetic field. Next, the two current loops collide and coalesce, while at the same time a helical structure grows further. Subsequently, the perturbed current, which is due to these helically bunched electrons, can drive a whistler instability. It should be noted in this case that the whistler wave is excited by the kinetic kink instability and not a beam instability. After the coalescence of two helical loops, tilting motions can be observed in the direction of left-hand rotation, and the helical structure will relax resulting in strong plasma heating mostly in the direction perpendicular to the ambient magnetic field. It is also shown that high-frequency electromagnetic waves can be emitted from the region where the two loops coalesce and propagate strongly in the direction of the electron drift velocity. These processes may be important in understanding heating mechansims for coronal loops as well as radio wave emission mechanisms from active regions of solar plasmas.

Nishikawa, K.-I.↗
NASA NTRS1986

Electron-cyclotron maser instability in relativistic plasmas

The electron-cyclotron maser instability is studied for the case of an anisotropic electron velocity distribution in the regime where the relativistic corrections to the wave dispersion are significant. Solution of the linear dispersion relation reveals that when the plasma frequency-gyrofrequency ratio is less than v(te)/c, the instability is localized just below k(perpendicular)c/Omega(e) = 1. The growth rate is then strongly peaked for emission at 90 deg to the magnetic field and is considerably larger than would be the case if the cold-plasma dispersion theory were valid. These features are confirmed by EM particle simulations.

Pritchett, P. L.↗
DOE OSTI2025

Accelerating detector simulations with Celeritas: profiling and performance optimizations

Celeritas is a GPU-optimized MC particle transport code designed to meet the growing computational demands of next-generation HEP experiments. It provides efficient simulation of EM physics processes in complex geometries with magnetic fields, detector hit scoring, and seamless integration into Geant4-driven applications to offload EM physics to GPUs. Recent efforts have focused on performance optimizations and expanding profiling capabilities. This paper presents some key advancements, including the integration of the Perfetto system profiling tool for detailed performance analysis and the development of track-sorting methods to improve computational efficiency.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
NASA NTRS1993

Interaction with the lower ionosphere of electromagnetic pulses from lightning - Heating, attachment, and ionization

A Boltzmann formulation of the electron distribution function and Maxwell's equations for the EM fields are used to simulate the interaction of lightning radiated EM pulses with the lower ionosphere. Ionization and dissociative attachment induced by the heated electrons cause significant changes in the local electron density, N(e). Due to 'slow' field changes of typical lightning EM pulses over time scales of tens of microsec, the distribution function follows the quasi-equilibrium solution of the Boltzmann equation in the altitude range of interest (70 to 100 km). The EM pulse is simulated as a planar 100 microsec long single period oscillation of a 10 kHz wave injected at 70 km. Under nighttime conditions, individual pulses of intensity 10-20 V/m (normalized to 100 km horizontal distance) produce changes in N(e) of 1-30 percent while a sequence of pulses leads to strong modification of N(e) at altitudes less than 95 km. The N(e) changes produce a 'sharpening' of the lower ionospheric boundary by causing a reduction in electron density at 75-85 km (due to attachment) and a substantial increase at 85-95 km (due to ionization) (e.g., the scale height decreases by a factor of about 2 at about 85 km for a single 20 V/m EM pulse). No substantial N(e) changes occur during daytime.

Taranenko, Y. N.↗
DOE OSTI2025

Evaluating Variable-Impedance Magnetically-Insulated Transmission Lines as a Risk-Mitigation Measure for Next-Generation Pulsed Power

This project has produced the first detailed characterizations of power flow resulting from applying the “variable-impedance MITL” concept to real-life systems in Sandia’s pulsed power program (Z and next-generation pulsed power (NGPP)). We present simulation results and analyses for constant-impedance versions of both Z and NGPP and survey the operational viability of several variable-impedance re-designs in the parameter space of linear tapers. Circuit modeling (SCREAMER/Bertha) was used to pinpoint promising candidate designs, and EM-PIC (Empire) simulations were used to evaluate these candidates more rigorously. This approach was particularly successful in the Z regime which resulted in the identification of several viable variable-impedance MITL designs for each level. The approach was more challenged in the operating space NGPP occupies, producing data points that speak to a more restrictive design space due to anode plasma turn-on. In the end, we were able to converge on one viable variable-impedance design for the highest inductance line (level “F”) and one for the highest current line (level “A”). Altogether, the body of simulation evidence presented in this report suggest there does exist flexibility in operating space for magnetically-insulated transmission lines (MITLs) having variable geometric impedance to be a potential enabling technology for safely increasing current delivery (and potentially lowering stack voltage) in pulsed-power drivers by manipulating electron losses; however, operating points for a particular design must be carefully screened. Circuit and EM-PIC modeling provided consistent verdicts in safe operating regimes for operational viability, but additional physics such as anode plasma turn-on which is included in Empire but not in SCREAMER/Bertha was found to be a critical factor affecting power flow that lead to different assessments between the codes. It is not always the case that the occurrence of anode plasma caused a design to fail (some designs turned on anode plasma yet still delivered load currents meeting design targets); the details matter such as how early in the pulse anode surfaces break down (and how large a region). However, in every case that it did fail it was found that the feedback from anode plasma was the cause (i.e., turning off the anode plasma model in Empire restored agreement with the circuit model prediction). As circuit simulations represent an efficient and practical means of surveying design space compared to more computationally-expensive approaches such as EM-PIC, it could be prudent to invest in the research and development of models to include the effects of anode plasma such as ion emission in circuit codes. The variable-impedance MITL design is a new concept that enables controlled manipulation of the initial electron losses in the outer MITL and can be tested on Z today. We encourage follow-on work to explore further optimization (including alternative variable-impedance profiles, e.g., having constant dZ/dR), and to confirm the major findings presented in this report by fielding test hardware on actual Z shots.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Accelerating detector simulations with Celeritas: Profiling and performance optimizations

Celeritas is a GPU-optimized Monte Carlo (MC) particle transport code designed to meet the growing computational demands of next-generation high energy physics (HEP) experiments. It provides efficient simulation of electromagnetic (EM) physics processes in complex geometries with magnetic fields, detector hit scoring, and seamless integration into Geant4-driven applications to offload EM physics to GPUs. Recent efforts have focused on performance optimizations and expanding profiling capabilities. This paper presents some key advancements, including the integration of the Perfetto system profiling tool for detailed performance analysis and the development of track-sorting methods to improve computational efficiency.

Lund, Amanda [Argonne National Laboratory (ANL)]↗
NASA NTRS1994

Superconducting axisymmetric finite elements based on a gauged potential variational principle. Part 1: Formulation

The present work is part of a research program for the numerical simulation of electromagnetic (EM) fields within conventional Ginzburg-Landau (GL) superconductors. The final goal of this research is to formulate, develop and validate finite element (FE) models that can accurately capture electromagnetic thermal and material phase changes in a superconductor. The formulations presented here are for a time-independent Ginzburg-Landau superconductor and are derived from a potential-based variational principle. We develop an appropriate variational formulation of time-independent supercontivity for the general three-dimensional case and specialize it to the one-dimensional case. Also developed are expressions for the material-dependent parameters alpha and beta of GL theory and their dependence upon the temperature T. The one-dimensional formulation is then discretized for finite element purposes and the first variation of these equations is obtained. The resultant Euler equations contain nonlinear terms in the primary variables. To solve these equations, an incremental-iterative solution method is used. Expressions for the internal force vector, external force vector, loading vector and tangent stiffness matrix are therefore developed for use with the solution procedure.

Schuler, James J.↗
DOE OSTI2026

Component-Level Inverse Design of Transmon Qubits Using Neural Networks

Designing a superconducting qubit to realize specific Hamiltonian parameters typically requires iterating through a time and compute-intensive forward loop in which the designer chooses a layout geometry, simulates it, extracts circuit parameters such as capacitances, and refines the geometry. We study the inverse version of this task using a neural-network workflow that maps target Hamiltonian parameters directly to component-level layout parameters, which we subsequently demonstrate on a planar transmon layout. During training, we pair the inverse model with a frozen forward surrogate model and evaluate the loss in Hamiltonian space rather than in layout-parameter space. In validation against a conventional EM solver, 97% of generated designs produce usable geometries, and the inverse-plus-surrogate pipeline reaches mean percent errors of 0.73% for qubit frequency and 1.58% for anharmonicity, comparable to or below the fabrication and simulation-to-measurement uncertainty expected for academic-process transmon devices of this type. A single pipeline query takes ~60 ms on CPU, versus ~2 min for a conventional EM capacitance extraction on the same hardware, a speedup of approximately 2,000x. Batching minimizes the AI model inference overhead, reducing the runtime to 3.1 microseconds per sample on CPU and 2.6 microseconds per sample on GPU at a batch size of 2048, resulting in speedups of 3.9 x 10^7 and 4.6 x 10^7, respectively, relative to a single conventional CPU EM extraction. Our results indicate that component-level inverse design usefully extends and complements conventional EM simulation, including for small datasets on the order of 1,000 samples.

Seidel, Olivia [Fermilab; Texas U., Arlington]↗
DOE OSTI2025

Magnetosphere Evolution and Precursor-driven Electromagnetic Signals in Merging Binary Neutron Stars

We detail new force-free simulations to investigate magnetosphere evolution and precursor electromagnetic (EM) signals from binary neutron stars. Our simulations fully follow a representative inspiral motion, capturing the intricate magnetospheric dynamics and their impact on EM outflows. We explore a range of stellar magnetic moment orientations and relative strengths, finding that the magnetospheres and Poynting flux evolution are strongly configuration dependent. The Poynting flux exhibits pulsations at twice the orbital frequency, 2Ω, and is highly anisotropic, following a power-law dependence on orbital frequency. The index ranges from 1 to 6, shaped by the intricate dynamics of the magnetospheres. Furthermore, we present the first computation of (1) the EM forces acting on the star surfaces, revealing the presence of torques that, for highly magnetized stars, could influence the orbital dynamics or break the crust; (2) the high-energy emission signals from these systems by adopting the established isolated pulsar theory. Assuming curvature radiation in the radiation-reaction limit, we find that photons could reach TeV–PeV energies in the last ∼ms for magnetic field strengths 10 10 –10 15 G. However, our analysis of single photon magnetic pair production suggests that these photons are unlikely to escape, with the MeV band emerging as a promising observational window for precursor high-energy emission. In this framework, we construct high-energy emission skymaps and light curves, exploring observational implications. Finally, we propose potential precursor radio emission and delayed afterglow echoes from magnetized outflows, which may contribute to late-time rebrightening in short gamma-ray bursts or to orphan afterglows.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Heterogeneous energetic material damage simulator (HEDS): A deep learning approach to simulate damage–sensitivity linkages

Damage in the microstructures of energetic materials (EMs), such as propellants and plastic bonded explosives (PBXs), can significantly alter their response to external loads. Both sensitization and desensitization can occur, causing concerns with safety and performance in the field; predictive models that connect damage and the sensitivity of EMs can enable design and provide confidence in their robustness and reliability. However, modeling of damage evolution is challenging for real microstructures of EMs; samples of damaged EMs are difficult to obtain, thereby hindering experiments and direct numerical simulations to determine the sensitivity of EMs at various stages of damage. Here, we develop an approach to generate synthetic, i.e., in silico produced, damaged microstructures for use in simulations to connect damage levels to sensitivity. The development of the present workflow to generate and impose varying levels of damage in microstructures, known as HEDS (Heterogeneous Energetic Material Damage Simulator), begins with a small set of images of damaged PBXs and combines a collection of deep neural network techniques to generate microstructures with varying levels of damage. By making the synthetic microstructures conform closely to those observed in available real, imaged microstructures, we develop an ensemble of damaged microstructures that can be used for in silico shock experiments. HEDS develops these microstructure ensembles as level set fields, which are directly employed in a sharp interface Eulerian hydrocode where shock simulations are performed to quantify the energy release rate from hotspot fields generated in the microstructure. These capabilities can be useful for the analysis and assessment of changes in the sensitivity of EMs and to design formulations that are less susceptible to damage-induced changes in sensitivity and performance.

Fang, Irene (ORCID:0009000844557122)↗
DOE OSTI2025

stor4build

The EnergyPlus simulation engine supports modeling and simulation of thermal energy storage (TES) systems in several ways, including using the Python-EMS feature, which extends the operation of the engine with custom code written in Python. Creation of models using this feature can be tedious and error prone, with the connection of the model components to the Python code a particularly troublesome area. The stor4build Python package simplifies this process by modifying an input model to add a selected TES technology (implemented with the Python-EMS feature) and runs the simulation. The package leverages the OpenStudio middleware software development kit to automate this process as much as possible, eliminating potential errors and simplifying usage of EnergyPlus. The package provides objects, functions, and OpenStudio measures that implement the necessary operations to automate the creation of EnergyPlus models that integrate TES technologies with building systems. In addition, two user interfaces are provided: a command line interface and a web application programming interface. The automated process implemented by the package greatly simplifies the modeling and simulation process, allowing for parametric studies to be executed much more efficiently and effectively. The OpenStudio-based workflow is also very flexible and will allow for future additions of new technologies.

DeGraw, JasonWilliam [Oak Ridge National Laborator↗