Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational structural mechanics”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Establishing model credibility for process-microstructure-property relationships in additive manufacturing using exascale computing

Additive Manufacturing (AM) of alloys holds significant promise as a disruptive technology in various industries, yet its adoption is often hindered by challenges in achieving consistent part quality. These issues are primarily due to the complex process-microstructure-property (PSP) relationships inherent to AM. Computational models can greatly aid in understanding these relationships, but their widespread impact and adoption has been limited by a lack of validated, open-source, and computationally efficient PSP modeling frameworks and hardware limitations. Here, this study leverages the ExaAM software suite and data from the AMBench-2018 series of laser powder bed fusion (LPBF) benchmark experiments to perform a comprehensive model assessment, including verification, validation, sensitivity analysis, and uncertainty quantification. The RADICAL-EnTK workflow manager was used to perform an ensemble of heat transport, solidification, and mechanical response simulations on the exascale computer Frontier, considering uncertainties in critical model inputs such as laser spot size and nucleation parameters, and consisting of 125 explicit grain structure simulations and 7875 crystal plasticity simulations. For a selected location within the Inconel 625 AMBench-2018 test artifact, sensitivity analysis and uncertainty quantification were performed using the predicted distributions of grain structure and mechanical properties. Qualitative agreement was found between the predicted grain size and texture and the observed AMBench-2018 microstructure, the mean predicted yield stress was within 5% of the experimental measurement mean, and the mean predicted engineering stress at 5% strain was within 10% of the experimental measurement mean. The insights gained from development and validation of the ExaAM PSP modeling framework will help guide future directions for enhancing the credibility and reliability of PSP models in AM, thereby accelerating the adoption of AM technologies in various industries.

Additive manufacturing↗
DOE OSTI2025

Impact of processing humidity on ionomer film structure and performance in hydroxide exchange membrane electrolyzers

Hydroxide exchange membrane electrolyzers (HEMELs) enable hydrogen production using low-cost, earth-abundant materials. Improving electrode fabrication is integral to enhancing device performance, and ionomer-responsible for transporting hydroxide and mechanically supporting the catalyst-is a major component. Here, we use experiments and computation to study the effects of relative humidity (RH) during the drying process of poly(aryl piperidinium) ionomer films on HEMEL electrodes. Broadly, the drying environments determine the physical structure and electrochemical traits of the ionomer network. High RH drying yields a highly porous network with excessive water uptake, structural defects, washout, and 64% reduction in hydroxide conductivity. Extremely low RH drying produces an overly compact pore network that hinders hydroxide mobility. In contrast, moderately low RH drying (9% RH) creates an ionomer film with well-balanced traits: excellent mechanical stability and connectivity needed for catalyst retention and hydroxide transport, which improves HEMEL performance by 40% at 1.8 V compared to suboptimal RHs. This research advances HEMEL manufacturing by providing a simple, scalable, and low-cost approach to optimize electrode ionomer films.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Extending PETSc's Composable Hierarchical Solvers (Final Technical Report)

This report documents research activities conducted at CU Boulder as part of Extending PETSc’s Composable Hierarchical Solvers, which has been part of a collaboration with Argonne National Laboratory (separate award). Our work has focused on performance-portable end-to-end GPU solvers demonstrated via exemplary applications in nonlinear fluid and structural mechanics. We describe advances in algorithmic composition and analysis in the context of these applications, but the implementations are fully documented and decoupled, and in use by other projects. We believe the vertical integration achieved through collaboration with ECP’s CEED and the PSAAP center at CU was necessary to take risks with data structures and algorithms.

42 ENGINEERING↗
DOE OSTI2024

Models for Single–Site Heterogeneous Catalysts on Carbon: MoO 2 Epoxidation Catalyst Anchored to a Fullerene

