Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Methods in electromagnetism”

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

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

At least 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

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↗
DOE OSTI2025

A Participation Factor-Based Approach for Defining the EMT Model Boundary for Power System Simulations with Inverter-Based Resources

The increasing penetration of inverter-based resources (IBRs) introduces new challenges to power system simulations, particularly with the emergence of fast electromagnetic transient (EMT) dynamics and sub-synchronous oscillations (SSO) that require time-consuming EMT simulations. To reduce the time cost for simulating a large-scale power grid with IBRs, this paper proposes a novel participation factor-based approach for defining a critical zone for detailed EMT modeling and simulations, which includes the IBRs, synchronous generators, and the network components participating significantly in simulated contingencies. Both model-based and response-based methods are introduced for the estimation of participation factors (PFs). The case study on the 240-bus Western Electricity Coordinating Council (WECC) system demonstrates that the EMT zone determined by the proposed approach can effectively capture power system dynamics involving IBRs.

EMT simulation↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Optimal Transport for $e/\pi^0$ Particle Classification in LArTPC Neutrino Experiments

The efficient classification of electromagnetic activity from $\pi^0$ and electrons is a notoriously challenging problem in the reconstruction of neutrino interactions in Liquid Argon Time Projection Chamber (LArTPC) detectors. We address this problem using the mathematical framework of Optimal Transport (OT), which has been successfully employed for event classification in other HEP contexts and is ideally suited to the high-resolution calorimetry of LArTPCs. Using a publicly available simulated dataset from the MicroBooNE collaboration, we show that OT methods achieve state-of-the-art reconstruction performance in $e/\pi^0$ classification. The success of this first application indicates the broader promise of OT methods for LArTPC-based neutrino experiments. This work motivates integrating OT in the reconstruction frameworks of LArTPC experiments such as SBN and DUNE more broadly. Since $\pi^0$s are a significant background for both oscillation experiments and BSM searches, OT can lead to sizeable improvements in the selection efficiency for these analyses by introducing a novel method with which to achieve $\pi^0$ rejection.

Caratelli, David [UC, Santa Barbara]↗
DOE OSTI2021

Constructing a new predictive scaling formula for ITER's divertor heat-load width informed by a simulation-anchored machine learning

