Search NASA⌕ Search

SEARCH · Search NASA

Results for “hybrid particle-in-cell code”

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.

40 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM SciDAC) (Technical Final Report)

Runaway electrons can severely damage the plasma facing components on ITER during a major disruption and pose a major risk for tokamak fusion. It has been recognized that an adequate disruption mitigation system (DMS) is essential for the safe operation of ITER. The United States is responsible for the design and implementation of the disruption mitigation system on ITER, and in July 2016 the Simulation Center for Runaway Electron Avoidance and Mitigation (SCREAM) was launched by DOE, in a joint Fusion Energy Sciences (FES) and Advanced Scientific Computing Research (ASCR) collaboration. SCREAM was a comprehensive theory and simulation SciDAC center that provided physics guidance in the avoidance and mitigation of runaway electrons, and in tandem with domestic and international experiments, helped establish the qualitative and quantitative bases for safe operational scenarios and viable mitigation techniques. The SCREAM center assembled a national team of experts in runaway electron physics, tokamak disruptions, magnetohydrodynamic (MHD) simulation, and advanced algorithms and computing. The team combined advanced simulation and analysis capability facilitated by direct participation of ASCR SciDAC institutes with theoretical models and code development by FES scientists to focus on the runaway risk for ITER and tokamaks in general. The research scope was focussed on integrated simulations of kinetic runaway electrons, including MHD and fluid models of impurity transport, within a research plan guided by theory. The specific research tasks were (1) establish the fundamental physics of runaway generation, saturation, and dynamical evolution in a tokamak; (2) examine the critical path toward runaway avoidance; and (3) investigate the viability and effectiveness of the leading candidate schemes for runaway mitigation. In all three areas, members of the team carried out scoping studies that established the readiness for rapid and critical advances, especially in the deployment and further development of large-to extreme-scale simulation tools. Our multi-pronged computational approach included (1) relativistic Fokker-Planck solvers with discretization in phase space, (2) self-consistent particle-in-cell techniques, (3) particle-based Monte-Carlo, and (4) MHD-particle hybrid simulations. Cross-check between these different methods provided an additional means for verification and further bolstered the fidelity of our physics prediction. Validation against experimental results brings confidence to the predictive capability for ITER and frequently leads to new ideas for understanding and mitigating the thermal quench driven runaway electron phenomenon.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Lower hybrid resonance acceleration of electrons and ions in solar flares and the associated microwave emission

The particle acceleration processes here studied are driven by the relaxation of unstable ion ring distributions; these produce strong wave activity at the lower hybrid resonance frequency which collapses, and forms energetic electron and ion tails. The results obtained are applied to the problem posed by the production of energetic particles by solar flares. The numerical simulation results thus obtained by a 2 1/2-dimensional particle-in-cell code show a simultaneous acceleration of electrons to 10-500 keV energies, and of ions to as much as the 1 MeV range; the energy of the latter is still insufficient to account for gamma-ray emission in the 4-6 MeV range, but furnish a seed population for further acceleration.

Mcclements, K. G.↗
NASA NTRS1997

Physics of Rotation Powered Pulsars and Their Nebulae

This letter is my progress report for the Astrophysics Theory grant. The first goal of the research supported by this grant is to produce a time dependent theory of the unsteady relativistic collisionless shock wave terminating the relativistic wind from a pulsar, and compare the predicted surface brightness fluctuations to Hubble Space Telescope observations of the wisps in the Crab Nebula. The second goal is to model the production of electron-positron pairs over the polar caps of rotation powered pulsars, and use the results to predict the heating of the surface due to particle trapping and bombardment of the atmosphere at the base of the polar field lines. We have succeeded in creating a one dimensional hybrid code, in which the electron-positron pairs incident on the shock structure are modeled as a relativistic, adiabatic, ideal MHD fluid, while the heavy ions are treated as particles using a particle-in-cell algorithm. The electromagnetic fields are evaluated from the currents and charge densities in the pairs and the ions, while the particles and the fluid accelerate in response to the computed self-consistent electromagnetic fields. The results are promising, in that the underlying ion cyclotron instability generates finite amplitude, propagating magnetosonic waves in the pairs, whose wavelengths and frequencies, when translated into physical units, are comparable to the observed running waves of brightness observed by HST near the Crab pulsar. The code is undergoing a number of tests, to assure us that this preliminary correspondence is not an artifact. In the coming year, the observational appearance of the models will be computed and compared to the HST observations of the Crab now in hand, and used to predict the HST results which will be obtained the year after next. WE also developed a one dimensional cascade theory for pair creation over pulsars' polar caps. A linear integral equation describing the synchrotron cascade has been derived and solved by iterative techniques, in the case when a high energy electron moving parallel to a star centered dipole magnetic field initiates the cascade through curvature gamma ray emission.

Arons, Jonathan↗