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 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Variational, stable, and self-consistent coupling of 3D electromagnetics to 1D transmission lines in the time domain

This work presents a new multiscale method for coupling the 3D Maxwell's equations to the 1D telegrapher's equations. While Maxwell's equations are appropriate for modeling complex electromagnetics in arbitrary-geometry domains, simulation cost for many applications (e.g. pulsed power) can be dramatically reduced by representing less complex transmission line regions of the domain with a 1D model. By assuming a transverse electromagnetic (TEM) ansatz for the solution in a transmission line region, we reduce the Maxwell's equations to the telegrapher's equations. Here, we propose a self-consistent finite element formulation of the fully coupled system that uses boundary integrals to couple between the 3D and 1D domains and supports arbitrary unstructured 3D meshes. Additionally, by using a Lagrange multiplier to enforce continuity at the coupling interface, we allow for an absorbing boundary condition to also be applied to non-TEM modes on this boundary. We demonstrate that this feature reduces non-physical reflection and ringing of non-TEM modes off of the coupling boundary. By employing implicit time integration, we ensure a stable coupling, and we introduce an efficient method for solving the resulting linear systems. We demonstrate the accuracy of the new method on two verification problems, a transient O-wave in a rectilinear prism and a steady-state problem in a coaxial geometry, and show the efficiency and weak scalability of our implementation on a cold test of the Z-machine MITL and post-hole convolute.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Microwave-Absorbing Foams with Adjustable Absorption Frequency and Structural Coloration

Microwave-absorbing materials with regulatable absorption frequency and optical camouflage hold great significance in intelligent electronic devices and advanced stealth technology. Herein, we present an innovative microwave-absorbing foam that can dynamically tune microwave absorption frequencies via a simple mechanical compression while in parallel enabling optical camouflage over broad spectral ranges by adjusting the structural colors. In this study, the vivid colors spanning different color categories generated from thin-film interference can be precisely regulated by adjusting the thickness of the conformal TiO 2 coatings on Ni/melamine foam. Enhanced interfacial and defect-induced polarizations resulting from the introduction of TiO 2 coating synergistically contribute to the dielectric attenuation performance. Consequently, such a foam exhibits exceptional microwave absorption capabilities, and the absorption frequency can be dynamically tuned from the S band to the Ku band by manipulating its compression ratio. Additionally, simulation calculations validate the adjustable electromagnetic wave loss behavior, offering valuable insights for the development of next-generation intelligent electromagnetic devices across diverse fields.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A curvature-based mechanism for the spontaneous emission of electromagnetic radiation at curved plasma–vacuum interfaces

In this work, we derive from first principles the dispersion relation for electromagnetic radiation that is spontaneously emitted at a 2D curved plasma–vacuum interface in the absence of both background and external electric and magnetic fields. This plasma geometry is motivated by plasmas generated through hypervelocity impact projectiles—an area that is of importance to satellite safety in the aerospace community, and may be of importance for the developing study of projectile impact-driven inertial confinement fusion. We find spontaneous radiation occurs when thermal surface waves propagate along the edge of a circular plasma, with a frequency lower than the plasma frequency and spatial decay proportional to 1/r. This is distinct from more well-known mechanisms of radiation, in which the radiation frequency is greater than or equal to the plasma frequency. This theory is validated against 2D, electromagnetic particle-in-cell simulations and shows reasonable agreement. Finally, we discuss the effects of plasma geometry and attribute the radiation generation mechanism to the acceleration of charged surface waves as they traverse azimuthally along a curved path.

Lau, Raymond (ORCID:0000000299318101)↗
DOE OSTI2024

Per-phase Phasor Modeling of GFL and GFM Inverters for Distribution System Dynamic Studies

A significant increase in inverter-based distributed energy resources (DERs) is expected in the near future. Existing electromagnetic models and sequential electromechanical models of DERs supporting dynamic analysis cannot be extended to largescale three-phase distribution systems usually with unbalanced construction. As such, there is a need for appropriate industry models that can support distribution system dynamic studies. This paper describes the latest phasor models of grid-following (GFL) and grid-forming (GFM) inverters that can capture distribution system dynamics under grid disturbances. The overall modeling performance is examined and compared with results from other platforms such as controller hardware-in-the-loop (CHIL), MATLAB/Simulink and PSCAD with electromagnetic transient (EMT) and/or phasor simulation capabilities under balanced and unbalanced faults.

