Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “siting and sizing”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Geothermal district energy systems coupled with seasonal underground thermal energy storage: a U.S. techno-economic screening by climate and geology

In the United States, cooling-dominated commercial building loads can cause geothermal heat pump-based district energy systems to accumulate a long-term subsurface thermal imbalance, motivating the incorporation of seasonal underground thermal energy storage. We developed a transferable workflow to evaluate geothermal district systems that pair ground heat exchangers with seasonal underground thermal energy storage. Using standardized hourly loads for seven commercial buildings and a uniform cost framework, we simulated ten U.S. cities with a physics-based ground heat exchanger model, subsurface storage simulations, and economic assessment to isolate the roles of climate and hydrogeology. In cooling-dominated cities, underground thermal energy storage supplied the majority of annual cooling, cutting electricity use and summer peaks substantially while achieving levelized costs comparable to or below conventional chiller-boiler plants. In cooler climates, the storage share shrunk, required borefield size and costs rose, and levelized cost of energy increased nearly linearly with declining underground thermal energy storage fraction, indicating storage fraction as the primary economic lever. Sensitivity analysis showed capital risk dominated by borefield drilling and surface heating, ventilation, and air-conditioning and piping, with underground thermal energy storage costs secondary. This workflow provides a transparent foundation for site-specific design and screening of next-generation geothermal district energy systems.

25 ENERGY STORAGE↗
DOE OSTI2026

Designing the Platinum Catalyst Layer for Improved Performance and Durability in Anion Exchange Membrane Water Electrolysis

To lower the cost of hydrogen produced by anion exchange membrane water electrolysis (AEMWE), it is critical to reduce the use of platinum group metal (PGM) catalysts within the device. While iridium has been successfully replaced with PGM-free catalysts at the anode, platinum-based (Pt) cathode catalysts are still required to meet the activity and durability targets. This study investigates the impact of commercial Pt/C catalyst loading, ionomer type and content, and electrode fabrication method on the cathode catalyst layer properties and AEMWE performance with the aim of determining the feasibility of reduced Pt loadings. While increased Pt loading is found to improve beginning-of-life performance, the effects are minimal above 0.6 mg/cm 2 . Ink characterization shows that ionomer type and content affect the ink stability, particle size, and percent of unbound ionomer, which further impact the homogeneity of the sprayed catalyst layers. The 5% PiperION cathode exhibited the highest performance, which may be attributed to a balance between the small particle size and the low proportion of unbound ionomer, minimizing kinetic and transport losses. Theoretical calculations show that the ionomers interact differently with the Pt surface, causing different surface charges and water adsorption strength and activating different mechanisms for hydrogen evolution. Pt-PiperION lowered the enthalpy of water-splitting by 0.1 eV compared to Pt alone and allowed for equal site access between adsorbed H* and OH* (both adsorbed at circa −2.2 eV). Although catalyst-coated membrane (CCM) fabrication techniques are desirable for scale-up, no performance enhancement is observed compared with the catalyst-coated substrate approach. Durability, as measured by degradation rates, Pt loss, and catalyst layer restructuring, was found to improve with increased Pt loadings, higher ionomer content, and CCM architectures. These findings provide important insight into the significant role of the cathode in AEMWE and strategies for maintaining the performance with low Pt loading or PGM-free catalysts.

08 HYDROGEN↗
DOE OSTI2024

Establishing robust correction schemes for improved and reliable ARM-AOS aerosol optical data products

Aerosol light absorption and scattering of solar radiation play an important role in the earth’s atmosphere in terms of direct and semi-direct radiative forcing. Optical parameters of importance to the US Department of Energy (DOE) climate models include absorption and scattering coefficients, single scattering albedo (SSA), absorption Angstrom exponents (AAE), and the asymmetry parameter (g). These parameters depend on aerosol size, shape and composition (refractive index), and are spectrally sensitive in the shortwave region. Additionally, these parameters have a complex dependency on the emission source, especially for carbonaceous aerosols. The DOE Atmospheric Radiation Measurement (ARM) user facility has deployed aerosol observing systems (AOS) containing several filter-based instruments to measure and constrain aerosol optical properties and related parameters at multiple sites worldwide. For measurement of aerosol light absorption, the AOS includes filter-based instruments (particle soot absorption photometer and tricolor absorption photometer) that infer particle-phase aerosol absorption coefficients at nominal red, green, and blue wavelength bands from the attenuation (ATN) of light passing through a particulate filter on which aerosols are deposited. Measurement of aerosol scattering is done in situ using nephelometers. By combining inferred absorption coefficients from filter-based ATN measurements and in situ scattering coefficients, value-added products (VAPs) such as SSA, AAE, and g are derived.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Case Study of the Glycoside Hydrolase Enzyme Mechanism Using an Automated QM-Cluster Model Building Toolkit

