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 235 records · Page 13

Quasielastic lepton scattering and back-to-back nucleons in the short-time approximation

Understanding quasielastic electron- and neutrino-scattering from nuclei has taken on new urgency with current and planned neutrino oscillation experiments, and with electron scattering experiments measuring specific final states, such as those involving nucleon pairs in “back-to-back” configurations. Accurate many-body methods are available for calculating the response of light (A ≤ 12) nuclei to electromagnetic and weak probes, but they are computationally intensive and only applicable to the inclusive response. In the present work we introduce a novel approach, based on realistic models of nuclear interactions and currents, to evaluate the short-time (high-energy) inclusive and exclusive response of nuclei. The approach accounts reliably for crucial two-nucleon dynamics, including correlations and currents, and provides information on back-to-back nucleons observed in electron and neutrino scattering experiments. We demonstrate that in the quasielastic regime and at moderate momentum transfers both initial- and final-state correlations, and two-nucleon currents are important for a quantitatively successful description of the inclusive response and final state nucleons. Finally, the approach can be extended to include relativistic—kinematical and dynamical—effects, at least approximately in the two-nucleon sector, and to describe the response in the resonance-excitation region.

6 ≤ A ≤ 19↗

Structural Analysis of a Serpentine Superconducting Magnet for the Interaction Region of the Electron Ion Collider

The electron ion collider (EIC) is under design to collide the high energy and highly polarized hadron with electron beams with luminosities up to 10 34 cm −2 s −1 . The superconducting magnet designs at the interaction region (IR) are quite challenging due to the close proximity of the hadron and electron beams. Several serpentine types of superconducting magnets have been designed due to the space restrictions at the IR. Here, to validate the designs in terms of the mechanical strength, structural integrity, as well as the allowed deformations, a detailed structural mechanical analysis for a 12-layer superconducting quadrupole serpentine magnet (Q1BpR) in the rear side hadron beams has been carried out using the finite element method (FEM) considering the pretension stress, shrinking due to cool down, and the electromagnetic force on the conductors. Anisotropy material properties and individual roving tension have been considered. The electro-magnetic (EM) simulations were performed using COMSOL, confirmed with RAT, the structural mechanical analysis were performed using COMSOL by considering the contact elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Criteria for Grid-Forming Inverters Transitioning Between Current Limiting Mode and Normal Operation

Here, this letter focuses on the fault recovery analysis of inverters using multi-loop droop based grid-forming control, taking into account different current limiting strategies. A criterion is proposed to determine when the inverter can exit the current limiting mode after fault clearance. The theoretical findings are verified through electromagnetic transient simulations using the direct (d)-axis priority-based current limiting method, circular current limiting method, and virtual impedance method.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Monitoring Offshore CO 2 Sequestration Using Marine CSEM Methods; Constraints Inferred from Field- and Laboratory-Based Gas Hydrate Studies

Offshore geological sequestration of CO 2 offers a viable approach for reducing greenhouse gas emissions into the atmosphere. Strategies include injection of CO 2 into the deep-ocean or ocean-floor sediments, whereby depending on pressure–temperature conditions, CO 2 can be trapped physically, gravitationally, or converted to CO 2 hydrate. Energy-driven research continues to also advance CO 2 -for-CH 4 replacement strategies in the gas hydrate stability zone (GHSZ), producing methane for natural gas needs while sequestering CO 2 . In all cases, safe storage of CO 2 requires reliable monitoring of the targeted CO 2 injection sites and the integrity of the repository over time, including possible leakage. Electromagnetic technologies used for oil and gas exploration, sensitive to electrical conductivity, have long been considered an optimal monitoring method, as CO 2 , similar to hydrocarbons, typically exhibits lower conductivity than the surrounding medium. We apply 3D controlled-source electromagnetic (CSEM) forward modeling code to simulate an evolving CO 2 reservoir in deep-ocean sediments, demonstrating sufficient sensitivity and resolution of CSEM data to detect reservoir changes even before sophisticated inversion of data. Laboratory measurements place further constraints on evaluating certain systems within the GHSZ; notably, CO 2 hydrate is measurably weaker than methane hydrate, and >1 order of magnitude more conductive, properties that may affect site selection, stability, and modeling considerations.

