Search NASA⌕ Search

SEARCH · Search NASA

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Implications of advanced collision operators for gyrokinetic simulation

In this paper, we explore both the potential improvements and pitfalls that arise when using advanced collision models in gyrokinetic simulations of plasma microinstabilities. Comparisons are made between the simple-but-standard electron Lorentz operator and specific variations of the advanced Sugama operator. The Sugama operator describes multi-species collisions including energy diffusion, momentum and energy conservation terms, and is valid for arbitrary wavelength. We report scans over collision frequency for both low and high ${k}_{\theta }{\rho }_{s}$ modes, with relevance for multiscale simulations that couple ion and electron scale physics. The influence of the ion–ion collision terms—not retained in the electron Lorentz model—on the damping of zonal flows is also explored. Collision frequency scans for linear and nonlinear simulations of ion-temperature-gradient instabilities including impurity ions are presented. Finally, implications for modeling turbulence in the highly collisional edge are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Laser-driven, ion-scale magnetospheres in laboratory plasmas. II. Particle-in-cell simulations

Ion-scale magnetospheres have been observed around comets, weakly magnetized asteroids, and localized regions on the Moon and provide a unique environment to study kinetic-scale plasma physics, in particular in the collision-less regime. In this work, we present the results of particle-in-cell simulations that replicate recent experiments on the large plasma device at the University of California, Los Angeles. Using high-repetition rate lasers, ion-scale magnetospheres were created to drive a plasma flow into a dipolar magnetic field embedded in a uniform background magnetic field. The simulations are employed to evolve idealized 2D configurations of the experiments, study highly resolved, volumetric datasets, and determine the magnetospheric structure, magnetopause location, and kinetic-scale structures of the plasma current distribution. We show the formation of a magnetic cavity and a magnetic compression in the magnetospheric region, and two main current structures in the dayside of the magnetic obstacle: the diamagnetic current, supported by the driver plasma flow, and the current associated with the magnetopause, supported by both the background and driver plasmas with some time-dependence. From multiple parameter scans, we show a reflection of the magnetic compression, bounded by the length of the driver plasma, and a higher separation of the main current structures for lower dipolar magnetic moments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Electrical and optical properties of warm dense beryllium along the principal Hugoniot

The electrical and optical properties of warm dense beryllium along the principal Hugoniot for temperatures from 0.95 eV to 10.65 eV and densities from 3.8 to 6.0 g/cm{sup 3} are investigated by using quantum molecular dynamics (QMD) simulations combined with the Kubo-Greenwood formulation. The dc conductivity σ{sub dc} and the ionization fraction are yielded by fitting the optical conductivity with the Drude-Smith model. The first-principles transport coefficients are compared with results of the Lee-More model and the Brysk model [Plasma Phys. 17, 473 (1975)]. Compared with the QMD result, the Lee-More model underestimates σ{sub dc} by 87% at low temperatures, approaches the QMD result gradually with the temperature rising, yet still underestimates σ{sub dc} by 49% corresponding to the temperature 10.65 eV. In the whole temperature range under investigation, the Brysk model overestimates the electronic thermal conductivity κ while the Lee-More model underestimates κ. The differences are reduced with the temperature increasing. At the temperature 10.65 eV, the Brysk κ is still around twice as large as the QMD result, and the Lee-More κ is smaller than the QMD data by about 40%. In addition, QMD Rosseland mean opacities are shown to be three orders of magnitude larger than results of the average-atom model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1998

Object-Oriented Scientific Programming with Fortran 90

Fortran 90 is a modern language that introduces many important new features beneficial for scientific programming. We discuss our experiences in plasma particle simulation and unstructured adaptive mesh refinement on supercomputers, illustrating the features of Fortran 90 that support the object-oriented methodology.

Fortran 90 Supercomputing Object-Oriented Adaptive↗
NASA NTRS1997

Boundary Layer Formation in the Magnetotail: Geotail Observations and Comparisons with a Global MHD Simulation

We present Geotail plasma and field observations from the middle magnetotail near X(sub GSE) = -46 R(sub E) for the time period 1400 to 1800 UT on December 14, 1994. During that period, the Wind satellite monitored the solar wind plasma and interplanetary magnetic field (IMF) upstream of the bow shock. The IMF was northward and the plasma parameters near average. Geotail observed slow tailward flows and a northward field. The plasma and field parameters indicate that Geotail is either in the plasma sheet or in a boundary layer. We used the Wind solar wind plasma and IMF data as input for a global simulation of that time interval. Comparison of the simulation results with the observational data show very good overall agreement of the magnitudes of the plasma and field parameters. In particular, the simulation reproduces the slow tailward flows and northward field found at Geotail. Small scale temporal, variations are less well reproduced. The simulation shows the formation of a broad boundary layer (which we call tail flank boundary layer, TFBL) that consists of closed flux which is formed by magnetic magnetic reconnection of IMF and lobe field lines. The simulation results indicate that Geotail is located very close to the TFBL and may have entered the TFBL proper. We show that the TFBL plays an important role in energy transport from the solar wind into the magnetosphere during northward IMF conditions.

