Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface operations analysis”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Computational Analysis of the Effect of Structured Packing Design on Absorption Column Hydrodynamics for Post-Combustion Carbon Capture Applications

Solvent based post-combustion carbon capture technologies have a potential for reducing carbon emissions from fossil-fuel-fired power plants and industrial sources where CO2 emissions are inherently harder to mitigate, such as steel or cement industries. A prominent technology to achieve this is by retrofitting absorption columns to the existing infrastructure. While these systems have been among the less costly alternatives for carbon capture, they still impose a considerable energy penalty to the operation of power plants or industrial facilities. The optimization of CO2 capture rate in solvent-based absorption process is complex as it depends on several factors including CO2 solubility, solvent reaction kinetics and temperature effects on the solubility, reaction rates, surface tension, and thermophysical properties of the solvent and the flue gas. The overall heat and mass transfer also depends on the hydrodynamics, which in turn, is affected by the packing geometry. In the current work, we systematically quantify the effects the design of the structured packing has on the column hydrodynamics, by performing detailed CFD simulations for different geometrical configurations and operating conditions. We then obtain relationships between the key hydrodynamic metrics, such as liquid holdup, interfacial area, wetted area, and pressure drop to the parameters defining the packing geometries and identify new more effective packing designs for the given operating conditions.

Shah, Yash Girish↗
DOE OSTI2023

Calendar life of lithium metal batteries: Accelerated aging and failure analysis

Lithium metal batteries (LMBs) are prime candidates for next-generation energy storage devices characterized by their remarkable energy storage capabilities. Despite the importance of understating the calendar aging impacts in LMBs, cycle life and calendar life have received inconsistent attention across various research and development stages. For acceptance into an application, especially electric vehicles, batteries are required to have sufficient calendar life to ensure performance over the life of vehicles which often experience extended periods of inactivity. In this study, an in-depth exploration into the calendar aging of LMB (Li/ Li[Ni 0.8 Mn 0.1 Co 0.1 ]O 2 in pouch cell format) is conducted. Specifically, the impact of calendar aging factors, state-of-charge (SOC), temperature, pressure, and operating conditions (open-circuit voltage, OCV and constant voltage, CV) on the surface morphology and thickness change rate in lithium metal anode is investigated. Notably, the results show that storing LMB pouch cells at an OCV, 80 % SOC, 25 °C, and 10 psi external pressure leads to a <1% reduction in capacity over 18 months. The results also suggest routes which can be utilized in the future for accelerated aging without altering the primary degradation modes. The results also show an accelerated calendar aging at high temperatures (45 °C), 100 % SOC, and CV suggesting they may serve as key accelerated aging methods for future studies. During calendar aging, pressure mitigates dendrite growth and extends calendar life. In conclusion, the insights from this study reinforce the viability of LMBs as a compelling solution for addressing increased energy storage needs in real-world applications.

25 ENERGY STORAGE↗
DOE OSTI2024

GCAM–GLORY v1.0: representing global reservoir water storage in a multi-sector human–Earth system model

