Results for “obtaining transport coefficients” — Page 4

DOE OSTI2021

Efficient calculation of carrier scattering rates from first principles

Abstract The electronic transport behaviour of materials determines their suitability for technological applications. We develop a computationally efficient method for calculating carrier scattering rates of solid-state semiconductors and insulators from first principles inputs. The present method extends existing polar and non-polar electron-phonon coupling, ionized impurity, and piezoelectric scattering mechanisms formulated for isotropic band structures to support highly anisotropic materials. We test the formalism by calculating the electronic transport properties of 23 semiconductors, including the large 48 atom CH 3 NH 3 PbI 3 hybrid perovskite, and comparing the results against experimental measurements and more detailed scattering simulations. The Spearman rank coefficient of mobility against experiment ( r s = 0.93) improves significantly on results obtained using a constant relaxation time approximation ( r s = 0.52). We find our approach offers similar accuracy to state-of-the art methods at approximately 1/500th the computational cost, thus enabling its use in high-throughput computational workflows for the accurate screening of carrier mobilities, lifetimes, and thermoelectric power.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Cycles-L: A Coupled, 3-D, Land Surface, Hydrologic, and Agroecosystem Landscape Model

Managing landscapes to increase agricultural productivity and environmental stewardship can be informed by spatially-distributed models that operate at spatial and temporal scales that are intervention-relevant. This paper presents Cycles-L, a landscape-scale agroecosystem and hydrologic modeling system, using as a test case a watershed in Pennsylvania. Cycles-L emerges from melding the landscape and hydrology structure of Flux-PIHM, a 3-D land surface hydrologic model, with the agroecosystem processes in the Cycles model. Consequently, Cycles-L can simulate processes affected by topography, soil heterogeneity, and management practices, owing to its physically-based hydrology that can simulate horizontal and vertical transport of solutes with water. The model was tested at a 730-ha experimental watershed within the Mahantango Creek watershed. Cycles-L simulated well stream water and mineral nitrogen discharge (Nash-Sutcliffe coefficient 0.55 and 0.60, respectively) and grain yield (root mean square error 1.2 Mg ha -1 ). Cycles-L outputs are as good or better than those obtained with the uncoupled Flux-PIHM (water discharge) and Cycles (grain yield) models. Modeled spatial patterns of nitrogen fluxes like denitrification illustrate the combined control of crop management and topography. For example, denitrification is almost twice as high when simulated with Cycles-L than when simulated with Cycles 1-D. Due to its spatial and temporal resolution, Cycles-L fills a gap in the availability of models that operate at a scale relevant to evaluate interventions in the landscape. Cycles-L can become a central component in tools for climate change scenario analysis, precision agriculture, precision conservation, and artificial intelligence-based decision support systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

A model for discrete fracture-clay rock interaction incorporating electrostatic effects on transport

Abstract A model based on the code CrunchClay is presented for a fracture-clay matrix system that takes electrostatic effects on transport into account. The electrostatic effects on transport include those associated with the development of a diffusion potential as captured by the Nernst-Planck equation, and the formation of a diffuse layer bordering negatively charged clay particles within which partial anion exclusion occurs. The model is based on a dual continuum formulation that accounts for diffuse layer and bulk water pore space, providing a more flexible framework than is found in the classical mean electrostatic potential models. The diffuse layer model is obtained by volume averaging ion concentrations in the Poisson-Boltzmann equation, but also includes the treatment of longitudinal transport within this continuum. The calculation of transport within the bulk and diffuse layer porosity is based on a new formulation for the Nernst-Planck equation that considers averaging of diffusion coefficients and accumulation factors at grid cell interfaces. Equations for function residuals and the associated Jacobian matrix are presented such that the system of nonlinear differential-algebraic equations can be solved with Newton’s method. As an example, we consider a 2D system with a single discrete fracture within which flow and advective transport occurs that is coupled to diffusion in the clay-rich matrix. The simulation results demonstrate the lack of retardation for anions (e.g., 36 Cl − ) of the contaminant plume within the fracture flow system because they are largely excluded from the charged clay rock, while the migration of cations (e.g., 90 Sr ++ ) is more strongly attenuated. The diffusive loss of divalent cations in particular from the fracture is accentuated by their accumulation in the diffuse layer within the clay-rich matrix.

58 GEOSCIENCES↗

DOE OSTI2022

Robust estimation of charge carrier diffusivity using transient photoluminescence microscopy

