SEARCH · Engineering Papers

Results for “surface complexation model”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress in modeling and experimental approaches towards the characterization of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in deep clay/shale/argillaceous rock. International collaboration activities such as heater tests and postmortem analysis of samples recovered from these have elucidated key information regarding changes in the engineered barrier system (EBS) material exposed to years of thermal loads. Chemical and structural analyses of sampled bentonite material from such tests has as well as experiments conducted on these are key to the characterization of thermal effects affecting bentonite clay barrier performance and the extent of sacrificial zones in the EBS during the thermal period. Thermal, hydrologic, and chemical data collected from heater tests and laboratory experiments has been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches (e.g., continuum vs. discrete) to tackle issues related to flow and transport at various scales of the host-rock and EBS design concept. Consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the EBS design. This is particularly important for backfilled repository concepts where temperature plays a key role in the EBS behavior and long-term performance. This report describes multiple R&D efforts on disposal in argillaceous geologic media through development and application of coupled THMC process models, experimental studies on clay/metal/cement barrier and host-rock (argillite) material interactions, molecular dynamic (MD) simulations of water transport during (swelling) clay dehydration, first-principles studies of metaschoepite (UO 2 corrosion product) stability, and advances in thermodynamic plus surface complexation database development. Drift-scale URL experiments provides key data for testing hydrological-chemical (HC) model involving strong couplings of fluid mixing and barrier material chemical interactions. The THM modeling focuses on heater test experiments in argillite rock and gas migration in bentonite as part of international collaboration activities at underground research laboratories (URLs). In addition, field testing at an URL involves in situ analysis of fault slip behavior and fault permeability. Pore-scale modeling of gas bubble migration is also being investigated within the gas migration modeling effort. Interaction experiments on bentonite samples from heater test under ambient and elevated temperatures permit the evaluation of ion exchange, phase stability, and mineral transformation changes that could impact clay swelling. Advances in the development, testing, and implementation of a spent nuclear fuel (SNF) degradation model coupled with canister corrosion focus on the effects of hydrogen gas generation and its integration with Geologic Disposal Safety Assessment (GDSA). GDSA integration activities includes evaluation of groundwater chemistries in shale formations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

A coupled equilibrium boundary layer model with stable water isotopes and its application to local water recycling

The contribution of evapotranspiration (ET) to regional precipitation, known as “local water recycling”, is a key process in the water cycle that can affect water management. However, the amount of planetary boundary layer (PBL) moisture that arises from ET is highly uncertain due to complex atmosphere and land surface conditions. In this study, an idealized two-layer equilibrium boundary layer model was coupled with a stable water isotope module including HDO and H 2 18 O to constrain PBL growing season water transport processes. The model was validated using turbulent heat fluxes and isotope ratios of water vapor (δ m ) and precipitation (δ P ) measured at a cropland site and a nearby tall tower in the Upper Midwest, United States. The results show that the PBL equilibrium features of δ m and δ P are well-constrained by thermal and moisture equilibrium in the PBL. For this study region, the summer values of rain evaporation fraction (f) and local water recycling ratio (LRR) are estimated to be 0.09 and 0.29 ± 0.12, respectively. The summer LRR values for the years 2006–2010 were 0.35, 0.36, 0.17, 0.29, and 0.29, respectively. The small value of LRR in 2008 corresponded to a drought condition with the lowest precipitation and second lowest ET among the five years. The summer magnitude of the amount effect is –2.8‰ (mm day –1 ) –1 and –0.8‰ (mm day –1 ) –1 for HDO and H 2 18 O, respectively. The local water recycling is identified as a significant factor influencing the continental isotope effect. Cropland has likely changed the regional LRR by –7.6 to 19.5% under different pre-agriculture land use scenarios. In conclusion, the feedback processes revealed here indicate that local water recycling is expected to be weakened under drought conditions, but it will be enhanced if irrigation is applied more intensely with more frequent drought events as the climate continues to warm.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Interpretable Machine Learning Models for Practical Antimonate Electrocatalyst Performance