Distribution system, Grid modelling, GridLAB-D, Po↗
DOE OSTI2022

Effects-Based Monitoring of Geomagnetically-Induced Current Using a Convolutional Neural Network

Geomagnetically-induced current (GIC) due to space weather can flow in the power grid causing undesirable effects such as transformer overheating, misoperation of protection devices, and potential blackouts. It is therefore important to monitor GIC in the power grid to improve online situational awareness and decision-making of system operators during a geomagnetic disturbance. To avoid the costly installation of GIC monitors at transformers’ neutrals, it is desirable to find correlations between GIC and already-monitored parameters. Hence, this work proposed the use of a convolutional neural network (CNN) to compute GIC amplitudes from learned patterns in the time-series data of transformer even harmonic currents. Using an electromagnetic transient program, GIC injection simulations were performed for a modeled Dominion Energy Virginia (DEV) substation with two 504 MVA, 500/230 kV transformers. Data collected from these offline simulations were used to train the CNN to provide online GIC monitoring. Testing the CNN performance involved using real GIC measurements from published literature and from a physical GIC monitor in the DEV area. Finally, the results showed that the proposed method was able to provide GIC readings with a root mean squared error of 1.56 A/phase (equivalent to an average accuracy of 94%) for these real GIC waveforms.

42 ENGINEERING↗
DOE OSTI2026

Nuclear recoil detection with color centers in bulk lithium fluoride

We present initial results on the detection of nuclear recoils in lithium fluoride (LiF) through the fluorescence of color centers created by particle interactions in the crystal lattice. Using light-sheet fluorescence microscopy, we image nuclear recoil tracks from both fast and thermal neutron interactions deep within a cubic-centimeter-scale sample. Automated three-dimensional feature extraction based on machine-learning tools enables the identification and classification of individual events. We observe that the fluorescence response of LiF to gamma irradiation is strongly suppressed, by a factor of 30–50 compared to neutron exposure, demonstrating intrinsic insensitivity to electromagnetic backgrounds. The observed and simulated event characteristics are consistent, including their number, size, and topology. These results establish the feasibility of LiF as a scalable detection medium for rare nuclear-recoil events and constitute a first step toward 10–1000 g scale detectors with single-event sensitivity for applications in reactor-neutrino detection, neutron spectroscopy, and dark matter searches.

Aroujo, G R [University of Zurich]↗
DOE OSTI2021

Redesign of the Jefferson Lab -300 kV DC Photo-Gun for High Bunch Charge Operations

Production of high bunch charge beams for the Electron-Ion Collider (EIC) is a challenging task. High bunch charge (a few nC) electron beam studies at Jefferson Lab using an inverted insulator DC high voltage photo-gun showed evidence of space charge limitations starting at 0.3 nC, limiting the maximum delivered bunch charge to 0.7 nC for beam at -225 kV, 75 ps (FWHM) pulse width, and 1.64 mm (rms) laser spot size. The low extracted charge is due to the modest longitudinal electric field (Ez) at the photocathode leading to beam loss at the anode and downstream beam pipe. To reach the few nC high bunch charge goal, and to correct the beam deflection exerted by the non-symmetric nature of the inverted insulator photo-gun the existing photo-gun was modified. This contribution discusses the electrostatic design of the modified photo-gun obtained using CST Studio Suite’s electromagnetic field solver. Beam dynamics simulations performed using General Particle Tracer (GPT) with the resulting electrostatic field map obtained from the modified electrodes confirmed the validity of the new design.

Wijethunga, S.A.K.↗
DOE OSTI2020

Cathode Modeling and Diode Performance in a Planar Geometry using Space-charge-limited Emission

