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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,261 records · Page 70

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Towards modelling AR Sco: generalized particle dynamics and strong radiation-reaction regimes

Numerical simulations of relativistic plasmas have become more feasible, popular, and crucial for various astrophysical sources with the availability of computational resources. The necessity for high-accuracy particle dynamics is especially highlighted in pulsar modelling due to the extreme associated electromagnetic fields and particle Lorentz factors. Including the radiation-reaction force in the particle dynamics adds even more complexity to the problem, but is crucial for such extreme astrophysical sources. We have also realized the need for such modelling concerning magnetic mirroring and particle injection models proposed for AR Sco, the first white dwarf pulsar. This paper demonstrates the benefits of using higher-order explicit numerical integrators with adaptive time-step methods to solve the full particle dynamics with radiation-reaction forces included. We show that for standard test scenarios, namely various combinations of uniform E- and B-fields and a static dipole B-field, the schemes we use are equivalent to and in extreme field cases outperform standard symplectic integrators in accuracy. We show that the higher-order schemes have massive computational time improvements due to the adaptive time-steps we implement, especially in non-uniform field scenarios and included radiation reaction where the particle gyro-radius rapidly changes. When balancing accuracy and computational time, we identified the adaptive Dormand–Prince eighth-order scheme to be ideal for our use cases. The schemes we use maintain accuracy and stability in describing the particle dynamics and we indicate how a charged particle enters radiation-reaction equilibrium and conforms to the analytical Aristotelian Electrodynamics expectations.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Adiabatic behavior of an elliptical vortex in a time-dependent external strain flow

The adiabatic behavior of two-dimensional (2D) fluid vortices subject to slowly varying, time-dependent external strain flows is studied theoretically, and experimentally using magnetized non-neutral electron plasmas. Here the E × B drift dynamics of the electrons perpendicular to the magnetic field are analogous to the motion of fluid vorticity under the 2D Euler equations describing ideal fluids. A low-dimensional elliptical vortex patch model is used along with a WKB approximation to derive a formula for the breaking of an adiabatic invariant due to the changing external flow. The invariant is interpreted as the amplitude of a perturbative oscillation about a stable fixed point corresponding to a steady elliptical vortex. Smooth, hyperbolic tangent and piecewise linear ramp functions are considered for the external strain time dependence. Standard exponential breaking is observed in the former case, whereas the latter exhibits a power-law breaking curve with periodic modulation. It is found that a driving term in the equations of motion contributes the majority of the breaking whereas the frequency variation plays a weaker, but significant role. The experimental data agree closely with the theoretical model. The most significant deviation is due to inviscid damping behavior associated with the smooth edges of the experimental vortices, which tends to reduce the amplitude of the oscillation. Finally, the results are compared and contrasted with other related experimental, numerical, and theoretical work.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1995

Kinetic structure of intermediate shocks: Implications for the magnetopause

A general study of the structure and stability of intermediate shocks (IS) in an isotropic plasma is presented using a hybrid as well as a resistive Hall MHD code. Special emphasis is put on the question of whether the rotational layers observed at the magnetopause can be intermediate shocks. The shocks are formed dynamically by the interaction between a flowing plasma and a stationary piston. Coplanar ISs (both strong and weak) are found to be stable in a collisionless plasma. The existence of slow shocks in a high beta plasma is also established for the first time. Noncoplanar ISs are found to be time-dependent, evolving toward a rotational discontinuity (RD) after some characteristic time tau which can be quite long (1000 Omega(exp -1), where Omega is the ion gyrofrequency). The value tau is larger the closer the rotation angle is to 180 deg. Rotations larger than 180 deg are found to be unstable, decaying into a state of minimum shear (i.e., rotation angle less than 180 deg). There are various length scales associated with an IS in the kinetic regime. The shortest scale is found to be the length scale over which rotation of the transverse component of the magnetic field takes place. This scale can have a half width as small as one ion inertial length (c/omega(sub p)) for electron sense rotations and 3c/omega(sub p) for ion sense rotations, for an upstream ion beta of unity. Both of these scales are consistent with the observed thickness at the magnetopause and identical to the corresponding RD scales. A detailed study of the mode conversion of the Alfven ion cyclotron waves (A/IC) waves across both slow and intermediate shocks and the resulting downstream wave spectrum are presented. The possibility that the large number of relfected ions observed at the magnetopause may be due to the presence of strong ISs is considered. The identification of strong ISs and their distinction from RDs should be possible in observations due to significant differences that exist between jump conditions and overall structure of the two discontinuities. The jumps in the plasma parameters across a weak IS are typically small. This together with the fact that the weak ISs and RDs have very similar thickness and other overall properties makes the distinction between weak ISs and RDs in the observations largely inconsequential. However, at large noncoplanarity angles the weak IS approaches the RD limit in a relatively short time (approximately less than 100 Omega (exp -1)). Thus, magnetopause rotations with large noncoplanarity angles are most likely either RDs or strong ISs. Finally, direct comparisons between fluid (resistive Hall MHD) and kinetic simulations show that fluid theory is not applicable to study of ISs in a collisionless plasma.

