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

Declining carbohydrate solubilization with increasing solids loading during fermentation of cellulosic feedstocks by Clostridium thermocellum: documentation and diagnostic tests

Abstract Background For economically viable 2nd generation biofuels, processing of high solid lignocellulosic substrate concentrations is a necessity. The cellulolytic thermophilic anaerobe Clostridium thermocellum is one of the most effective biocatalysts for solubilization of carbohydrate harbored in lignocellulose. This study aims to document the solubilization performance of Clostridium thermocellum at increasing solids concentrations for two lignocellulosic feedstocks, corn stover and switchgrass, and explore potential effectors of solubilization performance. Results Monocultures of Clostridium thermocellum demonstrated high levels of carbohydrate solubilization for both unpretreated corn stover and switchgrass. However, fractional carbohydrate solubilization decreases with increasing solid loadings. Fermentation of model insoluble substrate (cellulose) in the presence of high solids lignocellulosic spent broth is temporarily affected but not model soluble substrate (cellobiose) fermentations. Mid-fermentation addition of cells ( C. thermocellum ) or model substrates did not significantly enhance overall corn stover solubilization loaded at 80 g/L, however cultures utilized the model substrates in the presence of high concentrations of corn stover. An increase in corn stover solubilization was observed when water was added, effectively diluting the solids concentration mid-fermentation. Introduction of a hemicellulose-utilizing coculture partner, Thermoanaerobacterium thermosaccharolyticum, increased the fractional carbohydrate solubilization at both high and low solid loadings. Residual solubilized carbohydrates diminished significantly in the presence of T. thermosaccharolyticum compared to monocultures of C. thermocellum, yet a small fraction of solubilized oligosaccharides of both C 5 and C 6 sugars remained unutilized. Conclusion Diminishing fractional carbohydrate solubilization with increasing substrate loading was observed for C. thermocellum -mediated solubilization and fermentation of unpretreated lignocellulose feedstocks. Results of experiments involving spent broth addition do not support a major role for inhibitors present in the liquid phase. Mid-fermentation addition experiments confirm that C. thermocellum and its enzymes remain capable of converting model substrates during the middle of high solids lignocellulose fermentation. An increase in fractional carbohydrate solubilization was made possible by (1) mid-fermentation solid loading dilutions and (2) coculturing C. thermocellum with T. thermosaccharolyticum , which ferments solubilized hemicellulose. Incomplete utilization of solubilized carbohydrates suggests that a small fraction of the carbohydrates is unaffected by the extracellular carbohydrate-active enzymes present in the culture.

09 BIOMASS FUELS↗
DOE OSTI2022

Simulation of the High Temperature Test Facility (HTTF) Core Using the 2D Ring Model with SAM

A two-dimensional ring model is developed with SAM to model the core of the High Temperature Test Facility (HTTF) at the system level. The ring model simplifies the complex structure of the HTTF core by converting the hexagonal rows of heaters and flow channels into layers of concentric annular rings. The ring model is first compared against a 3D-1D model where the solid structures are fully resolved in three dimensions while the fluid structures are modeled as one-dimensional flows. Comparison between the 3D-1D and the ring models shows that the latter can predict major parameters reasonably well under steady-state normal operating condition, but the heater temperatures are underpredicted. Adjustment is made to the effective thermal conductivity of the ceramic core of the ring model to improve the heater temperature predictions. Furthermore, the ring model is also used to simulate a transient Pressurized Conduction Cooldown (PCC) condition and is benchmarked with the experimental data from the HTTF Test PG-27. Good agreement is obtained between the experimental data and the predictions by the ring model.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Semi-Annual Report for Horizontal Compact High Temperature Gas Reactor (HC-HTGR) Development during Performance Period April 2023 – September 2023