Space-charge-limited (SCL) emission parameters are varied to study the performance effects in a planar diode using an electromagnetic particle-in-cell simulation software suite, EMPIRE. Oscillations in the simulations are found and linked to the emission parameters, namely the breakdown threshold, the emission delay time, and the current density ramp time. The oscillations are suggested to be a transverse oscillator due to the perfect magnetic conductor boundary condition in steady-state operation and the formation of a virtual cathode in the diode driven by the SCL boundary condition.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Rapid Demonstration of Bremsstrahlung Diode Optimization

Optimization of the radiation pattern from a Bremsstrahlung target for a given application is possible by controlling the electron beam that impacts the high-atomic-number target. In this work, the electron beam is generated by a 13MV vacuum diode that terminates a coaxial magnetically insulted transmission line (MITL) on the HERMES-III machine at Sandia National Labs. Work by Sanford introduced a geometry for vacuum diodes that can control the flow within bounds. The "indented anode", as coined by Sanford, can straighten out the electron beam in a high-current diode that would otherwise be prone to beam pinching. A straighter beam will produce a more forwardly directed radiation pattern while a pinching electron beam will yield a focal point or hot spot on axis and a more diffuse radiation pattern. Either one of these may be desirable depending on the application. This work serves as a first attempt to optimize the radiation pattern in the former sense of collimating the radiation pattern given a limited parameter space. The optimization is attempted first using electromagnetic particle-in-cell simulations in the EMPIRE code suite. The setup of the models used in EMPIRE is discussed along with some basic theory behind some of the models used in the simulations such as anode heating and secondary ions. Theoretical work performed by Allen Garner and his students at Purdue is included here, which concerns the impact of collisions in these vacuum diodes. The EMPIRE simulations consider both an aggressive and a conservative design. The aggressive design is inherently riskier while the conservative design is chosen as something that, while still a risk, is more likely to perform as expected. The ultimate goal of this work was to validate the EMPIRE code results with experimental data. While the experiment that tested the diode designs proposed by the simulation results fell outside of the fiscal boundaries of this project (and for that reason the results of which are not included in this report), the hardware for the experiment was designed and drafted within those same fiscal boundaries, and is thus included in this report. However, there was yet another experiment performed in this project that tested a key feature of the diode: the hemispherical cathode. Those results are documented here as well, which show that the cathode tip is an important aspect to controlling the diode flow. A short series of simulations on this diode were also performed after the experiment in order to gain a better understanding of the effect of ions. on the flow pattern and faceplate dose profile.

43 PARTICLE ACCELERATORS↗
DOE OSTI2022

Simulating the filamentation of smoothed laser beams with three-dimensional nonlinear dynamics

In a plasma, the ponderomotive force of an inhomogeneous electromagnetic field expels plasma from regions of high intensity. When a laser propagates through a plasma, this force creates density wells that subsequently modify the index of refraction. The beam refracts and focuses into these wells and may filament. In extreme cases, the laser beam will spray due to increasing angular divergence of the beam. The threshold for ponderomotive self-focusing is well established for isolated laser hotspots or speckles. Here, we define a practical threshold for characterizing the filamentation of thousands of speckles that are found in the focal plane of high-power laser beams spatially smoothed with random phase plates as used at high energy and power laser facilities studying inertial confinement fusion. This threshold is tested against three-dimensional simulations of speckled laser light propagating through plasma. Four metrics are applied to assess filamentation: the fraction of power above five times the average intensity, an effective f-number, the mean-squared perpendicular wavenumber, and the fraction of rarefied density with deviation from the initial density exceeding | δn/ n| = 0.1. The speckled beams studied are generated by random phase plates, both with and without additional polarization smoothing, in a parameter regime of relevance to indirect drive experiments. While filamentation has been discussed extensively in the literature, we believe this to be the first published simulation study with three-dimensional nonlinear hydrodynamics that addresses the onset threshold of ponderomotive filamentation and establishes the lengths and time scales necessary to reach a statistical steady state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Frequency Emitter Geolocation Using Signal Strength Fingerprinting Informed By 3D Propagation Modeling

