SEARCH · Search NASA

Results for “critical zone”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Selective penetration behavior of microgels in superpermeable channels and reservoir matrices

Gel treatment is an effective way to attack excessive water production in many mature oilfields around the world. Selective penetration is desired for successful gel treatments. That is, gel materials should easily penetrate the target zones (i.e., channeling features such as superpermeable channels) without entering/damaging the nontarget zones (i.e., reservoir matrices or oil zones). This study revealed that presence of threshold penetration pressure (ΔP th ) was responsible for selective penetration behavior of tested microgels. The concept of ΔP th was utilized to figure out favorable working conditions for effective gel treatments. Microgel dispersions were injected into superpermeable (super-k) sandpacks (mimicking super-k channels in reservoirs, 60–221 darcies), heterogeneous models with super-k channels (79–230 darcies), and sandstone cores (mimicking reservoir matrices, 50–5000 md). The results demonstrated that a minimum differential driving pressure (i.e., threshold penetration pressure, ΔP th ) was required to push microgel particles to penetrate channels or matrices. The critical penetration behavior was closely related to the particle/pore size ratio. Low ΔP th at smaller particle/pore ratios was beneficial to allow easy penetration of gel materials into the channeling zones. On the contrary, high ΔPth at larger particle/pore ratios was desirable to prevent gel materials from massively invading and damaging the matrices. Instead, the gel particles accumulated at the inlet surface, and a gel cake was gradually formed. The cake further prevented the invasion of the gels. The cake could be removed by chemical breakers to resume the injectivity/productivity of the matrices. Correlations were developed to describe the relationship between ΔP th and particle/pore ratio. A distinct transition was identified at the particle/pore ratio of about 3. This research could help identify the favorable conditions to achieve successful gel treatments. In an effective conformance treatment, the particle/pore ratio in the channel should be sufficiently low to allow easy penetration of gel materials into the channel (e.g., particle/pore ratio<2 in this study). Meanwhile, the particle/pore ratio in the matrix should be large enough to support a high ΔP th and thus prevent massive gel invasion into the matrix. This study advances the current pore scale studies (a single particle passing through a single channel) to Darcy-scale characterization.

02 PETROLEUM↗

DOE OSTI2025

Chemical Reactor Network Modeling of Ammonia Rich-Quench-Lean Combustion Using a Partially Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However, low flammability and a propensity for high nitrogen oxide (NO x ) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NO x emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NO x concentrations below those of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing nonadiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and postflame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors' past work, this arrangement has demonstrated the significant sensitivity of exit NO x to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and nonpremixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NO x emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have been shown to create conditions that are conducive to increased N 2 O formation in the lean stage. In conclusion, the results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia.

Combustion↗

DOE OSTI2025

CRN Modeling of Ammonia RQL Combustion using a Partially-Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However low flammability and a propensity for high nitrogen oxide (NOx) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NOx emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel, while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NOx concentrations below that of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing non-adiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and post-flame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors’ past work, this arrangement has demonstrated the significant sensitivity of exit NOx to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and non-premixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NOx emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have shown to create conditions which are conducive to increased N2O formation in the lean stage. The results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia..

advanced gas turbines↗

DOE OSTI2022

Coastal Acoustic Buoy for Offshore Wind: Project Synthesis

