Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “system 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.

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

At least 1,153 records · Page 64

Wave energy resources assessment for the multi-modal sea state of Hawaii

This paper describes development and validation of a 32-year wave hindcast for Hawaii as part of the wave energy resources assessment for the Exclusive Economic Zones of the US and affiliated territories. The nested model system comprises structured WAVEWATCH III models for the entire globe and the Central North Pacific as well as an unstructured SWAN model around the Hawaiian Islands with wind forcing from the Global Forecast System Reanalysis and its regional downscaling. The development effort follows the standards introduced by the International Electrotechnical Commission (IEC). The hindcast dataset was firstly validated with buoy measurements in terms of the six IEC wave energy resources parameters including the omnidirectional wave power ????, significant wave height ????????0, energy period ????????, spectral width ????0, direction of maximum directionally resolved wave power ?????????? ????????, and directionality coefficient ????. The comparison shows good overall agreement but with greater uncertainties for the spectral parameters ????0, ???????????????????? , and ???? formulated for narrow band spectra. The multimodal sea state in Hawaii always includes a mix of short-period wind seas, long-period swells from the North and South Pacific as well as waves with moderate periods from multiple sources. We partition the hindcast spectra based on the respective period ranges and directions and validate the partitioned spectral energy and parameters with measured values for implementation in wave energy resources assessment. The partitioned wave parameters from the hindcast allow separate seasonal descriptions of wind seas, north and south swells, and moderate-period waves for in-depth understanding of the wave climate and resources in Hawaii.

Li, Ning↗

Interpretable Deep Learning for the Earth System with Fractal Nets

Focal Area 3: Explainable AI Our confidence in the projections made by Earth System Models (ESMs) depends on understanding them to be, in some important respects, faithful representations of the Earth system. Here we present an “explainable Artificial Intelligence (AI)” method that allows us to uncover the dynamical structure of the observed and modeled Earth system, discover hidden links across wide spatiotemporal scales, target model development efforts at poorly-represented dynamics, and optimize observed or modeled data collection to maximize predictive information. Science Challenge: Dynamical system science for the Earth system poses unique challenges given the large degree of internal climate variability. Thus, tools that help us understand how ESMs succeed and fail at representing these dynamics are crucial, particularly in relation to the observed system. Furthermore, the computational and memory constraints on ESM data output motivate in situ analysis of ESM dynamics, including automatic detection of dynamical shifts. Also, of key importance are procedures that leverage ESMs to optimize observational campaigns for improving process representation, reducing structural uncertainty and improving model skill.

54 ENVIRONMENTAL SCIENCES↗

High-fidelity multiphysics load following and accidental transient modeling of microreactors using NEAMS tools: Application of NEAMS codes to perform multiphysics modeling analyses of micro-reactor concepts

The feasibility of modeling microreactors using high-fidelity models with the Nuclear Energy Advanced Modeling and Simulation (NEAMS) tools is investigated in this report. Three overarching questions guided this research: can NEAMS tools readily be applied for high-fidelity multiphysics modeling of different types of transients in microreactor designs; how accurate are the results obtained; and are improvements needed in accuracy or user experience of NEAMS tools, especially considering newly developed capabilities? This work builds upon FY-2022 work, and two microreactor concepts considering heat pipe (HP-MR) and gas-cooled (GC-MR) technologies were further analyzed using high-fidelity multiphysics simulations. The NEAMS tools considered and coupled within the MultiApp environment are Griffin for neutronics, BISON for thermo-mechanics, Sockeye for heat pipe modeling (in HP-MR), SAM for 1D Fluid – 3D solid modeling of coolant channels and system modeling of balance of plant components (in GC-MR), and the SWIFT code for hydrogen redistribution in hydride moderator. The Heat Pipe MicroReactor (HP-MR) concept was further analyzed in FY-2023 to demonstrate the stochastic TRISO failure modeling capability in BISON to check operational limits of the TRISO fuel. A new full-core Gas-Cooled MicroReactor (GC-MR) model was developed based on the initial assembly-model used in Y-2022 and used for steady-state and accidental depressurization transient simulations. Accuracy of the simulations performed was assessed through 1) verification analyses completed on the different physics with code-to-code comparison, and 2) validation of the multiphysics simulations based on modeling of the Kilopower Reactor Using Stirling Technology (KRUSTY) experiment. In FY-2023, the mesh and model of KRUSTY was updated to closely match publicly available data, and the neutronic model was verified and validated against experimental control rod worth measurements. The multiphysics model of KRUSTY was developed and used for steady-state analysis and for modeling reactivity insertion transient. The calculated power increase and stabilization agrees well with experimental data following adjustment in fuel thermal expansion coefficient. As an important component of this project, the ANL team gathered experience with a wide range of NEAMS tools: the MOOSE Mesh System, Griffin, BISON, SWIFT, Sockeye, SAM, Workbench, and the MOOSE MultiApp System, and provided assessment of new capabilities. Noteworthy are the user assessment of the “vapor-only” flow model in Sockeye and development of a multiphysics startup transient in HP-MR unit cell for use as tutorial in Sockeye. The full-core GC-MR model was used for assessment of SAM for balance of plant modeling and for demonstrating the SWIFT code capability for hydrogen redistribution modeling in multiphysics transient analyses. In this process, several bugs/issues were identified and reported to developers. Finally, the assembly GC-MR model developed in FY-2022 coupling Griffin, BISON and SAM through flow blockage and rod ejection transients was published to the National Reactor Innovation Center (NRIC) Virtual Test Bed (VTB). The Heat Pipe MicroReactor (HP-MR) concept high-fidelity multiphysics coupling of Griffin/BISON/Sockeye in load-following and heat pipe failure transients was also published on the VTB. Those submissions are enabling thorough review of these models as well as wide distribution to industry, regulator, and university users. In this analysis, several new research questions were uncovered, and follow-up analyses are recommended to further improve some models, consider additional transients, and continue development of VTB models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Task 12 Sustainability - Methodological Guidelines on Net Energy Analysis of Photovoltaic Electricity (2nd Edition)

