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

Stability of compressible reduced magnetohydrodynamic equilibria-Analogy with magnetorotational instability

Stability analyses for equilibria of the compressible reduced magnetohydrodynamics (CRMHD) model are carried out by means of the Energy-Casimir (EC) method. Stability results are compared with those obtained for ideal magnetohydrodynamics (MHD) from the classical {delta}W criterion. An identification of the terms in the second variation of the free energy functional for CRMHD with those of {delta}W is made: two destabilizing effects present for CRMHD turn out to correspond to the kink and interchange instabilities in usual MHD, while the stabilizing roles of field line bending and compressibility are also identified in the reduced model. Also, using the EC method, stability conditions in the presence of toroidal flow are obtained. A formal analogy between CRMHD and a reduced incompressible model for magnetized rotating disks, due to Julien and Knobloch [EAS Pub. Series, 21, 81 (2006)], is discovered. In light of this analogy, energy stability analysis shows that the condition for magnetorotational instability (MRI) for the latter model corresponds to the condition for interchange instability in CRMHD, with the Coriolis term and shear velocity playing the roles of the curvature term and pressure gradient, respectively. Using the EC method, stability conditions for the rotating disk model, for a large class of equilibria with possible non-uniform magnetic fields, are obtained. In particular, this shows it is possible for the MRI system to undergo, in addition to the MRI, another instability that is analogous to the kink instability. For vanishing magnetic field, the Rayleigh hydrodynamical stability condition is recovered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Relationship between the ion drag and electric forces in dense dust clouds

It is shown that the relation between the ion drag and electric forces is strongly dependent on the dust number density in complex plasmas. The effect of the particle size and discharge parameters on the force balance is investigated. Examples are given for realistic complex plasma parameters and comparison with microgravity experiments is presented.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Central peaking of magnetized gas discharges

Partially ionized gas discharges used in industry are often driven by radiofrequency (rf) power applied at the periphery of a cylinder. It is found that the plasma density n is usually flat or peaked on axis even if the skin depth of the rf field is thin compared with the chamber radius a. Previous attempts at explaining this did not account for the finite length of the discharge and the boundary conditions at the endplates. A simple 1D model is used to focus on the basic mechanism: the short-circuit effect. It is found that a strong electric field (E-field) scaled to electron temperature T{sub e}, drives the ions inward. The resulting density profile is peaked on axis and has a shape independent of pressure or discharge radius. This “universal” profile is not affected by a dc magnetic field (B-field) as long as the ion Larmor radius is larger than a.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Laser-plasma source parameters for Kr, Gd, and Tb ions at 6.6 nm

There is increasing interest in extreme-ultraviolet (EUV) laser-based lamps for sub-10-nm lithography operating in the region of 6.6 nm. A collisional-radiative model is developed as a post-processor of a hydrodynamic code to investigate emission from resonance lines in Kr, Gd, and Tb ions under conditions typical for mass-limited EUV sources. The analysis reveals that maximum conversion efficiencies of Kr occur at 5 Multiplication-Sign 10{sup 10}W/cm{sup 2}, while for Gd and Tb it was Asymptotically-Equal-To 0.9%/2{pi}sr for laser intensities of (2-5) Multiplication-Sign 10{sup 12}W/cm{sup 2}.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Characterization and performance of a field aligned ion cyclotron range of frequency antenna in Alcator C-Mod