Single-site molybdenum dioxo catalysts, fullerenol/MoO 2 , are prepared via grafting precursor (DME)MoO 2 Cl 2 onto a highly polyhydroxylated fullerene (ful) and an isomerically-pure and well-defined fullerene (ful*). These catalyst structures are characterized by ICP-OES, XPS, XANES, EXAFS, DRIFT, Raman, and NMR spectroscopy, and DFT. Mo 3d 5/2 XPS and Mo K-edge XANES assign the oxidation state as Mo(VI). Mo EXAFS data fitting reveals two Mo=O double and two Mo–O single bonds at distances of 1.7 and 1.9 Å, respectively, while an Mo=O stretchingl mode is observed at ~950 cm –1 by DRIFT and Raman spectroscopy. These data align well with DFT computational results, supporting the proposed catalyst structure as Fullerene(-μ-O-) 2 M(=O) 2 . Additionally, DFT provides insight into the energetically favorable grafting sites for an isomerically pure fullerenol. The scope of fullerenol/MoO 2 mediated alkene epoxidation includes abiotic alkenes, natural occurring terpenes, and conjugated olefins. For cyclooctene the rate law is first-order in [Mo], near first order in [olefin] and zero-order in [t-butyl hydroperoxide]. A plausible reaction mechanism involves peroxide addition first and then cyclooctene addition directly across the peroxo bond forming the epoxide product, consistent with DFT computation. Overall, fullerenol/MoO 2 shows promise as a sustainable and structurally well-defined system with versatile catalytic activity and good epoxidation recyclability.

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

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes

Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH 2 -capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH 3 (CH 2 ) 5 CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young’s modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. Here, we construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a “Goldilocks” tube radius at which capillary wave fluctuations in the tube wall are minimized.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Mechanical behaviour of additively manufactured metals

Additive manufacturing is reshaping the production of engineering components in diverse industries, such as the automotive, aerospace, defense, and biomedical sectors, by offering unprecedented design flexibility. The non-equilibrium processing conditions of additive manufacturing generate materials with unique microstructures and tailored mechanical properties that are often unattainable through conventional routes. This review focuses on recent advances in additively manufactured metals that demonstrate distinctive mechanical behaviors, including strength-ductility synergy, microstresses and gradient plasticity, fracture and fatigue resistance, and high-temperature creep performance. Here, we examine the mechanisms and micromechanical effects arising from the heterogeneous microstructures fabricated by additive manufacturing, to guide the design of a wide range of high-performance structural materials. Furthermore, we discuss critical research needs and emerging opportunities in process control, alloy design, advanced characterization, high-fidelity computational modeling, and machine learning aimed at achieving exceptional mechanical properties in additively manufactured metals.

Additive Manufacturing↗
DOE OSTI2026

Gas-driven short disconnection mitigates thermal runaway in Li-ion batteries under mechanical abuse

Mechanical abuse poses a critical safety risk to Li-ion batteries by inducing internal short circuits that initiate thermal runaway. Remarkably, voltage recovery frequently emerges during thermal runaway initiation, a phenomenon that conventional theories fail to explain. Here, our study develops a new multiphysics mechanism to explain voltage recovery, termed gas-driven short disconnection, whereby internal gas pressure mechanically disengages short-circuit contacts and causes the cell voltage to rebound. This mechanism incorporates gas generation and its structural impact on the short circuit. Real-time optical and thermal imaging and X-ray computed tomography reveal fluid-structure interaction between internal gas flow and adjacent shorting contacts. We establish a mechanistic framework linking gas-driven short disconnection to cell-level voltage and temperature responses, elucidating the extension–truncation pattern of voltage recovery. Furthermore, thermal regime maps show that a voltage recovery duration exceeding 5 s correlates with limited temperature rise below 150 °C, indicating that sustained short-circuit disconnection suppresses Joule heating. Additionally, a dimensionless criterion is deduced from scaling analysis for physical plausibility of gas-driven short disconnection in mechanically abused cells. This finding inspires smart venting control, which regulates gas release to maintain the internal pressure while dissipating gas enthalpy, thereby providing a device-level strategy for thermal runaway mitigation.

Gas-driven short disconnection↗
DOE OSTI2020

The mechanism of carboxylative cyclization of propargylamine by N-heterocyclic carbene complexes of Au(I)

Here, the complex [Au(IPr)(2-oxazolidinone)] (1 = IKa) was prepared from reaction of [Au(IPr)Cl] (2), K 2 CO 3 , and propargyl amine (PPA). Kinetic studies have been performed for acid cleavage of 1 to yield the oxazolidinone product and an [Au(IPr)(X)] adduct. The fastest rates of cleavage were found to occur for the hydrogen chloride salt of PPA (PPA-HCl) and for the CO 2 adduct of PPA, PPA-CO 2 = the carbamic acid (CA). This transformation was studied as a function of [CA], pressure of CO 2 as well as temperature. Detailed computational studies support the formation of a key intermediate and are also in agreement with a rapid carbonylation/decarbonylation reaction. The computed reactions mechanisms for addition of PPA-HCl and CA are also presented as well as the crystal structure of 1.

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

Modulating CO 2 Electrocatalytic Conversion to the Organics Pathway by the Catalytic Site Dimension

