Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Averages”

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.

At least 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Bioenergy pathways within United States net-zero CO 2 emissions scenarios in the Energy Modeling Forum 37 study

The Energy Modeling Forum 37 study is organized around carbon dioxide (CO 2 ) mitigation scenarios reaching net-zero CO 2 emissions by 2050 in the United States. Here, this paper summarizes the potential contribution of bioenergy use in the electric power, transportation, industrial, and buildings sectors toward meeting that target based on model results. Thirteen modeling teams reported bioenergy consumption in the Reference and Net Zero scenarios. Consumption of bioenergy increased over time in the Reference scenario, from an average across models of 3.2 exajoules (EJ) in 2020 to 3.8 EJ in 2050. Average bioenergy consumption in 2050 increased further to 7.3 EJ in the Net Zero scenario. All scenarios that reach net-zero emissions required some form of carbon dioxide removal to offset emissions that are difficult to reduce. Carbon dioxide removal using bioenergy with CO 2 capture and storage (BECCS) varies widely across models, up to 1000 Mt CO 2 in 2050. Some models rely instead on direct air carbon capture and storage (DACCS), up to 2200 Mt CO 2 , and others use a combination of BECCS and DACCS. Model results show a strong inverse relationship between the amounts of BECCS and DACCS deployed. All modeling teams assumed a carbon sink from land use, land use change, and forestry, further offsetting a portion of emissions from fossil fuels and industry that are expensive to eliminate. Bioenergy consumption in 2050 decreased by an average of 1.5 EJ across eight models in a Net Zero+ scenario relative to the Net Zero scenario, due in part to a lower equilibrium carbon price resulting from optimistic cost assumptions for all energy technologies.

09 BIOMASS FUELS↗
DOE OSTI2026

Analyzing historical snow trends in interior Alaska

Study region The Chena River watershed in Interior Alaska, USA Study focus This study examines 40 years (water years 1982–2021) of snowpack characteristics to consider its hydrological implications in the 5350 km² Chena River basin. Using observations and a fine-scale physics model, we analyzed trends of snow water equivalent (SWE), snow onset and disappearance, and snow cover duration (SCD). New hydrological insights for the region Results indicate a decline in SWE across the modeled domain, averaging a decrease of 3 mm per decade, with larger decreases (up to 10 mm per decade) at lower elevations. While domain-averaged SWE trends were not statistically significant, observed SCD showed statistically significant decreases: −5.2, −5.0, and −4.4 days per decade at Teuchet Creek, Fairbanks F.O., and Little Chena Ridge, respectively. Notably, observations at SNOTEL stations and modeling revealed no statistically significant change in domain-averaged Rain-on-Snow (ROS) events over the 40-year period, contrasting some regional future estimates of increased ROS frequency. Peak streamflow did not consistently correlate with peak SWE levels, suggesting that other environmental factors such as ROS events and rapid temperature increases (e.g., a 10°C spike observed in 1992) are key drivers of hydrological outcomes. These findings improve understanding of complex subarctic hydrological processes impacting permafrost and highlight the need for adaptive water resource management to mitigate multi-factor risks like flooding and wildfire, requiring proactive planning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Reduced-Order Modeling Method for Phase-Change Thermal Energy Storage Heat Exchangers

