Search NASA⌕ Search

SEARCH · Search NASA

Results for “charged 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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Description of a Generalized Analytical Model for the Micro-dosimeter Response

An analytical prediction capability for space radiation in Low Earth Orbit (LEO), correlated with the Space Transportation System (STS) Shuttle Tissue Equivalent Proportional Counter (TEPC) measurements, is presented. The model takes into consideration the energy loss straggling and chord length distribution of the TEPC detector, and is capable of predicting energy deposition fluctuations in a micro-volume by incoming ions through both direct and indirect ionic events. The charged particle transport calculations correlated with STS 56, 51, 110 and 114 flights are accomplished by utilizing the most recent version (2005) of the Langley Research Center (LaRC) deterministic ionized particle transport code High charge (Z) and Energy TRaNsport WZETRN), which has been extensively validated with laboratory beam measurements and available space flight data. The agreement between the TEPC model prediction (response function) and the TEPC measured differential and integral spectra in lineal energy (y) domain is promising.

Badavi, Francis F.↗
NASA NTRS1997

Angular neutron transport investigation in the HZETRN free-space ion and nucleon transport and shielding computer program

Extension of the high charge and energy (HZE) transport computer program HZETRN for angular transport of neutrons is considered. For this paper, only light ion transport, He4 and lighter, will be analyzed using a pure solar proton source. The angular transport calculator is the ANISN/PC program which is being controlled by the HZETRN program. The neutron flux values are compared for straight-ahead transport and angular transport in one dimension. The shield material is aluminum and the target material is water. The thickness of these materials is varied; however, only the largest model calculated is reported which is 50 gm/sq cm of aluminum and 100 gm/sq cm of water. The flux from the ANISN/PC calculation is about two orders of magnitude lower than the flux from HZETRN for very low energy neutrons. It is only a magnitude lower for the neutrons in the 10 to 20 MeV range in the aluminum and two orders lower in the water. The major reason for this difference is in the transport modes: straight-ahead versus angular. The angular treatment allows a longer path length than the straight-ahead approximation. Another reason is the different cross section sets used by the ANISN/PC-BUGLE-80 mode and the HZETRN mode. The next step is to investigate further the differences between the two codes and isolate the differences to just the angular versus straight-ahead transport mode. Then, create a better coupling between the angular neutron transport and the charged particle transport.

Singleterry, R. C., Jr.↗
NASA NTRS1973

Methods of treating complex space vehicle geometry for charged particle radiation transport

Current methods of treating complex geometry models for space radiation transport calculations are reviewed. The geometric techniques used in three computer codes are outlined. Evaluations of geometric capability and speed are provided for these codes. Although no code development work is included several suggestions for significantly improving complex geometry codes are offered.

Hill, C. W.↗
DOE OSTI2025

Studies of particle transport in high-energy-density plasma in the presence of a megagauss magnetic field

Charged particle transport in a magnetic field is among the most fundamental phenomena in plasma physics. Magnetic fields are often introduced in conventional magnetized fusion devices to suppress particle and heat transport, a phenomenon that has been extensively studied. In contrast, particle and heat transport in magnetized high-energy-density (HED) plasmas remain relatively unexplored due to the challenges of generating an external magnetic field strong enough to significantly alter particle dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Resolution requirements for energy conservation in kinetic plasma simulations

The kinetic codes used to model the coupled dynamics of electromagnetic fields and charged particle transport have requirements for spatial, temporal, and charge resolution. These requirements may vary by the solution technique and scope of the problem. In this report, we investigate the resolution limits in the energy-conserving implicit particle-in-cell code CHICAGO. This report has the narrow aim of determining the maximum acceptable grid spacing for the dense plasmas generated in models of z-pinch target gases and power-flow electrode plasmas. In the 2D sample problem, the plasma drifts without external forces with velocity of 10 cm/µs. Simulations are scaled by plasma density to maintain uniform strides across the plasma and from the plasma to the boundaries. Additionally, the cloud-in-cell technique is used with 400 particles per cell and Δ t = 0.85× the Courant limit. For the linear cloud distribution, the criterion for conserving energy is Δ E/E tot < 0.01 for 50,000 time steps. The grid resolution is determined to crudely be Δx ≲ 3 l s , where l s is the electron collisionless skin depth. For the second-order cloud distribution the criterion is Δ E/E tot < 0.005 yielding Δx ≤ 15 l s . These scalings are functions of the chosen v d , Δ t , particles-per-cell, and number of steps.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1985