This work focuses on the problem of RF geolocation in complex multipath environments. Using 3D electromagnetic propagation modeling to characterize environments of interest will enable more accurate RF geolocation. Specifically, a path-loss radio map can be generated in simulation for use in received signal strength indicators (RSSI) fingerprinting, or pattern matching. RSSI fingerprinting is an example of data-based method that takes site specific information into account which should allow for better performance than other model-based methods that use a generalized model of electromagnetic propagation. This modeling capability will also be used to evaluate the relative performance of RSSI fingerprinting, pathloss model based RSSI methods such as differential received signal strength circles (DRSS), RSSI joint gaussian estimation, and time-difference of arrival (TDOA). New methods using this simulation derived electromagnetic characterization could improve the efficacy of currently deployed and future RF spectral monitoring solutions. Wireless InSite developed by Remcom is used as the simulation tool of choice in this work. An indoor location is simulated with a grid of fixed receivers and a grid of transmit locations. Using the output of the Wireless InSite simulation the response from a given transmit location to a given receive location can be generated. During the first year of the project various geolocation methods evaluated on purely synthetic, but realistic, data. The second year focused on testing and validating the efficacy of simulation informed RF geolocation using two physical testbeds. This work has shown that data-based approaches are more accurate than model-based ones at the expensive of requiring measured or simulated site-specific training data.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

Characterization of the Redesigned Extended MITL for HERMES III

Modeling and simulation of the legacy HERMES III Magnetically Insulated Transmission Line (MITL) has been performed using EMPHASIS, an unstructured time-domain electromagnetic (UTDEM) particle-in-cell (PIC) simulation software. This design when used lost roughly half of its current before the pulse reached the load. The cause of the current loss in the MITL was found to be the vacuum impedance changes along the MITL. The MITL was then redesigned to maintain constant impedance and simulated in EMPHASIS once again. Following predicting simulation results, the new MITL was then built, installed, and tested, showing minimal current loss and good agreement with simulation and theoretical results, all of which are reported here. Additionally, an analysis of experimental voltage calculation techniques using cathode and anode currents is performed and compared to simulation results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Accelerating particle-in-cell kinetic plasma simulations via reduced-order modeling of space-charge dynamics using dynamic mode decomposition

We present a data-driven reduced-order modeling of the space-charge dynamics for electromagnetic particle-in-cell (EMPIC) plasma simulations based on dynamic mode decomposition (DMD). The dynamics of the charged particles in kinetic plasma simulations such as EMPIC is manifested through the plasma current density defined along the edges of the spatial mesh. We showcase the efficacy of DMD in modeling the time evolution of current density through a low-dimensional feature space. Not only do such DMD based predictive reduced-order models help accelerate EMPIC simulations, they also have the potential to facilitate investigative analysis and control applications. Here, we demonstrate the proposed DMD-EMPIC scheme for reduced-order modeling of current density and speedup in EMPIC simulations involving electron beam under the influence of magnetic field, virtual cathode oscillations, and backward wave oscillator.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Identifying the microtearing modes in the pedestal of DIII-D H-modes using gyrokinetic simulations

Recent evidence points toward the microtearing mode (MTM) as an important fluctuation in the H-mode pedestal for anomalous electron heat transport. A study of the instabilities in the pedestal region carried out using gyrokinetic simulations to model an ELMy H-mode DIII-D discharge (USN configuration, 1.4 MA plasma current, and 3 MW heating power) is presented. The simulations produce MTMs, identified by predominantly electromagnetic heat flux, small particle flux, and a substantial degree of tearing parity. The magnetic spectrogram from Mirnov coils exhibits three distinct frequency bands—two narrow bands at lower frequency (~35–55 kHz and ~70–105 kHz) and a broader band at higher frequency (~300–500 kHz). Global linear GENE simulations produce MTMs that are centered at the peak of the ω* profile and correspond closely with the bands in the spectrogram. The three distinctive frequency bands can be understood from the basic physical mechanisms underlying the instabilities. For example (i) instability of certain toroidal mode numbers (n) is controlled by the alignment of their rational surfaces with the peak in the ω* profile, and (ii) MTM instabilities in the lower n bands are the conventional collisional slab MTM, whereas the higher n band depends on curvature drive. While many features of the modes can be captured with the local approximation, a global treatment is necessary to quantitatively reproduce the detailed band gaps of the low-n fluctuations. Notably, the transport signatures of the MTM are consistent with careful edge modeling by SOLPS.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Calomplification — the power of generative calorimeter models