Thermal energy storage can facilitate the effective utilization of renewable energy. To speed up the design process of thermal energy storage devices, it is critical to develop fast and accurate modeling methods for phase change material embedded heat exchangers (PCM HXs). This study developed and compared two approximation-assisted reduced-order PCM HX models for the simulation of thermal storage components and systems, which were verified against a validated finite-volume model. They are a pure black-box model, and a grey-box model based on the Number of Transfer Units (or effectiveness-NTU) approach. We used the reduced-order models to predict the performance of standalone PCM HX for 1000 cases with different fluid inlet conditions and HX designs. We also integrated the models with a vapor compression system to predict the compressor energy consumption and total charging time under various conditions. The results show that overall, the black-box model gave more accurate results than the grey-box model. On average, the mean absolute deviation in the PCM HX fluid outlet temperature was 0.05 K and 0.1 K for the black-box model and grey-box model, respectively. The grey-box model gave larger temperature deviation toward the end of the phase change process, due to the simplified two-node PCM representation. On the system level, the mean absolute deviation in compressor energy consumption was 0.2% and 0.3% for the black-box model and grey-box model, respectively. Their mean absolute deviation in total charging time was 1.1% and 2.6%, respectively. In terms of computation efficiency, the system simulation speedup ratio gained by using the reduced-order models was 11 to 57. System simulation time was decreased from an average of 1465 s to an average of 59 s. This shows that the proposed reduced-order modeling methods can be used to predict system performance of PCM HXs with less than 3% accuracy penalty, and 25 times less computational time than finite-volume models, enabling faster design and evaluation of PCM thermal storage devices.

25 ENERGY STORAGE↗
DOE OSTI2025

An interregional optimization approach for time series aggregation in continent-scale electricity system models

Modeling electric power systems with high shares of weather-dependent resources requires tradeoffs between temporal, spatial, and operational resolution. Many studies perform time series aggregation using clustering algorithms to reduce the temporal dimension, but when modeling continent-scale electricity systems that are large enough to contain multiple independent weather systems, this approach requires large numbers of representative periods to minimize errors in regional wind and solar capacity factors. Here, a new optimization-based approach for representative period selection and weighting is introduced that minimizes regional errors in average renewable capacity factors and electricity demand. The method delivers higher regional fidelity with fewer representative periods than alternative clustering methods when applied to wind, solar, and demand profiles for the contiguous United States. When representative periods are selected from multiple weather years, the optimized method reproduces regional averages with lower error than a complete 365-day time series from any single weather year. The method identifies only representative (as opposed to outlying) periods but can be combined with an iterative "stress period" identification approach to guide efficient decision-making considering both average and high-risk weather conditions.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Light-absorbing impurities accelerating glacial melting in southeastern Tibetan Plateau

Deposition of light-absorbing particles on glacier surfaces poses a series of adverse impacts on the cryospheric environment, climate and human health. Broad attention of the scientific community has been paid on insoluble light-absorbing impurities (ILAIs) in snow and ice on glaciers over the Tibetan Plateau (TP). However, systematic investigation of ILAIs in snowpack of glaciers on the TP is scarce. In this study, the properties and darkening effect of ILAIs in snowpack on glaciers are extensively investigated in the southeast of TP. Results show that ILAIs concentrations in multiple types of snow and ice samples were significantly different. Snowpit depths varied substantially from one profile to another during May and June 2016. The average concentrations of ILAIs in snowpits increase as snow melting progresses. Black carbon (BC) and dust cause snow albedo reduction more in snow with larger grain size Re. Based on a radiative transfer model calculation, the average albedo reduction induced by BC in the snowpack was 0.141 ± 0.02, and associated daily maximum radiative forcing (RF) was 72.97 ± 12.7 W m-2. BC is a controlling light-absorbing factor in snowpack and causes substantial albedo reduction and thus the associated daily maximum RF. The maximum reduction of snow cover duration was 4.56 ± 0.71 days caused by BC and dust in snowpack in southeastern TP. The average mass absorption cross-section (MAC) of BC from multiple snowpits was 3.26 ± 0.46 m2 g-1, which represents a typical value of MAC in snow on glaciers, but it is type-dependent of snow/ice samples. Tropospheric aerosols vertically extended up to 8 km over the TP and its surrounding areas, which indicates the transport of aerosols from remote sources through elevated pathways. A large amount of carbon stored in the brittle glaciers can be potentially released with meltwater runoff under a warming climate. This study provides a new insight for investigating carbonaceous and light-absorbing particles in glacierization areas.

Niu, Hewen↗
DOE OSTI2021

Initial plant community responses to hardwood control treatments in restoration of remnant longleaf pine (Pinus palustris) woodlands

