Search NASA⌕ Search

SEARCH · Search NASA

Results for “PLASMA DYNAMICS”

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 721 records · Page 40
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Laser Plasma Interaction in the MG Magnetic Field (Final Report)

Strong magnetic fields change properties of plasma. Employment of the magnetic fields can improve conditions for inertial confinement fusion. Compression and heating of the magnetized plasma is a base of the MagLIF pulsed power approach to the controlled fusion. Study of expansion and heating of plasma in the external magnetic fields, development of plasma instabilities, parametric and other effects are important for basic physics and applications. The main goal of the research for the grant was experimental investigation of laser produced plasma in the well characterized and controlled MG magnetic fields. We developed a platform for investigation of high intensity laser-plasma interaction (LPI) in the MG fields using a university-scale pulsed power machine. A Zebra machine at UNR generated transverse magnetic fields of 3 MG on the surface of the rod load and longitudinal fields of 1.4 MG in the coil loads. The magnetic fields were characterized with two-color Faraday rotation diagnostics. The pulsed power machine provided almost static magnetic fields for LPI compared to the much faster laser driven magnetic fields. We found that eddy currents in the rising magnetic field generate plasma on the metal targets. Eddy currents did not impact Si and CH dielectric targets. A Leopard laser with intensity in the focus up to 10 19 W/cm 2 produced plasma in the MG magnetic fields. The modified reflected Faraday diagnostics was used to study laser driven B-field at the MTW laser at Laboratory for Laser Energetics, UR. We demonstrated a slow dependence of the magnetic field in the coil on the applied laser intensity. The magnetic field increased by a factor of 2 if laser intensity increased by a factor of 30. This scaling is important for comparison of experiments at different laser facilities. Expansion of the plasma in the 2-3 MG external B-field was studied. Laser-produced plasma in the azimuthal magnetic field took the unique form of a thin disc expanding radially with a velocity of 250 km/s and confined in the vertical direction. A HYDRA MHD program at the University of Rochester was applied for simulations of LPI in the magnetic field B >3 MG. Simulations for the plasma disc parameters and the expansion dynamics were in a good agreement with experiments. Generation of narrow dense plasma jets in the longitudinal magnetic field of 0.6-0.8 MG was studied. Narrow plasma jets reached a length of 3-4 mm with the electron plasma density of (0.2-1.2)x10 20 cm -3 . A jet tip propagated with the velocity of 160-200 km/s. MHD simulations showed a good agreement of the dynamics of the formation of plasma jets with experiments. These jets are relevant to astrophysical jets. Simulations of K-shell spectra of Si plasma was performed for LPI in the B-field. The laser prepulse generated plasma near the laser target. A PrismSPECT model with a MeV electron beam produced by the laser was used. The spectral modelling showed the increased plasma density by a factor by 2-3 in the magnetic field. This was in agreement with the observed dynamics of the plasma plume confined by the axial magnetic field. The two-plasmon decay was studied in the 2.5-3 MG transversal magnetic fields. Strong 2-3.5 nm spectral widening and a 2-4 nm shifts of “red” and “blue” 3/2ω 0 spectral components were observed. The large shift and widening exceed the expected temperature and magnetic shift. PIC simulations are performing to clarify physics of this effect. The experimental research program for the grant is completed. New plasma effects in the MG fields were studied at the university-scale pulsed power generator. Two graduate and two undergraduate UNR and UR students carried experiments at the Zebra generator, supported and developed plasma diagnostics, provided data processing, and performed MHD simulations of LPI in the strong magnetic field. Research results are published in 9 referred papers and presented in 12 conferences. Program manager: Kramer U. Akli; Collaborators: R. C. Mancini and H. Sawada, University of Nevada, Reno (UNR); R. Betti and A. V. Maximov, University of Rochester (UR).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulating non-native cubic interactions on noisy quantum machines

As a milestone for general-purpose computing machines, we demonstrate that quantum processors can be programed to efficiently simulate dynamics that are not native to the hardware. Moreover, on noisy devices without error correction, we show that simulation results are significantly improved when the quantum program is compiled using modular gates instead of a restricted set of standard gates. We demonstrate the general methodology by solving a cubic interaction problem, which appears in nonlinear optics, gauge theories, as well as plasma and fluid dynamics. To encode the non-native Hamiltonian evolution, we decompose the Hilbert space into a direct sum of invariant subspaces in which the nonlinear problem is mapped to a finite-dimensional Hamiltonian simulation problem. Furthermore, in a three-states example, the resultant unitary evolution is realized by a product of approximately 20 standard gates, using which approximately ten simulation steps can be carried out on state-of-the-art quantum hardware before results are corrupted by decoherence. In comparison, the simulation depth is improved by more than an order of magnitude when the unitary evolution is realized as a single cubic gate, which is compiled directly using optimal control. Alternatively, parametric gates may also be compiled by interpolating control pulses. Modular gates thus obtained provide high-fidelity building blocks for quantum Hamiltonian simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

