Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,783 records · Page 99

Cleaning graphene: A first quantum/classical molecular dynamics approach

Graphene outstanding properties created a huge interest in the condensed matter community and unprecedented fundings at the international scale in the hope of application developments. Recently, there have been several reports of incomplete removal of the polymer resists used to transfer as-grown graphene from one substrate to another, resulting in altered graphene transport properties. Finding a large-scale solution to clean graphene from adsorbed residues is highly desirable and one promising possibility would be to use hydrogen plasmas. In this spirit, we couple here quantum and classical molecular dynamics simulations to explore the kinetic energy ranges required by atomic hydrogen to selectively etch a simple residue—a CH{sub 3} group—without irreversibly damaging the graphene. For incident energies in the 2–15 eV range, the CH{sub 3} radical can be etched by forming a volatile CH{sub 4} compound which leaves the surface, either in the CH{sub 4} form or breaking into CH{sub 3} + H fragments, without further defect formation. At this energy, adsorption of H atoms on graphene is possible and further annealing will be required to recover pristine graphene.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Mixed stochastic-deterministic density functional theoretic decomposition of Kubo–Greenwood conductivities in the projector augmented wave formalism

Pairing the accuracy of the Kohn–Sham density-functional framework with the efficiency of a stochastic algorithmic approach, mixed stochastic-deterministic density functional theory (mDFT) achieves a favorable computational scaling with system sizes and electronic temperatures. We employ the recently developed mDFT formalism to investigate the dynamic charge-transport properties of systems in the warm dense matter regime. The optical conductivity spectra are computed for single- and multi- component mixtures of carbon, hydrogen, and beryllium using two complementary approaches: Kubo–Greenwood in the mDFT picture and real-time time-dependent mDFT. We further devise a decomposition of the Onsager coefficients leading up to the Kubo–Greenwood spectra to exhibit contributions from the deterministic, stochastic, and mixed electronic state transitions at different incident photon energies.

36 MATERIALS SCIENCE↗

Multi-dimensional dynamics of stimulated Brillouin scattering in a laser speckle: Ion acoustic wave bowing, breakup, and laser-seeded two-ion-wave decay

Two- and three-dimensional particle-in-cell simulations of stimulated Brillouin scattering (SBS) in laser speckle geometry have been analyzed to evaluate the relative importance of competing nonlinear processes in the evolution and saturation of SBS. It is found that ion-trapping-induced wavefront bowing and breakup of ion acoustic waves (IAW) and the associated side-loss of trapped ions dominate electron-trapping-induced IAW wavefront bowing and breakup, as well as the two-ion-wave decay instability over a range of ZTe/Ti conditions and incident laser intensities. In the simulations, the latter instability does not govern the nonlinear saturation of SBS; however, evidence of two-ion-wave decay is seen, appearing as a modulation of the ion acoustic wavefronts. This modulation is periodic in the laser polarization plane, anti-symmetric across the speckle axis, and of a wavenumber matching that of the incident laser pulse. A simple analytic model is provided for how spatial “imprinting” from a high frequency inhomogeneity (in this case, the density modulation from the laser) in an unstable system with continuum eigenmodes can selectively amplify modes with wavenumbers that match that of the inhomogeneity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Instability of an electron-plasma shear layer in an externally imposed strain flow

The E x B shear instability of a two-dimensional (2D) filament (i.e., a thin, rectangular strip perpendicular to the magnetic field) of magnetized pure electron plasma is investigated experimentally in the presence of an externally imposed strain flow. Data are acquired using a specialized Penning–Malmberg trap in which strain flows can be applied in 2D by biasing segmented electrodes surrounding the plasma. The E x B drift dynamics are well-described by the Drift-Poisson equations, which are isomorphic to the 2D Euler equations describing ideal fluids. Thus, the experimental results correspond to the Rayleigh instability of a shear layer in a 2D ideal fluid, where the electron density is analogous to the fluid vorticity. Shear layers are prepared by stretching initially axisymmetric electron vortices using a strong, applied strain flow. The data at early times are in quantitative agreement with a linear model which extends Rayleigh’s work to account for the influence of an external strain flow. In the presence of weak strain, the system approximately maintains a phase relationship that corresponds to an instantaneous Rayleigh eigenmode. The instability develops into the nonlinear regime later in time and at smaller spatial scales as the strain rate is increased. A secondary vortex pairing instability is observed, but it is suppressed when the strain-to-vorticity ratio exceeds roughly 0.025. In this way, vorticity transport perpendicular to the filament is diminished due to the applied strain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fully plasma-based electron injector for a linear collider or XFEL

We demonstrate through high-fidelity particle-in-cell (PIC) simulations a simple approach for efficiently generating 20 + GeV electron beams with the necessary charge, energy spread, and emittance for use as an injector in a future linear collider or a next generation XFEL. A high quality injected bunch is generated by self-focusing an unmatched electron driver in a nonlinear plasma wakefield. Over pump depletion distances, the drive beam dynamics and self-loading effects lead to high energy, low-energy spread output beams. For plasma densities of 10 18 c⁢m −3 , PIC simulation results indicate that self-injected beams with 0.52 n⁢C charge can be accelerated to 20 GeV with projected core energy spreads of ≲ 1%, normalized slice emittances of 110n⁢m, peak normalized brightness of ≳ 10 19 A/m 2 /rad 2 , and transfer efficiencies of ≳ 44%.

Particle acceleration in plasmas↗

Color-based tracking of plasma dust particles

