Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Electromagnetic Simulation”

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

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

At least 379 records · Page 21

Physics-constrained machine learning for electrodynamics without gauge ambiguity based on Fourier transformed Maxwell’s equations

We utilize a Fourier transformation-based representation of Maxwell’s equations to develop physics-constrained neural networks for electrodynamics without gauge ambiguity, which we label the Fourier–Helmholtz–Maxwell neural operator method. In this approach, both of Gauss’s laws and Faraday’s law are built in as hard constraints, as well as the longitudinal component of Ampère–Maxwell in Fourier space, assuming the continuity equation. An encoder–decoder network acts as a solution operator for the transverse components of the Fourier transformed vector potential, $\hat{A}_⟂(k,t)$, whose two degrees of freedom are used to predict the electromagnetic fields. This method was tested on two electron beam simulations. Among the models investigated, it was found that a U-Net architecture exhibited the best performance as it trained quicker, was more accurate and generalized better than the other architectures examined. We demonstrate that our approach is useful for solving Maxwell’s equations for the electromagnetic fields generated by intense relativistic charged particle beams and that it generalizes well to unseen test data, while being orders of magnitude quicker than conventional simulations. We show that the model can be re-trained to make highly accurate predictions in as few as 20 epochs on a previously unseen data set.

97 MATHEMATICS AND COMPUTING↗

Gyrokinetic simulation of low-frequency Alfvénic modes in DIII-D tokamak

Global gyrokinetic simulations find that a beta-induced Alfvén eigenmode (BAE) and a Low frequency mode (LFM) coexist in the DIII-D tokamak experiments. The LFM can be excited without fast ions and has a frequency inside the gap of the beta-induced Alfvén-acoustic eigenmode (BAAE). The mode structure of the LFM from simulations and many of its parametric dependencies are consistent with experimental observations. For the BAE, the simulated mode structure is consistent with experiment but the frequency is lower than experiment. An antenna scan shows that the LFM is NOT the conventioanl BAAE; rather, it is an interchange-like electromagnetic mode excited by non-resonant drive of pressure gradients. The compressible magnetic perturbations significantly increase the growth rates of the BAE and LFM. On the other hand, effects of trapped electrons and equilibrium current have modest effects on the BAE and LFM. Finally, radial transfer of wave energy is different between the BAE and LFM, resulting in their respective triangular and elliptic mode shapes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modal Field Reconstruction in Resonant Cavities in the Fundamental and Undermoded Frequency Regimes

Theory, simulations, and experiments are presented that demonstrate reconstruction of electromagnetic fields in a cavity from sparse probe measurements. Such techniques are often referred to as virtual sensing, allowing fields at unobserved locations to be predicted. These methods are appropriate for the fundamental and undermoded regimes, providing the ability to estimate fields (and shielding effectiveness) throughout an arbitrarily shaped cavity from a few judiciously spaced probes. A modal simulation method is implemented that allows the response of arbitrarily shaped cavities to be rapidly computed with respect to varying probe locations and slot parameters, enabling statistical analysis of probe placement on reconstruction performance. A cylindrical vessel with numerous probe holes is developed for experiments, referred to as Perforated Vessel 2 (PV2). Experiments are performed on the vessel with and without a steel box inside, where transmit power is delivered into the vessel either through probes (probe injection) or through slots using an external antenna (slot excitation). Simulations and experiments illustrate that when the number of probes is minimal (equal to the number of mode coefficients to be estimated at each frequency), probe placement is critical to avoid missed peaks and to have acceptable reconstruction error. Probe placement becomes less important as the number of probes is increased, but care is still required to avoid probe locations giving poor performance.

42 ENGINEERING↗

Dissertation Proposal: Innovations in Pulsed Power and Plasma Science Theory Simulation and Experiment

