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

An unstructured mesh based neutronics optimization workflow

We have developed a fully automated workflow to optimize the neutronics performance of the Second Target Station (STS) at the Oak Ridge National Laboratory’s Spallation Neutron Source. The optimization workflow starts with the parametrized solid CAD engineering models and converts them into the unstructured mesh (UM) models for the neutronics calculations with MCNP6.2. Calculations are executed and their results are loaded into the Dakota optimization toolkit. Dakota analyzes the results and proposes new geometry parameters for the next design iteration. The cycle repeats until the optimal parameters are found. The automated CAD to MCNP conversion, the use of high-fidelity UM models, and the use of modern optimizer are the key elements that advance the entire optimization workflow in comparison with the original workflow. The original workflow was based on a simplified constructive solid geometry (CSG) modeling with MCNPX, mcnp_pstudy tool, and an in-house optimizer. Herein to demonstrate the new workflow, we present a case of neutronics optimization of the moderator–reflector assembly (MRA). Apart from the MRA, the workflow can optimize other major STS components, such as the spallation target, neutron beamlines, radiation shielding, and various accelerator components. Importantly, the new workflow opens the door to the advanced multi-physics multi-parameter optimization and has the potential for use in other nuclear physics and accelerator applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A hybrid-kinetic simulation tool for non-thermal warm x-ray z-pinch sources, with gas-puff and wire array exemplars

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Machine-learning interatomic potentials for interfaces in all-solid-state batteries: Perspectives on training data, model selection, and validation

Interfaces play a pivotal role in dictating the performance and reliability of all-solid-state batteries (ASSBs), where complex electro-chemo-mechanical phenomena at grain boundaries (GBs) and interfaces can lead to degradation and failure. Traditional atomistic simulation methods, such as first-principles calculations and classical molecular dynamics, face limitations in modeling these interfaces due to either high computational cost or insufficient transferability to the diverse atomic environments evolving at interfaces. Machine-learning interatomic potentials (MLIPs) have emerged as a transformative approach, enabling large-scale, high-accuracy simulations of disordered and chemically complex systems by leveraging the predictability of machine learning models trained on first-principles data. Recent applications of MLIPs have demonstrated their ability to capture intricate behaviors at ASSB interfaces, including ion transport, interfacial evolution, and degradation mechanisms, with accuracy and efficiency unattainable by conventional methods. This prospective paper presents comprehensive analysis and practical guidance for MLIP development for GBs and interfaces in ASSBs, with a focus on three key pillars: data generation, model selection, and validation. Here, we review the current state of MLIP applications for GBs and interfaces in both general and ASSB-specific materials, highlighting best practices and challenges in constructing diverse and representative datasets, choosing appropriate machine learning architectures, and rigorously validating model performance. We also discuss emerging strategies and opportunities for improved reliability and efficiency of MLIPs to simulate realistic interfaces in ASSBs.

Energy - Storage↗
DOE OSTI2022

Release of Contaminants from Burning Liquids and Solids

This report details model development, theory, and a literature review focusing on the emission of contaminants on solid substrates in fires. This is the final report from a 2-year Nuclear Safety Research and Development (NSRD) project. The work represents progress towards a goal of having modeling and simulation capabilities that are sufficiently mature and accurate that they can be utilized in place of physical tests for determining safe handling practices. At present, the guidelines for safety are largely empirically based, derived from a survey of existing datasets. This particular report details the development, verification and calibration of a number of code improvements that have been implemented in the SIERRA suite of codes, and the application of those codes to three different experimental scenarios that have been subject of prior tests. The first scenario involves a contaminated PMMA slab, which is exposed to heat. The modeling involved a novel method for simulating the viscous diffusion of the particles in the slab. The second scenario involved a small pool fire of contaminated combustible liquid mimicking historical tests and finds that the release of contaminants has a high functionality with the height of the liquid in the container. The third scenario involves the burning of a contaminated tray of shredded cellulose. A novel release mechanism was formulated based on predicted progress of the decomposition of the cellulose, and while the model was found to result in release that can be tuned to match the experiments, some modifications to the model are desirable to achieve quantitative accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Generating MCNP Input Files for Unstructured Mesh Geometries

