Search NASA⌕ Search

SEARCH · Search NASA

Results for “magnetic modeling and 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

Simulation Models for Exploring Magnetic Reconnection

Simulations have played a critical role in the advancement of our knowledge of magnetic reconnection. However, due to the inherently multiscale nature of reconnection, it is impossible to simulate all physics at all scales. For this reason, a wide range of simulation methods have been crafted to study particular aspects and consequences of magnetic reconnection. This article reviews many of these methods, laying out critical assumptions, numerical techniques, and giving examples of scientific results. Plasma models described include magnetohydrodynamics (MHD), Hall MHD, Hybrid, kinetic particle-in-cell (PIC), kinetic Vlasov, Fluid models with embedded PIC, Fluid models with direct feedback from energetic populations, and the Rice Convection Model (RCM).

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Machine-learning force-field models for dynamical simulations of metallic magnets

We review recent advances in machine-learning (ML) force-field methods for Landau–Lifshitz–Gilbert simulations of itinerant electron magnets, focusing on their scalability and transferability. Built on the principle of locality, a deep neural-network model is developed to efficiently and accurately predict electron-mediated forces governing spin dynamics. Symmetry-aware descriptors constructed through a group-theoretical approach ensure rigorous incorporation of both lattice and spin-rotation symmetries. The framework is demonstrated using the prototypical s-d exchange model widely employed in spintronics. ML-enabled large-scale simulations reveal novel nonequilibrium phenomena, including anomalous coarsening of tetrahedral spin order on the triangular lattice and the freezing of phase-separation dynamics in lightly hole-doped, strong-coupling square-lattice systems. These results establish ML force-field frameworks as scalable, accurate, and versatile tools for modeling nonequilibrium spin dynamics in itinerant magnets.

Artificial neural networks↗
DOE OSTI2024

First-principles study of the magneto-Raman effect in van der Waals layered magnets

Magneto-Raman spectroscopy has been used to study spin-phonon coupling in two-dimensional (2D) magnets. Raman spectra of CrI 3 show a strong dependence on the magnetic order within a layer and between the layers. Here we carry out the first systematic theoretical investigation of the magneto-Raman effect in 2D magnets by performing density functional theory calculations and developing a generalized polarizability model. Our first-principles simulations well reproduce experimental Raman spectra of CrI 3 with different magnetic states. The model reveals how the change of spin orientation in each layer is coupled to the layer’s vibration to induce or eliminate the spin-dependent anti-symmetric off-diagonal terms in the Raman tensor for altering the selection rules. We also uncover that the correlation between phonon modes and magnetic orders is a universal phenomenon, which should exist in other phonon modes and 2D magnets. Our predictive simulations and modeling are expected to guide the research in 2D magnets.

36 MATERIALS SCIENCE↗
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 OSTI2025

Machine learning force field model for kinetic Monte Carlo simulations of itinerant Ising magnets

Here, we present a scalable machine learning (ML) framework for large-scale kinetic Monte Carlo (kMC) simulations of itinerant electron Ising systems. As the effective interactions between Ising spins in such itinerant magnets are mediated by conducting electrons, the calculation of energy change due to a local spin update requires solving an electronic structure problem. Such repeated electronic structure calculations could be overwhelmingly prohibitive for large systems. Assuming the locality principle, a convolutional neural network (CNN) model is developed to directly predict the effective local field and the corresponding energy change associated with a given spin update based on Ising configuration in a finite neighborhood. As the kernel size of the CNN is fixed at a constant, the model can be directly scalable to kMC simulations of large lattices. Our approach is reminiscent of the ML force field models widely used in first-principles molecular dynamics simulations. Applying our ML framework to a square-lattice double-exchange Ising model, we uncover unusual coarsening of ferromagnetic domains at low temperatures. Our work highlights the potential of ML methods for large-scale modeling of similar itinerant systems with discrete dynamical variables.

machine learning↗
DOE OSTI2025

