Search NASA⌕ Search

SEARCH · Search NASA

Results for “PICS”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Collaborative Research: Particle Simulation of Collisionless Magnetic Reconnection Under Finite Guide Field (Final Technical Report)

The project investigated 3-D linear instabilities of Harris current sheet under a realistic ion-to-electron mass ratio and a finite guide field. The excitation and 3-D properties of lower hybrid drift instability (LHDI), along with two new instabilities were demonstrated. The report presents the nonlinear properties of current sheet instabilities by employing Gyrokinetic electron and Fully kinetic ion (GeFi) particle simulation code. The GeFi simulation results are compared with fully kinetic PIC simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Final Report for Center for Tokamak Transient Simulations at USU

Providing plasma fluid codes like NIMROD with continuum drift kinetic (CDK) physics that is quantitatively valid and computationally feasible throughout the spatial domain is difficult. Work at Utah State University (USU), in collaboration with the Center for Tokamak Transient Simulations (CTTS), focused on applying CDK closures in disruption-related calculations. Three examples where kinetic physics is paramount are (1) the electron stress tensor closure in Ohms law for accurately describing neoclassical tearing mode (NTM) evolution, (2) runaway electron (RE) density (nRE) and current (jRE) moments in NIMROD’s extended MHD model for self-consistent evolution of RE populations during disruptions and, (3) energetic ion effects on a myriad of MHD instabilities. While NTM simulations and continuum and PIC approaches to energetic ions in NIMROD have been a major goals of USU’s closure work for several years, the development of self-consistent CDK RE capability in NIMROD was started and extended considerably during the CTTS effort. Some goals of CDK RE in NIMROD are to explore the effects of the 2D relativistic phase space in 4D simulations and compare with NIMROD’s fluid RE model. Four publications and two PhD theses came out of the USU CTTS effort.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Developing nonlinear laser-plasma instability models for high-fidelity, multi-physics simulation capability for ICF/HED (IC Report) [Slides]

Modeling nonlinear effects of laser plasma instabilities (LPI) is critically important for ICF/HED experiments. Our project has the following goals: • Develop physics-based nonlinear LPI models using PIC simulations • Couple the nonlinear LPI effects to macroscopic modeling of ICF/HED exp. through laser ray-tracing (LRT).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Ion Channel Laser Based on Direct Laser Acceleration of a Shaped Beam Driver (Final Technical Report)

The main focus of our research was to understand the synergies between an electron bunch and a laser pulse when the two co-propagate through the plasma. We have clearly demonstrated using theoretical and computational modeling that the propagation distance of both the bunch and the laser pulse could be extended. The ability of the combined bunch/pulse system to extend the propagation distance and the size of the plasma bubble enables highly-efficient sources of relativistic electrons. As pointed out by the Plasma Decadal Study, such electron sources can be used for generating extremely bright x-rays for a variety of applications can serve as probes and diagnostics for other plasma experiments: high energy density (HED) sciences, inertial confinement fusion (ICF), and potentially Fusion Materials and Technology (FM&T). Development of advanced diagnostics of plasma-based accelerators is yet another key area identified by recent reports. Even broader security and medical applications of compact accelerator-based radiation sources, such as very high energy electron (VHEE) sources for FLASH radiobiology, have been identified by a recent multi-agency panel. A number of key technical issues were considered and successfully resolved during the course of the grant. Those include: beam loading effect of both the driver and witness bunches on the wake, laser channeling by the bunch, and the direct laser acceleration (DLA) of the driver bunch by the laser pulse. To demonstrate clear synergy, we were able to ascertain that the total energy gain of a witness bunch in the wake of the combined bunch/laser complex is large than the sum of the energy gains in the wake of the laser pulse alone, and the driver bunch alone. We have also demonstrated that not only the energy gain is improved, but the energy spread is not sacrificed. To carry out these simulations, we have to develop a range of in-house computational tools, including a fully-3D code WAND-PIC.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Investigation on laser absorption and x-ray radiation in microstructured titanium targets heated by short-pulse relativistic laser pulses