Changes in land use over the past century have contributed to substantial losses of longleaf pine (Pinus palustris) woodlands in the southeastern USA and replacement with higher density, mixed pine and hardwood stands that suppress understory development and limit application of prescribed fire. To increase understanding of limiting factors and identify potential approaches for restoring longleaf pine woodlands, we studied in this work initial availability of light and soil water and 2-year (2018–2019) plant community responses after controlling overstory hardwoods in five remnant longleaf pine woodlands having no evidence of previous agriculture at the Savannah River Site near Aiken, SC, USA. Seven hardwood control treatments and a non-treated check were compared in a randomized compete block experiment: cutting, cutting + shredding of logging residues, stem injection with imazapyr herbicide, cutting + basal spray with imazapyr herbicide, cutting + basal spray with triclopyr herbicide, cutting + directed foliar spray with a mixture of glyphosate and imazapyr herbicides, and cutting + broadcast foliar spray with the same herbicide mixture. In the year prior to hardwood cutting (2016), understory light availability averaged 23% of full sun. Throughout the year of treatment installation (2017), average soil water content (SWC) was below permanent wilting point (5.5% SWC) where overstory hardwoods were retained; whereas, it was above wilting point (7.5%) where they had been cut. Thus, combined effects of shade and root competition from overstory hardwoods probably limited cover of herbaceous species in the non-treated check. In the 2 years following treatment, hardwood survival averaged > 99%, 62%, 42% and < 1% for the non-treated check, cutting, shredding, and herbicide treatments, respectively. Relative to the non-treated check, herbaceous species richness was increased by shredding, stem injection, basal spray, or directed foliar treatments; herbaceous cover was increased by basal spray or directed foliar treatments; and woody cover was decreased by directed or broadcast foliar treatments. The directed foliar spray was the most effective treatment for achieving many of the desired understory characteristics of a longleaf pine woodland, including a diverse understory dominated by herbaceous vegetation capable of supporting periodic prescribed fires. Although the imazapyr basal spray treatment increased cover of remnant woodland indicator species, none of the treatments increased richness of this species group, suggesting that restoration of native species composition will be expedited by enrichment seedings or plantings.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Recovery of DNA from fired and unfired cartridge casings: comparison of two DNA collection methods

For over 10 years, various studies have attempted to increase the recovery of DNA from ammunition by modifying the DNA collection, extraction, purification, and amplification procedures, with varying levels of success. This study focused on the “soaking” method of Montpetit & O’Donnell [1] and the “rinse-and-swab" method of Bille et al. [2]. First, testing for the presence of exogenous DNA, 210 boxed cartridges (brass, steel, and nickel-plated) from nine manufacturers were swabbed and DNA was extracted, concentrated, and quantified. Extracts that quantified > 0 ng/µL (44 of 210) were amplified and genotyped with GlobalFiler™. Of those, only one extract yielded two alleles indicating that the manufacturing and packaging of ammunition was virtually DNA free. Next, to obtain a baseline comparison of two DNA collection methods on a non-metallic substrate and identify a suitable number of cells to spot on cartridges, different DNA input amounts of primary human adult epidermal keratinocytes (HEKa) were tested. Thereafter, 300 brass and 300 nickel-plated, cartridges were spotted with HEKa cells containing ~5 ng of DNA, fired or unfired, and processed with either method. Finally, five methods representing hybrids of the soaking and rinse-and-swab methods were tested to determine if variations of those methods could be used to increase DNA yield and recovery. The results show that the soaking method consistently yielded more DNA than the rinse-and-swab method from a non-metallic substrate. However, the comparison study demonstrated that both methods performed comparably for cartridges. On average, the soaking method recovered 0.25 ng of DNA (5.1% recovery) and the rinse-and-swab method recovered 0.28 ng (5.8% recovery). However, average recoveries were significantly different among three analysts and considerable variation in yields were observed, possibly due to storage time. Furthermore, consistent with prior reports, the DNA recovered from brass casings was only 16% of that recovered from nickel-plated casings and the average yield of DNA from fired casings was reduced to 67% of unfired casings. Moreover, DNA extracts from brass or nickel-plated casings did not appear to contain amplification inhibitors and only 30/596 appeared severely degraded. Finally, both the published rinse-and-swab and soaking methods yielded more DNA than all modifications of the two methods. Overall, both methods yielded equivalent DNA quantities. Additionally, recovery of DNA from any given cartridge casing may be dependent on storage time as well as the skill, proficiency, and experience of the analyst and may reflect stochastic effects, particularly for casings containing low copy and/or degraded DNA.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Titanium isotope systematics of refractory inclusions: Echoes of molecular cloud heterogeneity