Net Energy Analysis (NEA) is a structured, comprehensive method of quantifying the extent to which a given energy source is able to provide a net energy gain (i.e., an energy surplus) to the end user, after accounting for all the energy losses occurring along the chain of processes that are required to exploit it (i.e., for its extraction, processing and transformation into a usable energy carrier, and delivery to the end user), as well as for all the additional energy 'investments' that are required in order to carry out the same chain of processes. However, this general framework leaves the individual practitioner with a range of choices that can affect the results and thus, the conclusions of a NEA study. The current IEA PVPS guidelines were developed to provide guidance on assuring consistency, balance, and quality to enhance the credibility and reliability of the results from photovoltaic (PV) NEAs. The guidelines represent a consensus among the authors - PV NEA experts in North America and Europe - for assumptions made on PV performance, process inputs and outputs, methods of analysis, and reporting of the results. Guidance is given on photovoltaic-specific parameters used as inputs in NEA and on choices and assumptions in inventory data analysis and on implementation of modelling approaches. A consistent approach towards system modelling, the functional unit, the system boundaries and allocation aspects enhance the credibility of PV electricity NEA studies and enables balanced NEA-based comparisons. Specifically, "apples-to-oranges" comparisons of different energy carriers (e.g., fuels vs. electricity) are not methodologically sound and are to be avoided in all cases; also, any comparison across renewable and non-renewable electricity generation technologies must clearly point out the intrinsically short-term nature of the NEA viewpoint, which does not capture the long-term sustainability implications of renewable vs. non-renewable primary energy harvesting and use: non-renewable primary energy resources are depleted and finally exhausted (irrespective of the size of the EROI), while renewable primary energy resources are not. This document provides an in-depth discussion of a common metric of NEA, namely the energy return on investment (EROI), and how this is to be interpreted vis-a-vis the deceptively similar-sounding metrics in the field of Life Cycle Assessment (LCA): cumulative energy demand (CED) and non-renewable cumulative energy demand (nr-CED) per unit output. Specifically, a number of key differences are highlighted between these metrics as applied to electricity production systems, which are listed in Table S-1.

14 SOLAR ENERGY↗

Process Interactions Can Change Process Ranking in a Coupled Complex System Under Process Model and Parametric Uncertainty

