Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “PLASMA 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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

The dynamics of low-beta plasma clouds as simulated by a three-dimensional, electromagnetic particle code

A 3D electromagnetic and relativistic particle simulation code is used to investigate the dynamics of low-beta plasma clouds moving perpendicular to an ambient magnetic field in vacuum and in a background plasma. The simulations show the formation of the space charge sheaths at the sides of the cloud with the associated polarization electric field which facilitate the cross-field propagation, as well as the sheaths at the front and rear end of the cloud caused by the larger ion Larmor radius, which allows ions to move ahead and lag behind the electrons as they gyrate. It is found that, in a background plasma, electron and ion sheaths expand along the magnetic field at the same rate, whereas in vacuum the electron sheath expands much faster than the ion sheath. Sheath clouds are accelerated up to relativistic energies, which indicates that artificial plasma clouds released in the ionosphere or magnetosphere may generate optical emissions as energetic sheath electrons scatter in the upper atmosphere.

Neubert, T.↗
DOE OSTI2013

Multi-moment advection scheme in three dimension for Vlasov simulations of magnetized plasma

We present an extension of the multi-moment advection scheme [T. Minoshima, Y. Matsumoto, T. Amano, Multi-moment advection scheme for Vlasov simulations, Journal of Computational Physics 230 (2011) 6800–6823] to the three-dimensional case, for full electromagnetic Vlasov simulations of magnetized plasma. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, and advances them on the basis of their governing equations. Similar to the two-dimensional scheme, the three-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile with little numerical dispersion or diffusion. This is a very important property for Vlasov simulations of magnetized plasma. We apply the scheme to electromagnetic Vlasov simulations. Propagation of linear waves and nonlinear evolution of the electron temperature anisotropy instability are successfully simulated with a good accuracy of the energy conservation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Spatial decay of electromagnetic waves from hypervelocity impact plasmas

Due to the high kinetic energy of impactors in space, hypervelocity impacts on spacecraft can ablate spacecraft material and produce plasmas. These plasmas have been observed in ground-based experiments to produce electromagnetic pulses (EMPs) and are suspected to be a cause of spacecraft anomalies. The initial EMP that is produced after plasma formation is replicated using a 2D electromagnetic plasma simulation, consisting of a discontinuous Galerkin, particle-in-cell (DG-PIC) framework. The spontaneous formation of EMPs from an initial quasineutral plasma is consistent with previous works and show behavior similar to that of waves emitted by antennas. The EMP emits from the plasma with a frequency less than the plasma frequency associated with the point of maximum density, showing the importance of the density gradient, and the far-field region of the “plasma antenna” is estimated. From the proportionality between plasma frequency and EMP frequency, there is a direct correlation between plasma density and radiation. Finally, as the plasma density decreases, the EMP frequency decreases, which suggests that several observable phenomena, such as optical flashes, radio frequency emission, and microwave emission, all result from the same plasma-related mechanism within a single hypervelocity impact.

42 ENGINEERING↗
DOE OSTI2014

Radio imaging of a type IVM radio burst on the 14th of August 2010

Propagating coronal mass ejections (CMEs) are often accompanied by burst signatures in radio spectrogram data. We present Nançay Radioheliograph observations of a moving source of broadband radio emission, commonly referred to as a type IV radio burst (type IVM), which occurred in association with a CME on the 14th of August 2010. The event was well observed at extreme ultraviolet (EUV) wavelengths by SDO/AIA and PROBA2/SWAP, and by the STEREO SECCHI and SOHO LASCO white light (WL) coronagraphs. The EUV and WL observations show the type IVM source to be cospatial with the CME core. The observed spectra is well fitted by a power law with a negative slope, which is consistent with optically thin gyrosynchrotron emission. The spectrum shows no turn over at the lowest Nançay frequencies. By comparing simulated gyrosynchrotron spectra with Nançay Radioheliograph observations, and performing a rigorous parameter search we are able to constrain several key parameters of the underlying plasma. Simulated spectra found to fit the data suggest a nonthermal electron distribution with a low energy cutoff of several tens to 100 keV, with a nonthermal electron density in the range 10{sup 0}-10{sup 2} cm{sup –3}, in a magnetic field of a few Gauss. The nonthermal energy content of the source is found to contain 0.001%-0.1% of the sources thermal energy content. Furthermore, the energy loss timescale for this distribution equates to several hours, suggesting that the electrons could be accelerated during the CME initiation or early propagation phase and become trapped in the magnetic structure of the CME core without the need to be replenished.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Numerical simulation of a plasma actuator based on ion transport