Calcium-aluminum-rich inclusions (CAIs) are highly refractory objects found in different chondrite groups and represent some of the oldest known solids of the Solar System. As such, CAIs provide key information regarding the conditions prevailing in the solar protoplanetary disk as well as subsequent mixing and transport processes. Many studies have investigated CAIs for their isotopic compositions and reported nucleosynthetic isotope anomalies in numerous elements, which are typically explained by the variable incorporation of isotopically highly anomalous presolar phases. However, with the exception of 54 Cr-enriched nanospinels, the exact presolar phases responsible for the isotopic heterogeneities are yet to be identified. To address this issue, we here present in-situ Ti isotopic analyses obtained on a diverse set of CAIs from various CV3 chondrites. In this work, the in-situ measurements were performed by targeting individual mineral phases of 15 CAIs with laser-ablation mass spectrometry and indicate significant inter- and intra-CAI isotopic heterogeneity in the neutron-rich isotope 50 Ti. This is particularly pronounced for primitive fine-grained CAIs, whereas coarse-grained CAIs, which have been subject to melting, exhibit smaller degrees of Ti isotopic heterogeneity. To further investigate this Ti isotopic heterogeneity, we additionally obtained Ti isotopic compositions of sequential acid leachates from two fine-grained and two coarse-grained CAIs derived from CV3 chondrites. In contrast to potential expectations from the first part of the study, we do not observe any significant intra-CAI Ti isotopic heterogeneity between the different leaching steps. The lack of intra-CAI Ti isotopic heterogeneity in the acid leachate samples of this study likely reflects that the leaching procedure is unable to efficiently separate the carriers of isotopically anomalous Ti in CAIs. By comparing the bulk CAI Ti isotope compositions with Ti isotope data for hibonite-rich objects from the literature, we find that the range of Ti isotope compositions recorded by CAIs from various chondrite groups can be accounted for by the averaging of hibonite grains. In turn, the variable Ti isotope compositions of hibonite grains can be explained by the averaging of isotopically diverse presolar grains present in the Sun’s parental molecular cloud. This effect of averaging is statistically supported by the central limit theorem, and the concept has the potential to be useful for other isotopic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Functional variability in specific root respiration translates to autotrophic differences in soil respiration in a temperate deciduous forest

