Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Methods in 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Conditional Point Sampling: A Monte Carlo Method for Radiation Transport in Stochastic Media.

Current methods for stochastic media transport are either computationally expensive or, by nature, approximate. Moreover, none of the well-developed, benchmarked approximate methods can compute the variance caused by the stochastic mixing, a quantity especially important to safety calculations. Therefore, we derive and apply a new conditional probability function (CPF) for use in the recently developed stochastic media transport algorithm Conditional Point Sampling (CoPS), which 1) leverages the full intra-particle memory of CoPS to yield errorless computation of stochastic media outputs in 1D, binary, Markovian-mixed media, and 2) leverages the full inter-particle memory of CoPS and the recently developed Embedded Variance Deconvolution method to yield computation of the variance in transport outputs caused by stochastic material mixing. Numerical results demonstrate errorless stochastic media transport as compared to reference benchmark solutions with the new CPF for this class of stochastic mixing as well as the ability to compute the variance caused by the stochastic mixing via CoPS. Using previously derived, non-errorless CPFs, CoPS is further found to be more accurate than the atomic mix approximation, Chord Length Sampling (CLS), and most of memory-enhanced versions of CLS surveyed. In addition, we study the compounding behavior of CPF error as a function of cohort size (where a cohort is a group of histories that share intra-particle memory) and recommend that small cohorts be used when computing the variance in transport outputs caused by stochastic mixing.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2020

Oak Ridge National Laboratory Shutdown Dose Rate Code Suite

The Oak Ridge National Laboratory (ORNL) shutdown dose rate (SDDR) code suite calculates gamma dose rates at a location of interest due to activation of materials. The code suite is based on the rigorous two-step (R2S) method that involves two radiation transport calculations: (1) neutron radiation transport to determine neutron fluxes in the material for which activation calculations will be performed and (2) gamma radiation transport to determine dose rates at a location of interest, due to material activation. The ORNL SDDR code suite is used to determine an importance function that characterizes the neutron’s importance to the final SDDR when the Monte Carlo radiation transport method is used. This report demonstrates how the ORNL SDDR code suite is implemented through the following codes: Monte Carlo N-Particle (MCNP) code Version 5-1.60, Oak Ridge Isotope Generator (ORIGEN), MSX suite of utilities, Neutron Activation Gamma Source Sampler (NAGSS), and Automated Variance Reduction Generator (ADVANTG). For this demonstration, the International Thermonuclear Experimental Reactor (ITER) SDDR benchmark problem is used and the dose rate result is compared to previously published work.

61 RADIATION PROTECTION AND DOSIMETRY↗
NASA NTRS1982

A comparative study of computational methods in cosmic gas dynamics

Many theoretical investigations of fluid flows in astrophysics require extensive numerical calculations. The selection of an appropriate computational method is, therefore, important for the astronomer who has to solve an astrophysical flow problem. The present investigation has the objective to provide an informational basis for such a selection by comparing a variety of numerical methods with the aid of a test problem. The test problem involves a simple, one-dimensional model of the gas flow in a spiral galaxy. The numerical methods considered include the beam scheme, Godunov's method (G), the second-order flux-splitting method (FS2), MacCormack's method, and the flux corrected transport methods of Boris and Book (1973). It is found that the best second-order method (FS2) outperforms the best first-order method (G) by a huge margin.

Van Albada, G. D.↗
NASA NTRS2008

Neutron Transport Models and Methods for HZETRN and Coupling to Low Energy Light Ion Transport

Exposure estimates inside space vehicles, surface habitats, and high altitude aircraft exposed to space radiation are highly influenced by secondary neutron production. The deterministic transport code HZETRN has been identified as a reliable and efficient tool for such studies, but improvements to the underlying transport models and numerical methods are still necessary. In this paper, the forward-backward (FB) and directionally coupled forward-backward (DC) neutron transport models are derived, numerical methods for the FB model are reviewed, and a computationally efficient numerical solution is presented for the DC model. Both models are compared to the Monte Carlo codes HETCHEDS and FLUKA, and the DC model is shown to agree closely with the Monte Carlo results. Finally, it is found in the development of either model that the decoupling of low energy neutrons from the light ion (A<4) transport procedure adversely affects low energy light ion fluence spectra and exposure quantities. A first order correction is presented to resolve the problem, and it is shown to be both accurate and efficient.

Blattnig, S.R.↗
DOE OSTI2025

Automated Hybrid Variance Reduction on Advanced Architectures in the Shift Monte Carlo Code