Ion cyclotron range of frequency (ICRF) heating is expected to provide auxiliary heating for ITER and future fusion reactors where high Z metallic plasma facing components (PFCs) are being considered. Impurity contamination linked to ICRF antenna operation remains a major challenge particularly for devices with high Z metallic PFCs. Here, we report on an experimental investigation to test whether a field aligned (FA) antenna can reduce impurity contamination and impurity sources. We compare the modification of the scrape of layer (SOL) plasma potential of the FA antenna to a conventional, toroidally aligned (TA) antenna, in order to explore the underlying physics governing impurity contamination linked to ICRF heating. The FA antenna is a 4-strap ICRF antenna where the current straps and antenna enclosure sides are perpendicular to the total magnetic field while the Faraday screen rods are parallel to the total magnetic field. In principle, alignment with respect to the total magnetic field minimizes integrated E|| (electric field along a magnetic field line) via symmetry. A finite element method RF antenna model coupled to a cold plasma model verifies that the integrated E|| should be reduced for all antenna phases. Monopole phasing in particular is expected to have the lowest integrated E||. Consistent with expectations, we observed that the impurity contamination and impurity source at the FA antenna are reduced compared to the TA antenna. In both L and H-mode discharges, the radiated power is 20%–30% lower for a FA-antenna heated discharge than a discharge heated with the TA-antennas. However, inconsistent with expectations, we observe RF induced plasma potentials (via gas-puff imaging and emissive probes to be nearly identical for FA and TA antennas when operated in dipole phasing). Moreover, the highest levels of RF-induced plasma potentials are observed using monopole phasing with the FA antenna. Thus, while impurity contamination and sources are indeed reduced with the FA antenna configuration, the mechanism determining the SOL plasma potential in the presence of ICRF and its impact on impurity contamination and sources remains to be understood.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Bifurcations of nonlinear ion acoustic travelling waves in the frame of a Zakharov-Kuznetsov equation in magnetized plasma with a kappa distributed electron

Bifurcations of nonlinear propagation of ion acoustic waves (IAWs) in a magnetized plasma whose constituents are cold ions and kappa distributed electron are investigated using a two component plasma model. The standard reductive perturbation technique is used to derive the Zakharov-Kuznetsov (ZK) equation for IAWs. By using the bifurcation theory of planar dynamical systems to this ZK equation, the existence of solitary wave solutions and periodic travelling wave solutions is established. All exact explicit solutions of these travelling waves are determined. The results may have relevance in dense space plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Charge-to-mass-ratio-dependent ion heating during magnetic reconnection in the MST RFP

Temperature evolution during magnetic reconnection has been spectroscopically measured for various ion species in a toroidal magnetized plasma. Measurements are made predominantly in the direction parallel to the equilibrium magnetic field. It is found that the increase in parallel ion temperature during magnetic reconnection events increases with the charge-to-mass ratio of the ion species. This trend can be understood if the heating mechanism is anisotropic, favoring heating in the perpendicular degree of freedom, with collisional relaxation of multiple ion species. The charge-to-mass ratio trend for the parallel temperature derives from collisional isotropization. This result emphasizes that collisional isotropization and energy transfer must be carefully modeled when analyzing ion heating measurements and comparing to theoretical predictions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Temporal and spatial structure of a runaway electron beam in air at atmospheric pressure

The time- and spatial structure of a runaway electron beam generated in air at atmospheric pressure by a high-voltage pulse with a rise time of {approx}300 ps is studied experimentally and numerically. It is obtained that the duration of the runaway electron current is a few tens of picoseconds, and it can consist of two or many peaks. It is shown that the many-peak temporal structure of the beam is caused by the non-simultaneous appearance of several emission centers on the cathode edge.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Modelling of the beam-plasma interaction in a strongly inhomogeneous plasma

The interaction of a beam propagating in a inhomogeneous solar wind plasma is considered. The properties of the plasma waves are described by the one dimensional Zakharov's equation, the beam is modelled by means of particles moving in the electric fields of the Langmuir waves. We take into account the presence of high level density fluctuations that are known to be present in the solar wind. It is shown that when the level of density fluctuations is low, {delta}n/n{sub 0} Much-Less-Than 3k{sup 2}{lambda}{sub d}{sup 2}, the regime of beam relaxation is very similar to that occurring in a homogeneous plasma and can be described by the quasilinear (QL) equations. In this case, the relaxation length is very short and corresponds to that obtained using the QL approximation. On the contrary, when the level of density fluctuations overcomes some limit, i.e. {delta}n/n{sub 0}{>=}3k{sup 2}{lambda}{sub d}{sup 2}, where {delta}n/n{sub 0} Much-Less-Than 1, the plasma inhomogeneities crucially influence the process of relaxation. First, the linear wave growth becomes localized and clearly identifiable wave packets/clumps dominate the wave spectrum; this is associated with the kinematic properties of the waves' propagation and the wave-particle resonant interactions. Most of the wave packets grow in the regions of density gradients. The second important feature revealed by the beam relaxation consists in the generation of a tail of accelerated electrons with velocities V > V{sub b} exceeding the beam drift. The beam widens in both directions, toward lower velocities as well as higher velocities; the density of the accelerated electrons can reach more than 10 Division-Sign 20 % of the beam density and the energy flux carried by this population can become as large as 40% of the initial energy flux of the beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Comparison of kinetic and extended magnetohydrodynamics computational models for the linear ion temperature gradient instability in slab geometry

