Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle transport”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Review of the first charged-particle transport coefficient comparison workshop

We present the results of the first Charged-Particle Transport Coefficient Code Comparison Workshop, which was held in Albuquerque, NM October 4–6, 2016. In this first workshop, scientists from eight institutions and four countries gathered to compare calculations of transport coefficients including thermal and electrical conduction, electron–ion coupling, inter-ion diffusion, ion viscosity, and charged particle stopping powers. In this paper, we give general background on Coulomb coupling and computational expense, review where some transport coefficients appear in hydrodynamic equations, and present the submitted data. Large variations are found when either the relevant Coulomb coupling parameter is large or computational expense causes difficulties. Understanding the general accuracy and uncertainty associated with such transport coefficients is important for quantifying errors in hydrodynamic simulations of inertial confinement fusion and high-energy density experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Inference of main ion particle transport coefficients with experimentally constrained neutral ionization during edge localized mode recovery on DIII-D

Abstract The plasma and neutral density dynamics after an edge localized mode are investigated and utilized to infer the plasma transport coefficients for the density pedestal. The Lyman-Alpha Measurement Apparatus (LLAMA) diagnostic provides sub-millisecond profile measurements of the ionization and neutral density and shows significant poloidal asymmetries in both. Exploiting the absolute calibration of the LLAMA diagnostic allows quantitative comparison to the electron and main ion density profiles determined by charge-exchange recombination, Thomson scattering and interferometry. Separation of diffusion and convection contributions to the density pedestal transport are investigated through flux gradient methods and time-dependent forward modeling with Bayesian inference by adaptation of the Aurora transport code and IMPRAD framework to main ion particle transport. Both methods suggest time-dependent transport coefficients and are consistent with an inward particle pinch on the order of 1 m s −1 and diffusion coefficient of 0.05 m 2 s −1 in the steep density gradient region of the pedestal. While it is possible to recreate the experimentally observed phenomena with no pinch in the pedestal, low diffusion in the core and high outward convection in the near scrape-off layer are required without an inward pedestal pinch.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Particle transport in Markov stochastic mixtures with spatial gradients

Markov geometries are a prototype class of stochastic media that are widely used to model complex disordered systems. In a series of recent works, we have shown that spatially homogeneous (but possibly non-isotropic) Poisson tessellations can be conveniently used in order to generate three-dimensional realizations of Markov media and thus characterize the features of particle transport e.g. in fragmented fuel elements following severe nuclear accidents or Rayleigh-Taylor turbulent mixing for inertial confinement fusion. In this paper we show that Poisson tessellations can be extended to take into account the presence of spatial gradients, which commonly occur in real-world applications due to material stratification. For this purpose, we will provide an explicit construction of spatially-non-homogeneous Poisson tessellations, fully characterize their statistical properties, and finally illustrate the behaviour of particle transport in such random media. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Isotope physics of heat and particle transport with tritium in JET-ILW type-I ELMy H-mode plasmas

As part the DTE2 campaign in the JET tokamak, we conducted a parameter scan in T and D-T complementing existing pulses in H and D. For the different main ion masses, type-I ELMy H-modes at fixed plasma current and magnetic field can have the pedestal pressure varying by a factor of 4 and the total pressure changing from <?CDATA $\beta_{\mathrm{N}} = 1.0$?> β N = 1.0 to 3.0. We investigated the pedestal and core isotope mass dependencies using this extensive data set. The pedestal shows a strong mass dependence on the density, which influences the core due to the strong coupling between both plasma regions. To better understand the causes for the observed isotope mass dependence in the pedestal, we analysed the interplay between heat and particle transport and the edge localised mode (ELM) stability. For this purpose, we developed a dynamic ELM cycle model with basic transport assumptions and a realistic neutral penetration. The temporal evolution and resulting ELM frequency introduce an additional experimental constraint that conventional quasi-stationary transport analysis cannot provide. Our model shows that a mass dependence in the ELM stability or in the transport alone cannot explain the observations. One requires a mass dependence in the ELM stability as well as one in the particle sources. The core confinement time increases with pedestal pressure for all isotope masses due to profile stiffness and electromagnetic turbulence stabilisation. Interestingly, T and D-T plasmas show an improved core confinement time compared to H and D plasmas even for matched pedestal pressures. For T, this improvement is largely due to the unique pedestal composition of higher densities and lower temperatures than H and D. With a reduced gyroBohm factor at lower temperatures, more turbulent drive in the form of steeper gradients is required to transport the same amount of heat. This picture is supported by quasilinear flux-driven modelling using TGLF -SAT2 within Astra . With the experimental boundary condition TGLF -SAT2 predicts the core profiles well for gyroBohm heat fluxes <?CDATA ${\gt}15$?> > 15 , however, overestimates the heat and particle transport closer to the turbulent threshold.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Extension of the energetic particle transport kick model in TRANSP to multiple fast ion species