Computationally predicting the performance of catalysts under reaction conditions is a challenging task due to the complexity of catalytic surfaces and their evolution in situ, different reaction paths, and the presence of solid-liquid interfaces in the case of electrochemistry. We demonstrate here how relatively simple machine learning models can be found that enable prediction of experimentally observed onset potentials. Inputs to our model are comprised of data from the oxygen reduction reaction on non-precious transition-metal antimony oxide nanoparticulate catalysts with a combination of experimental conditions and computationally affordable bulk atomic and electronic structural descriptors from density functional theory simulations. From human-interpretable genetic programming models, we identify key experimental descriptors and key supplemental bulk electronic and atomic structural descriptors that govern trends in onset potentials for these oxides and deduce how these descriptors should be tuned to increase onset potentials. Here, we finally validate these machine learning predictions by experimentally confirming that scandium as a dopant in nickel antimony oxide leads to a desired onset potential increase. Macroscopic experimental factors are found to be crucially important descriptors to be considered for models of catalytic performance, highlighting the important role machine learning can play here even in the presence of small datasets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Effects of Size and Shape on the Tolerances for Misalignment and Probabilities for Successful Oriented Attachment of Nanoparticles

Oriented attachment (OA) of nanoparticles is an important pathway of crystal growth, but tools for quantitatively modeling OA are lacking. Here we present several simple models that relate the probability of achieving OA to basic geometric parameters such as particle size, shape, and lattice periodicity. A Moiré-domain model is applied to understand twist-misorientations between parallel surfaces, and it predicts that the range of twist angles yielding perfect OA is inversely related to the width of the contact area. This is confirmed and further developed using a surface functional model, which predicts how crystallographic registration forces drive the emergence of complex orientational energy landscapes. The energy landscapes are predicted to possess local minima that can trap particles in imperfect alignments, and these local minima become deeper and more numerous as the contact area increases, making OA more challenging for large particles. Further, a second set of models is presented to understand the sequence of events by which two crystallographic faces become co-planer after collision. We use a ‘central force approximation’ to quantitatively predict the odds of attaining coalignment between various faces when particles collide with random misalignments, and we show that in the absence of biasing forces, the probability of attaining alignment on a given face is roughly proportional to its solid angle as viewed from the center of the particle. The model predicts that OA is most favorable between well-faceted particles and becomes exceedingly unlikely for large spherical particles that express many microfacets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Perspective—Mass Conservation in Models for Electrodeposition/Stripping in Lithium Metal Batteries

25 ENERGY STORAGE↗

DOE OSTI2022

BattPhase—A Convergent, Non-Oscillatory, Efficient Algorithm and Code for Predicting Shape Changes in Lithium Metal Batteries Using Phase-Field Models: Part I. Secondary Current Distribution

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. Thermodynamics-based equations have been used to study phase changes in lithium batteries using phase-field approaches. However, grid convergence studies and the effect of additional parameters needed to simulate these models are not well-documented in the literature. In this paper, using a motivating example of a moving boundary model in one- and two-dimensions, we show how one can formulate phase-field models, implement algorithms for the same and analyze the results. An open-access code with no restrictions is provided as well. This article concludes with some thoughts on the computational efficiency of phase-field models for simulating dendritic growth.

25 ENERGY STORAGE↗

DOE OSTI2020

Next-Generation Ecosystem Experiments–Arctic Flyer

An important challenge for Earth system models (ESMs) is to accurately represent land surface and subsurface processes and their complex interactions. This is true for all regions of the world, but is especially important for high-latitude Arctic ecosystems that are characterized by ice-rich landscapes where topography, hydrology, vegetation, and biogeochemistry are inexplicitly linked. To address this challenge, the Terrestrial Ecosystem Science (TES) program within the Department of Energy’s (DOE) Office of Biological and Environmental Research (BER) is supporting a next-generation ecosystem experiments project in the Arctic (NGEE–Arctic).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Mechanisms of dissolution from gibbsite step edges elucidated by ab initio molecular dynamics with enhanced sampling

