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,297 records · Page 72

On the small divergence of laser-driven ion beams from nanometer thick foils

We report on experimental studies of divergence of proton beams from nanometer thick diamond-like carbon foils irradiated by a linearly polarized intense laser with high contrast. Proton beams with extremely small divergence (half angle) of 2° are observed in addition with a remarkably well-collimated feature over the whole energy range, showing one order of magnitude reduction of the divergence angle in comparison to the results from μm thick targets. Similar features are reproduced in two-dimensional particle-in-cell simulations with parameters representing our experiments, indicating a strong influence from the electron density distribution on the divergence of protons. Our comprehensive experimental study reveals grand opportunities for using nm foils in experiments that require high ion flux and small divergence.

43 PARTICLE ACCELERATORS↗

Current sheets and pressure anisotropy in the reconnection exhaust

A particle-in-cell simulation shows that the exhaust during anti-parallel reconnection in the collisionless regime contains a current sheet extending 100 inertial lengths from the X line. The current sheet is supported by electron pressure anisotropy near the X line and ion anisotropy farther downstream. Field-aligned electron currents flowing outside the magnetic separatrices feed the exhaust current sheet and generate the out-of-plane, or Hall, magnetic field. Existing models based on different mechanisms for each particle species provide good estimates for the levels of pressure anisotropy. The ion anisotropy, which is strong enough to reach the firehose instability threshold, is also important for overall force balance. It reduces the outflow speed of the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Impact of Solar Array Designs on High Voltage Operations

As power levels of advanced spacecraft climb above 25 kW, higher solar array operating voltages become attractive. Even in today s satellites, operating spacecraft buses at 100 V and above has led to arcing in GEO communications satellites, so the issue of spacecraft charging and solar array arcing remains a design problem. In addition, micrometeoroid impacts on all of these arrays can also lead to arcing if the spacecraft is at an elevated potential. For example, tests on space station hardware disclosed arcing at 75V on anodized A1 structures that were struck with hypervelocity particles in Low Earth Orbit (LEO) plasmas. Thus an understanding of these effects is necessary to design reliable high voltage solar arrays of the future, especially in light of the Vision for Space Exploration of NASA. In the future, large GEO communication satellites, lunar bases, solar electric propulsion missions, high power communication systems around Mars can lead to power levels well above 100 kW. As noted above, it will be essential to increase operating voltages of the solar arrays well above 80 V to keep the mass of cabling needed to carry the high currents to an acceptable level. Thus, the purpose of this paper is to discuss various solar array approaches, to discuss the results of testing them at high voltages, in the presence of simulated space plasma and under hypervelocity impact. Three different types of arrays will be considered. One will be a planar array using thin film cells, the second will use planar single or multijunction cells and the last will use the Stretched Lens Array (SLA - 8-fold concentration). Each of these has different approaches for protection from the space environment. The thin film cell based arrays have minimal covering due to their inherent radiation tolerance, conventional GaAs and multijunction cells have the traditional cerium-doped microsheet glasses (of appropriate thickness) that are usually attached with Dow Corning DC 93-500 silicone adhesive. In practice, these cover glasses and adhesive do not cover the cell edges. Finally, in the SLA, the entire cell and cell edges are fully encapsulated by a cover glass that overhangs the cell perimeter and the silicone adhesive covers the cell edges providing a sealed environment. These three types of blanket technology have been tested at GRC and Auburn. The results of these tests will be described. For example, 15 modules composed of four state-of-the-art 2x4 cm GaAs solar cells with 150 pm cover glasses connected in two-cell series strings were tested at high voltage, in plasma under hypervelocity impact. A picture of one of the modules is shown in figure 1. These were prepared by standard industry practice from a major supplier and had efficiencies above 18%. The test results and other fabrication factors that influenced the tests will be presented. In addition, results for SLA segments tested under the same conditions will be presented. Testing of thin film blankets at GRC will also be presented. Figure 1 : Typical GaAs Solar Cell Module These results will show significant differences in resistance to arcing that are directly related to array design and manufacturing procedures. Finally, the approaches for mitigating the problems uncovered by these tests will be described. These will lay the foundation for future higher voltage array operation, even including voltages above 300-600 V for direct drive SEP applications.

Brandhorst, Henry W., Jr.↗

Retarding potential analyzer for the Pioneer-Venus Orbiter Mission