The enhancement effect of a microstructured surface on laser absorption and characteristic Kα emission has been investigated by measuring K-shell emission from titanium (Ti) targets irradiated with high-intensity (~10 20 W cm -2 ), subpicosecond (500 fs) laser pulses. The experimental results indicate a modest enhancement (1.6×) of Kα emission from microstructured targets compared to flat foils, but with a similar intensity and profile of Heα and Li-like satellites. Particle-in-cell (PIC) simulations are implemented to further understand the underlying physical processes in the laser interaction with both targets, interpreting the mechanisms responsible for the Kα enhancement. The reasons for the lower-than-expected enhancement of Kα emission are discussed. The rapid heating of the bulk plasma might result in the premature shutdown of Kα emission before the thermalization of hot electrons or even the end of laser pulses, suggesting that the use of Kα emission as a diagnostic of the hot-electron yield or relaxation could lead to a misinterpretation.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Linac_Gen: Integrating Machine Learning and Particle-in-Cell Methods for Enhanced Beam Dynamics at Fermilab

Here, we introduce Linac_Gen, a tool developed at Fermilab, which combines machine learning algorithms with Particle-in-Cell methods to advance beam dynamics in linacs. Linac_Gen employs techniques such as Random Forest, Genetic Algorithms, Support Vector Machines, and Neural Networks, achieving a tenfold increase in speed for phase-space matching in Linacs over traditional methods, through the use of genetic algorithms. Crucially, Linac_Gen's adept handling of 3D field maps elevates the precision and realism in simulating beam instabilities and resonances, marking a key advancement in the field. Benchmarked against established codes, Linac_Gen demonstrates not only improved efficiency and precision in beam dynamics studies but also in the design and optimization of Linac systems, as evidenced in its application to Fermilab's PIP-II Linac project. This work represents a notable advancement in accelerator physics, marrying ML with PIC methods to set new standards for efficiency and accuracy in accelerator design and research. Linac_Gen exemplifies a novel approach in accelerator technology, offering substantial improvements in both theoretical and practical aspects of beam dynamics.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

MFiX Development Updates

This presentation discusses recent developments of the Multiphase Flow with Interphase eXchanges (MFiX) software. A brief review of all modeling approaches along with their cost versus accuracy is provided to guide users when selecting a model for a given application. The main new features of the past five official releases of MFiX are described. Improvement in chemistry management allow for easier simulation setup, and faster simulation speed. Progress in the implementation of a thin-wall boundary condition is presented. The major new model development over the past year is the release of the Glued Sphere Particle model (GSP), where component spheres are combined together to represent non-spherical particles. The integration of Machine Learning (ML) workflow in the CFD process is discussed with two applications: a surrogate model for stiff chemistry and the development of a PIC stress model using ML.

Dietiker, Jeff↗
DOE OSTI2025

New model for the ion collection by cylindrical probes over a wide range of collisionality

Langmuir probes remain one of the most important diagnostic tools for plasma processing applications. Modern probe analysis usually relies on the electron current part of the Langmuir probe characteristic using the Druyvesteyn method. However, for electronegative plasmas or for discharges containing dust the analysis of the ion current attracted by the probe can be desirable to determine the ion density. But, even at low pressures of a few Pa, the ion current is affected by collisions due to the large cross section for charge exchange. Available theories for collisional or collision-enhanced ion currents onto probes are complex and not well validated. Thus, in this contribution, we compare available collisional probe theories for the ion current to results of particle-in-cell (PIC) simulations. To this end, the probe surrounded by a semi-infinite plasma is simulated using a modified version of the open-source code EDIPIC. A dataset of currents for different neutral gas pressures is obtained and compared to the different theories from the literature. Based on these results, we propose a simpler and more intuitive model for the ion current collected by the probe, based on the model of Gatti and Kortshagen (Phys. Rev. E 78, 046402, 2008), developed for the charging of dust particles.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Digital Twin Framework for PIP-II Linac: AI-Driven Multi-Scale Modeling from Ion Source to 800 MeV

