Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fields simulating electromagnetic fields”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4

Full-wave simulations on helicon and parasitic excitation of slow waves near the edge plasma

Helicon waves are thought to be promising in various tokamaks, such as DIII-D, because they can penetrate reactor-grade high-density cores and drive the off-axis current with higher efficiency. In the frequency regime ~476 MHz, both slow electrostatic and fast electromagnetic helicon waves can coexist in DIII-D. If the antenna parasitically excites the slow mode, these waves can propagate along the magnetic field line into the scrape-off layer (SOL). Although the importance of the misalignment of the Faraday screen and the electron density in the SOL on the excitation and propagation of slow modes is well known, the conditions for minimizing slow mode excitation have yet to be optimized. Using the Petra-M simulation code in the 2D domain, we analyze the effects of the misalignment of the antenna in the poloidal direction, the misalignment of the Faraday screen in the toroidal direction, and the density in front of the antenna on slow mode generation. Our results suggest that the misalignment of the Faraday screen is a critical factor in reducing the slow mode and that the misalignment angle should be below ~5° to minimize the slow wave excitation. When the electron density is higher than 3.5 x 10 18 m –3 in the SOL, the generation of the slow mode from the antenna is minimized and unaffected by the misalignment of the Faraday screen.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

HiTEM

SAND2025-03438O HiTEM3d v1.0 is a user-friendly software designed to simulate how electromagnetic signals move through large geological environments, particularly around wellbores. It helps geophysicists and resource explorers understand subsurface conditions more effectively and affordably. By simplifying complex geological models, HiTEM3d allows users to visualize and analyze electromagnetic fields, making it easier to identify valuable resources like minerals or water. Its efficient processing capabilities ensure quick results, making it an essential tool for anyone involved in geophysical research or exploration. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Weiss, Chester [Sandia National Lab. (SNL-CA), Liv↗

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)↗

Higher-form symmetry and chiral transport in real-time Abelian lattice gauge theory

We study classical lattice simulations of theories of electrodynamics coupled to charged matter at finite temperature, interpreting them using the higher-form symmetry formulation of magnetohydrodynamics (MHD). We compute transport coefficients using classical Kubo formulas on the lattice and show that the properties of the simulated plasma are in complete agreement with the predictions from effective field theories. In particular, the higher-form formulation allows us to understand from hydrodynamic considerations the relaxation rate of axial charge in the chiral plasma observed in previous simulations. A key point is that the resistivity of the plasma – defined in terms of Kubo formulas for the electric field in the 1-form formulation of MHD – remains a well-defined and predictive quantity at strong electromagnetic coupling. However, the Kubo formulas used to define the conventional conductivity vanish at low frequencies due to electrodynamic fluctuations, and thus the concept of the conductivity of a gauged electric current must be interpreted with care.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Grid Resiliency with a 100% Renewable Microgrid

San Diego Gas & Electric Company (SDG&E) installed America’s first and largest utility-scale microgrid in Borrego Springs in 2013. The first generation Borrego Springs Microgrid utilized diesel generators to form and stabilize the microgrid island, with support from grid-scale batteries and local solar photovoltaic (PV) generation. In this project, SDG&E in partnership with National Renewable Energy Laboratory (NREL) demonstrated through modeling, simulation and utility field testing that blackstart and islanding of the microgrid can be led with 100% renewable, inverter based resources (IBRs), to help reduce community reliance on conventional generation resources. Through equipment upgrades, grid-forming island leader capability was transitioned to a battery IBR instead of the Borrego Springs Microgrid diesel generators. A new microgrid controller was integrated to the microgrid and programmed to control and manage multiple energy storage systems. Synchrophasor and other power quality data verified autonomous, high-speed response of the IBRs through blackstart, islanding, and load step testing. Results of project field evaluations provide distribution systems operators (DSO) with increased confidence that renewable, IBR can replace traditional generators to blackstart and island microgrids and rapidly establish stable island frequency with rapid changes in peak power demand. Importantly, the project validated the integration feasibility of a distributed energy resource management system (DERMS) controller that manages multiple grid-forming and grid-following IBRs, establishing a standard design interface to reduce the complexity of integrating new DERs in the future and supporting replication by the industry. As a result of learnings in this project, SDG&E has implemented the microgrid controller strategy at multiple other microgrid sites, thereby validating the replicability of the solution. Hardware-in-the-loop (HIL) simulations including power and controller HIL hardware — along with electromagnetic transient (EMT) simulations of Borrego Springs Microgrid —informed adjustments to inverter parameters and were important to characterize the performance of the IBRs in relevant operating conditions before deployment. The EMT and HIL simulations of islanding the entire community are important contributions in providing confidence in IBR performance prior to future islanding of the community in the field. High-fidelity EMT and/or HIL simulation of IBRs can de-risk field operations, and its relevance and importance as a tool is increasing as distribution grids and microgrids become more complex and dynamic with an increasing proportion of renewable generation, distributed energy storage, and two-way power and energy flows.