Pulsed power and plasma physics are topics of great study at both Sandia National Laboratories (SNL or Sandia) and the University of New Mexico (UNM). The goal of this research is to further knowledge and understanding of these fields using the resources of both SNL and UNM in three ways. The first way is through the comprehension, application, and testing of theory. Reading and analytically deriving theoretical solutions of problems both real-world and simplified will allow for a fresh perspective and the furthering of the theory. One such theory is Ottinger's generalized theory for voltage measurement in magnetically insulated transmission lines (MITLs). By working through the math, a deeper understanding of the theory is gained from which one may add more physically accurate and/or more detailed physics into the theory. Additionally, understanding the theory lays a good foundation from which one can analyze, test, and compare results to the theory in the following two ways that will advance the fields of pulsed power and plasma physics. The second way is through the modeling and simulation of real-world and simplified problems that utilize and test the afore mentioned theories. Theory can be applied to a simulation domain by using the unstructured time-domain electromagnetic (UTDEM) codes EMPHASIS and EMPIRE as well as the physical modeling software CUBIT, all of which were developed at SNL. Problems such as the modeling and design of the extended MITL on HERMES III, the understanding of space-charge-limited emission from vacuum cathodes, and the interaction between a relativistic electron beam and an ideal gas can all be modeled, simulated, and analyzed with this set of codes. Here the advantage is three-fold. Firstly, theory that describes our understanding of these problems can be put to the test and advanced through iterative simulation and analysis. Secondly, the understanding of these problems will have a positive impact on national security through the advancement of the technological capability of the United States of America. Thirdly, and not unrelated to the prior advantage, is the validation and verification of EMPIRE and EMPHASIS. This segues into the third way, which is through experiment and the comparison of experiment to simulated and theoretical results. Performing experimental comparisons completes the scientific method and grounds all of the work in reality. Being able to physically test theory and simulation is necessary for any real conclusions to be drawn. Another advantage for carrying out experimental work is to advance the physical testing capabilities of SNL. Several systems will be developed and tested through the course of this work that positively impact technological advancement of Sandia National Labs. Lastly, all of the above work will converge to yield a well-rounded perspective that ties the three categories of research together.

42 ENGINEERING↗

Simulations of heat fluxes in an ELMy H-mode discharge on HL-2A

In order to study the distribution and evolution of the transient heat flux on HL-2A during edge-localized-mode (ELM) bursts, the BOUT++ electromagnetic six-field two-fluid model is used to simulate the pedestal collapse under the lower single-null divertor geometry. The equilibrium profiles of HL-2A ELMy H-mode discharge No. 24 953 are adopted as the initial condition in the original case. In this instance, linear analysis shows that the resistive ballooning mode (RBM) and drift-Alfven wave are unstable to this equilibrium, and RBM is the dominant instability. The evolutions of the radial heat fluxes at the outer mid-plane and heat fluxes to the inner and outer targets during the ELM event are presented. Six more equilibria are constructed based on the original case to find out the influence of the pedestal profiles on the peak electron heat flux. The results indicate that the heat flux increases with temperature and/or density, and the theoretical analysis and simulation results consistently show that the heat flux q ∥e is proportional to n e0,SEP $T^{\frac32}_{e0,SEP}$.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of the EDGES Low-band Antenna Beam Model

The response of the antenna is a source of uncertainty in measurements with the Experiment to Detect the Global Epoch of Reionization Signature (EDGES). We aim to validate the electromagnetic beam model of the low-band (50–100 MHz) dipole antenna with comparisons between models and against data. We find that simulations of a simplified model of the antenna over an infinite perfectly conducting ground plane are, with one exception, robust to changes in the numerical electromagnetic solver code or algorithm. For simulations of the antenna with the actual finite ground plane and realistic soil properties, we find that two out of three numerical solvers agree well. Applying our analysis pipeline to a simulated drift-scan observation from an early EDGES low-band instrument that had a 10 m × 10 m ground plane, we find residual levels after fitting and removing a five-term foreground model from the simulated data binned in local sidereal time (LST) average about 250 mK with ±40 mK variation between numerical solvers. A similar analysis of the primary 30 m × 30 m sawtooth ground plane reduced the LST-averaged residuals to about 90 mK with ±10 mK between the two viable solvers. More broadly we show that larger ground planes generally perform better than smaller ground planes. Simulated data have a power that is within 4% of real observations, a limitation of net accuracy of the sky and beam models. We observe that residual spectral structures after foreground model fits match qualitatively between simulated data and observations, suggesting that the frequency dependence of the beam is reasonably represented by the models. We find that a soil conductivity of 0.02 S m{sup -1} and relative permittivity of 3.5 yield good agreement between simulated spectra and observations. This is consistent with the soil properties reported by Sutinjo et al. for the Murchison Radio-astronomy Observatory, where EDGES is located.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Report on G4-Med, a Geant4 benchmarking system for medical physics applications developed by the Geant4 Medical Simulation Benchmarking Group

