Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multicomponent”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nuclear–electronic orbital approach to quantization of protons in periodic electronic structure calculations

The nuclear–electronic orbital (NEO) method is a well-established approach for treating nuclei quantum mechanically in molecular systems beyond the usual Born–Oppenheimer approximation. In this work, we present a strategy to implement the NEO method for periodic electronic structure calculations, particularly focused on multicomponent density functional theory (DFT). The NEO-DFT method is implemented in an all-electron electronic structure code, FHI-aims, using a combination of analytical and numerical integration techniques as well as a resolution of the identity scheme to enhance computational efficiency. After validating this implementation, proof-of-concept applications are presented to illustrate the effects of quantized protons on the physical properties of extended systems, such as two-dimensional materials and liquid–semiconductor interfaces. Specifically, periodic NEO-DFT calculations are performed for a trans-polyacetylene chain, a hydrogen boride sheet, and a titanium oxide–water interface. The zero-point energy effects of the protons as well as electron–proton correlation are shown to noticeably impact the density of states and band structures for these systems. These developments provide a foundation for the application of multicomponent DFT to a wide range of other extended condensed matter systems.

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

First-principles calculation of lattice distortion, electronic structure, and bonding properties of GeTe-based and PbSe-based high-entropy chalcogenides

The massive amount of wasted heat energy from industry has pushed the development of thermoelectric (TE) materials that directly convert heat into electricity to a new level of concern. Recently, multicomponent alloys such as GeTe-based and PbSe-based high-entropy (HE) chalcogenides have attracted a great deal of attention due to their potential application as TE materials. The nature of the interatomic bonding, lattice distortion (LD), and the electronic structure in this class of materials is not fully understood. Herein, we report a comprehensive computational investigation of nine GeTe-based HE alloys with eight metallic elements (Ag, Pb, Sb, Bi, Cu, Cd, Mn, and Sn) with large supercells of 1080 atoms each; seven PbSe-based HE solid solutions: Pb0.99−ySb0.012SnySe1−2xTexSx (x = 0.1, 0.2, 0.25, 0.3, 0.35, 0.4, 0.45, with y = 0) with supercells of 1000 atoms each; and five Pb0.99−ySb0.012SnySe1−2xTexSx (y = 0.05, 0.1, 0.15, 0.2, 0.25 with x = 0.25) solid solutions with supercells of 1000 atoms each. All these HE models are theoretically investigated for the first time. The electronic structure, interatomic bonding, charge transfer, and lattice distortion (LD) are investigated by first-principles calculations based on density functional theory. Multicomponent HE alloys can cause a significant LD, which affects their mechanical, thermal, and TE properties. The calculations for the GeTe-based HE chalcogenides showed that they are semiconductors with a narrow bandgap, except for m8, which has a semi-metallic characteristic, and this makes them good candidates for TE applications. For most of these models, the Fermi level shifts upward and locates deeply in the conduction bands, resulting in the enhancement of the electrical conductivity (σ). The bonding properties showed that most bonds in m5 are more dispersed, indicating highest LD and lower lattice thermal conductivity. For PbSe-based HE solid solutions, the LD calculations showed that the models Pb0.99Sb0.012Se0.5Te0.25S0.25 and Pb0.89Sb0.012Sn0.1Se0.5Te0.25S0.25 have the higher LD, and thus a lower lattice thermal conductivity. Such investigations are in high demand since it enables us to design new HE chalcogenides for TE applications. We use the novel concept of total bond order density as a single quantum mechanical metric to characterize the internal cohesion of these HE alloys and correlate with calculated properties, especially the mechanical properties. This work provides a solid database for HE chalcogenides and a road map for many potential applications. Moreover, the computational procedure we developed can be used to design new HE chalcogenides for specific TE applications.

Physics↗
DOE OSTI2022

Low-Pressure Electrochemical Synthesis of Complex High-Pressure Superconducting Superhydrides