Monte Carlo transport methods are the most accurate schemes for solving problems with complex energy and spatial features, but they come with a high computational cost. Although hybrid methods have enabled the use of Monte Carlo transport for a large class of problems, they still require significant computing resources. Modern multicore CPUs with large numbers of compute cores and graphical processing units (GPUs) provide opportunities to optimize the memory and run-time costs of hybrid Monte Carlo methods. This paper documents the development and analysis of three Monte Carlo transport algorithms that support hybrid transport using the consistent adjoint-driven importance sampling (CADIS) and forward-weighted CADIS methods in the Shift Monte Carlo code: history-based transport using static and dynamic threading on multicore CPUs and event-based transport enabling weight window tracking on GPUs. The results are shown for two challenging hybrid problems on the Frontier supercomputer at the Oak Ridge Leadership Computing Facility. The results show that all three methods yield good performance and enable solutions of difficult fixed-source transport problems in less than 2 min on 20 nodes of Frontier. Dynamic threading was observed to give up to 20% better scaling behavior than static threading. Moreover, the AMD Instinct 250X GPU was found to give 9 to 11 times greater throughput per graphics compute die than the best CPU performance. In conclusion, additional opportunities for optimization of hybrid transport on GPUs are discussed.

Denovo↗
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 OSTI2020

Comparison of Generation of Higher-Order Neutron Scattering Cross Sections

The generation of high-order neutron scattering cross sections consistent with high-fidelity simulations remains an area of active research. Popular options include generating cross sections from continuous energy Monte Carlo calculations or generating cross sections from a deterministic neutron transport calculation with high-fidelity tabulated cross sections. Both options present challenges. Monte Carlo simulations can naturally process continuous energy cross section data and allow for the general description of anisotropic neutron scattering given a scattering law. However, Monte Carlo simulations are inherently related to particle weighting and it has been suggested that this may be unacceptable for generating high-order neutron scattering cross sections. Deterministic neutron transport calculations can easily calculate high-order moments of the flux to appropriately calculate higher-order neutron scattering cross sections but are generally limited by discretization of space, energy, and angle. In this work, the trade-offs between generation of high-order neutron scattering cross sections via Monte Carlo and deterministic neutron transport methods are investigated. The methods implemented in the Monte Carlo computer program Serpent and the deterministic fast reactor neutron cross section generator MC2-3 are compared. Cross sections resulting from these methods are used in Rattlesnake, a deterministic neutron transport code developed by Idaho National Laboratory (INL), and results are compared to a reference continuous energy Monte Carlo calculation. Whereas previous work investigating the effects of anisotropic neutron scattering has focused on light water reactor (LWR) simulations, this work focuses on high-order neutron scattering cross sections as they relate to fast reactor simulations. To investigate the consequences of the Serpent and MC2-3 methodologies, a test problem is developed. The test problem is a one-dimensional geometry with fast reactor materials designed to demonstrate deep penetration and exacerbate the effects of high-order neutron scattering. Based on the results of the deep penetration test problem, it is concluded that P3 neutron scattering cross sections are sufficient to describe anisotropic scattering in fast reactor materials. Any neutron scattering of order higher than P3 offers negligible change in the eigenvalue of the test problem. Additionally, it is determined that the methodology as implemented in Serpent is applicable for generating high-order neutron scattering cross sections through at least P3 in fast reactor materials.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Early Career Award: Single Molecule Fluorescence Imaging for a Background-Free Neutrinoless Double Beta Decay Search

This project aimed to develop a suite of novel technologies targeted at detecting individual barium ions produced in neutrinoless double beta decays of 136 Xe, as part of the NEXT program of gaseous xenon detectors. We realized novel molecules which could sense Ba 2+ ions in dry conditions; new microscopy techniques operable in dry noble gases, including the required operating conditions of xenon at 10 bar pressure; and new radiofrequency ion transport methods (developed both theoretically and experimentally) that can transport ions in dense xenon gas. Together, these developments represent a major advance toward a background-free neutrinoless double beta decay experiment based on single molecule fluorescence imaging of Ba 2+ daughter ions.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Vapor phase transport system and method for depositing perovskite semiconductors

Vapor phase transport systems and methods of depositing perovskite films are described. In an embodiment, a deposition method includes feeding a perovskite solution or constituent powder to a vaporizer, followed by vaporization and depositing the constituent vapor as a perovskite film. In an embodiment, a deposition system and method includes vaporizing different perovskite precursors in different vaporization zones at different temperatures, followed by mixing the vaporized precursors to form a constituent vapor, and depositing the constituent vapor as a perovskite film.

Bush, Kevin M. Alexander↗
DOE OSTI2022

Vapor transport growth of MnBi 2 Te 4 and related compounds

Motivated by fine tuning of the magnetic and topological properties of MnBi 2 Te 4 via defect engineering, in this work, we report the crystal growth of MnBi 2 Te 4 and related compounds using vapor transport method and crystal characterization by measuring elemental ratio, magnetic and transport properties, and scanning tunneling microscopy. For the growth of MnBi 2 Te 4 single crystals, I 2 , MnI 2 , MnCl 2 , TeCl 4 , and MoCl 5 are all effective transport agents; chemical transportation occurs faster in the presence of iodides than chlorides. We further successfully grow MnSb 2 Te 4 , MnBi 2-x Sb x Te 4 , and Sb-doped MnBi 4 Te 7 crystals. A small temperature gradient< 20°C between the hot and cold ends of the growth ampoule is critical for the successful crystal growth of MnBi 2 Te 4 and related compounds. Compared to flux grown crystals, vapor transported crystals tend to be Mn stoichiometric, and Sb-bearing compositions have more Mn/Sb site mixing. The vapor transport growth provides a new materials synthesis approach to fine tune the magnetic and topological properties of these intrinsic magnetic topological insulators where controlling defects is vital.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Phonon transport governed by intrinsic scattering in short-period AlN/GaN superlattices

