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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Computer simulations of cosmic-ray diffusion near supernova remnant shock waves

A plasma simulation model was used to study the resonant interactions between streaming cosmic-ray ions and a self-consistent spectrum of Alfven waves, such as might exist in the interstellar medium upstream of a supernova remnant shock wave. The computational model is a hybrid one, in which the background interstellar medium is an MHD fluid and the cosmic-rays are discrete kinetic particles. The particle sources for the electromagnetic fields are obtained by averaging over the fast cyclotron motions. When the perturbed magnetic field is larger than 10 percent of the background field, the macro- and microphysics are no longer correctly predicted by quasi-linear theory. The particles are trapped by the waves and show sharp jumps in their pitch-angles relative to the background magnetic field, and the effective ninety-degree scattering time for diffusion parallel to the background magnetic field is reduced to between 5 and 30 cyclotron periods. Simulation results suggest that Type 1 supernova remnants may be the principal sites of cosmic ray acceleration.

Max, C. E.↗
NASA NTRS1996

Three-Dimensional Electromagnetic Monte Carlo Particle-in-Cell Simulations of Critical Ionization Velocity Experiments in Space

Although the existence of the critical ionization velocity (CIV) is known from laboratory experiments, no agreement has been reached as to whether CIV exists in the natural space environment. In this paper we move towards more realistic models of CIV and present the first fully three-dimensional, electromagnetic particle-in-cell Monte-Carlo collision (PIC-MCC) simulations of typical space-based CIV experiments. In our model, the released neutral gas is taken to be a spherical cloud traveling across a magnetized ambient plasma. Simulations are performed for neutral clouds with various sizes and densities. The effects of the cloud parameters on ionization yield, wave energy growth, electron heating, momentum coupling, and the three-dimensional structure of the newly ionized plasma are discussed. The simulations suggest that the quantitative characteristics of momentum transfers among the ion beam, neutral cloud, and plasma waves is the key indicator of whether CIV can occur in space. The missing factors in space-based CIV experiments may be the conditions necessary for a continuous enhancement of the beam ion momentum. For a typical shaped charge release experiment, favorable CIV conditions may exist only in a very narrow, intermediate spatial region some distance from the release point due to the effects of the cloud density and size. When CIV does occur, the newly ionized plasma from the cloud forms a very complex structure due to the combined forces from the geomagnetic field, the motion induced emf, and the polarization. Hence the detection of CIV also critically depends on the sensor location.

Wang, J.↗
DOE OSTI2024

Plasma edge and scrape-off layer turbulence in gyrokinetic simulations of negative triangularity plasmas

Gyrokinetic simulations in the long-wavelength or drift-kinetic limit are carried out of DIII-D inner-wall-limited (IWL) plasmas to investigate the effect of triangularity on edge and scrape-off layer (SOL) turbulence. The effect of neutral interactions and triangularity on plasma blobs is explored due to the impact blobs can have in setting the SOL width or introducing impurities through interactions with plasma-facing components. Seeded blob simulations with neutrals in shaped SOL scenarios demonstrate that increasing elongation, triangularity, or Shafranov shift decreases radial blob velocities, but neutral interactions have a minor effect. Fully turbulent simulations of DIII-D IWL plasmas include both open- and closed-field-line regions. The negative triangularity (NT) simulation has lower average core T e , lower normalized T e fluctuations, and lower fluxes, but a greater number of coherent structures (blobs) identified with increased size and velocity, on average. Density and electron temperature profiles are within a factor of 2 of experimental values. In conclusion, the increased trapped electron particle fraction in NT simulations is consistent with previous studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Computer simulation of Alfven waves and double layers along auroral magnetic field lines

A plasma simulation has been developed to model interactions between inertial Alfven waves and double layers and to investigate their relative contributions at auroral particle acceleration. We use a novel one-dimensional particle-in-cell code, with periodic boundary conditions, to model the nonlinear excitation of current-driven weak double layers via free energy supplied by an inertial Alfven wave. Analysis of the simulation output shows that double layers are not the agent primarily responsible for electron acceleration. Rather, the inertial Alfven wave accelerates groups of electrons into a steepening beam as it encounters them. As the beam electrons reenter the main distribution, decelerated by anomalous resistive effects, they are replaced by electrons farther downstream. Hence, the particles do not free-stream over the length of the channel. Furthermore, this wave action persists even when the system is linearly stable to ion-acoustic modes, precluding the possibility that this behavior is brought about by the formation of ion-acoustic double layers.

Silberstein, M.↗
DOE OSTI2024

Intrinsic Toroidal Rotation Driven by Turbulent and Neoclassical Processes in Tokamak Plasmas from Global Gyrokinetic Simulations

