SEARCH · Search NASA

Results for “Numerical electromagnetics code”

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.

SourceResearch connectionFrom yearThrough year

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

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRS1994

Acceleration of electrons and ions by strong lower-hybrid turbulence in solar flares

One of the outstanding problems in solar flare theory is how to explain the 10-20 keV and greater hard x-ray emissions by a thick target bremsstrahlung model. The model requires the acceleration mechanism to accelerate approximately 10(exp 35) electrons sec(exp -l) with comparable energies, without producing a large return current which persists for long time scales after the beam ceases to exist due to Lenz's law, thereby, producing a self-magnetic field of order a few mega-Gauss. In this paper, we investigate particle acceleration resulting from the relaxation of unstable ion ring distributions, producing strong wave activity at the lower hybrid frequency. It is shown that strong lower hybrid wave turbulence collapses in configuration space producing density cavities containing intense electrostatic lower hybrid wave activity. The collapse of these intense nonlinear wave packets saturate by particle acceleration producing energetic electron and ion tails. There are several mechanisms whereby unstable ion distributions could be formed in the solar atmosphere, including reflection at perpendicular shocks, tearing modes, and loss cone depletion. Numerical simulations of ion ring relaxation processes, obtained using a 2 1/2-D fully electromagnetic, relativistic particle in cell code are discussed. We apply the results to the problem of explaining energetic particle production in solar flares. The results show the simultaneous acceleration of both electrons and ions to very high energies: electrons are accelerated to energies in the range 10-500 keV, while ions are accelerated to energies of the order of MeVs, giving rise to x-ray emission and gamma-ray emission respectively. Our simulations also show wave generation at the electron cyclotron frequency. We suggest that these waves are the solar millisecond radio spikes. The strong turbulence collapse process leads to a highly filamented plasma producing many localized regions for particle acceleration and resulting in approximately 10(exp 17) electron 'beamlets' of width approximately equal to 10 lambda sub De which eliminates the production of large magnetic fields. In this paper, we demonstrate that the model produces an energetic electron spectrum with the right flux to account for the hard x-ray observations.

Spicer, D. S.↗

NASA NTRS1991

FD-TD numerical simulation of an entire lightning strike on the C160 aircraft

Experimental transient electromagnetic field measurements were performed on a Transall C160 aircraft during in-flight lightning strikes. The data allow a test of the predictive capabilities of a three dimensional time domain finite difference code (ALICE) developed at ONERA in order to investigate lightning-aircraft interactions. Using a transfer function technique in the 3D code, it is shown that a bi-leader attached to an aircraft can be simulated by a linear model, and so the electromagnetic fields can be calculated anywhere on the vehicle. Comparison of experimental and numerical results were made for several lightning strikes. Skin current density and electromagnetic field distributions are discussed in detail.

Alliot, J. C.↗

NASA NTRS2010

The Final Merger of Massive Black Holes: Recoils, Gravitational Waves, and Electromagnetic Signatures

The final merger of two massive black holes produces a powerful burst of gravitational radiation, emitting more energy than all the stars in the observable universe combined. The resulting gravitational waveforms will be easily detectable by the space-based LISA out to redshifts z > 10, revealing the masses and spins ofthe black holes to high precision. If the merging black holes have unequal masses, or asymmetric spins, the final black hole that forms can recoil with a velocity exceeding 1000 km/s. And, when the black holes merge in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. For more than 30 years, scientists have tried to compute black hole mergers using the methods of numerical relativity. The resulting computer codes have been plagued by instabilities, causing them to crash well before the black holes in the binary could complete even a single orbit. Within the past few years, however, this situation has changed dramatically, with a series of remarkable breakthroughs. This talk will focus on new results that are revealing the dynamics and waveforms of binary black hole mergers, recoil velocities, and the possibility of accompanying electromagnetic outbursts.

Centrella, Joan↗

NASA NTRS2010

The Final Merger of Massive Black Holes: Recoils, Gravitational Waves, and Electromagnetic Signatures