58 GEOSCIENCES↗

Quasi-Helmholtz decomposition, Gauss' laws and charge conservation for finite element particle-in-cell

Development of particle-in-cell (PIC) methods using finite element based methods (FEMs) have been a topic of renewed interest; this has largely been driven by (a) the ability of finite element methods to better model geometry, (b) better understanding of function spaces that are necessary to represent all Maxwell quantities, and (c) more recently, the fundamental rubrics that should be obeyed in space and time so as to satisfy Gauss' laws and the equation of continuity. In that vein, methods have been developed recently that satisfy these equations and are agnostic to time stepping methods. While this development is indeed a significant advance, it should be noted that implicit FEM transient solvers support an underlying null space that corresponds to a gradient of a scalar potential ∇Φ(r) (or t∇Φ(r) in the case of wave equation solvers). While explicit schemes do not suffer from this drawback, they are only conditionally stable with time step sizes that are mesh dependent and very small. Furthermore, the null space produces spurious charge that can corrupt the desired physics of a PIC simulation. The way to overcome this bottleneck, and indeed, satisfy all four Maxwell's equation is to use a quasi-Helmholtz formulation on a tessellation. In the re-formulation presented, we strictly satisfy the equation of continuity and Gauss' laws for both the electric and magnetic flux densities. Results illustrating the efficacy of this scheme will be demonstrated by analyzing non-neutral beams with space-charge and the adiabatic expansion of a neutral plasma with realistic parameters (Debye length and real mass ratios).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Deep Learning-Assisted Real-Time Forward Modeling of Electromagnetic Logging in Complex Formations

Higher dimensional (i.e., 2-D and 3-D) modeling is indispensable to correctly evaluate the responses of electromagnetic logging tools in complex formation environments. However, limited by the high computational cost of rigorous modeling such as finite difference method and the finite element method, the real-time applications in the well logging industry primarily rely on the 1-D forward solver, which would result in erroneous formation evaluation for complex scenarios. As a result, aiming at realizing fast modeling for electromagnetic logging tools in complex formations, this paper proposes a general framework assisted by deep neural networks (DNNs). The framework consists of three modules: earth model classification, parameter extraction, and surrogate construction. Separate DNNs are trained and tested for different modules. The accuracy and efficiency of the DNN assisted fast modeling are validated by several experiments. Here, this study finds that the fast modeling assisted by DNNs is able to calculate the tool responses and reconstruct the subsurface formations in real-time.

97 MATHEMATICS AND COMPUTING↗

A Particle-in-Cell Method for Plasmas with a Generalized Momentum Formulation, Part I: Model Formulation

Here, this paper formulates a new particle-in-cell method for the Vlasov–Maxwell system. Under the Lorenz gauge condition, Maxwell’s equations for the electromagnetic fields can be written as a collection of scalar and vector wave equations. The use of potentials for the fields motivates the adoption of a Hamiltonian formulation for particles that employs the generalized (conjugate) momentum. A notable advantage offered by the Hamiltonian formulation is the elimination of time derivatives in the Lorenz gauge formulation that are required by the standard Newton–Lorentz treatment of the particles. This allows the fields to retain the full time-accuracy guaranteed by the field solver. The resulting updates for particles require only knowledge of the fields and their spatial derivatives. An analytical method for constructing these spatial derivatives is presented that exploits the underlying integral solution used in the field solver for the wave equations. Moreover, these derivatives are demonstrated to converge at the same rate as the fields in both time and space. The Method of Lines Transpose field solver we consider in this work is globally first-order accurate in time and high-order accurate in space (e.g., fourth- and fifth-order) and belongs to a larger class of methods which are unconditionally stable, can address geometry, and leverage $\mathcal {O}(N)$ fast summation methods for efficiency. We demonstrate the method on several well-established benchmark problems on bounded domains, including a plasma sheath as well as a relativistic particle beam. The efficacy of the proposed formulation is established by comparing with a second-order accurate finite-difference time-domain method that employs a leapfrog time advance for particles and a charge conserving map suitable for bounded domains. The new method shows mesh-independent numerical heating properties even in cases where the plasma Debye length is smaller than the grid spacing. This is an important feature of the new method for problems defined on bounded domains, because it permits the use of coarser grids in space in the representation of the fields. Such a capability has significant implications for the simulation of plasmas in bounded domains with complex geometry, where the ratio between the largest and smallest cells can vary significantly. The use of high-order spatial approximations in the new method also means that fewer grid points are required in order to achieve a fixed accuracy. Our results also suggest that the new method can be used with fewer simulation particles per cell compared to the benchmark explicit method, which permits further computational savings.