The PIP-II linac will enable >1.2 MW beam power for DUNE, requiring unprecedented operational reliability across its warm front-end (RFQ, MEBT) and five distinct SRF sections operating at 162.5/325/650 MHz. We present a comprehensive digital twin framework uniquely combining a fully differentiable fast beam transport code with neural network surrogates trained on high-fidelity PIC simulations, capturing space charge and nonlinear dynamics beyond traditional envelope codes while achieving 10⁴ speedup at <1% accuracy. End-to-end differentiability enables gradient-based optimization across 500+ parameters simultaneously previously impossible with conventional tools while the model incorporates static/dynamic errors and serves as a virtual commissioning platform for diverse hardware integration. The framework facilitates reinforcement learning for pulsed/CW mode transitions, predictive maintenance through anomaly detection, and autonomous tuning algorithm development with real-time execution capability. Validation against physics simulations shows excellent agreement for the front-end, with initial results demonstrating potential for 30% commissioning time reduction and proactive fault mitigation, providing a scalable blueprint for operating next-generation high-intensity accelerators.

Pathak, Abhishek [Fermilab] (ORCID:000000021704208↗
DOE OSTI2026

Collaborative Research: Enabling multi-scale studies of magnetic reconnection with interpretable data-driven models

The development of accurate reduced descriptions and improved closures for magnetic reconnection is an important and a long‐standing challenge in plasma physics. The four‐fluid approach, and associated closures, that were investigated have the potential to improve the accuracy of plasma fluid models, capturing physical effects which would otherwise require a kinetic description. If successful, this approach could have an important impact for the modeling of laboratory and space plasmas. The major goals of this project were to develop new machine learning (ML) tools based on sparse and symbolic regression techniques, and to extract interpretable and generalizable reduced models (e.g., in the form of partial differential equations - PDEs) from data generated by first principles plasma simulations. Preserving interpretability of such data‐driven models is key to addressing the long‐standing theoretical and numerical challenges. Prior proof‐of‐principle studies have demonstrated the enormous potential of this approach, by recovering the well‐established hierarchy of plasma equations (from Vlasov to MHD) from data produced by particle‐in‐cell (PIC) simulations. Our goal in this project was to extend and apply these new tools to construct better kinetic closures for magnetic reconnection; to derive better models of particle injection and acceleration by this fundamental plasma process; and to use this understanding to accelerate the development of multi‐scale plasma algorithms. While our immediate focus was on the problem of magnetic reconnection, the tools that were will developed are general and applicable to other areas of plasma physics, and more broadly to many‐body phenomena. We anticipate that the development of these multi‐scale models will have a significant impact across different areas of plasma science, from fusion to space and astrophysical plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

A kinetic study of fusion burn waves in compressed deuterium–tritium and proton–boron plasmas

We present particle-in-cell simulations with Monte Carlo collisions of fusion burn waves in compressed deuterium–tritium and proton–boron plasmas. We study the energy balance in the one-dimensional expansion of a hot-spot by simulating Coulomb collisions, fusion reactions, and bremsstrahlung emission with a Monte Carlo model and inverse bremsstrahlung absorption using a new PIC model. This allows us to self-consistently capture the alpha particle heating and radiative losses in the expanding hot-spot and surrounding cold fuel. After verifying our model in a code-to-code comparison with both kinetic and fluid codes for the case of a deuterium–tritium hot-spot, we simulate the expansion of a proton–boron hot-spot initialized at 200 keV and 1,000 g/cm 3 . Our model predicts that energy radiated by the hot-spot is recaptured by the surrounding high-density opaque fuel reducing the expansion work done by the propagating burn wave. As a result, we find the net fusion energy produced over the course of $20$~ps is twice the initial hot-spot energy independent of whether radiation physics is included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Refining the Census of the Upper Scorpius Association with Gaia

We have refined the census of stars and brown dwarfs in the Upper Sco association (∼10 Myr, ∼145 pc) by (1) updating the selection of candidate members from our previous survey to include the high-precision astrometry from the second data release of Gaia, (2) obtaining spectra of a few hundred candidate members to measure their spectral types and verify their youth, and (3) assessing the membership (largely with Gaia astrometry) of 2020 stars toward Upper Sco that show evidence of youth in this work and previous studies. We arrive at a catalog of 1761 objects that are adopted as members of Upper Sco. The distribution of spectral types among the adopted members is similar to those in other nearby star-forming regions, indicating a similar initial mass function. In previous studies, we have compiled mid-infrared photometry from the Wide-field Infrared Survey Explorer and the Spitzer Space Telescope for members of Upper Sco and used those data to identify the stars that show evidence of circumstellar disks; we present the same analysis for our new catalog of members. As in earlier work, we find that the fraction of members with disks increases with lower stellar masses, ranging from ≲10% for >1 M {sub ⊙} to ∼22% for 0.01–0.3 M {sub ⊙}. Finally, we have estimated the relative ages of Upper Sco and other young associations using their sequences of low-mass stars in M{sub G{sub R{sub P}}} versus G {sub BP} − G {sub RP}. This comparison indicates that Upper Sco is a factor of two younger than the β Pic association (21–24 Myr) according to both nonmagnetic and magnetic evolutionary models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Spectroscopic Orbits of Subsystems in Multiple Stars. VII

Spectroscopic orbits of main-sequence stars HIP 3150A, 6873B, 11537A, 22531A, 22534B, 31089B, 49336A, 104833C, and 107731A belonging to eight multiple systems are determined from high-resolution spectra taken with CHIRON. Two of those are twins with mass ratios above 0.95. HIP 11537 is a young three-tier quadruple system with inner periods of 22.3 and 1146 days and the outer period of 3 kyr. HIP 22531 (ι Pic) is the brightest star in a hierarchical system with six components. It is a spectroscopic binary with periods of 1.56 days and 2.75 yr, as well as a γ Dor variable with a period of 0.67 day, possibly in a 7:3 resonance with the inner orbit. HIP 22534, also member of this system, is a double-lined binary with a period of 208 days. For HIP 31089, both the spectro-interferometric 32 yr outer orbit and the 213 day orbit of the subsystem are determined. HIP 107731 is a triple system with an inner period of 470 days and a fast spatial motion, likely metal-poor. New orbits contribute to the statistics of hierarchical multiplicity in the solar neighborhood.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

A Census of the 32 Ori Association with Gaia*

I have used high-precision photometry and astrometry from the third data release of Gaia (DR3) to identify candidate members of the 32 Ori association. Spectral types and radial velocities have been measured for subsets of the candidates using new and archival spectra. For the candidates that have radial velocity measurements, I have used UVW velocities to further constrain their membership, arriving at a final catalog of 169 candidates. I estimate that the completeness of the survey is ~90% for spectral types of ≲M7 (≳0.06 M ⊙ ). The histogram of spectral types for the 32 Ori candidates exhibits a maximum at M5 (~0.15 M ⊙ ), resembling the distributions measured for other young clusters and associations in the solar neighborhood. The available UVW velocities indicate that the association is expanding, but they do not produce a well-defined kinematic age. Based on their sequences of low-mass stars in color–magnitude diagrams, the 32 Ori association and Upper Centaurus-Lupus/Lower Centaurus-Crux (UCL/LCC) are coeval to within ±1.2 Myr, and they are younger than the β Pic moving group by ~3 Myr, which agrees with results from previous analysis based on the second data release of Gaia. Finally, I have used mid-IR photometry from the Wide-field Infrared Survey Explorer to check for excess emission from circumstellar disks among the 32 Ori candidates. Disks are detected for 18 candidates, half of which are reported for the first time in this work. The fraction of candidates at ≤M6 that have full, transitional, or evolved disks is 10/149 = 0.07$^{+0.03}_{-0.02'}$ which is consistent with the value for UCL/LCC.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Ion Acceleration and the Development of a Power-law Energy Spectrum in Magnetic Reconnection

How charged particles are accelerated efficiently and form a power-law energy spectrum in magnetic reconnection is a problem that is not well understood. In a previous paper, it was shown that the electron Kelvin–Helmholtz instability (EKHI) in force-free magnetic reconnection generates fast-expanding vortices that can accelerate electrons in a few tens of ion gyroperiods (less than 1 ms in the solar corona) to form a power-law energy distribution. In this paper, we present a particle-in-cell (PIC) simulation study of ion acceleration in force-free magnetic reconnection in the presence of the EKHI-induced turbulence. We find that ions are not significantly accelerated by the EKHI-induced stochastic electric field until the magnetic vortices expand to sizes comparable to the ion gyroradius. The Alfvén waves generated by the EKHI couple with the magnetic vortices, leading to resonance between the ions inside the magnetic vortices and Alfvén waves and enhanced ion heating. The induced Alfvén wave resonance results in a broken power-law energy spectrum with a breakpoint at $\sim {m}_{i}{v}_{A}^{2}$, where v A is the Alfvén velocity. Here, we show that the process that forms the nonthermal tail is a second-order Fermi mechanism and the mean spectral index is α = (1 + 4 a 2 D/R) /2, where D is the spatial scale of the inductive electric field, R is that of vortices, and a = B g /B 0 , with ratio of guide field B g and asymptotic B 0 .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

On the Short-scale Spatial Variability of Electron Inflows in Electron-only Magnetic Reconnection in the Turbulent Magnetosheath Observed by MMS

We investigate the detailed properties of electron inflow in an electron-only reconnection event observed by the four Magnetospheric Multiscale (MMS) spacecraft in the Earth's turbulent magnetosheath downstream of the quasi-parallel bow shock. The lack of ion coupling was attributed to the small-scale sizes of the current sheets, and the observed bidirectional super-Alfvénic electron jets indicate that the MMS spacecraft crossed the reconnecting current sheet on both sides of an active X-line. Remarkably, the MMS spacecraft observed the presence of large asymmetries in the two electron inflows, with the inflows (normal to the current sheet) on the two sides of the reconnecting current layer differing by as much as a factor of four. Furthermore, even though the four MMS spacecraft were separated by less than seven electron skin depths, the degree of inflow asymmetry was significantly different at the different spacecraft. The asymmetry in the inflow speeds was larger with increasing distances downstream from the reconnection site, and the asymmetry was opposite on the two sides of the X-line. We compare the MMS observations with a 2D kinetic particle-in-cell (PIC) simulation and find that the asymmetry in the inflow speeds stems from in-plane currents generated via the combination of reconnection-mediated inflows and parallel flows along the magnetic separatrices in the presence of a large guide field.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Planar Collisionless Shock Simulations with the Semi-implicit Particle-in-cell Model FLEKS

This study investigates the applicability of the semi-implicit particle-in-cell code FLexible Exascale Kinetic Simulator (FLEKS) to heliospheric shock simulations. We examine one- and two-dimensional local planar shock simulations, initialized using MHD states with upstream conditions representative of plasmas in the hypersonic, β ∼ 1 regime, for both quasi-perpendicular and quasi-parallel configurations. The refined algorithm in FLEKS proves robust, enabling accurate shock simulations with a grid resolution on the order of the electron inertial length d e . Our simulations successfully capture key shock features, including shock structures (foot, ramp, overshoot, and undershoot), upstream and downstream waves (fast magnetosonic, whistler, Alfvén ion-cyclotron, and mirror modes), and non-Maxwellian particle distributions. Crucially, we find that at least two spatial dimensions are critical for accurately reproducing downstream-wave physics in quasi-perpendicular shocks and capturing the complex dynamics of quasi-parallel shocks, including surface rippling, shocklets, short, large-amplitude magnetic structures, magnetic reconnection, and jets. Furthermore, our parameter studies demonstrate the impact of mass ratio and grid resolution on shock physics. This work provides valuable guidance for selecting appropriate physical and numerical parameters for shock simulations using a semi-implicit PIC method, paving the way for incorporating kinetic shock processes into large-scale collisionless plasma simulations with the MHD-AEPIC model.

plasma astrophysics↗
DOE OSTI2025

On operation and control of CW magnetrons for superconducting accelerators

CW magnetrons, developed for industrial RF heaters, were suggested to feed RF cavities of superconducting accelerators due to higher efficiency and lower cost of RF power than provide traditionally used klystrons, IOTs or solid-state amplifiers. RF amplifiers driven by a master oscillator serve as coherent RF sources. CW magnetrons are regenerative RF auto generators with a huge regenerative gain. This causes regenerative instability with a notable noise when a magnetron operates as an auto generator i.e., with the anode voltage above the threshold of self-excitation. Traditionally, an injection locking by a small signal is used for phase stabilization of magnetrons. In this case CW magnetrons with the injection-locked (coherent) oscillations generate a notable level of noise. This may preclude use of CW magnetrons in this mode in the Superconducting RF (SRF) accelerators. This paper reviews PIC modeling and previously obtained experimental results for the operation and control of CW magnetrons, which led to the development of techniques most suitable for various SRF accelerators using forced RF oscillation of magnetrons, when the magnetron is launched by an injected forcing signal, and regenerative noise is suppressed.

Kazakevich, G. [MUONS Inc., Batavia] (ORCID:000000↗