Raeder, J.↗
DOE OSTI2024

Turbulence spreading effects on the ELM size and SOL width

BOUT++ turbulence simulations were performed to investigate the impact of turbulence spreading on the edge localized mode (ELM) size and divertor heat flux width (λ q ) broadening in small ELM regimes. Here, this study is motivated by EAST experiments. BOUT++ linear simulations of a pedestal radial electric field (E r ) scan show that the dominant toroidal number mode (n) shifts from high-n to low-n, with a narrow mode spectrum, and the maximum linear growth rate increases as the pedestal E r well deepens. The nonlinear simulations show that as the net E × B pedestal flow increases, the pressure fluctuation level and its inward penetration beyond the top of the pedestal both increase. This leads to a transition from small ELMs to large ELMs. Both inward and outward turbulence spreading are sensitive to the scrape-off-layer (SOL) plasma profiles. The inward turbulence spreading increases for the steep SOL profiles, leading to increasing pedestal energy loss in the small ELM regime. The SOL width (λ q ) is significantly broadened progressing from the ELM-free to small ELM regime, due to the onset of strong radial turbulent transport. The extent of the SOL width (λ q ) broadening depends strongly on outward turbulence spreading. The fluctuation energy intensity flux Γ ε at the separatrix can be enhanced by increasing either pedestal E r flow shear or local SOL pressure gradient. The λ q is broadened as the fluctuation energy intensity flux Γ ε at the last close flux surface (LCFS) increases. Local SOL E × B flow shear will restrain outward turbulence spreading and the associated heat flux width broadening. Operating in H-mode with small ELMs has the potential to solve two critical problems: reducing the ELM size and broadening the SOL width.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

Target compressions by working fluids driven with solid liner implosions

Compression by a spherical solid liner of a gold target surrounded by a hydrogen plasma is simulated. Two-dimensional simulations that treat only a subset of the physics included in the one-dimensional code were performed in an attempt to assess multidimensional effects. A one-dimensional numerical code has been developed to study the effects of thermal radiation and conduction. Results of pressure, density, and energy deposited for different initial plasma conditions are presented and discussed. Results from both one- and two-dimensional codes show that the average target density at peak compression is 39-43 g/cu cm, using the SHIVA Star facility at 90 kV discharge.

Chiang, P.-R.↗
DOE OSTI2017

Invited Review Article: Gas puff imaging diagnostics of edge plasma turbulence in magnetic fusion devices

Gas puff imaging (GPI) is a diagnostic of plasma turbulence which uses a puff of neutral gas at the plasma edge to increase the local visible light emission for improved space-time resolution of plasma fluctuations. This paper reviews gas puff imaging diagnostics of edge plasma turbulence in magnetic fusion research, with a focus on the instrumentation, diagnostic cross-checks, and interpretation issues. The gas puff imaging hardware, optics, and detectors are described for about 10 GPI systems implemented over the past similar to 15 years. Comparison of GPI results with other edge turbulence diagnostic results is described, and many common features are observed. Here, several issues in the interpretation of GPI measurements are discussed, and potential improvements in hardware and modeling are suggested.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Experimental quantification of the impact of heterogeneous mix on thermonuclear burn

In inertial confinement fusion, deuterium–tritium (DT) fuel is brought to densities and temperatures where fusion ignition occurs. However, mixing of the ablator material into the fuel may prevent ignition by diluting and cooling the fuel. MARBLE experiments at the National Ignition Facility provide new insight into how mixing affects thermonuclear burn. These experiments use laser-driven capsules containing deuterated plastic foam and tritium gas. Embedded within the foam are voids of known sizes and locations, which control the degree of heterogeneity of the fuel. Initially, the reactants are separated, with tritium concentrated in the voids and deuterium in the foam. During the implosion, mixing occurs between the foam and gas materials, leading to DT fusion reactions in the mixed region. Here, it is shown that by measuring the ratios of DT and deuterium–deuterium neutron yields for different macropore sizes and gas compositions, the effects of mix heterogeneity on thermonuclear burn may be quantified, supporting an improved understanding of these effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Alpha channeling with high-field launch of lower hybrid waves