Karimabadi, H.↗

NASA NTRS1992

Prospects for optical and radio-wave remote imaging of high latitude exospheric processes

The paper examines the use of optical, radio-wave, and extreme UV imaging to study the distribution of magnetospheric plasma continuously and thereby study magnetospheric global dynamics. Optical imaging of the terrestrial exosphere at FUV and EUV wavelengths is discussed in the light of evidence from recent observational campaigns. Radio imaging of the earth's magnetosphere is shown to provide good instantaneous images of the plasma from beyond 4 earth radii. Optical imaging is examined with respect to the study of other planetary exospheres, and the technique is shown to depend on the ability to do normal incidence optics in the EUV. The Cassini radio-uplink experiment is examined as an example of radio-wave tomography which can provide geometry-based tomographic montages to image the heliosphere and the planetary magnetospheres.

Chiu, Y. T.↗

NASA NTRS1994

High Altitude Plasma Instrument (HAPI) data analysis

The objectives of the Dynamics Explorer mission are to investigate the coupling of energy, mass, and momentum among the earth's magnetosphere, ionosphere, and upper atmosphere. At launch, on August 3, 1981, DE-1 was placed into an elliptical polar orbit having an apogee of 23,130 km to allow global auroral imaging and crossings of auroral field lines at altitudes of several thousand kilometers. At the same time DE-2 was placed into a polar orbit, coplanar with that of DE-1 but with a perigee altitude low enough (309 km) for neutral measurements and an apogee altitude of 1012 km. The DE-1 High Altitude Plasma Instrument (HAPI) provided data on low and medium energy electrons and ions from August 13, 1981 until December 1, 1981, when a high-voltage failure occured. Analysis of HAPI data for the time period of this contract has produced new results on the source mechanisms for electron conical distributions, particle acceleration phenomena in auroral acceleration regions, Birkeland currents throughout the nightside auroral regions, the source region for auroral kilometric radiation (AKR), and plasma injection phenomena in the polar cusp.

Burch, J. L.↗

DOE OSTI2021

The Consequences of Rayleigh-Taylor Instabilities on Implosion Dynamics in Gas-puff Z-pinch Experiments (Final Report)

The Laboratory of Plasma Studies (LPS), Cornell University, proposed to carry out a comprehensive series of gas-puff z-pinch experiments on a 1 MA pulsed power machine for the purpose of understanding the dynamics of imploding magnetized high energy density plasmas in the presence of Rayleigh-Taylor Instabilities (RTI). Cornell’s pulsed power facilities, including the associated diagnostic instruments, are largely supported under the auspices of the Center for Pulsed-power-driven High-energy-density Plasmas (the Center), an NNSA-sponsored Center of Excellence. The proposed experiments were carried out predominantly by LPS graduate student, Sophia Rocco (the “primary grad student”), in collaboration with other graduate students and the scientific and technical staff of the Center using Cornell’s 1 MA COBRA pulsed power facility and all of its associated diagnostic instruments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Plasma Surface Interactions: Bridging from the Surface to the Micron Frontier through Leadership Class Computing (Final Technical Progress Report)