Abstract. Reservoirs play a significant role in modifying the spatiotemporal availability of surface water to meet multi-sector human demands, despite representing a relatively small fraction of the global water budget. Yet the integrated modeling frameworks that explore the interactions among climate, land, energy, water, and socioeconomic systems at a global scale often contain limited representations of water storage dynamics that incorporate feedbacks from other systems. In this study, we implement a representation of water storage in the Global Change Analysis Model (GCAM) to enable the exploration of the future role (e.g., expansion) of reservoir water storage globally in meeting demands for, and evolving in response to interactions with, the climate, land, and energy systems. GCAM represents 235 global water basins, operates at 5-year time steps, and uses supply curves to capture economic competition among renewable water (now including reservoirs), non-renewable groundwater, and desalination. Our approach consists of developing the GLObal Reservoir Yield (GLORY) model, which uses a linear programming (LP)-based optimization algorithm and dynamically linking GLORY with GCAM. The new coupled GCAM–GLORY approach improves the representation of reservoir water storage in GCAM in several ways. First, the GLORY model identifies the cost of supplying increasing levels of water supply from reservoir storage by considering regional physical and economic factors, such as evolving monthly reservoir inflows and demands, and the leveled cost of constructing additional reservoir storage capacity. Second, by passing those costs to GCAM, GLORY enables the exploration of future regional reservoir expansion pathways and their response to climate and socioeconomic drivers. To guide the model toward reasonable reservoir expansion pathways, GLORY applies a diverse array of feasibility constraints related to protected land, population, water sources, and cropland. Finally, the GLORY–GCAM feedback loop allows evolving water demands from GCAM to inform GLORY, resulting in an updated supply curve at each time step, thus enabling GCAM to establish a more meaningful economic value of water. This study improves our understanding of the sensitivity of reservoir water supply to multiple physical and economic dimensions, such as sub-annual variations in climate conditions and human water demands, especially for basins experiencing socioeconomic droughts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Kinetic Modeling Analysis of Ar Addition to Atmospheric Pressure N 2 –H 2 Plasma for Plasma-Assisted Catalytic Synthesis of NH 3

Zero-dimensional kinetic modeling of atmospheric pressure Ar–N 2 –H 2 nonthermal plasma was carried out to gain mechanistic insights into plasma-assisted catalytic synthesis of ammonia. Ar dilution is a common technique for tailoring plasma discharge properties and has been shown to enhance NH 3 formation when added to N 2 –H 2 plasma. The kinetic model was developed for a coaxial dielectric barrier discharge quartz wool-packed bed reactor operating at near room temperature using a kHz-frequency plasma source. With 30% Ar mixed in a 1:1 N 2 –H 2 plasma at 760 Torr, we find that NH 3 production is dominated by Eley–Rideal (E-R) surface reactions, which heavily involve surface NH x species derived from N and H radicals in the gas phase, while the influence of excited N 2 molecules is negligible. This is contrary to the commonly proposed mechanism that excited N 2 molecules created by Penning excitation of N 2 by Ar(4s) and Ar(4p) play a significant role in assisting NH 3 formation. Our model shows that the enhanced NH 3 formation upon Ar dilution is unlikely due to the interactions between Ar and H species, as excited Ar atoms have a weak effect on H radical formation through H 2 dissociation compared to electrons. We find that excited Ar atoms contribute to 28% of the N radical production in the gas phase via N 2 dissociation, while the rest are dominated by electron-impact dissociation. Furthermore, Ar species play a negligible role in the product NH 3 dissociation. N 2 conversion sensitivity analyses were carried out for electron number density (ne) and reduced electric field (E/N), and contributions from Ar to gas-phase N radical production were quantified. Finally, the model can provide guidance on potential reasons for observing enhanced NH 3 formation upon Ar dilution in N 2 –H 2 plasma beyond changes in the discharge characteristics.

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

An Innovative Insight into Performance Degradation of NCM111 Cathode Induced by Suspension of Operation

The lifespan of lithium-ion batteries varies enormously from fundamental study to practical applications. This big difference has been typically ascribed to the high degree of uncertainty in unpredictable and complicated operation conditions in real-life applications. Here, in this paper, we report that the pause of the charging–discharging process, which is frequently operated in practice but rarely studied in academics, is an important reason for the performance degradation of the NCM111 cathode. It is found that the pause during cycling could trigger a remarkable drop in capacity, giving rise to ~30% more capacity decay compared with the continuously cycled sample. In situ synchrotron X-ray diffraction analysis reveals that the harmful H1–H2 phase transition, which typically appears in the initial cycle but disappears in subsequent cycles, is reactivated by the pausing process. The anisotropic lattice strains that occur during the H1–H2 transition result in mechanical fractures that terminate with an inert NiO-type rock-salt phase on the surface of particles. The present study indicates that the discontinuous usage of rechargeable batteries is also a key factor for cycle life, which might provide a distinct perspective on the performance decay in practical applications.

25 ENERGY STORAGE↗
DOE OSTI2022