Understanding and predicting divertor heat-load width λq is a critically important problem for an easier and more robust operation of ITER with high fusion gain. Previous predictive simulation data for λ q using the extreme-scale edge gyrokinetic code XGC1 [S. Ku et al., Phys. Plasmas 25, 056107 (2018)] in the electrostatic limit under attached divertor plasma conditions in three major US tokamaks [C. S. Chang et al., Nucl. Fusion 57, 116023 (2017)] reproduced the Eich and Goldston attached-divertor formula results [formula #14 in T. Eich et al., Nucl. Fusion 53, 093031 (2013) and R. J. Goldston, Nucl. Fusion 52, 013009 (2012)] and furthermore predicted over six times wider λ q than the maximal Eich and Goldston formula predictions on a full-power (Q = 10) scenario ITER plasma. After adding data from further predictive simulations on a highest current JET and highest-current Alcator C-Mod, a machine learning program is used to identify a new scaling formula for λ q as a simple modification to the Eich formula #14, which reproduces the Eich scaling formula for the present tokamaks and which embraces the wide λ q XGC for the full-current Q = 10 ITER plasma. Additionally, the new formula is then successfully tested on three more ITER plasmas: two corresponding to long burning scenarios with Q = 5 and one at low plasma current to be explored in the initial phases of ITER operation. The new physics that gives rise to the wider λ q XGC is identified to be the weakly collisional, trapped-electron-mode turbulence across the magnetic separatrix, which is known to be an efficient transporter of the electron heat and mass. Electromagnetic turbulence and high-collisionality effects on the new formula are the next study topics for XGC1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Isotope effect and critical magnetic fields of superconducting YH 6 : A Migdal-Eliashberg theory approach

The emergence of near-ambient temperature superconductivity under pressure in metal hydride systems has motivated a desire to further understand such remarkable properties, specifically critical magnetic fields. Here, YH 6 is suggested to be a departure from conventional superconductivity, due to apparent anomalous behavior. Using density functional calculations in conjunction with Migdal-Eliashberg theory we show that in YH 6 the critical temperature and the isotope effect under pressure, as well as the high critical fields, are consistent with strong-coupling conventional superconductivity, a property anticipated to extend to other related systems. Furthermore, strong-coupling corrections occur to the expected BCS values for the isotope effect coefficient (α), Ginzburg-Landau parameter [κ 1 ⁡(T)], London penetration depth [λ L (T)], electromagnetic coherence length [ξ⁡(T)], and the energy gap (Δ 0 ).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Optimal transport for 𝑒/𝜋 0 particle classification in LArTPC neutrino experiments

The efficient classification of electromagnetic activity from 𝜋 0 and electrons remains an open problem in the reconstruction of neutrino interactions in liquid argon time projection chamber (LArTPC) detectors. We address this problem using the mathematical framework of optimal transport (OT), which has been successfully employed for event classification in other high energy physics contexts and is ideally suited to the high-resolution calorimetry of LArTPCs. Using a publicly available simulated dataset from the MicroBooNE Collaboration, we show that OT methods achieve state-of-the-art reconstruction performance in 𝑒/𝜋 0 classification. The success of this first application indicates the broader promise of OT methods for LArTPC-based neutrino experiments.

Neutrino detection↗
DOE OSTI2021

Computational Fluid Dynamics Simulations to Support Efficiency Improvements in Aluminum Smelting Process

Smelting is broadly described as the extraction of a metal from its ore. In the United States, aluminum is commonly produced by smelting alumina in bauxite using the Hall-Héroult process. Optimization of equipment and processes in conventional smelting is crucial to enhancing process efficiency and productivity, is necessary for improving the techno-economic feasibility, which directly manifests as the growth of the American economy. To achieve optima, insightful data on the multiphysics phenomena that are inherent to the process must be obtained through physical investigation or high-fidelity numerical simulations. The resolution of relevant scales in time and space for smelting operations requires intensive, high-performance computing (HPC) simulations. Hostile operating conditions limit physical data acquisition to specific techniques; therefore, these data do not describe the multiscale interaction of simultaneous effects. Fortunately, in recent decades, significant advancements in computing hardware and computational methods have made the numerical resolution of such a complex process possible. In this study, a high-fidelity simulation of aluminum smelting was performed using an open-source tool, OpenFOAM, which analyzed many parameters characteristic to underlying phenomena. A multiphysics model based on the Eulerian-Eulerian multifluid approach was adopted. This model can resolve critical issues in the electrolytic smelting of aluminum, such as bubbling of carbon dioxide from the anode(s), magnetohydrodynamics from electromagnetic effects, ionic dissolution of the alumina in the electrolyte, and the evolution of thermal profiles. This study provides valuable connectivity for characteristic data that can direct the future designs of efficient smelters. A basic framework to model and simulate the smelting process using OpenFOAM is presented for user modification in keeping with process development. Of relevance to the flow field, a detailed investigation of vortices produced by bubble motion and electromagnetics is discussed, along with their impact on the evolution of thermal profiles. The predictions show small-scale vortices in the clearance between the anode and cathode caused by magnetic forces. Predictions also indicate relatively large-scale vortices in the inter-anode space resulting from carbon dioxide rising through the electrolytic flow field. The formation of vortices at the edges of anodes was shown to direct alumina charged by the feeder to the bottom of the anodes, thus preventing the entrapment of gas bubbles in the periphery of the bottom of the anode. Symmetry was observed in the location of cold spots in the electrolytic mixture in the vicinity of the feeder. Cold spots were also observed in the clearance between the anode and cathode due to the flow’s transmission of unconverted alumina to this region.

36 MATERIALS SCIENCE↗
DOE OSTI2021

A Fabrication Strategy for Reconfigurable Millimeter–Scale Metamaterials

Rather than depending on material composition to primarily dictate performance metrics, metamaterials can leverage geometry to achieve specific properties of interest. For example, reconfigurable metamaterials have enabled programmable shape transformations, tunable mechanical properties, and energy absorption. While several methods exist to fabricate such structures, they often place severe restrictions on manufacturing materials, or require significant manual assembly. Moreover, these arrays are typically composed of unit cells that are either macro–scale or micro–scale in dimension. Here, the fabrication gap is bridged, and laminate manufacturing is used to develop a method for designing reconfigurable metamaterials at the millimeter–scale, that is compatible with a wide range of materials, and that requires minimal manual assembly. In addition to showing the versatility of this fabrication method, how the use of laminate manufacturing affects the behavior of these multi–component arrays is also characterized. To this end, a numerical model that captures the deformations exhibited by the structures is developed, and an analytic model that predicts the strain of the structure under compressive stress is built. Overall, this approach can be leveraged to develop millimeter–scale metamaterials for applications that require reconfigurable materials, such as in the design of tunable acoustics, photonic waveguides, and electromagnetic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

DQ Admittance Extraction for Inverter-Based Resources

The power grid industry is pushing for electromagnetic transient (EMT)-based studies for generation interconnection and planning process due to high penetrations of inverter-based resource (IBRs). Vendor-specific and site-specific black-box IBR models are preferred in those simulation studies. For small-signal analysis, measurement-based admittance models are necessary. In this paper, we demonstrate the extraction of frequency-domain dq frame IBR models. These linear models are specific to operating conditions. We demonstrate two extraction methods: frequency scan and a step response-based method. The latter relies on converting time-domain responses to Laplace-domain expressions via eigensystem realization analysis (ERA). Both can lead to dq admittance representation and the latter is time saving.

admittance model↗
DOE OSTI2020

Plasmon-Mediated Intramolecular Methyl Migration with Nanoscale Spatial Control

Plasmonic materials interact strongly with light to focus and enhance electromagnetic radiation down to nanoscale volumes. Due to this localized confinement, materials that support localized surface plasmon resonances are capable of driving energetically unfavorable chemical reactions. In certain cases, the plasmonic nanostructures are able to preferentially catalyze the formation of specific photoproducts, which offers an opportunity for the development of solar-driven chemical synthesis. Here, using plasmonic environments, we report inducing an intramolecular methyl migration reaction, forming 4-methylpyridine from N-methylpyridinium. Using both experimental and computational methods, we were able to confirm the identity of the N-methylpyridinium by making spectral comparisons against possible photoproducts. This reaction involves breaking a C–N bond and forming a new C–C bond, highlighting the ability of plasmonic materials to drive complex and selective reactions. Additionally, we observe that the product yield depends strongly on optical illumination conditions. This is likely due to steric hindrance in specific regions on the nanostructured plasmonic substrate, providing an optical handle for driving plasmonic catalysis with spatial specificity. Furthermore, this work adds yet another class of reactions accessible by surface plasmon excitation to the ever-growing library of plasmon-mediated chemical reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

First characterisation of the MAGO cavity, a superconducting RF detector for kHz–MHz gravitational waves

Heterodyne detection using microwave cavities is a promising method for detecting high-frequency gravitational waves (GWs) or ultralight axion dark matter. In this work, we report on studies conducted on a spherical 2-cell cavity developed by the MAGO collaboration for high-frequency GWs detection. Although fabricated around 20 years ago, the cavity had not been used since. Due to deviations from the nominal geometry, we conducted a mechanical survey and performed room-temperature plastic tuning. Measurements and simulations of the mechanical resonances and electromagnetic properties were carried out, as these are critical for estimating the cavity’s GW coupling potential. Based on these results, we plan further studies in a cryogenic environment. The cavity characterisation does not only provide valuable experience for a planned physics run but also informs the future development of improved cavity designs.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Laser-plasma ion beam booster based on hollow-channel magnetic vortex acceleration

Laser-driven ion acceleration provides ultrashort, high-charge, low-emittance beams, which are desirable for a wide range of high-impact applications. Yet after decades of research, a significant increase in maximum ion energy is still needed. This paper introduces a quality-preserving staging concept for ultraintense ion bunches that is seamlessly applicable from the nonrelativistic plasma source to the relativistic regime. Full three-dimensional particle-in-cell simulations prove robustness and capture of a high-charge proton bunch, suitable for readily available and near-term laser facilities.

43 PARTICLE ACCELERATORS↗
DOE OSTI2025

Multi-frequency progressive refinement for learned inverse scattering

Interpreting scattered acoustic and electromagnetic wave patterns is a computational task that enables remote imaging in a number of important applications, including medical imaging, geophysical exploration, sonar and radar detection, and nondestructive testing of materials. However, accurately and stably recovering an inhomogeneous medium from far-field scattered wave measurements is a computationally difficult problem, due to the nonlinear and non-local nature of the forward scattering process. We design a neural network, called Multi-Frequency Inverse Scattering Network (MFISNet), and a training method to approximate the inverse map from far-field scattered wave measurements at multiple frequencies. We consider three variants of MFISNet, with the strongest performing variant inspired by the recursive linearization method — a commonly used technique for stably inverting scattered wavefield data — that progressively refines the estimate with higher frequency content. MFISNet outperforms past methods in regimes with high-contrast, heterogeneous large objects, and inhomogeneous unknown backgrounds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Penetration Through Slots in Overmoded Cavities

A resonant cavity undergoes three distinct behaviors with increasing frequency: 1) fundamental modes, localized in frequency with well defined modal distribution; 2) undermoded region, where modes are still separated, but are sufficiently perturbed by small imperfections that their spectral positions (and distributions) are statistical in nature; and 3) overmoded region, where modes overlap, field distributions follow stochastic distributions, and the slot acts as if in free space. Understanding the penetration through slots in the overmoded region is of great interest, and is the focus of this article. Since full-wave solvers may not be able to provide a timely answer for very high frequencies due to a lack of memory and/or computation resources, we develop bounding methods to estimate worst-case average and maximum fields within the cavity. Finally, after discussing the bounding formulation, we compare its results to full-wave simulations at the first, second, and third resonance supported by the slot in the case of a cylindrical cavity. Note that the bounding formulation indicates that results are nearly independent of cavity shape: only the cavity volume, frequency, and cavity quality factor affect the overmoded region, making this formulation a powerful tool to assess electromagnetic interference and electromagnetic compatibility effects within cavities.