24 POWER TRANSMISSION AND DISTRIBUTION↗

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

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

79 ASTRONOMY AND ASTROPHYSICS↗

FEM Analysis of Hybrid LTS/HTS Cos-Theta Dipole Magnet With Heterogeneous Cable Model

The Lawrence Berkeley National Laboratory (LBNL) and the National High Magnetic Field Laboratory (NHMFL) have published results on Bi-2212 superconductive magnets realized and tested in the canted cosine-theta and solenoid designs, respectively. Fermilab is now preparing for the assembly of the first Bi-2212 stress-managed cosine-theta insert magnet. The insert will be part of the first hybrid cosine-theta magnet made of Nb$_3$Sn outer layers within the US-MDP effort to reach a 20 T bore field. This paper presents the analytical analysis of the cosine-theta Nb$_3$Sn/Bi-2212 hybrid magnet. We report the parameters, logic, and implementation method of the 2D electromagnetic and mechanical FEM analysis of the LTS/HTS hybrid magnet. Results from a detailed heterogeneous model are compared to the homogeneous model implemented in the past. A Python code has been developed to simulate the current degradation due to stresses in the detail-modeled conductor areas. The current degradation has been introduced in the simulation dynamics for the HTS conductor as an iteration process, updating the input load of Lorentz forces of the energization at each step. The magnetic and mechanical analysis results of the 2D cosine-theta LTS/HTS dipole magnet have been described and analyzed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Emergent electric field induced by dissipative sliding dynamics of domain walls in a Weyl magnet

The dynamic motion of topological defects in magnets induces an emergent electric field, as exemplified by the continuous flow of skyrmion vortices. However, the electrodynamics underlying this emergent field remains poorly understood. In this context, magnetic domain walls—one-dimensional topological defects with two collective modes, sliding and spin-tilt—offer a promising platform for exploration. Here we demonstrate that the dissipative motion of domain walls under oscillatory current excitation generates an emergent electric field. We image domain patterns and quantify the domain-wall length under applied magnetic fields in mesoscopic devices based on the magnetic Weyl semimetal NdAlSi. These devices exhibit exceptionally strong domain-wall scattering and a pronounced emergent electric field, as observed in the imaginary component of the complex impedance. Spin dynamics simulations reveal that domain-wall sliding dominates over spin-tilting, in which the phase delay of the domain-wall motion with respect to the driving force impacts the emergent electric field. Our findings establish domain-wall dynamics as a platform for studying emergent electromagnetic fields and motivate further investigations of the coupled motion of magnetic solitons and conduction electrons.

Yamada, Rinsuke [The University of Tokyo, Japan]↗

Geant4 Event Biasing and Fast Simulation

Geant4 offers advanced event biasing techniques to significantly accelerate simulations involving rare events. Various biasing methods, such as leading particle selection, cross-section biasing, radioactive decay enhancement, and bremsstrahlung splitting, enable efficient event sampling, though they require careful handling. Additionally, Geant4 provides a Fast Simulation Interface, allowing the replacement of standard processes in specific region and for selected particles, enabling faster execution or external code integration. Applications of fast simulation include electromagnetic shower modeling in calorimeters, machine learning inference, and offloading tasks to specialized hardware like GPUs, making Geant4 a powerful tool for computationally demanding simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

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↗

Assessing time-dependent temperature profile predictions using reduced transport models for high performing NSTX plasmas

Time-dependent, predictive simulations were performed with the 1.5D tokamak integrated modeling code TRANSP on a large set of well-analyzed, high performing discharges from the National Spherical Torus Experiment (NSTX) in order to evaluate how well modern reduced transport models can reproduce experimentally observed temperature profiles in spherical tokamaks. Overall, it is found that simulations using the Multi-Mode Model (MMM) more consistently agree with the NSTX observations than those using the Trapped Gyro-Landau Fluid (TGLF) model, despite TGLF requiring orders of magnitude greater computational cost. When considering all examined discharges, MMM has median overpredictions of electron temperature (T e ) and ion temperature (T i ) profiles of 28% and 27%, respectively, relative to the experiment. TGLF overpredicts T e by 46%, with much larger variance than MMM, and underpredicts T i by 25%. As the ratio of kinetic to magnetic field pressure (β) is increased across NSTX discharges, TGLF predicts lower T e and significant flattening of the T i profile, conflicting with NSTX observations. When using an electrostatic version of TGLF, both T e and T i are substantially overpredicted, underscoring the importance of electromagnetic turbulence in the high β spherical tokamak regime. Additionally, calculations with neural net surrogate models for TGLF were performed outside of TRANSP with a time slice flux matching transport solver, finding better agreement with experiment than the TRANSP simulations, highlighting the impact of different transport solvers and simulation techniques. Furthermore, the reasonable agreement with experiment of temperature profiles predicted by MMM motivates a more detailed examination of the sensitivities of the TRANSP simulations with MMM to different NSTX plasma regimes in a companion paper (Lestz et al Plasma Phys. Control. Fusion 67 105030), in preparation for self-consistent, time-dependent predictive modeling of NSTX-U scenarios.