Geant4 is a Monte Carlo code extensively used in medical physics for a wide range of applications, such as dosimetry, micro- and nanodosimetry, imaging, radiation protection, and nuclear medicine. Geant4 is continuously evolving, so it is crucial to have a system that benchmarks this Monte Carlo code for medical physics against reference data and to perform regression testing. In this work, to respond to these needs, we developed G4-Med, a benchmarking and regression testing system of Geant4 for medical physics. G4-Med currently includes 18 tests. They range from the benchmarking of fundamental physics quantities to the testing of Monte Carlo simulation setups typical of medical physics applications. Both electromagnetic and hadronic physics processes and models within the prebuilt Geant4 physics lists are tested. The tests included in G4-Med are executed on the CERN computing infrastructure via the use of the geant-val web application, developed at CERN for Geant4 testing. The physical observables can be compared to reference data for benchmarking and to results of previous Geant4 versions for regression testing purposes. This paper describes the tests included in G4-Med and shows the results derived from the benchmarking of Geant4 10.5 against reference data.

60 APPLIED LIFE SCIENCES↗

Microwave-induced selective decomposition of cellulose: Computational and experimental mechanistic study

Understanding microwave-material interactions will help facilitate the utilization of microwave technology in gasification and renewable energy production. In this study, cellulose was used as a model compound to simulate the organic matter in biomass, and its catalytic decomposition under a microwave (MW) field was studied to identify structural changes from the reaction. The study was conducted using a MW source coupled to a fixed-bed gas-flow reactor, mass spectroscopic, and Fourier transform infrared spectroscopy post-reaction analysis. Zeolite 13X was chosen as a microwave absorber to study the catalytic enhancement of the decomposition of cellulose. Density functional theory (DFT) was used to gain insights into the molecular transformations occurring in the presence of a static electric field, which was used to simulate the electric field component of the microwave electromagnetic radiation. Theoretical calculations demonstrated that both the positive and negative portion of the electric field interact with the permanent dipoles of the cellulose effecting the decomposition through the glycosidic bond breaking mechanism. The theoretical result was verified using infrared spectroscopic analysis of the pure cellulose during microwave heating. The theoretical calculations suggest only the positive electric field component may be active for the cellulose decomposition in the presence of Zeolite 13X. Physically mixing Zeolite 13X with the cellulose led to a significant enhancement in the decomposition rate of the glycosidic bond. Zeolite 13X enhanced the glycosidic O-C decomposition at lower MW power (lower temperatures), whereas the O-H functional group required higher MW power (higher temperature) for its decomposition. The DFT study coupled with the reaction studies revealed that the electric field polarizability is dependent upon both the direction and the orientation of the cellulose. Gas products revealed that applying 250 W of MW power led to the production of CO, H 2 , along with some CO 2 , CH 4 , and benzene at 305 °C. Finally, reaction under 500, 750, and 1000 W of power at constant temperature (305 °C) revealed that higher power led to the complete decomposition of cellulose to mostly CO, H 2 , and CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Groundwater Circulation Within the Mountain Block: Combining Flow and Transport Models With Magnetotelluric Observations to Untangle Its Nested Nature

Abstract Mountains are vital water sources for humans and ecosystems, continuously replenishing lowland aquifers through surface runoff and mountain recharge. Quantifying these fluxes and their relative importance is essential for sustainable water resource management. However, our mechanistic understanding of the flow and transport processes determining the connection between the mountain block and the basin aquifer remains limited. Traditional conceptualizations assume groundwater circulation within the mountain block is predominantly shallow. This view neglects the role of deep groundwater flowpaths significantly contributing to the water, solute, and energy budgets. Overcoming these limitations requires a holistic characterization of the multiscale nature of groundwater flow along the mountain‐to‐valley continuum. As a proof‐of‐concept, we use a coupled groundwater flow and transport model to design a series of numerical experiments that explore the role of geology, topography, and weathering rates in groundwater circulation and their resulting resistivity patterns. Our results show that accumulating solutes near stagnation zones create contrasting electrical resistivity patterns that separate local, intermediate, and regional flow cells, presenting a target for magnetotelluric observations. To demonstrate the sensitivity of magnetotelluric data to features in our resistivity models, we use the MARE2DEM electromagnetic modeling code to perform forward and inverse simulations. This study highlights the potential of magnetotelluric surveys to image the resistivity structure resulting from multiscale groundwater circulation through relatively impervious crystalline basement rocks in mountainous terrains. This capability could change our understanding of the critical zone, offering a holistic perspective that includes deep groundwater circulation and its role in conveying solutes and energy.

54 ENVIRONMENTAL SCIENCES↗

Substation Cable Layouts for EMP Coupling Analysis