CO 2 release from forest soils (R s ) is a prominent flux in the global carbon cycle. Rs is derived from roots (autotrophic respiration, R a ) and microbial (heterotrophic) respiration and is highly dynamic, as it depends on edaphic and environmental conditions as well as root functional traits and microbial community composition. It is unclear how root functional traits affect root and microbial respiration rates; however, their consideration may help parse out the relative contributions of root and microbial respiration to R s . At a temperate forest site, root systems of 3–4 functional root orders and their surrounding surface soil were carefully excavated and placed into custom trays designed to repeatedly measure R s in situ on eight temperate tree species that varied in their root functional strategies and mycorrhizal affinity. R s was measured bi-weekly to monthly for nearly one year using a custom chamber attached to a gas exchange system. R s varied over time, ranging from 0.3 to 12 µmol m -2 s -1 . Comparable root systems of the same species were excised from the soil and specific root respiration rates (R r ) were measured. Rr ranged from 2.5 to 9.0 nmol g -1 s -1 and was negatively correlated with root tissue density and positively related to root tissue nitrogen concentration. Using R r to estimate R a , we estimate that R a accounts for <10%, on average 2–3%, of R s for individual root systems (averaging 1.2 g dry biomass) housed in surrounding soil (average 1.3 kg dry mass) in situ; thus, Ra was roughly 20 times greater than Rh per unit mass. The contribution of R a peaked in the fall and coincided with leaf senescence of the forest canopy. A soil-sterilizing experimental treatment designed to help isolate R a in situ reduced bacterial biomass and shifted fungal community composition, but there was no reduction in Rs of the in-situ root-soil tray systems. The relative R a to R s ratio increased with root functional strategies characterized by greater specific root length and tip abundance, but also to greater root tissue density. The ratio of R a to R s also increased with warmer soil temperatures and decreased slightly with increasing soil moisture. We discuss how incorporating root functional traits as modulators of the autotrophic contribution to R s could be considered when modeling total soil CO 2 efflux from forests.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Multi-site evaluation of stratified and balanced sampling of soil organic carbon stocks in agricultural fields

Estimating soil organic carbon (SOC) stocks in agricultural fields is essential for environmental and agronomic research, management, and policy. Stratified sampling is a classic strategy for estimating mean soil properties, and has recently been codified in SOC monitoring protocols. However, for the specific task of estimating the SOC stock of an agricultural field, concrete guidance is needed for which covariates to stratify on and how much stratification can improve estimation efficiency. It is also unknown how stratified sampling of SOC stocks compares to modern alternatives, notably doubly balanced sampling. To address these gaps, we collected high-density (average of 7 samples ha -1 ) and deep (average of 75 cm) measurements of SOC stocks at eight commercial fields under maize-soybean production in two US Midwestern states. We combined these measurements with a Bayesian geostatistical model to evaluate stratified and balanced sampling strategies that use a set of readily-available geographic, topographic, spectroscopic, and soil survey data. We examined the number of samples needed to achieve a given level of SOC stock estimation accuracy. While stratified sampling using these variables enables an average sample size reduction of 17% (95% CI, 11% to 23%) compared to simple random sampling, doubly balanced sampling is consistently more efficient, reducing sample sizes by 32% (95% CI, 25% to 37%). The data most important to these efficiency gains are a remotely-sensed SOC index, SSURGO estimates of SOC stocks, and the topographic wetness index. We conclude that in order to meet the urgent challenge of climate change, SOC stocks in agricultural fields could be more efficiently estimated by taking advantage of this readily-available data, especially with doubly balanced sampling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Regional-scale soil carbon predictions can be enhanced by transferring global-scale soil–environment relationships

Accurate modelling and mapping soil organic carbon are crucial for supporting soil health restoration and climate change mitigation at both regional and global scales. However, regional soil predictions often suffer from data scarcity and high prediction uncertainty. Utilizing a pre-trained global-to-regional soil carbon predictive model can be a potential solution to address this challenge. Despite its promise, how to construct and apply the global-scale model to enhance regional-scale soil carbon mapping remains largely unexplored. Here, we propose the Global Soil Carbon Pre-trained Model (GSoilCPM), a deep-learning-based domain adaptative model, to enhance regional-scale soil carbon predictions. Based on large amount of environmental covariate data and 106,167 soil samples across the globe, we verify our hypothesis of the effectiveness of this 'global-to-regional' modelling strategy. The pre-trained model can be then transferred and fine-tuned to bridge the regional- and global-scale soil–environment relationships. We applied and validated this modelling strategy in four regional-scale study areas, three in the Northern Hemisphere and one in the Southern Hemisphere, each with distinct environmental background. Compared to traditional modelling approaches as a baseline, four case studies all demonstrated significant improvement in prediction accuracy across diverse environments and varying data availabilities. The average percentage improvement across all regions is 10.93% (absolute values decreased by 1.20 g kg−1 averagely) in MAE and 29.04% (absolute values increased by 0.10 averagely) in CCC. The applicability and future horizons of using GSoilCPM were further discussed. We further reveal that regions with fewer soil samples or lower baseline accuracy benefit more from the pre-trained global model. Our findings highlight the advantages of leveraging the generalized knowledge from global models to enhance specifically localized soil modelling, positioning a potential paradigm shift in digital soil mapping, and far-reaching implications for soil monitoring and land management.