The final merger of two massive black holes produces a powerful burst of gravitational radiation, emitting more energy than all the stars in the observable universe combined. The resulting gravitational waveforms will be easily detectable by the space-based LISA out to redshifts z greater than 10, revealing the masses and spins of the black holes to high precision. If the merging black holes have unequal masses, or asymmetric spins, the final black hole that forms can recoil with a velocity exceeding 1000 km/s. And, when the black holes merge in the presence of gas and magnetic fields, various types of electromagnetic signals may also be produced. For more than 30 years, scientists have tried to compute black hole mergers using the methods of numerical relativity. The resulting computer codes have been plagued by instabilities, causing them to crash well before the black holes in the binary could complete even a single orbit. Within the past few years, however, this situation has changed dramatically, with a series of remarkable breakthroughs. This talk will focus on new results that are revealing the dynamics and waveforms of binary black hole mergers, recoil velocities, and the possibility of accompanying electromagnetic outbursts.

Centrella, Joan M.↗

NASA NTRS1995

A Numerical Analysis of Electromagnetic Scattering from Two-Dimensional Edge Terminations

Several techniques that influence the low frequency scattering from penetrable edge terminations are evaluated using a hybrid finite element and boundary element method code. The edge terminations consist of a dielectric skin forming an exterior shape with an internal conducting bulkhead. Some of the techniques considered are bulkhead shaping, internal material loading, placement of resistive cards, and the placement of lossy dielectric material rods. The intent of the various treatments is to find a combination or combinations that influence(s) the backscattered field to acceptable levels over a range of frequencies for both transverse magnetic and transverse electric polarizations.

Gray, B. E.↗

NASA NTRS1998

Numerical Differentiation for Adaptively Refined Finite Element Meshes

Postprocessing of point-wise data is a fundamental process in many fields of research. Numerical differentiation is a key operation in computational electromagnetics. In the case of data obtained from a finite element method with automatic mesh refinement much work needs still to be done. This paper addresses some issues in differentiating data obtained from a finite element electromagnetic code with adaptive mesh refinement, and it proposes a methodology for deriving the electric field given the magnetic field on a mesh of linear triangular elements. The procedure itself is nevertheless more general and might be extended for numerically differentiating any point-wise solution based on triangular meshes.

Borgioli, Andrea↗

DOE OSTI2021

Higher-order particle representation for particle-in-cell simulations

In this paper we present an alternative approach to the representation of simulation particles for unstructured electrostatic and electromagnetic PIC simulations. In our modified PIC algorithm we represent particles as having a smooth shape function limited by some specified finite radius, r 0 . A unique feature of our approach is the representation of this shape by surrounding simulation particles with a set of virtual particles with delta shape, with fixed offsets and weights derived from Gaussian quadrature rules and the value of r 0 . As the virtual particles are purely computational, they provide the additional benefit of increasing the arithmetic intensity of traditionally memory bound particle kernels. The modified algorithm is implemented within Sandia National Laboratories' unstructured EMPIRE-PIC code, for electrostatic and electromagnetic simulations, using periodic boundary conditions. We show results for a representative set of benchmark problems, including electron orbit, a transverse electromagnetic wave propagating through a plasma, numerical heating, and a plasma slab expansion. In this work, good error reduction across all of the chosen problems is achieved as the particles are made progressively smoother, with the optimal particle radius appearing to be problem-dependent.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1987

Hypercube matrix computation task

The Hypercube Matrix Computation (Year 1986-1987) task investigated the applicability of a parallel computing architecture to the solution of large scale electromagnetic scattering problems. Two existing electromagnetic scattering codes were selected for conversion to the Mark III Hypercube concurrent computing environment. They were selected so that the underlying numerical algorithms utilized would be different thereby providing a more thorough evaluation of the appropriateness of the parallel environment for these types of problems. The first code was a frequency domain method of moments solution, NEC-2, developed at Lawrence Livermore National Laboratory. The second code was a time domain finite difference solution of Maxwell's equations to solve for the scattered fields. Once the codes were implemented on the hypercube and verified to obtain correct solutions by comparing the results with those from sequential runs, several measures were used to evaluate the performance of the two codes. First, a comparison was provided of the problem size possible on the hypercube with 128 megabytes of memory for a 32-node configuration with that available in a typical sequential user environment of 4 to 8 megabytes. Then, the performance of the codes was anlyzed for the computational speedup attained by the parallel architecture.

Calalo, R.↗

NASA NTRS1993

On the stability of nongyrotropic ion populations - A first (analytic and simulation) assessment

