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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Conservative discontinuous Galerkin interpolation: Sheared boundary conditions

Local studies of accretion disks and laboratory magnetized plasmas employ analytical coordinate mappings that introduce sheared boundary conditions (BCs). Here, we present a discontinuous Galerkin (DG) algorithm to apply such BCs based on projections and quadrature-free integration. The procedure is high-order accurate, preserves moments exactly and works in multiple dimensions. Tests of increasing complexity are provided, beginning with translations of one and two dimensional fields, followed by 3D and 5D simulations with sheared (twist-shift) BCs. Results show that the algorithm is (p + 1)-order accurate in the DG representation and (p + 2) -order accurate in the cell averages, with p being the order of the polynomial basis. Quantification of the algorithm's hyperdiffusion and discussion of aliasing errors are given. This technique enables conservative local simulations of plasma turbulence with DG, not possible until now.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Transport and losses of fusion-born alpha particles in the presence of tearing modes using the new Toroidal Accelerated Particle Simulator (TAPaS)

Here, the transport and losses of fusion-born alpha particles are studied in the presence of a single-helicity tearing mode, characterized by (m = 2, n = 1). The analysis is performed by means of the recently developed Toroidal Accelerated Particle Simulator (TAPaS). Although such modes have usually been believed to result only in a local flattening of the radial profiles, it is shown that the density profile can exhibit a global modification leading to significant losses of alpha particles. This is due to the fact that, although the magnetic field does not exhibit any chaotic behaviour, the trajectories of alpha particles do, as revealed by their Poincaré maps. Such results are in qualitative agreement with past observations and simulations of energetic particles generated by neutral beam injection in TFTR, DIII-D and AUG tokamaks. In-depth analysis is carried out to characterize the impact of the tearing mode on the transport and losses of fusion-born alpha-particles with a realistic density profile. The impact of the amplitude is evidenced. Moreover, the effect of the island rotation frequency is assessed based on a detailed analysis of the linear resonances in phase-space, in agreement with the simulation results. Finally, the probability density function of the exit time has been computed and the transport of alpha particles has been found to be anomalous.

anomalous transport↗
NASA NTRS1980

Some features of inverted-V events as seen from simulated double layers

Results from a numerical simulation of a double layer show some features similar to those of inverted-V events. The strong heating of thermal and precipitating electrons is observed along with extremely low frequency fluctuations found during inverted-V events. It is suggested that after the acceleration of auroral electrons by the double layer, the precipitating free electrons are heated by the nonlinear effects of the electron beam plasma instability. Fluctuations and pulsations of auroral electron fluxes during auroral events are caused by a relaxation type of oscillation. The finite extent of one dimensional plasma is simulated by solving the Vlasov and Poisson equations as an initial and boundary value problem.

Singh, N.↗
NASA NTRS2000

Magnetic Flux Compression Concept for Nuclear Pulse Propulsion and Power

The desire for fast, efficient interplanetary transport requires propulsion systems having short acceleration times and very high specific impulse attributes. Unfortunately, most highly efficient propulsion systems which are within the capabilities of present day technologies are either very heavy or yield very low impulse such that the acceleration time to final velocity is too long to be of lasting interest, One exception, the nuclear thermal thruster, could achieve the desired acceleration but it would require inordinately large mass ratios to reach the range of desired final velocities. An alternative approach, among several competing concepts that are beyond our modern technical capabilities, is a pulsed thermonuclear device utilizing microfusion detonations. In this paper, we examine the feasibility of an innovative magnetic flux compression concept for utilizing microfusion detonations, assuming that such low yield nuclear bursts can be realized in practice. In this concept, a magnetic field is compressed between an expanding detonation driven diamagnetic plasma and a stationary structure formed from a high temperature superconductor (HTSC). In general, we are interested in accomplishing two important functions: (1) collimation of a hot diamagnetic plasma for direct thrust production; and (2) pulse power generation for dense plasma ignition. For the purposes of this research, it is assumed that rnicrofusion detonation technology may become available within a few decades, and that this approach could capitalize on recent advances in inertial confinement fusion ICF) technologies including magnetized target concepts and antimatter initiated nuclear detonations. The charged particle expansion velocity in these detonations can be on the order of 10 (exp 6)- 10 (exp 7) meters per second, and, if effectively collimated by a magnetic nozzle, can yield the Isp and the acceleration levels needed for practical interplanetary spaceflight. The ability to ignite pure fusion micro-bursts with reasonable levels of input energy is an equally challenging scientific problem. It remains to be seen, however, whether an effective ignition driver can be developed which meets the requirements for practical spaceflight application (namely high power density, compactness, low weight, and low cost). In this paper, system level performance and design issues are examined including generator performance, magnetic flux compression processes, magnetic diffusion processes, high temperature superconductor (HTSC) material properties, plasmadynamic processes, detonation plasma expansion processes, magnetohydrodynamic instabilities, magnetic nozzle performance, and thrust production performance. Representative generator performance calculations based on a simplified skin layer formulation are presented as well as the results of exploratory small-scale laboratory experiments on magnetic flux diffusion in HTSC materials. In addition, planned follow-on scientific feasibility experiments are described which utilize high explosive detonations and high energy gas discharges to simulate the plasma conditions associated with thermonuclear micro-detonations.