The population of North Atlantic right whales is critically endangered and their habitat overlaps with offshore windfarm leases. It is therefore imperative that effective mitigation strategies be used to avoid impacts on right whales during the construction of offshore windfarms. The Department of Energy issued FOA Number DE-FOA-0001924 to encourage the development of technology that could monitor large exclusion zones for right whales in order to mitigate potential impact of construction noise on right whales. This report summarizes the past two years of the development and evaluation of the Coastal Acoustic Buoy for Offshore Wind (CABOW) project which aimed to develop technology to monitor large exclusion zones for North Atlantic right whales. Over the course of the project SMRU Consulting have implemented a rigorous design process including comparison of different approaches (e.g., single sensor vs multiple sensors), as well as consideration of placement and timing of acoustic monitoring. We have evaluated critical components of the CABOW system including, reliability, detection range, and bearing accuracy in areas adjacent to offshore windfarm leases in Maryland by conducting 3,536 playbacks of simulated right whale upcalls. The maximum call detection range was 7.5 km when noise was 99 dB re 1µPa rms (50-225 Hz), but this reduced to < 1 km when ambient noise levels were high. Our detection probability in the field was measured as a function of range as well as the source-to-noise level ratio allowing us to build a model to predict the probability of detection under various scenarios (sample size: 3,536 calls x 5 buoys = 17,680). The median bearing error was -0.25° but this is likely an underestimate of error due to experimental design. Using the published recall and precision of the two detectors we implemented in the CABOW system, we estimate that at a recall of 80%, our precision was > 80%, within the range of what we were aiming for in this project. To estimate our exclusion zone false negative and false positive rates, we built a simulation model using the empirical data from our field trial. We modelled three to nine CABOW units placed on the 10 km exclusion zone and estimated our false negative rate to be 1% or less (which was our project goal) and our false positive rate to be between seven and nine percent, slightly above our goal of 5%. However, we also modelled an equivalent PAM system that does not have bearing capabilities and found the false positive rates for that system to be six to eight times higher than the CABOW rate. This higher false positive rate of PAM systems with low spatial information could have significant cost repercussions for offshore wind developers by adding work shutdowns or delays without providing additional protection for right whales. The model we built allows us to explore the placement of PAM systems under various scenarios and will thus help facilitate planning of PAM mitigation systems to meet NOAA Incidental Harassment Authorizations for specific windfarms. We achieved an average system uptime of 98.3%, just below our goal of 99%. The issues that caused these short losses of data have been identified and fixed. Right whale detections and audio clips were typically transferred via radio from the buoys to the base station in two to four seconds. We therefore believe we have developed a highly robust real-time PAM system. Based on the above, we feel we have achieved the stated funding goal of developing a cost-effective and robust real-time PAM system that enables the monitoring of large exclusion zones for right whales during the constructions of offshore windfarms. This should lead to decreased costs and risks for the offshore wind sector while providing robust mitigation for right whales. It is important to state that PAM mitigation will need to be implemented with other mitigation strategies (e.g., visual observers) to provide a complete mitigation strategy to ensure that any effects on right whales from offshore windfarm construction is minimized.

17 WIND ENERGY↗

DOE OSTI2026

Physics vs structure: A systematic benchmark of learning strategies for multi-zone building thermal dynamics

Recent advances in physics-informed and data-driven machine learning promise improved thermal models for advanced building control, yet there is limited quantitative evidence on when added physics structure and architectural complexity are beneficial. Here, this work presents a systematic benchmark of five representative system identification methods for modeling multi-zone building thermal dynamics: linear state-space models, multi-layer perceptrons, neural state-space models, neural ordinary differential equations, and physically-consistent neural networks. The methods are evaluated across multiple data regimes and zone coupling strategies. Using a high-fidelity multi-zone commercial building emulator, we examine short-term and long-term prediction accuracy, computational efficiency, and ease of development. Our results reveal critical trade-offs between prediction performance, model complexity, and physical consistency. We demonstrate that decoupled, nonlinear black-box models consistently outperform coupled physics-constrained architectures in both predictive accuracy and out-of-distribution robustness in majority of the test cases for the building type considered in the study. Our findings quantify the cost of complexity in building thermal modeling and provide concrete, actionable, scenario-based guidelines for selecting model classes for control-oriented applications.

Building thermal modeling↗

DOE OSTI2022

Interactions between molecular-scale processes and hyporheic exchange for understanding Fe-S-C cycling in riparian wetlands (Final Report)