Horizontal Compact High Temperature Gas Reactor (HC-HTGR) is being designed by a multi-disciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed HC-HTGR in 3 years and support its commercialization as a safe, low-cost HTGR. Argonne National Laboratory (Argonne) is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. Additionally, Argonne is providing analysis of the primary coolant system to ensure temperatures within the core remain below safety margins during steady-state and potential accident scenarios. This fourth semi-annual report summarized the progress made at Argonne on the two tasks during the second half of FY23. As a part of the RCCS design task, recent efforts have been made to complete a conceptual design of the RCCS for the HC-HTGR, including the design update of the water panel and system configuration favorable in point of view of fabrication and system operation. Design calculations were conducted under various heat load conditions to validate the system design. Transient simulations using RELAP5-3D were conducted to investigate system dynamics under transients of interest and to evaluate the system performance in the design condition. The results demonstrated the overall system feasibility that the RCCS design maintains structures temperatures lower than maximum allowable temperature with sufficient system inventory without any active heat sink. In the primary system thermal hydraulics task, the preliminary analysis was performed for a long term pressurized conduction cooldown (PCC) transient. This analysis used a combination of a fully resolved and coarse homogenized mesh to predict the temperature distribution for the steady-state initial condition and the PCC transient. The steady-state initial condition was determined using a fully resolved full core model with 3D solid to 1D fluid coupling. The fully resolved mesh was also used to model the first 20 seconds of the PCC. The temperature difference between fuel pins and the graphite matrix becomes minimal and the dominant heat transfer shifts to a larger scale radially towards the RCCS. After 20 seconds, a homogenized coarse mesh is used, greatly reducing the computational costs of the model. These results demonstrate that the core is designed to passively remove enough decay heat in a protected loss of primary coolant flow to prevent an unsafe rise of core temperatures.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Understanding the Solid-State Electrode–Electrolyte Interface of a Model System Using First-Principles Statistical Mechanics and Thin-Film X-ray Characterization

Intermixing of atomic species at the electrode-electrolyte boundaries can impact the properties of the interfaces in solid-state batteries. Herein, this work uses first-principles statistical mechanics along with experimental characterization to understand intermixing at the electrode-electrolyte interface. For the model presented in this work, lithium manganese oxide (LiMn 2 O 4 , LMO) and lithium lanthanum titanate (Li 3x La 2/3-x TiO 3 , LLTO) are employed as the cathode and electrolyte, respectively. The results of the computational work show that Ti-Mn intermixing at the interface is significant at synthesis temperatures. The experimental results in this work find that, at some critical temperatures between 600 and 700 degrees C for material preparation, the interface of LLTO-LMO becomes blurred. Calculations predict that the interface is unstable with regard to Ti-Mn intermixing starting at 0 K, suggesting that the critical temperature found in the experiment is related to kinetics. The work overall suggests that, in designing a solid-state battery, the fundamental reactions such as intermixing need to be considered.

25 ENERGY STORAGE↗
DOE OSTI2024

Modeling and Analysis of Polarization Losses in Solid Oxide Fuel Cells with Siloxane Contamination

In this study, the degradation of the solid-oxide fuel cell (SOFC) nickel-yttria stabilized zirconia anode under decamethyltetrasiloxane (L4) contamination is examined with experiments and modeling. A model is developed for the polarization losses based on the charge transfer coefficient,α, and diffusion layer thickness,δ, and fitted to the experimental data to understand how the siloxane degrades the SOFC performance with time. The results of the model indicate that the total polarization losses increase approximately 44% over the course of the 180 min experiment at 350 mA cm −2 . Activation losses dominate the polarization losses initially but decrease in their total contribution while concentration losses increase. Scanning electron microscopy (SEM) with wavelength dispersive X-ray spectroscopy (WDS) elemental mapping indicates that silicon deposition is highest at the outer edge of the anode and forms a barrier layer to fuel diffusion, increasing concentration losses. When the model is applied to other previous D4 and L4 siloxane experiments conducted over a period of 40 h, similar trends in polarization losses are observed. Polarization losses increase more rapidly with D4 compared to L4 siloxane contamination, with concentration losses increasing the fastest with both types of siloxane.

Electrochemistry↗
DOE OSTI2022

Explicit modeling of pebble temperature in the porous-media model for pebble-bed reactors

In this study, we developed a multiscale model to include an explicit pebble-temperature model nested in the porous-media model for pebble-bed reactor applications. The multiscale solid-phase energy balance model, including the pebble surface energy balance equation and an explicit modeling of pebble temperature, can predict the macroscopic (pebble bed) and microscopic (pebble) temperature distributions under both steady-state and transient conditions. The proposed multiscale model is solved in a fully coupled manner using the Newton- Krylov method, and therefore iterations between the macroscopic (pebble-bed-scale) and microscopic (pebble- scale) model are avoided. Extensive code verifications, validation, and demonstrations have been performed for this newly developed model. By explicitly modeling pebble temperatures, this new model addresses a major deficiency of the basic porous-media model, which assumes homogeneous solid-phase temperature and is not appropriate for pebble-bed reactor design and safety analyses.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

Electronic coarse-graining of long conjugated molecules: Case study of non-fullerene acceptors