Gyrokinetic tokamak plasmas can exhibit intrinsic toroidal rotation driven by the residual stress. While most studies have attributed the residual stress to the parallel-momentum flux from the turbulent E × B motion, the parallel-momentum flux from the drift-orbit motion (denoted $Π^D_\parallel$) and the E × B-momentum flux from the E × B motion (denoted $Π_{E×B}$) are often neglected. Here, we use the global total-f gyrokinetic code XGC to study the residual stress in the core and the edge of a DIII-D H-mode plasma. Numerical results show that both $Π^D_\parallel$ and $Π_{E×B}$ make up a significant portion of the residual stress. In particular, $Π^D_\parallel$ in the core is higher than the collisional neoclassical level in the presence of turbulence, while in the edge it represents an outflux of countercurrent momentum even without turbulence. Using a recently developed “orbit-flux” formulation, we show that the higher-than-neoclassical-level $Π^D_\parallel$ in the core is driven by turbulence, while the outflux of countercurrent momentum from the edge is mainly due to collisional ion orbit loss. In conclusion, these results suggest that $Π^D_\parallel$ and $Π_{E×B}$ can be important for the study of intrinsic toroidal rotation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRSDate not supplied

Nanoflare Heating of an X-Ray Bright Point

Nanoflares are thought to be one of the prime candidates that can keep the solar corona to its multimillion kelvin temperature. Individual nanoflares are difficult to detect with the present generation instruments, however their presence can be inferred by comparing the nanoflare heated simulated plasma emissions with the observed emission. Here, we present a simulation of emission from an X-ray Bright Point (XBP) that was observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS), along with concurrent observations from SDO/AIA and Hinode/XRT. We use EBTEL hydrodynamic code to simulate the XBP loops. Length and magnetic field strength of these loops are derived from the potential field extrapolation of the observed photospheric magnetogram by HMI/SDO. Each loop is assumed to be heated by random nanoflares, whose magnitude and frequency are determined by the looplength and magnetic field strength. The simulated outputs are used to predict the intensity of spectrally pure map of Fe-18, Fe-17, Ne-9 ,O-8, O-9, Ne-9 etc, which are then compared with the derived intensity from MaGIXS observation. Further we have predicted the intensity map as observed by AIA and XRT and compared them with the observation. We also estimated the temperature distribution of the XBP from the simulation and found a good agreement with the derived distribution from MaGIXS observation.

coronal heating↗
DOE OSTI2020

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

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

42 ENGINEERING↗
DOE OSTI2024

Modeling a divertor with mid-leg pumping for high-power H-mode scenarios in DIII-D considering E x B drift flows

Edge-plasma simulations of a baffled, long-legged divertor in DIII-D, performed using the multi-fluid code UEDGE, indicate that the position of the detachment front is constrained to the location of the pump duct along the low-field side (LFS) baffle. Simulations including magnetic and E x B drifts were performed for 12.5 MW deuterium plasmas including intrinsic carbon and seeded neon to assess the optimal location of the LFS divertor pump to create a stable detachment front between the target and the X-point. The radiation front position in the simulations, taken to be indicative of the detachment front, can be controlled between the pump and X-point in the favorable magnetic field direction for H-mode access by moving the pump duct location upstream of the target along the LFS baffle. In the unfavorable magnetic field direction, the radial E θ x B drift flows are directed towards the pumping surface, efficiently removing the injected deuterium gas and limiting the sensitivity of the radiation front location to the gas injection rate. The role of pumping rate and drift direction on the pumping efficiency are also found to affect the divertor plasma conditions and detachment front location in UEDGE simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1990

Particle acceleration by intense auroral VLF turbulence

Broadband turbulence in the lower-hybrid to plasma frequency range is found in a variety of forms in the suprauroral region, most notably as auroral hiss and VLF saucers. When the turbulence is intense, it is observed to be associated with ion conics (ions heated transverse to the geomagnetic field) and 'counter-streaming' electron fluxes (heated in both directions parallel to the field). A review is presented of the dispersion and propagation characteristics of whistler resonance-cone waves, which comprise the turbulence. Plasma simulation and mesoscale (Monte Carlo) simulation techniques are used to illustrate the interaction of the ambient plasma with the turbulence. These calculations demonstrate how this interaction results in transverse heating of the ions and parallel heating of the electrons of the plasma, leading to the formation of the observed heated and accelerated particle fluxes.

Retterer, John M.↗
NASA NTRS1990

Simulation study of ion two-stream instability in the auroral acceleration region

Particle in cell plasma simulation technique has been used to study the effect of the ion two-stream instability on hydrogen and oxygen distributions at small relative drifts. Whether heating due to parallel modes which are unstable at low drifts might be sufficient to quench oblique modes which are unstable at higher drifts is investigated. It is found that this is not the case and that the oblique modes contribute significantly to heating.

Gray, Perry C.↗
DOE OSTI2024

