Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solid modeling”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Thermodynamic description of molten salt systems: KCl-LiCl-NaCl and KCl-LiCl-NdCl 3

The mixture of KCl and LiCl has been used as electrolyte in the electrorefining process to recover uranium from used nuclear fuels due to the low melting point. However, lanthanides and sodium in the reactor waste continuously dissolve into it and thus alter its thermodynamic properties. To understand the thermodynamic behavior evolution of the electrolyte with the accumulation of impurities, thermodynamic modeling for KCl-LiCl-NaCl and KCl-LiCl-NdCl 3 and four constituent binary systems in the entire composition space was performed using the CALPHAD (CALculation of PHAse Diagrams) approach. The ionic liquid was described by the two-sublattice model, where neutral species were introduced to consider short-range ordering (SRO) within the melt, whereas the solid solution was modeled based on the Compound Energy Formalism. Literature data on phase equilibria and thermochemical properties were critically evaluated and used during the optimization of thermodynamic parameters for KCl-LiCl-NaCl and KCl-LiCl-NdCl 3 and their subsystems. The calculated phase diagrams and mixing enthalpies are in good agreement with the experimental data. The thermodynamic modeling for the KCl-LiCl-NdCl 3 system was carried out for the first time. To fill the gap in experimental measurement, enthalpy of mixing for the KCl-LiCl-NdCl 3 melt was estimated using the surrounded-ion model. These data then served as critical inputs for thermodynamic optimization. Furthermore, the present study can provide insights into thermodynamic property evolution of the electrolyte and solubility limit of various impurities during the electrorefining process.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Physics-Based Model to Represent Membrane-Electrode Assemblies of Solid-Oxide Fuel Cells Based on Gadolinium-Doped Ceria

This paper reports a physics-based model that predicts membrane-electrode assembly (MEA) performance of solid-oxide fuel cells (SOFCs) with Ce 0.9 Gd 0.1 O 2− δ (GDC10) electrolyte membranes. The paper derives self-consistent thermodynamic and transport properties for GDC1o mobile charged defects (oxide vacancies and reduced-ceria small polarons) by fitting published measurements of oxygen non-stoichiometry and conductivity over ranges of temperature and O 2 partial pressures. The button-cell model is applied to evaluate how mixed ionic-electronic conductivity influences the performance of an SOFC MEA with a GDC10 electrolyte sandwiched between a porous, composite Ni-GDC10 anode and a porous, composite cathode of Sm 0.5 Sr 0.5 CoO 3− δ (i.e., SSC) and GDC10. SSC properties are also derived by fitting published conductivity and oxygen non-stoichiometry measurements. Mixed conductivity of GDC10 and competing charge transfer reactions at both electrodes reduce open circuit voltages due to leakage current and buildup of defect concentrations at electrode-electrolyte interfaces. To fit polarization data, the button-cell model includes heterogeneous reaction rates for defect incorporation on the GDC10 surface along with Butler–Volmer expressions derived for competing charge transfer reaction rates from rigorous analyses assuming rate-limiting, elementary charge transfer reactions for each electrode. The calibrated MEA model can support rigorous SOFC modeling with GDC10 electrolytes over the range of conditions within a fully operating cell.

Electrochemistry↗
DOE OSTI2020

Precautions of Using Three-Electrode Configuration to Measure Electrode Overpotential in Solid Oxide Electrochemical Cells: Insights from Finite Element Modeling

Accurate determination of electrode overpotentials is essential to assess the performance of the electrode and understand the rate-limiting steps involved. Three-electrode configuration with the use of a reference electrode at a fixed potential is a standard way to measure overpotential of a specific electrode in liquid electrochemical systems. However, application of such three-electrode configuration to solid electrochemical cells for overpotential determination is not straightforward and requires extra caution. Here we report a theoretical Finite Element Analysis on the geometrical requirements for which the reference electrode can be applied to anode- or electrolyte-supported solid oxide button cells. The modeling results suggest that the symmetry of the working and counter electrodes is the key factor determining if a reference electrode is suitable to use. For anode-supported fuel cells with asymmetrical working and counter electrodes, reference electrode cannot be used under all circumstances. To use reference electrode for overpotential measurements, electrolyte-supported cells with symmetrical semicircular-shaped electrodes are preferred. A data processing procedure has also been presented to obtain the electrode overpotential from the measured potential using the three-electrode scheme in solid oxide electrochemical cells.