The Los Alamos National Laboratory’s (LANL) Monte Carlo N-Particle (MCNP)1 transport code version 6.3 (also known as MCNP6.3) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. This feature has been developed for performing calculations of complex geometry models because manually creating CSG models is time-consuming and error-prone as the complexities of geometries increase. A UM geometry model is a collection of finite elements representing a solid geometry. The first step of the MCNP UM calculation is using other software packages to create a finite element mesh representation of a solid 3D geometry because the MCNP code cannot be used to generate a UM model. Computer-aided design (CAD) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. Some mesh generation software packages may also be used to create solid geometries and thus CAD files are not needed. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software suite. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. Starting with a 6.3 version, the MCNP code can process HDF5 mesh input files. We only focus on the UM models formatted as Abaqus input files in this report since currently no external software can be used to generate HDF5 mesh input files for MCNP UM calculations. The MCNP code version 6.3 can be used to convert the Abaqus mesh input files into the HDF5 mesh input files, but this option is typically used by the MCNP code development team to test the HDF5 mesh input file feature. Several software packages (such as Abaqus, Attila4MC, or Cubit) can be used to create the Abaqus input files for MCNP UM calculations. An MCNP UM calculation using an Abaqus model requires two input file types: MCNP and Abaqus input files. The Abaqus input files needed for MCNP UM calcu lations must have the correct Abaqus syntax and meet the additional requirements by the MCNP code. The MCNP code can process only Abaqus input files that make use of part and assembly definitions, where elements in each part must be grouped into one or more element sets (i.e., elset) using *Elset keyword lines with specified naming formats. The MCNP and Abaqus input files required for MCNP UM simulations must be related; pseudo-cells in an MCNP input file must be constructed from mesh model data from an Abaqus input file. For large complex UM models, it is tedious to manually create MCNP UM input files. The um pre op (unstructured mesh pre operations) program with the -m option can be used to create a skeleton MCNP input file from an Abaqus input file [6]. Since the um pre op program was written in Fortran and was not written for optimized performance, this program is a deprecated feature in the MCNP code version 6.3 and may be removed in the next release of the code. To improve calculation flow of multiphysics calculations, a Python3 code called write mcnp um input has been developed to generate an MCNP input file instead of using the um_pre_op -m option. This Python code was initially released to the public in 2020. We have updated this Python code for MCNP6.3 and it was used to generate the MCNP input files used to verify the MCNP6.3 code. The write_mcnp_um_input code is included with the MCNP6.3 code package which will be released to the public through the Radiation Safety Information Computational Center (RSICC) at Oak Ridge National Laboratory. This report is a revision of LA-UR-20-27139 report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Microstructure-Based Modeling of Inner Oxygen Pressure in Solid Oxide Electrolysis Cells

One major degradation mechanism during long-term operation of solid oxide electrolysis cells (SOECs) is delamination of oxygen electrodes (OEs). The driving force for the electrode delamination could be the generated high inner oxygen pressure near the electrode-electrolyte interface during operation. However, the effects of transport properties and electrode thickness on the inner oxygen partial pressure are not well understood. Here, a microstructure-based electrochemical model which includes the conduction of electrons and oxygen ions coupled with Butler-Volmer-type chemical reactions at triple-phase-boundaries (TPBs), is employed to investigate the oxygen pressure in lanthanum strontium manganate (LSM)-based SOECs. The model is applied to both two-dimensional (2D) prototype microstructures and three-dimensional (3D) realistic microstructures, and the oxygen pressure is analyzed as a function of transport properties and electrode thickness under both potentiostatic and galvanostatic operations. The simulation results suggest strategies to suppress electrode delamination. The simulation results are compared to an analytical solution, and the discrepancies are attributed to the Butler-Volmer-type kinetics included in the microstructure-based model.

