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,495 records · Page 83
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Influence of plasma surface interactions on tokamak startup

The startup phase of a tokamak is a complex phenomenon involving burnthrough of the low-Z impurities and rampup of I{sub p}, the plasma current. The design considerations of a tokamak are closely connected with the startup modeling. Plasma evolution is analysed using a zero-dimensional model. The particle and energy balance is considered of two subclasses of plasmas which are penetrable by neutral gas, together with another component, neutrals trapped in the wall. The first subclass includes plasmas being penetrated by slow neutrals of (∼few eV) temperature. The second includes plasmas being penetrated only by fast neutrals having a temperature comparable to that of the ions. The impact of impurities on energy balance is considered through their generation by ion induced desorption of adsorbed oxygen on the first wall and physical and chemical sputtering of carbon. The paper demonstrates self-consistently that the evolution of initial phase of the discharge is intimately linked to the condition of the plasma facing components (PFCs) and the resultant plasma surface interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2013

Effect of volume and surface charges on discharge structure of glow dielectric barrier discharge

The effect of volume and surface charges on the structure of glow dielectric barrier discharge (DBD) has been investigated numerically by using two-dimensional (2D) fluid modeling. The local increase of volume or surface charges induces a kind of activation-inhibition effect, which enhances the local volume discharge and inhibits the discharge in neighborhoods, resulting in non-uniform discharge. The activation-inhibition effect due to the non-uniform volume and/or surface charges depends on the non-uniformity itself and the applied voltage. The activation-inhibition of non-uniform charges has different effects on the volume charges and the accumulated surface charges. The distribution of remaining free charges (seed electrons) in volume at the beginning of voltage pulse plays a key role for the glow DBD structure, resulting in a patterned DBD, when the seed electrons are non-uniform at higher frequency and moderate voltage or uniform DBD, when the seed electrons are uniform at lower frequency or high voltage. The distribution of surface charges is not the determining factor but a result of the formed DBD structure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Monte Carlo implementation of a guiding-center Fokker-Planck kinetic equation

A Monte Carlo method for the collisional guiding-center Fokker-Planck kinetic equation is derived in the five-dimensional guiding-center phase space, where the effects of magnetic drifts due to the background magnetic field nonuniformity are included. It is shown that, in the limit of a homogeneous magnetic field, our guiding-center Monte Carlo collision operator reduces to the guiding-center Monte Carlo Coulomb operator previously derived by Xu and Rosenbluth [Phys. Fluids B 3, 627 (1991)]. Applications of the present work will focus on the collisional transport of energetic ions in complex nonuniform magnetized plasmas in the large mean-free-path (collisionless) limit, where magnetic drifts must be retained.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Propagation of ion acoustic shock waves in negative ion plasmas with nonextensive electrons

Nonlinear ion acoustic shocks (monotonic as well as oscillatory) waves in negative ion plasmas are investigated. The inertialess electron species are assumed to be nonthermal and follow Tsallis distribution. The dissipation in the plasma is considered via kinematic viscosities of both positive and negative ion species. The Korteweg-de Vries Burgers (KdVB) equation is derived using small amplitude reductive perturbation technique and its analytical solution is presented. The effects of variation of density and temperature of negative ions and nonthermal parameter q of electrons on the strength of the shock structures are plotted for illustration. The numerical solutions of KdVB equation using Runge Kutta method are obtained, and transition from oscillatory to monotonic shock structures is also discussed in detail for negative ions nonthermal plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Gyro-induced acceleration of magnetic reconnection

The linear and nonlinear evolution of magnetic reconnection in collisionless high-temperature plasmas with a strong guide field is analyzed on the basis of a two-dimensional gyrofluid model. The linear growth rate of the reconnecting instability is compared to analytical calculations over the whole spectrum of linearly unstable wave numbers. In the strongly unstable regime (large Δ′), the nonlinear evolution of the reconnecting instability is found to undergo two distinctive acceleration phases separated by a stall phase in which the instantaneous growth rate decreases. The first acceleration phase is caused by the formation of strong electric fields close to the X-point due to ion gyration, while the second acceleration phase is driven by the development of an open Petschek-like configuration due to both ion and electron temperature effects. Furthermore, the maximum instantaneous growth rate is found to increase dramatically over its linear value for decreasing diffusion layers. This is a consequence of the fact that the peak instantaneous growth rate becomes weakly dependent on the microscopic plasma parameters if the diffusion region thickness is sufficiently smaller than the equilibrium magnetic field scale length. When this condition is satisfied, the peak reconnection rate asymptotes to a constant value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Formation of current filaments and magnetic field generation in a quantum current-carrying plasma