By considering only one electronic state per molecule, charge transport models of molecular solids neglect intramolecular charge transfer. This approximation excludes materials with quasi-degenerate spatially separated frontier orbitals, such as non-fullerene acceptors (NFAs) and symmetric thermally activated delayed fluorescence emitters. By analyzing the electronic structure of room-temperature molecular conformers of a prototypical NFA, ITIC-4F, we conclude that the electron is localized on one of the two acceptor blocks with the mean intramolecular transfer integral of 120 meV, which is comparable with intermolecular couplings. Therefore, the minimal basis for acceptor–donor–acceptor (A–D–A) molecules consists of two molecular orbitals localized on the acceptor blocks. This basis is robust even with respect to geometry distortions in an amorphous solid, in contrast to the basis of two lowest unoccupied canonical molecular orbitals withstanding only thermal fluctuations in a crystal. The charge carrier mobility can be underestimated by a factor of two when using single site approximation for A–D–A molecules in their typical crystalline packings.

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

Combined Data and Deep Learning Model Uncertainties: An Application to the Measurement of Solid Fuel Regression Rate

In complex physical process characterization, such as the measurement of the regression rate for solid hybrid rocket fuels, where both the observation data and the model used have uncertainties originating from multiple sources, combining these in a systematic way for quantities of interest (QoI) remains a challenge. In this paper, we present a forward propagation uncertainty quantification (UQ) process to produce a probabilistic distribution for the observed regression rate r. We characterized two input data uncertainty sources from the experiment (the distortion from the camera U c and the non-zero-angle fuel placement U Y ), the prediction and model form uncertainty from the deep neural network (U m ), as well as the variability from the manually segmented images used for training it (U s ). Here, we conducted seven case studies on combinations of these uncertainty sources with the model form uncertainty. The main contribution of this paper is the investigation and inclusion of the experimental image data uncertainties involved, and how to include them in a workflow when the QoI is the result of multiple sequential processes.

42 ENGINEERING↗
DOE OSTI2025

Thermal Hydraulics Analysis of a Divertor Monoblock Using SALAMANDER

Divertors are critical components in magnetic confinement fusion devices. One of the Divertor's crictial roles is absorbing the highest heat flux from the plasma. The Divertor accomplishes this through the use of divertor monoblocks and cooling channels. This work presents a thermal hydraulic analysis of a divertor monoblock using the Software for Advanced Large-scale Analysis of MAgnetic confinement for Numerical Design, Engineering and Research (SALAMANDER). SALAMANDER is an open source tool developed at Idaho National Laboratory for conducting multiphysics and multscale analysis of Plasma Facing Components (PFCs). For the monoblock, solid heat conduction is modeled in the block with a constant heat flux of 1E7 W/m^2 applied to the top of the block to represent peak heat flux from the plasma. The cooling channel (which is pressurized water at a Reynolds number of ~1,000,000) is modeled using a k-epsilon turbulence model with standard wall functions. The heat transfer between the solid and fluid domain is modeled using the Dittus-Boelter correlation for the convective heat transfer coefficient. This work aims to show the importance of a multiphysics modeling approach when performing simulations on PFCs.

70 - PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulation of flame acceleration and deflagration-to-detonation transition with heat transfer in HE product and fractured HE

A simple heat transfer model has been developed for high explosives. For the gaseous high explosive (HE) product, the Smargorinsky eddy viscosity model has been adopted for estimating the turbulent thermal conductivity. Convective heat transfer in the fractured solid HE is modeled as the effective thermal conduction in the porous media, which is equivalent to the enthalpy flux carried by the velocity determined by Darcy’s law with Forchheimer term.Other features of the model include a multi-reaction chain, equation-of-state for each species involved, and a two-phase hotspot model. Test calculations are performed with various initial disturbances. Onset and progress of flame acceleration (FA) and chemical reactions are strongly affected by the heat transfer and initial disturbance. Staged development of FA also affirms the necessity of modeling based on a multi-reaction chain. When conditions warrant,e.g., stronger initial disturbance, FA evolves to deflagration-to-detonation transition (DDT),demonstrating automatic capture of the DDT.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗
DOE OSTI2024

Enabling the broader adoption of fusion simulation on complex geometry