Color-based tracking to observe agglomeration of deposited particles inside a compact planar magnetron during plasma discharge was done by creating high dynamic range (HDR) images of photos captured by a Pentax K10D digital camera. Carbon erosion and redeposition was also monitored using the technique. The HDR images were subjected to a chromaticity-based constraint discoloration inside the plasma chamber indicating film formation or carbon redeposition. Results show that dust deposition occurs first near the evacuation pumps due to the pressure gradient and then accumulates at the positively charged walls of the chamber. This method can be applied to monitor dust formation during dusty plasma experiments without major modification of plasma devices, useful especially for large fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Criterion for the sign of wave energy

Criterion for the sign of wave energy are developed by using the symmetry properties of the plasma equilibrium and the fact that Vlasov dynamics is an incompressible flow in phase space, rather than the usual and more difficult procedure of calculating the value of the wave energy directly. Applications are made to the case of waves excited on a nonneutral plasma in a Malmberg-Penning trap and to waves excited on an infinitely long nonneutral beam.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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↗

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↗

grim: A Flexible, Conservative Scheme for Relativistic Fluid Theories

Hot, diffuse, relativistic plasmas such as sub-Eddington black-hole accretion flows are expected to be collisionless, yet are commonly modeled as a fluid using ideal general relativistic magnetohydrodynamics (GRMHD). Dissipative effects such as heat conduction and viscosity can be important in a collisionless plasma and will potentially alter the dynamics and radiative properties of the flow from that in ideal fluid models; we refer to models that include these processes as Extended GRMHD. Here we describe a new conservative code, grim, that enables all of the above and additional physics to be efficiently incorporated. grim combines time evolution and primitive variable inversion needed for conservative schemes into a single step using an algorithm that only requires the residuals of the governing equations as inputs. This algorithm enables the code to be physics agnostic as well as flexibility regarding time-stepping schemes. grim runs on CPUs, as well as on GPUs, using the same code. We formulate a performance model and use it to show that our implementation runs optimally on both architectures. grim correctly captures classical GRMHD test problems as well as a new suite of linear and nonlinear test problems with anisotropic conduction and viscosity in special and general relativity. As tests and example applications, we resolve the shock substructure due to the presence of dissipation, and report on relativistic versions of the magneto-thermal instability and heat flux driven buoyancy instability, which arise due to anisotropic heat conduction, and of the firehose instability, which occurs due to anisotropic pressure (i.e., viscosity). Finally, we show an example integration of an accretion flow around a Kerr black hole, using Extended GRMHD.

79 ASTRONOMY AND ASTROPHYSICS↗

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.↗

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.↗

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.↗

Proton imaging of stochastic magnetic fields

Recent laser-plasma experiments report the existence of dynamically significant magnetic fields, whose statistical characterization is essential for a complete understanding of the physical processes these experiments are attempting to investigate. In this paper, we show how a proton-imaging diagnostic can be used to determine a range of relevant magnetic-field statistics, including the magnetic-energy spectrum. To achieve this goal, we explore the properties of an analytic relation between a stochastic magnetic field and the image-flux distribution created upon imaging that field. This ‘Kugland image-flux relation’ was previously derived under simplifying assumptions typically valid in actual proton-imaging set-ups. We conclude that, as with regular electromagnetic fields, features of the beam’s final image-flux distribution often display a universal character determined by a single, field-scale dependent parameter – the contrast parameter $\unicode[STIX]{x1D707}\equiv d_{s}/{\mathcal{M}}l_{B}$ – which quantifies the relative size of the correlation length $l_{B}$ of the stochastic field, proton displacements $d_{s}$ due to magnetic deflections and the image magnification ${\mathcal{M}}$ . For stochastic magnetic fields, we establish the existence of four contrast regimes, under which proton-flux images relate to their parent fields in a qualitatively distinct manner. These are linear, nonlinear injective, caustic and diffusive. The diffusive regime is newly identified and characterized. The nonlinear injective regime is distinguished from the caustic regime in manifesting nonlinear behaviour, but as in the linear regime, the path-integrated magnetic field experienced by the beam can be extracted uniquely. Thus, in the linear and nonlinear injective regimes we show that the magnetic-energy spectrum can be obtained under a further statistical assumption of isotropy. This is not the case in the caustic or diffusive regimes. We discuss complications to the contrast-regime characterization arising for inhomogeneous, multi-scale stochastic fields, which can encompass many contrast regimes, as well as limitations currently placed by experimental capabilities on one’s ability to extract magnetic-field statistics. Furthermore, the results presented in this paper are of consequence in providing a comprehensive description of proton images of stochastic magnetic fields, with applications for improved analysis of proton-flux images.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A plasma deflagration accelerator as a platform for laboratory astrophysics

The replication of astrophysical flows in the laboratory is critical for isolating particular phenomena and dynamics that appear in complex, highly-coupled natural systems. Particularly, plasma jets are observed in astrophysical contexts at a variety of scales, typically at high magnetic Reynolds number and driven by internal currents. In this paper, we present detailed measurements of the plasma parameters within deflagration-produced plasma jets, the scaling of these parameters against both machine operating conditions and the corresponding astrophysical phenomena. Using optical and spectroscopic diagnostics, including Schlieren cinematography, we demonstrate the production of current-driven plasma jets of »100 km/s and magnetic Reynolds numbers of »100, and discuss the dynamics of their acceleration into vacuum. Finally, the results of this study will contribute to the reproduction of various types of astrophysical jets in the laboratory and indicate the ability to further probe active research areas such as jet collimation, stability, and interaction.

79 ASTRONOMY AND ASTROPHYSICS↗

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↗

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↗

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.↗