Two-dimensional numerical simulation of ion transport and flow around a single dielectric barrier discharge plasma actuator (PA) is performed. Spatial distributions of ions and electrons as well as their time evolution are obtained by solving the transport equations of monovalent positive ions, monovalent negative ions, and electrons. Voltage and frequency of the driving alternating-current signal are assumed to be 8 kV and 5 kHz, respectively. Special focus is laid upon the effect of voltage gradient dV/dt on the magnitude of the body force. The validity of steady force models often used in flow simulation is also examined. The simulation results show that the magnitude of the body force induced by the PA increases as the voltage gradient dV/dt increases and its increase rate becomes milder at higher voltage. The mechanism of body force generation is explained from the time evolution of number density fields of ions and electrons. A comparison between flow simulations using a time-resolved body force and its time-averaged counterpart demonstrates that the time-averaged model gives sufficiently accurate results when the time scale of the flow is more than 30 times greater than that of the PA.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

The generation of random variates from a relativistic Maxwellian distribution

A procedure for generating random variates from a relativistic Maxwellian distribution with arbitrary temperature and drift velocity is presented. The algorithm is based on the rejection method and can be used to initialize particle velocities in kinetic simulations of plasmas and gases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Development of a LIBS system for in situ surface measurements during plasma exposure in PISCES-A

A laser-induced breakdown spectroscopy (LIBS) system has been developed with a Q-switched Nd:YAG laser (wavelength = 1064 nm and pulse width ~5 ns) to conduct in situ surface measurements during plasma exposure in the PISCES-A linear divertor plasma simulator. The LIBS signal enhancement is obtained with both the magnetic field normal to the surface of a target and steady-state background plasma. Migration of sputtered Ta impurities onto the neighboring Cr surface is identified during He plasma exposure, only when cone structures are formed on the Cr surface. D retention in W during D plasma exposure is observed to decrease with increasing the sample temperature. The temporal evolution of D outgassing from W is measured in a time range of ~10-420 s right after D plasma exposure. A power law fit, t - α , yields α ~ 0.34 ± 0.09, which is nearly consistent with calculated and measured values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Experimental Measurements of Ion Heating in Collisional Plasma Shocks and Interpenetrating Supersonic Plasma Flows

Here, we present time-resolved measurements of ion heating due to collisional plasma shocks and interpenetrating supersonic plasma flows, which are formed by the oblique merging of two coaxial-gun-formed plasma jets. Our study is repeated using four jet species: N, Ar, Kr, and Xe. In conditions with small interpenetration between jets, the observed peak ion temperature T i is consistent with the predictions of collisional plasma-shock theory showing a substantial elevation of T i above the electron temperature T e and also the subsequent decrease of T i on the classical ion-electron temperature-equilibration timescale. In conditions of significant interpenetration between jets, such that shocks do not apparently form, the observed peak T i is still appreciable and greater than T e but much lower than that predicted by collisional plasma-shock theory. Experimental results are compared with multifluid plasma simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Path Integral Monte Carlo Simulations of Iron Plasmas (Final Technical Report)

