Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupled 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

FENIX: Towards a Fully Integrated Multiphysics Framework for Plasma Facing Component Modeling

Computational tools have a crucial role to play in accelerating the deployment of fusion as a clean, reliable, abundant, and sustainable energy source. Multiphysics, high-fidelity simulation capabilities can help model, study, and predict intricate interactions between materials performance, plasma exposure, neutron irradiation, and engineering processes. As such, they can assist in the resolution of scientific and engineering challenges underpinning design, construction, and commission of fusion power plants. To address these needs, ongoing efforts are leveraging the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework and delivering new computational tools for the fusion community. These tools inherit crucial attributes from MOOSE. They are open-source, modular, integrated with nuclear industry-standard software quality assurance processes, and enable multiphysics, multi-fidelity, fully integrated, zero- to three-dimensional, and massively parallel simulations. After a short overview of these capabilities, we will present the development of Fusion ENergy Integrated multiphys-X (FENIX), a MOOSE-based application designed to enable plasma facing component design and performance evaluation. Throughout their lifetime, plasma facing components are exposed to extreme thermal loads, repeated thermal shocks, and irradiation by plasma ions, neutral particles, and high-energy neutrons. Consequently, designing a plasma facing component with acceptable lifetime degradation is extremely challenging. FENIX aims to model the multiphysics environment in which plasma facing components evolve to accelerate their design studies. To that end, FENIX couples existing MOOSE capabilities such as heat transfer, thermomechanics, and thermal hydraulics, with tritium transport via the MOOSE-based Tritium Migration Analysis Program, Version 8 (TMAP8), with neutronics via the MOOSE-based high-fidelity neutron-photon transport and fluid dynamics code Cardinal, and finally with Particle-in-Cell plasma simulation capabilities being developed in this project. In this study, we present the current FENIX capabilities and preliminary results of its application to model the Tritium Plasma Experiment set up at Idaho National Laboratory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Fully generalized, turbulent trace impurity transport with Gkeyll and Flan in the DIII-D far-SOL

The Monte Carlo trace impurity turbulent transport code Flan is introduced for the first time. Flan follows impurities in a turbulent background plasma from Gkeyll using the Lorentz force to resolve the full particle gyro orbit. Collisions are handled using the Nanbu collision algorithm (Nanbu 1997 Phys. Rev. E 55 4642–52), and ionization/recombination is handled via ADAS coupling. The far-SOL of a generic DIII-D L-mode is simulated with and without the collision model to show how collisions affect radial tungsten transport. Anomalous diffusion coefficient (D r ) and pinch velocity (v p ) profiles are extracted from fits to the results. With collisions, D r and v p are between 0–1.0 m 2 s −1 and −100–100 m s −1 , respectively. Without collisions, D r and v p are between 0–0.3 m 2 s −1 and −50–50 m s −1 , respectively. Exponential fits to the radial W density profiles and experimental data from W deposition along a collector probe are in good agreement, demonstrating Flan as a useful interpretive modeling tool. Additional simulations show that impurity transport away from the wall increases with atomic number, though it is not clear why. Flan has the potential to better interpret existing data and improve reactor scale predictions of core contamination because the underlying physics model is very general and does not rely on arbitrary user-defined transport coefficients.

DIII-D↗
DOE OSTI2020

Depletion capabilities in the OpenMC Monte Carlo particle transport code

A depletion solver has been implemented in OpenMC and is described herein. The depletion solver is implemented in Python and interfaces with OpenMC’s transport solver through a C++ application programming interface, which enables an in-memory transport-depletion coupling. Multiple integration methods for advancing in time have been implemented and exhibit tradeoffs in cost, accuracy, and memory use. For all time integration methods, evaluation of the matrix exponential is performed by using the incomplete partial fraction form of the Chebyshev rational approximation method. Simulations of a pressurized water reactor (PWR) pincell and a sodium-cooled fast reactor (SFR) assembly were carried out with OpenMC and Serpent. For both problems, the use of a high-fidelity depletion chain results in predictions of k eff that agree within 20–30 pcm between OpenMC and Serpent. Predicted actinide concentrations were found to agree to a fraction of a percent, and most fission product concentrations were found to agree within 1%. Here, the few cases where larger differences were observed can be attributed either to differences in how the energy dependence of fission product yields is handled or deficiencies in the nuclear data used. OpenMC simulations of the PWR and SFR problems using a simplified 228-nuclide depletion chain demonstrate that it achieves accuracy close to that of the full, high-fidelity depletion chain with respect to the studied responses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Photon detector response function methodology using MCNP and shift hybrid radiation transport code for wide-area contamination assay applications