Glycoside hydrolase enzymes are important for hydrolyzing the β-1,4 glycosidic bond in polysaccharides for deconstruction of carbohydrates. The two-step retaining reaction mechanism of Glycoside Hydrolase Family 7 (GH7) was explored with different sized QM-cluster models built by the Residue Interaction Network ResidUe Selector (RINRUS) software using both the wild-type protein and its E217Q mutant. The first step is the glycosylation, in which the acidic residue 217 donates a proton to the glycosidic oxygen leading to bond cleavage. In the subsequent deglycosylation step, one water molecule migrates into the active site and attacks the anomeric carbon. Residue interaction-based QM-cluster models lead to reliable structural and energetic results for proposed glycoside hydrolase mechanisms. The free energies of activation for glycosylation in the largest QM-cluster models were predicted to be 19.5 and 31.4 kcal mol −1 for the wild-type protein and its E217Q mutant, which agree with experimental trends that mutation of the acidic residue Glu217 to Gln will slow down the reaction; and are higher in free energy than the deglycosylation transition states (13.8 and 25.5 kcal mol −1 for the wild-type protein and its mutant, respectively). For the mutated protein, glycosylation led to a low-energy product. This thermodynamic sink may correspond to the intermediate state which was isolated in the X-ray crystal structure. Hence, the glycosylation is validated to be the rate-limiting step in both the wild-type and mutated enzyme.

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

Heterogeneous multiphase flow properties of volcanic rocks and implications for noble gas transport from underground nuclear explosions

Of interest to the Underground Nuclear Explosion Signatures Experiment are patterns and timing of explosion-generated noble gases that reach the land surface. The impact of potentially simultaneous flow of water and gas on noble gas transport in heterogeneous fractured rock is a current scientific knowledge gap. This article presents field and laboratory data to constrain and justify a triple continua conceptual model with multimodal multiphase fluid flow constitutive equations that represents host rock matrix, natural fractures, and induced fractures from past underground nuclear explosions (UNEs) at Aqueduct and Pahute Mesas, Nevada National Security Site, Nevada, USA. Capillary pressure from mercury intrusion and direct air–water measurements on volcanic tuff core samples exhibit extreme spatial heterogeneity (i.e., variation over multiple orders of magnitude). Petrographic observations indicate that heterogeneity derives from multimodal pore structures in ash-flow tuff components and post-depositional alteration processes. Comparisons of pre- and post-UNE samples reveal different pore size distributions that are due in part to microfractures. Capillary pressure relationships require a multimodal van Genuchten (VG) constitutive model to best fit the data. Relative permeability estimations based on unimodal VG fits to capillary pressure can be different from those based on bimodal VG fits, implying the choice of unimodal vs. bimodal fits may greatly affect flow and transport predictions of noble gas signatures. The range in measured capillary pressure and predicted relative permeability curves for a given lithology and between lithologies highlights the need for future modeling to consider spatially distributed properties.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2026

Sparsity Applications for Gradient‐Based Optimization of Wind Farms

Optimizing wind farms is essential for designing efficient energy systems, especially as farms grow larger and span multiple sites. However, this optimization becomes increasingly challenging due to the rising computational cost associated with more turbines. Gradient‐based optimization methods scale better than gradient‐free approaches for large problems, but the most computationally expensive component remains the calculation of gradients for the objective function and constraint Jacobians. To address this, we propose leveraging sparsity to accelerate gradient evaluations and reduce the size of the constraint Jacobian. Wind farms naturally exhibit sparsity—many turbines do not influence each other under certain wind directions. However, unlike traditional sparse problems with fixed patterns, wind farm sparsity is dynamic, requiring new strategies to handle changing interactions efficiently. This paper presents a study of sparsity in wind farm optimization and introduces several methods to exploit it. These strategies are tested on multiple farms using the analytic Cumulative Curl model, with gradients computed via automatic differentiation (AD). The same sparsity‐aware techniques are also applicable to finite difference (FD) methods, where they can yield even greater speedups due to the high cost of directional evaluations. Results show that sparse methods achieve up to a 10x speedup with less than ± 5% variance in optimized wake losses compared to traditional methods. These findings suggest that sparsity‐aware optimization not only maintains solution quality but also scales efficiently with farm size, enabling more comprehensive design exploration at reduced computational cost.