The wave and dispersion equations for perturbations propagating parallel to an ambient magnetic field in magnetoplasmas with nongyrotropic ion populations show, in general, the occurrence of coupling between the parallel (left- and right-hand circularly polarized electromagnetic and longitudinal electrostatic) eigenmodes of the associated gyrotropic medium. These interactions provide a means to driving linearly one mode with free-energy sources of other modes in homogeneous media. Different types of nongyrotropy bring about distinct classes of coupling. The stability of a hydrogen magnetoplasma with anisotropic, nongyrotropic protons that only couple the electromagnetic modes to each other is investigated analytically (via solution of the derived dispersion equation) and numerically (via simulation with a hybrid code). Nongyrotropy enhances growth and enlarges the unstable spectral range relative to the corresponding gyrotropic situation. The relevance of the properties of nongyrotropic populations to space plasma environments is also discussed.

Brinca, A. L.↗

DOE OSTI2025

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗

DOE OSTI2022

Global gyrokinetic simulations of ASDEX Upgrade up to the transport timescale with GENE–Tango

An accurate description of turbulence up to the transport timescale is essential for predicting core plasma profiles and enabling reliable calculations for designing advanced scenarios and future devices. Here, we exploit the gap separation between turbulence and transport timescales and couple the global gyrokinetic code GENE to the transport-solver Tango, including kinetic electrons, collisions, realistic geometries, toroidal rotation and electromagnetic effects for the first time. This approach overcomes gyrokinetic codes' limitations and enables high-fidelity profile calculations in experimentally relevant plasma conditions, significantly reducing the computational cost. We present numerical results of GENE–Tango for two ASDEX Upgrade discharges, one of which exhibits a pronounced peaking of the ion temperature profile not reproduced by TGLF–ASTRA. We show that GENE–Tango can correctly capture the ion temperature peaking observed in the experiment. By retaining different physical effects in the GENE simulations, e.g., collisions, toroidal rotation and electromagnetic effects, we show that the ion temperature profile's peaking can be linked to electromagnetic effects of submarginal (stable) KBM modes. Furthermore, based on these results, the expected GENE–Tango speedup for the ITER standard scenario is larger than two orders of magnitude compared to a single gyrokinetic simulation up to the transport timescale, possibly making first-principles ITER simulations feasible on current computing resources.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1991

Numerical studies of electron dynamics in oblique quasi-perpendicular collisionless shock waves

Linear and nonlinear electron damping of the whistler precursor wave train to low Mach number quasi-perpendicular oblique shocks is studied using a one-dimensional electromagnetic plasma simulation code with particle electrons and ions. In some parameter regimes, electrons are observed to trap along the magnetic field lines in the potential of the whistler precursor wave train. This trapping can lead to significant electron heating in front of the shock for low beta(e). Use of a 64-processor hypercube concurrent computer has enabled long runs using realistic mass ratios in the full particle in-cell code and thus simulate shock parameter regimes and phenomena not previously studied numerically.

Liewer, P. C.↗

DOE OSTI2021

Numerical and Analytical Modeling of the Effect of Cracks on the Self-Inductance of a COTS YJ-41003-TC Toroid

COTS inductors and transformers often contain partial cracks whose effect on inductance, a key performance parameter, have not been carefully studied. In this report, the effects of both partial and complete cracks on the self-inductance of a 100 turn square cross section COTS YJ-41003-TC toroid comprised of J Material was comprehensively investigated using both analytically derived closed form expressions and 3D computational techniques employing commercial codes. Both partial (half-penny) and complete (air gap) cracks of 10 and 25 μm were investigated. The crack is defined as the physical distance between two faces of the toroid's magnetic core, such that the surface normal of either face is along the Φ-direction, in alignment with the B-field. For the purposes of validation, two different approaches were incorporated for both the analytical and numerical models. The two analytical methods are comprised of a first principles approach based on the physics of electromagnetics, as well as linear circuit theory. The former directly utilizes the integral form of Maxwell's equations while the latter exploits the interchangeable relationship between electric and magnetic circuits. Validation within the computational scheme is realized through a code-to-code comparison between commercial solvers, COMSOL Multiphysics and CST, with the former employing the Finite Element Method (FEM) and the latter the Finite Difference Time Domain (FDTD) technique. Sound agreement between all four methods (ie., two analytical and two numerical) is observed, with results indicating that only a perturbation in self-inductance occurs for the half-penny cracks, while a substantial reduction takes place for the case of complete cracks. It is important to note that even though a static μ r is applied, representing the linear region of the BH curve (based on manufacturer specifications), the complete crack results still place a lower conservative bound on the inductance. This follows from the fact that even in the case of a half-penny crack, if the magnetic core portion of the crack approaches saturation, the crack begins to behave like an air gap, or complete crack. When an air gap is introduced into a magnetic core, a substantial reduction in inductance can occur due to the significant difference in permeabilities between the two mediums (ie., μ core >> μ air ). The once intact bulk magnetic core of the toroid essentially begins to behave like an air core.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Efficient numerical algorithm for multi-level ionization of high-atomic-number gases