Here, radiation transport modeling using the Monte Carlo N-Particle (MCNP) radiation transport code and Monte Carlo code, Shift, were employed to model detector responses for a variety of wide-area photon contamination scenarios. In this study, 2" × 2" and 3" × 3" cylindrical NaI(Tl) scintillation detector configurations at source detector-distances of 0.5 cm, 1 cm, 2.54 cm, 10 cm, and 30 cm were modeled. Media of soil, concrete, and steel were evaluated for contamination depths ranging from surface to a depth of an infinite thickness in each medium for photon energies ranging from 20 keV to 3 MeV, which correspond to the energies that current detectors can discern. Monoenergetic photon surface contamination detector responses for each of the media, source–detector distances, and detectors were estimated using MCNP v6.2. Shift was harnessed for improved variance reduction of particle transport in highly attenuating media to obtain average cell fluxes in the two MCNP NaI(Tl) scintillation detector configurations. Average cell flux values in Shift were coupled with detector responses from MCNP to convert average cell flux in a void to energy distribution of pulses in the NaI(Tl) scintillation detector crystal of interest. An optimized detector response function methodology was developed by coupling the MCNP radiation transport method with the Consistent Adjoint Driven Importance Sampling (CADIS) hybrid radiation transport method built into Shift to significantly decrease the runtime of thousands of MCNP pulse height simulations. The methodology may be utilized to quickly and accurately facilitate the assessment of a broad range of wide-area environmental contamination assay and decommissioning cleanup applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Quasiparticle Couplings in Transport of Heat, Charge, and Spin for Novel Energy Materials. Final report

Understanding the microscopic processes involved in the transport and conversion of energy from the atomic-scale to the meso-scale is critical for the development of next-generation materials for energy conversion and information processing. The goal of this program was to elucidate the couplings between microscopic degrees-of-freedom of atomic vibrations, spins, and electrons, which underpin heat, charge, and spin transport. This project systematically investigated these quasi-particle interactions, using state-of-the-art neutron scattering techniques, combined with first-principles computer simulations. This research program uncovered new information and understanding, which helps to rationalize materials properties based on the interplay of their atomistic degrees of freedom.

36 MATERIALS SCIENCE↗
NASA NTRS2008

Sensor-Web Operations Explorer

Understanding the atmospheric state and its impact on air quality requires observations of trace gases, aerosols, clouds, and physical parameters across temporal and spatial scales that range from minutes to days and from meters to more than 10,000 kilometers. Observations include continuous local monitoring for particle formation; field campaigns for emissions, local transport, and chemistry; and periodic global measurements for continental transport and chemistry. Understanding includes global data assimilation framework capable of hierarchical coupling, dynamic integration of chemical data and atmospheric models, and feedback loops between models and observations. The objective of the sensor-web system is to observe trace gases, aerosols, clouds, and physical parameters, an integrated observation infrastructure composed of space-borne, air-borne, and in-situ sensors will be simulated based on their measurement physics properties. The objective of the sensor-web operation is to optimally plan for heterogeneous multiple sensors, the sampling strategies will be explored and science impact will be analyzed based on comprehensive modeling of atmospheric phenomena including convection, transport, and chemical process. Topics include system architecture, software architecture, hardware architecture, process flow, technology infusion, challenges, and future direction.

observation↗
DOE OSTI2025

Integrated modeling of alpha particle losses in JET DT plasmas