Wetlands represent some of the most productive ecosystems on the planet and critically influence global environmental health. Specifically, wetlands promote water quality by transforming nutrients and organic compounds and sequestering metals and contaminants. Riparian wetland hyporheic zones, where toxic surface water and anoxic groundwater mix, exhibit dynamic conditions that drive steep redox gradients and promote hotspots of diverse and fluctuating microbial activity. Changes in climate, water quality, and water quantity can disturb hydrologic flow and biogeochemical processing in these environments. Understanding how sulfate loading, impacted by hydrologic flux and anthropogenic inputs, influences iron and carbon cycling in wetlands will be crucial for predicting water quality issues driven by iron mineral precipitation and sorption, such as the release of heavy metals and other toxic elements. We used a fully integrated multi-scale and multi-method approach to develop a mechanistic understanding of how hydrologic flow influences coupled iron and sulfur cycles in riparian wetlands. This entailed hydrological, geochemical, and microbial observations at two locations: an anthropogenic sulfate-impacted riparian wetland in northern Minnesota, and a low-sulfate Fe-rich riparian wetland in Tims Branch at the Savannah River Site (SRS). These sites were characterized by hydrologically dynamic conditions where the stream and wetland systems oscillated between gaining (upward flow) and losing (downward flow) conditions that recharged the system occasionally with oxidants that fueled a variety of biogeochemical reactions. Aqueous geochemical measurements of surface water, groundwater, and porewater samples were made alongside solid-phase geochemical analyses of sediment gravity cores. Bulk X-ray absorption spectroscopy at the Advanced Photon Source (APS; Argonne) interrogated the speciation and distribution of Fe and S mineral phases of the sediments. Interestingly, it was discovered that, despite strongly reducing conditions, Fe(III) compounds and a variety of intermediate valence S compounds were stable in the subsurface. Indeed, compounds like thiosulfate, S(0), and intermediate valence organosulfur compounds were more prevalent than FeS and pyrite. The composition of the sediments did change with changes in hydrologic flow, showing their reactivity in changing redox conditions. The abundance of these intermediate S compounds was likely formed as a result of anaerobic oxidation by aqueous and solid-phase Fe(III) compounds, fueling a cryptic S cycle that is driving the breakdown of organic matter in the hyporheic zone. Microbiome surveys showed a core community that seemed stable across the landscape, but changed with increasing depth into the sediment. The community composition did not seem to change dramatically with changes in season or hydrologic flow, except for some organisms that have the potential to contribute to S cycling. More work is needed to confirm their functional activity. These fine process-scale analyses were placed within a dynamic field context using physical flow parameters from surface water and groundwater level measurements. These data helped shed light on sulfur-driven biogeochemical processes in hydrologically dynamic riparian wetlands, addressing a gap in our understanding about the impacts of pollution and other anthropogenic changes on ecologically sensitive environments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Application of New Los Alamos OPLIB Opacities in Solar Modeling Using the Mesa Code

The Sun provides the only local laboratory to test our understanding of stellar physics. The production, then, of valid dynamical models of the Sun is critical to how we think about other stars. Despite improvements in modeling methodologies and abundance measurements, solar models continue to reproduce errors in convection zone depth and surface helium abundance compared to helioseismic observations. Using the 1D stellar evolution code MESAstar, this work first evaluates the effect of changes to timestep resolution and mixing length theory on model outputs. A standard test case is then constructed with both GS98 and AGSS09 abundances using OP and OPAL opacities before implementing the new Los Alamos OPLIB opacities. It is shown that, when using GS98 abundances, OPLIB opacities improve agreement with helioseismic inference of the convection zone base radius, R czb , by 1.68-σ and 1.19-σ compared to OP and OPAL, respectively, while worsening agreement with the inferred envelope helium mass fraction, Y surf , by 0.63-σ and 0.98-σ for OP and OPAL.

14 SOLAR ENERGY↗

DOE OSTI2023

Worth More Dead than Alive? Quantifying Necromass Persistence for Terrestrial Carbon Storage