17 WIND ENERGY↗
DOE OSTI2021

Polycarboxylate ester adsorption on cement grains: Influence of polydispersity

Although in academic and industrial practice polymers are typically described by their average molecular parameters, most often they display a polydispersity both in size and in properties. In this paper, we focus on the influence of polydispersity on the adsorption of polycarboxylate ester on cement grains and on its influence on the definition of a so-called affinity of the polymer to the surface. We first show that a constant fraction of a polydisperse polycarboxylate ester always adsorbs no matter the availability of adsorption sites. We moreover suggest from our measurements and from scaling laws from literature that the degree of polydispersity in the structure of polycarboxylate ester should affect adsorption at both low and intermediate dosages and hydrodynamic radii in solution. However, it shall not have a major influence on polymer layer thickness or saturation dosage.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Arp 58 and Arp 68: two M 51-type systems

ABSTRACT We study two M 51-type systems Arp 68 and Arp 58, which strongly differ by their stellar masses, gas content, and environment. Long-slit spectral observations obtained at the Russian 6-m telescope were used to trace the distributions of a line-of-sight (LOS) velocity and a gas-phase oxygen abundance along the spectral cuts. Two systems are compared by their observed properties. We found a very strong large-scale non-circular motion of gas in both systems and a kpc-size saw-edged velocity profile along the tidal spiral arm of Arp 68, probably caused by the gas outflow due to the stellar feedback. A deep decrease of LOS velocity is also found in the ‘hinge’ region in Arp 58, where the inner spiral arm transforms into the tidal spiral arm, which was predicted earlier for M 51-type galaxies. Local sites of star formation and the satellites are compared with the evolutionary models at the colour–colour diagrams. Unlike the spiral galaxy Arp 58, the main galaxy in Arp 68 system is experiencing an ongoing burst of star formation. Gas-phase metallicity estimates show that Arp 58 has a higher metal abundance and reveals a shallow negative radial gradient of the gas-phase oxygen abundance. The emission gas in Arp 68 has noticeably lower metallicity than it is expected for a given luminosity of this galaxy, which may be connected with its space position in the local void.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Biexciton Condensation in Electron-Hole-Doped Hubbard Bilayers: A Sign-Problem-Free Quantum Monte Carlo Study

The bilayer Hubbard model with electron-hole doping is an ideal platform to study excitonic orders due to suppressed recombination via spatial separation of electrons and holes. However, suffering from the sign problem, previous quantum Monte Carlo studies could not arrive at an unequivocal conclusion regarding the presence of phases with clear signatures of excitonic condensation in bilayer Hubbard models. Here, we develop a determinant quantum Monte Carlo algorithm for the bilayer Hubbard model that is sign-problem-free for equal and opposite doping in the two layers and study excitonic order and charge and spin density modulations as a function of chemical potential difference between the two layers, on-site Coulomb repulsion, and interlayer interaction. In the intermediate coupling regime and in proximity to the SU(4)-symmetric point, we find a biexcitonic condensate phase at finite electron-hole doping, as well as a competing (π, π) charge density wave state. In this study, we extract the Berezinskii-Kosterlitz-Thouless transition temperature from superfluid density and a finite-size scaling analysis of the correlation functions and explain our results in terms of an effective biexcitonic hard-core boson model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Observing Profiles of Derived Kinematic Field Quantities Using a Network of Profiling Sites

Abstract Observations of thermodynamic and kinematic parameters associated with derivatives of the thermodynamics and wind fields, namely, advection, vorticity, divergence, and deformation, can be obtained by applying Green’s theorem to a network of observing sites. The five nodes that comprise the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) profiling network, spaced 50–80 km apart, are used to obtain measurements of these parameters over a finite region. To demonstrate the applicability of this technique at this location, it is first applied to gridded model output from the High-Resolution Rapid Refresh (HRRR) numerical weather prediction model, using profiles from the locations of ARM network sites, so that values calculated from this method can be directly compared to finite difference calculations. Good agreement is found between both approaches as well as between the model and values calculated from the observations. Uncertainties for the observations are obtained via a Monte Carlo process in which the profiles are randomly perturbed in accordance with their known error characteristics. The existing size of the ARM network is well suited to capturing these parameters, with strong correlations to model values and smaller uncertainties than a more closely spaced network, yet it is small enough that it avoids the tendency for advection to go to zero over a large area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Characterizing warm atmospheric boundary layer over land by combining Raman and Doppler lidar measurements