Transient microscopy has emerged as a powerful tool for imaging the diffusion of excitons and free charge carriers in optoelectronic materials. In many excitonic materials, extraction of diffusion coefficients can be simplified because of the linear relationship between signal intensity and local excited state population. However, in materials where transport is dominated by free charge carriers, extracting diffusivities accurately from multidimensional data is complicated by the nonlinear dependence of the measured signal on the local charge carrier density. To obtain accurate estimates of charge carrier diffusivity from transient microscopy data, statistically robust fitting algorithms coupled to efficient 3D numerical solvers that faithfully relate local carrier dynamics to raw experimental measurables are sometimes needed. Here, we provide a detailed numerical framework for modeling the spatiotemporal dynamics of free charge carriers in bulk semiconductors with significant solving speed reduction and for simulating the corresponding transient photoluminescence microscopy data. To demonstrate the utility of this approach, we apply a fitting algorithm using a Markov chain Monte Carlo sampler to experimental data on bulk CdS and methylammonium lead bromide (MAPbBr 3 ) crystals. Parameter analyses reveal that transient photoluminescence microscopy can be used to obtain robust estimates of charge carrier diffusivities in optoelectronic materials of interest, but that other experimental approaches should be used for obtaining carrier recombination constants. Additionally, simplifications can be made to the fitting model depending on the experimental conditions and material systems studied. Our open-source simulation code and fitting algorithm are made freely available to the scientific community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Global Sensitivity Analysis of a Reactive Transport Model for Mineral Scale Formation During Hydraulic Fracturing

Injection of water-based hydraulic fracturing fluid (HFF) into tight shale gas/oil formations can increase formation permeability and enhance production rates, but this process frequently causes mineral scale formation that can occlude pore space and hinder flow. To identify the most important factors that control the formation of mineral scales, we applied a novel global sensitivity analysis method—distance-based generalized sensitivity analysis (DGSA)—to a reactive transport model (RTM) that was previously built and calibrated to simulate precipitation of barite [BaSO4] and iron (hydr)oxide [Fe(OH) 3 ] in shale matrices and on fracture surfaces. Reactive transport simulations were run with model parameters randomly sampled based on assigned uncertainties. Modeling results for barite and Fe(OH)3 formation were clustered using machine-learning algorithms. A list of ranked critical input parameters was obtained after statistical quantification of cumulative distribution functions of input parameters. We found that barite formation is most sensitive to the rate of sulfate ion generation, which is determined by the pyrite dissolution rate coefficient and oxidant availability. In addition, barite formation is sensitive to the initial amounts of barite in HFF and shale, followed by barite thermodynamics/kinetics. For Fe(OH) 3 formation, the ranked factors are Fe(OH)3 precipitation rate coefficients, initial HFF pH, initial Fe(OH) 3 amount in HFF, and oxidant availability. Overall, our results provide insights into managing mineral scale formation during hydraulic fracturing to enhance production. Meanwhile, this study serves as an example of global sensitivity analysis of RTMs using the efficient, straightforward, and open-source DGSA method.

58 GEOSCIENCES↗

DOE OSTI2024

Illuminating early-stage dynamics of heavy-ion collisions through photons at RHIC BES energies

Heavy-ion collisions at $\sqrt {^SNN}$~10 GeV probe the QCD phase diagram at large baryon densities. Because the longitudinal Lorentz contraction is small at these collision energies, understanding the dynamics during the early phase of the collision is essential for the subsequent modeling of the system evolution and for constraining the QGP transport properties at finite baryon densities. Direct photons provide undistorted information on early-stage dynamics. We model relativistic heavy-ion collisions at RHIC Beam Energy Scan energies with a hybrid dynamical approach consisting of a 3D-Glauber initial state followed by viscous hydrodynamics and hadronic transport (MUSIC + UrQMD). The implemented thermal photon emission takes into account the enhancement from finite baryon chemical potentials. We show that direct photon spectra and their anisotropic flow coefficients have a strong sensitivity to the early stage of heavy-ion collisions. Thus, they provide constraints on QGP dynamics complementary to those obtained from hadronic observables.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Modeling Neutron Star Matter in the Age of Multimessenger Astrophysics

Abstract The interpretation of available and forthcoming data obtained from multimessenger astrophysical observations—potentially providing unprecedented access to neutron star properties—will require the development of novel, accurate theoretical models of dense matter. Of great importance, in this context, will be the capability for devising a description of thermal effects applicable to the study of quantities other than the equation of state, such as the transport coefficients and the neutrino mean free path in the nuclear medium. A formalism based on correlated basis states and the cluster expansion technique has been previously employed to derive a well-behaved effective interaction—suitable for use in standard perturbation theory—from a state-of-the-art nuclear Hamiltonian, including phenomenological two- and three-nucleon potentials. Here, we provide a comprehensive and self-contained account of the extension of this approach to the treatment of finite-temperature effects, and report the results of numerical calculations of a number of properties of nuclear matter with arbitrary neutron excess and temperature up to 50 MeV.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