The continuous cycle of plant-associated microbial biomass growth, death, and decay provides a consistent supply of necromass-C throughout the rooting zone of the soil profile. These inputs are obviously important for resource allocation and soil quality, but microbial necromass-C may play a critical, but currently unknown, role in mitigating the impacts of climate change through atmospheric decarbonization. Microbial necromass is the largest terrestrial sink of persistent carbon in soils, thus its fate and preservation will have a major impact on global C budgets. Our current understanding of the ecosystem controls on necromass generation, biogeochemical transformation, and stabilization is lacking. The objective of this project is to explore the specific influence of rhizosphere and necromass inputs on soil carbon pools.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Results of ISO/TS 6336-22 Evaluating Full Contact Zone

ISO/TS 6336-22 (Calculation of load capacity of spur and helical gears - Part 22: Calculation of micropitting load capacity) is the ISO technical specification containing a proposal for calculations of the risk of micropitting in gear sets. Micropitting is a Hertzian fatigue phenomenon that appears as ultra-fine cracking and pitting on the flanks of gear teeth. Since progressive micropitting can lead to macropitting and flank damage, critical applications such as wind turbines, marine drives, and high-speed gear drives seek to accurately predict whether their designs are susceptible to this damage. ISO/TS 6336-22 assesses micropitting risk through a safety factor which is calculated as the minimum specific film thickness in the contact zone divided by a permissible specific film thickness. In the previous paper, the calculations were performed using the simplified method (Method B) that evaluates points on the path of contact. This was done for three gear sets that experienced micropitting in operation. The minimum specific film thickness for the two field cases was very high, which indicates that the gears were operating in the full elastohydrodynamic lubrication (EHL) regime. A more accurate calculation for these cases (Method A) calculates the specific film thickness across the entire contact zone. This paper applies this Method to the case study from the previous paper. The results are compared to micropitting observed in operation. Results are also compared to the results of the previous paper. Conclusions are made regarding the accuracy of both Methods compared to the field cases and relative to each other.

ENGINEERING,WIND ENERGY↗

DOE OSTI2022

Results of ISO/TS 6336-22 Evaluating Full Contact Zone

gearbox↗

DOE OSTI2022

Results of ISO/TS 6336-22 Evaluating Full Contact Zone: Preprint

ISO/TS 6336-22 (Calculation of load capacity of spur and helical gears — Part 22: Calculation of micropitting load capacity) is the ISO technical specification containing a proposal for calculations of the risk of micropitting in gear sets. Micropitting is a Hertzian fatigue phenomenon that appears as ultra-fine cracking and pitting on the flanks of gear teeth. Since progressive micropitting can lead to macropitting and flank damage, critical applications such as wind turbines, marine drives, and high-speed gear drives seek to accurately predict whether their designs are susceptible to this damage. ISO/TS 6336-22 assesses micropitting risk through a safety factor which is calculated as the minimum specific film thickness in the contact zone divided by a permissible specific film thickness. In the previous paper, the calculations were performed using the simplified method (Method B) that evaluates points on the path of contact. This was done for three gear sets that experienced micropitting in operation. The minimum specific film thickness for the two field cases was very high, which indicates that the gears were operating in the full elastohydrodynamic lubrication (EHL) regime. A more accurate calculation for these cases (Method A) calculates the specific film thickness across the entire contact zone. This paper applies this Method to the case study from the previous paper. The results are compared to micropitting observed in operation. Results are also compared to the results of the previous paper. Conclusions are made regarding the accuracy of both Methods compared to the field cases and relative to each other.

gearbox↗

DOE OSTI2022

Direct visualization of polaron formation in the thermoelectric SnSe