25 ENERGY STORAGE↗
DOE OSTI2021

A graph-based network for predicting chemical reaction pathways in solid-state materials synthesis

Abstract Accelerated inorganic synthesis remains a significant challenge in the search for novel, functional materials. Many of the principles which enable “synthesis by design” in synthetic organic chemistry do not exist in solid-state chemistry, despite the availability of extensive computed/experimental thermochemistry data. In this work, we present a chemical reaction network model for solid-state synthesis constructed from available thermochemistry data and devise a computationally tractable approach for suggesting likely reaction pathways via the application of pathfinding algorithms and linear combination of lowest-cost paths in the network. We demonstrate initial success of the network in predicting complex reaction pathways comparable to those reported in the literature for YMnO 3 , Y 2 Mn 2 O 7 , Fe 2 SiS 4 , and YBa 2 Cu 3 O 6.5 . The reaction network presents opportunities for enabling reaction pathway prediction, rapid iteration between experimental/theoretical results, and ultimately, control of the synthesis of solid-state materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Eight decades of research on the long-term health effects of radiation in atomic bomb survivors and their offspring

Abstract This year marks the 80th anniversary of the atomic bombings of Hiroshima and Nagasaki. Over the past eight decades, large-scale cohort studies of atomic bomb survivors and their offspring conducted by the Radiation Effects Research Foundation and its predecessor, the Atomic Bomb Casualty Commission, have provided critical insights into the long-term health effects of radiation exposure. Key findings include early identification of radiation-associated leukemia, as well as excess risks of all solid cancers combined, and most individual cancer sites. Observed radiation dose–response relationships have generally followed a linear-quadratic model for leukemia and a linear model for all solid cancers. Recent findings indicating possible upward curvature in the dose–response for all solid cancers may reflect underlying heterogeneity in factors related to individual cancer sites and should be explored further. Generally, younger age at exposure, lower attained age, and female sex appear to show greater radiation sensitivity for all solid cancers combined but results differ by individual cancer site. Recent studies have also identified potential radiation-related excesses for non-cancer diseases such as cataracts, various circulatory diseases, and others. Studies of heritable effects on the offspring of exposed atomic bomb survivors, in contrast, have shown no elevated risk to date in offspring from parental radiation exposure, either at the molecular or disease level. With the cooperation of the atomic bomb survivors and their families, Radiation Effects Research Foundation’s research will continue to play a crucial role in informing the health of survivors, their families, and global radiation protection in the decades to come.

Oncology↗
DOE OSTI2021

Study the thermal radiation effects in gas-solid flows with gray and non-gray P1 models implemented in MFiX

Thermal radiation is a dominant mode of heat transfer in combustion/gasification, packed/circulating bed reactors, and energy storage/conversion devices. In hydrocarbon fuel combustion, the high absorption and emission of product gases like CO 2 and H 2 O significantly affect the heat transfer characteristics. It is essential to understand the radiative energy propagation in such applications. In the present study, we developed a framework to implement the P-1 radiation model in MFiX for Eulerian-Eulerian (Two fluid-TFM) and Eulerian-Lagrangian (Discrete element-DEM and Particle in cell-PIC) models. The P-1 is a simplification of the spherical harmonics method. The framework has been verified on various flue gas conditions. The verified solver results are validated against the benchmark results available in the literature. The verified and validated gray and non-gray weighted sum of gray gas models is evaluated for dry and wet flue gas conditions. Further, they are also used to study gas-solid flows in a fluidized bed/riser.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2022