Cosmological Neutron Stars Produce Diffuse Axion X-Ray Signatures

Axion-like particles can be abundantly produced through scattering processes in the cores of neutron stars (NSs). If they are ultralight ($m_a \lesssim 10^{-4}$ eV), then they can efficiently convert to detectable photons in the external NS magnetospheres, and if they are heavy ($m_a \gtrsim 1$ eV), then they can decay into photons before reaching Earth. In this work, we search for the resulting X-ray signatures from both of these channels summing over the $\textit{cosmological}$ NS population. We compare the predicted axion-induced X-ray signal to the cosmic X-ray background today as measured by a number of instruments such as NuSTAR, HEAO, Swift, and INTEGRAL. We model the axion-induced signal using NS cooling simulations and magnetic field evolution models. We find no evidence for axions and derive strong constraints for both ultralight and heavy axion scenarios, covering new parameter space for the axion-photon and axion-nucleon couplings. Our results rule out the axion-explanation of the Magnificent Seven X-ray excess from nearby isolated NSs.

Cosmology and Nongalactic Astrophysics (astro-ph.C↗
DOE OSTI2024

Transport coefficient sensitivities in a semi-analytic model for magnetized liner inertial fusion

Performance of magnetized liner inertial fusion (MagLIF) experiments is highly dependent on transport processes including magnetized heat flows and magnetic flux losses. Magnetohydrodynamic simulations used to model these experiments require a choice of model for the transport coefficients, which are the constants of proportionality relating driving terms, such as temperature gradients and currents, to the associated heat and magnetic field transport. The coefficients have been the subject of repeated recalculation using various methods throughout the years. Using a semi-analytic MagLIF model, we compare models for the transport coefficients. The choice of model modifies magnetic-flux losses caused by the Nernst thermoelectric effect and thermal conduction losses. We present simulated results from parameter scans conducted in order to compare the effects of the different models on parameters of interest in MagLIF. In some regions of parameter space, discrepancies of up to 38% are found in integrated quantities like the fusion yield. These results may serve as a guide for experimental validation of the various models, particularly as laser preheat energies and initial axial field strengths are increased on MagLIF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

The Verification and Validation of a Magnetic Plasma Fluid Model Utilizing the MOOSE (Multiphysics Object Oriented Simulation Environment) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification, benchmarking, and validation studies have been conducted. Verification studies involved utilizing the method of manufactured solutions and comparing the convergence slope of a known solution to the theoretical slope. Benchmarking consists of comparisons to existing edge codes, namely BOUT++ and SOLEDGE3X. Validation efforts focused on comparisons against open-source data from the TCV tokamak.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Detecting shearless phase-space transport barriers in global gyrokinetic turbulence simulations with test particle map models

In magnetically confined fusion plasmas, the role played by zonal 𝐸 x 𝐵 flow shear layers in the suppression of turbulent transport is relatively well understood. However, less is understood about the role played by the weak shear regions that arise in the non-monotonic radial electric field profiles often associated with these shear layers. In electrostatic simulations from the global total-𝑓 gyrokinetic particle-in-cell code XGC, we demonstrate how shearless regions with non-zero flow curvature form zonal ‘jets’ that, in conjunction with neighbouring regions of shear, can act as robust barriers to particle transport and turbulence spreading. By isolating quasi-coherent fluctuations radially localised to the zonal jets, we construct a map model for the Lagrangian dynamics of gyrokinetic test particles in the presence of drift waves. We identify the presence of shearless invariant tori in this model and verify that these tori act as partial phase-space transport barriers in the simulations. We also demonstrate how avalanches impinging on these shearless tori cause eddy detachment events that form ‘cold/warm core ring’ structures analogous to those found in oceanic jets, facilitating transport across the barriers without destroying them completely. We discuss how shearless tori may generically arise from tertiary instabilities or other types of discrete eigenmodes, suggesting their potential relevance to broader classes of turbulent fluctuations.

fusion plasma↗
DOE OSTI2025

MLSPICE: Machine Learning based SPICE Modeling Platform for Power Magnetics

Electrical power converters are critical to a wide range of applications ranging from renewable integration to transportation electrification, and can be a key factor determining the size, weight, and efficiency of energy conversion systems. Magnetic components are typically the largest and least efficient components in power electronics. While there have been major strides in the modeling and analysis of power semiconductor devices and circuit simulations, the necessary advances in the design of power magnetics have lagged. In this project, we have transformed the modeling and design of power magnetics with machine learning enabled methods and catalyze simultaneous disruptive improvements for ML-based power electronics design tools. A fully automated open-source machine learning based magnetics modeling platform – the MagNet project - with innovations in full stack have been developed to greatly accelerate the design process and provide new insights to magnetic material and geometry design. The ARPA-E funded MagNet platform contains three major building blocks: 1) a ML-Integrated Data Acquisition System (MIDAS): a highly automated data acquisition testbed which is capable of measuring a large number of magnetic cores with a wide range of electrical circuit excitations; 2) a ML-integrated Core Loss Model (MICLM): a machine-learning trained modeling method for modeling the core loss and saturation effects of magnetic materials for arbitrary excitation waveforms; 3) ML-guided Magnetics SPICE Simulation Tool (PMSPICE): a fully integrated CAD tool which can simulate the magnetics in SPICE. It can help the designers to quickly model the linear and non-linear characteristics of magnetic components and evaluate their behavior in SPICE simulations. The developed MagNet system has fully demonstrated the proposed performance target and has been open sourced to the entire power electronics community to advance the modeling and design of power magnetics from many different angles.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Effect of self-generated magnetic fields on x-ray emission in Kr-filled targets at the National Ignition Facility