There is great current interest in multicomponent superhydrides due to their unique quantum properties under pressure. A remarkable example is the ternary superhydride Li 2 MgH 16 computationally identified to have an unprecedented high superconducting critical temperature T c of ~470 K at 250 GPa. However, the very high synthesis pressures required remains a significant hurdle for detailed study and potential applications. In this Letter, we evaluate the feasibility of synthesizing ternary Li-Mg superhydrides by the recently proposed pressure-potential (P 2 ) method that uniquely combines electrochemistry and applied pressure to control synthesis and stability. Furthermore, the results indicate that it is possible to synthesize Li-Mg superhydrides at modest pressures by applying suitable electrode potentials. Using pressure alone, no Li-Mg ternary hydrides are predicted to be thermodynamically stable, but in the presence of electrode potentials, both Li 2 MgH 16 and Li 4 MgH 24 can be stabilized at modest pressures. Three polymorphs are predicted as ground states of Li 2 MgH 16 below 300 GPa, with transitions at 33 and 160 GPa. The highest pressure phase is superconducting, while the two at lower pressures are not. Our findings point out the potentially important role of the P 2 method in controlling phase stability of complex multicomponent superhydrides.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Hydration and reaction mechanisms on sodium silicate glass surfaces from molecular dynamics simulations with reactive force fields

Recent development of reactive force fields have enabled molecular dynamics simulations of interactions between silicate glasses and water at the atomistic scale. While multicomponent silicate glasses encompass a wide variety of compositions and properties, one common structural feature in these glasses is the combination of the network structure that is made up of silica tetrahedra linked through corner sharing interspersed with network modifiers like alkali and alkaline-earth ions that break up the Si–O–Si linkages by forming nonbridging oxygen. In reactions with water, ion exchange between alkali ions in the glass and proton or hydronium in the solution, as well as hydrolysis reaction of the Si–O–Si linkages and subsequent silanol formation, is observed and well documented. Furthermore, we have used a set of recently developed reactive force field to investigate the reactions between water and the surfaces of silica and sodium silicate glasses of different compositions for reactions up to 8 nanoseconds. Our results indicate sodium leaching into water and diffusion of water molecules up to 25 Å into the glass surface. We examined the structural and compositional changes inside the glass and around the diffused ions and use these to explain the rates of silanol formation at the surface. We also observed proton transport in the glass which has an indirect influence on the silanol formation rates. While the surface of the glass was rough to start with, it undergoes further modification into a hydrated gel-like structure in the glass for up to 5 Å in the higher alkali containing glasses. It was found that the leached sodium ions remain close to the interface and that fragments of silicate network from the surface is capable of dislodging from the bulk glass and enter the aqueous solution. These simulations thus provide insights into the formation and structure of an alteration layers commonly observed in multicomponent silicate glasses corroded in aqueous solutions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Cocrystals combining order and correlated disorder via colloidal crystal engineering with DNA

Colloidal cocrystallization enables the formation of multicomponent materials with unique physicochemical properties, yet the role of nanoparticle (NP) shape and specific ligand interactions to cocrystallize anisotropic and isotropic NPs, with order and correlated disorder, remains underexplored. Here, geometry-inspired strategies along with programmable DNA interactions are combined to achieve structural control of colloidal cocrystal assemblies. Coassembling polyhedral and spherical NPs with complementary DNA yields two classes of cocrystals: one where both components order, and another where polyhedral NPs form a periodic lattice, while spherical NPs remain disordered but spatially correlated with polyhedral edges and corners. The size ratio of the building blocks can be used to control the ordering of spherical NPs—smaller octahedral-to-sphere size ratios favor fully ordered cocrystals. Molecular dynamics simulations further elucidate the role of NP shapes and dimensions in the structural outcome of the cocrystal. This work provides a framework for deliberately targeting and accessing crystals with exotic multicomponent structures.

