SEARCH · Engineering Papers

Results for “Transport theory”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Intrinsic thermal conductivity of ZrC from low to ultrahigh temperatures: A critical revisit

Current phonon transport theory based on ground-state calculations has been successful in predicting thermal conductivity at room and medium temperatures but may misrepresent behavior at high temperatures. Here, in this work, we predict the thermal conductivity (κ) of ZrC including electronic and phonon contributions from 300 to 3500 K, by including high-order phonon scattering; lattice expansion; temperature-dependent (TD) second-, third-, and fourth-order force constants (2FC, 3FC, and 4FC); and interband phonon conduction by using first principles. For the phonon transport, we find that four-phonon scattering (4ph) significantly reduces the phonon thermal conductivity (κ ph ), by as much as ~75% at 3500 K. After including 4ph scattering and all other factors, κ ph shows a ~T –1.5 rather than ~T –1 dependence. TD 2FC decreases three-phonon scattering rates but increases 4ph rates by decreasing and increasing the scattering phase spaces, respectively. For 4ph phase space, the TD 2FC flattens phonon bands, and allows more redistribution-4ph processes (1 + 2 → 3 + 4) to happen. The combination effect of TD 2FC and TD 4FC reduces 4ph rates of acoustic modes but increases those of optical modes. The TD 3FC and 4FC decrease the phonon scattering cross section and increase the κ ph significantly (by 52% at 3500 K). The contribution from interband (Wigner) phonon conduction is small, even at ultrahigh temperatures. For electronic thermal transport, we find that it is sensitive to and can be changed by 20% by the TD lattice constants. The Lorenz number varies from 1.6 to 3.3 × 10 –8 WΩK –2 at different temperatures. The theoretical prediction in the literature overpredicts κph (e.g., ~28%) and underpredicts the κ el (e.g., ~38%), resulting in an overall underprediction of κ (~26% at 1500 K). The impacts of grain size and defects are found to be strong, leading to the lower observed thermal conductivity in experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Band structure dependent electronic localization in macroscopic films of single-chirality single-wall carbon nanotubes

Significant understanding has been achieved over the last few decades regarding chirality-dependent properties of single-wall carbon nanotubes (SWCNTs), primarily through single-tube studies. However, macroscopic manifestations of chirality dependence have been limited, especially in electronic transport, despite the fact that such distinct behaviors are needed for many applications of SWCNT-based devices. In addition, developing reliable transport theory is challenging since a description of localization phenomena in an assembly of nanoobjects requires precise knowledge of disorder on multiple spatial scales, particularly if the ensemble is heterogeneous. Here, we report an observation of pronounced chirality-dependent electronic localization in temperature and magnetic field dependent conductivity measurements on macroscopic films of single-chirality SWCNTs. The samples included large-gap semiconducting (6,5) and (10,3) films, narrow-gap semiconducting (7,4) and (8,5) films, and armchair metallic (6,6) films. Experimental data and theoretical calculations revealed Mott variable-range-hopping dominated transport in all samples, while localization lengths fall into three distinct categories depending on their band gaps. Armchair films have the largest localization length. Here our detailed analyses on electronic transport properties of single-chirality SWCNT films provide significant new insight into electronic transport in ensembles of nanoobjects, offering foundations for designing and deploying macroscopic SWCNT solid-state devices.

36 MATERIALS SCIENCE↗

DOE OSTI2025

An Atomistic Study of Reactivity in Solid-State Electrolyte Interphase Formation for Li/Li7P3S11