COMPASS-U Global Heat Balance Calculations

COMPASS-U is a medium size, high magnetic field experimental tokamak (R = 0.9 m, B t = 5 T, and I p = 2 MA), built at the Institute of Plasma Physics, Czech Academy of Sciences (IPP-CAS). This global heat balance calculation for COMPASS-U was done at Princeton Plasma Physics Laboratory (PPPL). Based on our previous experience of building global thermal model for National Spherical Torus Experiment-Upgrade (NSTX-U) at PPPL, this 2-D global thermal model geometry represents a typical cross section of COMPASS-U machine. The model includes thermal radiation, conduction, and convection among components and also between the machine and outer environment. Helium gas heating/cooling was modeled with fluid element and surface element. This model was used to calculate component temperatures and heat distribution during heat up, cryogenic cool down, normal operation, and fault operation scenarios. Definition of 14 main thermal scenarios was provided by IPP-CAS. First nine are normal operation scenarios. Scenarios #10–#14 are fault scenarios. Results of thermal scenario #5 will be given and discussed in this article, including peak temperatures, temperature ratcheting, energy distribution, and required cooling power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Hydropower Geotechnical Foundations: Current Practice and Innovation Opportunities for Low-Head Applications

Hydropower is a renewable energy resource that produces electricity from flowing water under pressure. Engineered hydropower structures, such as dams, are used to create a hydraulic head, enabling a turbinegenerator unit to convert pressurized flow into electricity. While hydropower has been a source of renewable energy since antiquity, new development in the United States has slowed in recent decades. Based on recent resource assessments, the largest opportunity to expand hydropower in the United States is from new stream-reach development (i.e., new hydropower development along stream-reaches that do not currently have hydroelectric facilities or other forms of infrastructure, such as dams). Roughly 75% of identified new stream-reach development potential is from low-head sites (less than 30 ft of head), which typically suffer from smaller power densities and higher normalized costs, given economies of scale. Hydropower developers and other stakeholders are thus interested in strategies to reduce initial capital costs while practicing sustainable development to maximize environmental compatibility with minimal disruption to natural aquatic life, sediment, and water flows. Historically, civil works have represented a significant cost driver for new hydropower development, with the foundation system representing a major cost component and source of uncertainty. The foundation system is the collection of engineered structural features (e.g. cutoff trenches, walls, grouting, anchors) constructed at or below the preconstruction ground surface that interfaces between the overlying structures (superstructures) and the bed material below (subsurface). Development of a hydropower foundation system must consider the various characteristics of the surrounding stream environment and subsurface while adhering to the engineering requirements of the superstructure that it supports. The care of water, excavation, and other construction activities are important features of foundation design and construction. The design and construction cost of the foundation system is largely dependent on the site geology and riverbed composition and is influenced by the level of geotechnical assessment required and conducted. Thus, a hydropower facility’s geotechnical foundation is often highly site-specific, with proper site selection and assessment being important to project success. The foundation system is designed to provide structural stability (of the foundation and dam), limit seepage, ensure public safety, and maintain functionality for the project life, during both construction and facility operations. Inadequate foundation or dam design can result in dam failure and the uncontrolled release of significant volumes of water, which could cause a high number of casualties and extensive property damage downstream of the failure. According to the Association of State Dam Safety Officials, approximately 30% of all historical dam failures in the United States are attributed to foundation or abutment defects, and another 20% are attributed to piping or seepage through the embankment, foundation, or abutment. To ameliorate these safety considerations, foundations often require massive amounts of construction material (e.g., grout, concrete, engineered dam fill) and long construction times. Foundation design also requires significant analysis prior to construction because the initial in-stream and abutment subsurface conditions are site-specific, and sufficient data for them often are lacking. Current practice requires on-site assessment, using expensive drilling and invasive and non-invasive investigation methods, to determine the expected cost of foundation material and treatment. Additionally, foundation construction often requires site dewatering (and other care of water activities), which involves constructing temporary diversion structures upstream and often downstream, called cofferdams, and water diversion systems that route water around the construction site. Cofferdams and water diversion systems can drastically increase construction costs and contribute to environmental disruption, including modification of flow patterns and benthic habitats. Given the technical, economic, and environmental challenges associated with hydropower foundations, opportunities exist to improve the current state of practice and to develop new and innovative solutions to Hydropower is a renewable energy resource that produces electricity from flowing water under pressure. Engineered hydropower structures, such as dams, are used to create a hydraulic head, enabling a turbinegenerator unit to convert pressurized flow into electricity. While hydropower has been a source of renewable energy since antiquity, new development in the United States has slowed in recent decades. Based on recent resource assessments, the largest opportunity to expand hydropower in the United States is from new stream-reach development (i.e., new hydropower development along stream-reaches that do not currently have hydroelectric facilities or other forms of infrastructure, such as dams). Roughly 75% of identified new stream-reach development potential is from low-head sites (less than 30 ft of head), which typically suffer from smaller power densities and higher normalized costs, given economies of scale. Hydropower developers and other stakeholders are thus interested in strategies to reduce initial capital costs while practicing sustainable development to maximize environmental compatibility with minimal disruption to natural aquatic life, sediment, and water flows. Historically, civil works have represented a significant cost driver for new hydropower development, with the foundation system representing a major cost component and source of uncertainty. The foundation system is the collection of engineered structural features (e.g. cutoff trenches, walls, grouting, anchors) constructed at or below the preconstruction ground surface that interfaces between the overlying structures (superstructures) and the bed material below (subsurface). Development of a hydropower foundation system must consider the various characteristics of the surrounding stream environment and subsurface while adhering to the engineering requirements of the superstructure that it supports. The care of water, excavation, and other construction activities are important features of foundation design and construction. The design and construction cost of the foundation system is largely dependent on the site geology and riverbed composition and is influenced by the level of geotechnical assessment required and conducted. Thus, a hydropower facility’s geotechnical foundation is often highly site-specific, with proper site selection and assessment being important to project success. The foundation system is designed to provide structural stability (of the foundation and dam), limit seepage, ensure public safety, and maintain functionality for the project life, during both construction and facility operations. Inadequate foundation or dam design can result in dam failure and the uncontrolled release of significant volumes of water, which could cause a high number of casualties and extensive property damage downstream of the failure. According to the Association of State Dam Safety Officials, approximately 30% of all historical dam failures in the United States are attributed to foundation or abutment defects, and another 20% are attributed to piping or seepage through the embankment, foundation, or abutment. To ameliorate these safety considerations, foundations often require massive amounts of construction material (e.g., grout, concrete, engineered dam fill) and long construction times. Foundation design also requires significant analysis prior to construction because the initial in-stream and abutment subsurface conditions are site-specific, and sufficient data for them often are lacking. Current practice requires on-site assessment, using expensive drilling and invasive and non-invasive investigation methods, to determine the expected cost of foundation material and treatment. Additionally, foundation construction often requires site dewatering (and other care of water activities), which involves constructing temporary diversion structures upstream and often downstream, called cofferdams, and water diversion systems that route water around the construction site. Cofferdams and water diversion systems can drastically increase construction costs and contribute to environmental disruption, including modification of flow patterns and benthic habitats. Given the technical, economic, and environmental challenges associated with hydropower foundations, opportunities exist to improve the current state of practice and to develop new and innovative solutions to challenges frequently encountered with traditional approaches. With this understanding, it is critically important to understand and document the current state of practice for hydropower geotechnical foundations, identify key challenges, and define opportunities for innovative solutions. To this end, this report documents the current state of practice across the three main phases of geotechnical foundation development: (1) geotechnical site assessment, (2) design, and (3) construction for hydropower systems. It also describes the major challenges with conventional approaches and identifies opportunities for innovation to reduce hydropower foundations costs, timelines, and risks. Key takeaways from this report include the following: Approximately 80% of available low-head sites are expected to have foundations on soil beds rather than rock beds, suggesting that rockfill and earthfill dams may be the most cost-effective conventional dam type for new projects.; Geotechnical and geologic investigation activities are time-consuming and expensive but are essential to define the parameters and criteria needed for foundation design.; Certain riverbed soil and bedrock types present significant technical challenges or require expensive foundation construction, which can prove financially prohibitive for low-head project development.; Modular hydropower design and prefabricated modular foundations represent a promising but unproven paradigm for new hydropower development. Design and construction approaches using optimized and highly repeatable, reliable components would benefit project cost, time, and risk but require additional research and development.; Temporary construction features for foundations, including cofferdams, water diversion, and water control systems, can prove costly and have inherent construction risk.; For economically viable development, hydropower geotechnical foundations should be limited to 4 to 15% of the project’s total initial capital costs. Many proposed projects have experienced cost overruns attributable to foundation difficulties or surprises during construction. These overruns may have been due to inadequate investigations, lack of adequate engineering effort to tailor the structures to site geology and topography, and/or contractual terms, among other considerations.; Challenges for hydropower foundations and opportunities for innovative technology solutions are identified in the following areas (consistent with the three main phases of foundation development): Geotechnical site assessment, Foundation design and materials, Construction methods and technology. Ultimately, this report aims to provide information about geotechnical foundations for low-head hydropower and to motivate transformative technologies to support hydropower growth.