An efficient numerical algorithm for laser driven multi-level ionization of high-atomic-number gases is proposed and implemented in an electromagnetic particle-in-cell code SPACE. The algorithm is based on analytical solutions to the system of differential equations describing ionization evolution. Using analytical solutions resolves the multiscale issue of ionization due to different characteristic time scales of ionization processes and the main code time step. Algorithm efficiency and memory requirements are significantly improved by using a locally reduced system of differential equations. The algorithm also assigns proper orbital quantum numbers and their projections to ionization states. The algorithm is verified and validated using experimental data.

Cheng, A. (ORCID:000000021945282X)↗

DOE OSTI2020

Designing stellarators using perpendicular permanent magnets

We have developed a fast method to design perpendicular permanent magnets for simplifying stellarator coils based on existing codes. Coil complexity is one of the main challenges for stellarators. To date, only electromagnetic coils have been used to generate 3D fields for stellarators. Permanent magnets provide an alternative way to produce the desired magnetic field for optimized stellarators. In this paper, we revisit the concept of representing surface current using magnetic dipoles and carry out numerical validations. A surface magnetization is proven to be equivalent to the surface current that can be linearly solved by existing coil design codes. An incremental multi-layer method has been developed to obtain a practical solution that is attainable with present permanent magnets. With this method, we can reproduce a half-Tesla NCSX configuration using specially designed neodymium magnets together with simple planar coils. It shows that stellarator coils could be substantially simplified by adopting permanent magnets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1997

AXAF Detector Backgrounds Produced By Cosmic Ray Protons

One of the science instruments on the Advanced X-ray Astrophysics Facility (AXAF), planned for launch in 1998 into a highly elliptical (10,000 km x 140,000 km) orbit, is a microchannel plate High Resolution Camera (HRC). This detector is designed to provide imaging and spectroscopic observations of x-rays emitted by stellar sources in the 0.1 to 10 keV energy range. Described here are analyses made to determine the expected time-dependent detector background from prompt and delayed (activation) radiation initiated by galactic cosmic-ray (GCR) proton interactions in the spacecraft and payload. Numerical simulations were made using the coupled set of Monte Carlo radiation transport codes, analysis software, and data bases shown. The major codes are HETC for nucleon-meson transport, EGS for simulating electromagnetic cascades, and MORSE for low-energy (less than 15 MeV) neutron transport. The simulation follows the transport history of photons in the energy range from - 100 GeV down to approx. 0.1 keV due to gamma-ray sources from neutral pion decay, high-energy (spallation) collisions, and low-energy neutron inelastic scattering and capture reactions. Also included is radioisotope production and the tracking of gamma-rays, electrons, and positrons from induced radioactivity.

Armstrong, T. W.↗

NASA NTRS1991

Advanced electromagnetic methods for aerospace vehicles

The Advanced Helicopter Electromagnetics is centered on issues that advance technology related to helicopter electromagnetics. Progress was made on three major topics: composite materials; precipitation static corona discharge; and antenna technology. In composite materials, the research has focused on the measurements of their electrical properties, and the modeling of material discontinuities and their effect on the radiation pattern of antennas mounted on or near material surfaces. The electrical properties were used to model antenna performance when mounted on composite materials. Since helicopter platforms include several antenna systems at VHF and UHF bands, measuring techniques are being explored that can be used to measure the properties at these bands. The effort on corona discharge and precipitation static was directed toward the development of a new two dimensional Voltage Finite Difference Time Domain computer program. Results indicate the feasibility of using potentials for simulating electromagnetic problems in the cases where potentials become primary sources. In antenna technology the focus was on Polarization Diverse Conformal Microstrip Antennas, Cavity Backed Slot Antennas, and Varactor Tuned Circular Patch Antennas. Numerical codes were developed for the analysis of two probe fed rectangular and circular microstrip patch antennas fed by resistive and reactive power divider networks.

Balanis, Constantine A.↗