Motivated by the high computational costs of classical simulations, machine-learned generative models can be extremely useful in particle physics and elsewhere. They become especially attractive when surrogate models can efficiently learn the underlying distribution, such that a generated sample outperforms a training sample of limited size. This kind of GANplification has been observed for simple Gaussian models. We show the same effect for a physics simulation, specifically photon showers in an electromagnetic calorimeter.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

EMT Simulation of Large PV Plant & Power Grid for Disturbance Analysis

One of the challenges faced during the analysis of various disturbances with inverter-based resources is the ability to replicate the same disturbance in simulations. The main issue arises from lack of electromagnetic transient (EMT) models and adequate fidelity of models. To overcome this problem, in this paper, EMT simulation of power grid and one of the affected photovoltaic (PV) plants is performed to replicate the observations during the Angeles Forest disturbance in 2018. The processes to (i) develop the EMT model of the power grid from traditional transient stability (TS) data, and (ii) develop high-fidelity EMT model of PV plant from data collected on a PV plant are described in the paper. The requirements in the high-fidelity EMT model of the PV plant are also described. Finally, simulation of the Angeles Forest disturbance is performed and similar simulation results (as compared to the observed disturbances) are observed.

Debnath, Suman↗
DOE OSTI2023

Neutrino transport in general relativistic neutron star merger simulations

Abstract Numerical simulations of neutron star–neutron star and neutron star–black hole binaries play an important role in our ability to model gravitational-wave and electromagnetic signals powered by these systems. These simulations have to take into account a wide range of physical processes including general relativity, magnetohydrodynamics, and neutrino radiation transport. The latter is particularly important in order to understand the properties of the matter ejected by many mergers, the optical/infrared signals powered by nuclear reactions in the ejecta, and the contribution of that ejecta to astrophysical nucleosynthesis. However, accurate evolutions of the neutrino transport equations that include all relevant physical processes remain beyond our current reach. In this review, I will discuss the current state of neutrino modeling in general relativistic simulations of neutron star mergers and of their post-merger remnants. I will focus on the three main types of algorithms used in simulations so far: leakage, moments, and Monte-Carlo scheme. I will review the advantages and limitations of each scheme, as well as the various neutrino–matter interactions that should be included in simulations. We will see that the quality of the treatment of neutrinos in merger simulations has greatly increased over the last decade, but also that many potentially important interactions remain difficult to take into account in simulations (pair annihilation, oscillations, inelastic scattering).

Foucart, Francois (ORCID:0000000346174738)↗
DOE OSTI2020

Unitary quantum lattice simulations for Maxwell equations in vacuum and in dielectric media

Utilizing the similarity between the spinor representation of the Dirac and the Maxwell equations that has been recognized since the early days of relativistic quantum mechanics, a quantum lattice algorithm (QLA) representation of unitary collision-stream operators of Maxwell's equations is derived for both homogeneous and inhomogeneous media. A second-order accurate 4-spinor scheme is developed and tested successfully for two-dimensional (2-D) propagation of a Gaussian pulse in a uniform medium whereas for normal (1-D) incidence of an electromagnetic Gaussian wave packet onto a dielectric interface requires 8-component spinors because of the coupling between the two electromagnetic polarizations. In particular, the well-known phase change, field amplitudes and profile widths are recovered by the QLA asymptotic profiles without the imposition of electromagnetic boundary conditions at the interface. The QLA simulations yield the time-dependent electromagnetic fields as the wave packet enters and straddles the dielectric boundary. QLA involves unitary interleaved non-commuting collision and streaming operators that can be coded onto a quantum computer: the non-commutation being the very reason why one perturbatively recovers the Maxwell equations.

Physics↗