13 HYDRO ENERGY↗
DOE OSTI2022

Robust palladium hydride catalyst for electrocatalytic formate formation with high CO tolerance

Palladium is unique in the electrocatalytic reduction of CO 2 to formate because of its low or even near-equilibrium onset potential. However, the inevitably produced CO molecules poison and deactivate the catalyst surface, resulting in an insufficient operating lifetime (<30 min) for conventional and optimized Pd catalysts. In this work, we present a hydrogen-rich Palladium hydride catalyst (PdH0.5/C) derived from a one-step solvothermal synthesis. This catalyst showed a 93.1 % faradaic efficiency towards formate at – 0.4 V (vs RHE). The working lifetime reached a record of 4 h, which was ~15 times longer than a commercial Pd catalyst and outperformed all previous Pd-based electrocatalysts for CO 2 -HCOO – conversion. The high CO tolerance was attributed to the selectivity improvement induced by lattice hydrogen and the weak CO adsorption strength on diverse active sites (i.e. kink, step, and terrace). Isotopic analysis revealed a direct participation of the lattice hydrogen in the protonation of CO 2 molecules in formate formation.

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

Observation of molecular hydrogen produced from bridging hydroxyls on anatase TiO2(101)

Anatase TiO2 is used extensively in a wide range of catalytic and photocatalytic processes and is a promising catalyst for hydrogen production. Here, we show that molecular hydrogen was produced from bridging hydroxyls, HOb, on the (101) surface of single crystal anatase (a-TiO2(101)). Electron bombardment at 30 K produced bridging oxygen vacancies in the surface. Deuterated bridging hydroxyls, DOb, were subsequently formed via dissociation of adsorbed D2O and confirmed by infrared reflection-absorption spectroscopy. During temperature programmed desorption (TPD) spectroscopy, D2 desorption was observed at 520 K. Density functional theory calculations show that both H2 and H2O production from HOb are endothermic at 0 K on a-TiO2(101), but H2 (H2O) desorption is entropically driven above 230 K (800 K). The calculated activation barrier for H2 desorption is 1.40 eV, which is similar to the desorption energy obtained from analysis of the D2 TPD spectra. The H2 desorption likely proceeds in two steps: H atom diffusion on the surface and then recombination. The calculations indicate that diffusion of the hydrogen atoms into the subsurface could also be an important process. We thank Drs. Z. Dohnalek and R. Rousseau for stimulating discussions. This work was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, Division of Chemical Sciences, Geosciences, & Biosciences and performed in EMSL, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at Pacific Northwest National Laboratory (PNNL). PNNL is a multiprogram national laboratory operated for DOE by Battelle. We also acknowledge the WPI computational facilities for providing resources to perform the DFT studies.