Numerical and analytic descriptions of cosmic-ray transport

It is not trivial to solve the equations that describe charged particle transport with the aid of computers, for instabilities, inaccuracies, and subtle artifacts are well known afflictions of numerical analysis. Two specific points are discussed. First to avoid inaccuracies, pitch angle scattering must be treated with great care. In particular, slightly inappropriate numerical formulations give rise to mean free paths that are in error by large factors. Secondly, A previously unrecognized artifact, numerical dispersion, is very similar to the physical phenomenon of dispersion. To avoid misinterpretations arising from this similarity, the spatial increment of the finite difference grid must be a small fraction of the mean free path. These points are illustrated by calculations based upon finite difference approximations to the transport equation.

Earl, J. A.↗
DOE OSTI2024

Optical and transport properties of plasma mixtures from ab initio molecular dynamics

Predicting the charged particle transport properties of warm dense matter/hot dense plasma mixtures is a challenge for analytical models. High accuracy ab initio methods are more computationally expensive, but can provide critical insight by explicitly simulating mixtures. In this work, we investigate the transport properties and optical response of warm dense carbon–hydrogen mixtures at varying concentrations under either conserved electronic pressure or mass density at a constant temperature. We compare options for mixing the calculated pure species properties to estimate the results of the mixtures. We find that a combination of the Drude model with the Matthiessen's rule works well for DC electron transport and low-frequency optical response. This breaks down at higher frequencies, where a volumetric mix of pure-species AC conductivities works better.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1994

A unifying comparison of nearly scatter free transport models