Lithium metal batteries offer superior volumetric and gravimetric specific capacities compared to those based on traditional graphite anodes. Although advancements in solid-state electrolytes address safety concerns, challenges remain, particularly regarding interphase formation in lithium metal anodes. This work presents a computational framework based on high-throughput first-principles density functional theory and machine-learning interatomic potentials (MLIPs) including automated iterative, active learning to enable robust computational exploration of interphase formation between lithium metal anodes and an inorganic solid-state electrolyte. As a demonstration, we apply the framework to a Li/Li7P3S11 interface and find that it accurately identifies the experimentally observed, thermodynamically stable interphase products as well as their overall spatial arrangement within a heterogeneous, amorphous layered structure, with Li2S domains of nanocrystallinity. Our simulations show two stages, a fast and slow diffusion reaction regime, that corroborate the relative phase formation rate of Li x P, Li2S, and Li3P. Using the Onsager transport theory, we capture time-dependent ionic diffusion within the reacting interface, including cross-correlation effects. We found that cross-correlation effects between Li-P and P-S ionic motion significantly influence P-ion diffusion, making it highly sensitive to the local environment and potentially leading to "kinetic trapping" of Li-P phases. The passivation of the interface is shown as the ionic fluxes all approach zero, effectively halting interphase growth.

Diffusion↗

DOE OSTI2024

Parallel Diffusion Coefficient of Energetic Charged Particles in the Inner Heliosphere from the Turbulent Magnetic Fields Measured by Parker Solar Probe

Diffusion coefficients of energetic charged particles in turbulent magnetic fields are a fundamental aspect of diffusive transport theory but remain incompletely understood. In this work, we use quasi-linear theory to evaluate the spatial variation of the parallel diffusion coefficient κ ∥ from the measured magnetic turbulence power spectra in the inner heliosphere. We consider the magnetic field and plasma velocity measurements from Parker Solar Probe made during Orbits 5–13. The parallel diffusion coefficient is calculated as a function of radial distance from 0.062 to 0.8 au, and the particle energy from 100 keV to 1 GeV. We find that κ ∥ increases exponentially with both heliocentric distance and energy of particles. The fluctuations in κ ∥ are related to the episodes of large-scale magnetic structures in the solar wind. By fitting the results, we also provide an empirical formula of κ ∥ = (5.16 ± 1.22) × 10 18 r 1.17 ± 0.08 E 0.71 ± 0.02 (cm 2 s -1 ) in the inner heliosphere, which can be used as a reference in studying the transport and acceleration of solar energetic particles as well as the modulation of cosmic rays.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Refining the two-band model for highly compensated semimetals using thermoelectric coefficients

In studying compensated semimetals, the two-band model has proven extremely useful in capturing electrical conductivity under magnetic field, as a function of density and mobility of electronlike and holelike carriers. However, it rarely offers practical insight into magnetothermoelectric properties. Here, in this work, we report the field dependence of thermoelectric (TE) coefficients in a highly compensated semimetal NbSb 2 , where we find the Seebeck (𝑆 𝑥⁢𝑥 ) and Nernst (𝑆 𝑥⁢𝑦 ) coefficients increase quadratically and linearly with applied magnetic field, respectively. Such field dependence was predicted in previous work that studied a system of two parabolic bands, within semiclassical Boltzmann transport theory when the following two conditions are simultaneously met: 𝜔 𝑐 ⁢𝜏 ≫ 1 and tan ⁡𝜃 𝐻 ≪ 1, where 𝜔 𝑐 and 𝜏 refer to the cyclotron frequency and relaxation time, respectively, and 𝜃 𝐻 is the Hall angle. Under these conditions, we find the field dependence of the TE coefficients directly provides a relation between the electronlike (𝑛 𝑒 ) and holelike (𝑛 ℎ ) carrier densities, which in turn can be used to refine two-band model fitting. With this, we find the compensation factor ($\frac{|Δ⁢𝑛|}{𝑛_{𝑒}}$, where Δ⁢𝑛 = 𝑛 𝑒 −𝑛 ℎ ) of NbSb 2 is two orders of magnitude smaller than what was found in unrestricted fitting, resulting in a larger saturation-field scale for magnetoresistance. Within the same framework of the semiclassical theory, we also deduce that the thermoelectric Hall angle tan⁡ 𝜃 𝛾 = $\frac{𝑆_{𝑥⁢𝑦}}{𝑆_{𝑥⁢𝑥}}$ can be expressed as ($\frac{Δ⁢𝑛}{𝑛_{𝑒}}$ ×𝜔 𝑐 ⁢𝜏) −1 , which serves as a parameter to predict the degree of compensation. Our findings offer crucial insights into identifying empirical conditions for field-induced enhancement of TE performance and into engineering-efficient thermoelectric devices based on semimetallic materials.