Deskins, N A.↗
DOE OSTI2023

A global perspective on bacterial diversity in the terrestrial deep subsurface

While recent efforts to catalogue Earth’s microbial diversity have focused upon surface and marine habitats, 12–20% of Earth’s biomass is suggested to exist in the terrestrial deep subsurface, compared to ~1.8% in the deep subseafloor. Metagenomic studies of the terrestrial deep subsurface have yielded a trove of divergent and functionally important microbiomes from a range of localities. However, a wider perspective of microbial diversity and its relationship to environmental conditions within the terrestrial deep subsurface is still required. Our meta-analysis reveals that terrestrial deep subsurface microbiota are dominated by Betaproteobacteria, Gammaproteobacteria and Firmicutes, probably as a function of the diverse metabolic strategies of these taxa. Evidence was also found for a common small consortium of prevalent Betaproteobacteria and Gammaproteobacteria operational taxonomic units across the localities. This implies a core terrestrial deep subsurface community, irrespective of aquifer lithology, depth and other variables, that may play an important role in colonizing and sustaining microbial habitats in the deep terrestrial subsurface. An in silico contamination-aware approach to analysing this dataset underscores the importance of downstream methods for assuring that robust conclusions can be reached from deep subsurface-derived sequencing data. Understanding the global panorama of microbial diversity and ecological dynamics in the deep terrestrial subsurface provides a first step towards understanding the role of microbes in global subsurface element and nutrient cycling.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Local-Distance Seismic Event Relocation and Relative Magnitude Estimation, Applications to Mining Related Seismicity in the Powder River Basin, Wyoming