PBL plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, for PBL characterization and process study, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary. Here, new approaches are developed to provide reliable PBLH and MLH to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) user facility. Compared with widely used lidar aerosol measurements for PBLH, WVMR is a better trace for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to re-construct 2-D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm seasons’ RL and ML measurements at the SGP site for five years to study warm season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable PBLH and MLH dataset will be valuable for PBL process study, model evolution, and PBL parameterization improvement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Connecting particle interactions to agglomerate morphology and rheology of boehmite nanocrystal suspensions

Rheology imposes significant challenges on processing of complex suspensions such as nuclear waste slurries at the Hanford and Savannah River sites. Understanding rheology connecting to microstructures and underlying particle interactions in complex slurries is therefore important for both fundamental knowledge and practical applications. Here, we use suspensions of aluminum oxyhydroxide minerals in the form of boehmite as an analog of the radioactive waste slurry to gain physical insights on the correlation between particle interactions, microstructures, and slurry rheology. Specifically, we use a combination of Couette rheometry and small-angle scattering techniques (independently and simultaneously) to understand how the slurry microstructure changes under flow and how these structural changes manifest themselves in the bulk rheology of the suspensions. Our experiments show that the boehmite slurries are thixotropic, with the rheology and structure of the suspensions changing with increasing exposure to flow. In the slurries, particle aggregates begin as loose, system-spanning clusters, but exposure to moderate shear rates causes the aggregates to irreversibly consolidate into denser clusters of finite size. The microstructural changes directly influence the rheological properties of the slurries such as viscosity and viscoelasticity. More importantly, our study shows that solution pH affects the amount of structural rearrangement and the kinetics of the rearrangement process, with an increase in pH leading to faster and more dramatic changes in the bulk rheology. Such dynamic microstructural changes and resultant rheology are understood via correlations between particle interactions and strength of particle network, coupled with surface chemistry and anisotropic nature of particle interactions. Nearly identical structural changes are also observed in Poiseuille flow geometries, implying that the observed changes are relevant in the pipe flow conditions present during waste processing.

Weston, Javen S.↗
DOE OSTI2024

Catalytic Consequences of Hierarchical Pore Architectures within MFI and FAU Zeolites for Polyethylene Conversion

The benefits of hierarchical zeolites for the conversion of bulky molecules like polymeric waste have been reported in the literature; however, the impact of mesopore sizes and connectivities on rates, product selectivities, and catalyst deactivation in the context of plastic upcycling has not been systematically probed. Here, in this study, we synthesized a suite of hierarchical MFI and FAU zeolites via desilication under varying conditions for metal-free polyethylene conversion reactions under batch and flow conditions (473–523 K). Polyethylene (solid) conversion rates (normalized by Bro̷nsted acid site density) were higher on hierarchical than parent microporous MFI regardless of mesopore connectivities, i.e., open or constricted, suggesting that the incorporation of mesopores facilitates diffusion of intermediate products to access medium-pore protons for successive scission events. Furthermore, higher branched:linear gaseous product ratios were produced on hierarchical than parent MFI, since mesopores allow for egress of bulkier molecules without undergoing further secondary events, e.g., isomerization back to linear alkanes/alkenes or beta scission. Solid conversion rates on hierarchical FAU synthesized via desilication with cetyltrimethylammonium bromide (CTABr), however, were not higher than parent FAU, likely because the presence of CTABr facilitates recrystallization of leached species to form composites (hierarchical FAU and ordered mesoporous materials) with more isolated mesopores. The stagnation in rates, despite increased mesopore volumes (>0.22 cm 3 g –1 ), highlights the importance of confinement effects provided by micropores for cleaving C–C bonds at modest reaction conditions. In situ 1 H MAS NMR performed on polyethylene with MFI zeolite show that PE isomerizes (and potentially deconstructs) at temperatures near 450 K, highlighting the role of Bro̷nsted acid sites in activating C–C bonds under mild reaction conditions. Catalyst recyclability studies showed that all catalysts undergo deactivation during plastic upcycling reactions, but to varying extents. Overall, hierarchical materials have better catalyst stability than parent materials, although the differences in stability between hierarchical and parent FAU are smaller than those for MFI. Taken together, these findings demonstrate how rates, selectivities, and catalyst deactivation from plastic upcycling reactions can be controlled via fine-tuning the identity and connectivity of mesopores.

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