High-power tunable laser driven THz generation in corrugated plasma waveguides

The excitation of Terahertz (THz) radiation by the interaction of an ultrashort laser pulse with the modes of a miniature corrugated plasma waveguide is considered. Here, the axially corrugated waveguide supports the electromagnetic modes with appropriate polarization and subluminal phase velocities that can be phase matched to the ponderomotive potential associated with the laser pulse, making significant THz generation possible. This process is studied via full format Particle-in-Cell simulations that, for the first time, model the nonlinear dynamics of the plasma and the self-consistent evolution of the laser pulse in the case where the laser pulse energy is entirely depleted. It is found that the generated THz is characterized by lateral emission from the channel, with a spectrum that may be narrow or broad depending on the laser intensity. A range of realistic laser pulse and plasma parameters is considered with the goal of maximizing the conversion efficiency of optical energy to THz radiation. As an example, a fixed drive pulse (0.55 J) with a spot size of 15 μm and a duration of 15 fs produces a THz radiation of 37.8 mJ of in a 1.5 cm corrugated plasma waveguide with an on axis average density of 1.4 × 10 18 cm –3 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

The locking and unlocking thresholds for tearing modes in a cylindrical tokamak

The locking and unlocking thresholds for tearing modes are in general different. In this work, the physics origin for this difference is illustrated from theory analysis, and a numerical procedure is developed to find both locking and unlocking thresholds. In particular, a new scaling law for the unlocking threshold that is valid in both weak and strong rotation regimes has been derived from the lowest amplitude of the RMP (resonant magnetic perturbation) allowed for the locked-mode solution. Above the unlocking threshold, the criterion for the phase-flip instability is extended to identify the entire locked-mode states. Two different regimes of the RMP amplitude in terms of the accessibility of the locked-mode states have been found. In the first regime, the locked-mode state may or may not be accessible depending on the initial conditions of an evolving island. In the second regime, the locked-mode state can always be reached regardless of the initial conditions of the tearing mode. The lowest RMP amplitude for the second regime is determined to be the mode-locking threshold. The different characteristics of the two regimes above the unlocking threshold reveal the underlying physics for the gap between the locking and unlocking thresholds and provide an explanation for the closely related and widely observed hysteresis phenomena in island evolution during the sweeping process of the RMP amplitude up and down across that threshold gap.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

How do magnetized HED flows transition to collimated plasma jets? (Final Report)

We present results from a pulsed-power experiment designed to generate accretion-like flows in the laboratory that may be used to study the dynamics of astrophysical accretion-driven outflows. To emulate systems like Young Stellar Object (YSO) outflows, we create radially converging plasma flows in a short disk that transition into much longer bipolar outflows. Our experiment uses a ~1 MA pulsed-power driver to generate plasma from an array of thin aluminum wires mounted to a 3D-printed current path that controls the angular momentum and magnetic field topology. In this paper, we show that experimental data scale favorably to observed properties of YSO jets. This set of experiments tests the simplest version of our possible plasma dynamics, where outflows are generated without rotation or axial magnetic field present. The temperature, velocity and density of the disk and outflows are characterized using interferometry, gated optical and ultraviolet imaging, and Thomson scattering diagnostics. We discuss the Reynolds number, magnetic Reynolds number and Mach numbers calculated in three different stages of the experiment, and compare our results to 3D numerical simulations in the PERSEUS magnetohydrodynamics code. This report is a summary of the work done by the graduate students who were supported by the project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Plasma–liquid interactions: a review and roadmap

We present that plasma-liquid interactions represent a growing interdisciplinary area of research involving plasma science, fluid dynamics, heat and mass transfer, photolysis, multiphase chemistry and aerosol science. This review provides an assessment of the state-of-the-art of this multidisciplinary area and identifies the key research challenges. The developments in diagnostics, modeling and further extensions of cross section and reaction rate databases that are necessary to address these challenges are discussed. Finally, the review focuses on non-equilibrium plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Helium atmospheric pressure plasma jets touching dielectric and metal surfaces