Li, Yuanwei [Northwestern Univ., Evanston, IL (Uni↗
DOE OSTI2026

The scientific case for concurrent neutron and X-ray scattering and spectroscopy

The interrogation of materials with X-rays or neutrons to determine structure, energetics, and dynamics is fundamental to advancing physical and chemical materials science and enabling innovative material technologies. A persistent challenge in materials development is that progress depends on understanding structure and dynamics across multiple length and time scales in increasingly complex, multicomponent systems featuring interfaces, heterogeneity, and hierarchical organization. Despite rapidly growing demands on materials characterization, current experimental approaches are almost exclusively based on isolated X-ray or neutron scattering and spectroscopy, reflecting a paradigm largely unchanged for decades. To assess the scientific need for a new experimental paradigm, a 3-day workshop sponsored by the U.S. National Science Foundation (NSF) was held at the SpringHill Suites, San Jose, California, from June 2 to 4, 2022. The workshop brought together 70 national and international experts who critically evaluated opportunities enabled by concurrent neutron and X-ray (NeX) scattering, spectroscopy, and imaging experiments. The participants reached a clear consensus that establishing NeX capabilities is crucial for advancing the science of complex materials in the United States. This report illustrates the scientific drivers for NeX experiments through representative examples spanning biomaterials, energy materials, soft matter, nanomaterials, quantum materials, geoscience, and applied materials research. The complementarity of neutrons and X-rays is essential for robust model development and refinement, particularly in multiphase and multicomponent systems. While joint refinement of data from separate experiments is valuable, concurrent measurements uniquely eliminate uncertainties arising from sample evolution, environmental drift, and irreproducibility associated with experiments performed at different locations and times. Realizing NeX capabilities will require the development of new instrumentation, data analysis frameworks, and robust sample environments compatible with both neutron and X-ray probes. Addressing these challenges will enable unambiguous interpretation of complex materials behavior and open new frontiers in materials research.

X-ray↗
DOE OSTI2022

The CALPHAD approach for HEAs: Challenges and opportunities

Phases are key microstructural features that determine material properties. The most intuitive way of representing phase stability in a material is by its phase diagram. The vastness of the composition space is the most compelling reason for investigating high-entropy alloys (HEAs), but it also creates major challenges in controlling phases and microstructures of HEAs due to lack of experimental phase diagrams in highly concentrated multicomponent space. The CALPHAD (Calculation of Phase Diagram) approach, although it needs improvements, is the only viable way to calculate multicomponent phase diagrams. The two major property databases in CALPHAD, Gibbs energy and mobility, when coupled, can provide not only the equilibrium phase information, but also information about diffusion and kinetics. In this article, we will review the use of the CALPHAD approach in understanding phase stability and diffusion in high-entropy alloys. Finally, we will discuss the notable trend of coupling the high-throughput CALPHAD approach with machine learning in HEA design.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Equilipy: a python package for calculating phase equilibria

The CALPHAD (CALculation of PHAse Diagram) approach (Nigel Saunders & Miodownik, 1998) provides predictions for thermodynamically stable phases in multicomponent-multiphase materials across a wide range of temperatures. Consequently, the CALPHAD calculations became an essential tool in materials and process design (Luo, 2015). Such design tasks frequently require navigating a high-dimensional space due to multiple components involved in the system. This increasing complexity demands high-throughput CALPHAD calculations, especially in the rapidly evolving field of alloy design. In response to the need, we developed Equilipy an open-source Python package designed for calculating phase equilibria of multicomponent-multiphase systems. Equilipy is specifically tailored for high-throughput CALPHAD calculations, offering parallel computations across multiple processors and nodes with the given NPT input conditions namely elemental compositions (N), pressure (P), and temperature (T). Equilipy utilizes the program structure and Gibbs energy functions from the Fortran-based program, Thermochimica (Piro et al., 2013), with incorporating a new Gibbs energy minimization algorithm. This algorithm, originally developed by Capitani and Brown in 1987 (Capitani & Brown, 1987), has been revised and implemented to enhance the stability and performance of calculations. The Fortran codes are precompiled and interfaced with Python via F2PY, ensuring high computation speed. Benchmark tests shown in Figure 1 demonstrate that Equilipy’s computation speed is comparable to those of established commercial software, TC-Python and PanPython. This result highlights its efficiency and potential applications in various scientific and industrial fields.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

High Performance non-PGM Transition Metal Oxide ORR Catalysts of PEMFCs

This project was designed to develop acid-stable PGM-free transition metal oxide oxygen reduction reaction (ORR) electrocatalysts to meet or exceed the performance and durability of the DOE 2020 Technical Targets for platinum-group metal (PGM) free electrocatalysts from first-principles to incorporation into membrane-electrode assemblies (MEAs) for polymer electrolyte membrane (PEM) fuel cells. The planned project was to accomplish this goal with a multi-step approach: 1) materials modeling and experimental screening to identify acid-stable oxides with high ORR activity, 2) optimization of catalyst particle size and catalyst/carbon/ionomer catalyst layer composition, 3) fabrication of MEAs for performance and durability testing. During the course of this project: 1) acid-stability descriptors were developed for manganese oxides; 2) a family of acid-stable multicomponent oxides based on antimony were developed; 3) a flexible synthesis for the formation of nanocrystalline nonstoichiometric oxides was developed; 4) limited ORR activity but modest and improving oxygen evolution reaction (OER) activity was measured for multicomponent antimony oxides. The project was programmed into five technical tasks: The development of acid-stable ORR electrocatalytic oxides through 1) identification and optimization of acid-stable oxide compositions, and 2) electrochemical characterization; 3) optimization of catalyst layer composition for MEAs using identified ORR electrocatalysts; 4) MEA fabrication and performance testing to result in performance of 44 mA-cm-2 at 0.9 V vs. RHE; and 5) accelerated-stress testing of optimized MEAs. Because no acid-stable oxide was identified with the requisite activity for ORR (4.4 µA-cm-2oxide intrinsic activity at 0.9 V vs. RHE) within the time and budget allotted in Tasks 1 and 2, the project was halted at the end of phase 1, with no activity in Tasks 3-5. The research output, while not succeeding in developing ORR electrocatalysts, advanced the development of acid-stable oxide materials, showing potential for further improvement as OER electrocatalysts.