The nonlinear dynamics of filamentation instability and magnetic field in a current-carrying plasma is investigated in the presence of quantum effects using the quantum hydrodynamic model. A new nonlinear partial differential equation is obtained for the spatiotemporal evolution of the magnetic field in the diffusion regime. This equation is solved by applying the Adomian decomposition method, and then the profiles of magnetic field and electron density are plotted. It is shown that the saturation time of filamentation instability increases and, consequently, the instability growth rate and the magnetic field amplitude decrease in the presence of quantum effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Dynamics of ion acoustic double layers in a magnetized two-population electrons plasma

The obliquely propagating ion acoustic (IA) double-layers are investigated in a magnetized two population electron plasmas. The extended Korteweg–de Vries equation is derived by using the reductive perturbation technique. The effect of obliqueness (l{sub z}) and magnitude of the external magnetic field (ω{sub ci}), as well as the electron number density (β) on the double-layer profile, is studied, and then the ranges of parameters for which the double-layers exist are investigated in detail. We found that the combined effects of l{sub z}, ω{sub ci}, and β significantly modify the basic properties (viz. amplitude and width) of the IA double-layers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Thrust and efficiency model for electron-driven magnetic nozzles

A performance model is presented for magnetic nozzle plasmas driven by electron thermal expansion to investigate how the thrust coefficient and beam divergence efficiency scale with the incoming plasma flow and magnetic field geometry. Using a transformation from cylindrical to magnetic coordinates, an approximate analytical solution is derived to the axisymmetric two-fluid equations for a collisionless plasma flow along an applied magnetic field. This solution yields an expression for the half-width at half-maximum of the plasma density profile in the far-downstream region, from which simple scaling relations for the thrust coefficient and beam divergence efficiency are derived. It is found that the beam divergence efficiency is most sensitive to the density profile of the flow into the nozzle throat, with the highest efficiencies occurring for plasmas concentrated along the nozzle axis. Increasing the expansion ratio of the magnetic field leads to efficiency improvements that are more pronounced for incoming plasmas that are not concentrated along the axis. This implies that the additional magnet required to increase the expansion ratio may be worth the added complexity for plasma sources that exhibit poor confinement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Filament velocity scaling laws for warm ions

The dynamics of filaments or blobs in the scrape-off layer of magnetic fusion devices are studied by magnitude estimates of a comprehensive drift-interchange-Alfvén fluid model. The standard blob models are reproduced in the cold ion case. Even though usually neglected, in the scrape-off layer, the ion temperature can exceed the electron temperature by an order of magnitude. The ion pressure affects the dynamics of filaments amongst others by adding up to the interchange drive and the polarisation current. It is shown how both effects modify the scaling laws for filament velocity in dependence of its size. Simplifications for experimentally relevant limit regimes are given. These are the sheath dissipation, collisional, and electromagnetic regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Negative ion extraction from hydrogen plasma bulk

A two-dimensional particle-in-cell/Monte Carlo collision model has been developed and used to study low electronegative magnetized hydrogen plasma. A configuration characterized by four electrodes is used: the left electrode is biased at V{sub l} = −100 V, the right electrode is grounded, while the upper and lower transversal electrodes are biased at an intermediate voltage V{sub ud} between 0 and −100 V. A constant and homogeneous magnetic field is applied parallel to the lateral (left/right) electrodes. It is shown that in the magnetized case, the bulk plasma potential is close to the transversal electrodes bias inducing then a reversed sheath in front of the right electrode. The potential drop within the reversed sheath is controlled by the transversal electrodes bias allowing extraction of negative ions with a significant reduction of co-extracted electron current. Furthermore, introducing plasma electrodes, between the transversal electrodes and the right electrode, biased with a voltage just above the plasma bulk potential, increases the negative ion extracted current and decreases significantly the co-extracted electron current. The physical mechanism on basis of this phenomenon has been discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Nonlinear evolution of the magnetized Kelvin-Helmholtz instability: From fluid to kinetic modeling

The nonlinear evolution of collisionless plasmas is typically a multi-scale process, where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic (MHD), Hall-MHD, two-fluid, hybrid kinetic, and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligible in the nonlinear regime. This study shows that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the collisionless regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Force-free Jacobian equilibria for Vlasov-Maxwell plasmas

New analytic force-free Vlasov-Maxwell equilibria for thin current sheets are presented. The magnetic flux densities are expressed in terms of Jacobian elliptic functions of one Cartesian spatial coordinate. The magnetic flux densities reduce to previously reported hyperbolic functions in one limit and sinusoidal functions in another limit of the modulus k. A much wider class of nonlinear force-free Vlasov-Maxwell equilibria open expanded possibilities for modeling of solar system, astrophysical and laboratory plasmas. Modified Maxwellian distribution functions are determined explicitly in terms of Jacobian elliptic functions. Conditions for double peaked distribution functions that could be unstable are developed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