Surface reactions underpin our collective understanding of mineral dissolution. This impacts a range of predictive geochemical models (for weathering, the fate and transport of metals) as well as industrial utilization of minerals. Yet atomistic details are rarely known due to the complex mineral/fluid interfacial environment. There is significant need to understand the mechanistic details of dissolution reactions and how they depend on surface morphology and solution conditions. This work utilizes surface pit models in conjunction with changes to solution composition that mimic pH to explore surface detachment of aluminate from gibbsite, which is a primary source of Al in soils and within the industrial processing of aluminum. Ab initio molecular dynamics simulations with enhanced sampling has been used to explore the detailed process of the detachment, the results of which indicate two potential pathways that are differentiated based upon the extent of water hydration. Here, the heights of the energy barriers depend upon the local morphology which influence the number of bridges (quasi-)simultaneously broken (1 or 2) or the Al-O coordination of the neighboring aluminum atoms (5 or 6) at the armchair edge. pH effects are significant, with a nearly 50% reduction in barrier height under alkaline conditions that are relevant to geothermal fluids and Al extraction from minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Dynamic PRA-Based Estimation of PWR Coping Time Using a Surrogate Model for Accident Tolerant Fuel

In this study, we propose an interpolation-based response surface surrogate methodology to manage a large number of scenarios in dynamic probabilistic risk assessment. It adopts the shape Dynamic Time Warping algorithm to cluster the interpolation neighborhood from time series sample data. The interpolation method was adapted from Taylor Kriging to allow a reduced-order model of the Taylor series. In order to demonstrate its applicability to complex issues in risk assessment for nuclear engineering, an example risk response surface to estimate emergency core cooling system (ECCS) criteria for triplex silicon carbide (SiC) accident-tolerant fuel was constructed. The response surface was exploited to estimate the cumulative failure probability of the fuel cladding structure due to the uncertainties in operator actions and safety systems. The functional failures were assessed based on a combination of individual layer failures computed by coupling Risk Analysis Virtual Environment software with a pressurized water reactor 1000-MW(electric) RELAP5 model and the in-house fuel performance assessment module. Results showed that SiC cladding failure probability spiked less than 1 min after a large-break loss-of- coolant accident whenever the current ECCS criteria for Zircaloy-4 (Zr-4) cladding was used. However, it still provides an increased safety margin of three orders of magnitude compared to Zr-4. This positive margin could be utilized to relax active ECCS requirements by allowing deviations of up to 450 s in its actuation time. The proposed surrogate methodology generated a response surface of SiC cladding failure probability reasonably well, with a significant savings of computation time. This methodology is expected to be useful in the analysis of system response with complex uncertainty sources.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Guidelines to Achieving High Selectivity for the Hydrogenation of α,β-Unsaturated Aldehydes with Bimetallic and Dilute Alloy Catalysts: A Review

Selective hydrogenation of unsaturated aldehydes to unsaturated alcohols is a challenging class of reactions, yielding valuable intermediates for the production of pharmaceuticals, perfumes, and flavorings. On monometallic heterogeneous catalysts, the formation of the unsaturated alcohols is thermodynamically disfavored over the saturated aldehydes. Hence, new catalysts are required to achieve the desired selectivity. In this work, the literature of three major research areas in catalysis is integrated as a step toward establishing the guidelines for enhancing the selectivity: reactor studies of complex catalyst materials at operating temperature and pressure, surface science studies of crystalline surfaces in ultrahigh vacuum, and first-principles modeling using density functional theory calculations. Aggregate analysis shows that bimetallic and dilute alloy catalysts significantly enhance the selectivity to the unsaturated alcohols compared to monometallic catalysts. This comprehensive review focuses primarily on the role of different metal surfaces as well as the factors that promote the adsorption of the unsaturated aldehyde via its C=O bond, most notably by electronic modification of the surface and formation of the electrophilic sites. Furthermore, challenges, gaps, and opportunities are identified to advance the rational design of efficient catalysts for this class of reactions, including the need for systematic studies of catalytic processes, theoretical modeling of complex materials, and model studies under ambient pressure and temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Molecular Dynamics Simulations of a Catalytic Multivalent Peptide–Nanoparticle Complex

Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide–nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure–function relationship, for technological and nanomedical applications.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Global fire modelling and control attributions based on the ensemble machine learning and satellite observations