Utah FORGE 4-2492: Design and Implementation of Innovative Stimulation Treatments to Maximize Energy Recovery - Workshop Presentation

This is a presentation on the Design and Implementation of Innovative Stimulation Treatments to Maximize Energy Recovery Efficiency at the Utah FORGE Site project by The University of Texas at Austin, presented by Professor Mukul M. Sharma. The project's objectives were to place fractures uniformly in a horizontal well (improve cluster efficiency) to ensure uniform distribution of flow into all hydraulic fractures, maximize the area of the created fracture network, ensure connectivity of the fractures from the injector to the producer, and optimize fracture size. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 7, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗
DOE OSTI2024

Photovoltaic and Cost Analysis for Winston-Salem, North Carolina

This study assesses the feasibility of solar installations at various sites in Winston-Salem, focusing on factors such as solar resource availability, electricity costs, and rooftop area for photovoltaic systems. The analysis begins by estimating daily energy requirements based on annual electricity use, followed by adjusting for seasonal and operational fluctuations. Using regional solar data from the National Solar Radiation Database (NSRDB), we calculate average peak sun hours to determine effective system sizing. A parametric approach using the System Advisor Model (SAM) refines this sizing process, incorporating a safety margin of 1.2 to address demand peaks. Each PV system is designed to fit available rooftop space, and a coverage threshold of 70% is identified as optimal for maximizing cost savings and sustainability. This threshold allows installations to meet substantial energy demands, supporting energy resilience and enhancing economic returns.

14 SOLAR ENERGY↗
DOE OSTI2025

ARM SGP PBLH and MLH datasets from Raman lidar and Doppler lidar

The planetary boundary layer (PBL) plays a critical role in the atmosphere by transferring heat, moisture, and momentum. The warm PBL has a distinct diurnal cycle including the daytime convective mixing layer (ML) and nighttime residual layer developments. Thus, simultaneous determinations of PBL height (PBLH) and ML height (MLH) are necessary for studying PBL characterization and processes. Here, new approaches are developed to provide reliable PBLH and MLH estimates to characterize warm PBL evolution. The approaches use Raman lidar (RL) water vapor mixing ratio (WVMR) and Doppler lidar (DL) vertical velocity measurements at the Southern Great Plains (SGP) atmospheric observatory, which was established by the Atmospheric Radiation Measurement (ARM) User Facility. Compared to widely used lidar aerosol measurements for PBLH, WVMR is a better tracer for PBL vertical mixing. For PBLH, the approach classifies PBL water vapor structures into a few general patterns, then uses a slope method and dynamic threshold method to determine PBLH. For MLH, wavelet analysis is used to reconstruct 2D variance from DL vertical wind velocity measurements according to the turbulence eddy size to minimize the impacts of gravity wave and eddy size on variance calculations; then, a dynamic threshold method is used to determine MLH. Remotely-sensed PBLHs and MLHs are compared with radiosonde measurements based on the Richardson number method. Good agreements between them confirm that the proposed new algorithms are reliable for PBLH and MLH characterization. The algorithms are applied to warm-season RL and ML measurements at the SGP site for five years to study warm-season PBL structure and processes. The weekly composited diurnal evolutions of PBLHs and MLHs in a warm climate were provided to illustrate diurnal and seasonal PBL evolutions. This reliable data set of PBLH and MLH values will be valuable for studying PBL processes, model evolution, and PBL parameterization improvements. The MLH dataset includes the MLH in values of km above ground level. The PBLH dataset includes the PBLH in values of km above ground level, along with a flag ("situation_PBLH") to determine the state of the PBL (1 = Cloudy Condition, 2 = Stable Layer, 3 = Multi-layer WVMR structure, 4 = Well-Mixed PBL, 5 = A de-coupled layer, 6 = Other).

mixing layer height↗
DOE OSTI2020