The retarding potential analyzer on the Pioneer-Venus Orbiter Mission has been designed to measure most of the thermal plasma parameters within and near the Venusian ionosphere. Parameters include total ion concentration, concentrations of the more abundant ions, ion temperatures, ion drift velocity, electron temperature, and low-energy (0-50 eV) electron distribution function. To accomplish these measurements on a spinning vehicle with a small telemetry bit rate, several functions, including decision functions not previously used in RPA's, have been developed and incorporated into this instrument. The more significant functions include automatic electrometer ranging with background current compensation; digital, quadratic retarding potential step generation for the ion and low-energy electron scans; a current sampling interval of 2 ms throughout all scans; digital logic inflection point detection and data selection; and automatic ram direction detection. Extensive numerical simulation and plasma chamber tests have been conducted to verify adequacy of the design for the Pioneer Mission.

Knudsen, W. C.↗

Comparison for non-local hydrodynamic thermal conduction models

Inertial confinement fusion and specifically shock ignition involve temperatures and temperature gradients for which the classical Spitzer-Haerm thermal conduction breaks down and a non-local operator is required. In this article, two non-local electron thermal conduction models are tested against kinetic Vlasov-Fokker-Planck simulations. Both models are shown to reproduce the main features of thermal heat front propagation at kinetic timescales. The reduction of the thermal conductivity as a function of the scalelength of the temperature gradient is also recovered. Comparisons at nanosecond timescales show that the models agree on the propagation velocity of the heat front, but major differences appear in the thermal precursor.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Cascaded target normal sheath acceleration

A cascaded target normal sheath acceleration (TNSA) scheme is proposed to simultaneously increase energy and improve energy spread of a laser-produced mono-energetic proton beam. An optimum condition that uses the maximum sheath field to accelerate the center of the proton beam is theoretically found and verified by two-dimensional particle-in-cell simulations. An initial 10 MeV proton beam is accelerated to 21 MeV with energy spread decreased from 5% to 2% under the optimum condition during the process of the cascaded TNSA. The scheme opens a way to scale proton energy lineally with laser energy.

43 PARTICLE ACCELERATORS↗

Acceleration and stability of a high-current ion beam in induction fields

A one-dimensional nonlinear analytic theory of the filamentation instability of a high-current ion beam is formulated. The results of 2.5-dimensional numerical particle-in-cell simulations of acceleration and stability of an annular compensated ion beam (CIB) in a linear induction particle accelerator are presented. It is shown that additional transverse injection of electron beams in magnetically insulated gaps (cusps) improves the quality of the ion-beam distribution function and provides uniform beam acceleration along the accelerator. The CIB filamentation instability in both the presence and the absence of an external magnetic field is considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Optimization of laser-target interaction for proton acceleration

The aim of this study is to optimize the characteristics, in particular, the maximum energy and the number of high energy protons accelerated by a high intensity laser from different targets by using two- and three-dimensional particle-in-cell simulations. Two principal ways are considered to increase the maximum proton energy: (1) the use of an optimized target design (by choosing its atomic composition, density, structure, thickness, and transverse size); (2) variation of the laser pulse parameters (duration, power, intensity, focal spot size, polarization, pulse shape, etc.). Our analysis demonstrates a possibility to accelerate protons to the energies exceeding 200 MeV with 20 J laser pulses with appropriately chosen focusing and duration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Influence of axial self-magnetic field component on arcing behavior of spiral-shaped contacts

The transverse magnetic field (TMF) contact design is commonly used in vacuum interrupters. When arcing occurs between the TMF contacts, the contact structure can create a self-induced magnetic field that drives the arc to move and rotate on the contact, and thus local overheating and severe erosion can be avoided. However, TMF contacts could also create an axial self-magnetic component, and the influence of this component on the arc behavior has not been considered to date. In this paper, five different types of Cu-Cr spiral-shaped TMF contacts with three different structures are investigated in a demountable vacuum chamber that contains a high-speed charge-coupled device video camera. It was found that the contact structure greatly influenced the arc behavior, especially in terms of arc rotation and the effective contact area, while contacts with the same slot structure but different diameters showed similar arc behavior and arc motion. The magnetic field distribution and the Lorentz force of each of the three different contact structures are simulated, and the axial self-magnetic field was first taken into consideration for investigation of the TMF contact design. It was found that contact designs that have higher axial self-magnetic field components tend to have arc columns with larger diameters and show poorer arc motion and rotation performance in the experiments.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Multicomponent kinetic simulation of Bernstein–Greene–Kruskal modes associated with ion acoustic and dust-ion acoustic excitations in electron-ion and dusty plasmas