Atmospheric pressure plasma jets (APPJs) are being investigated in the context plasma medicine and biotechnology applications, and surface functionalization. The composition of the surface being treated ranges from plastics, liquids, and biological tissue, to metals. The dielectric constant of these materials ranges from as low as 1.5 for plastics to near 80 for liquids, and essentially infinite for metals. The electrical properties of the surface are not independent variables as the permittivity of the material being treated has an effect on the dynamics of the incident APPJ. In this paper, results are discussed from a computational investigation of the interaction of an APPJ incident onto materials of varying permittivity, and their impact on the discharge dynamics of the plasma jet. The computer model used in this investigation solves Poisson's equation, transport equations for charged and neutral species, the electron energy equation, and the Navier-Stokes equations for the neutral gas flow. The APPJ is sustained in He/O{sub 2} = 99.8/0.2 flowing into humid air, and is directed onto dielectric surfaces in contact with ground with dielectric constants ranging from 2 to 80, and a grounded metal surface. Low values of relative permittivity encourage propagation of the electric field into the treated material and formation and propagation of a surface ionization wave. High values of relative permittivity promote the restrike of the ionization wave and the formation of a conduction channel between the plasma discharge and the treated surface. The distribution of space charge surrounding the APPJ is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Molecular-Dynamics Analysis of the Mechanical Behavior of Plasma-Facing Tungsten

We report a systematic computational analysis of the mechanical behavior of plasma-facing component (PFC) tungsten focusing on the impact of void and helium (He) bubble defects on the mechanical response beyond the elastic regime. Specifically, we explore the effects of porosity and He atomic fraction on the mechanical properties and structural response of PFC tungsten, at varying temperature and bubble size. We find that the Young modulus of defective tungsten undergoes substantial softening that follows an exponential scaling relation as a function of matrix porosity and He atomic content. Beyond the elastic regime, our high strain rate simulations reveal that the presence of nanoscale spherical defects (empty voids and He bubbles) reduces the yield strength of tungsten in a monotonically decreasing fashion, obeying an exponential scaling relation as a function of tungsten matrix porosity and He concentration. Our detailed analysis of the structural response of PFC tungsten near the yield point reveals that yielding is initiated by emission of dislocation loops from bubble/matrix interfaces, mainly 1/2<111> shear loops, followed by gliding and growth of these loops and reactions to form <100> dislocations. Furthermore, dislocation gliding on the <111>{211} twin systems nucleates 1/6<111> twin regions in the tungsten matrix. These dynamical processes reduce the stress in the matrix substantially. Subsequent dislocation interactions and depletion of the twin phases via nucleation and propagation of detwinning partials lead the tungsten matrix to a next deformation stage characterized by stress increase during applied straining. Furthermore, our structural analysis reveals that the depletion of twin boundaries (areal defects) is strongly impacted by the density of He bubbles at higher porosities. After the initial stress relief upon yielding, increase in the dislocation density in conjunction with decrease in the areal defect density facilitates the initiation of dislocation-driven deformation mechanisms in the PFC crystal.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Energy transfer in reconnection and turbulence

Reconnection and turbulence are two of the most commonly observed dynamical processes in plasmas, but their relationship is still not fully understood. Using 2.5D kinetic particle-in-cell simulations of both strong turbulence and reconnection, we compare the cross-scale transfer of energy in the two systems by analyzing the generalization of the von Kármán Howarth equations for Hall magnetohydrodynamics, a formulation that subsumes the third-order law for steady energy transfer rates. Even though the large scale features are quite different, the finding is that the decomposition of the energy transfer is structurally very similar in the two cases. In the reconnection case, the time evolution of the energy transfer also exhibits a correlation with the reconnection rate. Finally, these results provide explicit evidence that reconnection dynamics fundamentally involves turbulence-like energy transfer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Holographic entropy and real-time dynamics of quarkonium dissociation in non-Abelian plasma

The peak of the heavy quark pair entropy at the deconfinement transition, observed in lattice QCD, suggests that the transition is effectively driven by the increase of the entropy of bound states. The growth of the entropy with the interquark distance leads to the emergent entropic force that induces dissociation of quarkonium states. Since the quark-gluon plasma around the transition point is a strongly coupled system, we use the gauge-gravity duality to study the entropy of heavy quarkonium and the real-time dynamics of its dissociation. In particular, we employ the improved holographic QCD model as a dual description of large N c Yang-Mills theory. Studying the dynamics of the fundamental string between the quarks placed on the boundary, we find that the entropy peaks at the transition point. We also study the real-time dynamics of the system by considering the holographic string falling in the black hole horizon where it equilibrates. As a result, in the vicinity of the deconfinement transition, the dissociation time is found to be less than a fermi, suggesting that the entropic destruction is the dominant dissociation mechanism in this temperature region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2016