08 HYDROGEN↗
DOE OSTI2022

Machine learning pipeline to predict defect behavior in metallic alloy systems

The interaction between defect and solute atoms is critical to the thermodynamic and kinetic behavior of metallic alloys under exposure to high-energy radiation, causing irradiation damage in materials. Radiation can generate non-equilibrium concentrations of point defects such as vacancies and interstitials. The excess point defects not only accelerate diffusional processes such as precipitation that cause radiation embrittlement, but also change the pathway of phase transformations, including nucleation processes. Understanding these defect behaviors is complicated by the challenge and complexity of addressing each possible local and discrete distribution of environments and chemical interactions around targeted defects-solute or solute-solute complexes. To resolve the challenge, machine learning regression techniques have emerged as powerful tools that can train and construct an energy model to accurately describe the chemical interactions of solutes and defects. In Fiscal Year 2022, the work focused on the workflow development and demonstration using machine learning regression, density functional theory, cluster expansion, and Monte Carlo simulation to predict the effects of ternary solute elements (e.g., aluminum and molybdenum) and point defects on the Cr-rich $\alpha^{\prime}$ precipitation in multicomponent FeCr model alloys. The computational outcomes include the prediction of the ternary phase diagram, vacancy formation energy for different compositions, and the effect of vacancies on the nucleation of Cr-rich clusters. The simulations predict a pronounced change of Cr solubility in bcc Fe by the addition of Al and the rejection of Al atoms from $\alpha^{\prime}$ precipitates. Additionally, the simulations show the formation of Cr-vacancy clusters as the initial nuclei for stable nucleation and growth of $\alpha^{\prime}$ particles. The results demonstrate important outcomes and applications of using machine learning pipeline to study model or commercial alloys with multicomponent solute species and point defects.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

FY24 Task 4: Studies of Phosphate and Fluoride Solubility for Dissolution of High Phosphate Tank Waste

Knowledge gaps have been identified in phosphate solubility in almost every single- and multi-component system, and particularly for aluminate in phosphate/hydroxide, where uncertainties in predicted values of aluminum solubility are greater than 50%. Solubility data of multicomponent, aqueous electrolytes containing sodium hydroxide, sodium fluoride, sodium phosphate, sodium nitrite, sodium nitrate, and dissolved gibbsite are also sparse. For example, for solutions of (i) sodium nitrate and sodium phosphate, data are only available at 30 and 50 °C, and for (ii) dissolved gibbsite in sodium hydroxide and sodium phosphate, only two data points are available at 20 and 40 °C. There is no data available on mixtures of sodium hydroxide, sodium phosphate, sodium fluoride, and dissolved gibbsite. These knowledge gaps were identified in a technical review of waste solubility data and the impact of dilution on solution stabilities and will be addressed in this work to predict aluminum hydroxide and sodium phosphate solubility in multicomponent electrolytes upon dilution, and upon variation of temperature. Results will provide the technical basis to develop accurate models for (i) gibbsite solubility and mass transfer of aluminum between solid and liquid forms following the sluicing of sludge and saltcake with water; and (ii) further blending of these suspensions with bismuth from bismuth phosphate waste, and zirconium and uranium left from the fuel decladding. This work will be essential to developing a disposition path for retrieval solutions from bismuth phosphate wastes. This effort would also support sludge washing to further reduce phosphate concentration if that process were to be added back to the flowsheet in the future.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2021