A series of numerical simulations based on a recurrence-free Vlasov kinetic algorithm presented earlier [Abbasi et al., Phys. Rev. E 84, 036702 (2011)] are reported. Electron-ion plasmas and three-component (electron-ion-dust) dusty, or complex, plasmas are considered, via independent simulations. Considering all plasma components modeled through a kinetic approach, the nonlinear behavior of ionic scale acoustic excitations is investigated. The focus is on Bernstein–Greene–Kruskal (BGK) modes generated during the simulations. In particular, we aim at investigating the parametric dependence of the characteristics of BGK structures, namely of their time periodicity (τ{sub trap}) and their amplitude, on the electron-to-ion temperature ratio and on the dust concentration. In electron-ion plasma, an exponential relation between τ{sub trap} and the amplitude of BGK modes and the electron-to-ion temperature ratio is observed. It is argued that both characteristics, namely, the periodicity τ{sub trap} and amplitude, are also related to the size of the phase-space vortex which is associated with BGK mode creation. In dusty plasmas, BGK modes characteristics appear to depend on the dust particle density linearly.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Intermittency in Hall-magnetohydrodynamics with a strong guide field

We present a detailed study of intermittency in the velocity and magnetic field fluctuations of compressible Hall-magnetohydrodynamic turbulence with an external guide field. To solve the equations numerically, a reduced model valid when a strong guide field is present is used. Different values for the ion skin depth are considered in the simulations. The resulting data are analyzed computing field increments in several directions perpendicular to the guide field, and building structure functions and probability density functions. In the magnetohydrodynamic limit, we recover the usual results with the magnetic field being more intermittent than the velocity field. In the presence of the Hall effect, field fluctuations at scales smaller than the ion skin depth show a substantial decrease in the level of intermittency, with close to monofractal scaling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Integrated System and Application Continuous Performance Monitoring and Analysis Capability

Scientific applications run on high-performance computing (HPC) systems are critical for many national security missions within Sandia and the NNSA complex. However, these applications often face performance degradation and even failures that are challenging to diagnose. To provide unprecedented insight into these issues, the HPC Development, HPC Systems, Computational Science, and Plasma Theory & Simulation departments at Sandia crafted and completed their FY21 ASC Level 2 milestone entitled "Integrated System and Application Continuous Performance Monitoring and Analysis Capability." The milestone created a novel integrated HPC system and application monitoring and analysis capability by extending Sandia's Kokkos application portability framework, Lightweight Distributed Metric Service (LDMS) monitoring tool, and scalable storage, analysis, and visualization pipeline. The extensions to Kokkos and LDMS enable collection and storage of application data during run time, as it is generated, with negligible overhead. This data is combined with HPC system data within the extended analysis pipeline to present relevant visualizations of derived system and application metrics that can be viewed at run time or post run. This new capability was evaluated using several week-long, 290-node runs of Sandia's ElectroMagnetic Plasma In Realistic Environments ( EMPIRE ) modeling and design tool and resulted in 1TB of application data and 50TB of system data. EMPIRE developers remarked this capability was incredibly helpful for quickly assessing application health and performance alongside system state. In short, this milestone work built the foundation for expansive HPC system and application data collection, storage, analysis, visualization, and feedback framework that will increase total scientific output of Sandia's HPC users.

97 MATHEMATICS AND COMPUTING↗

Integrated System and Application Continuous Performance Monitoring and Analysis Capability (Final)

Scientific applications run on high-performance computing (HPC) systems are critical for many national security missions within Sandia and the NNSA complex. However, these applications often face performance degradation and even failures that are challenging to diagnose. To provide unprecedented insight into these issues, the HPC Development, HPC Systems, Computational Science, and Plasma Theory & Simulation departments at Sandia crafted and completed their FY21 ASC Level 2 milestone entitled "Integrated System and Application Continuous Performance Monitoring and Analysis Capability." The milestone created a novel integrated HPC system and application monitoring and analysis capability by extending Sandia’s Kokkos application portability framework, Lightweight Distributed Metric Service (LDMS) monitoring tool, and scalable storage, analysis, and visualization pipeline. The extensions to Kokkos and LDMS enable collection and storage of application data during run time, as it is generated, with negligible overhead. This data is combined with HPC system data within the extended analysis pipeline to present relevant visualizations of derived system and application metrics that can be viewed at run time or post run. This new capability was evaluated using several week-long, 290-node runs of Sandia’s ElectroMagnetic Plasma In Realistic Environments (EMPIRE) modeling and design tool and resulted in 1TB of application data and 50TB of system data. EMPIRE developers remarked this capability was incredibly helpful for quickly assessing application health and performance alongside system state. In short, this milestone work built the foundation for expansive HPC system and application data collection, storage, analysis, visualization, and feedback framework that will increase total scientific output of Sandia’s HPC users.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-dimensional particle-in-cell modeling of terahertz gyrotrons with cylindrical and planar configurations of the interaction space