Recent efforts to characterize small ( M w < 3 ) seismic events at local distances have become more important because of the increased observation of human-triggered and induced seismicity and the need to advance nuclear explosion monitoring capabilities. The signals generated by low-magnitude seismic sources necessitate the use of nearby short-period observations, which are sensitive to local geological heterogeneity. Local to near-regional distance ( < 300 km ) surface and shear waves can dominate short-period observations from small, shallow seismic sources. In this work, we utilize these observations to estimate precise, relative locations and magnitudes of ~ 700 industrial mining events in Wyoming, using nearly 360,000 observations. The precise, relative location estimates (with formal location uncertainty estimates of less than 1 km) collapse a diffuse collection of mining events into discrete clusters associated with individual blasting operations. We also invert the cross-correlation amplitudes to estimate precise, relative moment magnitude estimates, which help validate and identify disparities in the event sizes reported by regional network catalogs. Joint use of multiple phases allows for the inclusion of more seismic events due to the increase in the number of observations. In some cases, using a single phase allowed us to relocate only 50% of the original reported seismic events within a cluster. Combining shear- and surface-wave phases increased the number of events to above 90% of the original events, allowing us to characterize a broader range of event sizes, source to station distances, and event distributions. This analysis takes a step toward making a fuller characterization of small industrial seismic events observed at local distances.

58 GEOSCIENCES↗
DOE OSTI2025

Simulation insights into wetting properties of hydrogen-brine-clay for hydrogen geo-storage