Contemporary fire dynamics is one of the most complex and least understood land surface phenomena. Global fire controls related to climate, vegetation, and anthropogenic activity are usually intertwined, and difficult to disentangle in a quantitative way. Here, we leveraged an ensemble of five machine learning (ML) models and multiple satellite-based observations to conduct global fire modeling for three fire metrics (burned area, fire number, and fire size), and quantified driving mechanisms underlying annual fire changes in a spatially resolved manner for the period 2003–2019. Ensemble learning is a meta-approach that combines multiple ML predictions to improve accuracy, robustness, and generalization performance. We found that the optimized ensemble ML well reproduced annual dynamics of global burned area (R 2 = 0.90, P < 0.001), total fire numbers (R 2 = 0.86, P < 0.001), and averaged fire size (R 2 = 0.70, P < 0.001). Additionally, the ensemble ML captured key spatial patterns of multi-year mean magnitudes, annual variabilities, anomalies, and trends for different fire metrics. Our ML-based fire attributions further highlighted the dominant role of enhanced anthropogenic activity in reducing global burned area (–1.9 Mha/yr, P < 0.01), followed by climate control (–1.3 Mha/yr, P < 0.01) and insignificant positive vegetation control (0.4 Mha/yr, P = 0.60). Spatially, climate dominated a much larger burned area (53.7%) than human (23.4%) or vegetation control (22.9%); however, the counteracting effects from regional wetting and drying trends weakened the net climate impacts on global burned area. The fire number and fire size exhibited similar spatial control patterns with burned area; globally, however, fire number tended to be more affected by climate while fire size more influenced by human activities. Overall, our study confirmed the feasibility and efficiency of ensemble ML in global fire modeling and subsequent control attributions, providing a better understanding of contemporary fire regimes and contributing to robust fire projections in a changing environment.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Modeling Profiles of Micrometeorological Variables in a Tropical Premontane Rainforest Using Multi-Layered CLM (CLM-ML)

This study updates the multi-layered Community Land Model (CLM-ml) for hillslopes and compares predictions from against observations collected in tropical montane rainforest, Costa Rica. Modifications are made in order to capture a wider array of vertical leaf area distributions, predict CO 2 profiles, account for soil respiration, and adjust wind forcings for difficult topographic settings. Test results indicate that the modified multi-layer CLM model can successfully replicate the shape of various micrometeorological profiles (humidity, CO 2 , temperature, and wind speed) under the canopy. In the single-layer models (CLM4.5 and CLM5), excessive day-to-night differences in leaf temperature and leaf wetness were originally noted, but CLM-ml significantly improved these issues, decreasing the amplitudes of diurnal cycles by 67% and 47%. Sub-canopy considerations, such as canopy shapes and turbulent transfer parameters, also played a significant role in model performance. More importantly, unlike single layer models, the results that CLM-ml produces can be compared to variables measured within the canopy to provide far more detailed diagnostic information. Further observations and model developments, aimed at reflecting surface heterogeneity, will be necessary to adequately capture the complexity and the features of the tropical montane rainforest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Inverse design of pore wall chemistry and topology through active learning of surface group interactions

Design of next-generation membranes requires a nanoscopic understanding of the effect of biologically inspired heterogeneous surface chemistries and topologies (roughness) on local water and solute behavior. In particular, the rejection of small, neutral solutes, such as boric acid, poses a heretofore unsolved challenge. In prior work, a computational inverse design technique using an evolutionary optimization successfully uncovered new surface design strategies for optimized transport of water over solutes in smooth, model pores consisting of two surface chemistries. However, extending such an approach to more complex (and realistic) scenarios involving many surface chemistries as well as surface roughness is challenging due to the expanded design space. In this work, we develop a new approach that uses active learning to optimize in a reduced feature space of surface group interactions, finding parameters that lead to their assembly into ordered, optimal patterns. This approach rapidly identifies novel surface functionalizations that maximize the difference in water and boric acid transport through the nanopore. Moreover, we find that the roughness of the nanopore wall, independent of its chemistry, can be leveraged to enhance transport selectivity: oscillations in the pore wall diameter optimally inhibit boric acid transport by creating energetic wells from which the solute must escape to transport down the pore. Furthermore, this proof-of-concept demonstrates the potential for active learning strategies, in concert with molecular simulations, to rapidly navigate complex design spaces of aqueous interfaces and is promising as a tool for engineering water-mediated surface interactions for a broad range of applications.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Radionuclide Interaction with Hydrothermally Altered Repository Materials (M4SF-23LL010302052)