Abstract Alfvénic instabilities (AEs) are well known to cause enhanced transport of energetic particles (EPs) in fusion devices. Most studies until now have focused on characterizing and understanding AE stability in single-species plasmas heated by neutral beams (NB), where deuterium is typically used as both main plasma species and NB fuel. As the fusion community moves toward fusion reactors that target burning plasma conditions, such as ITER, the single-species picture breaks down. Burning plasmas, which will use a mix of deuterium and tritium (DT) as main fuel, also feature the presence of several supra-thermal fusion products such as alpha particles, protons, helium isotopes and high-energy tritium ions. This work presents the extension of the EP transport kick model implemented in the TRANSP time-dependent tokamak transport code to study the combined effect of multiple EP species on AE stability and, in turn, the response of different EP species to plasma instabilities in terms of their redistribution and losses. Further validation of the enhanced model is planned based on experimental results expected from the JET DT campaign scheduled for 2021, in preparation for ITER plasmas and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Density profiles in stellarators: an overview of particle transport, fuelling and profile shaping studies at TJ-II

We provide an overview of activities carried out at the TJ-II stellarator for improving our understanding of- and developing plasma physics models for particle density profiles in stellarators. Namely, we report on recent progress in turbulent particle transport simulation, validation of pellet deposition models, density profile shaping for performance control and new experimental techniques for edge turbulence and plasma-neutral interaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Lagrangian dynamics of particle transport in oral and nasal breathing

Here, we present a large-eddy simulation (LES) of saliva particle transport during normal human breathing through the nose and mouth. The flow of the air–saliva mixture is modeled using an Eulerian LES that is coupled with a Lagrangian particle tracking module to obtain trajectories of saliva particles in a room with stagnant air conditions. The coupled Eulerian–Lagrangian simulation yields novel insights into the intricate dynamics of Lagrangian coherent structures (LCS) and fundamental material lines that emerge from the saliva particles' trajectories during several breathing cycles. Further, we systematically compare the quantitative LCS diagnostics of mouth breathing with those of mouth and nose normal breathing. Analyzing the simulation results of human breathing from the mouth and nose, we show that, soon after the first breathing cycle, saliva particles form a series of roll-up vortex rings that propagate forward. The forward propagation of these vortex rings leads to the formation of an asymmetrical primary forefront vortex. The individual vortex rings continuously propagate forward, merging with the forefront vortex, and ascending along the limb of the leading vortex.

42 ENGINEERING↗
DOE OSTI2025

A charged particle transport approximation for thick and thin plasmas

In this work, we will construct a simple method to be utilized in the interpretation of the properties of a thin plasma from a study of the transport of charged particles through it. We do so by first demonstrating how charged particle fluxes encode local plasma information in their spectra for sufficiently thick plasmas and then propose how to extend that analysis to the thin plasma limit where the extensive geometry cannot be neglected through arguments from scale separation. We provide a numerical treatment of the transport problem and demonstrate that it is in good agreement with the approximation we present. Finally, we utilize both the numerical solution and our novel approximation to study the impact of the extensive scale of the plasma on the yield and flux normalized high-energy neutron spectra resulting from the upscattering of charged fuel ions. Using this analysis, we show that (after controlling for fusion yield) for the same uniform densities and temperatures, larger plasmas have higher magnitude but softer reaction-in-flight (RIF) neutron spectra relative to smaller plasmas. This is because larger plasmas retain more knocked-on suprathermal ions within their bulk and can downscatter them to lower average energies, while smaller plasmas allow a larger fraction of high energy particles to “range out” of the system prior to substantial downscattering or inducing an RIF reaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Machine Learning Approach to Modeling of Neutral Particles Transport in Plasma

A propagator‐based approach is investigated for Monte‐Carlo (MC) modeling of neutral particle transport in fusion boundary plasmas. The propagator is based on a Green's function for the neutral kinetic equation, which depends on the plasma profiles. A neural network (NN)‐based model for the propagator provides a fast and accurate solution for the neutral distribution function in plasma. Preliminary results from a small 1D test problem look encouraging. The proposed approach, a propagator‐based NN model for neutral transport in plasma, has potential for generalization to higher dimensions and efficient coupling with plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Comparing particle transport in JET and DIII-D plasmas: gyrokinetic and gyrofluid modelling