Hydrogen geo-storage is attracting substantial interdisciplinary interest as a cost-effective and sustainable option for medium- and long-term storage. Hydrogen can be stored underground in diverse formations, including aquifers, salt caverns, and depleted oil and gas reservoirs. The wetting dynamics of the hydrogen-brine-rock system are critical for assessing both structural and residual storage capacities, and ensuring containment safety. Through molecular dynamics simulations, we explore how varying concentrations of cushion gases (CO 2 or CH 4 ) influence the wetting properties of hydrogen-brine-clay systems under geological conditions (15 MPa and 333 K). We employed models of talc and the hydroxylated basal face of kaolinite (kaoOH) as clay substrates. Our findings reveal that the effect of cushion gases on hydrogen-brine-clay wettability is strongly dependent on the clay-brine interactions. Notably, CO 2 and CH 4 reduce the water wettability of talc in hydrogen-brine-talc systems, while exerting no influence on the wettability of hydrogen-brine-kaoOH systems. Detailed analysis of free energy of cavity formation near clay surfaces, clay-brine interfacial tensions, and the Willard-Chandler surface for gas-brine interfaces elucidate the molecular mechanisms underlying wettability changes. Our simulations identify empirical correlations between wetting properties and the average free energy required to perturb a flat interface when clay-brine interactions are less dominant. Here, our thorough thermodynamic analysis of rock-fluid and fluid-fluid interactions, aligning with key experimental observations, underscores the utility of simulated interfacial properties in refining contact angle measurements and predicting experimentally relevant properties. These insights significantly enhance the assessment of gas geo-storage potential. Prospectively, the approaches and findings obtained from this study could form a basis for more advanced multiscale simulations that consider a range of geological and operational variables, potentially guiding the development and improvement of geo-storage systems in general, with a particular focus on hydrogen storage.

25 ENERGY STORAGE↗
DOE OSTI2022

A three-dimensional pore-scale model for redox flow battery electrode design analysis

A three-dimensional (3-D) pore-scale model has been developed to construct multiscale fibrous electrodes for redox flow batteries (RFB). New designs, such as biporous electrodes modify electrode structures by creating secondary pores on single carbon fiber to reduce internal battery resistance. Existing pore-scale models only resolve electrodes to the single carbon fiber scale and cannot incorporate recent multiscale electrode designs into numerical models. Our new model aims to bridge this gap and provide a tool to rapidly screen new electrode configurations. Two multiscale electrodes, laser-perforated and biporous electrodes, were investigated at varying operation conditions with the proposed 3-D pore-scale models. The laser-perforated electrode exhibits a reduced pressure drop, but follows the same permeability correlation compared to the corresponding pristine electrode. For the biporous electrode, the added specific surface area and faster reaction kinetics from the secondary pores are the most influential factors leading to improved battery efficiency. However, operating the biporous electrode in limiting current density conditions should be avoided due to the decreased mass transfer efficiency and a more significant voltage loss. Finally, we believe that our 3-D pore-scale model can accelerate the flow battery electrode design process and provide new insights into electrode geometry optimizations.

25 ENERGY STORAGE↗
DOE OSTI2025

Studying the Durability of Molybdenum Carbide Aerogel as an Electrocatalyst Support for Proton Exchange Membrane Fuel Cells

Carbon-supported platinum nanoparticles (Pt/C) are currently the state-of-the-art catalyst in proton exchange membrane fuel cells (PEMFCs). Unfortunately, the carbon support lacks the ability to stabilize the metal catalyst, as platinum tends to dissolve and agglomerate, significantly compromising the durability. Herein, we synthesized a ceramic material, molybdenum carbide aerogel (MCAG), and utilized it as a Pt support for PEMFCs, as an alternative for conventional carbon supports. N 2 adsorption and XRD analysis showed that the MCAG possesses a combination of high porosity and a well-defined ceramic crystalline structure. The Pt/MCAG system was studied for its electrocatalytic activity toward ORR in a half-cell and demonstrated satisfactory reaction kinetics and electrochemical active surface area, comparable to the commercial Pt/C. For durability examination, an accelerated stress test (AST) in a single cell was conducted with the Pt/MCAG catalyst at the cathode, following the U.S. DOE electrocatalyst AST protocol. Intermediate polarization curves and cyclic voltammograms were recorded over the course of the AST to monitor the aging process. The results indicated that the MCAG support exhibited superior durability in long-term fuel cell operation compared to a carbon support. The Pt/MCAG system demonstrated stable behavior during the progressive stages of the AST, retaining more than 50% of its initial performance by the end of the test, whereas carbon-supported Pt preserved only 30%.

25 ENERGY STORAGE↗
DOE OSTI2026

Memory-Aware External Facelist Calculation: A Data-Parallel Atomic Hash Counting Approach

