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,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Implementation and assessment of an extended hydrogenic molecular model in UEDGE

Abstract Two sets of pure deuterium plasmas are simulated using UEDGE, one including atoms and molecules and the other including atoms only as separate fluid species. Simulation results in one and two dimensions are reported to assess the role of molecules in tokamak plasmas. It is shown that thermal coupling of the molecules to the plasma can be relevant for the onset of detachment under highly collisional conditions in simple geometries, but play a modest role at low collisionality. Ion‐molecule equipartition presents an additional ion energy loss channel, dissipating ion energy that would otherwise heat the electrons by equipartition. The resulting reduced electron temperature increases the ion‐electron recombination by an order of magnitude as the plasma density is increased, yielding earlier and deeper detachment.

Holm, A.↗
DOE OSTI2013

Effects of the non-uniform initial environment and the guide field on the plasmoid instability

Effects of non-uniform initial mass density and temperature on the plasmoid instability are studied via 2.5-dimensional resistive magnetohydrodynamic (MHD) simulations. Our results indicate that the development of the plasmoid instability is apparently prevented when the initial plasma density at the center of the current sheet is higher than that in the upstream region. As a result, the higher the plasma density at the center and the lower the plasma β in the upstream region, the higher the critical Lundquist number needed for triggering secondary instabilities. When β=0.2, the critical Lundquist number is higher than 10{sup 4}. For the same Lundquist number, the magnetic reconnection rate is lower for the lower plasma β case. Oppositely, when the initial mass density is uniform and the Lundquist number is low, the magnetic reconnection rate turns out to be higher for the lower plasma β case. For the high Lundquist number case (>10{sup 4}) with uniform initial mass density, the magnetic reconnection is not affected by the initial plasma β and the temperature distribution. Our results indicate that the guide field has a limited impact on the plasmoid instability in resistive MHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

One-dimensional full wave simulation on XB mode conversion in electron cyclotron heating

The XB mode conversion in electron cyclotron resonance frequency heating has been studied in detail through 1D full wave simulation. The field pattern depends on the density scale length, and the wave absorption near upper hybrid resonance is maximized beyond the R(X) mode cutoff density for optimized density scale length. The simulated mode conversion efficiency has been compared with that of an analytic formula, showing good agreements except for the phase dependent term of the X wave. The mode conversion efficiency is calculated for oblique injections as well, and it is found that the efficiency decreases as the injection angles increases. Short magnetic field scale length is confirmed to relax the short density scale length condition maximizing the XB mode conversion efficiency. Finally, the simulation code is used to analyze the mode conversion and power absorption of a pre-ionization plasma in versatile experiment spherical torus.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Generalization of the Child-Langmuir law to the alternate extraction of positive and negative ions

Using a combined analytical and simulation approach, we investigate positive and negative ion extraction between two electrodes from an ion-ion plasma source. With a square voltage waveform applied to the electrodes, we obtain approximate analytical solutions for the time-averaged extracted current densities, which are given simply by: J{sub p}{sup ac}=[α−fL√((M{sub p})/(q{sub p}V{sub 0}) )]J{sub p}{sup dc}, and J{sub n}{sup ac}=[(1−α)−fL√((M{sub n})/(q{sub n}V{sub 0}) )]J{sub n}{sup dc}, where J{sup ac} is the time-averaged current density, α is the square waveform duty cycle, f is the frequency, L is the electrode gap length, M is the ion mass, q is the ion charge, V{sub 0} is the applied voltage amplitude, J{sup dc} is the dc extracted current density, and the subscripts p and n refer to positive and negative ions, respectively. In particular, if J{sup dc} is the dc space-charge limited current density, then these equations describe the square waveform generalization of the Child-Langmuir law.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Generation of laser plasma bunches with a high efficiency of energy concentration for laboratory simulation of collisionless shock waves in magnetised cosmic plasma