Xue, Fei↗
DOE OSTI2021

Coupling particle scale model and SuperDEM–CFD for multiscale simulation of biomass pyrolysis in a packed bed pyrolyzer

An efficient biomass pyrolysis process requires a comprehensive understanding of the chemical and physical phenomena that occur at multi-length and time scales. In this study, a multiscale computational approach was developed and validated for biomass pyrolysis in a packed-bed reactor by integrating pyrolysis kinetics, a particle scale model, and SuperDEM-CFD in open-source code MFiX. A one-dimensional particle-scale model that discretizes the characteristic length of biomass particle into layers was developed to predict the intraparticle phenomena inside a single particle. The 1D model was validated by comparing it with a single biomass particle pyrolysis experiment. Here, a recently developed SuperDEM-CFD model was employed to simulate the non-spherical particle-particle contact and fluid-particle interaction. The coupled model was applied to simulate the pyrolysis of cubic biomass particles in a packed bed and validated by comparing with experimental data. Simulation with and without particle-scale model was compared, and the effect of the gas-solid heat transfer models was also investigated.

09 BIOMASS FUELS↗
DOE OSTI2023

IDAES-PSE 2.3.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications.. New Features and Models • New diagnostics toolboxes and examples o Tutorial for using the DiagnosticsToolbox o Methods to identify potential evaluation errors in models o SVDToolbox for performing singular value decomposition analysis on models to identify possible degeneracies and poor scaling o DegeneracyHunter for identifying irreducible degenerate sets in models • Model for solid-liquid separation which takes two inputs (solid and liquid streams) and produces three outlets (solids, liquid entrained with solids and pure liquid) • Example of temperature swing adsorption in models_extra Deprecation Warnings • With the update to Pyomo 6.7.0, the plate heat exchanger model has begun failing tests on some CI platforms. If this is not addressed by the February release, this model may be removed due to failing tests. See #1294. Offers to help identify the cause and fix this issue are welcome. Improvements and bug fixes • Fixed oversight which precluded using Mixer and Separator models when inherent reactions were present in property package • Added reporting methods to the MSContactor model • Minor corrections in some examples where values were being fixed outside of declared variable bounds

DiagnosticsToolbox,IDAES,PSE,Process Systems Engin↗
DOE OSTI2021

Modeling the performance and faradaic efficiency of solid oxide electrolysis cells using doped barium zirconate perovskite electrolytes

Y-doped BaZrO 3 (BaZr 1 x Y x O 3 δ , or “BZY”), a proton-conducting ceramic featuring high bulk conductivity and good chemical stability, is a promising electrolyte material for solid oxide electrolysis cells. Further doping with Ce and/or Yb (creating materials “BCZY” and “BCZYYb”) can improve conductivity and sintering properties, but at significant penalty to cells’ faradaic efficiency (FE). Studies have proposed that reduction of lattice Ce can occur in the hydrogen electrode, which consumes some hydrogen produced by the hydrogen evolution reaction, leading to decreased FE. Despite studies suggesting this phenomenon, the mechanism is largely unknown. In this work, we developed a multiphysics model to study the transport of multiple defect species and the performance of BZY, BCZY, and BCZYYb, capturing the tradeoff between enhanced performance at the cost of FE for BCZY and BCZYYb electrolytes compared to BZY. We also found that increasing the water content of the anode gas supply lowers the current output of the cell but results in better FE. The model, which uses several parameters previously unavailable in the literature, was validated to experiments varying temperature, steam water content, and electrolyte material, as well as two performance metrics (performance curves and FE). Results verify and explain observed trends, informing future work on Ce-doped BZY electrolytes.

08 HYDROGEN↗
DOE OSTI2022