Direct coupling of early-time high-altitude electromagnetic pulse (HEMP) to substation control cables is simulated for cable layouts based on surveys of seven electrical substations in the United States. An analytic transmission line modeling code is used to estimate worst-case coupled current at the terminations of cable segments in or near the control shack. Where applicable, an induced voltage due to cable shield grounding is also estimated. Various configurations are simulated, including cables with different elevations, lengths, radii, and terminations. Plots of the coupled HEMP effects are given, and general relationships between these effects and the substations geometric and material parameters are highlighted and discussed.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Celeritas R&D Report: Accelerating Geant4

Celeritas is a new Monte Carlo (MC) detector simulation code designed for computationally intensive applications on high-performance heterogeneous architectures. In the past two years Celeritas has advanced from prototyping a Graphics Processing Unit (GPU)-based single physics model in infinite medium to implementing a full set of electromagnetic (EM) physics processes in complex geometries. The current release of Celeritas, version 0.4, has incorporated full device-based navigation, an event loop in the presence of magnetic fields, and detector hit scoring. New functionality incorporates a scheduler to offload electromagnetic physics to the GPU within a Geant4-driven simulation, enabling straightforward integration of Celeritas into the high energy physics (HEP) experimental frameworks CMSSW and ATLAS FullSimLight. On the Perlmutter supercomputer, Celeritas performs EM physics between 3× and 18× faster using the machine’s Nvidia GPUs compared to using only CPUs, corresponding to an electrical power efficiency up to a factor of 5. When running a multithreaded Geant4 ATLAS test beam application with full hadronic physics, using Celeritas to accelerate the EM physics results in an overall simulation speedup of 1.7–2.2× on GPU and 1.2× on CPU. In a CMS test application using tt¯ events and the prototype Run 4 configuration, compared to Geant4 CPU, Celeritas with a Nvidia A100 improves overall throughput up to a factor of 2.7× but cannot be efficiently shared with more than 8 cores.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Understanding the electromagnetic topology during the ohmic breakdown in tokamaks considering self-generated electric fields

We report that the physical mechanisms of the ohmic breakdown in a tokamak have been understood based on the classical Townsend avalanche theory. However, a new systematic theory (Yoo et al 2018 Nat. Commun.9 3523) recently demonstrated that electron avalanches during the ohmic breakdown are completely different from the Townsend avalanche due to strong self-generated electric fields. In this study, we elucidate the multi-dimensional effects of the self-generated electric field on plasma dynamics during the ohmic breakdown. We also propose a novel electromagnetic topology analysis method that can easily predict the overall plasma behavior and where the main plasma is generated. The topology analysis method is validated by a state-of-art particle simulation for various magnetic configurations. New physical insights into the complex electromagnetic topology would facilitate designing more reliable and optimized ohmic breakdown scenarios in future tokamaks, such as ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Robustness of neutron star merger simulations to changes in neutrino transport and neutrino-matter interactions

Binary neutron star mergers play an important role in nuclear astrophysics: Their gravitational wave and electromagnetic signals carry information about the equation of state of cold matter above nuclear saturation density, and they may be one of the main sources of r-process elements in the Universe. Neutrino-matter interactions during and after merger impact the properties of these electromagnetic signals and the relative abundances of the produced r-process elements. Existing merger simulations are, however, limited in their ability to realistically model neutrino transport and neutrino-matter interactions. Here, we perform a comparison of the impact of the use of state-of-the art two-moment or Monte Carlo transport schemes on the outcome of merger simulations, for a single binary neutron star system with a short-lived neutron star remnant [(5–10) ms]. We also investigate the use of different reaction rates in the simulations. While the best transport schemes generally agree well on the qualitative impact of neutrinos on the system, differences in the behavior of the high-density regions can significantly impact the collapse time and the properties of the hot tidal arms in this metastable merger remnant. The chosen interaction rates and transport algorithm as well as recent improvements by Radice et al. to the two-moment algorithms can all contribute to changes at the (10–30)% level in the global properties of the merger remnant and outflows. Furthermore, the limitations of previous moment schemes fixed by Radice et al. also appear sufficient to explain the large difference that we observed in the production of heavy-lepton neutrinos in a previous comparison of Monte Carlo and moment schemes in the context of a low-mass binary neutron star system.

79 ASTRONOMY AND ASTROPHYSICS↗

Suitability of Altair FEKO's Computational Electromagnetic Techniques for Tree Model Scattering and Penetration at Ku-band

