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.
Being the “mother distributions” of all types of two-parton correlation functions, generalized transverse momentum dependent parton distributions (GTMDs) have attracted a lot of attention over the last years. We argue that exclusive double production of pseudoscalar quarkonia (η c or η b ) in nucleon-nucleon collisions gives access to GTMDs of gluons.
Triple-conducting materials have been proved to improve the performance of popular protonic ceramic electrolysis cells. However, partially because of the complexity of the water splitting reaction involving three charge carriers, that is, oxygen (O 2– ), proton (H + ), and electron (e – ), the triple-conducting reaction mechanism was not clear, and the reaction conducting pathways have seldom been addressed. In this study, the triple conducting Ruddlesden–Popper phase Pr 1.75 Ba 0.25 NiO 4+δ as an anode on the BaCe 0.7 Zr 0.1 Y 0.1 Yb 0.1 O 3–δ electrolyte was fabricated and its electroresponses were characterized by electrochemical impedance spectroscopy with various atmospheres and temperatures. The impedance spectra are deconvoluted by means of the distribution of the relaxation time method. The surface exchange rate and chemical diffusivity of H + and O 2– are characterized by electrical conductivity relaxation. The physical locations of electrochemical processes are also identified by atomic layer deposition with a surface inhibitor. A microkinetics model is proposed toward conductivities, triple-conducting pathways, reactant dependency, surface exchange and bulk diffusion capabilities, and other relevant properties. Lastly, the rate-limiting steps and suggestions for further improvement of electrode performance are presented.
Solvolytic extraction of lignin from biomass is a critical step in lignin-first biorefining, including the reductive catalytic fractionation (RCF) process. Key to optimal RCF processing is the ability to rapidly extract lignin from biomass at high delignification extents and transfer the lignin molecules to a catalyst surface in a time frame that minimizes lignin condensation reactions. Here, we use a flow-through reactor to study the effects of temperature (175-250 °C), residence time (9 to 36 min), and solvent composition (methanol and methanol-water) on lignin extraction and condensation. We evaluated three metrics at each condition: total delignification, delignification rate, and extent of condensation, the latter measured by a decrease in monomer yield for batch hydrogenolysis reactions of solvolysis liquor compared to batch RCF reactions. We observe that delignification is predominantly determined by temperature, while residence time dictates the lignin condensation extent. Moreover, the extent of both extraction and condensation increased in the methanol-water solvent system compared to that in the methanol system. Lignin extracted in methanol is stable up to 18-min residence times at or below 225 °C, while a majority of the lignin extracted in methanol-water is condensed with a 9-min residence time at 200 °C. These results can inform reactor designs and solvent selection for lignin-first biorefining processes that aim to physically separate the biomass and catalyst.
Hyporheic exchange induced by riverbed topography and roughness provides important ecosystem functions. We investigated how cobble clusters embedded in a finer-grained sand bed affect the near-bed channel flow and the exchange of surface and subsurface water. We tested how the spacing and embeddedness of cobbles altered hyporheic exchange through a three-dimensional fully coupled surface-subsurface model. The 3D modeling framework captured the full physics of the exchange process locating the lateral position of upwelling zones on the side of the cobble. As the cobbles protrude more into the channel, eddies appear downstream, with smaller vortices where cobbles are closer together. In our simulations, hyporheic exchange increased with the spacing and protrusion ratio. The travel time in the hyporheic zone also increased with the protrusion ratio because the hyporheic flow paths are longer when the cobble is more protruded. With different spacings at the highest protrusion ratio, the average slope of early time breakthrough curve was steeper for large spacing settings: the travel time increases as the spacing decreases.
Heavy precipitation, often associated with weather phenomena such as tropical cyclones, extratropical cyclones (ETCs), atmospheric rivers (ARs), and mesoscale convective systems (MCSs), can cause significant socio-economic loss. Here, in this study, we apply atmospheric feature trackers to quantify the contributions of these storm types in observational data sets and climate model short-range hindcasts. We generate a global hourly storm data set at 0.25° spatial resolution covering 2006–2020, based on the tracking results from TempestExtremes and Python FLEXible object TRacKeR. Our analyses show that these four storm types account for 67% of global annual mean precipitation and 82% of top 1% precipitation extremes, with MCSs mainly over the tropics, and ARs and ETCs over the midlatitudes. The percentage of precipitation contributions from these storms also show strong seasonality over many geographical locations. We further apply the tracking results to the Energy Exascale Earth System Model (E3SM) short-range hindcasts and evaluate how well these storms are simulated. The evaluation show that E3SM, with ∼1° resolution, significantly underestimates storm-associated precipitation totals and extremes, especially for MCSs in the tropics. Our analysis also suggests that model fails to capture the correct mean diurnal phases and amplitude of MCS precipitation. This phenomenon-based approach provides a better understanding of precipitation characteristics and can lead to enhanced model evaluation by revealing underlying problems in model physics related to precipitation processes associated with the heavy-precipitating storms.
Methane (CH 4 ) is globally the second most critical greenhouse gas after carbon dioxide, contributing to 16%–25% of the observed atmospheric warming. Wetlands are the primary natural source of methane emissions globally. However, wetland methane emission estimates from biogeochemistry models contain considerable uncertainty. One of the main sources of this uncertainty arises from the numerous uncertain model parameters within various physical, biological, and chemical processes that influence methane production, oxidation, and transport. Sensitivity Analysis (SA) can help identify critical parameters for methane emission and achieve reduced biases and uncertainties in future projections. This study performs SA for 19 selected parameters responsible for critical biogeochemical processes in the methane module of the Energy Exascale Earth System Model (E3SM) land model (ELM). The impact of these parameters on various CH 4 fluxes is examined at 14 FLUXNET- CH 4 sites with diverse vegetation types. Given the extensive number of model simulations needed for global variance-based SA, we employ a machine learning (ML) algorithm to emulate the complex behavior of ELM methane biogeochemistry. We found that parameters linked to CH 4 production and diffusion generally present the highest sensitivities despite apparent seasonal variation. Comparing simulated emissions from perturbed parameter sets against FLUXNET-CH 4 observations revealed that better performances can be achieved at each site compared to the default parameter values. This presents a scope for further improving simulated emissions using parameter calibration with advanced optimization techniques.
As a confined thin sheet crumples, it spontaneously segments into flat facets delimited by a network of ridges. Despite the apparent disorder of this process, statistical properties of crumpled sheets exhibit striking reproducibility. Experiments have shown that the total crease length accrues logarithmically when repeatedly compacting and unfolding a sheet of paper. Here, we offer insight to this unexpected result by exploring the correspondence between crumpling and fragmentation processes. We identify a physical model for the evolution of facet area and ridge length distributions of crumpled sheets, and propose a mechanism for re-fragmentation driven by geometric frustration. This mechanism establishes a feedback loop in which the facet size distribution informs the subsequent rate of fragmentation under repeated confinement, thereby producing a new size distribution. We then demonstrate the capacity of this model to reproduce the characteristic logarithmic scaling of total crease length, thereby supplying a missing physical basis for the observed phenomenon.
Abstract The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) was a yearlong expedition supported by the icebreaker R/V Polarstern , following the Transpolar Drift from October 2019 to October 2020. The campaign documented an annual cycle of physical, biological, and chemical processes impacting the atmosphere-ice-ocean system. Of central importance were measurements of the thermodynamic and dynamic evolution of the sea ice. A multi-agency international team led by the University of Colorado/CIRES and NOAA-PSL observed meteorology and surface-atmosphere energy exchanges, including radiation; turbulent momentum flux; turbulent latent and sensible heat flux; and snow conductive flux. There were four stations on the ice, a 10 m micrometeorological tower paired with a 23/30 m mast and radiation station and three autonomous Atmospheric Surface Flux Stations. Collectively, the four stations acquired ~928 days of data. This manuscript documents the acquisition and post-processing of those measurements and provides a guide for researchers to access and use the data products.
Lake trophic state is a key ecosystem property that integrates a lake’s physical, chemical, and biological processes. Despite the importance of trophic state as a gauge of lake water quality, standardized and machine-readable observations are uncommon. Remote sensing presents an opportunity to detect and analyze lake trophic state with reproducible, robust methods across time and space. We used Landsat surface reflectance data to create the first compendium of annual lake trophic state for 55,662 lakes of at least 10 ha in area throughout the contiguous United States from 1984 through 2020. The dataset was constructed with FAIR data principles (Findable, Accessible, Interoperable, and Reproducible) in mind, where data are publicly available, relational keys from parent datasets are retained, and all data wrangling and modeling routines are scripted for future reuse. Together, this resource offers critical data to address basic and applied research questions about lake water quality at a suite of spatial and temporal scales.
Explore the source record for details and available documents.
In this paper, we provide a numerical tool to study a material’s coherence from a set of 2D Lagrangian trajectories sampling a dynamical system, i.e., from the motion of passive tracers. We show that eigenvectors of the Burau representation of a topological braid derived from the trajectories have levelsets corresponding to components of the Nielsen–Thurston decomposition of the dynamical system. One can thus detect and identify clusters of space–time trajectories corresponding to coherent regions of the dynamical system by solving an eigenvalue problem. Unlike previous methods, the scalable computational complexity of our braid-based approach allows the analysis of large amounts of trajectories.
Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.
The wind-energy (WE) industry relies on numerical weather prediction (NWP) forecast models as foundational or base models for many purposes, including wind-resource assessment and wind-power forecasting. During the Second Wind Forecast Improvement Project (WFIP2) in the Columbia River Basin of Oregon and Washington, a significant effort was made to improve NWP forecasts through focused model development, to include experimental refinements to the High Resolution Rapid Refresh (HRRR) model physics and horizontal grid spacing. In this study, the performance of an experimental version of HRRR that includes these refinements is tested against a control version, which corresponds to that of the operational HRRR run by National Oceanic and Atmospheric Administration/National Centers for Environmental Protection at the outset of WFIP2. Furthermore, the effects of horizontal grid resolution were also tested by comparing wind forecasts from the HRRR (with 3-km grid spacing) with those from a finer-resolution HRRR nest with 750-m grid spacing. Model forecasts are validated against accurate wind-profile measurements by three scanning, pulsed Doppler lidars at sites separated by a total distance of 71 km. Model skill and improvements in model skill, attributable to physics refinements and improved horizontal grid resolution, varied by season, by site, and during periods of atmospheric phenomena relevant to WE. In general, model errors were the largest below 150 m above ground level (AGL). Experimental HRRR refinements tended to reduce the mean absolute error (MAE) and other error metrics for many conditions, but degradation in skill (increased MAE) was noted below 150 m AGL at the two lowest-elevation sites at night. Finer resolution was found to produce the most significant reductions in the error metrics.
The spatial resolution and imaging characteristics of plastic scintillators are characterized using laser-driven proton beams. Laser-driven proton beams typically have broad energy spectra and are accompanied by relativistic electrons and high-energy photons, both potentially contributing to background noise. Different types and thicknesses of Eljen Technology scintillators are compared to determine their intrinsic point spread function. Point-projection imaging of a mesh is used to compare the imaging resolution of the scintillator to the usual imaging detector, radiochromic film, and is found to be reasonably comparable and sufficient for many experimental applications.
Bulk metallic glasses (BMGs) have successfully been used to replicate molds that are structured at the nano- and even atomic scale through thermoplastic forming (TPF), an ability that was speculated to be rooted in the glass’ featureless atomic structure. These previous demonstrations of atomically precise imprinting, however, were performed under conditions where mold atomic feature dimensions coincided with the unit cell size of constituents in the BMG. In order to evaluate if accurate atomic-scale replication is possible in general, i.e., independent of the accidental presence of favorable constituent size/feature size relationships, we have used Pt 57.5 Cu 14.7 Ni 5.3 P 22.5 to replicate three different crystalline facets of LaAlO 3 single crystals, each exposing distinct atomic step heights. We find that in all cases, the terraced surface termination can be copied with remarkable fidelity, corroborating that BMGs when thermoplastic formed are capable of adapting to any externally imposed confinement with sub-angstrom precision without being limited by factors related to the specifics of their internal structure. This unprecedented capability of quasi-limitless replication fidelity reveals that the deformation mechanism in the supercooled liquid state of BMGs is essentially homogeneous and suggests TPF of BMGs to be a versatile toolbox for atomic and precision nanoscale imprinting.
Here, recent work shows that strong stability and dimensionality freedom are essential for robust numerical integration of thermostatted ringpolymer molecular dynamics (T-RPMD) and path-integral molecular dynamics, without which standard integrators exhibit non-ergodicity and other pathologies [R. Korol et al., J. Chem. Phys. 151, 124103 (2019) and R. Korol et al., J. Chem. Phys. 152, 104102 (2020)]. In particular, the BCOCB scheme, obtained via Cayley modification of the standard BAOAB scheme, features a simple reparametrization of the free ring-polymer sub-step that confers strong stability and dimensionality freedom and has been shown to yield excellent numerical accuracy in condensed-phase systems with large time steps. Here, we introduce a broader class of T-RPMD numerical integrators that exhibit strong stability and dimensionality freedom, irrespective of the Ornstein–Uhlenbeck friction schedule. In addition to considering equilibrium accuracy and time step stability as in previous work, we evaluate the integrators on the basis of their rates of convergence to equilibrium and their efficiency at evaluating equilibrium expectation values. Within the generalized class, we find BCOCB to be superior with respect to accuracy and efficiency for various configuration-dependent observables, although other integrators within the generalized class perform better for velocity-dependent quantities. Extensive numerical evidence indicates that the stated performance guarantees hold for the strongly anharmonic case of liquid water. Both analytical and numerical results indicate that BCOCB excels over other known integrators in terms of accuracy, efficiency, and stability with respect to time step for practical applications.
We report on systematic changes to the adsorption geometry of the dye N3 {[cis-bis(isothiocyanato)bis(2,2' -bipyridyl-4,4' -dicarboxylato ruthenium(II)]} on a gold substrate as the pH of the deposition environment is altered. The protonation states of the four --COOH groups of the N3 dye change according to the modified pH conditions, thus affecting the number of --COOH and --NCS functional groups that participate in the adsorption to gold. Here, we use heterodyne detected vibrational sum frequency generation (HD-VSFG) spectroscopy to obtain surface specific vibrational information on both --COOH and --NCS groups as a function of pH of the deposition conditions. Polarization-dependent HD-VSFG yields sets of complex χ (2) spectra, enabling us to perform a simultaneous fitting procedure to the polarization-dependent real and imaginary components and thus extract detailed structural information of the N3/gold interface. Overall, our results show that N3 preferentially adsorbs to gold either with two --COOH groups and one --NCS group in more acidic conditions or with one --COOH group and two --NCS groups in more basic conditions.
We report the recent boom in computational chemistry has enabled several projects aimed at discovering useful materials or catalysts. We acknowledge and address two recurring issues in the field of computational catalyst discovery. First, calculating macro-scale catalyst properties is not straightforward when using ensembles of atomic-scale calculations [e.g., density functional theory (DFT)]. We attempt to address this issue by creating a multi-scale model that estimates bulk catalyst activity using adsorption energy predictions from both DFT and machine learning models. The second issue is that many catalyst discovery efforts seek to optimize catalyst properties, but optimization is an inherently exploitative objective that is in tension with the explorative nature of early-stage discovery projects. In other words, why invest so much time finding a “best” catalyst when it is likely to fail for some other, unforeseen problem? We address this issue by relaxing the catalyst discovery goal into a classification problem: “What is the set of catalysts that is worth testing experimentally?” Here, we present a catalyst discovery method called myopic multiscale sampling, which combines multiscale modeling with automated selection of DFT calculations. It is an active classification strategy that seeks to classify catalysts as “worth investigating” or “not worth investigating” experimentally. Our results show an ~7–16 times speedup in catalyst classification relative to random sampling. These results were based on offline simulations of our algorithm on two different datasets: a larger, synthesized dataset and a smaller, real dataset.