A buoyancy–shear–drag–scalar-based turbulence model for power-law acceleration-driven Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing

A previously developed phenomenological turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz instability-induced mixing based on a general buoyancy–shear–drag model [O. Schilling, “A buoyancy–shear–drag-based turbulence model for Rayleigh–Taylor, reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing,” Physica D 402, 132238 (2020)] is extended to include active or passive scalar mixing and power-law acceleration-driven Rayleigh–Taylor mixing. The buoyancy–shear–drag equations are coupled to a scalar variance equation that is used to define the molecular mixing parameter θ m , and when the scalar is active, modifies the Rayleigh–Taylor and Kelvin–Helmholtz mixing layer growth parameters to depend on the asymptotic value of this parameter, θ mol . Here, the scalar variance equation is closed by algebraically or differentially modeling the scalar variance dissipation rate. Nonlinear analytical solutions of the model are obtained in the total and separate bubble and spike mixing layer width formulations with the algebraic scalar variance dissipation rate for each instability, which are then used to calibrate the mechanical and scalar equation coefficients to predict specific values of physical observables and molecular mixing parameters. Surrogate mechanical and scalar turbulent fields can be constructed by multiplying a presumed self-similar spatial profile by appropriate functions of the width and its time derivative, and of the scalar obtained by solving the ordinary differential model equations. The explicit modeling and solution of turbulent transport equations are not required. The bubble and spike mixing layer width and scalar variance equations are then solved numerically for constant-acceleration Rayleigh–Taylor, impulsively reshocked Richtmyer–Meshkov, and Kelvin–Helmholtz mixing, confirming that the prescribed level of molecular mixing is correctly predicted and illustrating the spatiotemporal evolution of the scalar fields.

Buoyancy–drag↗

DOE OSTI2024

Halogenated Dibenzo[ f , h ]quinoxaline Units Constructed 2D‐Conjugated Guest Acceptors for 19% Efficiency Organic Solar Cells

Abstract Halogenation of Y‐series small‐molecule acceptors (Y‐SMAs) is identified as an effective strategy to optimize photoelectric properties for achieving improved power‐conversion‐efficiencies (PCEs) in binary organic solar cells (OSCs). However, the effect of different halogenation in the 2D‐structured large π‐fused core of guest Y‐SMAs on ternary OSCs has not yet been systematically studied. Herein, four 2D‐conjugated Y‐SMAs (X‐QTP‐4F, including halogen‐free H‐QTP‐4F, chlorinated Cl‐QTP‐4F, brominated Br‐QTP‐4F, and iodinated I‐QTP‐4F) by attaching different halogens into 2D‐conjugation extended dibenzo[ f , h ]quinoxaline core are developed. Among these X‐QTP‐4F, Cl‐QTP‐4F has a higher absorption coefficient, optimized molecular crystallinity and packing, suitable cascade energy levels, and complementary absorption with PM6:L8‐BO host. Moreover, among ternary PM6:L8‐BO:X‐QTP‐4F blends, PM6:L8‐BO:Cl‐QTP‐4F obtains a more uniform and size‐suitable fibrillary network morphology, improved molecular crystallinity and packing, as well as optimized vertical phase distribution, thus boosting charge generation, transport, extraction, and suppressing energy loss of OSCs. Consequently, the PM6:L8‐BO:Cl‐QTP‐4F‐based OSCs achieve a 19.0% efficiency, which is among the state‐of‐the‐art OSCs based on 2D‐conjugated Y‐SMAs and superior to these devices based on PM6:L8‐BO host (17.70%) and with guests of H‐QTP‐4F (18.23%), Br‐QTP‐4F (18.39%), and I‐QTP‐4F (17.62%). The work indicates that halogenation in 2D‐structured dibenzo[ f , h ]quinoxaline core of Y‐SMAs guests is a promising strategy to gain efficient ternary OSCs.

14 SOLAR ENERGY↗

DOE OSTI2020

Interpretation of Lyman opacity measurements in JET with the ITER-like wall using a particle balance approach