Lestz, J. B. [General Atomics, San Diego, CA (Unit↗

Spindown of Pulsars Interacting with Companion Winds: The Impact of Magnetospheric Compression

Pulsars in binary systems with strong companion winds can have the magnetopause separating their magnetosphere from the wind located well within their light cylinder. This bow-like enclosure effectively creates a waveguide that confines the pulsar's electromagnetic fields and can significantly alter its spindown. In this paper, we study the spindown of compressed pulsar magnetospheres in such systems. We parameterize the confinement as the ratio between the equatorial position of the magnetopause (or standoff distance) R m and the pulsar's light cylinder R LC . Using particle-in-cell simulations, we quantify the pulsar spindown for a range of compressions, R m /R LC = 1/3–1, and inclination angles, χ = 0°...90°, between magnetic and rotation axes. Our strongly confined models (R m /R LC = 1/3) show two distinct limits. For χ = 0°, the spindown of a compressed pulsar magnetosphere is enhanced by approximately a factor of three compared to an isolated pulsar due to the increased number of open magnetic field lines. Conversely, for χ = 90°, the compressed pulsar spins down at less than 40% of the rate of an isolated reference pulsar due to the mismatch between the pulsar wind stripe wavelength and the waveguide size. We apply our analysis to the 2.77 s oblique rotator (χ = 60°) in the double-pulsar system PSR J0737-3039. With the numerically derived spindown estimate, we constrain its surface magnetic field to B* ≈ (7.3 ± 0.2) × 10 11 G. We discuss the time modulation of its period derivative, the effects of compression on its braking index, and implications for the radio eclipse in PSR J0737-3039.

79 ASTRONOMY AND ASTROPHYSICS↗

Spatially Resolved Near Field Spectroscopy of Vibrational Polaritons at the Small N Limit

Vibrational polaritons, which have been primarily studied in Fabry–Pérot cavities with a large number of molecules (N ~ 10 6 –10 10 ) coupled to the resonator mode, exhibit various experimentally observed effects on chemical reactions. However, the exact mechanism is elusively understood from the theoretical side, as the large number of molecules involved in an experimental strong coupling condition cannot be represented completely in simulations. This discrepancy between theory and experiment arises from computational descriptions of polariton systems typically being limited to only a few molecules, thus failing to represent the experimental conditions adequately. To address this mismatch, we used surface phonon polariton (SPhP) resonators as an alternative platform for vibrational strong coupling. SPhPs exhibit strong electromagnetic confinement on the surface and thus allow for coupling to a small number of molecules. As a result, this platform can enhance nonlinearity and slow down relaxation to the dark modes. In this study, we fabricated a pillar-shaped quartz resonator and then coated it with a thin layer of cobalt phthalocyanine (CoPc). By employing scattering-type scanning near-field optical microscopy (s-SNOM), we spatially investigated the dependency of vibrational strong coupling on the spatially varying electromagnetic field strength and demonstrated strong coupling with 38,000 molecules only–reaching to the small N limit. Through s-SNOM analysis, we found that strong coupling was observed primarily on the edge of the quartz pillar and the apex of the s-SNOM tip, where the maximum field enhancement occurs. In contrast, a weak resonance signal and lack of coupling were observed closer to the center of the pillar. This work demonstrates the importance of spatially resolved polariton systems in nanophotonic platforms and lays a foundation to explore polariton chemistry and chemical dynamics at the small N limit–one step closer to reconcile with high-level quantum calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Symplectic neural network and its application to charged particle dynamics in electromagnetic fields

Recently, machine learning models have shown many successes in various applications in science and technology. In this work, we focus on the charged particle dynamics, with the development of a class of symplectic neural networks, including a linear version, SympMat, and a nonlinear version, HénonNet. Both are designed to preserve the structure of Hamiltonian systems. We show that they can be used to model relevant Hamiltonian systems of interest in plasma physics and astrophysics, for linear and nonlinear charged particle dynamics, with the potential to bridge multi-scale simulations. These symplectic neural networks are adapted to the applications in plasma simulations and particle-wave interaction with parametric dependence and periodicity, where we have investigated their performance and accuracy. In particular, SympMat is shown to outperform the traditional Boris particle pusher down to the sub-gyroperiod scale in the case of charged particles in uniform magnetic fields. HénonNet successfully predicts the hot electron distribution, which is validated against theoretical results. These results highlight the potential of symplectic neural networks as a trajectory integrator for particle-in-cell simulations or a fast surrogate to replace conventional numerical schemes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Phase II Field Demonstration at Lansing Smith Generating Plant, Southport, Florida (Final Report)

The Final Technical Report: Field Test Design and Pressure Management Strategies for Phase II Field Demonstration of Optimal Pressure Control, Plume Management, and Produced Water Strategies presents the culmination of multi-year efforts under the U.S. Department of Energy’s Brine Extraction and Storage Test (BEST) program, led by EPRI in partnership with Gulf Power at the Plant Smith site near Panama City, Florida. The project was designed to evaluate and demonstrate the technical feasibility of managing subsurface pressures and fluid movement associated with large-scale CO₂ injection, using low-salinity water as a proxy fluid. Through a combination of field injection testing, reservoir modeling, and optimization studies, the research team developed and refined active and passive brine extraction strategies aimed at controlling injection-induced pressure buildup, mitigating risks of fault activation, and managing plume migration. The field demonstration incorporated a new injection well (TIW-2), a new monitoring/extraction well (TEMW-A), and an existing well (TIW-1) repurposed for passive pressure relief. Complementary geophysical monitoring was designed to track plume development and evaluate the effectiveness of pressure management strategies. The report details the integrated workflow encompassing wellfield development, regulatory permitting, model calibration, and survey design. It includes a comprehensive electromagnetic (EM) modeling and inversion study used to develop a cost-effective, time-lapse geophysical monitoring plan capable of imaging the evolving low-salinity plume within the highly saline Lower Tuscaloosa formation. Reservoir simulation results guided the design of the 17-month injection program and the timing of active extraction to maintain formation pressures below a hypothetical fault reactivation threshold. Supporting analyses evaluated potential injectivity risks related to clay fines migration and geochemical incompatibility, identifying practical mitigation measures such as salinity and pH control. The final design recommends a series of crosswell EM surveys—conducted before, during, and after injection—as the most effective and economical approach for plume imaging, supplemented by continuous downhole pressure and flow monitoring. Collectively, the study provides a field-ready framework for cost-effective pressure management and monitoring in support of future CO₂ storage projects.

01 COAL, LIGNITE, AND PEAT↗

Reduced-action-integral approach for photon-photon interactions in vacuum

Electromagnetic waves propagating through vacuum can polarize virtual electron–positron pairs; this polarization, in turn, nonlinearly modifies their propagation. A semi-classical nonlinear wave equation describing the propagation is derived from the Euler–Heisenberg Lagrangian density, which captures vacuum polarization effects up to the one-loop level. In this article, we present a reduced-actionintegral approach that enables rapid modeling of nonlinear phenomena arising from the Euler– Heisenberg Lagrangian. Application of the variational principle to the reduced action provides equations of motion for familiar light-pulse parameters, such as spot size, phase, polarization, and phase-front curvature, without requiring full-field simulations. Three examples demonstrate the utility of the approach: phase modulation, birefringence, and frequency mixing.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rank-Limiting Strategies for Optimizing Tensor-Train Finite-Difference Time-Domain Simulations

We introduce rank-limiting strategies to optimize tensor-train decompositions for three-dimensional finite-difference time-domain simulations using the relationship between the tensors and their specific dimensionality. These include the use of hard caps on the inner ranks of the tensor train decomposition and the use of a group rounding algorithm taking into account all field components simultaneously. Here, several numerical examples are considered to verify the efficacy of the proposed optimization strategies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Off-shell pion properties: Electromagnetic form factors and light-front wave functions

The off-shell pion electromagnetic form factors are explored with corresponding off-shell light-front wave functions modeled by constituent quark and antiquark. We apply the Mandelstam approach for the microscopic computation of the form factors relating the model parameters with the pion decay constant and charge radius. Analyzing the existing data on the cross sections for the Sullivan process, 1 H(e,e',π + )⁢n, Charged pion form factor between Q 2 = 0.60, and 2.45 GeV 2 . I. Measurements of the cross section for the 1 H⁡(e,e'⁢π + )⁢n reaction, we extract the off-shell pion form factor using the relation derived from the generalized Ward-Takahashi identity for the pion electromagnetic current. They are compared with our previous results from exactly solvable manifestly covariant model of a (3+1)-dimensional fermion field theory. We find that the adopted constituent quark model reproduces the extracted off-shell form factor F 1 ⁡(Q 2 ,t) from the experimental data within a few percent difference and matches well with our previous theoretical simulation which exhibits a variation of about 10% for the extracted off-shell pion form factor g(Q 2 ,t). We also identify the pion valence parton distribution function (PDF) and transverse momentum distribution (TMD) in terms of the light-front wave function and discuss their off-shell properties.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