Thermodynamic modeling of adsorption at the liquid-solid interface

Adsorptive separation techniques are significantly energy efficient in comparison to conventional thermal separation techniques such as distillation. Despite extensive research and development activities undertaken for mixed-gas adsorption, the use of adsorption techniques for the separation of multicomponent liquid mixtures is still limited. A major barrier is the lack of accurate adsorption thermodynamic models, which form the scientific foundation of process simulation of such systems, making the translation to industrial scale challenging. In this work, we have rigorously computed the surface excess of adsorption for six binary liquid mixtures on silica gel at 303 K using the frameworks of the generalized Langmuir isotherm (gL) and the adsorbed solution theory (AST). The six binary liquid mixtures were formed by the pair-wise combinations of four components: benzene, 1,2-dichloroethane, cyclohexane, and n-heptane. We have based our calculations by considering simultaneous equilibria of three phases: saturated vapor phase, bulk liquid phase, and adsorbed phase. Further, the composition of the corresponding saturated vapor phase was estimated by the Nonrandom Two-Liquid activity coefficient model and experimental vapor-liquid equilibria data. The activity coefficients of the adsorbed phase, the central issue of multicomponent adsorption thermodynamics, were calculated using the adsorption Nonrandom Two-Liquid activity coefficient model. Devoid of simplifying assumptions, gL and AST provide rigorous thermodynamic frameworks for adsorption equilibria of multicomponent liquid mixtures.

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

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory. The shear modulus model for the ambient phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress–strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston–Tonks–Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to ⁠10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid–liquid but never a fully liquid state, never fully lost strength, stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a 10 multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory (DFT). The shear modulus model for the ambient β phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the β phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress-strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston-Tonks-Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to 10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability 20 of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid-liquid but never fully liquid state, never fully lost strength, 25 stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Semi-Empirical Material Model for Hydrogen Uptake Kinetics by 1,4-bis(phenylethynyl)benzene (DEB)-based Getters

We report on development of a semi-empirical model describing the kinetics of H 2 /D 2 reaction with 1,4-bis(phenylethynyl)benzene (DEB)-based getters in pellet and powder form. The physical basis of the model is solid-state diffusion of DEB plus reaction of DEB with H atoms on the surface of nanometer-scale Pd catalyst particles. The purpose of this material model is to facilitate more accurate predictions of performance for DEB-based getter assemblies. These predictions will be accomplished by integrating the model presented here into finite element, reaction/diffusion simulations of getter assembly hydrogen uptake performance. As of Fall 2020, this integration is underway.

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

Development and Experimental Optimization of High-Temperature Modeling Tools and Methods for Concentrated Solar Power Particle - Systems

A novel, open-source radiative modeling toolset was developed to extend the functionality of particle-based modeling software (e.g. discrete element method (DEM)) to environmental conditions relevant to concentrated solar power applications. This toolset was optimized for deployment on desktop workstations instead of high-performance computing systems, to render such tools more accessible to the research community. Both particle-based modeling and radiative exchange modeling are computationally expensive and often require specialized programming expertise, making these methods cumbersome to use. Recent developments in DEM software by DCS Computing have greatly reduced these challenges, providing a graphical-user-interface based platform and modeling optimization for desktop workstations, HPCs, and cloud computing. The University of Dayton leveraged the experience of DCS Computing in developing a user-friendly, open-source radiative heat transfer expansion for DEM modeling. The University of Dayton DEM+ radiative modeling toolset was developed using a combination of fundamental experimental measurements, modeling, and simplified flow experiments over a range of temperatures and flow conditions. The toolset provides researchers with access to multiple radiative models including an accelerated Monte-Carlo Ray Tracing (application agnostic, highly computationally expensive), an expanded database of distance-based approximations (application limited, computationally light), and a weighted blending of the two methods capable of achieving over 90% reduction in computation time with equivalent accuracy compared to Monte-Carlo Ray Tracing. Through a graphical user interface, users can customize the radiative models to match their desired accuracy and available computational resources, improving access to particle based modeling for the research community. Ceramic sintered bauxite proppants were used in modeling and experimentally as a baseline. Both the radiative heat transfer and flow properties for particulate systems were investigated at elevated temperatures up to 800 °C. The major accomplishments for this work include a verified, open-source radiative modeling toolset to be distributed amongst the research community and the fabrication of three small-scale test facilities to investigate particle behavior and tune DEM flow properties for operation up to 800 °C. The findings have been shared with the research community via conference modeling workshops, deployment of the tools in DCS Computing Aspherix®, and open-source access to the developed radiative modeling tool. The development of next-generation CSP facilities and thermal energy storage systems based on ceramic particles requires providing access to computationally efficient and accurate modeling tools. Particles will experience a wide range of environments (20-800 °C) and handling conditions (dilute curtains or dense packing), requiring specially designed and optimized equipment. Optimizing solid particle physics models and establishing best-practices for particle modeling in CSP environments will assist researchers with designing optimized equipment, accelerating the deployment of more economically-competitive CSP facilities.

