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

Modeling Ni Coarsening and Migration in the Hydrogen Electrode of Solid oxide Cells under Operating Conditions

Ni coarsening and migration is the most important degradation in the hydrogen electrode of solid oxide cells. This presentation reviews our works on Ni coarsening and migration in NETL over the past years. We used simulation techniques including phase-field modeling and density function theory to investigate possible mechanisms of Ni coarsening and migration, including self-diffusion of Ni, diffusion of gaseous Ni(OH)2, surface diffusion of Ni(OH)x and Ni-YSZ wettability change, and examined the effect of operating conditions on these mechanisms. So far, none of the mechanisms can fully explain the experiments. Remaining questions and possible paths forward are summarized and discussed.

Lei, Yinkai↗
DOE OSTI2022

Scale-Up and Modeling Efforts Using an Omniphobic Surface Treatment for Mitigating Solids Deposition

Abstract Gas hydrates, waxes and asphaltenes represent some of the most significant flow assurance challenges, especially in subsea lines, where treatment options can be limited. Currently, complete avoidance is the primary strategy for hydrate management, while chemical or mechanical flushing/pigging may be utilized for other solids. Each system must also be treated with individualized solutions, as there have been no proved one-size-fits-all technologies demonstrated to date. As an alternative to constant chemical injection or thermodynamic controls such as insulation or heating, a robust omniphobic surface treatment material has been developed which has been shown in previous studies to significantly reduce the adhesion of flow assurance solids, resulting in lower risk for deposition and plugging of gas hydrates, waxes, and asphaltenes. As part of a Department of Energy study, laboratory scale tests were performed on a variety of apparatuses ranging from micromechanical force (micron scale) to laboratory flow loop testing (meter scale). While the results from these tests have been promising, there can be some disconnects between lab-scale observations and field-scale testing. Several factors can result in discrepancies between small-scale bench-top studies and deployment-scale efforts. These include flow conditions (Reynolds number, flow regime), chemistry (synthetic vs. crude oils, gas composition, and additives), physical differences (pipeline material, size, and geometry), etc. This work will describe the actions taken to better understand the deployment and performance of a material such as this omniphobic surface treatment in field scenarios, and to build confidence for deployment in a production scenario. These efforts include modeling field production scenarios using models developed from lab observations, application of the material to field-scale equipment to better understand application challenges, and expanded survivability and longevity testing. Understanding the bridge between laboratory testing and full-scale deployment allows for better technique and risk mitigation for field scale testing.

03 NATURAL GAS↗
DOE OSTI2025

A Perspective on Multiscale Modeling of Explicit Solvation-Enabled Simulations of Catalysis at Liquid–Solid Interfaces

Catalysis at liquid-solid interfaces is profoundly influenced by the interfacial solvent structure, which affects catalytic activity, selectivity, and reaction pathways. This perspective discusses state-of-the-art multiscale modeling methods that integrate quantum mechanics and molecular mechanics approaches to apply explicit solvent molecules to capture these interfacial phenomena. Specifically, the construction of multiscale models, the importance of capturing the interfacial solvent structure, and the computational strategies used to achieve this are explored, and the challenges in balancing chemical accuracy with computational expense are highlighted. Additionally, this perspective addresses the limitations of current methods. Opportunities for integrating machine learning are proposed. Here, by advancing the efficiency and user friendliness of multiscale modeling, it is argued that deeper insights into heterogeneous catalysis in liquid phases can be provided, which will ultimately contribute to the development of more efficient catalytic processes.

Ab initio molecular dynamics↗
DOE OSTI2021

Local Multimodal Electro-Chemical-Structural Characterization of Solid-Electrolyte Grain Boundaries