97 MATHEMATICS AND COMPUTING↗

Design of High Efficient Mid-Wavelength Infrared Polarizer on ORMOCHALC Polymer

While an organically modified chalcogenide (ORMOCHALC) can be used to fabricate a polymeric mid-wavelength infrared (MWIR) polarizer with competitive extinction ratio compared to the commercial wire-grid polarizers, which are composed of fragile inorganic materials, there is still a knowledge gap regarding the systematic design process to obtain high transmission efficiency and extinction ratio. As such, a computational parameter study for design optimization is conducted with the geometric parameters of the bilayer grating ORMOCHALC polarizer. The computational study shows that the Fabry–Pérot cavity is the primary mechanism that determines the transmission behaviors and the extinction ratio. A bilayer grating design, guided by the parameter study, is realized through the thermal nanoimprint and metal deposition processes. The extinction ratios measured with the Fourier-transform infrared are 245, 304, and 351 at the wavelength of 3, 4, and 5 μm, respectively. Compared to the state-of-the-art of the polymeric MWIR linear polarizers, the extinction ratio is improved by 1.4 times, and the transmission efficiency is increased by 2.5 times. Theoretical analysis with the multiple-layer model based on the transfer matrix method predicts a matched transmission behavior with the experiment and a full-wave electromagnetic simulation.

42 ENGINEERING↗

Implementing Faraday effect measurement constraints into the Grad–Shafranov equilibrium fitting code EFIT

A new tool for the exploration and diagnosis of the internal magnetic field of plasmas in the DIII-D tokamak in the form of a constraint on the EFIT (Equilibrium Fitting) Grad–Shafranov code based on the Faraday-effect Radial Interferometer-Polarimeter (RIP) diagnostic is presented, including description, verification, and sample application. The physics underlying the diagnostic and its implementation into EFIT are discussed, and the results showing the verification of the model are given, and the model’s limitations are discussed. The influence of the diagnostic’s input on the resulting equilibrium parameters is characterized. Furthermore, the effect of electron density profile refinement is evaluated and found to be negligible. A sample application of the diagnostic is shown, indicating that the RIP constraint has similar effects on the equilibrium as motional Stark effect constraints do.

47 OTHER INSTRUMENTATION↗

Review of Safety and Exposure Limits of Electromagnetic Fields (EMF) in Wireless Electric Vehicle Charging (WEVC) Applications

This study reviews the exposure limits and safety of intermediate frequency (IF) electromagnetic field (EMF) emissions for wireless electric vehicle charging (WEVC) applications. A review of the electromagnetic field exposure limits identified in international guidelines are presented. An overview of the electromagnetic field shielding technologies is provided including recommended geometries, materials, and performances of the methods available in the literature. Available laboratory results of EMF emissions are summarized considering several wireless power transfer studies in different power levels. Possible EMF reduction techniques are discussed with shielding practices and ORNL case studies. Also, living object detection (LOD) and foreign object detection (FOD) methods are reviewed from a safety aspect.

Asa, Erdem↗

Cosmogenic Background Rejection at the ICARUS [Poster]

ICARUS is on the surface of Earth and is exposed to huge cosmic activity and is the focus for several physics analyses. In summary, the overburden is very effective in reducing the hadrons while fully eliminating the electromagnetic cosmic ray components. Exploiting the information from three sub detectors and benchmarking the method with data is in progress.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Comparison of Electromagnetic Transient and Phasor-based Simulation for the Stability of Grid-Forming-Inverter-based Microgrids