Electrochemical reduction of carbon dioxide to organic chemicals provides a value-added route for mitigating greenhouse gas emissions. Here, we report a family of carbon-supported Sn electrocatalysts with the tin size varying from single atom, ultrasmall clusters to nanocrystallites. High single-product Faradaic efficiency (FE) and low onset potential of CO 2 conversion to acetate (FE = 90% @ –0.6 V), ethanol (FE = 92% @ –0.4 V), and formate (FE = 91% @ –0.6 V) were achieved over the catalysts of different active site dimensions. The CO 2 conversion mechanism behind these highly selective, size-modulated p-block element catalysts was elucidated by structural characterization and computational modeling, together with kinetic isotope effect investigation.

10 SYNTHETIC FUELS↗
DOE OSTI2026

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

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

First-principles study of an S = 1 quasi one-dimensional quantum molecular magnetic material

We use density functional theory to study the structural, magnetic, and electronic structures of the organometallic quantum magnet NiCl 2-4 SC (NH 2 ) 2 (DTN). Recent work has demonstrated the quasi one-dimensional nature of the molecular crystal and studied its quantum phase transitions at low temperatures. The system includes a magnetoelectric (ME) coupling and, when doped with Br, the presence of an exotic Bose-glass state. Using the generalized gradient approximation with inclusion of a van der Waals term to account for weak intermolecular forces and by introducing a Hubbard U term to the total energy, we systematically show that our calculations reproduce the magnetic anisotropy, the intermolecular exchange coupling strength, and the magnetoelectric effect in DTN, which have been observed in previous experiments. Further analysis of the electronic structure gives insight into the underlying magnetic interactions, including what mechanisms may be causing the ME effect. Using this computationally efficient model, we predict what effect applying an electric field might have on the magnetic properties of this quantum magnet.

1-dimensional spin chains↗
DOE OSTI2022

Aircraft Icing Severity Evaluation

Aircraft icing refers to the ice buildup on the surface of an aircraft flying in icing conditions. The ice accretion on the aircraft alters the original aerodynamic configuration and degrades the aerodynamic performances and may lead to unsafe flight conditions. Evaluating the flow structure, icing mechanism and consequences is of great importance to the development of an anti/deicing technique. Studies have shown computational fluid dynamics (CFD) and machine learning (ML) to be effective in predicting the ice shape and icing severity under different flight conditions. CFD solves a set of partial differential equations to obtain the air flow fields, water droplets trajectories and ice shape. ML is a branch of artificial intelligence and, based on the data, the self-improved computer algorithms can be effective in finding the nonlinear mapping relationship between the input flight conditions and the output aircraft icing severity features.

Li, Sibo (ORCID:000000021705844X)↗
DOE OSTI2020

Different bonding type along each crystallographic axis: Computational study of poly($\textit{p}$-phenylene terephthalamide)

Poly($\textit{p}$-phenylene terephthalamide) (PPTA) exhibits van der Waals (vdW) bonding along its $\textit{a}$ axis, hydrogen bonding along its $\textit{b}$ axis, and covalent bonding along its $\textit{c}$ axis. In this work, we explore the structural and mechanical properties of PPTA using density functional theory with various functionals including LDA, PBE, PBE+rVV10L, SCAN, and SCAN+rVV10, compared with available experiments. The hierarchy of nonempirical semilocal functionals LDA, PBE, and SCAN (not fitted to any multicenter bonded system) includes differing amounts of intermediate-range vdW interaction. rVV10 is the long-range vdW correction. (rVV10L differs from rVV10 only in the value of a range parameter.) Among the tested functionals, SCAN shows the best performance for the lattice parameters of PPTA along the two crystal directions involving vdW or hydrogen-bond interaction. The equilibrium lattice constants obtained by SCAN and PBE+rVV10L are closest to experimental data, while SCAN+rVV10 slightly overbinds the system. We study the mechanical response of PPTA by applying strain along three lattice directions. Due to the inclusion of vdW interaction, SCAN, PBE+rVV10L, and SCAN+rVV10 all exhibit correct bonding strain-energy curves. On the contrary, PBE strongly underestimates the vdW interaction needed to resist uniaxial stretching along the $\textit{a}$ axis. The Young's modulus and yield strength of PPTA are computed and compared with previous results. The experimental values are much smaller than the computed ones, mainly due to the fact that the PPTA fiber samples used for measurements are mechanically weaker than the perfect molecular crystals considered in the simulations. Interestingly, when a compressive uniaxial stress of 25 Kbar is applied along the $\textit{b}$ axis, a structural phase transition, in which the hydrogen bonds reform along one diagonal of the $\textit{ab}$ rectangle, is predicted by SCAN.

36 MATERIALS SCIENCE↗
DOE OSTI2022

4D evolution of Cr 23 C 6 precipitates in neutron-irradiated and annealed HT-UPS steel observed via synchrotron micro-computed tomography