Alpha particle confinement is crucial for sustaining burning plasmas and designing future reactor concepts. Along with classical/prompt losses, various magnetohydrodynamic instabilities can lead to wave-particle interactions which can transport alpha particles outward from the plasma. This can result in a reduction in plasma heating/performance, and, at worst, damage in-vessel components. Joint European Torus’s recent deuterium–tritium campaigns in 2021–2023 have produced numerous alpha particle loss measurements with its scintillator probe and Faraday cup array fast ion loss detectors as discussed in Bonofiglo et al (2024 Nucl. Fusion 64 096038). This paper will report on integrated energetic particle transport modeling in support of those measurements. The modeling is accomplished with the TRANSP and ORBIT-kick codes with the use of recently developed reduced models which calculate mode structure, amplitude, and the evolving dynamics. When possible, constraints and comparisons to experiment are conducted. Case studies are performed on a variety of magnetohydrodynamic activity, including: fishbones, tearing modes (TMs), and sawtooth crashes. Additionally, a special case of an alpha-driven toroidal Alfvén eigenmode is briefly discussed, where modeling showed marginally weak alpha losses and was unable to support experimental observations. Coupled effects between a TM and toroidal field ripple are presented and were unable to replicate the observations in lost particle pitch but did duplicate the localized flattening of the measured neutron profile. Additional modeling results compare the magnitude of losses and energy/velocity-space sensitivities against experimental observations/measurements for each scenario. This work corroborates numeric alpha transport modeling while also identifying model deficiencies. While this report details alpha transport, it also presents open issues for discussion in assessing the validity of our numerical models towards burning plasmas.