14 SOLAR ENERGY↗
DOE OSTI2022

Residual Stress Modeling and Advanced Diffraction Measurements of 347H Steel Weldments

The 347H stainless steel is a primary high-temperature material for many energy and power generation industries. Stress relief cracking (SRC) has been a particular concern in welding of this material. The residual stress induced by welding and its evolution during post-welding heat treatment (PWHT) and subsequent operating and service conditions is one of the primary factors contributing to SRC. The lifetime of welded structure components is also controlled by the precipitation kinetics that accompanies PWHT, stress relaxation process, and long-term aging and complex synergistic factors. Various theories have been proposed in the past to explain SRC. However, a widely accepted approach to predicting the entire damage evolution and the resulting performance reduction is still lacking. This study is to demonstrate a reliable solution for two critical issues that affect the predictions. First, the residual stress distribution obtained from both simulation and neutron diffraction is compared, which increases the accuracy of mechanical analysis simulation model and therefore builds a solid basis for the lifetime prediction model. Second, the generation, evolution, and annihilation of precipitates are monitored by the synchrotron diffraction experiment. Preliminary results demonstrate the critical importance of precipitation kinetics on the residual stress distribution/redistribution during heat treatments.

Yang, Yi↗
DOE OSTI2020

Arctic Coastal Erosion: Modeling and Experimentation

Increasing Arctic coastal erosion rates have put critical infrastructure and native communities at risk while also mobilizing ancient organic carbon into modern carbon cycles. Although the Arctic comprises one-third of the global coastline and has some of the fastest eroding coasts, current tools for quantifying permafrost erosion are unable to explain the episodic, storm-driven erosion events. Our approach, mechanistically coupling oceanographic predictions with a terrestrial model to capture the thermo-mechanical dynamics of erosion, enables this much needed treatment of transient erosion events. The Arctic Coastal Erosion Model consists of oceanographic and atmospheric boundary conditions that force a coastal terrestrial permafrost environment in Albany (a multi-physics based finite element model). An oceanographic modeling suite (consisting of WAVEWATCH III, Delft3D-FLOW, and Delft3D-WAVE) produced time-dependent surge and run-up boundary conditions for the terrestrial model. In the terrestrial model, a coupling framework unites the mechanical and thermal aspects of erosion. 3D stress/strain fields develop in response to a plasticity model of the permafrost that is controlled by the frozen water content determined by modeling 3D heat conduction and solid-liquid phase change. This modeling approach enables failure from any allowable deformation (block failure, slumping, etc.). Extensive experimental work has underpinned the ACE Model development including field campaigns to measure in situ ocean and erosion processes, strength properties derived from thermally driven geomechanical experiments, as well as extensive physical composition and geochemical analyses. Combined, this work offers the most comprehensive and physically grounded treatment of Arctic coastal erosion available in the literature. The ACE model and experimental results can be used to inform scientific understanding of coastal erosion processes, contribute to estimates of geochemical and sediment land-to-ocean fluxes, and facilitate infrastructure susceptibility assessments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Solid Oxide Cell Manufacturing Cost Tool