We present the results of first experiments on the formation of collisionless shock waves (CSWs) in background plasma by injecting laser plasma bunches transverse to the magnetic field (as a piston) with a maximum energy up to 100 J per unit of solid angle and with a high enough degree of ion magnetisation. With this aim in view, on a unique KI-1 facility at the Institute of Laser Physics, Siberian Branch of the Russian Academy of Sciences (ILP), a plastic (polyethylene) target irradiated by a CO{sub 2} laser in the most energy-efficient regime (near the plasma formation threshold) and a highly ionised hydrogen plasma with a high concentration in a large volume (not less than 1 m{sup 3}) have been employed. As a result of model experiments performed on the basis of a model of collisionless interaction of plasma flows, developed at the VNIIEF and being adequate to the problem under consideration, not only an intensive, background-induced, deceleration of a super-Alfven laser plasma flow, but also the formation in that flow of a strong perturbation having the properties of a subcritical CSW and propagating transverse to the magnetic field, have been first registered in the laboratory conditions. (extreme laser radiation: physics and fundamental applications)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Experiment and simulation of novel liquid crystal plasma mirrors for high contrast, intense laser pulses

Here, we describe the first demonstration of plasma mirrors made using freely suspended, ultra-thin films formed dynamically and in-situ. We also present novel particle-in-cell simulations that for the first time incorporate multiphoton ionization and dielectric models that are necessary for describing plasma mirrors. Dielectric plasma mirrors are a crucial component for high intensity laser applications such as ion acceleration and solid target high harmonic generation because they greatly improve pulse contrast. We use the liquid crystal 8CB and introduce an innovative dynamic film formation device that can tune the film thickness so that it acts as its own antireflection coating. Films can be formed at a prolonged, high repetition rate without the need for subsequent realignment. High intensity reflectance above 75% and low-field reflectance below 0.2% are demonstrated, as well as initial ion acceleration experimental results that demonstrate increased ion energy and yield on shots cleaned with these plasma mirrors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Plasma physics in strong-field regimes: Theories and simulations

In strong electromagnetic fields, unique plasma phenomena and applications emerge whose description requires recently developed theories and simulations [Y. Shi, Ph.D. thesis, Princeton University (2018)]. In the classical regime, to quantify effects of strong magnetic fields on three-wave interactions, a convenient formula is derived by solving the fluid model to the second order in general geometry. As an application, magnetic resonances are exploited to mediate laser pulse compression, whereby higher intensity pulses can be produced in wider frequency ranges, as confirmed by particle-in-cell simulations. In even stronger fields, relativistic-quantum effects become important, and a plasma model based on scalar quantum electrodynamics (QED) is developed which unveils observable corrections to Faraday rotation and cyclotron absorption in strongly magnetized plasmas. Beyond the perturbative regime, lattice QED is extended as a numerical tool for plasma physics, using which the transition from wakefield acceleration to electron-positron pair production is captured when laser intensity exceeds the Schwinger threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Global current profile effects on the evolution and saturation of magnetic islands

The nonlinear evolution of magnetic islands is investigated by means of numerical simulations. The growth and saturation of the island are shown to depend not only on asymptotic tearing mode stability parameter {Delta}{sup Prime} but also on the initial global current profile. Presence of the external current sheets leads to the formation of different island state for the same value of {Delta}{sup Prime }. It is found that the flow vorticity generated by the interacting current sheets is an important feature of nonlocal interactions and nonlocal effects in magnetic islands dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Improving cryogenic deuterium–tritium implosion performance on OMEGA