Litchford, Ronald J.↗
DOE OSTI2018

Alternative model of space-charge-limited thermionic current flow through a plasma

It is widely assumed that thermionic current flow through a plasma is limited by a “space-charge-limited” (SCL) cathode sheath that consumes the hot cathode's negative bias and accelerates upstream ions into the cathode. In this paper, we formulate a fundamentally different current-limited mode. In the “inverse” mode, the potentials of both electrodes are above the plasma potential, so that the plasma ions are confined. The bias is consumed by the anode sheath. There is no potential gradient in the neutral plasma region from resistivity or presheath. The inverse cathode sheath pulls some thermoelectrons back to the cathode, thereby limiting the circuit current. Thermoelectrons entering the zero-field plasma region that undergo collisions may also be sent back to the cathode, further attenuating the circuit current. In planar geometry, the plasma density is shown to vary linearly across the electrode gap. A continuum kinetic planar plasma diode simulation model is set up to compare the properties of current modes with classical, conventional SCL, and inverse cathode sheaths. SCL modes can exist only if charge-exchange collisions are turned off in the potential well of the virtual cathode to prevent ion trapping. With the collisions, the current-limited equilibrium must be inverse. Inverse operating modes should therefore be present or possible in many plasma devices that rely on hot cathodes. Evidence from past experiments is discussed. Finally, the inverse mode may offer opportunities to minimize sputtering and power consumption that were not previously explored due to the common assumption of SCL sheaths.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Gyrokinetic particle simulation of a field reversed configuration

Gyrokinetic particle simulation of the field-reversed configuration (FRC) has been developed using the gyrokinetic toroidal code (GTC). The magnetohydrodynamic equilibrium is mapped from cylindrical coordinates to Boozer coordinates for the FRC core and scrape-off layer (SOL), respectively. A field-aligned mesh is constructed for solving self-consistent electric fields using a semi-spectral solver in a partial torus FRC geometry. This new simulation capability has been successfully verified and driftwave instability in the FRC has been studied using the gyrokinetic simulation for the first time. Initial GTC simulations find that in the FRC core, the ion-scale driftwave is stabilized by the large ion gyroradius. In the SOL, the driftwave is unstable on both ion and electron scales.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Sparsely corrupted stimulated scattering signals recovery by iterative reweighted continuous basis pursuit

In this paper, we consider the problem of extracting the desired signals from noisy measurements. This is a classical problem of signal recovery which is of paramount importance in inertial confinement fusion. To accomplish this task, we develop a tractable algorithm based on continuous basis pursuit and reweighted ℓ{sub 1}-minimization. By modeling the observed signals as superposition of scale time-shifted copies of theoretical waveform, structured noise, and unstructured noise on a finite time interval, a sparse optimization problem is obtained. We propose to solve this problem through an iterative procedure that alternates between convex optimization to estimate the amplitude, and local optimization to estimate the dictionary. The performance of the method was evaluated both numerically and experimentally. Numerically, we recovered theoretical signals embedded in increasing amounts of unstructured noise and compared the results with those obtained through popular denoising methods. We also applied the proposed method to a set of actual experimental data acquired from the Shenguang-II laser whose energy was below the detector noise-equivalent energy. Both simulation and experiments show that the proposed method improves the signal recovery performance and extends the dynamic detection range of detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2013

Micro-sphere layered targets efficiency in laser driven proton acceleration