SnSe is a layered material that currently holds the record for bulk thermoelectric efficiency. The primary determinant of this high efficiency is thought to be the anomalously low thermal conductivity resulting from strong anharmonic coupling within the phonon system. Here we show that the nature of the carrier system in SnSe is also determined by strong coupling to phonons by directly visualizing polaron formation in the material. We employ ultrafast electron diffraction and diffuse scattering to track the response of phonons in both momentum and time to the photodoping of free carriers across the bandgap, observing the bimodal and anisotropic lattice distortions that drive carrier localization. Relatively large (18.7 Å), quasi-one-dimensional (1D) polarons are formed on the 300-fs timescale with smaller (4.2 Å) 3D polarons taking an order of magnitude longer (4 ps) to form. This difference appears to be a consequence of the profoundly anisotropic electron–phonon coupling in SnSe, with strong Fröhlich coupling only to zone-center polar optical phonons. These results demonstrate a high density of polarons in SnSe at optimal doping levels. Strong electron-phonon coupling is critical to the thermoelectric performance of this benchmark material and, potentially, high performance thermoelectrics more generally.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Vertical Characterization of Flow and Contaminant Concentration in New Hanford Extraction Wells

Effective groundwater remediation should focus on high-concentration, low-permeability zones within the aquifer, as these areas often present significant challenges to contaminant removal. Characterizing the vertical distribution of contamination in an aquifer is therefore critical for optimizing and monitoring pump-and-treat (P&T) systems to meet target cleanup timelines. Flow patterns driven by hydraulic gradients and variations in permeability are described to provide insight into contaminant mass removal from long-screened wells. The objective of this multi-phase work is to provide actionable recommendations for improving remedy performance monitoring and P&T operations, which ongoing evaluations will address. This report summarizes the preliminary results of vertical flow and contaminant distribution monitoring in new extraction wells at the Hanford Site.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Reconciling Contrasting Relationships Between Relative Dispersion and Volume-Mean Radius of Cloud Droplet Size Distributions

Cloud droplet spectral relative dispersion is critical to parameterizations of cloud radiative properties, warm-rain initiation, and aerosol-cloud interactions in models; however, there is no consistent relationship between relative dispersion and volume-mean radius in literature, which hinders improving relative dispersion parameterization and calls for physical explanation. Here we show, by analyzing aircraft observations of cumulus clouds during Routine AAF [Atmospheric Radiation Measurement (ARM) Aerial Facility] Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations, that the correlation between relative dispersion and volume-mean radius changes from positive to negative as volume-mean radius increases. With the new observation, we postulate that the sign of the correlation is determined by whether or not condensation (evaporation) occurs simultaneously with significant new activation (deactivation). The hypothesis is validated by simulations of both an adiabatic cloud parcel model and a parcel model accounting for entrainment-mixing. A new quantity, first bin strength, is introduced to quantify this new observation. Theoretical analysis of truncated gamma and modified gamma size distributions further supports the hypothesis and reconciles the contrasting relationships between relative dispersion and volume-mean radius, including the results in polluted fog observations. The results could shed new light on the so-called “twilight zone” between cloudy and cloud-free air, which in turn affects evaluation of aerosol-cloud interactions and retrieval of aerosol optical depth.

Lu, Chunsong↗

DOE OSTI2026

Predicting multiphase flow and tracer transport for an underground chemical explosive test

Detecting radionuclide gas seepage from clandestine underground nuclear tests is central to nonproliferation explosion monitoring research. Yet, early-time (<6 day) gas transport driven by the explosive pressure wave remains poorly constrained due to scarcity of field data. We simulate multi-phase gas transport in the vadose zone using pre-shot data from a recent chemical explosion in P-Tunnel at the Nevada National Security Site, USA. Despite using a simplified 2D-radial model, predictions of tracer arrival matched observations within one order-of-magnitude. Our results show how transient blast forcing rapidly mobilizes gases from the cavity into surrounding rock – critical for optimizing sensor placement and test planning. This unique integration of field data and modeling represents a significant improvement in our ability to predict gas migration from underground explosions. More broadly, it offers insights into the coupled dynamics of pressure waves and contaminant transport in the vadose zone, with implications for monitoring and hazard assessment.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Development of a Full-Scale Connected U-Net for Reflectivity Inpainting in Spaceborne Radar Blind Zones