The increasing interest in studying feeder-level, inverter-based microgrids requires appropriate simulation methods that can effectively simulate the dynamics of distribution feeders that typically have thousands of nodes. To examine the accuracy of phasor-based simulation for the stability of feeder-level, inverter-based microgrids, a comparative study between electromagnetic transient (EMT) and phasor-based simulation was conducted. Case studies of a two-inverter small microgrid and a modified version of the IEEE 34-Node Test System that has seven grid-forming inverters show that phasor-based simulation agrees well with EMT simulation when the P-f droop gains of grid-forming inverter controllers are within normal ranges. However, under extremely high P-f droop gains, phasor-based simulation can result in over-optimistic stability results. This study provides insight into the applicability of phasor-based simulation for the stability of feeder-level, inverter-based microgrids.

Du, Wei↗

Sensitivity of geophysical techniques for monitoring secondary CO 2 storage plumes

For geologic carbon storage, the ability to detect secondary CO 2 plumes—defined as those CO 2 plumes accumulating outside the intended storage reservoir—is fundamental to preventing unexpected CO 2 migration into groundwater resources and for risk and liability management. Understanding the sensitivity of various geophysical methods to secondary plumes is crucial for designing cost-effective monitoring schemes. We use several modeling scenarios to demonstrate the process of assessing sensitivities and detection thresholds of three primary geophysical techniques—surface seismic, borehole-to-surface electromagnetic (EM), and surface and borehole gravity—for early detection of secondary CO 2 plumes in the post-injection phase. While seismic reflection methods are often considered in monitoring strategies to track the evolution of CO 2 plumes, they are also the most expensive. Due to cost considerations, especially for long-term post-injection monitoring, other techniques complement seismic monitoring when designing an adaptive monitoring network. Borehole-to-surface EM or surface gravity surveys are feasible for time-lapse monitoring of deep secondary CO 2 plumes. Furthermore, these surveys could be carried at intervals defined by site-specific conditions. If time-lapse EM and/or gravity surveys detect any signal responses beyond the expected change, it would trigger a need for the higher resolution seismic survey.

58 GEOSCIENCES↗

Prototype tests of the electromagnetic particle injector-2 for fast time response disruption mitigation in tokamaks

Predicting and controlling disruptions is an important and urgent issue for ITER. Some disruptions with a short warning time may be unavoidable. For these cases, a fast time response disruption mitigation method is essential. Experimental tests on a prototype system of a novel, rapid time-response disruption mitigation system being developed for tokamak-based reactors and as a backup option for ITER, referred to as the electromagnetic particle injector (EPI), have been able to verify the primary advantages of the concept. These are its ability to meet short warning time scales of <10 ms while attaining the projected high velocities for deep radiative payload penetration in reactor-scale plasmas. The EPI relies on an electromagnetic propulsion system. A metallic sabot is accelerated electromagnetically to the required velocities (>1 km s -1 ) within 2 ms, at which point it releases a radiative payload consisting of a shell pellet or well-defined microspheres. Initial experimental tests from the prototype system show attainment of over 600 m s -1 in about 1 ms. Finally, essential aspects of payload separation from the sabot and sabot capture have also been demonstrated at 200 m s -1 , and the method can be extended to over 2 km s -1 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On numerical errors to the fields surrounding a relativistically moving particle in PIC codes

The particle-in-cell (PIC) method is widely used to model the self-consistent interaction between discrete particles and electromagnetic fields. It has been successfully applied to problems across plasma physics including plasma based acceleration, inertial confinement fusion, magnetically confined fusion, space physics, astrophysics, high energy density plasmas. In many cases the physics involves how relativistic particles (those with high relativistic γ factors) are generated and interact with plasmas. However, when relativistic particles stream across the grid, both in “vacuum” and in plasma, many numerical issues may arise which can lead to unphysical results. We present a detailed analysis of how discretized Maxwell solvers used in PIC codes can lead to numerical errors to the fields that surround particles that move at relativistic speeds across the grid. Expressions for the axial electric field as integrals in k space are presented that reveal two types of errors. The first arises from errors to the numerator of the integrand and leads to unphysical fields that are antisymmetric about the particle. Furthermore, the second arises from errors to the denominator of the integrand and lead to Cherenkov like radiation in “vacuum”. These fields are not anti-symmetric, extend behind the particle, and cause the particle to accelerate or decelerate depending on the solver and parameters. The unphysical fields are studied in detail for two representative solvers - the Yee solver and the FFT based solver. Although the Cherenkov fields are absent, the space charge fields are still present in the fundamental Brillouin zone for the FFT based solvers. In addition, the Cherenkov fields are present in higher order zones for the FFT based solvers. Comparison between the analytical solutions and PIC simulation results are presented. A solution for eliminating these unphysical fields by modifying the k operator in the axial direction is also presented. Using a customized finite difference solver, this solution was successfully implemented into OSIRIS [1]. Results from the customized solver are also presented. Additionally, this solution will be useful for a beam of particles that all move in one direction with a small angular divergence.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Classification of Intensity Distributions of Transmission Eigenchannels of Disordered Nanophotonic Structures Using Machine Learning