electrical conductivity↗

DOE OSTI2026

The Physics Imposed on a Streaming Operator by Spherical Transport Problems

The streaming operator, which generates a displacement of a particle on a straight line at a constant speed in transport theory, is derived algebraically from a spherical coordinate formulation of Newton’s second law. This derivation leads to an operator that has more partial derivatives than a Cartesian coordinate formulation of the operator. The additional partial derivatives, which are with respect to the normalized velocity variables of a particle, take into account the intrinsic curvature of a ball. Moreover, these partial derivatives mitigate ray effects, which arise when a finite number of normalized velocities (also called directions or discrete ordinates) are used to simulate a continuous S 2 sphere of directions, by rotating the polar axis of the S 2 sphere into the radial direction of the coordinate system. As a consequence of this rotation, the number of actual discrete ordinates is greatly amplified to an enormous number of effective discrete ordinates by a multiplier that is equal to the number of patches that partitions a spherical surface. In addition to the derivation of the streaming operator, we provide in closed form a solution to the system of characteristic equations that is equivalent to the streaming operator. Furthermore, the solution to the system of characteristic equations enables the construction of an integral operator that is the inverse to the streaming operator. Examples in which ray effects are immensely mitigated by spherical coordinates are presented.

integral operator↗

DOE OSTI2023

Synthesis and Observation of Emergent Phenomena in Epitaxial Heusler Compound Heterostructures

The proposal was on the synthesis and observation of emergent phenomena in epitaxial Heusler compound heterostructures. The large range of properties and number of Heusler compounds opens up a wide number of potential compounds that will exhibit emergent phenomena. The similarity, large range of relatively inexpensive, large area, high crystal quality, III-V bulk substrates, lattice parameters and the ability to tune the lattice parameters through ternary or quaternary III-V compound semiconductor epitaxial growth, makes III-V semiconductors an ideal choice for substrates for epitaxial growth of Heusler compounds. A number of half Heusler compounds have been predicted to exhibit band inversion, making them topological and are therefore expected to exhibit spin-momentum locked topological surface states with linear dispersion. Others are predicted to be semimetals with Weyl points and others semiconducting and magnetic. During the course of this grant, emphasis has been on investigating Heusler compounds with emergent phenomena and demonstrating the ability to tune their properties through alloying and strain. We have grown toplogical semimetal (PtLuSb, PtMnBi), Weyl (Co 2 MnAl, Co 2 TiGe), half metal (PtMnSb, Co 2 MnSi, Co 2 MnAl x Si 1-x , Co2FeAl), and semiconducting (CoTiSb, NiTiSn) and tuned their properties through alloying and epitaxial strain. We also investigated the closely related materials of rare-earth monopnictide, some of which have also been predicted to be topological. During the attempts to grow the PtMnBi, it was discovered that Bi, another predicted topological material when ultrathin, could be grown epitaxially on InSb, results for which are also reported here. The main focus for this effort has been on using variable photon energy and spin-dependent angle resolved photoemission (ARPES) to determine bulk band structure and surface states of pristine epitaxial films grown on III-V semiconductor and MgO substrates and correlate results with theory and transport measurements. Theory has been critical to interpretation of experimental results and has been essential in guiding experiments. The research benefited from several strong collaborations between the PIs and the beamline scientists at the Advanced Light Source at Lawrence Berkeley Laboratory, the Stanford Linear Accelerator Center (SLAC) at Stanford and at the Max Lab at Lund University in Sweden. The strong experiment - theory collaboration between the PI’s groups, the Palmstrøm group at UCSB and the Janotti group at the University of Delaware, has been critical for interpreting the experimental ARPES and magnetotransport measurements results and making predictions to guide experiments. Weekly interactive Zoom meetings made this work well. A collaboration between the Palmstrøm group and Dr. Alexei Fedorov at the Advanced Light Source (ALS) resulted in significant modifications to his end chamber to accommodate the vacuum suitcase that was designed and constructed in the Palmstrøm group at UCSB. In collaboration with beamline scientists, Drs. Makoto Hashimoto and Donghui Lu at SLAC, Palmstrøm made modifications to the vacuum suitcase and developed special sample holders that allowed samples to be grown in the Palmstrøm MBE systems at UCSB and transported in the UHV vacuum suitcase to SLAC for ARPES measurements. The development of the vacuum suitcase was essential for this grant as it has allowed variable photon energies to be used to identify surface versus bulk states on samples that could not be capped and decapped using As- or Sb-capping layers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Effects of aluminum content on thermoelectric performance of Al x CoCrFeNi high-entropy alloys