CloudSat’s Cloud Profiling Radar is a valuable tool for remotely monitoring high-latitude snowfall, but its ability to observe hydrometeor activity near the Earth’s surface is limited by a radar blind zone caused by ground clutter contamination. This study presents the development of a deeply supervised U-Net-style convolutional neural network to predict cold season reflectivity profiles within the blind zone at two Arctic locations. The network learns to predict the presence and intensity of near-surface hydrometeors by coupling latent features encoded in blind zone-aloft clouds with additional context from collocated atmospheric state variables (i.e., temperature, specific humidity, and wind speed). Results show that the U-Net predictions outperform traditional linear extrapolation methods, with low mean absolute error, a 38% higher Sørensen–Dice coefficient, and vertical reflectivity distributions 60% closer to observed values. The U-Net is also able to detect the presence of near-surface cloud with a critical success index (CSI) of 72% and cases of shallow cumuliform snowfall and virga with 18% higher CSI values compared to linear methods. An explainability analysis shows that reflectivity information throughout the scene, especially at cloud edges and at the 1.2-km blind zone threshold, along with atmospheric state variables near the tropopause, are the most significant contributors to model skill. This surface-trained generative inpainting technique has the potential to enhance current and future remote sensing precipitation missions by providing a better understanding of the nonlinear relationship between blind zone reflectivity values and the surrounding atmospheric state.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

A continuum model for heat and mass transfer in moving-bed reactors for thermochemical energy storage

In this work, a continuum heat and mass transfer model coupling transport phenomena and high-temperature thermochemical reactions is developed for stationary packed-bed and counter-flow moving-bed reactors. After presenting the general modeling framework, we focus on the 2D axisymmetric version of the model for which validation is conducted with experimental results for a packed-bed reactor in the literature for manganese-iron oxide reduction/oxidation and an in-house counter-flow moving-bed reactor for magnesium-manganese oxide reduction up to 1450 °C. Transient simulation results including the local distributions of gas/solid temperatures, oxygen concentration and the extent of reaction, as well as the various energy flow components and energy conversion efficiencies are reported. The results based on the 2D axisymmetric model are also compared with those obtained from a previous 1D model. The comparison shows that capturing the radial variation is critical in reactor modeling and the 2D results demonstrate improved agreement with experiments. Specifically, large temperature variations along the radial direction are observed especially in the reaction zone; this non-uniform radial temperature distribution has a significant effect on the chemical reaction extent due to its strong dependence on temperature; and the overall oxygen concentration at the reactor exit and the predicted system efficiency are slightly lower in the 2D model compared to the 1D model. Finally, the present heat and mass transfer model can provide valuable insights into reactor design, scale-up, and operating conditions selection to maximize system energy storage efficiency.

25 ENERGY STORAGE↗

DOE OSTI2025

EnergyPlus-MCP: A model-context-protocol server for ai-driven building energy modeling

Traditional building energy modeling with the EnergyPlus building performance simulation engine requires domain expertise, programming skills, and intensive manual efforts limiting its effective adoption. This paper introduces EnergyPlus-MCP, the first open-source Model Context Protocol (MCP) server specifically designed for EnergyPlus simulation workflows, establishing a new foundational infrastructure for AI-driven building energy modeling. The MCP server implements a layered architecture with 35 specialized tools spanning model management, editing and analysis, HVAC and other systems configuration inspection, and simulation execution, enabling Large Language Models to interact with EnergyPlus through conversational interfaces. The server addresses critical workflow barriers by automating model validation, streamlining energy efficiency measures modification, and providing intelligent output management with interactive visualization. Through practical demonstrations using a multi-zone building retrofit analysis, we show how the EnergyPlus-MCP server significantly reduces manual efforts while maintaining full simulation rigor. By providing accessible natural language interfaces to sophisticated building energy analysis, this approach enables scalable deployment of simulation expertise across public and private organizations, educational institutions, and research teams, fundamentally transforming traditional building energy modeling practices.

AI↗