Here we employ density functional theory based phonon transport methods to provide a rigorous understanding of the nature of thermal transport in coherent short-period AlN/GaN superlattices (SLs), with period lengths up to three unit cells of each, and compare these with properties of their bulk constituents. Increasing the period length leads to phonon band folding with frequency gaps and thus smaller phonon velocities and more phonon scattering than in bulk, both of which reduce lattice thermal conductivity (κ). Contrary to expectations, we find that velocity variations among larger-period AlN/GaN SLs play only a minor role in cross-plane κ reductions, while variations in intrinsic phonon scattering are strongly correlated with their transport behaviors. This work provides insights into the microscopic behaviors of technologically relevant nanostructured materials, which are likely applicable to a wider range of SL systems and other nanostructures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1992

Ground-based measurements of galactic cosmic ray fragmentation in shielding

The mean free path for nuclear interactions of galactic cosmic-rays is comparable to shielding and tissue thicknesses present in human interplanetary exploration, resulting in a significant fraction of nuclear reaction products at depth. In order to characterize the radiation field, the energy spectrum, the angular distribution, and the multiplicity of each type of secondary particles must also be known as a function of depth. Reactions can take place anywhere in a thick absorber; therefore, it is necessary to know these quantities as a function of particle energy for all particles produced. HZE transport methods are used to predict the radiation field; they are dependent on models of the interaction of man-made systems with the space environment to an even greater extent than methods used for other types of radiation. Hence, there is a major need to validate these transport codes by comparison with experimental data. The most cost-effective method of validation is a comparison with ground-based experimental measurements. A research program to provide such validation measurements using neon, iron and other accelerated heavy ion beams will be discussed and illustrated using results from ongoing experiments and their comparison with current transport codes. The extent to which physical measurements yield radiobiological predictions will be discussed.

NASA Discipline Radiation Health↗
NASA NTRS1981

Hot oxygen atoms in the upper atmosphere of Venus

Using data obtained by an ultraviolet spectrometer on the Pioneer-Venus Orbiter, energy and altitude distribution of nonthermal oxygen atoms in the Venus atmosphere for typical daytime and nighttime conditions are calculated using two different numerical methods. Agreement is found between results obtained by the two-stream transport method and the diffusion equation technique. Dissociative recombination of molecular oxygen ions and charge exchange of atomic oxygen ions with the neutral hydrogen and oxygen gas are the source terms, and calculations show that the dissociative recombination term is the dominant one. It is shown that measured hot atom concentrations exceed estimates of the thermal hydrogen atom concentrations derived from Brinton et al. (1980) over most of the daytime exosphere, and the calculated concentrations do so to an even greater extent. It is concluded that hot oxygen atoms are an important component of the dayside exosphere of Venus, their most important effect being the absorption of solar wind.

Nagy, A. F.↗
DOE OSTI2021

Constructing angular distributions of neutrinos in core-collapse supernovae from zeroth and first moments calibrated by full Boltzmann neutrino transport

Two-moment neutrino transport methods have been widely used for developing theoretical models of core-collapse supernovae (CCSN), since they substantially reduce the computational burden inherent in the multidimensional neutrino-radiation hydrodynamical simulations. The approximation, however, comes at a price; the detailed structure of angular distribution of neutrinos is sacrificed, that is the main drawback of this approach. In this paper, we develop a novel method by which to construct angular distributions of neutrinos from the zeroth and first angular moments. In our method, the angular distribution is expressed with two quadratic functions of the neutrino angle in a piecewise fashion. We determine the best parameters in the fitting function by comparing to the neutrino data in a spherically symmetric CCSN model with full Boltzmann neutrino transport. We demonstrate the capability of our method by using our recent 2D CCSN model. In this work, we find that the essential features of the angular distributions can be well reconstructed, whereas the angular distributions of incoming neutrinos tend to have large errors that increase with flux factor ($\textit{κ}$). This issue originates from the insensitiveness of incoming neutrinos to $\textit{κ}$, that is an intrinsic limitation in moment methods. Based on the results of the demonstration, we assess the reliability of electron- neutrinos lepton number (ELN)-crossing searches with the two-moment neutrino transport. This analysis is complementary to our [L. Johns and H. Nagakura arXiv:2104.04106] that scrutinizes the limitation of crossing searches with a few moments. We find that the systematic errors of angular distributions for incoming neutrinos lead to misjudgements of the crossing at $\textit{κ}$ ≳ 0.5. This casts doubt on the results of ELN-crossing searches based on two-moment methods in some previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