Introducing a non-regular distribution in the mass and bonding by including distinctly different elements can reduce the phonon transport even within structurally well-ordered materials. These distributions are a quality of all high-entropy alloys (HEAs), however the inclusion of aluminum in Al x CoCrFeNi is particularly impactful due to the large mismatch in atomic mass with other components. The resultant low phonon conductivity is a requirement for high thermoelectric performance, motivating the investigation of the effects of Al content on phonon transport as well as other thermoelectric properties. This work examines the phonon and electron transport and thermoelectric conversion properties with various Al contents (0 ≤ x Al ≤ 2) in this Cantor alloy using first-principles calculations, molecular dynamics, and semi-classical Boltzmann transport theory. The calculated phonon density of states and thermoelectric properties present reasonable agreements with experiments, including neutron scattering. A large reduction of phonon conductivity (k L ) is observed even with low x Al s, which we attribute to effective phonon scatterings by the large mass mismatch. However, its temperature dependence is not significant, demonstrating a minor contribution of interphonon scattering. In contrast, electrical conductivity (σ) and Seebeck coefficient (S) increase with temperature at higher x Al s with body-centered cubic structures. Therefore, the thermoelectric figure of merit (ZT) of Al x CoCrFeNi HEAs is enhanced by increasing the Al content mainly due to the increase of the thermoelectric power factor (σS 2 ) at high temperatures, while at low temperatures the phonon-scattering enhancement by mass mismatch is also important.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Connection Between Mass Flux Transport and Eddy Diffusivity in Convective Atmospheric Boundary Layers

Turbulence parameterizations for convective boundary layer in coarse-scale atmospheric models usually consider a combination of the eddy-diffusive transport and a non-local transport, typically in the form of a mass flux term, such as the widely adopted eddy-diffusivity mass-flux (EDMF) approach. These two types of turbulent transport are generally considered to be independent of each other. Using results from large-eddy simulations, here, we show that a Taylor series expansion of the updraft and downdraft mass-flux transport can be used to approximate the eddy-diffusivity transport in the atmospheric surface layer and the lower part of the mixed layer, connecting both eddy-diffusivity and mass-flux transport theories in convective conditions, which also quantifies departure from the Monin-Obukhov similarity (MOS) in the surface layer. This study provides a theoretical support for a unified EDMF parameterization applied to both the surface layer and mixed layer and highlights important correction required for surface models relying on MOS

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Collisionless cooling of perpendicular electron temperature in the thermal quench of a magnetized plasma

Thermal quench of a nearly collisionless plasma against an isolated cooling boundary or region is an undesirable off-normal event in magnetic fusion experiments, but an ubiquitous process of cosmological importance in astrophysical plasmas. Parallel transport theory of ambipolar-constrained tail electron loss is known to predict rapid cooling of the parallel electron temperature $T_{e\Vert}$ although $T_{e\Vert}$ is difficult to diagnose in actual experiments. Instead direct experimental measurements can readily track the perpendicular electron temperature $T_{e\bot}$ via electron cyclotron emission. The physics underlying the observed fast drop in $T_{e\bot}$ requires a resolution. Here two collisionless mechanisms, dilutional cooling by infalling cold electrons and wave-particle interaction by two families of whistler instabilities, are shown to enable fast $T_{e\bot}$ cooling that closely tracks the mostly collisionless crash of $T_{e\Vert}$. These findings motivate both experimental validation and reexamination of a broad class of plasma cooling problems in laboratory, space, and astrophysical settings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Selective permeation up a chemical potential gradient to enable an unusual solvent purification modality