A verification of the gyrokinetic microstability codes GEM, GYRO, and GS2

A previous publication [R. V. Bravenec et al., Phys. Plasmas 18, 122505 (2011)] presented favorable comparisons of linear frequencies and nonlinear fluxes from the Eulerian gyrokinetic codes gyro[J. Candy and R. E. Waltz, J. Comput. Phys. 186, 545 (2003)] and gs2[W. Dorland et al., Phys. Rev. Lett. 85, 5579 (2000)]. The motivation was to verify the codes, i.e., demonstrate that they correctly solve the gyrokinetic-Maxwell equations. The premise was that it is highly unlikely for both codes to yield the same incorrect results. In this work, we add the Lagrangian particle-in-cell code gem[Y. Chen and S. Parker, J. Comput. Phys. 220, 839 (2007)] to the comparisons, not simply to add another code, but also to demonstrate that the codes' algorithms do not matter. We find good agreement of gem with gyro and gs2 for the plasma conditions considered earlier, thus establishing confidence that the codes are verified and that ongoing validation efforts for these plasma parameters are warranted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

A phenomenological model on the kink mode threshold varying with the inclination of sheath boundary

In nature and many laboratory plasmas, a magnetic flux tube threaded by current or a flux rope has a footpoint at a boundary. The current driven kink mode is one of the fundamental ideal magnetohydrodynamic instabilities in plasmas. It has an instability threshold that has been found to strongly depend on boundary conditions (BCs). We provide a theoretical model to explain the transition of this threshold dependence between nonline tied and line tied boundary conditions. We evaluate model parameters using experimentally measured plasma data, explicitly verify several kink eigenfunctions, and validate the model predictions for boundary conditions BCs that span the range between NLT and LT BCs. Based on this model, one could estimate the kink threshold given knowledge of the displacement of a flux rope end, or conversely estimate flux rope end motion based on knowledge of it kink stability threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

On the atomic line profiles in high pressure plasmas

In a previous contribution to this journal [H. P. Stormberg, J. Appl. Phys. 51(4), 1963 (1980)], Stormberg presented an analytical expression for the convolution of Lorentz and Levy line profiles, which models atomic radiative transitions in high pressure plasmas. Unfortunately, the derivations are flawed with errors and the final expression, while correct, is accompanied by misguiding comments about the meaning of the symbols used therein, in particular the “complex error function.” In this paper, we discuss the broadening mechanisms that give rise to Stormberg's model and present a correct derivation of his final result. We will also provide an alternative expression, based on the Faddeeva function, which has decisive computational advantages and emphasizes the real-valuedness of the result. The MATLAB/Octave scripts of our implementation have been made available on the publisher's website for future reference.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Penetration of lower hybrid current drive waves in tokamaks

Lower hybrid (LH) ray propagation in toroidal plasma is shown to be controlled by combination of the azimuthal spectrum launched by the antenna, the poloidal variation of the magnetic field, and the scattering of the waves by the drift wave fluctuations. The width of the poloidal and radial radio frequency wave spectrum increases rapidly as the rays penetrate into higher density and scatter from the drift waves. The electron temperature gradient (ETG) spectrum is particularly effective in scattering the LH waves due to its comparable wavelengths and phase velocities. ETG turbulence is also driven by the radial gradient of the electron current profile giving rise to an anomalous viscosity spreading the LH driven plasma currents. The LH wave scattering is derived from a Fokker-Planck equation for the distribution of the ray trajectories with diffusivities derived from the drift wave fluctuations. The condition for chaotic diffusion for the rays is derived. The evolution of the poloidal and radial mode number spectrum of the lower hybrid waves are both on the antenna spectrum and the spectrum of the drift waves. Antennas launching higher poloidal mode number spectra drive off-axis current density profiles producing negative central shear [RS] plasmas with improved thermal confinement from ETG transport. Core plasma current drive requires antennas with low azimuthal mode spectra peaked at m = 0 azimuthal mode numbers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Modeling target bulk heating resulting from ultra-intense short pulse laser irradiation of solid density targets

Isochoric heating of solid-density matter up to a few tens of eV is of interest for investigating astrophysical or inertial fusion scenarios. Such ultra-fast heating can be achieved via the energy deposition of short-pulse laser generated electrons. Here, we report on experimental measurements of this process by means of time- and space-resolved optical interferometry. Our results are found in reasonable agreement with a simple numerical model of fast electron-induced heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