NETL has developed a detailed solid oxide cell (SOC) cost modeling tool. The tool provides a fully customizable spreadsheet that aids in cost estimation of large volume manufacturing of SOC (both fuel cell and electrolysis cell). The tool can be used to assist commercial developers and other end users to evaluate the costs of producing cells and stack components as SOC systems increase in scale. The manufacturing materials, labor, component costs, and specific equipment can be tailored to accurately reflect the desired manufacturing process. The tool is paired with a detailed user manual to assist in the understanding of the full functionality. This document is the spreadsheet tool. To access the associated user guide, please visit: <a href="https://netl.doe.gov/energy-analysis/details?id=f5823130-0d20-4a04-9c53-341682e58ac1" rel="noopener noreferrer">Energy Analysis | netl.doe.gov</a>

cost modeling↗
DOE OSTI2023

Multiple and nonlocal cation redox in Ca–Ce–Ti–Mn oxide perovskites for solar thermochemical applications

Modeling-driven design of redox-active off-stoichiometric oxides for solar thermochemical H 2 production (STCH) seldom has resulted in empirical demonstration of competitive materials. Here, we report the theoretical prediction and experimental evidence that the perovskite Ca 2/3 Ce 1/3 Ti 1/3 Mn 2/3 O 3 is synthesizable with high phase purity, stable, and has desirable redox thermodynamics for STCH, with a predicted average neutral oxygen vacancy (VO) formation energy, E v = 3.30 eV. Flow reactor experiments suggest potentially comparable or greater H 2 production capacity than recent promising Sr–La–Mn–Al and Ba–Ce–Mn metal oxide perovskites. Utilizing quantum-based modeling of a solid solution on both A and B sub-lattices, we predict the impact of nearest-neighbor composition on E v and determine that A-site Ce 4+ reduction dominates the redox-activity of Ca 2/3 Ce 1/3 Ti 1/3 Mn 2/3 O 3 . X-ray absorption spectroscopy measurements provide evidence that supports these predictions and reversible Ce 4+ -to-Ce 3+ reduction. Our models predict that Ce 4+ reduces even when it is not nearest-neighbor to the V O , suggesting that refinement of Ce stoichiometry has the possibility of further enhancing performance.

08 HYDROGEN↗
DOE OSTI2025

New developments and verification of fusion blanket simulation capabilities in the MOOSE framework

Multiphysics modeling capabilities have a crucial role to play in the accelerated deployment of fusion energy. To that end, we developed new multiphysics fusion blanket simulation capabilities in the Multiphysics Object-Oriented Simulation Environment (MOOSE). Firstly, we expanded on the existing capabilities of the previously published work, by coupling 3D tritium transport modeling capabilities using the Tritium Migration Analysis Program, version 8 (TMAP8) to an existing tool including thermal hydraulics, fully three-dimensional (3D) heat transfer, and loosely coupled neutronics analysis. Secondly, we performed a thorough verification of the new capabilities and increased testing code coverage to meet MOOSE’s software quality standards. The MOOSE framework follows a strict software quality assurance plan to be Nuclear Quality Assurance, Level 1 compliant. The new multiphysics fusion blanket simulation capabilities are now held to the same standard. Thirdly, to demonstrate MOOSE’s new fusion blanket modeling capabilities, we performed a fully integrated, multiphysics simulation of a 3D solid ceramic breeder blanket design. This proof-of-concept simulation provides the temperature and tritium distribution across the blanket. In conclusion, the combined efforts towards software quality and the development of multiphysics coupling capabilities provide an effective and reliable framework for modeling solid ceramic fusion blankets using MOOSE.

modeling and simulation↗