High-temperature-ultrafine precipitate strengthened (HT-UPS) steel is a potential structural material for advanced nuclear reactors; however, its irradiation response is not well understood. This research provides insight into irradiation-induced effects, such as precipitate evolution mechanisms and four-dimensional morphological evolution, in HT-UPS steel using synchrotron micro-computed tomography. Identical specimens were characterized pre-irradiation and post-irradiation following neutron exposure up to 0.3 displacements per atom at 600 °C. Irradiation effects were also differentiated from the annealing response of precipitates. Following neutron irradiation, the average Cr 23 C 6 precipitate size reduced, affected by the synergy of nucleation and growth, ballistic dissolution, and inverse coarsening, which was observed at fluences an order of magnitude lower than previously observed. Annealing at 600 °C for 32 h increased the average Cr 23 C 6 precipitate size and decreased the phase fraction, attributed to precipitate coarsening. The precipitate morphology evolution and resultant mechanisms can be utilized to parameterize and validate microstructural models simulating radiation damage or annealing.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Computational Modeling of Molten Salt Infiltration and Oxidation in Nuclear Graphite

Graphite is utilized as a moderator and reflector in advanced nuclear reactor designs due to its high thermal conductivity, neutron moderation properties, and resistance to radiation damage. However, its longterm performance and reliability are challenged by degradation mechanisms such as molten salt infiltration in molten salt reactors (MSRs) and oxidation in gas-cooled reactors (GCRs). These mechanisms can compromise the structural integrity and operational lifetime of graphite components, necessitating a more detailed assessment of their physical behavior. This report focuses on the development of computational models for molten salt infiltration and oxidation of graphite to aid the design and performance analysis of graphite components. For molten salt infiltration, a computational framework is developed that couples incompressible Navier-Stokes and phase-field model to simulate the penetration of molten salt into graphite?s interconnected pore structure. Initial model verification is performed using two-phase flows in two dimensions, demonstrating the models ability to capture fundamental physical behavior and agree with analytical solution. This framework is then applied to a realistic IG110 nuclear graphite , where a computed tomography extracted pore geometry is used to analyse the infiltration behavior of FLiNaK molten salt. This model provides insights into how the microstructure and other relevant parameters influence the transport pathways of molten salt into graphite, potentially offering a means to rapidly evaluate a graphite grade?s resistance to infiltration. For oxidation, the report details pore-scale mass and heat transport models, describing the diffusion of gases, reaction kinetics, and thermal effects. Additionally, this report highlights inconsistencies in the existing volume-averaged macroscopic model, particularly in upscaling of reaction kinetics and flux terms, and surface to volume transformations. These inconsistencies suggest that current formulations may not accurately capture the experimentally observed graphite oxidation process, highlighting the need for improved model development. This work advances the development of physics-based computational models for graphite degradation, contributing to improved predictive models for next-generation nuclear reactor designs. Future efforts will focus on refining the infiltration model to address non-physical behaviors and enhance its robustness. Additionally, for oxidation, further studies will employ the principles of volume averaging to rigorously derive the upscaled equations, potentially in collaboration with subject matter experts.

Computational Modeling of Molten Salt Infiltration↗
DOE OSTI2024

Machine-Learned Linear Structural Dynamics

The tension between accuracy and computational cost is a common thread throughout computational simulation. One such example arises in the modeling of mechanical joints. Joints are typically confined to a physically small domain and yet are computationally expensive to model with a high-resolution finite element representation. A common approach is to substitute reduced-order models that can capture important aspects of the joint response and enable the use of more computationally efficient techniques overall. Unfortunately, such reduced-order models are often difficult to use, error prone, and have a narrow range of application. In contrast, we propose a new type of reduced-order model, leveraging machine learning, that would be both user-friendly and extensible to a wide range of applications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Fused Deposition Modeling Additive Manufacturing of Carbonized Structures via Waste-Enhanced Filaments

This report details the development, characterization, and use of coal-enhanced composite materials in additive manufacturing applications. High coal loading formulations—containing up to 70 wt.% coal—were successfully extruded and processed using commercially available 3D printers. Extensive experimental testing was conducted to assess the mechanical, thermal, and microstructural properties of the composites. In parallel, multi-scale computational modeling was employed to elucidate atomistic interactions and evaluate the effects of printing-induced defects on structural performance. Large-scale printability trials demonstrated the feasibility of fabricating complex components for tooling and construction applications, including wind turbine blade molds and modular wall sections. Techno-economic analyses demonstrated the cost-effectiveness and scalability of coal-enhanced composites for large-scale additive manufacturing applications such as wind turbine blade tooling.

01 COAL, LIGNITE, AND PEAT↗