Recent measurements of the Ly β /D α ratio in the JET ITER-like wall (ILW) divertor show comparatively more pronounced Ly β reabsorption relative to previous results in JET with the carbon wall (JET-C). At the outer horizontal target ion current rollover point a Ly β reabsorption rate of 60% is measured, increasing to 80% in more pronounced detachment, implying a Ly α reabsorption rate of 90%–98%. The radially resolved Ly β /D α measurements are used to constrain Lyman opacity corrections to atomic rate coefficients using the population escape factor technique in order to capture the local changes to the excited state population structure in the high Lyman opacity regions at and outboard of the outer strike point. To check the self-consistency of the Lyman opacity measurements, a detailed spectroscopic interpretation of the outer divertor particle balance is presented, in which the impact of opacity corrections to the Ly α inverse photon efficiency coefficients is assessed. A five-fold deficit in the estimated D + source rate obtained with optically thin plasma assumptions is reconciled once the opacity corrections are factored into the Ly α photon rate to ionization rate conversion. The experiment results are reproduced in EDGE2D-EIRENE density scan simulations in which an imposed ad hoc Lyman reabsorption rate is recovered spectroscopically using synthetic measurements, and a similar shortfall in the D + source rate estimates is reconciled using Lyman opacity corrected atomic data. The model limitations prevent a more detailed self-consistent analysis of the Lyman opacity impact on divertor parameters and detachment evolution, and hence motivate renewed efforts to re-establish routine exploitation of the photon transport modelling capabilities in the EIRENE code package. Lyman opacity corrections to atomic data coefficients are necessary for spectroscopic interpretation of the JET-ILW divertor plasma, with significant influence on the divertor plasma also a likely consequence of the presence of strong Ly α reabsorption. Finally, a scan of auxiliary heating powers spanning L-mode and H-mode conditions reveals a strong correlation of the measured Ly β opacity with the outer target temperature, suggesting new possibilities for using Lyman opacity measurements in establishing detachment scalings.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2025

VARI3D & PERSENT: Perturbation and Sensitivity Analysis

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry based transport code. This manuscript serves as a single manual for two separate codes: VARI3D and PERSENT. The VARI3D code (VARIational 3D) is based upon the classic finite difference diffusion theory solver available in DIF3D. The PERSENT code (PERturbation and SENitivity for Transport) is based upon the variational nodal method employed in DIF3D termed VARIANT. The VARIANT solver was added to DIF3D in 1995 and has seen continued development and use for the last 18 years. Because VARI3D primarily uses deprecated coding practices, rather than incorporating the perturbation and sensitivity treatments for transport within VARI3D, a new coding development was built using modern Fortran coding. The primary purpose of this manual is to describe the theory behind PERSENT (and by convenience, that of VARI3D) and discuss the input and output of PERSENT along with giving potential users an idea of how to use it. While this manuscript does describe the input and output of VARI3D, the PERSENT code is intended to be the replacement capability of VARI3D as PERSENT can generate nearly identical (if not superior) diffusion theory results. In this manuscript, the relevant aspects of generalized perturbation theory and exact perturbation theory that apply to both VARI3D and PERSENT are covered. The input and output of VARI3D is displayed by excerpting several of the example problems. Similarly, the input and output of PERSENT is displayed along with tips on how best to use the code. Note that the input and output of the inhomogeneous solver wrapped around DIF3D (DIF3D_IFS) is also discussed as it is needed to carry out some of the sensitivities in PERSENT such as reaction rate ratios. This manuscript describes several perturbation and sensitivity problems, and the results computed using PERSENT. From these sections, potential users should find that PERSENT provides not only the typical tables of numbers desired in perturbation and sensitivity analysis work, but also can visually plot the result for a more thorough understanding of the space and energy distribution (Section 5). Overall, PERSENT is observed to produce accurate reactivity worths and sensitivities for the displayed set of test problems and clearly demonstrates the need to have a transport-based sensitivity capability as evident from the thousands of percent errors observed in the 21-group hexagonal fast reactor problem (covered in Section 7). The uncertainty calculation capability is described in Section 3 and demonstrated in Section 7.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2021

Using Computationally-Determined Properties for Machine Learning Prediction of Self-Diffusion Coefficients in Pure Liquids