We examine the effects of self-generated magnetic fields in a Kr gas pipe x-ray source platform. X-ray emission from Kr plasma is dependent on the plasma conditions, as the ionization state is largely a function of temperature. Magnetic fields are known to limit heat conduction, which increases temperature. We show that the emission in simulations of the gas pipe x-ray source is dependent on how self-generated magnetic fields are modeled. The inclusion of self-generated magnetic fields in simulations more accurately captures the emission of lower energy x-ray emission (L-shell), bringing results closer to experiments. The modeled x-ray emission and self-generated magnetic fields are shown to be particularly sensitive to the inclusion of the Nernst effect in simulations. Severely limiting the Nernst effect leads to a hotter Kr plasma, which can account for the discrepancy seen in earlier studies. By modifying the Nernst effect multiplier, we can achieve better experimental agreement in x-ray emission from gas pipes; the value of the multiplier that leads to the best agreement is dependent on the laser power of the drive. Currently, the suppression factor of the Nernst effect needed for high power drives (PL>200 TW) is more restrictive than what is currently put forward by non-local models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

An Integrated Assessment of a G3 GMD Event on Large-Scale Power Grids: From Magnetometer Data to Geomagnetically Induced Current Analysis

Solar activities can cause geomagnetic disturbances (GMDs) that give rise to geomagnetically induced currents (GICs) which may compromise the reliability of the power system. In order to build more reliable models representing GMD interactions with the power grid, the power system’s detailed electrical model must be considered along with fluctuations in the earth’s magnetic and induced surface electric fields. Here, this study investigates the impact of incorporating spatially varying magnetic fields into surface electric field models on GMD risk metrics. A spatially independent magnetic field model and a spatially varying model are compared through simulations. To perform this analysis, the earth’s magnetic field disturbances are transformed into surface electric fields using respective one-dimensional earth conductivity models. Then, the modeling impact of these electric fields is studied using a 2,000-bus grid for Texas and a 25,000-bus grid for the northeast and mid- Atlantic regions of the United States. Simulation results reveal that the inclusion of spatially varying magnetic fields results in considerable differences in GMD risk metrics, highlighting the importance of accounting for spatial variability when assessing GMD risks in the power system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Structure of Self-Generated Magnetic Fields in Laser-Solid Interaction from Proton Tomography