The need for energy-efficient recovery of organic solutes from aqueous streams is becoming more urgent as chemical manufacturing transitions toward nonconventional and bio-based feedstocks and processes. In addition to this, many aqueous waste streams contain recalcitrant organic contaminants, such as pharmaceuticals, industrial solvents, and personal care products, that must be removed prior to reuse. We observe that rigid carbon membrane materials can remove and concentrate organic contaminants via an unusual liquid-phase membrane permeation modality. Surprisingly, detailed thermodynamic calculations on the chemical potential of the organic contaminant reveal that the organic species has a higher chemical potential on the permeate side of the membrane than on the feed side of the membrane. This unusual observation challenges conventional membrane transport theory that posits that all permeating species move from high chemical potential states to lower chemical potential states. Based on experimental measurements, we hypothesize that the organic is concentrated in the membrane relative to water via favorable binding interactions between the organic and the carbon membrane. The concentrated organic is then swept through the membrane via the bulk flow of water in a modality known as “sorp-vection.” We highlight via simplified nonequilibrium thermodynamic models that this “uphill” chemical potential permeation of the organic does not result in second-law violations and can be deduced via measurements of the organic and water sorption and diffusion rates into the carbon membrane. Moreover, this work identifies the need to consider such nonidealities when incorporating unique, rigid materials for the separations of aqueous waste streams.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Space-Dependent Calculation of the Multiplicity Moments for Shells With the Inclusion of Scattering

In recent work, we extended the methodology of multiplicity counting in nuclear safeguards by elaborating the one-speed stochastic transport theory of the calculation of the so-called multiplicity moments, i.e., the factorial moments of the number of neutrons emitted from a fissile item, following a source event from an internal neutron source [spontaneous fission and (α; n) reactions]. Calculations were made for solid spheres and cylinders, with the source being homogeneously distributed within the item. Recent measurements of the Rocky Flats Shells during the Measurement of Uranium Subcritical and Critical (MUSIC) campaign conducted by Los Alamos National Laboratory and assisted by the University of Michigan inspired us to extend the model to spherical shell geometry with a point source in the middle of the central cavity. Comparison of the calculated results with the experimental ones indicated that accounting for fission as the only neutron reaction (the standard procedure in the point model, adapted also in our work so far) was not sufficient for reaching good agreement with measurements. The model was therefore extended to include elastic scattering into the one speed formalism, whereas the effect of inelastic scattering was accounted for in an empirical way. After these extensions, good agreement was found between the calculated and the measured values. The paper describes the extension of the theory and provides concrete quantitative results.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Prehydrodynamic evolution and its signatures in final-state heavy-ion observables

We investigate the effects of prehydrodynamic evolution on final-state observables in heavy-ion collisions using state-of-the art event simulations coupled to different prehydrodynamic scenarios, which include the recently developed effective kinetic transport theory evolution model KøMPøST. Differential flow observables are found to be mostly insensitive to the details of prehydrodynamic evolution. The main effect we observe is in the p T spectra, particularly the mean transverse momentum. However, at least part of this effect is a consequence of the underlying conformal invariance assumption currently present in such approaches, which is known to be violated in the temperature regime probed in heavy-ion collisions. This assumption of early time conformal invariance leads to an artificially large out-of-equilibrium bulk pressure when switching from (conformal) prehydrodynamic evolution to hydrodynamics (using the nonconformal QCD equation of state), which in turn increases the transverse momentum. Our study indicates that a consistent treatment of prehydrodynamic evolution in heavy-ion collisions requires the use of nonconformal models of early time dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Phase-space methods for neutrino oscillations: Extension to multibeams

