SEARCH · Engineering Papers
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.
The Orthogonal Diffusion Discretizations (ODD) library is a collection of low order (Diffusion, P1, Ray Tracing, and Variable Eddington Factor (VEF)) thermal radiation solvers for orthogonal structured grids driven by an application programing interface (API).
This is an hourly future weather dataset for energy modeling applications. The dataset is primarily based on the output of a regional climate model (RCM), i.e., the Weather Research and Forecasting (WRF) model version 3.3.1. The WRF simulations are driven by the output of a general circulation model (GCM), i.e., the Community Climate System Model version 4 (CCSM4). This dataset is in the EPW format, which can be read or translated by more than 25 building energy modeling programs (e.g., EnergyPlus, ESP-r, and IESVE), energy system modeling programs (e.g., System Advisor Model (SAM)), indoor air quality analysis programs (e.g., CONTAM), and hygrothermal analysis programs (e.g., WUFI). It contains 13 weather variables, which are the Dry-Bulb Temperature, Dew Point Temperature, Relative Humidity, Atmospheric Pressure, Horizontal Infrared Radiation Intensity from Sky, Global Horizontal Irradiation, Direct Normal Irradiation, Diffuse Horizontal Irradiation, Wind Speed, Wind Direction, Sky Cover, Albedo, and Liquid Precipitation Depth. The weather data is created for two emissions scenarios: RCP4.5 and RCP8.5 and spans two 10-year time slices in the future: 2045 - 2054 and 2085 - 2094. It offers a spatial resolution of 12 km by 12 km with extensive coverage across most of North America. Due to the enormous size of the entire dataset, in the first stage of its distribution, we provide 20 years of future weather data for the centroid of each Public Use Microdata Area (PUMA), excluding Hawaii. PUMAs are non-overlapping, statistical geographic areas that partition each state or equivalent entity into geographic areas containing no fewer than 100,000 people each. The 2,378 PUMAs as a whole cover the entirety of the U.S. The weather data can be utilized alongside the large-scale energy analysis tools, ResStock and ComStock, developed by National Renewable Energy Laboratory, whose smallest resolution is at the PUMA scale. The data for RCP4.5 is still being processed and will be published soon.
We report significant advances in understanding and reducing nonradiative Shockley-Read-Hall recombination in polycrystalline CdSe x Te 1-x , leading to microsecond charge carrier lifetimes. In undoped Al 2 O 3 -passivated heterostructures we find external radiative efficiency 0.2%, quasi-Fermi level splitting 950 mV, mobility 100 cm 2 /(Vs), and diffusion length 14 µm. In solar cells measured lifetimes can exceed 1 µs. We interpret this data to indicate MgZnO/CdSeTe interface recombination velocity <; 100 cm/s. Based on our results, it appears CdTe PV technology has potentially overcome longstanding “recombination lifetime” limitation and in the near future will transition to improving other aspects of device design.
The discrete ordinates discontinuous Galerkin (SN-DG) method is a well-established and practical approach for solving the radiative transport equation. In this paper, we study a low-memory variation of the upwind SN-DG method. The proposed method uses a smaller finite element space that is constructed by coupling spatial unknowns across collocation angles, thereby yielding an approximation with fewer degrees of freedom than the standard method. Like the original SN-DG method, the low-memory variation still preserves the asymptotic diffusion limit and maintains the characteristic structure needed for mesh sweeping algorithms. While we observe second-order convergence in the scattering dominated, diffusive regime, the low-memory method is in general only first-order accurate. To address this issue, we use upwind reconstruction to recover second-order accuracy. Finally, for both methods, numerical procedures based on upwind sweeps are proposed to reduce the system dimension in the underlying Krylov solver strategy.
Plasma nitridation was conducted to modify the surfaces of Zircaloy-4. Scanning electron microscopy (SEM), transmission electron microscopy (TEM), and Raman analysis were used to characterize microstructures and phases. Surface indentation and cross-sectional indentation were performed to evaluate mechanical property changes. Nitridation forms a thin layer of ZrN phase, followed by a much deeper layer affected by nitrogen diffusion. The ZrN phase is confirmed by both TEM and Raman characterization. The Raman peaks of ZrN phase show a temperature dependence. The intensity increases with increasing nitridation temperatures, reaches a maximum at 700 °C, and then decreases at higher temperatures. The ZrN layer appears as continuous small columnar grains. The surface polycrystalline ZrN phase is harder than the bulk by a factor of ~8, and the nitrogen diffusion layer is harder by a factor of ~2–5. The activation energy of nitrogen diffusion was measured to be 2.88 eV. The thickness of the nitrogen-hardened layer is controllable by changing the nitridation temperature and duration.
Ni-Cr alloys are widely employed as structural materials in fast nuclear reactors but are vulnerable to high-temperature helium (He) embrittlement (HTHE) under fast neutron irradiation. A comprehensive understanding of He diffusion in Ni-Cr alloys, which governs the kinetics of HTHE, is therefore essential for developing resilient materials and preventing failure. In this work, we reveal the underlying mechanisms of He diffusion in pure Ni and Ni-Cr alloys by integrating density functional theory (DFT) with atomic kinetic Monte Carlo (AKMC) simulations. Our findings uncover a non-monotonic dependence of He diffusivity on Cr concentration, contradicting the monotonic trends predicted by DFT-parameterized theories. At low Cr concentrations, He diffusion is dominated by short-range trapping, characterized by multiple trapping sites and a distinct mechanism within the first nearest neighbor of Cr, differing from that in pure Ni. At high Cr concentrations, these local traps become interconnected, forming long-range fast diffusion channels that enhance He mobility. The competition between localized trapping and extended channeling results in a diffusivity that first decreases, then increases with rising Cr content. These atomic-scale insights offer critical guidance for the design of radiation-tolerant Ni-based alloys. Moreover, the combined DFT-AKMC methodology and the concept of random walker diffusion through interconnected energy basins present a broadly applicable framework for studying transport phenomena in disordered systems.
Starting in 2019, the Lagrangian Applications Project (LAP) and the Transport Project set out to develop a new ALE/Lagrangian radiation-hydrodynamics (RH) capability in a new code product named LUMOS. This work was done under the guidance of the Advanced Simulation and Computing (ASC) program with the goal of producing software capable of leveraging the high-order thermal radiative transfer (TRT) solvers provided by the Jayenne and Capsaicin software projects. This capability supplements the existing gray diffusion solver that is currently available in LAP’s FLAG code [1, 2, 3] for RH. The remainder of this memo describes the coupling of the radiation and hydrodynamics solvers within LAP’s LUMOS code. Initially delivered in 2020 as part of an L2 milestone [4], this capability continues to mature in FY21 with more efficient and robust algorithms, new support for ALE, and support for mixed materials per cell. Today, LUMOS is provided as a standard end-user product in the suite of LAP tools provided in each release cycle. Code access can be requested at https://asc.lanl.gov.
The study, quantification, and monitoring of iodine in nuclear power generation systems have been paramount due to its environmental and health implications. Thorough speciation and transport studies strengthened the development and deployment of water reactors, which are now required for the new generation of Molten Salt Reactors (MSR), where the fuel and coolant are a complex liquid system of suitable halides. The speciation, radiation-induced chemistry, and transport of iodine (I-127, I-129 and I-131) under extreme conditions of temperature and ionizing radiation are among the challenges for the future implementation of MSR’s. In aqueous systems, the triiodide-iodine-iodide (I3--I2-I-) equilibrium is strongly dependent on the conditions of the medium. In molten iodides and chlorides, the triiodide detection (and its disproportionation) poses challenges due to the extreme conditions of temperature and radiation, but studies have shown a dependence on temperature, absorbed dose and salt composition. Through absorption spectroscopy techniques (diffuse reflectance and high-temperature UV-Vis), as well as resonance (Electron Paramagnetic Resonance) it is possible to identify primary transients induced by ionizing radiation trapped in the solid phase (Solid-state radiolysis), and steady-state main absorption components with a fast-acquisition real-time spectroscopy. The presence of heterogeneous phase (colloids), viscosity changes, and gas evolution are among the observed processes. A specific case of study is speciation in the presence of Ni. A corrosion product with magnetic properties due to unpaired electrons in its electronic structure. Ferromagnetic materials are monitored in solidified salts to identify aggregation and phase separation. Comparative studies of eutectic mixtures of LiI-KI and LiCl-KCl with known addition of dopants will be presented, with a perspective of radiation chemistry and its possible implications in complex systems.
Herein, to understand the microstructural stability of candidate plasma-facing materials under fusion-relevant environments, neutron irradiation of W and W-3%Re alloys with and without K and La dopants was performed in the mixed-spectrum High Flux Isotope Reactor at nominal temperatures of ~850 °C and ~1100 °C to calculated doses between 0.42 and 0.47 dpa. To the best of our knowledge, this study presents the first experimental evidence of radiation-enhanced recrystallization in W and undoped W–Re alloys at ~850 °C, conditions where thermal annealing does not cause any grain growth in a similar timescale. Potassium- or lanthanum-doped tungsten alloys showed more resistance to radiation-enhanced grain growth. We explain the acceleration of grain growth by analyzing the self-diffusion constant under atomic displacement environments. The microstructural observations of the studied W variants suggest that La doping is more effective than K doping for mitigating recrystallization. This study also found that radiation-enhanced recrystallization is an important consideration when designing and applying W to plasma-facing components in future nuclear fusion reactors.
We have irradiated high-density polyethylene ribbons and cylinders to ~200 kGy at a dose rate of ~5 kGy day -1 , using a 60 Co source. The radiolytically-generated gas consists almost entirely of H 2 . Based on the time- and dose-dependence of the outgassing, we have determined the radiation- chemical yield for the evolved H 2 , as well as the diffusion coefficient of H 2 in HDPE. We find that the evolved gas yield G evolved depends on the interplay of sample geometry, experimental timescale, and the diffusion coefficient of H 2 in HDPE. When these factors facilitate completion of diffusion prior to gas quantification, such as for sub-mm-thickness ribbons, we measure a time- and dose-invariant value of G evolved = 0.0076 cm 3 std. H 2 J -1 that is in agreement with literature reports. When sample geometry and experiment timescale do not permit completion of diffusion, we measure a time- and dose-dependent value of G evolved . We have developed a model for the transient outgassing behavior of the cylinders, and fitted it to the experimental data. Based on these fits, we estimate a room-temperature diffusivity of D H 2 /HDPE ≈ 1.67 ± 0.12 x 10 -6 cm 2 s -1 for H 2 in HDPE, in good agreement with literature values. At present, the effects of dose rate and γ-photon energy on HDPE gas yield are not known and require further investigation.
We have developed a first-passage kinetic Monte Carlo approach for materials aging to investigate the sensitivity of void swelling to model parameters, including helium bubble density and size distribution. In addition to explicitly accounting for the spatial distribution of individual point defects, bubbles, and voids, our approach can simulate total doses equivalent to 100 years of natural aging on statistically representative volumes of materials. This technique enables us to study the effects on swelling and radiation damage evolution due to temperature and dose rate (as altered in artificially aged experiments), differences in effective interaction radii between vacancies and interstitials, and varying defect diffusion activation energies, while providing more detailed information than previous rate-equation based approaches. In conclusion, our results indicate that spatial effects that are not modeled in mean-field rate theories could play a significant role in void swelling initiation and growth for certain regimes of model parameters.
Vacancies are missing atoms in a crystalline material, and occur both at equilibrium (varying with temperature) and out of equilibrium such as when crystalline materials are damaged with radiation or corrosion. While we know of their importance, particularly regarding diffusive mechanisms, it is not straightforward to experimentally measure their concentration or to visualize them directly. Traditionally, measurements such as positron annihilation spectroscopy and to some degree X-ray diffraction can measure average concentrations, but generally lack the ability to visualize or quantify a heterogenous concentration of vacancies that can occur at the level of individual defects and microstructural features. Here, we present a method to map vacancy concentrations and their distribution using local lattice parameter measurements with a high-resolution electron microscope. Our method utilizes a Au thin film as a model to demonstrate the method via four-dimensional scanning transmission electron microscopy (4D-STEM) by correlating the differences between changes in lattice parameter and the volumetric thermal expansion during in situ heating experiments. The vacancy mapping methodology is also applied to non-equilibrium defects accumulated in pure Al via knock-on electron beam irradiation. Our method demonstrates the ability to map point defect concentrations in heterogeneous systems in situ with nanometer spatial resolution. The result is a technique that can provide direct measurements of vacancy concentrations at the level of individual defects in studies of materials in and out of equilibrium.
We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.
Atomic defects in the solid state are a key component of quantum repeater networks for long-distance quantum communication1. Recently, there has been significant interest in rare earth ions, in particular Er 3+ for its telecom band optical transition that allows long-distance transmission in optical fibres. However, the development of repeater nodes based on rare earth ions has been hampered by optical spectral diffusion, precluding indistinguishable single-photon generation. Here, we implant Er 3+ into CaWO 4 , a material that combines a non-polar site symmetry, low decoherence from nuclear spins and is free of background rare earth ions, to realize significantly reduced optical spectral diffusion. For shallow implanted ions coupled to nanophotonic cavities with large Purcell factor, we observe single-scan optical linewidths of 150 kHz and long-term spectral diffusion of 63 kHz, both close to the Purcell-enhanced radiative linewidth of 21 kHz. This enables the observation of Hong–Ou–Mandel interference between successively emitted photons with a visibility of V = 80(4)%, measured after a 36 km delay line. We also observe spin relaxation times T 1,s = 3.7 s and T 2,s > 200 μs, with the latter limited by paramagnetic impurities in the crystal instead of nuclear spins. Furthermore, this represents a notable step towards the construction of telecom band quantum repeater networks with single Er 3+ ions.
We numerically investigate the internal evolution of multiphase clouds, which are at rest with respect to an ambient, highly ionized medium (HIM) representing the hot component of the circumgalactic medium. Time-dependent saturated thermal conduction and its implications like condensation rates and mixing efficiency are assessed in multiphase clouds. Our simulations are carried out by using the adaptive mesh refinement code FLASH. The model clouds are initially in both hydrostatic and thermal equilibrium and are in pressure balance with the HIM. Thus, they have steep gradients in both temperature and density at the interface to HIM leading to non-negligible thermal conduction. Several physical processes are considered numerically or semi-analytically: thermal conduction, radiative cooling and external heating of gas, self-gravity, mass diffusion, and dissociation of molecules and ionization of atoms. It turns out that saturated thermal conduction triggers a continuous condensation irrespective of cloud mass. Dynamical interactions with ambient HIM all relate to the radial density gradient in the clouds: (1) mass flux due to condensation is the higher the more homogeneous the clouds are; (2) mixing of condensed gas with cloud gas is easier in low-mass clouds, because of their shallower radial density gradient; and thus (3) accreted gas is distributed more efficiently. A distinct and sub-structured transition zone forms at the interface between cloud and HIM, which starts at smaller radii and is much narrower as deduced from analytical theory.
ABSTRACT We present magnetohydrodynamic (MHD) simulations of the star-forming multiphase interstellar medium (ISM) in stratified galactic patches with gas surface densities Σgas = 10, 30, 50, and 100 $\mathrm{M_\odot \, pc^{-2}}$. The silcc project simulation framework accounts for non-equilibrium thermal and chemical processes in the warm and cold ISM. The sink-based star formation and feedback model includes stellar winds, hydrogen-ionizing UV radiation, core-collapse supernovae, and cosmic ray (CR) injection and diffusion. The simulations follow the observed relation between Σgas and the star formation rate surface density ΣSFR. CRs qualitatively change the outflow phase structure. Without CRs, the outflows transition from a two-phase (warm and hot at 1 kpc) to a single-phase (hot at 2 kpc) structure. With CRs, the outflow always has three phases (cold, warm, and hot), dominated in mass by the warm phase. The impact of CRs on mass loading decreases for higher Σgas and the mass loading factors of the CR-supported outflows are of order unity independent of ΣSFR. Similar to observations, vertical velocity dispersions of the warm ionized medium (WIM) and the cold neutral medium (CNM) correlate with the star formation rate as $\sigma _\mathrm{z} \propto \Sigma _\mathrm{SFR}^a$, with a ∼ 0.20. In the absence of stellar feedback, we find no correlation. The velocity dispersion of the WIM is a factor ∼2.2 higher than that of the CNM, in agreement with local observations. For $\Sigma _\mathrm{SFR} \gtrsim 1.5 \times 10^{-2}\, \mathrm{M}_\odot \, \mathrm{yr}^{-1}\, \mathrm{kpc}^{-2}$ the WIM motions become supersonic.
Group-V doping can enhance the efficiency of CdTe photovoltaic devices by increasing the open-circuit voltage (Voc), if long lifetimes can be retained and radiative voltage loss can be avoided. In this preliminary study, ex-situ Sb doping via drive-in diffusion improved grain sizes compared to Sb-free films without CdCl2 treatment. Sb doping did not reduce the carrier lifetime for uniform CdSeTe absorbers. Additionally, the CST-air surface recombination decreased with Sb. No potential fluctuations were observed by PL at this doping level. The [Sb] and carrier concentration for this study are not yet measured so the conclusions drawn are all preliminary
Light water is the most important neutron moderator in many reactor applications. Thermal neutron scattering kernels for H bound in H 2 O are represented by ENDF-format thermal scattering laws (TSLs) evaluated at discrete temperatures T . Nuclear data evaluations are conventionally validated utilizing a combination of critical benchmarks and experimental cross section data. Existing public critical benchmarks may be inadequate for testing water TSLs over the full range of T of interest in reactor applications, and experimental thermal scattering cross section data for water is sparse at elevated T . In this work, MC21 is used to simulate the decay of the equilibrium thermal neutron flux in light-water spheres of several radii. The fundamental-mode time eigenvalue $\mathcal{a}$ is calculated at 22 °C and 227 °C (at saturation pressure) for each sphere using three different H-H 2 O TSLs. Fitting $\mathcal{a}$ to a polynomial function of geometric buckling allows calculation of the thermal neutron diffusion length L . Extrapolation lengths are treated as a function of the transport mean free path and geometry. Experimental L data from 25 publications at 49 temperatures (from 10 °C to 295 °C), computed by several different time-dependent and space-dependent decay methods, is used to develop an empirical fit of L vs. T . The MC21-calculated L for the TSLs tested are within ±1% of the predicted values at 22 °C and 227 °C. This validation approach, which may be repeated at arbitrary T , constitutes an integral benchmark specific to and characterizing the detailed physics of the thermal scattering kernel applied.