Deep learning↗
DOE OSTI2024

Investigation of precooling unit design options in hydrogen refueling station for heavy-duty fuel-cell electric vehicles

Precooling gaseous hydrogen fuel to a cold temperature before refueling a heavy-duty (HD) hydrogen fuel cell electric vehicle (HFCEV) is essential to avoid overheating the vehicle tank, as well as achieving a high state of charge (SOC). Because a large volume of hydrogen is dispensed during each fill, the need for a shorter fill time amplifies the need to precool each load for refueling a HFCEV. Thus, the design and operation of a precooling unit (PCU), as well as the associated capital and operating costs, plays a pivotal role in any plans for heavy-duty hydrogen refueling stations. Here, in this paper, we present a thermodynamic and technoeconomic analysis of a PCU in a gaseous hydrogen refueling station (HRS) for HD HFCEVs. By employing Argonne National Laboratory's hydrogen station cost optimization and performance evaluation (H2SCOPE) model, the refueling of 50 kg of hydrogen on-board type IV tank at ambient temperatures of 15–45 °C and varying fill rates is simulated. The required degree of precooling temperature to obtain either 100% or the maximum possible SOC is obtained from the simulation. Additionally, the simulation results demonstrate that the average flow rate of hydrogen is approximately 40% lower than the maximum flow rate during a typical fill; which motivates further evaluation of the instantaneous hydrogen mass flow rate profile and suggests the scope of improving precooling unit design. Accordingly, a hybrid strategy of precooling hydrogen has been proposed to address the cooling load by sizing the refrigeration unit for the average flow rate of hydrogen, while supplementing the above average peak hydrogen flow cooling load through thermal buffering. The combined technique enables the downsizing of the original PCU capacity by 25–40% and demonstrates a potential cost reduction of the PCU by approximately 30%, which translates to an installed cost reduction of ∼$125,000 per dispenser.

25 ENERGY STORAGE↗
DOE OSTI2023

Energy-efficient cooperative resource allocation and task scheduling for Internet of Things environments

Offloading Internet of Things (IoT) tasks to the cloud for further processing might not always lead to an optimal execution time, particularly in situations such as resource contention, under-provisioning, over-provisioning, and fragmentation. In addition, dynamically optimizing the number of Virtual Machines (VMs) for resource scheduling in order to meet application requirements remains a major research challenge. Further, existing resource scheduling algorithms focus primarily on minimizing operational costs while maximizing resource sharing and utilization. Considering energy utilization as part of the resource allocation and scheduling process as an optimization objective for maintaining load balancing has often been neglected. To address these challenges and more, we propose a cooperative energy-aware resource allocation and scheduling strategy based on a Technique for Order of Preference by Similarity to Ideal Solution (TOPSIS) multi-criteria decision-making method. Here we used the Grid Workloads Archive dataset to evaluate our proposed approach named TOPREAL. Experimental results with respect to the allocation of VM resources when considering processing a large segment of tasks indicate that TOPREAL outperforms existing algorithms in terms of energy savings, with an average improvement of 40.25%, while maintaining an average improvement of 16.21% when it comes to execution time. Results also demonstrate that our method can save an average of 78.06 processing hours and 63,215kJ of energy when compared to existing scheduling algorithms. These results demonstrate the effectiveness of our proposed model and the viability of using multi-criteria decision-making techniques such as TOPSIS to solve the resource allocation and scheduling problem in edge environments.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

The impact of Uber and Lyft on vehicle ownership, fuel economy, and transit across U.S. cities