This project addressed a key barrier to advanced fusion and nuclear simulation: the difficulty of performing high-fidelity Monte Carlo neutronics directly on complex, real-world CAD geometry. Traditional workflows require engineers to rebuild CAD models as simplified constructive solid geometry, a time-consuming and error-prone process that limits design iteration and broader adoption of simulation tools. The goal of this Phase I SBIR was to make CAD-based neutronics practical, accessible, and robust for industrial and research users. During the project, Coreform significantly enhanced the Direct Accelerated Geometry Monte Carlo (DAGMC) workflow and fully integrated it into Coreform Cubit as a first-class capability. Major achievements include optimized material assignment and surface meshing workflows, substantial performance improvements to geometry imprinting and preparation, native export of DAGMC models, and new visualization tools to support OpenMC source definition and lost-particle debugging. Coreform also expanded Cubit’s capabilities as a full OpenMC preprocessor, including the ability to convert OpenMC constructive solid geometry models back into CAD for visualization, multiphysics coupling, and debugging. In collaboration with Argonne National Laboratory, the project delivered comprehensive new DAGMC documentation and training materials, transforming DAGMC from a research-oriented tool into a production-ready workflow. Results were disseminated through tutorials, conference training, and multiple well-attended webinars demonstrating integrated CAD-based neutronics and multiphysics workflows. Overall, this project demonstrated that high-fidelity Monte Carlo simulations can be performed directly on complex CAD geometry, reducing setup time, improving usability, and enabling faster, more informed design decisions for fusion and nuclear energy systems.

42 ENGINEERING↗
DOE OSTI2026

Native Architecture of Wheat Straw Cell Walls: A Unified Model from X-ray Scattering and Solid-State NMR

Plant secondary cell walls constitute the dominant reservoir of renewable biomass, comprising tightly packed cellulose, hemicellulose, and lignin at the nanoscale. Recent advances in solid-state NMR spectroscopy and the availability of small-angle X-ray scattering for biomass characterization have led to an accumulation of experimental data on cell wall organization, yet no explicit structure model has simultaneously satisfied both Xray and NMR observations. Using wheat straw as a model system, we propose a structural framework consistent with current knowledge of cellulose biosynthesis, X-ray scattering data, and one- and two-dimensional 13 C solid-state NMR spectra. In this model, 18-chain elementary fibrils align in parallel and populate the cross-section at random. Arabinose-substituted xylan shows no conformational dependence for cellulose-binding in wheat, and only a minor fraction of 2-fold xylan appears in close proximity to cellulose, unlike in Arabidopsis, where xylan is more tightly attached to the cellulose surface. While NMR data cannot unambiguously resolve the internal arrangement of the 18 glucan chains, X-ray scattering profiles uniquely constrain the fibril size and exclude the possibility of tight bundling in the intact walls. The specific interaction between the matrix polymers and the cellulose elementary fibrils must be reconsidered in light of the small interfibril spaces, which bring the matrix components into spatial proximity with cellulose even in the absence of attractive interactions. These findings provide fundamental molecular-level insight into cellulose fibril architecture and matrix−polymer interactions, resolving longstanding discrepancies between spectroscopic and scattering data and advancing our understanding of biopolymer assembly into structurally and functionally versatile lignocellulosic biomaterials.

Carbohydrates↗
DOE OSTI2022

Correction to: Investigating the rheology of fluidized and non-fluidized gas-particle beds: implications for the dynamics of geophysical flows and substrate entrainment

We report natural geophysical mass flows are among the most complex granular systems and their dynamics are often modified by the presence of an interstitial fluid. Prediction of their runout requires the development of models estimating the solid stresses in these hazardous currents wherein excess pore-fluid pressure can develop. We use discrete element modelling (DEM-CFD) with a Coarse-Graining post-processing technique (CG) to investigate the rheology of unsteady gas-particle fluidized to non-fluidized granular beds placed on horizontal and inclined planes. Similar to fluidized beds immersed in viscous fluids, the effective friction coefficient of air-fluidized beds can be defined as a function of the classic μ(I)-rheology and the non-dimensional fluid or solid pressure to explain the failure and dynamics of granular flows with excess pore pressure on inclines. However, dilation imposed by fluid drag and particle collisions in gas-particle fluidized beds can drastically change its effective frictional properties. In contrast with the common assumption in water-particle flows that granular temperature is negligible, in our gas-particle simulations, the contribution of the velocity fluctuations to the stress tensor is significant. Hence, the shear stress is found to be non-zero even when the flow is fully fluidized in the inertial regime. These results suggest the need to better understand velocity fluctuations to predict the effective viscosity of sheared fluidized mixtures and are relevant for many applications. Notably, a unified approach is useful for many geophysical flows that encompass a range of fluidization conditions in a single flow such as pyroclastic density currents and snow avalanches.

58 GEOSCIENCES↗
DOE OSTI2021

Deposition control of model glasses with surface-mediated orientational order