We perform linear stability studies of the ion temperature gradient (ITG) instability in unsheared slab geometry using kinetic and extended magnetohydrodynamics (MHD) models, in the regime k{sub ∥}/k{sub ⊥}≪1. The ITG is a parallel (to B) sound wave that may be destabilized by finite ion Larmor radius (FLR) effects in the presence of a gradient in the equilibrium ion temperature. The ITG is stable in both ideal and resistive MHD; for a given temperature scale length L{sub Ti0}, instability requires that either k{sub ⊥}ρ{sub i} or ρ{sub i}/L{sub Ti0} be sufficiently large. Kinetic models capture FLR effects to all orders in either parameter. In the extended MHD model, these effects are captured only to lowest order by means of the Braginskii ion gyro-viscous stress tensor and the ion diamagnetic heat flux. We present the linear electrostatic dispersion relations for the ITG for both kinetic Vlasov and extended MHD (two-fluid) models in the local approximation. In the low frequency fluid regime, these reduce to the same cubic equation for the complex eigenvalue ω=ω{sub r}+iγ. An explicit solution is derived for the growth rate and real frequency in this regime. These are found to depend on a single non-dimensional parameter. We also compute the eigenvalues and the eigenfunctions with the extended MHD code NIMROD, and a hybrid kinetic δf code that assumes six-dimensional Vlasov ions and isothermal fluid electrons, as functions of k{sub ⊥}ρ{sub i} and ρ{sub i}/L{sub Ti0} using a spatially dependent equilibrium. These solutions are compared with each other, and with the predictions of the local kinetic and fluid dispersion relations. Kinetic and fluid calculations agree well at and near the marginal stability point, but diverge as k{sub ⊥}ρ{sub i} or ρ{sub i}/L{sub Ti0} increases. There is good qualitative agreement between the models for the shape of the unstable global eigenfunction for L{sub Ti0}/ρ{sub i}=30 and 20. The results quantify how far fluid calculations can be extended accurately into the kinetic regime. We conclude that for the linear ITG problem in slab geometry with unsheared magnetic field when k{sub ∥}/k{sub ⊥}≪1, the extended MHD model may be a reliable physical model for this problem when ρ{sub i}/L{sub Ti0}<10{sup −2} and k{sub ⊥}ρ{sub i}<0.2.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Ideal magnetohydrodynamic stability of configurations without nested flux surfaces

Existing numerical tools for calculating the MHD stability of magnetically confined plasmas generally assume the existence of nested flux surfaces. These tools are therefore not immediately applicable to configurations with magnetic islands or regions with an ergodic magnetic field. However, in practice, these islands or ergodic regions are often small, and their effect on MHD stability can then be evaluated using a perturbation theory developed in the present paper. This procedure allows the effect of the broken magnetic topology on the stability of each eigenmode to be calculated without requiring any knowledge about the perturbed eigenfunctions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

The transition mechanism from a symmetric single period discharge to a period-doubling discharge in atmospheric helium dielectric-barrier discharge

Period-doubling and chaos phenomenon have been frequently observed in atmospheric-pressure dielectric-barrier discharges. However, how a normal single period discharge bifurcates into period-doubling state is still unclear. In this paper, by changing the driving frequency, we study numerically the transition mechanisms from a normal single period discharge to a period-doubling state using a one-dimensional self-consistent fluid model. The results show that before a discharge bifurcates into a period-doubling state, it first deviates from its normal operation and transforms into an asymmetric single period discharge mode. Then the weaker discharge in this asymmetric discharge will be enhanced gradually with increasing of the frequency until it makes the subsequent discharge weaken and results in the discharge entering a period-doubling state. In the whole transition process, the spatial distribution of the charged particle density and the electric field plays a definitive role. The conclusions are further confirmed by changing the gap width and the amplitude of the applied voltage.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Electrical studies and plasma characterization of an atmospheric pressure plasma jet operated at low frequency