The performance demands on plasma-facing components (PFCs) in future fusion power plants are beyond the capability of current materials. In defining the plasma-surface interaction (PSI), it is evident that three coupled spatial regions influence PFC materials evolution and performance. These regions consist of (1) the edge and scrape-off layer region of the plasma, (2) the near-surface material response to extreme thermal and particle fluxes under the influence of, and feedback to, the plasma sheath, and (3) the structural material’s response to an intense, 14 MeV peaked neutron spectrum, which produces very high concentrations of transmuted elements through (n,p) and (n,α) reactions and structural material property degradation. The objective of this project is to develop and deploy validated, high-performance simulation tools capable of predicting the performance of tungsten-based plasma-facing components (PFCs) in a burning fusion plasma environment, which includes modeling surface morphology evolution in either erosion or re-deposition regimes, and the recycling of hydrogen species. This requires the development of a leadership-scale computational code to predict PFC behavior that is well integrated to a suite of multiscale modeling techniques to bridge the scales needed to address complex physical and computational issues at the plasma surface interface and the transition region below the surface where neutron damage processes in the bulk dominate material behavior. Successful completion of this project will provide PFC simulation tools to evaluate steady-state performance of tungsten-based PFC and divertor components in burning plasma environments. This will enable the identification of critical experiments to confirm whether practical PFC solutions exist for magnetic fusion energy beyond ITER. More specifically, the research activities within this project focus on two broadly defined research thrusts: 1) Developing a new simulation code, Xolotl-PSI to predict PFC operating lifetime and performance, that is specifically designed to take advantage of leadership class computing facilities; and 2) Integrating and applying discrete particle-based, as well as continuum-based, multiscale modeling techniques to provide scientific discovery of the mechanisms controlling PFC and bulk materials evolution under fusion plasma and 14-MeV neutron exposure; in this thrust, the challenge of scale-bridging between atomistic/microstructural modeling and the continuum-scale of PFC device performance is achieved through multiscale modeling techniques. Throughout the grant period, research at UMass Amherst focused on the atomic-scale analysis of the dynamics and kinetics (transport and reactions) of small mobile helium clusters near surfaces and grain boundaries in plasma-exposed tungsten, as well as the proper parameterization and incorporation of these processes into our SciDAC team’s continuum-level cluster dynamics code (Xolotl). Our work also included characterization of large-scale molecular-dynamics (MD) simulation results of helium implantation into tungsten based on the findings of our analysis and comparisons of Xolotl simulation predictions with the large-scale MD simulation results for Xolotl validation purposes; emphasis was placed on the impact of helium cluster dynamics/kinetics on tungsten surface morphology, near-surface tungsten structure, and helium retention. Furthermore, our work included initiation of a research program on evaluating the impact of structural defects in tungsten due to plasma exposure (such as helium nanobubbles) on tungsten’s thermophysical properties, focusing on lattice thermal conductivity. Our studies contributed to a fundamental understanding of the effects of plasma-surface interactions on the dynamical response of plasma-facing materials in nuclear fusion devices and the development of predictive hierarchical multiscale computational tools for the quantitative description of the material’s dynamical response to plasma exposure.

36 MATERIALS SCIENCE↗

NASA NTRS1995

Semikinetic modeling and observations of high-latitude plasma outflow

In one fo the most exciting areas of progress, we have now developed a dynamic semikinetic model for examining the synergistic effects of waves and magnetospheric hot plasma populations on the outflowing ionospheric plasma. We have done this by imposing hot biMaxwellian ion and electron distributions at the top of our auroral simulation flux tube (4 R(sub e)), as well as a spectrum of waves with altitude which perpendicularly heats the ionospheric ions. We have also addressed the quasi-statistical properties of out flowing O(+) through bulk parameter analysis of DE-1/RIMS observations when DE-1 was in the midaltitude polar cap magnetosphere. A very exciting paper which was both submitted and appeared during this period in JGR concerned the centrifugal acceleration effect on the polar wind. During this year, we also developed the extension of our GSK code to include the lower ionosphere and associated processes. Finally, perhaps our most exciting project currently is a study of F-region upflows using a modification of the FLIP ionospheric fluid dynamics model of Richards and Torr to allow for the effects of soft electron precipitation ionization and convection-driven ion heating, and performed comparisons with satellite and radar data. A list of reference papers that have been published is listed.

Horwitz, James L.↗

DOE OSTI2018

Quantifying Platelet Margination in Diabetic Blood Flow