Typical models of polycrystalline ionic materials treat the grain boundary properties as single valued, without consideration of the full range of values that define the macroscopically measured average. Here a unique experimental platform suitable for local multimodal characterization of individual grain boundaries in bicrystal fibers is reported. A variation of three orders of magnitude in the grain boundary conductivity of ceria is observed, as measured across six individual bicrystals by both alternating current impedance spectroscopy and direct current (D.C.) linear sweep voltammetry. Nonlinear behavior of the D.C. measurements is consistent with resistance due to a space charge effect. Time-of-flight secondary ion beam spectroscopy reveals a correlation between grain boundary resistance and the concentration of impurities Si, Al, and Ca segregated at the grain boundaries, although the bulk concentrations of these impurities are negligible. Electron backscatter diffraction analysis of the crystal orientations suggests a correlation between the misorientation across the grain boundaries and grain boundary resistance. These correlations point towards a grain boundary resistance that arises from impurity-generated space charge effects and variations in impurity concentration and hence resistivity driven by the energetics of impurity segregation to grain boundaries of differing surface energies.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Numerically exact generalized Green's function cluster expansions for electron-phonon problems

We generalize the family of approximate momentum average methods to formulate a numerically exact, convergent hierarchy of equations whose solution provides an efficient algorithm to compute the Green's function of a particle dressed by bosons suitable in the entire parameter regime. We use this approach to extract ground-state properties and spectral functions. Our approximation-free framework, dubbed the generalized Green's function cluster expansion (GGCE), allows access to exact numerical results in the extreme adiabatic limit, where many standard methods struggle or completely fail. We showcase the performance of the method, specializing three important models of charge-boson coupling in solids and molecular complexes: the molecular Holstein model, which describes coupling between charge density and local distortions, the Peierls model, which describes modulation of charge hopping due to intersite distortions, and a more complex Holstein + Peierls system with couplings to two different phonon modes, paradigmatic of charge-lattice interactions in organic crystals. Furthermore, the GGCE serves as an efficient approach that can be systematically extended to different physical scenarios, thus providing a tool to model the frequency dependence of dressed particles in realistic settings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Cell and Stack Degradation Evaluation and Modeling

Presentation on NETL's solid oxide fuel cell (SOFC) work plan research (FWP 1022411) given by invitation to the 2024 Hydrogen Program Annual Merit Review meeting in Arlington, VA on May 7, 2024.

Abernathy, Harry↗
DOE OSTI2020

Evolution of solid electrolyte interphase and active material in the silicon wafer model system

Alloying anode materials for lithium-ion batteries, such as silicon (Si), are an important focus of materials research due to the demand for increased energy density for electric vehicles and stationary grid storage. However, progress towards next generation anode materials has been hindered by poor capacity retention due to the instability of the Si active material and the solid electrolyte interphase (SEI). Herein, the evolution of the active Si material and the SEI are simultaneously investigated from the perspective of chemical, structural, morphological, and electronic evolution in the Si wafer model system through its cycling life, using time-of-flight secondary ion mass spectrometry (TOF-SIMS), scanning transmission electron microscopy (STEM), atomic force microscopy (AFM), and scanning spreading resistance microscopy (SSRM). The results illustrate a dynamic evolution of the SEI and active Si material through cycling. Through an improved understanding of evolution of SEI and Si active material within the Si wafer model system, mitigation strategies may be designed to extend the lifetime of Si anodes.

25 ENERGY STORAGE↗
DOE OSTI2024

A layered solid finite element formulation with interlaminar enhanced displacements for the modeling of laminated composite structures

Accurate modeling of layered composite structures often requires the use of detailed finite element models which can sufficiently resolve the kinematics and material behavior within each layer of the composite. However, individually discretizing each material layer into finite elements presents a prohibitive computational expensive given the large number of thin layers comprising some laminated composites. To address these challenges, an 8-node layered solid hexahedral finite element is formulated with the aim of striking an appropriate balance between efficiency and fidelity. The element is discretized into an arbitrary number of distinct material layers, and employs reduced in-plane integration within each layer. The chosen reduced integration scheme is supplemented by a novel physical stabilization approach which includes layerwise enhancements to mitigate various forms of locking phenomena. The proposed framework additionally supports the inclusion of interlaminar enhanced displacements to better represent the kinematics of general layered composite materials. Finally, the described element formulation has been implemented in the ParaDyn finite element code, and its efficacy for modeling laminated composite structures is demonstrated on a variety of verification problems.