Low-temperature, high-pressure plasma jets have an extensive use in medical and biological applications. Much work has been devoted to study these applications while comparatively fewer studies appear to be directed to the discharge itself. In this work, in order to better understand the kind of electrical discharge and the plasma states existing in those devices, a study of the electrical characteristics of a typical plasma jet, operated at atmospheric pressure, using either air or argon, is reported. It is found that the experimentally determined electrical characteristics are consistent with the model of a thermal arc discharge, with a highly collisional cathode sheet. The only exception is the case of argon at the smallest electrode separation studied, around 1 mm in which case the discharge is better modeled as either a non-thermal arc or a high-pressure glow. Also, variations of the electrical behavior at different gas flow rates are interpreted, consistently with the arc model, in terms of the development of fluid turbulence in the external jet.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Spatio-temporal evolution of the H → L back transition

Since ITER will operate close to threshold and with limited control, the H → L back transition is a topic important for machine operations as well as physics. Using a reduced mesoscale model [Miki et al., Phys. Plasmas 19, 092306 (2012)], we investigate ELM-free H → L back transition dynamics in order to isolate transport physics effects. Model studies indicate that turbulence spreading is the key process which triggers the back transition. The transition involves a feedback loop linking turbulence and profiles. The I-phase appears during the back transition following a slow power ramp down, while fast ramp-downs reveal a single burst of zonal flow during the back transition. The I-phase nucleates at the pedestal shoulder, as this is the site of the residual turbulence in H-mode. Hysteresis in the profile gradient scale length is characterized by the Nusselt number, where Nu=χ{sub i,turb}/χ{sub i,neo}. Relative hysteresis of temperature gradient vs density gradient is sensitive to the pedestal Prandtl number, where Pr{sub ped}=D{sub ped}/χ{sub i,neo}. We expect the H-mode to be somewhat more resilient in density than in temperature.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Self-similar expansion of a warm dense plasma

The properties of an expanding plasma composed of degenerate electron fluid and non-degenerate ions are studied. For our purposes, we use fluid equations for ions together with the electron momentum equation that include quantum forces (e.g., the quantum statistical pressure, forces due to the electron-exchange and electron correlations effects) and the quasi-neutrality condition. The governing equation is written in a tractable form by using a self-similar transformation. Numerical results for typical beryllium plasma parameters revealed that, during the expansion, the ion acoustic speed decreases for both isothermal and adiabatic ion pressure. When compared with classical hydrodynamic plasma expansion model, the electrons and ions are found to initially escape faster in vacuum creating thus an intense electric field that accelerates most of the particles into the vacuum ahead of the plasma expansion. The relevancy of the present model to beryllium plasma produced by a femto-second laser is highlighted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Automation of the guiding center expansion

We report on the use of the recently developed Mathematica package VEST (Vector Einstein Summation Tools) to automatically derive the guiding center transformation. Our Mathematica code employs a recursive procedure to derive the transformation order-by-order. This procedure has several novel features. (1) It is designed to allow the user to easily explore the guiding center transformation's numerous non-unique forms or representations. (2) The procedure proceeds entirely in cartesian position and velocity coordinates, thereby producing manifestly gyrogauge invariant results; the commonly used perpendicular unit vector fields e{sub 1},e{sub 2} are never even introduced. (3) It is easy to apply in the derivation of higher-order contributions to the guiding center transformation without fear of human error. Our code therefore stands as a useful tool for exploring subtle issues related to the physics of toroidal momentum conservation in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Kelvin-Helmholtz instability in non-Newtonian complex plasma

The Kelvin-Helmholtz (KH) instability is studied in a non-Newtonian dusty plasma with an experimentally verified model [Phys. Rev. Lett. 98, 145003 (2007)] of shear flow rate dependent viscosity. The shear flow profile used here is a parabolic type bounded flow. Both the shear thinning and shear thickening properties are investigated in compressible as well as incompressible limits using a linear stability analysis. Like the stabilizing effect of compressibility on the KH instability, the non-Newtonian effect in shear thickening regime could also suppress the instability but on the contrary, shear thinning property enhances it. A detailed study is reported on the role of non-Newtonian effect on KH instability with conventional dust fluid equations using standard eigenvalue analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