Effect of Polymer Removal on the Morphology and Phase of the Nanoparticles in All-Inorganic Heterostructures Synthesized via Two-Step Polymer Infiltration

Polymer templates play an essential role in the robust infiltration-based synthesis of functional multicomponent heterostructures with controlled structure, porosity, and composition. Such heterostructures are be used as hybrid organic–inorganic composites or as all-inorganic systems once the polymer templates are removed. Using iron oxide/alumina heterostructures formed by two-step infiltration of polystyrene-block-polyvinyl pyridine block copolymer with iron and aluminum precursors from the solution and vapor-phases, respectively, we show that the phase and morphology of iron oxide nanoparticles dramatically depend on the approach used to remove the polymer. We demonstrate that thermal and plasma oxidative treatments result in iron oxide nanoparticles with either solid or hollow morphologies, respectively, that lead to different magnetic properties of the resulting materials. Our study extends the boundaries of structure manipulations in multicomponent heterostructures synthesized using polymer infiltration synthesis, and hence their properties.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Experiments and numerical modeling of surrogate droplet combustion for a 93 octane gasoline calibration fuel

Gasoline is the world’s most prolific liquid transportation fuel for ground transportation vehicles. Yet, a fundamental understanding of its burning characteristics is still lacking. Complicating this matter is the highly multicomponent nature of gasoline (comprised of hundreds of chemical species) and the seasonable variability of its composition due to variations in the crude oil feedstocks used in refining. The need to bring order to this complex problem has led to development of calibration fuels (Tier II in presented study) which on average will provide repeatable metrics that replicate performance of refined gasoline for various engine designs. This study presents results of experiments using the configuration of an isolated droplet to study the combustion dynamics of calibration fuel and three different surrogates (S4, S5, S6) for the calibration fuel. Droplet burning is under microgravity conditions, where gas transport is determined entirely by fuel evaporation in order to promote spherical symmetry in the burning process. The initial droplet diameters ranged between 0.55 mm and 0.60 mm. The droplets were spark-ignited and their burning histories were recorded by digital video cameras to obtain quantitative data. For each fuel, the evolutions of droplet, flame, and soot-shell diameters were measured to calculate the burn rates, flame-droplet diameter ratios (FSR) and soot-droplet diameter ratios. Additionally, a detailed numerical modeling is presented. This simulation includes multicomponent phase equilibrium, unsteady gas and liquid transport, radiative heat transfer, internal liquid diffusion and detailed combustion chemistry of different fuels during their burning process. Experiment results show that the calibration fuel and surrogates have a close burn rates and FSR, and S5 presents a relatively similar results with calibration fuel. Predictions from the numerical simulation using the surrogate are compared with experiment measurements, which showed a good agreement with the experiment results, including the droplet burning rate and FSR.

Guo, Songtao↗
DOE OSTI2023

Thermo-Hydrological Modeling of Thermal Energy Storage in a Depleted Oil Reservoir

Thermal energy storage in oil and gas reservoirs leverages the existing surface and subsurface infrastructure, which can pave the way for economic production of geothermal energy. Existing studies on geothermal energy storage are focused mostly on the use of aquifers with more homogeneous rock and fluid properties. Coupling of heat and fluid flow in a multiphase-multicomponent system, such as an oil reservoir, is imperative especially if existing oil field assets need to be repurposed as required for a sustainable energy transition. The objective is to model the subsurface thermo-hydrological processes associated with reservoir performance and operational sustainability. The model evaluates formation pressure and temperature within the reservoir and at the injection/production wells during multiple charge and discharge cycles. Hot water (approximately 200 degrees C) heated by Concentrating Solar Power (CSP) at high pressure is injected into the existing oil reservoir for storage and produced as thermal energy for power generation, which will be accompanied by enhanced oil recovery. To demonstrate the coupled fluid and heat flow during the injection/production cycle in the subsurface reservoir, TOUGH3 (developed by Berkeley Lab) is used to simulate the thermo-hydrological (TH) processes in a multiphase, multicomponent system. Two well geometries are considered within the reservoir grid: 1) a single-well huff-n-puff system (same well is used for injection and production), and 2) an isolated injection-production well doublet. Seasonal charge and discharge cycling are implemented based on the scheduling specified in the model input file. The model reports pressure, temperature, enthalpy, liquid fluxes, heat fluxes, pore velocities, and changes in porosity & permeability due to temperature and pressure variations during the cyclic Reservoir Thermal Energy Storage (RTES) operations. The results from the simulations can be used to optimize the operational parameters (such as well spacing and injection/production rates) and round-trip efficiency for surface power-plants coupled with thermal energy storage over time. They can also serve as important inputs for levelized cost of storage estimations. The research will help to design and integrate surface renewable energy sources, such as concentrating solar power (CSP), with RTES to help balance out power supply and demand on the grid.