A flexible direct-drive target platform is used to implode cryogenic deuterium–tritium (DT) capsules on the OMEGA laser [Boehly et al., Opt. Commun. 133, 495 (1997)]. The goal of these experiments is to demonstrate ignition hydrodynamically equivalent performance where the laser drive intensity, the implosion velocity, the fuel adiabat, and the in-flight aspect ratio (IFAR) are the same as those for a 1.5-MJ target [Goncharov et al., Phys. Rev. Lett. 104, 165001 (2010)] designed to ignite on the National Ignition Facility [Hogan et al., Nucl. Fusion 41, 567 (2001)]. The results from a series of 29 cryogenic DT implosions are presented. The implosions were designed to span a broad region of design space to study target performance as a function of shell stability (adiabat) and implosion velocity. Ablation-front perturbation growth appears to limit target performance at high implosion velocities. Target outer-surface defects associated with contaminant gases in the DT fuel are identified as the dominant perturbation source at the ablation surface; performance degradation is confirmed by 2D hydrodynamic simulations that include these defects. A trend in the value of the Lawson criterion [Betti et al., Phys. Plasmas 17, 058102 (2010)] for each of the implosions in adiabat–IFAR space suggests the existence of a stability boundary that leads to ablator mixing into the hot spot for the most ignition-equivalent designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

Sub-picosecond and nanometer scale dynamics of aluminum target surface heated by ultrashort laser pulse

The possibility of measuring surface hydrodynamics at picosecond and nanometer scales by using a time resolved optical interfero-polarimeter is demonstrated. This method is applied to observe the surface position of aluminum target heated by ultrashort laser pulse. Numerical simulations of the surface relaxation during a few tens of picoseconds after irradiation are in very good agreement with the measurements. Results show the important contribution of oxide layer in the surface dynamic evolution and the great accuracy of our experimental method.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Magnetic helicity signature produced by cross-field 2D turbulence

Hybrid numerical simulations of freely decaying 2D turbulence are presented. The background magnetic field is perpendicular to the simulation plane, which eliminates linear kinetic Alfven waves from the system. The normalized magnetic helicity of the initial large-scale fluctuations is zero, while the normalized cross-helicity is not. As the turbulence evolves, it develops nonzero magnetic helicity at smaller scales, in the proton kinetic range. In the quasi-steady state of evolution, the magnetic helicity spectrum has a peak consistent with the solar wind observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Finite Mach number spherical shock wave, application to shock ignition

A converging and diverging spherical shock wave with a finite initial Mach number M{sub s0} is described by using a perturbative approach over a small parameter M{sub s}{sup −2}. The zeroth order solution is the Guderley's self-similar solution. The first order correction to this solution accounts for the effects of the shock strength. Whereas it was constant in the Guderley's asymptotic solution, the amplification factor of the finite amplitude shock Λ(t)∝dU{sub s}/dR{sub s} now varies in time. The coefficients present in its series form are iteratively calculated so that the solution does not undergo any singular behavior apart from the position of the shock. The analytical form of the corrected solution in the vicinity of singular points provides a better physical understanding of the finite shock Mach number effects. The correction affects mainly the flow density and the pressure after the shock rebound. In application to the shock ignition scheme, it is shown that the ignition criterion is modified by more than 20% if the fuel pressure prior to the final shock is taken into account. A good agreement is obtained with hydrodynamic simulations using a Lagrangian code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Study of detached plasma profile in the divertor simulation experimental module of tandem mirror GAMMA 10/PDX

A divertor simulation experimental module (D-module) in the tandem mirror GAMMA 10/PDX was used for the study of plasma detachment. In previous studies, it was difficult to measure far-upstream plasma parameters in the D-module, and only electrostatic probes on the target plate were used to perform electron temperature and density measurements. Furthermore, to study the detached plasma structure, a Thomson scattering (TS) system and a microwave interferometer system have been installed to measure the inside plasma parameters of the D-module, and a movable electrostatic probe has been placed at the inlet of the D-module to measure the inlet plasma density and temperature. The TS system in the central cell observed the electron temperature and density of the core plasma simultaneously. These measurements revealed the entire density and temperature structure from the core plasma to the divertor plate. The line average electron density measured by the microwave interferometer showed a rollover behavior during detachment. The results indicated that the ionization region was located around the center of the D-module, and it appears to move upstream along the axis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2013