Center for Tokamak Transient Simulations (RPI Unstructured Mesh Developments for FES SciDAC4 Partnerships) (Final Report)

The goal of this project was the development of unstructured mesh technologies for fusion simulation codes” for FES SciDAC partnerships and to integrate those technologies into the simulation workflows of those partnerships. Specific developments were executed in support of the following FES SciDAC4 partnerships: Partnership Center for High‐fidelity Boundary Plasma Simulation (HBPS), Center for Integrated Simulation of Fusion Relevant RF (RF‐SciDAC), Center for Plasma Surface Interactions: Predicting the Performance and Impact of Dynamic PFC Surfaces (PSI2), and Center for Tokamak Transient Simulations (CTTS). The key unstructured mesh development areas addressed in this project include (i) methods to most effectively perform PIC calculations on unstructured meshes; (ii) creating meshes for fusion systems accounting for any desired level of geometric complexity and providing physics aware mesh configurations, (iii) adapting unstructured meshes to control the discretization errors, (iv) executing unstructured mesh calculations on GPU accelerated systems, (v) supporting physics/application‐specific PIC operations including surface/wall interactions of particles, (vi) providing infrastructure tools to support the interactions of solvers with unstructured meshes, and (vii) providing advanced methods for coupling plasma physics codes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1983

Numerical simulation of a disk-shaped electron accelerating electrostatic probe

This paper presents a method utilizing the numerical plasma simulation particle-in-cell technique for studying electron accelerating probe characteristics. A probe configuration other than the commonly used plane, cylindrical, or spherical case will be studied; the probe configuration used is that of a disk. Cylindrical geometry is used in the simulation model. One of the advantages of a numerical simulation is the ease with which the details of the system can be examined. Therefore, in addition to the current-voltage characteristics of the static probe, other results are presented such as the time evolution of the plasma properties and studies of the sheath region surrounding the probe.

Nonnast, J. H.↗
NASA NTRS1993

Hybrid simulations of the effects of interstellar pickup hydrogen on the solar wind termination shock

Hybrid (kinetic ions/fluid electrons) plasma simulations are used to study the effects of a population of energetic interstellar pickup hydrogen ions on the solar wind termination shock. The pickup hydrogen is treated as a second ion species in the simulations, and thus the effects of the pick-ups on the shock, as well as the effects of the shock on the pickups, are treated in a fully self-consistent manner. For quasi-perpendicular shocks with 10-20 percent pickup hydrogen the pickup ions manifest themselves in a small foot ahead of the shock ramp caused by pickup ion reflection. For oblique shocks with smaller angles between the field and the shock normal, a large fraction of the pickup ions are reflected and move back upstream where they excite large amplitude magnetosonic waves which steepen into shocklets. These backstreaming pickup ions may provide advance warning of a spacecraft encounter with the termination shock.

Liewer, P. C.↗
DOE OSTI2025

Electron Cyclotron Emission–Based Separatrix Identification in ITER with OMFIT Synthetic Modeling

Accurate determination of the separatrix location is essential for understanding edge plasma behavior and optimizing confinement in tokamaks, especially in next-generation devices such as ITER. In this study, a synthetic microwave diagnostics module was developed and implemented in the OMFIT framework to assess the feasibility of an electron cyclotron emission–based separatrix detection method in ITER plasmas. Simulations were carried out using ITER H-mode equilibrium scenarios with different plasma density profiles and different pedestal widths. Here, the results show that the electron emission temperature profiles consistently exhibit an inversion pattern near the edge, with a well-defined minimum point that could serve as a proxy for the separatrix location. However, unlike in the DIII-D, the minimum point in ITER is systematically offset by approximately 2 cm into the scrape-off layer, independent of density or pedestal width, which is within the radial resolution range (2 to 5 cm) determined by the 500-MHz channel spacing. While the method does not provide the exact separatrix location, it offers a reliable indicator of the boundary region and has potential applications for real-time boundary monitoring in ITER and other future fusion devices.

Electron cyclotron emission↗
NASA NTRS1974

Computational study of nonlinear plasma waves

A low-noise plasma simulation model is developed and applied to a series of linear and nonlinear problems associated with electrostatic wave propagation in a one-dimensional, collisionless, Maxwellian plasma, in the absence of magnetic field. It is demonstrated that use of the hybrid simulation model allows economical studies to be carried out in both the linear and nonlinear regimes with better quantitative results, for comparable computing time, than can be obtained by conventional particle simulation models, or direct solution of the Vlasov equation. The characteristics of the hybrid simulation model itself are first investigated, and it is shown to be capable of verifying the theoretical linear dispersion relation at wave energy levels as low as .000001 of the plasma thermal energy. Having established the validity of the hybrid simulation model, it is then used to study the nonlinear dynamics of monochromatic wave, sideband instability due to trapped particles, and satellite growth.

Matsuda, Y.↗