CSP↗
DOE OSTI2023

Integrated membrane material design and system synthesis

In designing membrane systems, the synergy between membrane materials and the process design is often overlooked. In this paper, we present a mixed-integer nonlinear programming (MINLP) model for synthesizing membrane systems while simultaneously designing the respective membrane materials for multicomponent gas separation. The approach considers superstructure representations for systems with: (1) same, (2) potentially different, and (3) property-targeting membrane materials. In the first two systems, the selection of membrane material is a decision, while in the final type, membrane permeances are subject to optimization. Physics-based surrogate models are used to describe permeation in crossflow and countercurrent flow permeators. We show that, through a case study of biogas upgrading, our approach obtains high quality solutions. Furthermore, we use the proposed approach while considering permeance-based production cost to find the optimal membrane.

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

A unified non-equilibrium phase change model for injection flow modeling

The homogenous relaxation model (HRM) is one of the most widely used models to describe the liquid- gas phase transition. However, in its original formulation, it is unable to handle multispecies vapor-liquid equilibrium (VLE), which limits its applicability to single-component fluids. In this work, a unified non-equilibrium phase change model that considers the VLE of multicomponent mixtures is proposed building upon the HRM's structure. A time factor is introduced to mimic the effect of different phase change timescales due to different mechanisms, e.g., cavitation, flash-boiling, and evaporation. Here to assess the model's performance, computational fluid dynamics simulations of the internal and near-nozzle injection flow of the Engine Combustion Network's Spray G injector were performed using the nine-component PACE-20 fuel with both the unified model and the original HRM. The predicted fuel density in the near-nozzle region matched well with X-ray tomography measurements. The simulation results indicated that, whereas the HRM failed to capture the vaporization due to convective mixing between the fuel and ambient gas, the unified model performed well in predicting the mixing-driven vaporization and the corresponding evaporative cooling. Further comparisons using the nine-component fuel formula and a single-component fuel surrogate demonstrated the unified model's ability to predict preferential vaporization, which affects the predictions of local mixture composition and rate of vaporization. Finally, it is shown that the unified model is capable of representing multiple phase change mechanisms, and the relaxation time factor plays an important role in determining the degree of phase change due to the different mechanisms.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Fundamental study of gas species transport in the oxygen electrode of solid oxide fuel and electrolysis cells

A fundamental analysis of multicomponent gas transport models was performed in application to the oxygen electrodes of solid oxide cells. It is common practice to neglect the effect of pressure gradients within oxygen electrodes, even though a net molar flux at the electrolyte surface implies that a pressure gradient must exist. The influence of both Darcy velocity and Knudsen flux are considered in the context of ordinary (Fickian) diffusion, the dusty gas model, and the binary friction model. Comparisons between the models and different sets of assumptions are made via parametric studies on operating load, oxygen partial pressure, microstructural properties, and electrode thickness. Results show that the pressure gradient will have a significant impact on the oxygen concentration distribution and therefore the concentration overpotential. In electrolysis mode, pressure increases up to 1 atm are predicted, indicating that pressure at the electrode/electrolyte interface could contribute to electrode delamination. Additionally, it is found that Darcy's law is insufficient for calculating the pressure distribution without accounting for the flux due to Knudsen diffusion. Additionally, it is found that for the range of properties typical of oxygen electrodes, there is a negligibly small difference between the dusty gas model and binary friction model from a practical standpoint.

08 HYDROGEN↗