Bonofiglo, Phillip J. [Princeton Plasma Physics La↗
NASA NTRS2008

An Analytical Model for the Prediction of a Micro-Dosimeter Response Function

A rapid analytical procedure for the prediction of a micro-dosimeter response function in low Earth orbit (LEO), correlated with the Space Transportation System (STS, shuttle) Tissue Equivalent Proportional Counter (TEPC) measurements is presented. The analytical model takes into consideration the energy loss straggling and chord length distribution of the detector, and is capable of predicting energy deposition fluctuations in a cylindrical micro-volume of arbitrary aspect ratio (height/diameter) by incoming ions through both direct and indirect (ray) events. At any designated (ray traced) target point within the vehicle, the model accepts the differential flux spectrum of Galactic Cosmic Rays (GCR) and/or trapped protons at LEO as input. On a desktop PC, the response function of TEPC for each ion in the GCR/trapped field is computed at the average rate of 30 seconds/ion. The ionizing radiation environment at LEO is represented by O'Neill fs GCR model (2004), covering charged particles in the 1 less than or equal to Z less than or equal to 28. O'Neill's free space GCR model is coupled with the Langley Research Center (LaRC) angular dependent geomagnetic cutoff model to compute the transmission coefficient in LEO. The trapped proton environment is represented by a LaRC developed time dependent procedure which couples the AP8MIN/AP8MAX, Deep River Neutron Monitor (DRNM) and F10.7 solar radio frequency measurements. The albedo neutron environment is represented by the extrapolation of the Atmospheric Ionizing Radiation (AIR) measurements. The charged particle transport calculations correlated with STS 51 and 114 flights are accomplished by using the most recent version (2005) of the LaRC deterministic High charge (Z) and Energy TRaNsport (HZETRN) code. We present the correlations between the TEPC model predictions (response function) and TEPC measured differential/integral spectra in the lineal energy (y) domain for both GCR and trapped protons, with the conclusion that the model correctly accounts for the increase in flux at low y values where energetic ions are the primary contributor. We further discuss that, even with the incorporation of angular dependency in the cutoffs, comparison of the GCR differential/integral flux between STS 51 and 114 TEPC measured data and current calculations indicates that there still exists an underestimation by the simulations at low to mid range y values. This underestimation is partly related the exclusion of the secondary pion particle production from the current version of HZETRN.

Badavi, Francis F.↗
NASA NTRS2010

Faster and More Accurate Transport Procedures for HZETRN

Several aspects of code verification are examined for HZETRN. First, a detailed derivation of the numerical marching algorithms is given. Next, a new numerical method for light particle transport is presented, and improvements to the heavy ion transport algorithm are discussed. A summary of various coding errors is also given, and the impact of these errors on exposure quantities is shown. Finally, a coupled convergence study is conducted. From this study, it is shown that past efforts in quantifying the numerical error in HZETRN were hindered by single precision calculations and computational resources. It is also determined that almost all of the discretization error in HZETRN is caused by charged target fragments below 50 AMeV. Total discretization errors are given for the old and new algorithms, and the improved accuracy of the new numerical methods is demonstrated. Run time comparisons are given for three applications in which HZETRN is commonly used. The new algorithms are found to be almost 100 times faster for solar particle event simulations and almost 10 times faster for galactic cosmic ray simulations.

Slaba, Tony C.↗
DOE OSTI2021

Modeling Coupled Chemo-Mechanical Behavior of Randomly Oriented NMC811 Polycrystalline Li-Ion Battery Cathodes

This paper develops a three-dimensional, transient, chemo-mechanical model that predicts the performance of single secondary particle Li-ion battery cathodes. The secondary particles are composed of numerous (approximately 60) randomly oriented single-crystal primary particles. The model incorporates concentration-dependent and anisotropic material properties. As much as possible, electrochemical, transport, and structural properties for crystalline NMC811 (Li x Ni 0.8 Mn 0.1 Co 0.1 O 2 ) are taken from prior publications. Weak Van der Waals bonding between primary particles is modeled empirically using a spring analogy, which enables local primary-particle separations (disintegration) and subsequent reattachments. The model fully couples Li diffusion and the mechanical response. Results include predictions of local Li-concentrations and stresses. High stresses are found near grain boundaries, especially when the lattice orientations are greatly misaligned. Particle separations are characterized in terms of a damage parameter. The model is used to predict the effects of design and operating conditions, including charge/discharge rates, cycling scenarios, and particle sizes.

25 ENERGY STORAGE↗
DOE OSTI2021

Sustained edge-localized-modes suppression and radiative divertor with an impurity-driven instability in tokamak plasmas

We report simultaneous control of the large edge localized modes (ELMs) and divertor heat fluxes in a metal wall environment is a critical issue for steady-state operation of a tokamak fusion reactors. Here we report a sustained ELM suppression scenario achieved in the EAST tokamak compatible with radiative divertor using different seeding impurity species over a wide range of conditions. A low-n mode appears, as manifested by the oscillations of a radiation front near the X-point. This mode appears to drive strong particle transport and tungsten exhaust, which is essential to the maintenance of the ELM-stable state. We have developed a model to explain the mode excitation, by coupling the impurity radiative condensation instability to drift waves, which could explain some characteristics of the low-n mode well. The low-n mode may offer a new ELM-stable scenario compatible with radiative divertor for future fusion reactors.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1982

Method of making formulated plastic separators for soluble electrode cells

A method making a membrane comprised of a hydrochloric acid-insoluble sheet of a mixture of a rubber and a powdered ion transport material is disclosed. The sheet can be present as a coating upon a flexible and porous substrate. These membranes can be used in oxidation-reduction electrical accumulator cells wherein the reduction of one member of a couple is accompained by the oxidation of the other member of the couple on the other side of the cell and this must be accompained by a change in chloride ion concentration in both sides. The method comprises preparing a mixture of fine rubber particles, a solvent for the rubber and a powdered ion transport material. The mixture is formed into a sheet and dried to produce a microporous sheet. The ion transport material includes particles ranging from about 0.01 to 10 microns in size and comprises from 20 to 50 volume percent of the microporous sheet.

Sheibley, D. W.↗
DOE OSTI2024

Cross-scale interaction between microturbulence and meso-scale reversed shear Alfvén eigenmodes in DIII-D plasmas

Abstract This paper reports global nonlinear gyrokinetic simulations that couple meso-scale reversed shear Alfvén eigenmodes (RSAEs) driven by energetic particles (EPs) and ion temperature gradient (ITG) microturbulence driven by thermal plasma, using equilibrium and profiles from DIII-D discharge #159243. In simulations focusing only on the ITG, electrostatic ITG drives a huge thermal ion heat transport, which is reduced by a factor of 10 to a level close to the experimental value in electromagnetic simulation due to finite β effect. In the simulations coupling the RSAE and ITG, ITG can scatter the resonant EP nonlinearly trapped by the RSAE and damp the zonal flows generated by the RSAE. The regulation of the RSAE by the ITG greatly reduces the initial saturation amplitude of the RSAE but increases the RSAE amplitude and associated EP transport to experimental levels in the quasi-steady state. The RSAE effects on the ITG, specifically the stronger zonal flows generated by the RSAE and the RSAE frequency modulation of the ITG-induced thermal ion heat transport, in turn, leads to a reduction of the thermal ion heat transport by more than a factor of 2 . For a stronger background ITG, the regulation of the RSAE by the ITG is stronger, while the RSAE effects on the ITG are weaker. This work highlights the importance of cross-scale coupling in the dynamics of the AE turbulence and EP transport.

Physics↗
DOE OSTI2025

Verification of a Fluid-Based Plasma-Edge Model Within the Multiphysics Object-Oriented Simulation Environment (MOOSE) Framework

As the goal of achieving fusion power on the grid comes closer to fruition, fully coupled multiphysics models of fusion devices will be crucial. These models must incorporate the interconnected phenomena of these devices, including plasma physics, neutronics, first wall interactions, and tritium transport. Currently, there are two main approaches to developing these platforms: (1) loosely coupled, where one couples existing codes and solvers together through input and output parameters and data, and (2) tightly coupled, where one develops the necessary models within a singular, integrated framework. This work focuses on the latter approach for magnetically confined fusion devices by developing a fluid-based plasma-edge model within the Multiphysics Object Oriented Simulation Environment (MOOSE) Framework. This effort is coordinated with other efforts to develop, test, demonstrate, and deploy fusion relevant multiphysics capabilities including electromagnetics, particle-in-cell plasma, tritium transport, and fusion blanket design. This new model is an expansion of the MOOSE-based plasma application, Zapdos, which was originally formulated to model low-temperature, non-magnetized plasma processes. Verification studies have been conducted using newly developed magnetic plasma capabilities. These involved convergence analyses utilizing the method of manufactured solutions to verify new operators and case studies. A modular approach was taken here to demonstrate increasingly complicated simulation scenarios, which included a singular fluid with uniform magnetic field case, a singular fluid with spatially varying magnetic field case, and a coupled multifluid case.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Utilization of ACE nuclear data file toolkit ACEtk to calculate relative sensitivity coefficients of point-kinetics parameters

Sensitivity and uncertainty methods are quintessential for nuclear criticality safety and experiment design. This type of analysis relies on calculations of sensitivity coefficients; sensitivity coefficients of the effective neutron multiplication factor with respect to some nuclear data are predominantly calculated and used. As a part of the Laboratory Directed Research & Development project EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) at Los Alamos National Laboratory, sensitivity coefficients of many radiation detector measurement responses with respect to nuclear data were investigated. Specifically, this paper outlines a method to calculate point-kinetics parameters relative sensitivity coefficients with respect to nuclear data. Point-kinetics parameters such as the prompt neutron decay constant, effective delayed neutron fraction, and neutron generation time are especially important to experimenters and reactor operators designing systems with dynamic neutron populations. This method couples capabilities of the ACE (A Compact ENDF) nuclear data file toolkit, ACEtk, with the ability to load cross sections into the radiation transport code Monte Carlo N-Particle (MCNP). In conclusion, key aspects of optimizing this method for a particular application and sensitivity profiles of the Jezebel criticality experiment are examined and discussed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2023