We introduce a minimal model of solid-forming anisotropic molecules that displays, in thermal equilibrium, surface orientational order without bulk orientational order. The model reproduces the nonequilibrium behavior of recent experiments in which a bulk nonequilibrium structure grown by deposition contains regions of orientational order characteristic of the surface equilibrium. This order is deposited, in general, in a nonuniform way because of the emergence of a growth-poisoning mechanism that causes equilibrated surfaces to grow slower than non-equilibrated surfaces. We use evolutionary methods to design oscillatory protocols able to grow nonequilibrium structures with uniform order, demonstrating the potential of protocol design for the fabrication of this class of materials.

Whitelam, Stephen (ORCID:0000000200866803)↗
DOE OSTI2025

Phase-field model of alloy solidification far from chemical equilibrium at the solid-liquid interface

We further develop a recently introduced phase-field model of far-from-equilibrium alloy solidification under additive manufacturing conditions [K. Ji et al., Phys. Rev. Lett. 130, 026203 (2023)]. This model utilizes enhanced solute diffusivity within the spatially diffuse interface region to quantitatively capture solute trapping with a larger interface width, thereby making simulations on experimentally relevant length and timescales computationally feasible. The main developments presented here include testing the robustness of different variational formulations, extending the model to concentrated alloys by incorporating solid and liquid free energies from thermodynamic databases, as illustrated for hypoeutectic Al-Ag alloys with CALPHAD, extending convergence tests as a function of interface width to 3D, and carrying out simulations in both 2D and 3D to examine existing theories of microstructure development. Our results indicate that the simplest variational formulation that interpolates the bulk free-energy density between its solid and liquid forms is the most robust. Remarkably, for hypoeutectic Al-Ag alloys, this formulation yields a high-velocity nonequilibrium phase diagram that is independent of interface width, thereby demonstrating that the framework of enhanced solute diffusivity can be nontrivially extended to concentrated alloys. Other variational formulations have restricted ranges of materials or processing parameters that can be reliably modeled. We use 2D simulations to construct high-velocity microstructure selection maps for dilute Al-Cu alloys. The results validate the important role of latent heat rejection at the interface and extend the limited predictions of linear stability analysis [A. Karma and A. Sarkissian, Phys. Rev. E 47, 513 (1993)] and sharp-interface 1D simulations to fully nonlinear regimes. Furthermore, 3D simulations, carried out using a computationally tractable axisymmetric cellular/dendritic interface shape, demonstrate a good convergence similar to that observed in 2D as a function of interface width. Full 3D simulations, in turn, reveal that the standard theory of absolute stability is a good predictor of the upper critical velocity beyond which steady-state growth becomes unstable, despite the different morphological manifestations of this instability in 2D and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Observing and Modeling the Sequential Pairwise Reactions that Drive Solid‐State Ceramic Synthesis

Abstract Solid‐state synthesis from powder precursors is the primary processing route to advanced multicomponent ceramic materials. Designing reaction conditions and precursors for ceramic synthesis can be a laborious, trial‐and‐error process, as heterogeneous mixtures of precursors often evolve through a complicated series of reaction intermediates. Here, ab initio thermodynamics is used to model which pair of precursors has the most reactive interface, enabling the understanding and anticipation of which non‐equilibrium intermediates form in the early stages of a solid‐state reaction. In situ X‐ray diffraction and in situ electron microscopy are then used to observe how these initial intermediates influence phase evolution in the synthesis of the classic high‐temperature superconductor YBa 2 Cu 3 O 6+ x (YBCO). The model developed herein rationalizes how the replacement of the traditional BaCO 3 precursor with BaO 2 redirects phase evolution through a low‐temperature eutectic melt, facilitating the formation of YBCO in 30 min instead of 12+ h. Precursor selection plays an important role in tuning the thermodynamics of interfacial reactions and emerges as an important design parameter in planning kinetically favorable synthesis pathways to complex ceramic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data-driven nonlocal model for fragmentation in the crushing of solids

A technique is proposed for reproducing particle size distributions in three-dimensional simulations of the crushing and comminution of solid materials. The method is designed to produce realistic distributions over a wide range of loading conditions, especially for small fragments. In contrast to most existing methods, the new model does not explicitly treat the small-scale process of fracture. Instead, it uses measured fragment distributions from laboratory tests as the basic material property that is incorporated into the algorithm, providing a data-driven approach. The algorithm is implemented within a nonlocal peridynamic solver, which simulates the underlying continuum mechanics and contact interactions between fragments after they are formed. Finally, the technique is illustrated in reproducing fragmentation data from drop weight testing on sandstone samples.

58 GEOSCIENCES↗