Self-generated magnetic fields in laser-solid interactions are experimentally characterized to reveal the 3D location and local field strength, rather than path-integrated quantities, using multi-view proton radiography and tomographic inversion. We infer magnetic fields that extend several millimeters off the target into the hot, rarefied corona, sufficient to strongly magnetize the plasma (Ω e τ e ≫ 1). The data are compared to MHD simulations incorporating recent improvements in modeling magnetic field generation and transport; the volume-averaged coronal field strength and magnetic flux agree to within 25% using a model with magnetic re-localization of transport, although the near-target morphology is not reproduced. This work demonstrates tomographic proton radiography as a valuable tool for investigating magnetic fields in laser-produced plasmas.

High-energy-density plasmas↗
DOE OSTI2025

Programmable exploration of magnetic states in Lieb-kagome interpolated lattices

We investigate a hybrid modeling framework in which a quantum annealer is used to simulate magnetic interactions in molecular qubit lattices inspired by experimentally realizable systems. Using phthalocyanine assemblies as a structurally constrained prototype, we model a continuous deformation from a Lieb to a kagome lattice, revealing frustration-driven disorder and magnetic field-induced reordering in the spin structure. Here, the goal is to show how a quantum annealer can operate as a physically instantiated, programmable platform to emulate experimentally relevant lattice deformations and produce observables in a manner analogous to an experimental measurement, enabling the characterization of magnetic arrangements beyond the reach of current molecular architectures. This surrogate modeling approach offers a pathway to explore and iteratively design tunable magnetic states in synthetic materials. The synthetic design, structural characterization, and quantum simulation framework established here defines a modular and scalable paradigm for probing the limits of engineered matter across chemistry, condensed matter, and quantum information science.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Analysis of the HAE activity in the TJ-II stellarator using a Landau closure model

The aim of this study is to analyze the stability of helical Alfvén eigenmodes (HAEs) in TJ-II discharges and the stabilizing effect of the energetic particles generated by the neutral beam injector (NBI) on pressure gradient-driven modes (PGDMs). HAE and PGDM stability is studied using the linear version of the gyro-fluid code FAR3d and the continuous structure by the STELLGAP code. First, Alfvén eigenmode (AE) and PGDM activity observed in the experiments is reproduced by the simulations, identifying unstable m/n = 4/7 − 2/3 and 7/12 − 5/8 HAEs triggered around ρ = 0.66 showing a frequency of 209 and 204 kHz, respectively, as well as 5/3 PGDM. Next, a parametric study is performed with respect to the thermal ion density and iota profile in the middle-outer plasma region to verify the robustness of the simulation results with respect to the uncertainty of experimental profiles. The analysis confirms that experimental uncertainty does not cause large deviations in the simulation results, showing the destabilization of the same HAEs for all the configurations tested. The simulations also indicate the decay of the 5/3 PGDM growth rate as the energetic particle (EP) population in the plasma increases, consistent with the experiment. Stability analysis of the n = 3, 7, 11, n = 5, 9, 13, n = 6, 10, 14, and n = 8, 12 helical families is performed with respect to the NBI operational regime for different EP energies, β as well as deposition profiles. The most unstable configuration is the radially localized on-axis NBI operation (stiff EP density profile gradients nearby the magnetic axis). Using the simulation model that reproduces the observed Alfvén activity, we extend the study to analyze NBI performance within a theoretical framework. It shows that increasing NBI voltage (which raises EP energy) leads to a degradation in NBI performance for a given power (related to EP β and their density). To achieve better NBI operation, higher voltage must be balanced with lower injection power, ensuring stable AEs while keeping the same EP β.

Plasma confinement↗
DOE OSTI2026

Comparing ICME simulations with scaled laboratory experiment