overmoded cavity↗
DOE OSTI2023

Nucleon Resonance Electroexcitation Amplitudes and Emergent Hadron Mass

Understanding the strong interaction dynamics that govern the emergence of hadron mass (EHM) represents a challenging open problem in the Standard Model. In this paper we describe new opportunities for gaining insight into EHM from results on nucleon resonance (N*) electroexcitation amplitudes (i.e., γvpN* electrocouplings) in the mass range up to 1.8 GeV for virtual photon four-momentum squared (i.e., photon virtualities Q2) up to 7.5 GeV2 available from exclusive meson electroproduction data acquired during the 6-GeV era of experiments at Jefferson Laboratory (JLab). These results, combined with achievements in the use of continuum Schwinger function methods (CSMs), offer new opportunities for charting the momentum dependence of the dressed quark mass from results on the Q2-evolution of the γvpN* electrocouplings. This mass function is one of the three pillars of EHM and its behavior expresses influences of the other two, viz. the running gluon mass and momentum-dependent effective charge. A successful description of the Δ(1232)3/2+ and N(1440)1/2+ electrocouplings has been achieved using CSMs with, in both cases, common momentum-dependent mass functions for the dressed quarks, for the gluons, and the same momentum-dependent strong coupling. The properties of these functions have been inferred from nonperturbative studies of QCD and confirmed, e.g., in the description of nucleon and pion elastic electromagnetic form factors. Parameter-free CSM predictions for the electrocouplings of the Δ(1600)3/2+ became available in 2019. The experimental results obtained in the first half of 2022 have confirmed the CSM predictions. We also discuss prospects for these studies during the 12-GeV era at JLab using the CLAS12 detector, with experiments that are currently in progress, and canvass the physics motivation for continued studies in this area with a possible increase of the JLab electron beam energy up to 22 GeV. Such an upgrade would finally enable mapping of the dressed quark mass over the full range of distances (i.e., quark momenta) where the dominant part of hadron mass and N* structure emerge in the transition from the strongly coupled to perturbative QCD regimes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Polarization‐Engineered Near‐Field Generation Using a Hybrid Tip–Antenna System