This documents is the final technical report for our grant entitled "Path Integral Monte Carlo Simulations of Iron Plasmas" that focused on developing path integral Monte Carlo (PIMC) computer simulations. This techniques will be developed to study plasmas composed of heavier elements including iron and other third row elements. Equations of state (EOS) and transport properties will be derived in the regime of warm dense matter (WDM) and dense plasmas where existing first-principles methods cannot be applied. While standard density functional theory (DFT) has been used to accurately predict the structure of many solids and liquids up to temperatures on the order of 100,000 K, this method is not applicable at much higher temperature because the number of partially occupied electronic orbitals reaches intractably large numbers or the use of finite-temperature free energy functionals in orbital-free DFT introduces an uncontrolled approximation. Here we focus on PIMC methods that become more and more efficient with increasing temperatures and still include all electronic correlation effects. In this approach, electronic excitations increase the efficiency rather than reduce it. While it had commonly been assumed this method could only be applied to elements without core electrons, we showed that PIMC with free-particle nodes works well for first-row elements (PRL 108 (2012) 115502). Most recently, we extended the applicability range of all-electron PIMC to second-row elements by adopting localized nodal surfaces (PRL 115 (2015) 176403). To simulate third-row elements efficiently under WDM conditions, we propose a new method to remove core electrons by introducing pseudo-nodes. We explain our approach step by step and present preliminary results. We focus our method development on getting PIMC simulations of iron to work because of its fundamental importance for WDM and astrophysics. Then we move on to krypton and copper-doped beryllium, a ICF ablator material. We plan to continue working on key second-row material such as Na, Mg, MgO, Al, silica, and silicon-doped plastic ablators. Our collaborators at LLNL, will use our PIMC EOS data both as comparisons to existing semi-empirical, EOS-generating schemes, and as input for continuum radiation hydrodynamics simulations. We will establish an efficient pipeline from PIMC to macroscopic continuum studies of materials response. An emphasis will be placed on benchmarking such methods for plasmas of heavy elements at the very high temperatures (~100 eV) and low densities that are generated when Hohlraum radiation heats the ablator material in indirect drive laser experiments. Results from changes to the EOS will be of immeasurable importance to the designers at the National Ignition Facility (NIF) and at other facilities. Starting with our EOS of Cu-doped Be, our second collaborator at LLE, will perform real-time simulations of laser fusion experiments at the Omega laser and at the NIF to determine how sensitive the compression path depends on the ablator EOS. Since our collaborator also has experience in performing orbital-free DFT calculations, we propose to compare predictions from this method with PIMC results. In joint publications, we plan to analyze the accuracy of different free-energy functionals in order to understand why existing orbital-free DFT calculations do not predict compression peaks along the shock Hugoniot curve that we see with PIMC. The peaks are caused by the ionization of various electron shells. Their accurate characterization is important to compare with experimental results. We will break new ground by developing PIMC techniques that can simulate iron and all other third row elements in the plasma and WDM regimes. We introduce the concept of pseudo-nodes for the efficient treatment of the core-electrons. The EOS and transport properties will be derived and published online in the form of a new WDM database. Our PIMC EOS calculations will benchmark and possibly replace semi-analytical EOS tables like QEOS or SESAME, which will impact the hydrocode simulation community and will affect the design of NIF targets. Our PIMC results will help to improve the accuracy of orbital-free DFT simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Injection of an overdense electron beam in space

A three-dimensional particle simulation model is used to study the injection and propagation of an overdense electron beam in the vicinity of a conductor in space. Beam electrons with a density of more than 100 times the ambient electron density are modeled using large-scale plasma simulations; in these simulations the surface area of the conductor is several thousands times that of the beam cross section at the injection point. The parameters of the simulations are chosen to allow the realistic simulation of active space experiments such as the Space Shuttle's Spacelab 2 electron beam mission. These simulations confirm space observations that an overdense electron beam can at least partially escape the near vicinity of the spacecraft, even in a fully ionized plasma. Once they have escaped from the vicinity of the spacecraft, these beam electrons should be able to propagate away freely until their energy is dissipated. It is suggested that such large-scale simulations could be useful in interpreting data from space experiments.

Okuda, Hideo↗
DOE OSTI2017

Multiscale Modeling of Deformation, Fracture and Failure of Fusion Materials and Structures Final Report