The ability to predict transport properties of liquids quickly and accurately will greatly improve our understanding of fluid properties both in bulk and complex mixtures, as well as in confined environments. Such information could then be used in the design of materials and processes for applications ranging from energy production and storage to manufacturing processes. As a first step, we consider the use of machine learning (ML) methods to predict the diffusion properties of pure liquids. Recent results have shown that Artificial Neural Networks (ANNs) can effectively predict the diffusion of pure compounds based on the use of experimental properties as the model inputs. In the current study, a similar ANN approach is applied to modeling diffusion of pure liquids using fluid properties obtained exclusively from molecular simulations. A diverse set of 102 pure liquids is considered, ranging from small polar molecules (e.g., water) to large nonpolar molecules (e.g., octane). Self-diffusion coefficients were obtained from classical molecular dynamics (MD) simulations. Since nearly all the molecules are organic compounds, a general set of force field parameters for organic molecules was used. The MD methods are validated by comparing physical and thermodynamic properties with experiment. Computational input features for the ANN include physical properties obtained from the MD simulations as well as molecular properties from quantum calculations of individual molecules. Furthermore, fluid properties describing the local liquid structure were obtained from center of mass radial distribution functions (COM-RDFs). Feature sensitivity analysis revealed that isothermal compressibility, heat of vaporization, and the thermal expansion coefficient were the most impactful properties used as input for the ANN model to predict the MD simulated self-diffusion coefficients. The MD-based ANN successfully predicts the MD self-diffusion coefficients with only a subset (2 to 3) of the available computationally determined input features required. A separate ANN model was developed using literature experimental self-diffusion coefficients as model targets. Although this second ML model was not as successful due to a limited number of data points, a good correlation is still observed between experimental and ML predicted self-diffusion coefficients.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Modeling impurity concentrations in liquid argon detectors

Impurities in noble liquid detectors used for neutrino and dark matter experiments can significantly impact the quality of data. We present an experimentally verified model for describing the dynamics of impurity distributions in liquid argon (LAr) detectors. The model considers sources, sinks, and transport of impurities within and between the gas and liquid argon phases. Measurements of oxygen concentrations in a 20-L LAr multi-purpose test stand are compared to calculations made with this model to show that an accurate description of the concentrations under various operational conditions can be obtained. A result of this analysis is a determination of Henry’s coefficient for oxygen in LAr. Additionally, these calculations also show that some processes have small effects on the impurity dynamics and excluding them yields a solution as a sum of two exponential terms. This solution provides a simple way to extract Henry’s coefficient with negligible approximation error. It is applied to the data and the Henry’s coefficient for oxygen in LAr is obtained as 0.84 +0.09 -0.05 , consistent with literature results. Finally, based on the analysis of the data with the model, we further suggest that, for a large liquid argon detector, barriers to flow (“baffles”) installed in the gas phase to restrict flow can help reduce the ultimate impurity concentration in the LAr.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

Hall effect in gated single-wall carbon nanotube films

The presence of hopping carriers and grain boundaries can sometimes lead to anomalous carrier types and density overestimation in Hall-effect measurements. Previous Hall-effect studies on carbon nanotube films reported unreasonably large carrier densities without independent assessments of the carrier types and densities. Here, we have systematically investigated the validity of Hall-effect results for a series of metallic, semiconducting, and metal–semiconductor-mixed single-wall carbon nanotube films. With carrier densities controlled through applied gate voltages, we were able to observe the Hall effect both in the n- and p-type regions, detecting opposite signs in the Hall coefficient. By comparing the obtained carrier types and densities against values derived from simultaneous field-effect-transistor measurements, we found that, while the Hall carrier types were always correct, the Hall carrier densities were overestimated by up to four orders of magnitude. This significant overestimation indicates that thin films of one-dimensional SWCNTs are quite different from conventional hopping transport systems.

42 ENGINEERING↗

DOE OSTI2024

VARI3D & PERSENT: Perturbation and Sensitivity Analysis (Revision 5)

The nodal diffusion method is one of the most widely used approaches in modern reactor analysis. In the nodal diffusion method, a coarse multi-group set of “homogenized” parameters is constructed such that the complex geometry of a reactor core along with the energy dependence of neutron and gamma ray cross sections in a nuclear reactor are conserved in the simpler geometry. The homogenization is typically done on a fuel assembly level as is the case in the DIF3D code developed at Argonne National Laboratory. The nodal methodology is used primarily to predict fuel cycle behavior of nuclear systems of which there is a substantial amount of validation in the literature. Another use of the nodal method is to obtain reactivity coefficients and kinetics parameters for use in a safety analysis of a given nuclear reactor. While there are many ways to obtain reactivity worth and kinetics parameters, the work presented in this manuscript is unique as it provides the user with the ability to compute reactivity worths, kinetics parameters, and cross section sensitivities with a Cartesian and hexagonal geometry-based transport code.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Viscous Properties of a Degenerate One-Dimensional Fermi Gas

In this paper, we study the viscous properties of a system of weakly interacting spin-1/2 fermions in one dimension. Accounting for the effect of interactions on the quasiparticle energy spectrum, we obtain the bulk viscosity of this system at low temperatures. Our result is valid for frequencies that are small compared with the rate of fermion backscattering. For frequencies larger than this exponentially small rate, the excitations of the system become decoupled from the center of mass motion, and the fluid is described by two-fluid hydrodynamics. We calculate the three transport coefficients required to describe viscous dissipation in this regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