In stellar physics and astrophysics, numerical simulations and laboratory experiments are often compared to observational data to support their representation of the real world. However, there is also merit in comparing numerical simulations to properly scaled experiments, especially when the experiment and the simulation are both emulating the solar phenomena. Confirming the credibility of scaled experiments and their scaling with well-validated models is important to expand our knowledge of the associated physical phenomena. This is significant because experiments and simulations can be performed frequently, whereas observations may be limited by location, field of view, and missing data. In this work, we use the Alfvén Wave Solar atmosphere Model, a well-validated magnetohydrodynamic model, to simulate an interplanetary coronal mass ejection (ICME) and compare it to an experiment which provides a scaled analog to a physical ICME. The experiment was performed on the Big Red Ball facility and scaled using dimensionless parameters such as plasma β and magnetosonic Mach number to reproduce the main structure of an ICME. We compare the model-simulated temperature, density, and magnetic field to those from the experiment, as well as the scaling parameters used in the experiment, to those calculated from the simulation. This comparison is performed to further justify the scaling arguments made by the experiment. Additionally, the comparison would lead to the development of stronger scaling arguments for future experiments.

Bryant, K. [University of Michigan, Ann Arbor, MI ↗
DOE OSTI2025

The Quench Protection of Direct Wind Magnets

The direct wind magnets have several unique features. The conductor layout is characterized by thin single conductors wound over long lengths and separated from the other turns with thick layer of insulation. Moreover, there are several layers of winding. The strong non-uniformity in the magnetic field distribution results in different current and thermal margins for quenching in different regions. The quench protection modeling involves multi-physics coupling between electrical, thermal and magnetic transient. The material properties vary over time and location. This coupled with the complexity in the geometry adds to long computation times. The objective of this study is to understand the quench propagation in an EIC high inductance direct wind magnet which enables one to design appropriate quench protection strategies to protect the magnet. An advanced modeling program has been developed to simulate direct wind magnet quench. In conclusion, we validate theoretical simulations with experimental data.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Modelling pulsed field magnetization of iron-based bulk superconductors

Abstract Bulk superconductors can be used as super-strength quasi-permanent magnets capable of providing magnetic flux densities considerably superior to conventional permanent magnets. This makes them attractive for several engineering applications that rely on strong magnetic fields like rotating machines, NMR/MRI and magnetic drug delivery systems. Recently, the authors reported a record trapped magnetic field in an iron-based bulk superconductor: 2.83 T was trapped in potassium-doped barium iron arsenide (Ba, K)Fe 2 As 2 (or Ba122) at 5 K. Of particular significance is that the strength and temporal stability of this magnetic field exceeds the requirements of MRI machines, indicating iron-based bulks can now perform at levels demanded by engineering applications. One crucial challenge for their practical use, however, is the need to apply and remove an external magnetic field to magnetize them. Pulsed field magnetization (PFM) shows great promise as a practical method of magnetizing bulks, but the process generates heat in the bulk that is detrimental to its superconducting performance and ability to act as a super-strength magnet. In this paper, coupled electromagnetic–thermal numerical models are used to simulate the PFM of iron-based bulk superconductors. Here we focus on the recent-record-breaking, fine-grain polycrystalline K-doped Ba122 bulks. The impact that the specific J c ( B ) characteristics and thermal properties of the Ba122 material—all of which have been experimentally measured from state-of-the-art samples—have on the magnetic flux dynamics and thermal behaviour during PFM, including the final trapped field, is investigated. We show that because the thermal properties are similar to those of REBa 2 Cu 3 O 7 −δ bulks, a similar response to pulsed fields is obtained. A maximum trapped field of ∼0.81 T (∼43.3% of the maximum trapped field capability under ideal, field-cooling conditions) was simulated at 5 K, with a magnetization efficiency of ∼54%. The modelling framework provides a fast and flexible tool for optimising the practical PFM process at different operating temperatures to maximise the trapped field in state-of-the-art Ba122 bulks and to guide the design of future experiments.

bulk superconductors↗