The main objective of the research performed at UCLA under grant number DE-FG02-03ER54708 was to develop a physical understanding of the underlying mechanisms that control high-temperature deformation and fracture of fusion structural materials. The key aspect is the synergistic interaction between thermo-mechanical damage phenomena (e.g. creep, fatigue, and fracture), and irradiation damage phenomena (helium and hydrogen bubbles, and irradiation-induced defects). Another thrust area of the research was on the development of constitutive equations for the mechanical behavior of irradiated structures in fusion components (mainly the First Wall and Plasma-Facing components), and in utilizing them in descriptions in large-scale Finite Element Modeling of fusion structures. A parallel experimental effort on thermo-mechanical and plasma-induced damage is designed to validate and confront the developed theoretical and computational predictions. The models have been experimentally validated through experimental tests in simulated plasma transients on specially-designed W samples. This grant started at UCLA on November 14, 2003 and has had successive continuations till November 14, 2017. Research performed under the grant was to develop a range of hierarchical models for the post-elastic deformation, fracture and failure of fusion structural materials. The developed multiscale modeling approach has been based on rigorous mathematical, physical, and computational methods at the forefront of computational materials science. At the fundamental and smallest length scale (nm-μm), we developed advanced rate theory and Monte Carlo approaches to model microstructure evolution, non-equilibrium phase transformations, and dislocation-defect interactions. Microscopic and mesoscopic models of radiation hardening, ductile-to-brittle-transition, post-elastic deformation, plastic instabilities, and fracture processes are based on Dislocation Dynamics (DD) and Grain Growth Dynamics (GGD). At the continuum level, we focused on modeling fracture and failure mechanisms for Virtual Integrated Testing (VISTA), and for understanding the limits of both materials-by-design versus structural component design through advancing the new approach of multiscale-multiphysics. The overall research has been strongly coupled with the U.S. national experimental program to ascertain and verify the range of investigated phenomena. In the following, we give a summary of research progress supported by this grant since its inception till November 14, 2017. Detailed research results are available in our publications, listed in the publications section of this report.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Informing ARC divertor design and plasma facing material selection through integrated modeling

This INFUSE 2023 project between UCLA and Commonwealth Fusion Systems used computer modeling to test whether tungsten materials can survive in CFS's ARC fusion reactor divertor. The team simulated plasma conditions and material responses, finding that tungsten-rhenium alloys resist grain growth better than pure tungsten, and that hydrogen buildup depends more on particle flux than temperature. The work helps CFS design durable plasma-facing components for their fusion power plant.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Verification of 5D continuum gyrokinetic code COGENT: Studies of kinetic drift wave instability

COGENT (Continuum Gyrokinetic Edge New Technology) is a kinetic plasma simulation code that is being developed by the edge simulation laboratory (ESL) collaboration. The original version of the code has been developed in a 4D phase space (2D configuration space and 2D velocity space) to address kinetic plasma phenomena in a complex magnetic field geometry including the core, magnetic separatrix, and scrape-off layer regions. This paper is focused on extending the original 4D phase space to a 5D phase space (3D2V) to address full kinetic turbulences in a tokamak edge region. Here, we report on the current status of 5D COGENT, which presently operates in a shear-less, simple slab geometry. As a verification study, we use the problems of collision-less drift wave instability (universal instability) including its modification in the presence of collisional effects. The simulation model includes the gyrokinetic equations for ion and electron species coupled to the long-wavelength limit of the gyro-Poisson equation. Collisional effects are represented by the Krook collision model. Linear analytical results for the drift mode growth rate and real frequency are recovered, and the non-linear stage is modelled and analysed as well. In addition, extensive 5D runs have been performed to address the effects of the drift wave instability on blob/filamentary structures characteristic of a tokamak edge. Finally, a helical-shaped potential perturbation is observed to grow exponentially in time while spinning around the filament axis with electron drift frequency.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Final Technical Report: Magnetic Reconnection in High-Energy Laser-Produced Plasmas

This report describes the final results from the DOE Grant DE-SC0007168, “Fast Magnetic Reconnection in HED Laser-Produced Plasmas.” The recent generation of laboratory high-energy-density physics facilities has opened significant physics opportunities for experimentally modeling astrophysical plasmas. The goal of this proposal is to use these new tools to study fundamental problems in plasma physics and plasma astrophysics. Fundamental topics in this area involve study of the generation, amplification, and fate of magnetic fields, which are observed to pervade the plasma universe and govern its evolution. This project combined experiments at DOE laser facilities with kinetic plasma simulation to study these processes. The primary original goal of the project was to study magnetic reconnection using a new experimental platform, colliding magnetized laser-produced plasmas. However through a series of fortuitous discoveries, the work broadened out to allow significant advancement on multiple topics in laboratory astrophysics, including magnetic reconnection, Weibel instability, and collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Dispersed ion structures at the poleward edge of the auroral oval: Low-altitude observations and numerical modeling