We estimate the effects of transportation network companies (TNCs) Uber and Lyft on vehicle ownership, fleet average fuel economy, and transit use in U.S. urban areas using a set of difference-in-difference propensity score-weighted regression models that exploit staggered market entry across the U.S. from 2011 to 2017. We find evidence that TNC entry into urban areas causes an average 0.7% increase in vehicle registrations with significant heterogeneity in these effects across urban areas: TNC entry produces larger vehicle ownership increases in urban areas with higher initial ownership (car-dependent cities) and in urban areas with lower population growth (where TNC-induced vehicle adoption outpaces population growth). We also find no statistically significant average effect of TNC entry on fuel economy or transit use but find evidence of heterogeneity in these effects across urban areas, including larger transit ridership reductions after TNC entry in areas with higher income and more childless households.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Analysis and prediction of reaction kinetics using the degree of rate control

“Degree of rate control” (DRC) analysis provides a quantitative approach for analysing the kinetics of multi-step reaction mechanisms that has been widely applied to both heterogeneous and homogeneous catalysis research, as well as electrocatalysis. The DRC of any given transition state or intermediate is defined as a partial derivative such that it approximately equals the fractional increase in net rate to the product of interest per differential decrease in its standard-state free energy for that species (÷RT), holding constant the standard-state free energies of all other transition states and intermediates. Even very complex mechanisms usually have only a few species with non-zero DRCs and are thus the species whose interactions with the catalyst most strongly affect the net rate. These key DRC values thus offer a simple and intuitive route to optimize catalyst materials, especially with the assistance of computational methods like density functional theory (DFT). These high-DRC species are also the species whose energetics must be most accurately measured or calculated to achieve an accurate kinetic model for any reaction mechanism. In simple cases with a single “rate-determining step”, the DRC for its transition state (TS) is + 1. Catalyst-bound intermediates, on the other hand, often have negative DRCs equal to a small integer times their fractional occupancy of catalyst sites. The apparent activation energy equals a weighted average of the standard-state enthalpies (relative to reactants) of all the species (intermediates, transition states and products) in the reaction mechanism, each weighted by its DRC (+RT). It has been shown that the apparent transfer coefficient in electrocatalysis, an inverted form of the Tafel slope, is a weighted average of the number of electrons transferred to generate each intermediate or product species in the mechanism, weighted again by the DRC. Quantitative analysis of kinetic isotope effects (KIEs, or the ratio of net rates for different reactant isotopes) in complex mechanisms has shown that the logarithm of the KIE equals the weighted average over all species in the mechanism of the difference between the two isotopes in their standard-state free energies (÷RT), again weighted by the DRC. The reaction orders with respect to fluid-phase concentrations of reactants, products and intermediates have also been proven to be directly related to DRCs. Thus, there are numerous experimental observables which equate to short linear combinations of DRCs, so that combinations of experimental measurements might provide access to DRC values. Since its invention for transition states in 1994 and its generalization to include intermediates in 2009, the DRC has thus far mainly been calculated for microkinetic models based either entirely on DFT or on DFT with the key energies (i.e., those for high-DRC species) being fine-tuned to match experiments. The relationships summarized in this work provide new opportunities for using experiments earlier in the development and optimization of microkinetic models that require input from computational catalysis.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Microdynamics of active particles in defect-rich colloidal crystals