For a complex hydrologic system with multiple processes and process interactions, global sensitivity analysis is often used to identify important or influential parameters for model development and improvement. The identification is complicated by process model uncertainty, when a system process can be represented by multiple process models. This study develops a new total-effect process sensitivity index to identify influential processes under model uncertainty. This is done by extending Sobol's total-effect parameter sensitivity index for one system model to total-effect process sensitivity index for multiple system models to account for uncertainty in process models and model parameters. The total-effect process sensitivity index includes not only the first-order process sensitivity index for measuring the importance of individual processes but also higher-order indices that account for process interactions. The total-effect process sensitivity index can identify an influential process that itself and its interactions with other processes influence a model output. Here, the total-effect process sensitivity index is applied to two numerical examples: (a) Sobol's G*-functions with analytical solutions of first-order and total-effect process sensitivity indices, and (b) groundwater flow models with interactions between recharge, geology, and snowmelt processes. The second evaluation shows that, due to second-order and higher-order process interactions, the first-order and total-effect process sensitivity indices give different process ranking. It is thus necessary to estimate both first-order and total-effect process sensitivity indices to appreciate the difference between the first-order impact of a process alone and the overall total-effect impact of the process itself and its interactions with other processes on a model output.

54 ENVIRONMENTAL SCIENCES↗

Modeling Global Carbon Costs of Plant Nitrogen and Phosphorus Acquisition

Abstract Most Earth system models (ESMs) do not explicitly represent the carbon (C) costs of plant nutrient acquisition, which leads to uncertainty in predictions of the current and future constraints to the land C sink. We integrate a plant productivity‐optimizing nitrogen (N) and phosphorus (P) acquisition model (fixation & uptake of nutrients, FUN) into the energy exascale Earth system (E3SM) land model (ELM). Global plant N and P uptake are dynamically simulated by ELM‐FUN based on the C costs of nutrient acquisition from mycorrhizae, direct root uptake, retranslocation from senescing leaves, and biological N fixation. We benchmarked ELM‐FUN with three classes of products: ILAMB, a remotely sensed nutrient limitation product, and CMIP6 models; we found significant improvements in C cycle variables, although the lack of more observed nutrient data prevents a comprehensive level of benchmarking. Overall, we found N and P co‐limitation for 80% of land area, with the remaining 20% being either predominantly N or P limited. Globally, the new model predicts that plants invested 4.1 Pg C yr −1 to acquire 841.8 Tg N yr −1 and 48.1 Tg P yr −1 (1994–2005), leading to significant downregulation of global net primary production (NPP). Global NPP is reduced by 20% with C costs of N and 50% with C costs of NP. Modeled and observed nutrient limitation agreement increases when N and P are considered together ( r 2 from 0.73 to 0.83).

54 ENVIRONMENTAL SCIENCES↗

Towards provision of regularly updated climate data from the Coupled Model Intercomparison Project

The Coupled Model Intercomparison Project (CMIP) is a flagship of the World Climate Research Programme (WCRP). CMIP has become a recognised ‘brand’ in climate circles evolving over the last thirty years from a targeted research activity by a small number of climate modelling centres intercomparing their Earth System Model (ESM) simulations to a broad international coordinated research effort (Durack et al, 2025). CMIP is organized as a research activity leveraging funded and in-kind contributions from experts within modelling centres and the broader scientific community supported more recently by a fully-funded International Project Office. Within CMIP, Model Intercomparison Projects (MIPs) are community-designed to understand past, present and future climate. CMIP data provides a valuable resource for climate research and is routinely used to assess model representation of climate processes and test scientific hypotheses in the context of model uncertainty and (forced and internal) variability as evident from its prolific use in scientific publications1 . The impact relies on enabling infrastructure (most prominently via the Earth System Grid Federation (ESGF)), which allows sharing of simulation output, provision of the boundary conditions used in each simulation, and definition of the data standards that are essential to facilitating wide use of the data. The impact is supplemented by the wide-ranging scrutiny to which model simulations are subjected. Beyond its use in research, CMIP data is a key resource for communities producing derived climate information from downscaling and impact studies, such as the Coordinated Regional Downscaling Experiment (CORDEX; Gutowski et al., 2016) and the Intersectoral Impacts MIP (ISIMIP; Frieler et al., 2024). Government, academic and commercial entities also increasingly rely on CMIP and its downstream data for climate risk assessments and climate services (for example, Copernicus Climate Change Service and World Bank portal). This means that, although CMIP is a research activity, it increasingly serves a secondary and very relevant role as a provider of climate data – a long-recognised dichotomy (Stevens, 2024). Research and applications have distinct needs, with the former requiring flexibility and generality and the latter consistency. Here we explain how the design of the research activity has been adapted to reduce the burdens imposed by applications and how the research infrastructure might evolve to further enable scientific inquiry. We propose one possible approach to consistently providing model information and projections for applications in the future.

Environmental sciences↗

Recent Code Developments/Improvements to SAM Multi-dimensional Flow Model