Although lower hybrid waves are effective at driving currents in present-day tokamaks, they are expected to interact strongly with high-energy particles in extrapolating to reactors. In the presence of a radial alpha particle birth gradient, this interaction can take the form of wave amplification rather than damping. While it is known that this amplification more easily occurs when launching from the tokamak high-field side, the extent of this amplification has not been made quantitative. Here, by tracing rays launched from the high-field-side of a tokamak, the required radial gradients to achieve amplification are calculated for a temperature and density regime consistent with a hot-ion-mode fusion reactor. These simulations, while valid only in the linear regime of wave amplification, nonetheless illustrate the possibilities for wave amplification using high-field launch of the lower hybrid wave.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Kinetic simulations of laser parametric amplification in magnetized plasmas

Laser pulse compression using magnetized resonance near the upper-hybrid frequency is promising for achieving higher output intensity in regimes previously thought impossible using unmagnetized plasmas. Using one dimensional particle-in-cell simulations, we verify that, by partially replacing plasma with an external transverse magnetic field of megagauss scale, the output pulse can be intensified by a factor of a few, due to the increased allowable amplification time despite a decreased growth rate. Further improvement is impeded by the generation of an electromagnetic wakefield, to which the amplified pulse loses more energy than it does in the unmagnetized case. This limitation can however be circumvented by the use of a stronger pump. In contrast to unmagnetized compression, the magnetized amplification remains efficient when the pump intensity is well above the wavebreaking threshold, until a higher phase-mixing threshold is exceeded. Here, this surprising resilience to wavebreaking in magnetized plasma is of great benefit for magnetized compression.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Evidence of electron microbunching in laser-driven modulated downramp injection and prospects for beam-driven implementation

Plasma accelerators can generate high-energy, high-brightness electron beams over centimeter-scale distances, offering novel pathways to compact x-ray free-electron lasers. Generating beams pre-bunched at the desired radiation wavelength would significantly enhance longitudinal coherence and reduce saturation length. Plasma density-modulated downramp injection offers an in-situ way to generate such beams with nanometer-scale bunching. Here we report the first experimental evidence of this mechanism in a laser-driven wakefield accelerator, showing that modulated density downramps generate modulated electron energy spectra absent in unmodulated cases. Particle-in-cell simulations reproduce these observations and reveal bunching factors of 0.05 at 0.4 μm, with a compression factor of approximately 7. Building on this demonstration, we propose a beam-driven implementation for FACET-II to generate multi-GeV beams pre-bunched at hundreds of nanometers wavelength, with sub-micrometer emittance, kiloampere peak current, and sub-percent slice energy spread. Two-stage magnetic compression enables tunable bunching from optical to extreme ultraviolet wavelengths while achieving peak currents exceeding 100 kA. Coherent transition radiation calculations confirm diagnostic feasibility. This approach offers a promising path towards compact, high-energy pre-bunched electron sources for advanced photon science applications.

electron microbunching↗
DOE OSTI2015

Numerical simulations of stable explosive-emission center plasma expansion in vacuum

We present the results of our numerical simulations of the stable cathode-flare expansion stage in the vacuum-breakdown switching period. We used the model of two-fluid hydrodynamics with corrections that allowed it to be applied to a consistent description of the electron and ion transport in plasma and the electron transport in a plasma-anode vacuum gap without imposing any joining conditions at the expanding plasma boundary. We have analyzed the dynamics of the density, directed velocity, and temperature profiles for both components and the influence of various factors on the plasma-cloud expansion velocity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2004

Complexity and Intermittent Turbulence in Space Plasmas

Sporadic and localized interactions of coherent structures arising from plasma resonances can be the origin of "complexity" of the coexistence of non- propagating spatiotemporal fluctuations and propagating modes in space plasmas. Numerical simulation results are presented to demonstrate the intermittent character of the non-propagating fluctuations. The technique of the dynamic renormalization-group is introduced and applied to the study of scale invariance of such type of multiscale fluctuations. We also demonstrate that the particle interactions with the intermittent turbulence can lead to the efficient energization of the plasma populations. An example related to the ion acceleration processes in the auroral zone is provided.

Chang, Tom↗
DOE OSTI2017

Monte Carlo simulation of excitation and ionization collisions with complexity reduction