Foliage penetration (FOPEN) radar at lower frequencies (VHF, UHF) is a well-studied area with many contributions. However, there is growing interest in using higher Ku-band frequencies (12-18 GHz) for FOPEN. Specifically, the reduced wavelength sizes provide some key saliencies for developing more optimized detection solutions. The disadvantage is that exploiting Ku-band for FOPEN is complicated because higher frequencies have pronounced scattering effects due to their smaller wavelengths. A methodology h as been developed to model and simulate FOPEN problems that characterize the phenomenology of Ku-band electromagnetic ( EM ) wave transmissions through moderate foliage. The details of this research (i.e. the realistic tree models, simulation setup and results) are documented in multiple reports. The main focus of this report is to describe the preliminary validation and verification of Altair FEKO, the computational EM (CEM) software used for this research, as well as present a simplified symmetrical tree model and an introductory CAD tree model.

42 ENGINEERING↗

Orbit-following simulations of fast-ion transport and losses due to the Alfvén eigenmode burst in the Large Helical Device

Orbit-following simulations of fast-ion transport and losses with time-dependent electromagnetic perturbations are performed to clarify the roles of Alfvén eigenmodes (AEs) and the low-frequency magnetohydrodynamic (MHD) mode observed in the kinetic-MHD hybrid simulation of AE bursts in the Large Helical Device. Fast-ion pressure profile flattening in the kinetic-MHD hybrid simulation can be reproduced by an orbit-following simulation with only the primary single AE of the time-dependent amplitude following the kinetic-MHD hybrid simulation result, while orbit-following simulations with constant AE amplitude of average level during AE burst cannot reproduce the fast-ion pressure profile flattening observed. The effects of other modes are negligible on the fast-ion pressure profile flattening. The fast-ion losses in kinetic-MHD hybrid simulation can be reproduced by an orbit-following simulation with time-dependent amplitude when the low-frequency MHD mode is considered in addition to multiple AEs. This indicates the synergetic effect of multiple AEs and the low-frequency MHD mode on fast-ion losses.

Seki, Ryohsuke↗

Novel Real- Time EMT- TS Modeling Architecture for Feeder Blackstart Simulations

This paper presents the development and benchmarking of a novel real-time electromagnetic-transient and transient-stability (EMT-TS) modeling architecture for distribution feeder restoration studies. The work presents for the first time in the literature a real-time EMT-TS testbed in which the grid-forming unit is simulated in EMT domain, operating as the slack bus of the phasor domain while including unbalanced voltage conditions. To evaluate the performance and limitations of the proposed model, an equivalent EMT testbed is developed to serve as a benchmark. First, the co-simulation framework and the domain coupling methodology are introduced. Next, the steady state operation of the EMT-TS and EMT models are matched. Then, tests are conducted to analyze the transient performance of the proposed EMT-TS model when compared to its EMT benchmark. Results reveal that when utilizing a battery energy storage system (BESS) as the grid-forming unit, the EMT-TS testbed can maintain high accuracy for typical load steps.

Paduani, Victor↗

Characterizing the Electromagnetic Response of the NOvA Detectors

The NOvA Test Beam Program was commissioned to complement the full-scale NOvA experiment in its goal of measuring neutrino oscillations. The NOvA Experiment consists of a 300-ton Near Detector and a 14-kton Far Detector, both exposed to the same beam, which can be operated in neutrino or antineutrino mode from Fermilab at an off-axis angle of 14.6~mrad. These detectors can discriminate between $\nu_\mu$-CC and $\nu_e$-CC neutrino (or corresponding antineutrino) events, enabling appearance and disappearance measurements of these flavors across the \SI{810}{\kilo\meter} near-to-far detector baseline. From these measurements, the underlying parameters of the three-flavor PMNS oscillations framework are measured. The NOvA Test Beam uses identical technology to the full-scale detectors, but scaled down to 30~tons. It is exposed to a beam of charged particles closely matched to the daughter products of NOvA (anti)neutrino interactions in both energy and species. An instrumented beamline determines both particle momentum and species upstream of the NOvA detector, allowing for a crisp measurement of the NOvA detector response to these known inputs. This dissertation presents an analysis of the electromagnetic response of the NOvA hardware and establishes a framework for incorporating these results to enhance the physics reach of the experiment. Prior to this work, the electromagnetic energy response of the NOvA detector relied primarily on simulation-driven calibration together with a globally applied energy scale uncertainty. The measurements presented here provide direct test beam constraints on the detector response, enabling particle-dependent and potentially energy-dependent uncertainty models derived from controlled electromagnetic inputs.

Myers, Dalton Grant [U. Texas, Austin (main)] (ORC↗