We have compared the AUREOL 3 (A3) observations of auroral ion precipitation, particularly ion beams, with the results from the global kinetic model of magnetotail plasma of Ashour-Abdalla et al. (1993). We have identified 101 energetic keV H(+) velocity dispersed precipitating ion structures (VDIS) with fluxes above 10(exp -3) ergs./sq cm./s in the A3 record between the end of 1981 and mid-1984. These beams display a systematic increase in energy with increasing latitude and were observed in a narrow region within less than 1 deg in latitude of the polar cap boundary. The VDIS are the most distinctive feature in the auroral zone of the plasma sheet boundary layer. We report first on a statistical analysis of the possible ralationships between magnetic activity or substorm phase and the VDIS properties. Our particle simulations of the precipitating ions have been extended by using a series of modified versions of the Tsyganenko (1989) magnetic field model and by varying the cross-magnetosphere electric field. In the simulations, plasma from a mantle source is subject to strong nonlinear acceleration, forming beams which flow along the PSBL. Only 3 to 4% of these beams precipitate into the ionosphere to form the VDIS while the majority return to the equatorial plane after mirroring and form the thermalized central plasma sheet. The final energy and the dispersion of the beams in the model depend on the amplitude of the cross-tail electric field. Two unsual observations of low-energy (less than 5 keV) O(+) VDIS, shifted by 4 deg 5 deg in invariant latitude equatorward of H(+) VDIS are analyzed in detail. The sparsity of such O(+) events and the absence of the changes in the flux and frequency of occurrence indicate a solar wind origin for the plasma. Finally, large-scale kinetic modeling, even with its simplifications and assumptions (e.g., static magnetic field, solar wind source), reproduces low-altitude auroral ion features fairly well; it may therefore be presented as an appropriate framework into which data on energization and transport of the hot plasma, obtained in the equatorial plane, could be inserted in the near future.

Bosqued, J. M.↗
DOE OSTI2023

Parallel Algebraic Multigrid for Fusion and Higher-Order PDEs

Multigrid methods play a key role in large-scale scientific simulation because they are among the fastest and most scalable approaches for solving the underlying sparse linear systems of equations that arise from a wide array of Partial Differential Equation (PDE) discretizations. Algebraic multigrid (AMG) is a special type of multigrid method that depends only on the description of the linear system, giving it better portability and broader applicability than geometric multigrid, as it requires no explicit knowledge of the problem geometry. Even though these methods are widely used today, there are still applications where further development is needed. In this report, we focus on PDEs with higher-order terms (e.g., fourth order), concentrating on a PDE that arises in tokamak edge plasma simulations (a tokamak is a machine that confines a plasma using magnetic fields and is believed to be the leading plasma confinement concept for future fusion power plants). General multigrid relaxes a linear system on coarser grids and reverses this process with interpolation, but standard AMG methods struggle with the aforementioned higher-order PDEs. We investigate cyclic coarsening and interpolation heuristics, as well as new iterative approximation methods of refining the solution at each grid to improve the existing multigrid approach. To this end, we ensure that these techniques are transferable to a parallelized setting with LLNL’s supercomputers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Fully Implict Magneto-hydrodynamics Simulations of Coaxial Plasma Accelerators

The resistive Magneto-Hydrodynamic (MHD) model describes the behavior of a strongly ionized plasma in the presence of external electric and magnetic fields. We developed a fully implicit MHD simulation tool to solve the resistive MHD governing equations in the context of a cell-centered finite-volume scheme. The primary objective of this study is to use the fully-implicit algorithm to obtain insights into the plasma acceleration and jet formation processes in Coaxial Plasma accelerators; electromagnetic acceleration devices that utilize self-induced magnetic fields to accelerate thermal plasmas to large velocities. We also carry out plasma-surface simulations in order to study the impact interactions when these high velocity plasma jets impinge on target material surfaces. Scaling studies are carried out to establish some basic functional relationships between the target-stagnation conditions and the current discharged between the coaxial electrodes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