The phase-space approach (PSA), which was originally introduced in Lacroix [] to describe neutrino flavor oscillations for interacting neutrinos emitted from stellar objects is extended to describe arbitrary numbers of neutrino beams. The PSA is based on mapping the quantum fluctuations into a statistical treatment by sampling initial conditions followed by independent mean-field evolution. A new method is proposed to perform this sampling that allows treating an arbitrary number of neutrinos in each neutrino beams. We validate the technique successfully and confirm its predictive power on several examples where a reference exact calculation is possible. We show that it can describe many-body effects, such as entanglement and dissipation induced by the interaction between neutrinos. Due to the complexity of the problem, exact solutions can only be calculated for rather limited cases, with a limited number of beams and/or neutrinos in each beam. The PSA approach considerably reduces the numerical cost and provides an efficient technique to accurately simulate arbitrary numbers of beams. Examples of PSA results are given here, including up to 200 beams with time-independent or time-dependent Hamiltonians. We anticipate that this approach will be useful to bridge exact microscopic techniques with more traditional transport theories used in neutrino oscillations. It will also provide important reference calculations for future quantum computer applications where other techniques are not applicable to classical computers. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2023

Shifting and Splitting of Resonance Lines due to Dynamical Friction in Plasmas

We report a quasilinear plasma transport theory that incorporates Fokker-Planck dynamical friction (drag) and pitch angle scattering is self-consistently derived from first principles for an isolated, marginally unstable mode resonating with an energetic minority species. It is found that drag fundamentally changes the structure of the wave-particle resonance, breaking its symmetry and leading to the shifting and splitting of resonance lines. In contrast, scattering broadens the resonance in a symmetric fashion. Comparison with fully nonlinear simulations shows that the proposed quasilinear system preserves the exact instability saturation amplitude and the corresponding particle redistribution of the fully nonlinear theory. Even in situations in which drag leads to a relatively small resonance shift, it still underpins major changes in the redistribution of resonant particles. This novel influence of drag is equally important in plasmas and gravitational systems. In fusion plasmas, the effects are especially pronounced for fast-ion-driven instabilities in tokamaks with low aspect ratio or negative triangularity, as evidenced by past observations. The same theory directly maps to the resonant dynamics of the rotating galactic bar and massive bodies in its orbit, providing new techniques for analyzing galactic dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Experimental Measurements of Ion Diffusion Coefficients and Heating in a Multi-Ion-Species Plasma Shock

Collisional plasma shocks generated from supersonic flows are an important feature in many astrophysical and laboratory high-energy-density plasmas. Compared to single-ion-species plasma shocks, plasma shock fronts with multiple ion species contain additional structure, including interspecies ion separation driven by gradients in species concentration, temperature, pressure, and electric potential. We present time-resolved density and temperature measurements of two ion species in collisional plasma shocks produced by head-on merging of supersonic plasma jets, allowing determination of the ion diffusion coefficients. Here our results provide the first experimental validation of the fundamental inter-ion-species transport theory. The temperature separation, a higher-order effect reported here, is valuable for advancements in modeling HED and ICF experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Invariant Forms of Dissolution Fingers

Dissolution of fractured and porous media introduces a positive feedback between fluid flow and reactant transport, leading to the emergence of pronounced, fingerlike channels. We investigate the formation of these structures using a microfluidic Hele-Shaw cell with a soluble bottom. Our experiments show that the shape of dissolution fingers is invariant and reveals itself over time as the fingers extend into the system. Here, by combining reactive-transport theory and conformal mapping techniques we derive these invariant forms. We relate these results to natural dissolution fingers in karst landscapes, and illustrate how to determine the groundwater flow rate responsible for their formation based on the finger shape.

Flows in porous media↗

DOE OSTI2020

LA17-V-InverseRadTransportProblem-PD2Kb. Final Project Report

In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.

07 ISOTOPE AND RADIATION SOURCES↗