Kinetic simulation of plasmas with detailed excitation and ionization collisions presents a significant computational challenge due to the multiscale feature of the collisional rates. In the present work, we propose a complexity reduction method based on atomic level grouping for modeling excitation and ionization collisions. High order of accuracy of the reduction method is realized by allowing an internal distribution within each group. We apply the reduction method to the standard Monte Carlo collision algorithm to model an atomic Hydrogen plasma. Numerical results suggest that the stiffness of the collisional kinetics can be significantly reduced with minimal loss in accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

Kinetic structure of rotational discontinuities: Implications for the magnetopause

Magnetic field rotations in the high ion beta magnetosheath that are part of the magnetopause structure are expected to have only a small normal component. We have studied the properties of rotational discontinuities (RDs) under these conditions, viewed as the limit of weak intermediate shocks (ISs), by performing hybrid simulations with a reflecting wall boundary condition (piston method). With this dynamic formation, the sense and size of rotation are not arbitrarily predetermined, but rather evolve from the given upstream (magnetosheath) and downstream (magnetospheric) boundary conditions, similar to what takes place at the magnetopause. This work focuses on several aspects: the observed minimum shear of RDs, their width, their internal signature, and their relation to ISs in isotropic plasmas. Our simulation results are in agreement with the minimum shear observations, that is, the RDs choose the sense of rotation that corresponds to the minimum angle between the upstream and downstream field vector. The RDs are stable, with a unique scale size. Typical gradient scale half widths are one to four ion inertial lengths with a total width up to ten times of that, in agreement with magnetopause observations. We develop a generalized fluid theory of RDs and discuss the characteristic internal signatures of the rotational layer, comparing the kinetic simulation results to predictions from the generalized fluid theory. The results show that ion inertia, anisotropic pressure, finite Larmor radius effects, nonzero ion heat flux, and reflected ions all contribute to the signatures of RDs on kinetic scales. The RDs may have upstream or downstream wave trains, which become weak for high ion beta and small normal components of the magnetic field. We explain the presence and direction of wave trains in terms of the kinetic properties of the Alfven/ion-cyclotron mode. Away from the RD limit there is a smooth transition to weak intermediate shocks, which have small jumps close to expected Rankine-Hugoniot values. Apart from that, there are few kinetic plasma signatures that distinguish RDs from their neighboring ISs. However, noncoplanar ISs evolve in time into thin RDs. Using the properties of RDs and ISs, we make specific suggestions how these discontinuities can be distinguished observationally in the case of an isotropic plasma.

Krauss-Varban, D.↗
DOE OSTI2019

Simulation study of CO 2 laser-plasma interactions and self-modulated wakefield acceleration

3D numerical simulations of the interaction of a powerful CO 2 laser with hydrogen jets demonstrating the role of ionization on the characteristics of induced wakes are presented. Simulations using SPACE, a parallel relativistic Particle-in-Cell code, are performed in support of the plasma wakefield accelerator experiments being conducted at the Brookhaven National Laboratory (BNL) Accelerator Test Facility (ATF). A novelty of the SPACE code is its set of efficient atomic physics algorithms that compute ionization and recombination rates on the grid and transfer them to particles. The influence of ionization on the spectrum of pump laser has been studied for a range of gas densities.Simulations reproduce both Stokes and anti-Stokes shifts in the spectrum of pump laser, similar to those observed in experiments in the spectrum of probe. A good agreement has been achieved with the experiments on the effect of variation in gas density on Stokes/anti-Stokes intensity. In addition, self-injection and trapping of electrons into the self-modulated wakes has been observed and analyzed. Here, the experimentally validated code SPACE will be used for predictive simulations to guide future experiments at BNL ATF.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Multi-camera imaging to characterize jet and liner uniformity on the Plasma Liner Experiment (PLX)

At Plasma Liner Experiment, a set of 36 coaxial plasma guns are deployed quasi-uniformly over a 9 ft diameter spherical chamber and are used to form a high-Z spherically compressive plasma liner. Simulations indicate that for the concept to ultimately achieve optimal target density and temperature, a high degree of timing uniformity is required between all guns. Therefore, to aid in quantifying and correcting gun-to-gun nonuniformities, a key diagnostic will consist of up to six fisheye-view CCD cameras positioned inside the main chamber such that each has all plasma guns within its view. The individual cameras can be triggered at different times to determine each plasma jet’s muzzle velocity and structure for different operating conditions. This camera array is currently under development, and the implementation needs and challenges for this camera array are discussed here. Additionally, we detail the analysis methodology for determining jet-to-jet uniformity deviations and how we can correct them, thereby improving overall liner uniformity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