Hypothesis. Because they are self-propulsive, active colloidal particles can interact with their environment in ways that differ from passive, Brownian particles. Here, we explore how interactions in different microstructural regions may contribute to colloidal crystal annealing. Experiments. We investigate active particles propagating in a quasi-2D colloidal crystal monolayer produced by alternating current electric fields (active-to-passive particle ratio ~ 1:720). The active particle is a platinum Janus sphere propelled by asymmetric decomposition of hydrogen peroxide. Crystals are characterized for changes in void properties. The mean-squared-displacement of Janus particles are measured to determine how active microdynamics depend on the local microstructure, which is comprised of void regions, void-adjacent regions (defined as within three particle diameters of a void), and interstitial regions. Findings. At active particle energy E A = 2.55 k B T, the average void size increases as much as three times and the average void anisotropy increases about 40% relative to the passive case. Further, the average microdynamical enhancement, <δ(t)>, of Janus particles in the crystal relative to an equivalent passive Janus particle is reduced compared to that of a free, active particle (<δ(t) > is 1.88 ± 0.04 and 2.66 ± 0.08, respectively). The concentration of active particles is enriched in void and void-adjacent regions. Active particles exhibit the greatest change in dynamics relative to the passive control in void-adjacent regions (<δ(t)> = 2.58 ± 0.06). The results support the conjecture that active particle microdynamical enhancement in crystal lattices is affected by local defect structure.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Explicit physics-informed neural networks for nonlinear closure: The case of transport in tissues

In upscaling methods, closures for nonlinear problems present a well-known challenge. While a number of theoretical methods have been proposed for handling such closures, nonlinearities still remain a significant obstacle for many problems. In this work, we use a combination of formal upscaling and data-driven machine learning for explicitly closing a nonlinear transport and reaction process in multiscale tissues. The classical effectiveness factor model is used to formulate the macroscale reaction kinetics. We train a multilayer perceptron network using training data generated by direct numerical simulations over microscale examples. Once trained, the network is used in an algorithm for numerically solving the upscaled (coarse-grained) differential equation describing mass transport and reaction in two example tissues. The network is described as being explicit in the sense that the network is trained using macroscale concentrations and gradients of concentration as components of the feature space rather than incorporating them as part of a constraint in the optimization process. Network training and solutions to the macroscale transport equations were computed for two different tissues. The two tissue types (brain and liver) exhibit markedly different geometrical complexity and spatial scale (cell size and sample size). The upscaled solutions for the average concentration are compared with numerical solutions derived from the microscale concentration fields by a posteriori averaging. There are three outcomes of this work of particular note. 1) Our overall approach results in an upscaled nonlinear PDE. The PDE is closed using a neural network, and our approach results in the definition of the classical effectiveness factor for effecting closure. 2) We identify particular source terms for the closure problem that are important for representing the structure of the closure. These source terms involve macroscale concentrations and their gradients. We adopt these source terms to use as explicit features in the learning algorithm. We find the trained networks that include the macroscale source terms generate models that are able to predict the correction factor with increased fidelity over those that do not. 3) We find that the trained network exhibits good generalizability, and it is able to predict the effectiveness factor with high fidelity for realistically-structured tissues despite the significantly different scale and geometrical complexity of the two example tissue types. This latter result emphasizes our purposeful connection between conventional averaging methods with the use of machine learning for closure; this contrasts with some machine learning methods for upscaling where the exact form of the macroscale equation remains unknown.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Local conservation of energy in fully implicit PIC algorithms

We consider the issue of strict, fully discrete local energy conservation for a whole class of fully implicit local-charge- and global-energy-conserving particle-in-cell (PIC) algorithms. Earlier studies demonstrated these algorithms feature strict global energy conservation. However, whether a local energy conservation theorem exists (in which the local energy update is governed by a flux balance equation at every mesh cell) for these schemes is unclear. In this study, we show that a local energy conservation theorem indeed exists. We begin our analysis with the 1D electrostatic PIC model without orbit-averaging, and then generalize our conclusions to account for orbit averaging, multiple dimensions, and electromagnetic models (Darwin). In all cases, a temporally, spatially, and particle-discrete local energy conservation theorem is shown to exist, proving that these formulations (as originally proposed in the literature), in addition to being locally charge conserving and globally energy conserving, are strictly locally energy conserving as well. In contrast to earlier proofs of local conservation in the literature, which only considered continuum time, our result is valid for the fully implicit time-discrete version of all models considered, including important features such as orbit averaging. We demonstrate the local-energy-conservation property numerically with a paradigmatic numerical example.

97 MATHEMATICS AND COMPUTING↗