BDEM (Discrete-element-simulator for high-solids granular flows) [SWR-22-72]

BDEM is a discrete element method based simulation tool developed specifically for modeling high-solids granular flows that include polydispersity, heat-transfer, moving boundaries and chemistry. Our solver provides facilities for simulating spherical/non-spherical particles with modified contact and friction models in complex dynamic geometries defined using level-sets or triangulated files. The solver is developed on top of NREL's open-source performance portable library, AMReX, providing parallel execution capabilities on current and upcoming high-performance-computing (HPC) architectures. Simulations at the scale of several millions to billion particles have been performed using this software on large scale computing resources. This software can be applied to non-reacting solids dominant flows in silos, hoppers and screw conveyors as well as in high temperature reacting systems such as screw kilns and auger reactors.

Sitaraman, Hariswaran↗
DOE OSTI2024

Assessing Critical Conditions for Scour Near Obstructions using Bed Shear, Particle Onset of Motion Balances, and CFD-DEM Modeling of Granular Beds

Computational Fluid Dynamics combined with a Discrete Element Method is one of the computational methods that can be used to model multiphase flows. In this method various phases, gas and liquid or solid, are present in the same computational domain. The local averaged Navier–Stokes equations determine the flow of the continuous phase fluid and are solved using the traditional CFD finite volume approach. DEM is based on a Lagrangian formulation, which solves the equations of motion, expressed in ordinary differential equations, for representative particles as they move in space and time. The interactions between the continuous fluid phase and discrete solid phase are modeled with the use of Newton’s laws of motion via drag force. The particles interact with each other and with the boundaries of the fluid continuum, and the resulting contact forces are included in the equations of motion. The properties of solid particles and boundaries are treated as elastic bodies, with specified density, elastic modulus, and Poisson’s ratio. Particle shapes may vary from single spherical particles to more complex-shaped composite particles. The particles may be introduced into the domain by random or structured injection at a point, surface, or volume, depending on the application. More details on the formulation can be found in the Simcenter STAR-CCM+ User’s Manual and OpenFOAM website.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Modeling thin layers of analytes on substrates for spectral analysis: use of solid/liquid n and k values to model reflectance spectra

Since solids are only sometimes seen en masse in a pure bulk form, and for liquids other than water almost never, a capability to model reflectance spectra from analytes deposited on various substrates would be highly advantageous. If available, the real, n ( ν ) , and imaginary, k ( ν ) , components of the complex refractive index, n n + i k , can be used to simulate infrared spectra, accounting for reflection, refraction, and absorption phenomena as a function of wavelength. We focus on using the Pacific Northwest National Laboratory (PNNL) derived n / k vectors for solid and liquid analytes deposited as thin layers on different types of substrates including conductors, such as aluminum, and inorganic dielectrics, such as glass. The model is an adaptation of the Monte Carlo ray trace modeling program, TracePro, extended through the use of its macrolanguage. The model is tested using thin films of organic liquids including silicone oil and no. 2 diesel fuel, as well as organic solids such as caffeine and acetaminophen on aluminum and glass. The predicted spectra for the solid films were compared to experimental hemispherical reflectance data measured using a Fourier transform spectrometer with an integrating sphere. The thickness of the calculated layer is a parameter for predicting the (transflectance) spectra and is obtained using the areal density measured from gravimetric methods to generate the thin-layer samples. Comparison of the calculated spectra with experimental hemispherical reflectance data shows excellent agreement, indicating promise for the use of measured n / k data to synthesize reference spectral data. In particular, accounting for the inhomogeneity of the deposits greatly improved the match with experimental data. Finally, the theoretical modeling shows that for thicker layers (ca. 20 to 100 μm) of typical organics possessing moderately strong k values, the longwave infrared features are often saturated and better spectral contrast is obtained from the overtone/combination bands in the shortwave infrared.