On the Impact of Mechanics on Electrochemistry of Lithium-Ion Battery Anodes

Abstract Models exploring electrochemistry-mechanics coupling in liquid electrolyte lithium-ion battery anodes have traditionally incorporated stress impact on thermodynamics, bulk diffusive transport, and fracture, while stress-kinetics coupling is more explored in the context of all solid-state batteries. Here, we showcase the existence of strong link between active particle surface pressure and reaction kinetics affecting performance even in liquid electrolyte systems. Traction-free and immobile particle surface mechanical boundary conditions are used to delineate the varying pressure magnitudes in graphite host during cycling. Both tensile and compressive stresses are generated in traction-free case, while a fixed surface subjects the entire particle to a compression state. Pressure magnitudes are nearly two to three orders of magnitude higher for the latter resulting in significant depression of open circuit potential and improvement of exchange current densities compared to stress-free state. The results demonstrate the need for incorporating stress-kinetics linkage in models and provide a rationale for putting battery electrodes under compression to improve kinetics.

25 ENERGY STORAGE↗
DOE OSTI2023

Automated bidirectional coupling of multiscale models of aerosol dosimetry: Validation with subject-specific deposition data

Assessing the toxicity of airborne particulate matter or the efficacy of inhaled drug depends upon accurate estimates of deposited fraction of inhaled materials. In silico approaches can provide important insights into site- or airway-specific deposition of inhaled aerosols in the respiratory system. In this study, we improved on our recently developed 3D/1D model that simulate aerosol transport and deposition in the whole lung over multiple breath cycles. A subject-specific multiscale lung model of a healthy male subject using computational fluid-particle dynamics (CFPD) in a 3D model of the oral cavity through the large bronchial airways entering each lobe was bidirectionally coupled with a recently improved Multiple Path Particle Dosimetry (MPPD) model to predict aerosol deposition over the entire respiratory tract over multiple breaths for four conditions matching experimental aerosol exposures in the same subject from which the model was developed. These include two particle sizes (1 and 2.9 μm) and two subject-specific breathing rates of ~300 ml/s (slow breathing) and ~750 ml/s (fast breathing) at a target tidal volume of 1 L. In silico predictions of retained fraction were 0.31 and 0.29 for 1 μm and 0.66 and 0.62 for 2.9 μm during slow and fast breathing, respectively, and compared well with experimental data (1 μm: 0.31 ± 0.01 (slow) and 0.27 ± 0.01 (fast), 2.9 μm: 0.63 ± 0.03 (slow) and 0.68 ± 0.02 (fast)). Furthermore, these results provide a great deal of confidence in the validity and reliability of our approach.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Clustering of inertial particles in turbulent flow through a porous unit cell