The System Analysis Module (SAM) is under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. In addition to its conventional one-dimensional flow network module, the multi-dimensional flow model of SAM offers significant benefits to its end-users, including the U.S. NRC, who uses it to develop reference models for advanced reactor concepts. This report provides a summary of the recent progress achieved under DOE-NE’s Nuclear Energy Advanced Modeling and Simulation program in the continuous code development and improvements of the SAM code, specifically its multi-dimensional flow model. The improvements include enhancements to the physical model, code usability, addressing user feedback, and compliance with the SQA program standards that aim to enhance and maintain the quality of the software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Boundary Layer Controls on the Shallow-to-Deep Cumulus Transition (Final Technical Report)

This project advanced our understanding of the processes governing cumulus cloud formation and provides an improved observational basis for validating earth system models. To be specific, the project used laser- and radar-remote sensors to examine the physical properties of updrafts that rise from earth’s surface and initiate clouds deeper in the atmosphere. These updrafts comprise “thermals” and “plumes” and occur at small spatial and temporal scales (e.g., 10s of minutes, 100s of meters). These small scales preclude explicit representation in most earth system and climate models, and thus necessitate “sub-grid-scale” parameterization of updraft processes. The innovation of this project was to directly measure the size, shape, strength and water vapor content of these updrafts with Doppler and Raman lidars, respectively, and to link these updraft properties to cloud processes using vertically pointed weather radars. The resulting data sets comprise 100s of thousands of updrafts and thousands of clouds, which far exceeds previous efforts, and thereby provides a robust statistical and physical representation of these processes. From these large datasets the project produced a sequence of scientific analyses that: (1) Elucidate how variations in the turbulent structure of the convective boundary layer control shallow cumulus convection, (2) Quantify the upward transport of water vapor to cloud base via thermals and plumes, (3) Validate large-eddy simulations of updrafts and shallow convective clouds, (4) Demonstrate a size-to-strength relationship between updraft width and updraft speed, and (5) Demonstrate how updrafts interact with the stability at the top of the convective boundary layer to modulate the depth and vigor of convective clouds. These results have been disseminated via several published journal articles, academic theses, and conference presentations. Collectively these results contribute to the Atmospheric System Research (ASR) program’s goal to “improve understanding of the key cloud, aerosol, precipitation, and radiation processes that affect the Earth’s radiative balance and hydrological cycle, particularly processes that limit the predictive ability of regional and global earth system models”.

54 ENVIRONMENTAL SCIENCES↗

Physics-informed regularization and structure preservation for learning stable reduced models from data with operator inference

Operator inference learns low-dimensional dynamical-system models with polynomial nonlinear terms from trajectories of high-dimensional physical systems (non-intrusive model reduction). Here, this work focuses on the large class of physical systems that can be well described by models with quadratic and cubic nonlinear terms and proposes a regularizer for operator inference that induces a stability bias onto learned models. The proposed regularizer is physics informed in the sense that it penalizes higher-order terms with large norms and so explicitly leverages the polynomial model form that is given by the underlying physics. This means that the proposed approach judiciously learns from data and physical insights combined, rather than from either data or physics alone. Additionally, a formulation of operator inference is proposed that enforces model constraints for preserving structure such as symmetry and definiteness in linear terms. Numerical results demonstrate that models learned with operator inference and the proposed regularizer and structure preservation are accurate and stable even in cases where using no regularization and Tikhonov regularization leads to models that are unstable.

97 MATHEMATICS AND COMPUTING↗

SERGHEI (SERGHEI-SWE) v1.0: a performance-portable high-performance parallel-computing shallow-water solver for hydrology and environmental hydraulics

The Simulation EnviRonment for Geomorphology, Hydrodynamics, and Ecohydrology in Integrated form (SERGHEI) is a multi-dimensional, multi-domain, and multi-physics model framework for environmental and landscape simulation, designed with an outlook towards Earth system modelling. At the core of SERGHEI's innovation is its performance-portable high-performance parallel-computing (HPC) implementation, built from scratch on the Kokkos portability layer, allowing SERGHEI to be deployed, in a performance-portable fashion, in graphics processing unit (GPU)-based heterogeneous systems. In this work, we explore combinations of MPI and Kokkos using OpenMP and CUDA backends. In this contribution, we introduce the SERGHEI model framework and present with detail its first operational module for solving shallow-water equations (SERGHEI-SWE) and its HPC implementation. This module is designed to be applicable to hydrological and environmental problems including flooding and runoff generation, with an outlook towards Earth system modelling. Its applicability is demonstrated by testing several well-known benchmarks and large-scale problems, for which SERGHEI-SWE achieves excellent results for the different types of shallow-water problems. Finally, SERGHEI-SWE scalability and performance portability is demonstrated and evaluated on several TOP500 HPC systems, with very good scaling in the range of over 20 000 CPUs and up to 256 state-of-the art GPUs.