DISTRIBUTION OF PLASMOIDS IN POST-CORONAL MASS EJECTION CURRENT SHEETS

Recently, the fragmentation of a current sheet in the high-Lundquist-number regime caused by the plasmoid instability has been proposed as a possible mechanism for fast reconnection. In this work, we investigate this scenario by comparing the distribution of plasmoids obtained from Large Angle and Spectrometric Coronagraph (LASCO) observational data of a coronal mass ejection event with a resistive magnetohydrodynamic simulation of a similar event. The LASCO/C2 data are analyzed using visual inspection, whereas the numerical data are analyzed using both visual inspection and a more precise topological method. Contrasting the observational data with numerical data analyzed with both methods, we identify a major limitation of the visual inspection method, due to the difficulty in resolving smaller plasmoids. This result raises questions about reports of log-normal distributions of plasmoids and other coherent features in the recent literature. Based on nonlinear scaling relations of the plasmoid instability, we infer a lower bound on the current sheet width, assuming the underlying mechanism of current sheet broadening is resistive diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Orientational order and formation of topological defects in two-dimensional systems

The formation dynamics of topological defects and the orientational order in nonideal systems of particles interacting with a screened Coulomb potential is numerically simulated. Calculations are performed in a broad range of parameters corresponding to experimental conditions in a laboratory dusty plasma. A relation is obtained between the number of defects, the shape of the orientational correlation function, and the coupling parameter of the system. New approximations are proposed for the orientational correlation function in the liquid and hexatic phases of structures analyzed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

Preplasma effects on the generation of high-energy protons in ultraintense laser interaction with foil targets

It is shown that the intense quasistatic electric and magnetic fields self-generated near the axis of the laser-driven channel in an appropriately profiled preplasma during ultraintense laser interaction with a thin target can create dense relativistic electron bunches. The latter easily penetrate through the target and can greatly enhance the sheath field at the rear, resulting in significant increase in the laser-to-ion energy conversion efficiency and the maximum energy of the target normal sheath accelerated ions. Particle-in-cell simulations show that with a hydrogen targets a proton beam of peak energy ∼38 MeV and energy conversion efficiency ≥6.5% can be produced by a linearly polarized 5 × 10{sup 19} W/cm{sup 2} laser. An analytical model is also proposed and its results agree well with those of the simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

A simulation study of radial expansion of an electron beam injected into an ionospheric plasma

Injections of nonrelativistic electron beams from a finite equipotential conductor into an ionospheric plasma have been simulated using a two-dimensional electrostatic particle code. The purpose of the study is to survey the simulation parameters for understanding the dependence of beam radius on physical variables. The conductor is charged to a high potential when the background plasma density is less than the beam density. Beam electrons attracted by the charged conductor are decelerated to zero velocity near the stagnation point, which is at a few Debye lengths from the conductor. The simulations suggest that the beam electrons at the stagnation point receive a large transverse kick and the beam expands radially thereafter. The buildup of beam electrons at the stagnation point produces a large electrostatic force responsible for the transverse kick. However, for the weak charging cases where the background plasma density is larger than the beam density, the radial expansion mechanism is different; the beam plasma instability is found to be responsible for the radial expansion. The simulations show that the electron beam radius for high spacecraft charging cases is of the order of the beam gyroradius, defined as the beam velocity divided by the gyrofrequency. In the weak charging cases, the beam radius is only a fraction of the beam gyroradius. The parameter survey indicates that the beam radius increases with beam density and decreases with magnetic field and beam velocity. The beam radius normalized by the beam gyroradius is found to scale according to the ratio of the beam electron Debye length to the ambient electron Debye length. The parameter dependence deduced would be useful for interpreting the beam radius and beam density of electron beam injection experiments conducted from rockets and the space shuttle.

Koga, J.↗