We report direct numerical simulation is used to investigate effects of turbulent flow in the confined geometry of a face-centred cubic porous unit cell on the transport, clustering and deposition of fine particles at different Stokes numbers (St = 0.01, 0.1, 0.5, 1, 2) and at a pore Reynolds number of 500. Particles are advanced using one-way coupling and the collision of particles with pore walls is modelled as perfectly elastic with specular reflection. Tools for studying inertial particle dynamics and clustering developed for homogeneous flows are adapted to take into account the embedded, curved geometry of the pore walls. The pattern and dynamics of clustering are investigated using the volume change of Voronoi tesselation in time to analyse the divergence and convergence of the particles. Similar to the case of homogeneous, isotropic turbulence, the cluster formation is present at large volumes, while cluster destruction is prominent at small volumes and these effects are amplified with the Stokes number. However, unlike homogeneous, isotropic turbulence, the formation of a large number of very small volumes was observed at all Stokes numbers and attributed to the collision of particles with the pore wall. Multiscale wavelet analysis of the particle number density indicates that the peak of the energy density spectrum, representative of enhanced particle clustering, shifts towards larger scales with an increase in the Stokes number. Scale-dependent skewness and flatness quantify the intermittent void and cluster distribution, with cluster formation observed at small scales for all Stokes numbers, and void regions at large scales for large Stokes numbers.

42 ENGINEERING↗