Second order gyrokinetic theory for particle-in-cell codes

The main idea of the gyrokinetic dynamical reduction consists in a systematical removal of the fast scale motion (the gyromotion) from the dynamics of the plasma, resulting in a considerable simplification and a significant gain of computational time. The gyrokinetic Maxwell–Vlasov equations are nowadays implemented in for modeling (both laboratory and astrophysical) strongly magnetized plasmas. Different versions of the reduced set of equations exist, depending on the construction of the gyrokinetic reduction procedure and the approximations performed in the derivation. The purpose of this article is to explicitly show the connection between the general second order gyrokinetic Maxwell–Vlasov system issued from the modern gyrokinetic theory and the model currently implemented in the global electromagnetic Particle-in-Cell code ORB5. Necessary information about the modern gyrokinetic formalism is given together with the consistent derivation of the gyrokinetic Maxwell–Vlasov equations from first principles. The variational formulation of the dynamics is used to obtain the corresponding energy conservation law, which in turn is used for the verification of energy conservation diagnostics currently implemented in ORB5. This work fits within the context of the code verification project VeriGyro currently run at IPP Max-Planck Institut in collaboration with others European institutions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Data-driven linear time advance operators for the acceleration of plasma physics simulation

In this study, we demonstrate the application of data-driven linear operator construction for time advance with a goal of accelerating plasma physics simulation. We apply dynamic mode decomposition (DMD) to data produced by the nonlinear SOLPS-ITER (Scrape-off Layer Plasma Simulator - International Thermonuclear Experimental Reactor) plasma boundary code suite in order to estimate a series of linear operators and monitor their predictive accuracy via online error analysis. We find that this approach defines when these dynamics can be represented by a sequence of approximate linear operators and is essential for providing consistent projections when compared to an unconstrained application. For linear diffusion and advection–diffusion fluid test problems, we construct and apply operators within explicit and implicit time advance schemes, demonstrating that stability can be robustly guaranteed in each case. We further investigate the use of the linear time advance operators within several integration methods including forward Euler, backward Euler, and the matrix exponential. The application of this method to simulation data from SOLPS-ITER, with varying levels of Markov chain Monte Carlo numerical noise, shows that constrained DMD operators yield a capability to identify, extract, and integrate a (slow) subset of the present timescales. Example applications show that for projected speedup factors of [Formula: see text], and [Formula: see text], a mean relative error of 3%, 5%, and 8% and maximum relative error less than 20% are achievable, which appears acceptable for typical SOLPS-ITER steady-state simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hamiltonian formulations for perturbed dissipationless plasma equations

The Hamiltonian formulations for the perturbed Vlasov–Maxwell equations and the perturbed ideal magnetohydrodynamics (MHD) equations are expressed in terms of the perturbation derivative ∂ F / ∂ ϵ ≡ [ F , S ] of an arbitrary functional F [ ψ ] of the Vlasov–Maxwell fields ψ = ( f , E , B ) or the ideal MHD fields ψ = ( ρ , u , s , B ) , which are assumed to depend continuously on the (dimensionless) perturbation parameter ϵ . In this study, [ , ] denotes the functional Poisson bracket for each set of plasma equations and the perturbation action functional S is said to generate dynamically accessible perturbations of the plasma fields. The new Hamiltonian perturbation formulation introduces a framework for functional perturbation methods in plasma physics and highlights the crucial roles played by polarization and magnetization in Vlasov–Maxwell and ideal MHD perturbation theories. One application considered in this paper is a formulation of plasma stability that guarantees dynamical accessibility and leads to a natural generalization to higher-order perturbation theory.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Generation of disc-like plasma from laser-matter interaction in the presence of a strong external magnetic field

