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.
Heat transfer behavior in sCO2 power cycles when firing solid fuels
As an asteroid descends toward Earth, it deposits energy in the atmosphere through aerodynamic drag and ablation. Asteroid impact risk assessments rely on energy deposition estimates to predict blast overpressures and ground damage that may result from an airburst, such as the one that occurred over Chelyabinsk, Russia in 2013. The rates and altitudes at which energy is deposited along the entry trajectory depend upon how the bolide fragments, which in turn depends upon its internal structure and composition. In this work, we have developed an analytic asteroid fragmentation model to assess the atmospheric energy deposition of asteroids with a range of structures and compositions. The modeling approach combines successive fragmentation of larger independent pieces with aggregate debris clouds released with each fragmentation event. The model can vary the number and masses of fragments produced, the amount of mass released as debris clouds, the size-strength scaling used to increase the robustness of smaller fragments, and other parameters. The initial asteroid body can be seeded with a distribution of independent fragment sizes amid a remaining debris mass to represent loose rubble pile conglomerations, can be given an outer regolith later, or can be defined as a coherent or fractured monolith. This approach enables the model to represent a range of breakup behaviors and reproduce detailed energy deposition features such as multiple flares due to successive burst events, high-altitude regolith blow-off, or initial disruption of rubble piles followed by more energetic breakup of the constituent boulders. These capabilities provide a means to investigate sensitivities of ground damage to potential variations in asteroid structure.
Previous Monte Carlo simulations provide a data base for properties of secondary electron emission (SEE) from insulators and metals. Incident primary electrons are considered at energies up to 1200 eV. The behavior of secondary electrons is characterized by (1) yield vs. primary energy E(sub p), (2) distribution vs. secondary energy E(sub s), and (3) distribution vs. angle of emission theta. Special attention is paid to the low energy range E(sub p) up to 50 eV, where the number and energy of secondary electrons is limited by the finite band gap of the insulator. For primary energies above 50 eV the SEE yield curve can be conveniently parameterized by a Haffner formula. The energy distribution of secondary electrons is described by an empirical formula with average energy about 8.0 eV. The angular distribution of secondaries is slightly more peaked in the forward direction than the customary cos theta distribution. Empirical formulas and parameters are given for all yield and distribution curves. Procedures and algorithms are described for using these results to find the SEE yield, and then to choose the energy and angle of emergence of each secondary electron. These procedures can readily be incorporated into numerical simulations of plasma-solid surface interactions in low earth orbit.
Abstract Electrochemical biomass oxidation coupled with hydrogen evolution offers a promising route to generate value‐added chemicals and clean energy. The complex adsorption behavior of 5‐hydroxymethyl furfural (HMF) and hydroxyl ions (OH − ) on the electrocatalyst surface during HMF electrooxidation reaction (HMFOR) necessitates an in‐depth understanding of active sites available for adsorption. Herein, oxygen vacancy (V O ) defects are introduced in NiFe layered double hydroxide (LDH) using Ce dopants to manipulate electronic structure. Synchrotron‐based HE‐XRD and XAS indicate negligible V O in La‐doped NiFe while Ce doping leads to V O defects due to flexible Ce redox (Ce 3+ ↔ Ce 4+ ). The V O ‐rich Ce‐NiFe exhibits higher Faradic efficiency of ≈90% to produce 2,5‐furan dicarboxylic acid (FDCA), far greater than ≈60% for NiFe V O in Ce‐NiFe act as alternative active sites for OH − adsorption, hence reducing adsorption competition for the same metal sites. DFT calculation results corroborate experimental findings by showcasing that the presence of V O in Ce‐NiFe manipulates the adsorption energies and facilitates the chemical adsorption OH − in V O to improve HMFOR. In situ HE‐XRD derived pair distribution function coupled to RMC simulations confirm OH − trapping in V O and HMF adsorption on metal centers as evident by interlayer distance evolution. Taken together, this work showcases routes for dual‐site electrocatalyst design for improved biomass electrooxidation.
Topological Weyl semimetal (WSM) is a solid-state realization of chiral Weyl fermions, whose phonon behaviors provide in-depth knowledge of their electronic properties. In this work, anisotropic Fano resonance is observed in a type-II WSM candidate LaAlSi by polarized Raman spectroscopy. The asymmetric line shape occurs for the $B^2_1$ phonon mode of LaAlSi only for 488- and 532-nm laser excitations but not for 364-, 633-, and 785-nm excitations, suggesting the excitation selectivity. The asymmetry, frequency, and linewidth of the $B^2_1$ phonon mode, along with the spectral background, all show fourfold rotational symmetry as a function of the polarization angle in the polarized Raman spectra. While the shift of Raman frequency in a metal or semimetal is typically attributed to Kohn anomaly, here we show that the anisotropic frequency shift in LaAlSi cannot be explained by the effect of Kohn anomaly, but potentially by the anisotropic scattering background of Fano resonance. Origins of the excitation-energy dependence and anisotropic behavior of the Fano resonance are discussed by the first-principles calculated electronic band structure and phonon dispersion.
Covalent functionalization of single-walled carbon nanotubes (SWCNTs) with organic molecules results in red-shifted emissive states associated with sp 3 -defects in the tube lattice, which facilitate their improved optical functionality, including single-photon emission. The energy of the defect-based electronic excitations (excitons) depends on the molecular adducts, the configuration of the defect, and concentration of defects. In this work, we model the interactions between two sp 3 -defects placed at various distances in the (6,5) SWCNT using time-dependent density functional theory. Calculations reveal that these interactions conform to the effective model of J-aggregates for well-spaced defects (>2 nm), leading to a red-shifted and optically allowed (bright) lowest energy exciton. H-aggregate behavior is not observed for any defect orientations, which is beneficial for emission. The splitting between the lowest energy bright and optically forbidden (dark) excitons and the pristine excitonic band are controlled by the single-defect configurations and their axial separation. These findings enable a synthetic design strategy for SWCNTs with tunable near-infrared emission.
Laser powder directed energy deposition was employed to fabricate Haynes® 282 Ni-based superalloy for high temperature creep testing. Material was heat treated using a single-step aging heat treatment of 800 °C for 4 h. Microstructural analysis was conducted in the as-deposited, solution-annealed, and age-hardened conditions. γ′ precipitates were absent in the as-deposited and annealed condition but formed during aging with an average size of 48.5 ± 12 nm. γ’ denuded zones ranging in thickness from 73 to 469 nm were observed adjacent to an effectively continuous grain boundary carbide film that covered 72.5% of the measured boundary length. The heat-treated material contained a high-volume fraction of TiC carbides ranging in size from 38 to 457 nm. Electron backscatter diffraction showed partial recrystallization, with 77-83% recrystallized area and a heterogeneous grain-size distribution. Creep testing at 750 °C at 300, 320, and 350 MPa revealed significantly inferior performance compared to wrought material, with higher strain rates, earlier tertiary onset, and shorter rupture lives. Fractography shows that failure was dominated by intergranular rupture. Post rupture microstructural analysis results suggest that the creep behavior was primarily governed by the non-equilibrium interface network rather than the bulk γ/γ’ matrix. In conclusion, the findings indicate that achieving wrought-like creep performance necessitates substantial improvements in the additive process, specifically focusing on enhanced homogenization, full recrystallization, and precise control over the morphology and distribution of carbides.
Advances in technology have begun to open new opportunities for behavior-based and technical approaches to managing residential energy use and meet sustainability-related objectives. Visions of the future predict homes with smart technologies delivering enhanced comfort and cost savings to residents; utility-partners who can remotely optimize energy resources to meet grid needs; and occupants who play more active roles in the energy system supported by advanced information communication technologies. Each of these scenarios implies augmented control over home energy use, yet uncertainties remain regarding which ones will deliver the greatest grid benefits and services to customers in a given situation. These scenarios also raise broader questions regarding customer agency and the relationship between customers and third parties moving forward. While both the provision of information to spur behavior change and automated technologies theoretically enhance control over energy use in the built environment, these strategies are not often studied from an integrated perspective. Seeking to address this gap and develop a deeper understanding of the evolving paradigm of control over home energy use, this paper presents a taxonomy to evaluate perspectives from public policy (ex. demand-side management), technological innovation (automated controls), and user-agency (ex. the role of behavior change) on approaches to managing home energy use. We draw on theoretical and empirical evidence from across disciplines to detail the dimensions and implications of deploying programs that incorporate various levels of control and anticipate such a taxonomy will help holistically map out and evaluate tradeoffs between different approaches to demand-side management moving forward.
The lack of completeness and systematic studies of cumulative fission product yield (FPY) from the three major actinides of 235 U, 238 U, and 239 Pu, over a broad energy range, was the primary motivation for the LLNL-LANL-TUNL collaboration. Given the focus on resolving FPY behavior around fission spectrum neutron energies, our 2016 data provides FPY measurements in fine energy steps in the range of 0.5–5.5 MeV and an additional anchor point at 14.8 MeV. Recently, our collaboration filled the gap between second chance fission and 14.8 MeV, providing a comprehensive study at eleven incident neutron energies from 0.5 to 14.8 MeV. However, during the course of the last measurements, it was found that the reported masses of the actinide deposits in the fission ionization chambers, which were used to measure the number of fissions, were not consistent with the complementary tests. Therefore, some corrections to the old data are needed. The purpose of this paper is to summarize all of the high-yield FPY data obtained for 235 U, 238 U, and 239 Pu isotopes using quasi-monoenergetic neutron beams at 11 incident energies ranging from E n = 0.5 to 14.8 MeV.
Numerous studies of the ground-state decay of the pygmy dipole resonance (PDR) have been carried out in the past. However, data on the decay of the PDR to low-lying excited states is still very scarce due to limitations of the sensitivity to weak branching transitions of experimental setups. Here, we present a detailed examination of the low-energy dipole response of 128 Te and 130 Te below their neutron separation thresholds of 8.8 and 8.5 MeV, respectively. Photonuclear reactions with the subsequent γ-ray spectroscopy of the decay channel with continuous-energy bremsstrahlung at varying endpoint energies and linearly polarized quasimonochromatic γ-ray beams with energies ranging from 2.7 to 8.9 MeV in steps of roughly 250 keV were used for probing the decay behavior of the low-energy dipole response in 128Te and 130Te. In addition, (γ,γ' γ") reactions were used to study the population of low-lying states of 128 Te. Spin-parity quantum numbers and reduced transition probabilities are determined for individual photo-excited states. The analysis of average decay properties for nuclear levels in narrow excitation-energy bins enable the extraction of photoabsorption cross sections, average branching ratios to the $2$$^{+}_{1}$ state, and the distinction between E1 and M1 transitions to the ground state and to the $2$$^{+}_{1}$ state accounting for resolved and unresolved transitions.
Vacancy and self-interstitial atomic diffusion coefficients in concentrated solid solution alloys can have a non-monotonic concentration dependence. Here, the kinetics of monovacancies and $\langle$100$\rangle$ dumbbell interstitials in Ni–Fe alloys are assessed using lattice kinetic Monte Carlo (kMC). The non-monotonicity is associated with superbasins, which impels using accelerated kMC methods. Detailed implementation prescriptions for first passage time analysis kMC (FPTA-kMC), mean rate method kMC (MRM-kMC), and accelerated superbasin kMC (AS-kMC) are given. The accelerated methods are benchmarked in the context of diffusion coefficient calculations. The benchmarks indicate that MRM-kMC underestimates diffusion coefficients, while AS-kMC overestimates them. In this application, MRM-kMC and AS-kMC are computationally more efficient than the more accurate FPTA-kMC. Our calculations indicate that composition dependence of migration energies is at the origin of the vacancy’s non-monotonic behavior. In contrast, the difference between formation energies of Ni–Ni, Ni–Fe, and Fe–Fe dumbbell interstitials is at the origin of their non-monotonic diffusion behavior. Additionally, the migration barrier crossover composition—based on the situation where Ni or Fe atom jumps have lower energy barrier than the other one—is introduced. KMC simulations indicate that the interplay between composition dependent crossover of migration energy and geometrical site percolation explains the non-monotonic concentration-dependence of atomic diffusion coefficients.
Direct numerical simulations are used to investigate the interaction of isotropic quasi-incompressible turbulence with a weak shock wave. A linear analysis of the interaction is conducted for comparison with the simulations. Both the simulations and the analysis show that turbulence is enhanced during the interaction. Turbulent kinetic energy and transverse vorticity components are amplified, and turbulent lengthscales are decreased. It is suggested that the amplification mechanism is primarily linear. Simulations also showed a rapid evolution of turbulent kinetic energy just downstream of the shock, a behavior not reproduced by the linear analysis. Analysis of the budget of the turbulent kinetic energy transport equation shows that this behavior can be attributed to the pressure transport term. Multiple compression peaks were found along the mean streamlines at locations where the local shock thickness had increased significantly.
We point out some common qualitative features of the Coon amplitude — a family of deformations of the Veneziano amplitude with logarithmic Regge trajectories — and the open string scattering amplitude for strings ending on a D-brane in AdS. Both reduce to the Veneziano amplitude at relatively low energies. Both systems have an accumulation point in their spectrum, with an infinite number of states below a certain energy. The approach to this point is very similar. Both have the same high-energy fixed-angle behavior. Nevertheless, we find some differences in the spectrum of states with highest angular momentum. These similarities suggest that there may exist a string background that realizes the Coon amplitude.
Abstract Grain boundaries in polycrystalline materials significantly affect their properties, such as ionic transport, corrosion, and chemical durability. The pyrochlore compound (Gd 2 Ti 2 O 7 ) is employed as a model for complex oxides and is known for its diverse applications, including nuclear waste immobilization. Density functional theory-based first-principles molecular dynamics simulations were performed at different temperatures on the hydrated grain boundary system. The results show extensive transformations within the grain boundaries among hydrous water species (OH − , H 2 O, and H 3 O + ). The temperature dependence of self-diffusion coefficients follows Arrhenius behavior, with an activation energy of 35.9 kJ/mol for hydrogen and 46.3 kJ/mol for oxygen. The lifetime of OH − is about three to four times longer than that of H 2 O at temperatures from 800 to 2100 K, suggesting the greater stability of OH − over H 2 O, a unique characteristic of the grain boundaries. The estimated lifetime of the hydrous species decreases as the temperature increases, with an activation energy of 9.9 kJ/mol for OH − and 13.4 kJ/mol for H 2 O. While Gd 3 + is more mobile than Ti 4+ , both the Gd 3 + and Ti 4+ cations are orders of magnitude less mobile than the water species. The results suggest that water species are much more mobile within grain boundaries than in the bulk crystal and have the potential to penetrate deep into polycrystalline materials through grain boundaries, leading to grain boundary degradation and dissolution. The different mobilities of cations in complex oxides can lead to leaching of certain cations and incongruent dissolution during the chemical weathering of Earth and industrial materials.
Refractory multi-principal element alloys (RMPEAs) combine a chemically disordered lattice with a structurally ordered, single‐phase body‐centered cubic (bcc) crystal structure. Chemical fluctuations in these alloys give rise to significant energy barriers that impede dislocation motion. In this study, we use phase field dislocation dynamics to examine how spatial temperature fluctuations compete with randomness in energy barriers to affect the motion of long screw dislocations in three equi-atomic MoNbTa‐based RMPEAs: MoNbTa, MoNbTaW, and MoNbTaVW. Over a wide range of homologous temperatures (T h ≈ 0–0.6), we determined a screw dislocation ‘flow stress’ as the minimum applied stress to sustain continuous motion over a long excursion distance within a fixed timeframe. All three RMPEAs exhibit temperature dependent flow stress with three characteristic glide regimes: at low homologous temperatures, flow stress drops sharply, and screw glide remains planar and rectilinear; at intermediate homologous temperatures, the flow stress levels off as glide becomes planar but wavy; and at high homologous temperatures, flow stress plateaus as screws exhibit nonplanar, three‐dimensional motion. Glide kinetics and transition temperatures are controlled by the chemically induced fluctuations in the energy landscape. The rectilinear to wavy transition temperature is controlled by statistically weakest local barriers in the glide plane, whereas the wavy to 3D transition temperature is governed by statistically strongest local barriers. At low and intermediate homologous temperatures, the relative spread in energy barriers governs glide behavior by controlling the local kinetics of kink pair formation and kink pinning. At high homologous temperatures, the average barrier height governs the glide behavior by controlling the number of out-of-plane excursions during 3D glide. These findings reveal how random chemical fluctuations determine screw‐driven plasticity in RMPEAs, providing critical insight for the design of high temperature structural alloys.
Highlights: • Provides an overview of the current understanding of the many-body localization phase transition. • Discusses the subtleties of finite-size scaling near this transition • Assesses the implications of these subtleties for numerical studies of spin chains. • Explores scaling of diagnostics in models with known localization transitions. • Presents suggestions for future numerical work. We re-examine attempts to study the many-body localization transition using measures that are physically natural on the ergodic/quantum chaotic regime of the phase diagram. Using simple scaling arguments and an analysis of various models for which rigorous results are available, we find that these measures can be particularly adversely affected by the strong finite-size effects observed in nearly all numerical studies of many-body localization. This severely impacts their utility in probing the transition and the localized phase. In light of this analysis, we discuss a recent study (Šuntajs et al., 2020) of the behaviour of the Thouless energy and level repulsion in disordered spin chains, and its implications for the question of whether MBL is a true phase of matter.
Recent years have witnessed the broad use of carbon electrodes for electric double layer capacitors (EDLCs) because of large surface area, high porosity and low cost. Whereas experimental investigations are mostly focused on the device performance, computational studies have been rarely concerned with electrochemical properties at conditions remote from equilibrium, limiting their direct applications to materials design. Through a comprehensive analysis of extensive experimental data with various machine-learning methods, we report herein quantitative correlations between the structural features of carbon electrodes and the in-operando behavior of EDLCs including energy and power density. Machine learning allows us to identify important characteristics of activated carbons useful to optimize their efficiency in energy storage.
Here, the impact of hydrogen charge states on the stability and transport characteristics of hydrogen interstitials in alumina polymorphs is evaluated by multiscale computational methods including density functional theory (DFT), ab initio molecular dynamics (AIMD) and machine learned force fields. Thermodynamic calculations show that the protonic H i +1 interstitial is the most stable defect species for most values of the electronic bandgap in both and amorphous alumina (Al 2 O 3 ). Further, active learned Gaussian approximation potentials (GAP) were developed using AIMD data to study temperature dependent long time proton diffusion in alumina. Diffusivity calculations from GAP-MD simulations are found to be comparable with of the AIMD data, while being ~340 times faster and scalable to larger systems. Comparisons with diffusivity values for other interstitial charge states (H i 0 and H i -1 ) and published experimental literature indicate that H i +1 diffusion is the likely mechanism of hydrogen transport. A good agreement is obtained between H i +1 diffusivity calculated in α-Al 2 O 3 from DFT: 5.05 10 -3 exp(-0.81 eV/k B /T) cm 2 /s and reported experiment: 9.7X10 -4 exp(-0.83 eV/k B /T) cm 2 /s. Computationally and experimentally calculated energy barriers (0.81 and 0.83 eV respectively) only differ by 2.5%. Similarly, the pre-exponential diffusion coefficients only differ by 0.5 orders of magnitude. Moreover, the diffusivity of H i +1 in amorphous Al 2 O 3 in the 1000–2000 K range is calculated to be 2.53X10 -2 exp(-0.89 eV/k B /T), just one order of magnitude higher than the corresponding value in α-Al 2 O 3 . This suggests that local structural disorder does not significantly affect the energy landscape and diffusion behavior of H i +1 in Al 2 O 3 . Overall, these results show promise for the application of alumina polymorphs as hydrogen permeation barriers.