Proton acceleration from the interaction of high contrast, 25 fs laser pulses at >10{sup 19} W/cm{sup 2} intensity with plastic foils covered with a single layer of regularly packed micro-spheres has been investigated experimentally. The proton cut-off energy has been measured as a function of the micro-sphere size and laser incidence angle for different substrate thickness, and for both P and S polarization. The presence of micro-spheres with a size comparable to the laser wavelength allows to increase the proton cut-off energy for both polarizations at small angles of incidence (10∘). For large angles of incidence, however, proton energy enhancement with respect to flat targets is absent. Analysis of electron trajectories in particle-in-cell simulations highlights the role of the surface geometry in the heating of electrons.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron Modulational Instability in the Strong Turbulent Regime for an Electron Beam Propagating in a Background Plasma

We study collective processes for an electron beam propagating through a background plasma using simulations and analytical theory. A new regime where the instability of a Langmuir wave packet can grow locally much faster than ion frequency is clearly identified. The key feature of this new regime is an electron modulational instability that rapidly creates a local Langmuir wave packet, which in its turn produces local charge separation and strong ion density perturbations because of the action of the ponderomotive force, such that the beam-plasma wave interaction stops being resonant. Three evolution stages of the process and observed periodic burst features are discussed. In conclusion, different physical regimes in the plasma and beam parameter space are demonstrated for the first time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

A phase-shift-periodic parallel boundary condition for low-magnetic-shear scenarios

Abstract We formulate a generalized periodic boundary condition as a limit of the standard twist-and-shift parallel boundary condition that is suitable for simulations of plasmas with low magnetic shear. This is done by applying a phase shift in the binormal direction when crossing the parallel boundary. While this phase shift can be set to zero without loss of generality in the local flux-tube limit when employing the twist-and-shift boundary condition, we show that this is not the most general case when employing periodic parallel boundaries, and may not even be the most desirable. A non-zero phase shift can be used to avoid the convective cells that plague simulations of the three-dimensional Hasegawa–Wakatani system, and is shown to have measurable effects in periodic low-magnetic-shear gyrokinetic simulations. We propose a numerical program where a sampling of periodic simulations at random pseudo-irrational flux surfaces are used to determine physical observables in a statistical sense. This approach can serve as an alternative to applying the twist-and-shift boundary condition to low-magnetic-shear scenarios, which, while more straightforward, can be computationally demanding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Fluid vs. kinetic magnetic reconnection with strong guide fields

The fast rates of magnetic reconnection found in both nature and experiments are important to understand theoretically. Recently, it was demonstrated that two-fluid magnetic reconnection remains fast in the strong guide field regime, regardless of the presence of fast-dispersive waves. This conclusion is in agreement with recent results from kinetic simulations, and is in contradiction to the findings in an earlier two-fluid study, where it was suggested that fast-dispersive waves are necessary for fast reconnection. In this paper, we give a more detailed derivation of the analytic model presented in a recent letter and present additional simulation results to support the conclusions that the magnetic reconnection rate in this regime is independent of both collisional dissipation and system-size. In particular, we present a detailed comparison between fluid and kinetic simulations, finding good agreement in both the reconnection rate and overall length of the current layer. Finally, we revisit the earlier two-fluid study, which arrived at different conclusions, and suggest an alternative interpretation for the numerical results presented therein.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Global gyrokinetic simulations of the H-mode tokamak edge pedestal

Global gyrokinetic simulations of DIII-D H-mode edge pedestal show two types of instabilities may exist approaching the onset of edge localized modes: an intermediate-n, high frequency mode which we identify as the “kinetic peeling ballooning mode (KPBM),” and a high-n, low frequency mode. Our previous study [W. Wan et al., Phys. Rev. Lett. 109, 185004 (2012)] has shown that when the safety factor profile is flattened around the steep pressure gradient region, the high-n mode is clearly kinetic ballooning mode and becomes the dominant instability. Otherwise, the KPBM dominates. Here, the properties of the two instabilities are studied by varying the density and temperature profiles. It is found that the KPBM is destabilized by density and ion temperature gradient, and the high-n mode is mostly destabilized by electron temperature gradient. Nonlinear simulations with the KPBM saturate at high levels. The equilibrium radial electric field (E{sub r}) reduces the transport. The effect of the parallel equilibrium current is found to be weak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Fast ignition: Dependence of the ignition energy on source and target parameters for particle-in-cell-modelled energy and angular distributions of the fast electrons