Gombosi et al. (1993) recently derived a modified telegrapher's equation for charged particle transport under the influence of isotropic scattering. This equation obeys causality and disallows upstream diffusion for particles with random velocities smaller than the bulk flow velocity. The acausal diffusion equation was obtained to lowest order in the expansion of smallness prameters. The paper by Gombosi et al. (1993) prompted responses from Pauls et al. (1993) and Earl (1993). This paper is written to explain the differences between the methods, assumptions, and results of Gombosi et al. (1993), Pauls et al. (1993), and Earl (1993) and presents a new method of obtaining approximate solutions. It is shown that the assumptions used by Gombosi et al. (1993) and Pauls et al. (1993) are physically equivalent. In our solution method, the solution of the modified telegrapher's equation is obtained as the casual limit of solutions accurate to second order in the smallness parameter expansion. In order to investigate the coherent velocity, we have also developed `wavenumber eigenfunctions' which account for all the pitch angle dependence in our Boltzmann equation. Using truncation, Earl (1993) obtains approximations for the wavenumber dependence of the lowest two frequency modes, which correspond to two of the wavenumber eigenmodes. We find that a consequence of including only two wavenumber eigenmodes is that one obtains solutions which disobey causality at sufficiently short times. Furthermore, the coherent velocity of the two eigenmodes is strongly dependent on wavenumber and approaches the particle velocity in the limit of large wavenumber for both isotropic and anisotropic scattering processes. We conclude that Earl's (1993) solutions and solutions obtained using the new solution method implicitly assume weak acausality and reasonable behavior in the temporal regime, t less than 4 tau. The solutions are not strictly consistent with the behavior of the lowest two frequency modes but have similar behavior in the regime of low wavenumber.

Schwadron, N. A.↗
DOE OSTI2025

Group Conductivity and Nonadiabatic Born Effective Charges of Disordered Metals, Warm Dense Matter, and Hot Dense Plasma

The average ionization state is a critical parameter in plasma models for charged particle transport, equations of state, and optical response. The dynamical or nonadiabatic Born effective charge (NBEC), calculated via first principles time-dependent density-functional theory, provides exact ionic partitioning of bulk electron response for both metallic and insulating materials. The NBEC can be transformed into a “group conductivity,” i.e., the electron conductivity ascribed to a subset of ions. We show that for disordered metallic systems, such as warm dense matter (WDM) and hot dense plasma, the static limit of the NBEC is different from the average ionization states, but that the ionization state can be extracted from the group conductivity even in mixed systems. Here we demonstrate this approach using a set of archetypical examples, including cold and warm aluminium, low- and high-density WDM carbon, and a WDM carbon-beryllium-hydrogen mixture.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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.↗
DOE OSTI2025

The SHMS 11 GeV/c spectrometer in Hall C at Jefferson Lab

The Super High Momentum Spectrometer (SHMS) has been built for Hall C at the Thomas Jefferson National Accelerator Facility (Jefferson Lab). With a momentum capability reaching 11 GeV/ c , the SHMS provides measurements of charged particles produced in electron-scattering experiments using the maximum available beam energy from the upgraded Jefferson Lab accelerator. The SHMS is an ion-optics magnetic spectrometer comprised of a series of new superconducting magnets which transport charged particles through an array of triggering, tracking, and particle-identification detectors that measure momentum, energy, angle and position in order to allow kinematic reconstruction of the events back to their origin at the scattering target. The detector system is protected from background radiation by a sophisticated shielding enclosure. The entire spectrometer is mounted on a rotating support structure which permits measurements to be taken with a large acceptance over laboratory scattering angles from 5.5° to 40°, thus allowing a wide range of low cross-section experiments to be conducted. Finally, these experiments complement and extend the previous Hall C research program to higher energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Multi-fidelity equations of state and transport coefficient datasets for pulsed-power applications

Reliably simulating experiments relevant to the National Nuclear Security Administration (NNSA) requires a detailed description of material properties across a wide range of conditions. Such properties include the equations of state, charged-particle transport coefficients, and optical properties like the opacity. Together, these properties make up the material models used in radiation-magnetohydrodynamic simulations of nuclear fusion experiments. Many of these models do not incorporate uncertainties in the data used to produce them. It is unknown whether these uncertainties significantly impact the interpretation of simulation results and diagnostics. The purpose of this work is to quantify how such uncertainties impact simulations of pulsed-power experiments. We accomplished this task by first assessing discrepancies between approaches used to generate the data. This included bringing together members of the high-energy-density community spanning the three NNSA laboratories and multiple universities. Then, using these data, we developed a general framework that systematically incorporates physical uncertainties within the material models suitable for uncertainty quantification analyses. The framework utilizes machine learning, Bayesian inference, and incorporates multi-fidelity datasets. We demonstrated the framework by quantifying the impact that material model uncertainties have on simulations of pulsed-power experiments underway on Z at Sandia National Laboratories. As a result of this work, we discovered that modest uncertainties in material models (roughly 20%) correspond to significant uncertainties in the outputs from simulations. Our framework has enabled rapid construction of material models through an automated procedure and allows for the generation of material models of interest to the NNSA.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Transport coefficients of warm dense matter from Kohn-Sham density functional theory

We present a comprehensive study of transport coefficients including DC electrical conductivity and related optical properties, electrical contribution to the thermal conductivity, and the shear viscosity via ab initio molecular dynamics and density functional theory calculations on the “priority 1” cases from the “Second Charged-Particle Transport Coefficient Workshop” [Stanek et al., Phys. Plasmas (to be published 2024)]. The purpose of this work is to carefully document the entire workflow used to generate our reported transport coefficients, up to and including our definitions of finite size and statistical convergence, extrapolation techniques, and choice of thermodynamic ensembles. In pursuit of accurate optical properties, we also present a novel, simple, and highly accurate algorithm for evaluating the Kramers–Kronig relations. These heuristics are often not discussed in the literature, and it is hoped that this work will facilitate the reproducibility of our data.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

A Primer on Nuclear "Recoil" Data

For many years the ENDF-6 format has existed to contain evaluated nuclear data. One primary driver for the format and the data is neutron transport calculations. Evaluated data can be processed by a code like NJOY into either continuous energy form (ACE format) for Monte Carlo codes like MCNP or into multi-group form (NDI tables) for deterministic codes like Partisn. The reaction cross sections are found in the MF 3 section of the ENDF-6 format. If the reaction produces one or more neutrons as outputs, then secondary neutron data must also be given in MF 4,5 or 6 sections of the format. Energy and angular distributions for the output neutrons must be given in one of several available formats. The most general formats are found in MF 6. MF 4 is for angular distributions, MF 5 is for Energy distributions, and MF 6 contains both. In recent years, more interest has developed in the other output particles (i.e., the “recoil” particles) from neutron induced reactions. This has been driven by interest in charged particle transport and in more specialized partial kermas. ( A separate total neutron kerma has been available for a long time.) Partial kermas are a breakdown of the Kinetic Energy Released into the MAterial by output particle or by neutron reaction. The sum of the partial kermas should be equal to the total kerma on a group-wise basis. The ENDF-6 format is general enough for these new data requirements. Ideally, evaluated data would exist for every output particle (including the secondary neutrons) from every neutron-induced reaction. In the MF 6 format section, data for multi-particle outputs may be entered using the LAW =1 option. This option explicitly allows energy and angular output distributions for each particle produced in the reaction. Such information should preserve the balance between partial kermas and the total kerma at the groupwise level as well as the individual particle averaged energies. For simpler 2-body reactions, LAW = 2 is available. This allows the evaluation to specify only enough data to specify the 2-body reaction fully. Such a simplification does not exist for 3 (or more) body breakups. A simplified form of LAW 2, i.e., LAW 3, also exists to generate approximate recoil output distributions in the absence of full data. However, evaluated data does not generally exist at this fine granularity for reactions involving more than 2 output particles ( e.g., the 3-body break-up reaction). When the multi-body detailed output distribution data is not available, LAW=6 may be employed in NJOY. LAW 6 produces approximate output distributions for all the particles produced in the reaction. Just like in the 2-body case, the smaller particles will generally carry off more energy. At any particular energy, the balance between the sum of the partial kermas and the total kerma will not necessarily be preserved. However, the average energy of each output particle from a reaction is preserved.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

OpenEdge: A collaborative, open-source, multi-purpose direct simulation Monte Carlo for plasma simulation in magnetic fusion environments

OpenEdge is a collaborative, open-source, object-oriented Direct Simulation Monte Carlo (DSMC) code, designed specifically for plasma simulations in magnetic fusion environments. Here, the code features include advanced structures, robust capabilities, and an effective parallelization strategy, all of which significantly enhance performance. It includes specialized modules for managing complex particle interactions, including collisions, ionization/recombination, and reflection/sputtering. Benchmarks and performance analyses have confirmed its efficiency and scalability. Versatile and adaptable, OpenEdge is applied across a broad spectrum of plasma-material interaction studies and charged particle transport in various fusion research settings.

Boundary plasma↗
DOE OSTI2021

Parameter space for magnetization effects in high-energy-density plasmas

Magnetic fields are well known to affect the evolution of fluids via the J × B force, where J is the current density and B is the magnetic field. This force leads to the influence of magnetic fields on hydrodynamics (magnetohydrodynamics). Magnetic fields are often neglected in modeling of high-energy-density plasmas, since J × B is very small compared with the plasma pressure gradients. However, many experiments lie in a separate part of parameter space where the plasma is indirectly affected via magnetization of the heat flux and charged particle transport. This is true even for initially unmagnetized plasmas, since misaligned density and temperature gradients can self-generate magnetic fields. By comparing terms in the induction equation, we go on to estimate the regions of parameter space where these self-generated fields are strong enough to affect the hydrodynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