According to this article, patients with type 2 diabetes mellitus (T2DM) develop thrombotic abnormalities strongly associated with cardiovascular diseases. In addition to the changes of numerous coagulation factors such as elevated levels of thrombin and fibrinogen, the abnormal rheological effects of red blood cells (RBCs) and platelets flowing in blood are crucial in platelet adhesion and thrombus formation in T2DM. An important process contributing to the latter is the platelet margination. We employ the dissipative particle dynamics method to seamlessly model cells, plasma, and vessel walls. We perform a systematic study on RBC and platelet transport in cylindrical vessels by considering different cell shapes, sizes, and RBC deformabilities in healthy and T2DM blood, as well as variable flowrates and hematocrit. In particular, we use cellular-level RBC and platelet models with parameters derived from patient-specific data and present a sensitivity study. We find T2DM RBCs, which are less deformable compared to normal RBCs, lower the transport of platelets toward the vessel walls, whereas platelets with higher mean volume (often observed in T2DM) lead to enhanced margination. Furthermore, increasing the flowrate or hematocrit enhances platelet margination. We also investigated the effect of platelet shape and observed a nonmonotonic variation with the highest near-wall concentration corresponding to platelets with a moderate aspect ratio of 0.38. We examine the role of white blood cells (WBCs), whose count is increased notably in T2DM patients. We find that WBC rolling or WBC adhesion tends to decrease platelet margination due to hydrodynamic effects. To the best of our knowledge, such simulations of blood including all blood cells have not been performed before, and our quantitative findings can help separate the effects of hydrodynamic interactions from adhesive interactions and potentially shed light on the associated pathological processes in T2DM such as increased inflammatory response, platelet activation and adhesion, and ultimately thrombus formation.

Biophysics↗

DOE OSTI2021

Non-conformal attractor in boost-invariant plasmas

We study the dissipative evolution of (0+1)-dimensionally expanding media with Bjorken symmetry using the Boltzmann equation for massive particles in relaxation-time approximation. Breaking conformal symmetry by a mass induces a non-zero bulk viscous pressure in the medium. It is shown that even a small mass (in units of the local temperature) drastically modifies the well-known attractor for the shear Reynolds number previously observed in massless systems. For generic nonzero particle mass, neither the shear nor the bulk viscous pressure relax quickly to a non-equilibrium attractor; they approach the hydrodynamic limit only late, at small values of the inverse Reynolds numbers. Only the longitudinal pressure, which is a combination of thermal, shear and bulk viscous pressures, continues to show early approach to a far-off-equilibrium attractor, driven by the rapid longitudinal expansion at early times. Second-order dissipative hydrodynamics based on a gradient expansion around locally isotropic thermal equilibrium fails to reproduce this attractor.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2019

Probing ion species separation and ion thermal decoupling in shock-driven implosions using multiple nuclear reaction histories

Simultaneously measured DD, DT, and D 3 He reaction histories are used to probe the impacts of multi-ion physics during the shock phase of inertial confinement fusion implosions. In these relatively hydrodynamiclike (burn-averaged Knudsen number &lang;N K &rang; ~0.3) shock-driven implosions, average-ion hydrodynamic DUED simulations are able to reasonably match burnwidths, nuclear yields, and ion temperatures. However, kinetic-ion FPION simulations are able to better simulate the timing differences and time-resolved reaction rate ratios between DD, DT, and D 3 He reactions. FPION simulations suggest that the D 3 He/DT reaction rate ratio is most directly impacted by ion species separation between the 3 He and T ions, whereas the D 3 He/DD reaction rate ratio is affected by both ion species separation and ion temperature decoupling effects.

Hydrodynamic codes↗

DOE OSTI2021

On the merit of hot ion mode for tearing mode stabilization

The stabilization of tearing modes with rf driven current benefits from the cooperative feedback loop between rf power deposition and electron temperature within the island. This effect, termed rf current condensation, can greatly enhance and localize the current driven within magnetic islands. It has previously been shown that the condensation effect opens the possibility of passive stabilization with broad rf profiles, as would be typical of LHCD for steady state operation. Here, we show that this self-healing effect can be dramatically amplified by operation in a hot ion mode, due to the additional electron heat source provided by the hotter ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Linear stability analysis via simulated annealing and accelerated relaxation

Simulated annealing (SA) is a kind of relaxation method for finding equilibria of Hamiltonian systems. A set of evolution equations is solved with SA, which is derived from the original Hamiltonian system so that the energy of the system changes monotonically while preserving Casimir invariants inherent to noncanonical Hamiltonian systems. The energy extremum reached by SA is an equilibrium. Since SA searches for an energy extremum, it can also be used for stability analysis when initiated from a state where a perturbation is added to an equilibrium. The procedure of the stability analysis is explained, and some examples are shown. Because the time evolution is computationally time consuming, efficient relaxation is necessary for SA to be practically useful. An acceleration method is developed by introducing time dependence in the symmetric kernel used in the double bracket, which is part of the SA formulation described here. An explicit formulation for low-beta reduced magnetohydrodynamics (MHD) in cylindrical geometry is presented. In conclusion, since SA for low-beta reduced MHD has two advection fields that relax, it is important to balance the orders of magnitude of these advection fields.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Ground state rotational and kinetic temperatures of C 2 molecules in a laser-produced plasma