42 ENGINEERING↗
DOE OSTI2021

OpenMC

OpenMC is an open source, community-developed framework for performing Monte Carlo particle transport calculations.It is capable of simulating neutron and photon transport below ~50 MeV and also includes a model for bremsstrahlung production of photons by electrons/positrons. Fixed source, k-eigenvalue, and subcritical multiplication problems can be solved. Problem geometries can be modeled either using constructive solid geometry or a CAD representation. A flexible and efficient tally system enables a wide variety of physical quantities to be tallied and analyzed. OpenMC can run in parallel using a hybrid MPI and OpenMP programming model and has been extensively tested on leadership class supercomputers. A built-in Bateman equation solver enables the modeling of material composition changes due to irradiation. One of the unique features of OpenMC is its rich, extensible Python and C/C++ programming interfaces that enable programming pre- and post-processing, multigroup cross section generation, workflow automation, multiphysics coupling, and the visualization of geometry and tally results. In addition to the core Monte Carlo transport solver and associated APIs, OpenMC includes a Python-based nuclear data interface that enables power users to inspect, modify,and perform various types of nuclear data processing on ENDF, ACE, and OpenMC¿s native HDF5 files

ROMANO, PaulK.↗
DOE OSTI2020

Evidence for a Solid-Electrolyte Inductive Effect in the Superionic Conductor Li 10 Ge 1– x Sn x P 2 S 12

Strategies to enhance ionic conductivities in solid electrolytes typically focus on the effects of modifying their crystal structures or of tuning mobile-ion stoichiometries. A less-explored approach is to modulate the chemical bonding interactions within a material to promote fast lithium-ion diffusion. Recently, the idea of a solid-electrolyte inductive effect has been proposed, whereby changes in bonding within the solid-electrolyte host framework modify the potential energy landscape for the mobile ions, resulting in an enhanced ionic conductivity. Direct evidence for a solid-electrolyte inductive effect, however, is lacking—in part because of the challenge of quantifying changes in local bonding interactions within a solid-electrolyte host framework. Here, we consider the evidence for a solid-electrolyte inductive effect in the archetypal superionic lithium-ion conductor Li 10 Ge 1–x Sn x P 2 S 12 . Substituting Ge for Sn weakens the {Ge,Sn}–S bonding interactions and increases the charge density associated with the S 2– ions. This charge redistribution modifies the Li + substructure causing Li + ions to bind more strongly to the host framework S 2– anions, which in turn modulates the Li + ion potential energy surface, increasing local barriers for Li + ion diffusion. Each of these effects is consistent with the predictions of the solid-electrolyte inductive effect model. Density functional theory calculations predict that this inductive effect occurs even in the absence of changes to the host framework geometry due to Ge → Sn substitution. These results provide direct evidence in support of a measurable solid–electrolyte inductive effect and demonstrate its application as a practical strategy for tuning ionic conductivities in superionic lithium-ion conductors.

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

Towards a NEAMS-based high-fidelity model of the MARVEL reactor