Unstructured volumetric meshes serve as fundamental data representations in various scientific simulations and analyses. They play a crucial role in representing complex computational domains and are essential for important numerical techniques, such as finite element analysis. Whenever such a mesh is read from a file, streamed in-situ, or generated by algorithms, scientific visualization libraries rely on calculating the external surface of a geometry, named “external facelist”, to produce a polygonal mesh for rendering. Consequently, external facelist calculation has become one of the most widely used algorithms in the scientific visualization domain, necessitating optimal performance. In this paper, we explore relevant work on external facelist calculation algorithms in two common visualization libraries, VTK and Viskores, assess their performance and memory constraints, and introduce a novel memory-aware external facelist calculation algorithm employing an atomic hash counting approach. This algorithm fully leverages Viskores' data-parallel primitive operations, facilitating its execution across diverse many-core architectures. Our algorithm features the lowest memory footprint on the GPU and the second-lowest on the CPU among all evaluated methods, and it also delivers the fastest performance on both CPU and GPU. It has been made available under an open-source license in the VTK and Viskores visualization systems.

Tsalikis, Spiros [Kitware] (ORCID:0000000151137195↗
DOE OSTI2025

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Simulation and Analysis on Reactor Pressure Vessel (RPV) subjected to Pressurized Thermal Shock (PTS) under SBLOCA scenario by using Cardinal to support the fracture mechanics analyses

The structural components that comprise nuclear reactors and their supporting structures are subjected to harsh operating environments that can challenge their integrity, especially after exposure for extended durations or under accident condition. As one of the most significant components of a Reactor, the Reactor Pressure Vessel (RPV) is exposed to an aggressive environment during the operation time (e.g. more than 40 years). Ageing degradation mechanisms (e.g. thermos-fatigue) could grow initial defects up to a critical size, increasing the susceptibility to failure in the RPV. The conventional methods are mostly based on simple crack and structure geometries. Very limited studies consider the real conditions of the RPV subjected to a thermal shock due to a Loss of Coolant Accident (LOCA). During a LOCA event, the most severe conditions take place when the emergency core cooling (ECC) water is injected inside the cold legs filled initially with hotter water and/or steam. The rapid cooling of the down-comer and the internal RPV surface followed probably by re-pressurization of the RPV causes large temperature gradients and variation of pressure which induces thermal-mechanical stresses. In order to develop the model for integrity assessment of a reactor pressure vessel (RPV) subjected to pressurized thermal shock (PTS), a multi-physics simulation, which includes the thermo-hydraulic, thermo-mechanical and fracture mechanics analyses is necessary. The multi-physics simulations are performed using Cardinal, a wrapping of the GPU-oriented spectral element Computational Fluid Dynamics (CFD) code NekRS and other multi-physics sub-modules within the MOOSE framework. Cardinal now fully supports MOOSE stochastic perturbations of NekRS models with varying boundary conditions, initial conditions, material properties, and any other quantity which is defined by a kernel (such as coefficients in a momentum source model). The implementation is designed in a flexible manner so that scalar values are sent from MOOSE into a user scratch space in NekRS, which can then be applied for any purpose within the NekRS case files (both on the host and device). When modeling PTS, several factors can impact the results significantly. In this report, the impacts of the geometry of the model, Reynolds number and buoyancy effect are investigated. Two geometry, i.e., a simplified model and a realistic RPV model, with both laminar and turbulent flow condition are adopted for the PTS simulation with and without buoyancy effect. The purpose of the investigation is to understand the impact of these factors on the prediction of temperature history of RPV. The accurate prediction on the temperature evolution, which will be exported to Grizzly code for further analyses on the progression of aging mechanisms and their effects on the integrity of RPV structures, is very crucial. Based on the understanding of these factors, a more sophisticated model is built to analysis the PTS under SBLOCA scenario. A literature survey is conducted to pick the SBLOCA scenario for the multi-physics simulation. The analysis helps to explain the form and the transformation of the cold plum when the ECC is activated under SBLOCA. This model can be can be applied to study the PTS effect for different RPV configurations. The results can help to assess structural component degradation for advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