This study presents an experimental comparison between rotational and kinetic temperatures of molecular species in a laser-produced plasma, utilizing tunable laser absorption spectroscopy. We analyzed high-resolution and time-resolved absorption spectra of 31 rotational transitions within the Phillips A 1 II u –X 1 Σ g 3-0 band of the C 2 molecule, measured from a laser-produced graphite plasma in ~2 Torr argon, within the 771–781 nm spectral range. Spectral fitting provided quantitative state column densities of C 2 rotational energy levels in the ground electronic and vibrational states. Boltzmann analysis revealed thermalized rotational distributions showing rotational temperatures decreasing from 3000 to 600 K within 5–500 μs following the onset of the laser plasma. The measured kinetic and rotational temperatures of C 2 molecules are comparable, suggesting that collisional processes in the laser-produced plasma are sufficient to maintain thermal equilibrium between the C 2 molecules and the Ar background gas.

Absorption spectroscopy↗

DOE OSTI2023

Precision Tests of the Nonlinear Mode Coupling of Anisotropic Flow via High-Energy Collisions of Isobars

Valuable information on dynamics of expanding fluids can be inferred from the response of such systems to perturbations in their initial geometry. We apply this technique in high-energy 96 Ru+ 96 Ru and 96 Zr+ 96 Zr collisions to scrutinize the expansion dynamics of the quark-gluon plasma, where the initial geometry perturbations are sourced by the differences in deformations and radial profiles between 96 Ru and 96 Zr, and the collective response is captured by the change in anisotropic flow V n between the two collision systems. Using a transport model, we analyze how the nonlinear coupling between lower-order flow harmonics V 2 and V 3 to the higher-order flow harmonics V 4 and V 5 , expected to scale as V 4 NL = χ 4 V 2 2 and V 5NL = χ 5 V 2 V 3 , gets modified as one moves from 96 Ru+ 96 Ru to 96 Zr+ 96 Zr systems. We find that these scaling relations are valid to high precision: variations of order 20% in V 4NL and V 5NL due to differences in quadrupole deformation, octupole deformation, and nuclear skin modify χ 4 and χ 5 by about 1–2%. Percent-level deviations are however larger than the expected experimental uncertainties and could be measured. Therefore, collisions of isobars with different nuclear structures are a unique tool to isolate subtle nonlinear effects in the expansion of the quark-gluon plasma that would be otherwise impossible to access in a single collision system.

Physics↗

DOE OSTI2025

Effects of size of an oxide layer on the surface properties of lithium irradiated by deuterium plasma

We investigate how growth of an oxide film will influence the deuterium recycling properties of Li on plasma facing surfaces. Lithium films on the walls or plasma-facing material surfaces of a fusion vacuum vessel improves plasma performance in part by removing residual impurity atoms from the plasma. Oxygen atoms, from residual water vapor or eroded oxide surfaces, are a prominent ubiquitous impurity that can create a layer of oxide on lithium films. The results of the work herein also indicate a possibility for determining the thickness of the oxide layer on such Li films by measuring the reflection probability of incident deuterium particles from the surface.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Impact of tracking resolution on 𝜙-meson spin alignment measurements

Measurements of global spin alignment of vector mesons in relativistic heavy-ion collisions can provide unique insights into spin-orbit interactions and vector meson dynamics in the quark-gluon plasma produced in those collisions. The global spin alignment is measured by the 00th coefficient of the spin density matrix, 𝜌 00 , via the polar angle (𝜃*) of the decay-kaon momentum in the parent rest frame with respect to the direction of the orbital angular momentum of the collision. Such measurements are affected by the angular and momentum resolutions of the reconstructed tracks in the experiment. Such effects are nontrivial because of kinematic complications caused by the boost to the parent rest frame, and could be important given that the global spin alignment signal is weak. In this paper, we investigate the effects of experimental tracking resolutions on measurements of the 𝜙 (1020) meson 𝜌 00 . We study these effects for two methods of 𝜌 00 measurements, the conventional method analyzing the 𝜙⁢−meson yield versus cos 2 ⁡𝜃* and the invariant mass (𝑚 inv ) method utilizing ⟨cos 2 ⁡𝜃*⟩ versus 𝑚 inv . Using typical resolution values from experiments, we find that the effect of track resolution on 𝜌 00 is small, well within typical measurement uncertainties.

Relativistic heavy-ion collisions↗