This report outlines the progress of Idaho National Laboratory in developing a high-fidelity and high-resolution model of the Microreactor Applications Research Validation and Evaluation reactor. The model was developed under the Nuclear Energy Advanced Modeling and Simulation microreactor application driver at Idaho National Laboratory. The overarching objective of this activity is the development of a high-fidelity multiphysics MARVEL model using NEAMS tools, and to verify and validate NEAMS tools against MARVEL reference simulation and experimental data, respectively. This is a unique opportunity to conduct multiphysics analysis on a soon-to-be-deployed microreactor. This multiphysics model developed under the NEAMS-funded INL microreactor application driver leverages three single-physics models coupled via the MOOSE’s MultiApp and Transfer systems. The latter systems enable in-memory data transfer between MOOSE-based and MOOSE-wrapped applications. The first single-physics model, that functions as main application, leverages Griffin to model the neutron transport in the core through the discontinuous finite element (DFEM) discrete ordinates solver (SN). Several optimization flags that were developed by the Griffin developer team were beta-tested to enhance the solver’s performance. These include the combined use of using_average_xs and update_averaged_xs_on that enable to avoid expensive on-the-fly cross sections evaluations at each linear iterations in favor of evaluations of the macroscopic cross sections at each Picard iteration. The second single-physics model uses BISON to handle solid heat transfer and asymptotic hydrogen redistribution analysis in the fuel. While the model returns consistent results for the temperature and hydrogen distribution in the fuel, a mismatch was noticed in the calculated temperature in the reflector due to the value of the gap conductance used in our model. Ongoing investigations are being performed to assess the origin of this discrepancy. Finally, the System Analysis Module (SAM) was used to model the flow of the sodium-potassium eutectic in the primary loop. A first verification was also performed showing good agreement in terms of mass flow rate and inlet temperature. All mesh files were generated using the MOOSE Reactor module, removing the need for external meshing tools. Notably, this workscope represents one of the initial applications of the MOOSE Reactor module for modeling highly irregular geometries. The use of the reactor module significantly streamlined the mesh generation process. The full multiphysics mode, that combines all the single physics models, was leveraged to conduct initial steady-state multiphysics simulations to compute power, and temperature distribution in the reactor. Initial testing was performed for transient simulations as well. In this case, the new checkpoint restart capability for eigenvalue calculations was tested showing the capability for streamlined restart of transient calculations. Future work will focus on improving the fidelity of the model by performing comprehensive code-to-code comparisons. For instance, the full-core Griffin neutronics model will be benchmarked against MCNP reference results, that were provided by the MARVEL design team. Additionally, the SAM T/H model will be verified against reference RELAP-5 results for selected accident scenarios. Besides code-to-code verification exercises, the model fidelity will be improved by replacing the single-channel SAM model with a more complex SAM-Pronghorn coupled model, in which the sub-channel capability is deployed to obtain radial temperature resolution in the coolant. This model will be developed in synergy with the NEAMS thermal hydraulics team.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Thermal Pressure in the Thermal Equation of State for Solid and a Proposed Substitute

Abstract The thermal equation of state (TEOS) for solids is a mathematic model among pressure, temperature and density, and is essential for geophysical, geochemical, and other high pressure–temperature (high P–T) researches. However, in the last few decades, there has been a growing concern about the accuracy of the pressure scales of the calibrants, and efforts have been made to improve it by either introducing a reference standard or building new thermal pressure models. The existing thermal equation of state, P ( V , T ) = P ( V , T 0 ) + P th ( V , T ), consists of an isothermal compression and an isochoric heating, while the thermal pressure is the pressure change in the isochoric heating. In this paper, we demonstrate that, for solids in a soft pressure medium in a diamond anvil cell, the thermal pressure can neither be determined from a single heating process, nor from the thermal pressure of its calibrant. To avoid the thermal pressure, we propose to replace the thermal pressure with a well-known thermal expansion model, and integrate it with the isothermal compression model to yields a Birch–Murnaghan-expansion TEOS model, called VPT TEOS. The predicted pressure of MgO and Au at ambient pressure from Birch–Murnaghan-expansion VPT TEOS model matches the experimental pressure of zero (0) GPa very well, while the pressure prediction from the approximated Anderson PVT TEOS exhibit a big deviation and a wrong trend.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Computational quantum chemistry of metal–organic frameworks

Metal–organic frameworks (MOFs) have premium exceptional properties for a variety of functions, such as gas separation and storage and catalysis. The large variety of possible inorganometallic nodes and organic linkers provide an almost unlimited number of combinations for assembling MOFs, which makes the experimental characterization and examination of all potentially useful combinations practically impossible. Furthermore, experimental studies of MOFs typically fall short in uncovering crucial details regarding their mechanisms of action or the molecular details responsible for their functional properties, such as the nature of adsorbate binding or the structures of transition states. Computational modeling has, therefore, become an efficient and important tool for strategizing the functionalization of MOFs and explicating the mechanisms of their functions. Furthermore, we review the computational methodologies used for computational studies of MOFs, especially Kohn–Sham density functional theory and combined quantum mechanical and molecular mechanical methods for calculating their structural, electronic, and magnetic properties, as well as for understanding the mechanisms of MOFs' applications to magetic devices, thermal conduction, gas adsorption, separation, storage, and sensing, thermal catalysis, photocatalysis, and electrocatalysis.

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