We perform 3D particle-in-cell simulations of terahertz gyrotrons with two different configurations of the interaction space. For a gyrotron with conventional cylindrical configuration of the interaction cavity, we demonstrate reasonable agreement between simulations and experimental results, including output frequency, structure of the higher-order operating mode (TE{sub 17,4}), output power, and ohmic losses. For a novel planar gyrotron scheme with transverse energy extraction, a possibility of further increasing the oversized factor with the single-mode operation regime retained is shown. Frequency detuning by mechanical variation of the gap between waveguide plates is also demonstrated.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Enhanced high harmonic generation and the phase effect in double-sided relativistic laser-foil interaction

High harmonic generation (HHG) from relativistic laser-foil interaction is investigated analytically and through particle-in-cell simulations. Previous work has shown that when two counter-propagating circularly polarized (CP) laser pulses interact with a thin foil, electrons can be well confined spatially to form a high density layer. The layer electrons oscillate in certain transversal direction and radiate intense high order harmonics. It is demonstrated here that there is a critical foil thickness, only below which can high harmonics be generated efficiently. Furthermore, to enhance the intensity in higher order region, the third linearly polarized (LP) short-wavelength laser pulse with much lower intensity is introduced. Analysis and simulations both show that the enhancement is determined by the relative phase {delta}{phi} between the driving CP laser pulses and LP pulse. The enhancement at high order is quite considerable and very sensitive to the relative phase {delta}{phi}, thus offering not only a way to efficiently produce HHG but also a new method to measure the phase of intense high-frequency laser pulses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Nonlinear evolution of the Kelvin-Helmholtz instability in the double current sheet configuration

The nonlinear evolution of the Kelvin-Helmholtz (KH) instability driven by a radially antisymmetric shear flow in the double current sheet configuration is numerically investigated based on a reduced magnetohydrodynamic model. Simulations reveal different nonlinear fate of the KH instability depending on the amplitude of the shear flow, which restricts the strength of the KH instability. For strong shear flows far above the KH instability threshold, the linear electrostatic-type KH instability saturates and achieves a vortex flow dominated quasi-steady state of the electromagnetic (EM) KH turbulence with large-amplitude zonal flows as well as zonal fields. The magnetic surfaces are twisted significantly due to strong vortices but without the formation of magnetic islands. However, for the shear flow just over the KH instability threshold, a weak EM-type KH instability is saturated and remarkably damped by zonal flows through modifying the equilibrium shear flow. Interestingly, a secondary double tearing mode (DTM) is excited subsequently in highly damped KH turbulence, behaving as a pure DTM in a flowing plasma as described in Mao et al. [Phys. Plasmas 21, 052304 (2014)]. However, the explosive growth phenomenon is replaced by a gradually growing oscillation due to the extremely twisted islands. As a result, the release of the magnetic energy becomes slow and the global magnetic reconnection tends to be gentle. A complex nonlinear interaction between the EM KH turbulence and the DTMs occurs for the medium shear flows above the KH instability threshold, turbulent EM fluctuations experience oscillatory nonlinear growth of the DTMs, finally achieves a quasi-steady state with the interplay of the fluctuations between the DTMs and the EM KH instability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PPPL-5081

This paper formulates the Tokamak Magneto-Hydrodynamics (TMHD), initially outlined by X. Li and L.E. Zakharov [Plasma Science and Technology, accepted, ID:2013-257 (2013)] for proper simulations of macroscopic plasma dynamics. The simplest set of magneto-hydrodynamics equations, sufficient for disruption modeling and extendable to more refined physics, is explained in detail. First, the TMHD introduces to 3-D simulations the Reference Magnetic Coordinates (RMC), which are aligned with the magnetic field in the best possible way. The numerical implementation of RMC is adaptive grids. Being consistent with the high anisotropy of the tokamak plasma, RMC allow simulations at realistic, very high plasma electric conductivity. Second, the TMHD splits the equation of motion into an equilibrium equation and the plasma advancing equation. This resolves the 4 decade old problem of Courant limitations of the time step in existing, plasma inertia driven numerical codes. The splitting allows disruption simulations on a relatively slow time scale in comparison with the fast time of ideal MHD instabilities. A new, efficient numerical scheme is proposed for TMHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