The paper models the impact of fueling versus particle transport in the DIII-D and JET tokamaks. Experimentally, it was found that in JET an increase in the density peaking at lower collisionality could be attributed to an increase in core fueling, however on DIII-D, the opposite was observed. In DIII-D, the increase in density peaking was the result of an increase in inward transport, not an increase in core fueling. The paper shows that the difference in temperature gradients and thus turbulence dynamics between the two experiments can explain the experimentally observed discrepancy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Transition in particle transport under resonant magnetic perturbations in a tokamak

Abstract Nonlinear 3D MHD simulations and validations reveal that the hybrid particle-MHD transport is a key process for driving the pump-out in the presence of Resonant Magnetic Perturbations (RMPs) in the KSTAR tokamak. Particle transport and the resulting density pump-out by RMPs are shown to be composed of not only the classical flow convection near magnetic islands due to polarization but also the neoclassical ion diffusion across perturbed magnetic surfaces. The latter is known as the Neoclassical Toroidal Viscosity (NTV) and is integrated into nonlinear MHD simulations here for the first time, revealing that the two-stage pump-outs observed in KSTAR experiments are reproduced only with such integrated nonlinear MHD and transport evolution. Near-resonant responses, which have received less attention than the resonant response, play distinct roles in the pump-out along with the island formation. In addition, this modeling is used to investigate the pump-outs in double-null-like plasmas and numerically capture the effect of the double-null shape on the pump-outs, which may explain the difficulty of Edge Localized Mode (ELM) suppression access in double-like plasmas. This reveals new aspects of the impact toroidal geometry and mode coupling have on 3D physics and reveals the importance of near-resonant components in suppressing ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

GeantV: Results from the Prototype of Concurrent Vector Particle Transport Simulation in HEP

Full detector simulation was among the largest CPU consumers in all CERN experiment software stacks for the first two runs of the Large Hadron Collider. In the early 2010s, it was projected that simulation demands would scale linearly with increasing luminosity, with only partial compensation from increasing computing resources. The extension of fast simulation approaches to cover more use cases that represent a larger fraction of the simulation budget is only part of the solution, because of intrinsic precision limitations. The remainder corresponds to speeding up the simulation software by several factors, which is not achievable by just applying simple optimizations to the current code base. In this context, the GeantV R&D project was launched, aiming to redesign the legacy particle transport code in order to benefit from features of fine-grained parallelism, including vectorization and increased locality of both instruction and data. This paper provides an extensive presentation of the results and achievements of this R&D project, as well as the conclusions and lessons learned from the beta version prototype.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Charged Particle Transport Libraries (Final Report FY2020 Level-2 Milestone #7127)

This document is the final report on the status and developments of the Charged-Particle Transport (CPT) Level-2 (“CP2020”) ASC/PEM Milestone work. This section gives an overview of the CP2020 work partitioned into two subsections: 1) a summary of the code tools used for generating CP (and other) data in Section 1.2; and 2) and library data and codes access information for end users in Section 1.3. Section 2 gives an overview of the status of previous work on charged-particle (CP) nuclear data, particularly focused on the work of the FY2011 CP Level-2 Milestone. New work for CP2020 is detailed in the following Section 3. Section 4 discusses testing and validation of the evaluated and processed data for the CP2020 library. Conclusions based on the experience of this Level-2 work and the outlook for future development of CPT-relevant data are given in the final Section 5. Appendices, detailed throughout the report, are provided for further information. This report, related documentation, codes, and the CPT nuclear data libraries are available via online repository.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Neutron (and other Particle) Transport at LANL: An Overview [Presentation]

For decades, Los Alamos National Laboratory has been at the forefront of neutron transport methods research and code development. One such code is PARTISN, the LANL parallel time-dependent discrete ordinate neutron transport code. In this presentation, we describe the various research efforts currently underway by the PARTISN and other code teams. Some examples of current research are a block automated mesh refinement scheme, the application of tensor trains to the discretized neutron transport equation, and GPU code porting. The block automated mesh refinement scheme uses cross section information to refine and coarsen the solution mesh to improve time to solution and reduce memory. The tensor train approach expresses discretized transport operators as tensor products of vectors and matrices to compress the size of linear systems being solved by transport codes. Rather than relying on matrix-free methods such as the transport sweep, we have access to an operator that can be inverted, reshaped, or manipulated algebraically. Finally, we describe how PARTISN is used, what problems we are looking to solve, and what the future holds for neutron transport at LANL. In addition to this, we briefly describe the various research efforts in other particle transport teams using both deterministic and Monte Carlo methods. In the presentation, we list possible opportunities for collaboration between the laboratory and faculty and students.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Implementation of the D1S Methodology for Shutdown Dose Rate Calculations in the OpenMC Monte Carlo Particle Transport Code