Behind the Meter Storage for Electric Vehicle Charging, Electrochemical and Thermal Energy Storage, and Solar Photovoltaic

In response to the potentially large and irregular demand from EVs, along with changing load profiles from buildings with on-site generation, utilities are evaluating multiple options for managing dynamic loads, including time-of-use pricing, demand charges, battery storage, and curtailment of variable generation. Buildings, as well as commercial, public, and workplace EV charging operations, can use a combination of electrochemical battery storage and thermal energy storage coupled with on-site generation to manage energy costs as well as provide resiliency and reliability for EV charging and building energy loads. We are completing a behind the meter storage analysis that focuses on determining the optimal system designs and energy flows for thermal and electrochemical behind the meter storage with on-site solar photovoltaic (PV) generation enabling electric vehicle charging in various climates, building types, and utility rate structures. In completing this analysis, we have developed a tool that combines existing battery models via the System Advisor Model (SAM) and building modeling software via EnergyPlus into a single interface. This tool allows us to simulate a building with a detailed battery model to properly size the battery, thermal energy storage, and solar PV systems to maximize profit for the system owner. This also allows us to assess how the battery degrades under various supervisory control dispatch algorithms to control charging/discharging; we can also see how thermal energy storage is created and used to complement the battery to reduce thermal loads in the building. With this project, we can analyze new batteries that are designed specifically for energy storage, rather than designed to be extremely energy dense for electric vehicle applications, using battery lifetime models from other national labs and the existing SAM battery model, which has detailed lifetime and degradation parameters. We can also assess novel thermal storage technologies by integrating them into the whole building energy simulation program EnergyPlus. Because the model calls both SAM and EnergyPlus, required inputs need to be compatible for both models. These inputs include, on a high-level, the following: weather files, building and electric vehicle load profiles, electricity rate tariff information, and system cost information for the stationary battery, solar PV, and thermal storage system. The various buildings we are studying for this analysis are retail big-box grocery store, commercial office building, fleet vehicle depot and operations facility, multi-family residential, and electric vehicle charging station. For these different applications, the battery and thermal storage will be dispatched differently, and the various technologies are sized differently to optimize cost.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Land use-land cover gradient demonstrates the importance of perennial grasslands with intact soils for building soil carbon in the fertile Mollisols of the North Central US

The impact of land use change and agricultural management on the cycling of soil organic carbon (SOC) is not well understood, limiting our ability to manage for, and accurately model, soil carbon changes at both local and regional scales. To address this issue, here we combined long-term soil incubations with acid-hydrolysis and dry combustion to parse total SOC (C t ) into three operationally defined SOC pools (active, slow, and recalcitrant) from 9 long-term sites with varying land uses on current and former tallgrass prairie soil. Land uses represented a gradient of soil disturbance histories including remnant prairie, restored prairie, grazed pasture, annual crop rotations, and continuous maize. Dry combustion was used to estimate total carbon (C t , physical), while acid hydrolysis of both the active (C a ) and slow (C s ) pools was used to estimate a recalcitrant carbon pool (C r , chemical). Non-linear modeling of CO 2 efflux data from the long-term incubations was then used to estimate C a , and the decomposition rates of both C a and C s (k a and k r , biological). The size of the slow pools C s was then defined mathematically as C t- (C a + C r ). Remnant prairie had the highest C t , while cool-season pasture and a 35-y-old restored prairie had higher C t than the other agricultural systems. All agricultural systems, including pasture, had the highest fraction of C t as C r (~50%), whose mean residence time (MRT) in these soils is ≥500 years (Paul et al., 2001a) demonstrating that this fraction persists, while the more labile fractions were lost over the course of a few months (C a ) to a few decades (C s ) as a result of tillage-intensive agriculture. The two- to four-decade MRT time of C s indicated a pool likely to be more responsive to the 20 to 40 years of land-use practices used at some of the sites. The C s pool was largest in the remnant- and 35-y-old prairies indicating significant C accrual and stabilization compared to the agricultural ecosystems. Interestingly, the remnant prairie maintained the highest C a pool as well, demonstrating the strong connection between the quantity of fresh C inputs and the potential for long-term C stabilization and accrual. The accumulation of C in active (≈labile) pools as a first step toward long-term stabilization highlights the tenuous nature of early carbon gains, which can be quickly lost in response to climate change or poor management.

54 ENVIRONMENTAL SCIENCES↗