Minimum conditions for accurate modeling of urea production via co-electrolysis

Co-electrolysis of carbon oxides and nitrogen oxides promise to simultaneously help restore the balance of the C and N cycles while producing valuable chemicals such as urea. However, co-electrolysis processes are still largely inefficient and numerous knowledge voids persist. Here, we provide a solid thermodynamic basis for modelling urea production via co-electrolysis. First, we determine the energetics of aqueous urea produced under electrochemical conditions based on experimental data, which enables an accurate assessment of equilibrium potentials and overpotentials. Next, we use density functional theory (DFT) calculations to model various co-electrolysis reactions producing urea. The calculated reaction free energies deviate significantly from experimental values for well-known GGA, meta-GGA and hybrid functionals. These deviations stem from errors in the DFT-calculated energies of molecular reactants and products. In particular, the error for urea is approximately -0.25 ± 0.10 eV. Finally, we show that all these errors introduce large inconsistencies in the calculated free-energy diagrams of urea production via co-electrolysis, such that gas-phase corrections are strongly advised.

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

Characterization and Accelerated Life Testing of a New Solid Oxide Electrolysis Cell

The project aims to develop mechanistic degradation models that realistically predict long-term Solid Oxide Electrolysis Cell (SOEC) durability, using input data from accelerated electrochemical life testing combined with quantitative microstructural and microchemical evaluation. A promising SOEC type will be further developed. The understanding achieved by combining experimental results and theory will be used to guide improvements in long-term SOEC durability, and validate that this technology can reach DOE H 2 production cost and durability targets. Planned outcomes include: (1) new experimental accelerated testing data on SOEC degradation at the oxygen electrode, fuel electrode, and electrolyte; (2) predictive theory of degradation developed based on the accelerated test data; (3) new electrode materials and processing methods that minimize degradation rate; (4) determination of the impact of SOEC electrolyte thickness and electrode composition/microstructure on performance and degradation; (5) determination of the operating conditions where degradation is minimized; and (6) validation of degradation models for extrapolating accelerated test data to predict long-term durability.

08 HYDROGEN↗
DOE OSTI2020

Crystallization histories of the group IIF iron meteorites and Eagle Station pallasites

The group IIF iron meteorites and Eagle Station pallasites (PES) have highly siderophile element abundances (HSE; Re, Os, Ir, Ru, Pt, and Pd) of metal that are consistent with formation in planetesimal cores by fractional crystallization with minor to major solid metal-liquid metal mixing. Modeling of HSE abundances of the IIF irons indicates a complex formation history that included the mixing of primitive and evolved solid and liquid metals. By contrast, modeling of HSE abundances of PES metal suggests these meteorites formed mainly as equilibrium solids from a common liquid. Abundances of some of the siderophile elements in the IIF irons and PES are permissive of a common core origin, however, the abundances of W and Ni indicate the PES ultimately formed on a more oxidized body. The PES most likely formed by the injection of olivine present at the core-mantle boundary into a metallic core liquid as a result of impact. The core then crystallized inward, trapping the olivine.

58 GEOSCIENCES↗
DOE OSTI2021

Peridynamic Model for Microballistic Perforation of Multilayer Graphene

The peridynamic theory of solid mechanics is applied to the continuum modeling of the impact of small, high-velocity silica spheres on multilayer graphene targets. The model treats the laminate as a brittle elastic membrane. The material model includes separate failure criteria for the initial rupture of the membrane and for propagating cracks. Material variability is incorporated by assigning random variations in elastic properties within Voronoi cells. Here, the computational model is shown to reproduce the primary aspects of the response observed in experiments, including the growth of a family of radial cracks from the point of impact.

42 ENGINEERING↗