Dynamics of laser produced plasma in a strong magnetic field was studied here using a 1 MA pulsed power generator coupled to an intense, high-energy laser. A 2–2.5 MG magnetic field was generated on the surface of a rod load 0.8–1.2 mm in diameter. A sub-nanosecond laser pulse with intensity of 3 × 10 15 W cm -2 was focused on the rod load surface. Side-on laser diagnostics showed the generation of two collimated jets 1–3 mm long on the front and rear sides of the load. End-on laser diagnostics reveal that the laser produced plasma in the MG magnetic field takes the form of a thin disc as the plasma propagates along the magnetic field lines. The disc-like plasma expands radially across the magnetic field with a velocity of 250 km s -1 . An electron temperature of 400 eV was measured in the laser-produced plasma on the rod load.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Investigation of shock-shock interaction and Mach reflection in laterally colliding laser-blow-off plasmas

Interactions of two Li plasma plumes and shock waves are investigated at various pressures (∼10{sup −5} to 3 mbar) in the argon gas ambient. Fast imaging and optical emission spectroscopy are used to study the plume dynamics and characteristic emission of plasmas. The plasma plumes are created in laser-blow-off geometry. The expansion of plasma plumes in the ambient gas leads to the formation of an interaction zone. The formation of interaction zone is dependent on the ambient pressure and below a certain pressure, no significant change is observed in the shape and size of the interaction plasma. In the higher pressure, formation of interaction zone and its shape are dependent on ambient pressure. Dynamics of seed plasmas and interaction zone are also affected by the shock-shock interactions. The shock-shock interaction depends on the angle of incidence (α) between two shock waves at the initial time of interaction but as the plumes expand, the shock-shock interaction does not follow α dependence.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Molecular dynamics simulations of Si etching in Cl- and Br-based plasmas: Cl{sup +} and Br{sup +} ion incidence in the presence of Cl and Br neutrals

Classical molecular dynamics (MD) simulations have been performed for Cl{sup +} and Br{sup +} ions incident on Si(100) surfaces with Cl and Br neutrals, respectively, to gain a better understanding of the ion-enhanced surface reaction kinetics during Si etching in Cl- and Br-based plasmas. The ions were incident normally on surfaces with translational energies in the range E{sub i} = 20–500 eV, and low-energy neutrals of E{sub n} = 0.01 eV were also incident normally thereon with the neutral-to-ion flux ratio in the range Γ{sub n}{sup 0}/Γ{sub i}{sup 0} = 0–100, where an improved Stillinger--Weber potential form was employed for the interatomic potential concerned. The etch yields and thresholds presently simulated were in agreement with the experimental results previously reported for Si etching in Cl{sub 2} and Br{sub 2} plasmas as well as in Cl{sup +}, Cl{sub 2}{sup +}, and Br{sup +} beams, and the product stoichiometry simulated was consistent with that observed during Ar{sup +} beam incidence on Si in Cl{sub 2}. Moreover, the surface coverage of halogen atoms, halogenated layer thickness, surface stoichiometry, and depth profile of surface products simulated for Γ{sub n}{sup 0}/Γ{sub i}{sup 0} = 100 were in excellent agreement with the observations depending on E{sub i} reported for Si etching in Cl{sub 2} plasmas. The MD also indicated that the yield, coverage, and surface layer thickness are smaller in Si/Br than in Si/Cl system, while the percentage of higher halogenated species in product and surface stoichiometries is larger in Si/Br. The MD further indicated that in both systems, the translational energy distributions of products and halogen adsorbates desorbed from surfaces are approximated by two Maxwellians of temperature T{sub 1} ≈ 2500 K and T{sub 2} ≈ 7000–40 000 K. These energy distributions are discussed in terms of the desorption or evaporation from hot spots formed through chemically enhanced physical sputtering and physically enhanced chemical sputtering, which have so far been speculated to both occur in the ion-enhanced surface reaction kinetics of plasma etching.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

A splitting Chebyshev collocation method for Schrödinger–Poisson system

We develop a splitting Chebyshev collocation (SCC) method for the time-dependent Schrödinger–Poisson (SP) system arising from theoretical analysis of quantum plasmas. By means of splitting technique in time, the time-dependant SP system is first reduced to uncoupled Schrödinger and Poisson equations at every time step. The space variables in Schrödinger and Poisson equations are next represented by high-order Chebyshev polynomials, and the resulting system are discretized by the spectral collocation method. Finally, matrix diagonalization technique is applied to solve the fully discretized system in one dimension, two dimensions and three dimensions, respectively. The newly proposed method not only achieves spectral accuracy in space but also reduces the computer-memory requirements and the computational time in comparison with conventional solver. Numerical results confirm the spectral accuracy and efficiency of this method, and indicate that the SCC method could be an efficient alternative method for simulating the dynamics of quantum plasmas.

97 MATHEMATICS AND COMPUTING↗