The energy and angular distributions of the fast electrons predicted by particle-in-cell (PIC) simulations differ from those historically assumed in ignition designs of the fast ignition scheme. Using a particular 3D PIC calculation, we show how the ignition energy varies as a function of source-fuel distance, source size, and density of the pre-compressed fuel. The large divergence of the electron beam implies that the ignition energy scales with density more weakly than the ρ{sup −2} scaling for an idealized beam [S. Atzeni, Phys. Plasmas 6, 3316 (1999)], for any realistic source that is at some distance from the dense deuterium-tritium fuel. Due to the strong dependence of ignition energy with source-fuel distance, the use of magnetic or electric fields seems essential for the purpose of decreasing the ignition energy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

A direct density modulation cathode in magnetron

A direct Density Modulation Cathode (DMC) in magnetron is proposed in this paper. By removing the velocity modulation process, electron spokes corresponding to the dominant operating mode can be quickly formed when the DMC is used. Based on theoretical analysis, particle-in-cell simulations and experimental investigations are carried out for demonstration. The final results show that compared with conventional solid cathode and transparent cathode, the DMC can increase 68% and even 146% of relative microwave widths, respectively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Stabilization of high-compression, indirect-drive inertial confinement fusion implosions using a 4-shock adiabat-shaped drive

Hydrodynamic instabilities and poor fuel compression are major factors for capsule performance degradation in ignition experiments on the National Ignition Facility. Using a recently developed laser drive profile with a decaying first shock to tune the ablative Richtmyer-Meshkov (ARM) instability and subsequent in-flight Rayleigh-Taylor growth, we have demonstrated reduced growth compared to the standard ignition pulse whilst maintaining conditions for a low fuel adiabat needed for increased compression. Using in-flight x-ray radiography of pre-machined modulations, the first growth measurements using this new ARM-tuned drive have demonstrated instability growth reduction of ∼4× compared to the original design at a convergence ratio of ∼2. Corresponding simulations give a fuel adiabat of ∼1.6, similar to the original goal and consistent with ignition requirements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Numerical study on Rayleigh-Taylor instabilities in the lightning return stroke

The Rayleigh-Taylor (R-T) instabilities are important hydrodynamics and magnetohydrodynamics (MHD) phenomena that are found in systems in high energy density physics and normal fluids. The formation and evolution of the R-T instability at channel boundary during back-flow of the lightning return stroke are analyzed using the linear perturbation theory and normal mode analysis methods, and the linear growth rate of the R-T instability in typical condition for lightning return stroke channel is obtained. Then, the R-T instability phenomena of lightning return stroke are simulated using a two-dimensional Eulerian finite volumes resistive radiation MHD code. The numerical results show that the evolution characteristics of the R-T instability in the early stage of back-flow are consistent with theoretical predictions obtained by linear analysis. The simulation also yields more evolution characteristics for the R-T instability beyond the linear theory. The results of this work apply to some observed features of the return stroke channel and further advance previous theoretical and experimental work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

A possible mechanism of the enhancement and maintenance of the shear magnetic field component in the current sheet of the Earth’s magnetotail

The influence of the shear magnetic field component, which is directed along the electric current in the current sheet (CS) of the Earth’s magnetotail and enhanced near the neutral plane of the CS, on the nonadiabatic dynamics of ions interacting with the CS is studied. The results of simulation of the nonadiabatic ion motion in the prescribed magnetic configuration similar to that observed in the magnetotail CS by the CLUSTER spacecraft demonstrated that, in the presence of some initial shear magnetic field, the north-south asymmetry in the ion reflection/refraction in the CS is observed. This asymmetry leads to the formation of an additional current system formed by the oppositely directed electric currents flowing in the northern and southern parts of the plasma sheet in the planes tangential to the CS plane and in the direction perpendicular to the direction of the electric current in the CS. The formation of this current system perhaps is responsible for the enhancement and further maintenance of the shear magnetic field near the neutral plane of the CS. The CS structure and ion dynamics observed in 17 intervals of the CS crossings by the CLUSTER spacecraft is analyzed. In these intervals, the shear magnetic field was increased near the neutral plane of the CS, so that the bell-shaped spatial distribution of this field across the CS plane was observed. The results of the present analysis confirm the suggested scenario of the enhancement of the shear magnetic field near the neutral plane of the CS due to the peculiarities of the nonadiabatic ion dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