This progress report (Level 4 Milestone Number M4SF-23LL010302052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline Activity Number SF-23LL01030205. The research is focused on actinide and radionuclide sequestration in hydrothermally altered repository materials. In FY23, a manuscript was in preparation for publication summarizing our analysis of radionuclide sorption and coprecipitation into Fe oxide phases and evaluation of radionuclide partitioning values across a range of radionuclides relevant to performance assessment. We demonstrated our approach in detail using Se sorption and coprecipitation with iron oxide minerals. These data were presented in our FY22 annual report and will not be repeated here. We also initiated experiments to identify radionuclide interaction with hydrothermally altered crystalline repository and backfill materials. Recent research performed at Los Alamos National Laboratory (LANL) and Sandia National Laboratory (SNL) has provided key insights regarding the hydrothermal alteration behavior of bentonite backfill in the presence of repository materials (steel, concrete, etc.). We are now examining how mineral alteration affects retardation behavior of plutonium and a suite of other radionuclides. These experiments also allow us to test the predictive ability of our component additivity approach to surface complexation and ion exchange. Our guiding hypothesis is that a robust surface complexation/ion exchange model and associated database, developed using our L-SCIE approach, can effectively predict changes in radionuclide sorption behavior resulting from the hydrothermal alteration of mineralogy in a repository near field. A short update of results to date is presented below. A manuscript was also in preparation describing a self-consistent model and approach to simulating Se(IV) and Se(VI) sorption to 5 iron oxide phases based on our L-SCIE community database. This will be the first implementation of a multi-mineral and multi-oxidation state sorption model using our new L-SCIE database and workflow. A short summary of these results is presented below.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

DOE OSTI2026

M3SF-25LL010302052 - Radionuclide Interaction with Hydrothermally Altered Repository Materials

This progress report (Level 3 Milestone Number M3SF-25LL010302052) summarizes research conducted at Lawrence Livermore National Laboratory (LLNL) within the Crystalline Host Rock Properties & Processes - LLNL Number SF-25LL01030205. Observed changes in radionuclide sorption after bentonite/clay heating have implications for radionuclide diffusive transport through engineered barriers and must be considered when designing waste disposal repositories. Recent research performed at Los Alamos National Laboratory (LANL) has provided key insights regarding the hydrothermal alteration behavior of bentonite backfill in the presence of repository materials (steel, concrete, etc.). We are examining how this mineral alteration affects retardation behavior of a suite of radionuclides of interest to repository performance assessment. Sorption experiments and data analysis for 233 U were initiated in FY24 following earlier experiments performed on 243 Am, 90 Sr, 137 Cs. In FY25, we completed the 233 U study and initiated and completed a study of 237 Np sorption. Below, we summarize the results and potential impacts of hydrothermal alteration on radionuclide retardation and assess the importance of this process to radionuclide migration from a GHRDC. We also use statistical tools (i.e. PCA) to help us determine the major drivers in affecting changes in measured Kd values induced by hydrothermal alteration. These experiments also allow us to test the predictive ability of our component additivity approach to surface complexation and ion exchange. Our guiding hypothesis is that a robust surface complexation/ion exchange model and associated database can effectively predict changes in radionuclide sorption behavior resulting from the hydrothermal alteration of mineralogy in a repository near field. In November 2024 the paper “Selenium interaction with iron minerals: Quantitative comparison of sorption and coprecipitation impacts on mobility” was published in Applied Geochemistry. A short update of results to date is presented below. In FY25, we also actively supported the DITUSC project, which is part of the EURADII initiative, as associate partners. We are executing the Migration2025 conference and support associated the NEA-TDB and Thermochimie workshops that will provide critical international engagements and develop consensus and synergy in thermodynamics as it relates to supporting the US nuclear waste repository program.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2023

Sensitivities of subgrid-scale physics schemes, meteorological forcing, and topographic radiation in atmosphere-through-bedrock integrated process models: a case study in the Upper Colorado River basin

Abstract. Mountain hydrology is controlled by interacting processes extending from the atmosphere through the bedrock. Integrated process models (IPMs), one of the main tools needed to interpret observations and refine conceptual models of the mountainous water cycle, require meteorological forcing that simulates the atmospheric process to predict hydroclimate then subsequently impacts surface–subsurface hydrology. Complex terrain and extreme spatial heterogeneity in mountainous environments drive uncertainty in several key considerations in IPM configurations and require further quantification and sensitivity analyses. Here, we present an IPM using the Weather Research and Forecasting (WRF) model which forces an integrated hydrologic model, ParFlow-CLM, implemented over a domain centered over the East River watershed (ERW), located in the Upper Colorado River basin (UCRB). The ERW is a heavily instrumented 300 km2 region in the headwaters of the UCRB near Crested Butte, CO, with a growing atmosphere-through-bedrock observation network. Through a series of experiments in the water year 2019 (WY19), we use four meteorological forcings derived from commonly used reanalysis datasets, three subgrid-scale physics scheme configurations in WRF, and two terrain shading options within WRF to test the relative importance of these experimental design choices for key hydrometeorological metrics including precipitation and snowpack, as well as evapotranspiration, groundwater storage, and discharge simulated by the ParFlow-CLM. Our hypothesis is that uncertainty from synoptic-scale forcings produces a much larger spread in surface–subsurface hydrologic fields than subgrid-scale physics scheme choice. Results reveal that the WRF subgrid-scale physics configuration leads to larger spatiotemporal variance in simulated hydrometeorological conditions, whereas variance across meteorological forcing with common subgrid-scale physics configurations is more spatiotemporally constrained. Despite reasonably simulating precipitation, a delay in simulated discharge peak is due to a systematic cold bias across WRF simulations, suggesting the need for bias correction. Discharge shows greater variance in response to the WRF simulations across subgrid-scale physics schemes (26 %) rather than meteorological forcing (6 %). The topographic radiation option has minor effects on the watershed-average hydrometeorological processes but adds profound spatial heterogeneity to local energy budgets (±30 W m−2 in shortwave radiation and 1 K air temperature differences in late summer). This is the first presentation of sensitivity analyses that provide support to help guide the scientific community to develop observational constraints on atmosphere-through-bedrock processes and their interactions.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Investigating the Influence of Ni, ZrO 2 , and Y 2 O 3 from SOFC Anodes on Siloxane Deposition

Siloxanes, as a type of impurity in biogas, can poison the Ni-YSZ anode of SOFCs. However, the influence of individual components of the anode, such as Ni, ZrO 2 , and Y 2 O 3 , on the siloxane deposition process has not been investigated extensively. In this study, Ni, ZrO 2 , and Y 2 O 3 pellets were exposed to H 2 + N 2 + H 2 O + D4 (octamethylcyclotetrasiloxane, 2.5 ppmv) and H 2 + N 2 + D4 (2.5 ppmv siloxane) gas mixtures at 750 °C to investigate their affinity and tolerance for siloxane degradation. Surface morphology analysis and electrochemical analysis including electrochemical impedance spectroscopy (EIS), related distribution of relaxation times (DRT) analysis and equivalent circuit modeling with complex nonlinear least square (CNLS) fitting were conducted. Here, a microstructure parameter—tortuosity factor to porosity ratio $\tau /\varepsilon $ calculated by diffusion polarization resistance was utilized for siloxane deposition evaluation. After comparing pellets surface morphology changes before and after experiments and $\tau /\varepsilon $ change following the contamination, Ni is considered as a major factor in siloxane deposition reactions in Ni-YSZ anode.

25 ENERGY STORAGE↗