Light-matter interaction optimization in complex nanophotonic structures is a critical step towards the tailored performance of photonic devices. The increasing complexity of such systems requires new optimization strategies beyond intuitive methods. For example, in disordered photonic structures, the spatial distribution of energy densities has large random fluctuations due to the interference of multiply scattered electromagnetic waves, even though the statistically averaged spatial profiles of the transmission eigenchannels are universal. Classification of these eigenchannels for a single configuration based on visualization of intensity distributions is difficult. However, successful classification could provide vital information about disordered nanophotonic structures. Emerging methods in machine learning have enabled new investigations into optimized photonic structures. In this work, we combine intensity distributions of the transmission eigenchannels and the transmitted speckle-like intensity patterns to classify the eigenchannels of a single configuration of disordered photonic structures using machine learning techniques. Specifically, we leverage supervised learning methods, such as decision trees and fully connected neural networks, to achieve classification of these transmission eigenchannels based on their intensity distributions with an accuracy greater than 99%, even with a dataset including photonic devices of various disorder strengths. Simultaneous classification of the transmission eigenchannels and the relative disorder strength of the nanophotonic structure is also possible. Our results open new directions for machine learning assisted speckle-based metrology and demonstrate a novel approach to classifying nanophotonic structures based on their electromagnetic field distributions. These insights can be of paramount importance for optimizing light-matter interactions at the nanoscale.

79 ASTRONOMY AND ASTROPHYSICS↗

Manufactured solutions for an electromagnetic slot model

The accurate modeling of electromagnetic penetration is an important topic in computational electromagnetics. Electromagnetic penetration occurs through intentional or inadvertent openings in an otherwise closed electromagnetic scatterer, which prevent the contents from being fully shielded from external fields. To efficiently model electromagnetic penetration, aperture or slot models can be used with surface integral equations to solve Maxwell's equations. A necessary step towards establishing the credibility of these models is to assess the correctness of the implementation of the underlying numerical methods through code verification. Surface integral equations and slot models yield multiple interacting sources of numerical error and other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the electric-field integral equation with a slot model. Finally, we demonstrate the effectiveness of these approaches for a variety of cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A review of high order strong stability preserving two-derivative explicit, implicit, and IMEX methods

High order strong stability preserving time discretizations ensure the nonlinear non-inner-product strong stability properties of spatial discretizations suited for the stable simulation of hyperbolic PDEs in a wide variety of application areas including fluid dynamics, magnetohydrodynamics, semiconductor devices, electromagnetics, and astrophysics. Over the past decade multiderivative time-stepping have been increasingly used for the time-evolution hyperbolic PDEs, so that the strong stability properties of these methods have become important. In this work we review sufficient conditions for a two-derivative multistage method to preserve the strong stability properties of spatial discretizations in a forward Euler and different conditions on the second derivative. In particular we present the strong stability preserving theory for explicit and implicit two-derivative Runge–Kutta schemes, including a special condition on the second derivative under which these implicit methods may be unconditionally strong stability preserving. This special condition is natural for the stiff component of wide range of plasma physics problems, and can be useful in the context of strong stability preserving implicit-explicit multi-derivative Runge–Kutta schemes, where the time-step restriction is then independent of the stiff term. Lastly, we present the strong stability preserving theory for implicit-explicit multi-derivative general linear methods, and some novel second and third order methods where the time-step restriction is independent of the stiff term.

97 MATHEMATICS AND COMPUTING↗