58 GEOSCIENCES↗

Inverter-Based Operation of Maui: Electromagnetic Transient Simulations

As larger and larger power systems approach and reach 100% inverter-based resource (IBR) operation during some hours of the year, questions arise regarding the stability of such extremely-high IBR power systems, the potential need for grid-forming (GFM) inverter technology, and the potential need for synchronous condensers. Relatedly, questions also arise about the ability of conventional positive sequence power system modeling tools to capture high-IBR system dynamics. This presentation introduces electromagnetic transient (EMT) simulations in PSCAD of the near-future (year 2023) Maui power system at and near 100% IBR operation and compares those simulations to positive sequence (PSSE) simulations. The Maui PSCAD model is parallelized on 30 cores and includes the entire transmission system (>200 three-phase buses), >170 individual and aggregate IBR models, four wind plants, and three synchronous generators plus six synchronous condensers at two locations. We investigate system stability with varying levels of inertia using conventional grid-following IBR controls, and then we investigate the impact of GFM controls on stability. Results suggest that: 1) positive sequence simulations can miss key dynamics in extremely high IBR cases; 2) EMT simulations can also miss key dynamics if inverter inner control dynamics are not modeled; 3) synchronous condensers can stabilize a system in which 100% of the energy is supplied by IBRs, even conventional grid-following IBRs; 4) GFM controls on just some of the IBRs can stabilize a 100% IBR power system, even if that system has zero inertia (i.e. no synchronous condensers, though synchronous condensers may be needed for other purposes such as protection system operations).

41 EE - Solar Energy Technologies Office (EE-4S)↗

SNL will be able to provide well-calibrated models to weapons systems analysts by integrating an engine for constitutive model calibration with the GRANTA materials database

Researchers at Sandia National Laboratories have integrated the GRANTA materials database with the MatCal calibration engine to calibrate material models from characterization data. GRANTA is gaining acceptance across the NNSA Tri-lab complex and is being populated with weapons-specific test data by Sandia experimentalists. To use that data to create material models for use by weapons systems analysts, MatCal has been enabled import calibration data and test conditions from GRANTA to quickly and reproducibly produce a calibrated set of parameters for a given constitutive model. The team is currently working to store the parameters characterizing material behavior in GRANTA to make them accessible by all weapons analysts.

36 MATERIALS SCIENCE↗

Network analysis of memristive device circuits: dynamics, stability and correlations

Abstract Networks with memristive devices are a potential basis for the next generation of computing devices. They are also an important model system for basic science, from modeling nanoscale conductivity to providing insight into the information-processing of neurons. The resistance in a memristive device depends on the history of the applied bias and thus displays a type of memory. The interplay of this memory with the dynamic properties of the network can give rise to new behavior, offering many fascinating theoretical challenges. But methods to analyze general memristive circuits are not well described in the literature. In this paper we develop a general circuit analysis for networks that combine memristive devices alongside resistors, capacitors and inductors and under various types of control. We derive equations of motion for the memory parameters of these circuits and describe the conditions for which a network should display properties characteristic of a resonator system. For the case of a purely memresistive network, we derive Lyapunov functions, which can be used to study the stability of the network dynamics. Surprisingly, analysis of the Lyapunov functions show that these circuits do not always have a stable equilibrium in the case of nonlinear resistance and window functions. The Lyapunov function allows us to study circuit invariances, wherein different circuits give rise to similar equations of motion, which manifest through a gauge freedom and node permutations. Finally, we identify the relation between the graph Laplacian and the operators governing the dynamics of memristor networks operators, and we use these tools to study the correlations between distant memristive devices through the effective resistance.

97 MATHEMATICS AND COMPUTING↗

The MOSAiC Distributed Network: Observing the coupled Arctic system with multidisciplinary, coordinated platforms