42 ENGINEERING↗
DOE OSTI2025

Thermodynamic modeling of countercurrent chemical looping reverse water gas shift process for redox material screening

The reverse water gas shift (RWGS) reaction is a key pathway for CO 2 utilization, particularly within Power-to-X process chains aimed at sustainable fuel and chemical production. Countercurrent chemical looping (CL-RWGS) using non-stoichiometric oxides can overcome equilibrium limitations of conventional RWGS reactors, enabling significantly higher CO 2 conversions. However, modeling the limiting performance of such systems is challenging due to their multiphase nature and coupled spatial and temporal variation in chemical composition. In this work, we present a discretized batch equilibrium model that simulates CL-RWGS reactors as a series of localized equilibrium exchanges between gas and solid elements. The model is numerically stable, computationally efficient, and free of kinetic source terms, making it well-suited for parametric studies and system-level integration. It is validated against established convection–diffusion models and shown to predict reasonable upper bounds on experimental results. Application of the model to a range of oxygen carrier materials identifies cerium–zirconium solid solutions, particularly Ce 0.80 Zr 0.20 O 2 , as a promising class offering superior oxygen storage characteristics compared to state-of-the-art La 0.6 Sr 0.4 FeO 3 . This framework provides a robust platform for materials screening, reactor sizing, and performance optimization in chemical looping systems. The model implementation is available as open-source software to support further research and development.

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

Efficient Reformulation of Linear and Nonlinear Solid-Phase Diffusion in Lithium-ion Battery Models using Symmetric Polynomials: Mass Conservation and Computational Efficiency

Lithium-ion batteries are typically modeled using porous electrode theory coupled with various transport and reaction mechanisms, along with suitable discretization or approximations for the solid-phase diffusion equation. The solid-phase diffusion equation represents the main computational burden for typical pseudo-2-dimensional (p2D) models since these equations in the pseudo r -dimension must be solved at each point in the computational grid. This substantially increases the complexity of the model as well as the computational time. Traditional approaches towards simplifying solid-phase diffusion possess certain significant limitations, especially in modeling emerging electrode materials which involve phase changes and variable diffusivities. A computationally efficient representation for solid-phase diffusion is discussed in this paper based on symmetric polynomials using Orthogonal Collocation and Galerkin formulation (weak form). A systematic approach is provided to increase the accuracy of the approximation (p form in finite element methods) to enable efficient simulation with a minimal number of semi-discretized equations, ensuring mass conservation even for non-linear diffusion problems involving variable diffusivities. These methods are then demonstrated by incorporation into the full p2D model, illustrating their advantages in simulating high C-rates and short-time dynamic operation of Lithium-ion batteries.

25 ENERGY STORAGE↗
DOE OSTI2024

Phase Field Modeling of Pressure Induced Densification in Solid Electrolytes

Adoption of dense and homogeneous solid electrolytes can possibly mitigate the propagation of lithium dendrites and enable lithium metal anodes. Application of external pressure helps to minimize the sintering temperature in oxide ceramics and can potentially densify softer sulfide electrolytes even under room temperature conditions. Here, a previously developed phase field-based computational scheme for predicting the high-temperature sintering-induced densification of oxide ceramic solid electrolytes is extended in the present context to capture the influence of external pressure for densifying solid electrolytes. Two different bulk deformation mechanisms, namely, "reorganization" and "creep deformation," are dominant under external pressure, which is different from the surface and grain-boundary diffusion-induced densification of solid electrolytes that occurs during high temperature sintering. External pressure also increases the points of contact between the particles, which further enhances the propensity of diffusion-induced sintering process. Results obtained from simulations indicate that densification under external pressure is independent of the solid electrolyte particle morphology. Finally, a phase map is generated between applied pressure and temperature for achieving complete densification of oxide ceramics, which can possibly guide the synthesis of thin and dense solid electrolyte separators.

25 ENERGY STORAGE↗