We present an implementation of the direct one-step (D1S) methodology for shutdown dose rate (SDR) calculations in the OpenMC Monte Carlo particle transport code. In addition to being the first fully open-source D1S implementation, it is also the first to require no ad hoc source code or nuclear data library modifications. The code can seamlessly switch between production of prompt and decay photons based on a user input parameter, and the decay data needed for decay photon generation are made available through a depletion chain file, which is already used for OpenMC’s built-in depletion/activation solver. A set of Python functions significantly eases the burden of computing and applying time correction factors needed to properly account for the time dependence of radionuclide activity. To assess the accuracy of the D1S implementation, SDR calculations have been carried out for three problems: a prism of iron irradiated by 14-MeV neutrons, the ITER port plug computational benchmark, and the Frascati Neutron Generator (FNG) ITER dose rate benchmark problem from the Shielding INtegral Benchmark Archive and Database (SINBAD). For each of these problems, comparisons were made to calculations using the rigorous two-step (R2S) method. The results on the iron prism problem illustrate how the D1S method achieves superior spatial resolution compared to the R2S method without the need for spatial discretization of the activation regions. The D1S and R2S results for the ITER port plug benchmark agree well with previously reported results in the literature. While the D1S results are 10% to 15% lower than the R2S results, this may be due to stochastic uncertainty and/or spatial discretization in the R2S calculations. On the FNG dose rate benchmark problem, the D1S method produces dose rate estimates that are within 4% of the dose rates predicted using a cell-based R2S workflow. The D1S estimates of the SDR are also in reasonable agreement with the experimental measurements and show the same basic trends that have been observed in previous works. A qualitative analysis of the execution time and uncertainty for the R2S and D1S workflows suggests that the D1S method would attain a higher figure of merit.

D1S method↗
DOE OSTI2022

Validation of energetic particle transport models for time-dependent integrated simulations of burning plasmas (Final ST Report)

The major goal of this project is to deliver a suite of validated energetic particle (EP) modeling tools for time-dependent integrated tokamak simulations (with the code TRANSP) that can be used for modeling plasma discharges featuring a substantial EP content, such as burning plasmas. During FY20, issues were resolved with an energetic particle transport model called TGLF-EP, with improvements made to the TGLF-EP workflow, generation of input files, and interpretation of runs on the NERSC supercomputing cluster. The recipient also gave input on simulations with the code TRANSP to reproduce the reference JET discharge, and assessed the accuracy and uncertainties of the simulated B field, current evolution, thermal plasma profiles, stored energy, etc. Input was also given on interpretation of fast-ion instability calculations and recommendations for modifications to JET afterglow discharge for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Internal erosion, particle transport, and channelization driven by fluid flow

The investigations summarized in this final technical report were accomplished under Department of Energy, Basic Energy Sciences research grant number DE-SC0010274 titled Internal erosion, particle transport, and channelization driven by fluid flow, $196,683 from 07/15/13 to 12/31/17, and $149,999 from 01/01/18 to 12/31/21. The work focused on developing physical models of fluid driven evolution of sediment beds as a result of dissolution and erosion of embedded particles that lead to channelization and fracture networks in porous rock. The grant work has resulted in seven peer reviewed publications, with another one currently under review, all of which are listed at the end of the report. These publications can be freely accessed through a technical library or directly from the publisher for a nominal fee. The main results obtained during the last grant cycle are described in brief in the following.

58 GEOSCIENCES↗
DOE OSTI2022

Celeritas: GPU-accelerated particle transport for detector simulation in High Energy Physics experiments

Within the next decade, experimental High Energy Physics (HEP) will enter a new era of scientific discovery through a set of targeted programs recommended by the Particle Physics Project Prioritization Panel (P5), including the upcoming High Luminosity Large Hadron Collider (LHC) HL-LHC upgrade and the Deep Underground Neutrino Experiment (DUNE). These efforts in the Energy and Intensity Frontiers will require an unprecedented amount of computational capacity on many fronts including Monte Carlo (MC) detector simulation. In order to alleviate this impending computational bottleneck, the Celeritas MC particle transport code is designed to leverage the new generation of heterogeneous computer architectures, including the exascale computing power of U.S. Department of Energy (DOE) Leadership Computing Facilities (LCFs), to model targeted HEP detector problems at the full fidelity of Geant4. This paper presents the planned roadmap for Celeritas, including its proposed code architecture, physics capabilities, and strategies for integrating it with existing and future experimental HEP computing workflows.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