Precise control of light polarization at the nanoscale is critical for accessing chiral optical responses and manipulating spin–photon interactions in advanced materials. Yet, conventional scattering-type near-field probes predominantly generate out-of-plane linear polarization and offer little control over phase or polarization state. Here, in this study, we introduce a polarization-engineered near-field methodology based on a combined metallic tip and planar dipole nanoantenna system. Using full-wave electromagnetic simulations, we show that the tip acts as a vertically oriented plasmonic resonator, while the antenna supports an in-plane dipolar mode. By tuning the tip–antenna geometry and tip height, the two orthogonal field components attain comparable amplitudes and a controllable ∼90° phase offset, producing circularly polarized nano-light in the antenna gap. The proposed system effectively functions as a nanoscale quarter-wave plate, converting linearly polarized illumination into circularly polarized hotspots without external polarization optics. This method establishes an experimentally accessible route toward polarization-programmable near-field nanoscopy, enabling chiral spectroscopy, selective excitation of spin/valley degrees of freedom, and quantum optical investigations at the nanoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Independent normalization for γ -ray strength functions: The shape method

Here, the shape method, a novel approach to obtain the functional form of the γ-ray strength function (γSF), is introduced. In connection with the Oslo method the slope of the nuclear level density (NLD) and γSF can be obtained simultaneously even in the absence of neutron resonance spacing data. The foundation of the shape method lies in the primary γ-ray transitions which preserve information on the functional form of the γSF. The shape method has been applied to 56 Fe, 92 Zr, and 164 Dy, which are representative cases for the variety of situations encountered in typical NLD and γSF studies. The comparisons of results from the shape method to those from the Oslo method demonstrate that the functional form of the γSF is retained regardless of nuclear structure details or J π values of the states fed by the primary transitions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