Central Arctic properties and processes are important to the regional and global coupled climate system. The Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) Distributed Network (DN) of autonomous ice-tethered systems aimed to bridge gaps in our understanding of temporal and spatial scales, in particular with respect to the resolution of Earth system models. By characterizing variability around local measurements made at a Central Observatory, the DN covers both the coupled system interactions involving the ocean-ice-atmosphere interfaces as well as three-dimensional processes in the ocean, sea ice, and atmosphere. The more than 200 autonomous instruments (“buoys”) were of varying complexity and set up at different sites mostly within 50 km of the Central Observatory. During an exemplary midwinter month, the DN observations captured the spatial variability of atmospheric processes on sub-monthly time scales, but less so for monthly means. They show significant variability in snow depth and ice thickness, and provide a temporally and spatially resolved characterization of ice motion and deformation, showing coherency at the DN scale but less at smaller spatial scales. Ocean data show the background gradient across the DN as well as spatially dependent time variability due to local mixed layer sub-mesoscale and mesoscale processes, influenced by a variable ice cover. The second case (May–June 2020) illustrates the utility of the DN during the absence of manually obtained data by providing continuity of physical and biological observations during this key transitional period. We show examples of synergies between the extensive MOSAiC remote sensing observations and numerical modeling, such as estimating the skill of ice drift forecasts and evaluating coupled system modeling. The MOSAiC DN has been proven to enable analysis of local to mesoscale processes in the coupled atmosphere-ice-ocean system and has the potential to improve model parameterizations of important, unresolved processes in the future.

54 ENVIRONMENTAL SCIENCES↗

Active and passive cooling approaches for a Southern California residential community

This study assesses cooling strategies in a low-income community in Southern California that lacks air conditioning and struggles with heat and air pollution. We used an urban building energy model and an electric distribution system model to evaluate active and passive cooling measures. The most effective space cooling measures were high-performance air-source heat pumps, cool coatings, window films, and harnessing the space cooling effect from heat pump water heaters. The results show that combining heat pump water heaters with window films and cool coatings reduces heat index hazard hours within buildings by 95 % to 99 % but increases total energy costs (equipment costs plus changes in utility bills) by 20 % to 60 % where the higher end includes building electrical upgrades. These measures also led to increased space heater use during colder months to avoid overcooling. Replacing conventional heaters with air source heat pumps eliminated unsafe indoor temperatures and reduced total energy use, but increased cost by 125 % to 150 %. In total, using heat pumps for space and water heating could reduce primary energy use by up to 57 %. The higher cost of active and passive cooling measures can be mitigated by existing and emerging incentive programs, especially those that support heat pumps. Electric distribution upgrades to support community electrification are estimated to increase utility costs by $\$$25 to $\$$40 per ratepayer per year. The results underscore the potential and challenges of adapting building infrastructure in communities at risk from climate and environmental stressors.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Precipitation of dopants on acceptor-doped LaMnO3±δ revealed by defect chemistry from first principles

Perovskite oxides degrade at elevated temperatures while precipitating dopant-rich particles on the surface. A knowledge-based improvement of surface stability requires a fundamental and quantitative understanding of the dopant precipitation mechanism on these materials. We propose that dopant precipitation is a consequence of the variation of dopant solubility between calcination and operating conditions in solid oxide fuel cells (SOFCs) and electrolyzer cells (SOECs). To study dopant precipitation, we use 20% (D = Ca, Sr, Ba)-doped LaMnO3+δ (LDM20) as a model system. We employ a defect model taking input from density functional theory calculations. The defect model considers the equilibration of LDM20 with a reservoir consisting of dopant oxide (DO), peroxide (DO2), and O2 in the gas phase. The equilibrated non-stoichiometry of the A-site and B-site as a function of temperature, T, and oxygen partial pressure, p(O2), reveals three regimes for LDM20: A-site deficient (oxidizing conditions), A-site rich (atmospheric conditions), and near-stoichiometric (reducing conditions). Assuming an initial A/B non-stoichiometry, we compute the dopant precipitation boundaries in a p–T phase diagram. Our model predicts precipitation both under reducing (DO) and under highly oxidizing conditions (DO2). We found precipitation under anodic, SOEC conditions to be promoted by large dopant size, while under cathodic, SOFC conditions precipitation is promoted by initial A-site excess. The main driving forces for precipitation are oxygen uptake by the condensed phase under oxidizing conditions and oxygen release assisted by B-site vacancies under reducing conditions. Possible strategies for mitigating dopant precipitation under in electrolytic and fuel cell conditions are discussed.

Hess, Franziska (ORCID:000000024460643X)↗