Search NASASearch

SEARCH · Search NASA

Results for “Computational solid 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.

1,011 records · Page 22

Evaluating Fracture Behavior of Bioinspired Alumina-YSZ Composites through Static and Dynamic Mechanical Testing

Ceramic materials are known for their high hardness and strength but are limited by low toughness and sudden failure. Inspired by the microstructure of dental enamel, which features undulating rods that promote crack deflection and energy absorption, architected specimens were designed and manufactured. A bespoke ball-on-ring (BoR) testing apparatus measured the static biaxial rupture strength of these bioinspired materials. Additionally, impact testing with spherical steel projectiles assessed their behavior under dynamic conditions. Monolithic alumina specimens were compared to architected alumina reinforced with yttria-stabilized zirconia (YSZ) rods produced via direct ink write 3D printing. BoR tests revealed that monolithic alumina had a Weibull modulus of 10.53 and a characteristic strength of 372.3 MPa, while the composite specimens showed a Weibull modulus of 5.46 and a characteristic strength of 213.3 MPa. Although the composite failed at lower stresses and projectile velocities, it exhibited notable crack deflection and fracture resistance, with cracks being effectively interrupted by the rod inclusions. The composite specimens also resulted in fewer and larger fragments upon impact. Finite element simulations confirmed the effectiveness of the rod structures in enhancing fracture resistance.

Bioinspired materials

Boundary Layer Analysis of Shock Tube Flows

Shock tubes offer a controlled environment to reproduce kinetic and radiative phenomena characteristic of atmospheric entry flows under ground-test conditions. The boundary layer developing behind the incident shock wave determines the available test time, influences particle residence times important for similarity scaling, and can affect radiative energy transport. In this study, we couple a quasi-1D space marcher with the compressible boundary-layer equations to numerically compute the post-shock flow in a shock frame of reference for various test gas mixtures representative of different planetary atmospheres. The solvers are individually verified against CFD simulations and analytical correlations available in the literature. Coupled solutions are computed for a finite-rate chemistry in the boundary layer and a non-catalytic isothermal wall. Results yield refined estimates of the maximum separation distance, as well as insights into concentration profiles of relevant species within the boundary layer.

Andrea Fagnani

Time-Resolved Neutron Imaging for Hydrogen Uptake in Subsurface Lithologies

Geologic hydrogen production and underground storage are increasingly important for meeting rising energy demands while providing clean-combustion advantages. However, hydrogen’s high diffusivity and propensity for leakage through porous media necessitate direct evaluation of its transport behavior in subsurface materials. Whereas X-ray microcomputed tomography (μCT) studies often employ contrast agents or surrogate gases, this study leverages neutron transmission radiography/CT to observe hydrogen migration in situ. This work represents the first demonstration of real-time neutron radiography of hydrogen migration in reservoir and caprock lithologies. Cylindrical cores of Indiana limestone, Amherst Gray sandstone, and Tumey shale were subjected to constant-pressure hydrogen charging and scanned in real time using high-resolution neutron radiography. Results indicate immediate hydrogen infiltration in sandstone and limestone, with homogeneous distribution detected throughout their pore structure. In contrast, hydrogen remained largely absent from fine-grained shale under the same pressure, except in an apparently localized fracture zone, where neutron signatures confirmed the presence of hydrogen. Subsequent neutron CT of the sandstone sample, using image subtraction against an uncharged reference, corroborated hydrogen distribution patterns. Even under lowpressure, single-phase conditions, distinct neutron imaging signatures of hydrogen were achieved. These preliminary findings underscore the potential of neutron imaging for advancing subsurface hydrogen migration research.

08 HYDROGEN

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation

Galaxy Zoo DESI: large-scale bars as a secular mechanism for triggering AGNs

ABSTRACT Despite the evidence that supermassive black holes (SMBHs) co-evolve with their host galaxy, and that most of the growth of these SMBHs occurs via merger-free processes, the underlying mechanisms which drive this secular co-evolution are poorly understood. We investigate the role that both strong and weak large-scale galactic bars play in mediating this relationship. Using 48 871 disc galaxies in a volume-limited sample from Galaxy Zoo DESI, we analyse the active galactic nucleus (AGN) fraction in strongly barred, weakly barred, and unbarred galaxies up to $z = 0.1$ over a range of stellar masses and colours. After controlling for stellar mass and colour, we find that the optically selected AGN fraction is $31.6 \pm 0.9$ per cent in strongly barred galaxies, $23.3 \pm 0.8$ per cent in weakly barred galaxies, and $14.2 \pm 0.6$ per cent in unbarred disc galaxies. These are highly statistically robust results, strengthening the tantalizing results in earlier works. Strongly barred galaxies have a higher fraction of AGNs than weakly barred galaxies, which in turn have a higher fraction than unbarred galaxies. Thus, while bars are not required in order to grow an SMBH in a disc galaxy, large-scale galactic bars appear to facilitate AGN fuelling, and the presence of a strong bar makes a disc galaxy more than twice as likely to host an AGN than an unbarred galaxy at all galaxy stellar masses and colours.

Astronomy & Astrophysics

Egg yolk as a model for gelation: From rheometry to flow physics

Egg yolks are an excellent model for studying sol-gel transitions, particularly the power law viscoelasticity that defines the critical point of gelation. However, prior studies lack comprehensive datasets and fail to visualize flow behavior linked to temperature and time-dependent linear and nonlinear rheology. Here, we present a detailed dataset characterizing egg yolk viscoelasticity across temperature, time, and forcing amplitude using oscillatory shear, step strain, step stress, and constant high strain rate. Novel protorheology visualizations link rheological properties with observable flow behavior. Our findings highlight the nuanced determination of the critical gel point, emphasizing observation timescale dependencies. We compare methods to identify critical temperatures for gelation, including power law viscoelasticity, moduli crossover, diverging zero-shear viscosity, and emerging equilibrium elastic modulus, while visualizing flow consequences near these transitions. Egg yolk is an accessible, realistic, and nontoxic material relevant to the physicist and the chef alike, making it ideal for understanding the rheology of critical gels. By integrating protorheology photos and videos with rigorous rheometric data, we deepen the understanding of critical gels, with broader impacts for studying other materials with sol-gel transitions.

Marsh, Maxwell C. [Department of Mechanical Scienc

Chevron Natural Gas Carbon Capture Technology Testing Project (Final Technical Report)

The objective of this project was to design, construct, commission, and operate an engineering scale post combustion carbon capture system to validate the technical maturity, operability, and scalability of Svante’s VeloxoTherm™ solid sorbent carbon capture technology under representative flue gas conditions. The project was executed at Chevron’s Kern River oil field in the San Joaquin Valley of California and evaluated carbon capture performance using slipstream flue gas from natural gas fired once through steam generators (OTSG), including indicative coal and natural gas combined cycle (NGCC) operating scenarios.

03 NATURAL GAS

All-Polyester Multilayer Plastics (‘All-Polyester MLPs’): A Redesign for Inherently Recyclable Plastics

Multilayer plastics (MLPs) are widely used in modern packaging because they combine multiple functions such as oxygen and moisture barriers, mechanical strength, puncture resistance, and heat sealability into lightweight and cost-effective packaging solutions. These attributes are essential for food, beverage, pharmaceutical, and consumer goods packaging. However, conventional commercial MLPs typically consist of five to twelve layers made from chemically incompatible materials, including polyolefins, polyethylene terephthalate, nylon, ethylene-vinyl alcohol (EVOH), adhesives, and tie layers. This complexity makes MLPs extremely difficult to recycle. As a result, the vast majority of multilayer plastics (MLPs) are disposed of through landfilling or incineration. This disposal pathway perpetuates demand for virgin material production, thereby driving high industrial energy consumption, increasing greenhouse gas (GHG) emissions, and contributing to the long-term accumulation of plastic waste. The scope of work included: (1) Design and fabrication of all-polyester multilayer structures using commercially relevant processing methods; (2) Experimental validation of barrier, mechanical, and sealing performance; (3) Demonstration of both mechanical and chemical recycling pathways; and (4) Comprehensive techno-economic analysis (TEA) and life-cycle assessment (LCA) to quantify cost, energy, and environmental impacts.

36 MATERIALS SCIENCE

QED corrections to the thermal neutrino interaction rate

Abstract Motivated by precision computations of neutrino decoupling at MeV temperatures, we show how QED corrections to the thermal neutrino interaction rate can be related to the electron-positron spectral function as well as an effective$$\overline{\nu }\nu \gamma $$vertex. The spectral function is needed both in a timelike and in a spacelike domain, and for both of its physical polarization states (transverse and longitudinal with respect to spatial momentum). Incorporating an NLO evaluation of this spectral function, an estimate of the$$\overline{\nu }\nu \gamma $$vertex, and HTL resummation of scatterings mediated by soft Bose-enhancedt-channel photons, we compute the interaction rate as a function of the neutrino momentum and flavour. Effects on the – (0. . .2)% level are found, noticeably smaller than a previous estimate of a related quantity.

Physics

32 examples of LLM applications in materials science and chemistry: towards automation, assistants, agents, and accelerated scientific discovery

Abstract Large language models (LLMs) are reshaping many aspects of materials science and chemistry research, enabling advances in molecular property prediction, materials design, scientific automation, knowledge extraction, and more. Recent developments demonstrate that the latest class of models are able to integrate structured and unstructured data, assist in hypothesis generation, and streamline research workflows. To explore the frontier of LLM capabilities across the research lifecycle, we review applications of LLMs through 32 total projects developed during the second annual LLM hackathon for applications in materials science and chemistry, a global hybrid event. These projects spanned seven key research areas: (1) molecular and material property prediction, (2) molecular and material design, (3) automation and novel interfaces, (4) scientific communication and education, (5) research data management and automation, (6) hypothesis generation and evaluation, and (7) knowledge extraction and reasoning from the scientific literature. Collectively, these applications illustrate how LLMs serve as versatile predictive models, platforms for rapid prototyping of domain-specific tools, and much more. In particular, improvements in both open source and proprietary LLM performance through the addition of reasoning, additional training data, and new techniques have expanded effectiveness, particularly in low-data environments and interdisciplinary research. As LLMs continue to improve, their integration into scientific workflows presents both new opportunities and new challenges, requiring ongoing exploration, continued refinement, and further research to address reliability, interpretability, and reproducibility.

Computer Science

Advances in Tactical Rocket Propulsion

An AGARD Colloquium on 'Progress in Tactical Rocket Propulsion' was held at La Jolla, California, on April 22 and 23, 1965. The unclassified papers, or unclassified versions of classified papers, presented at this meeting form the contents of the present AGARD conference proceedings, which has been identified by a new title that reflects more properly the scope of the remaining technical papers. The articles contained in this volume represent significant evaluations of the status of solid-propellant, liquid-propellant and hybrid-propellant rocket development. Together with the reviews, they should provide a useful reference source for research workers and students alike.

S S Penner

Uncertainty quantification for competing failure mechanisms in unidirectionally reinforced carbon–carbon composites

Microstructure-informed finite element models play a key role in the carbon–carbon composite design process. Variability in manufacturing process parameters and experimental limitations introduce model parameter uncertainty. This study quantifies the effect of model parameter uncertainty on transverse tensile fracture behavior and proposes a methodology to predict the failure mode based on competing microscale damage mechanisms. Finite element simulations incorporate fiber–matrix interface debonding with cohesive zones and matrix damage with a smeared crack band approach in a unidirectional carbon–carbon composite. Results from a variance-based global sensitivity analysis identifies interfacial and matrix damage parameters as the primary source of variability in fracture behavior. Sobol’ indices indicate that matrix and cohesive zone strengths contribute 94% of the variance in the effective ultimate stress. A local analysis elucidates the relationship between these constituent strength parameters and failure mode by estimating the probability of cohesive, matrix, and mixed-mode dominated failure. Based on the results for 4000 simulations, 93% exhibit mixed-mode or interfacial dominated failure, which underscores the crucial role of fiber–matrix interface debonding in the transverse tensile failure of carbon–carbon composites. These uncertainty quantification results facilitate more efficient model calibration and provide a framework for microstructure-informed failure predictions in the face of manufacturing-induced uncertainty.

36 MATERIALS SCIENCE

Crystal Growth Scale-Up and Stability of RbSr 2 X 5 :Eu Scintillators

The discovery and development of new scintillation materials support national security needs. In these applications, large volumes of scintillator crystals are needed to achieve efficient screening for contraband. Therefore, one important step in the discovery of new scintillator compositions is testing their feasibility for scale-up and their stability. In this work, high-quality Ø22 mm crystals of two new scintillators RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu were grown via the Vertical Bridgman method, and their scintillation properties were characterized. Here, the Ø22 mm RbSr 2 I 5 :Eu crystals could be grown with fast translation rates up to 3.5 mm/h. Both RbSr 2 Br 5 :Eu and RbSr 2 I 5 :Eu had high scintillation performance, including light yields of 46,000 and 61,000 ph/MeV, respectively, for Ø22 × 35 mm crystals. Additionally, properties related to physical stability were investigated, including the coefficients of thermal expansion via high-temperature X-ray diffraction (HTXRD) as well as moisture sensitivity. HTXRD confirmed the absence of solid–solid phase transitions and showed that RbSr 2 Br 5 had minimal thermal expansion anisotropy compared to RbSr 2 I 5 and some other inorganic metal halide scintillators, which favors the growth of large-sized crystals.

36 MATERIALS SCIENCE

Selectivity mechanisms of ion intercalation in Prussian blue analogs

Prussian blue analogs (PBAs) are a family of materials with facile, reversible, and selective ion transport capability for various ions via electrochemical intercalation, owing to their vacancy structure. The large tunable compositional space of PBAs allows for manipulation of intercalation behavior and selectivity by controlling structural vacancy level through choice of transition metal centers and modifications to the synthesis process. However, a lack of understanding of the mechanisms of ion selectivity hinders the material’s design process. Here, for this work, we investigated the origins of ion selectivity using a model PBA, copper hexacyanoferrate, and focused on eight technologically and biologically prominent ions, for which we determined a sequence of selectivity: Rb + > K + > Na + > Ba 2+ > Sr 2+ ≈ Ca 2+ > Mg 2+ > Li + . We provide electrochemical, structural, and redox evidence of strong correlation between the ion identity, the dominant charge-compensating redox, and preferred occupancy site. Specifically, using synchrotron anomalous X-ray diffraction (AXRD), we reveal that monovalent ions exhibit significant association with the corner sites of the unit cell and iron redox, whereas divalent ions display affinity toward the center site with higher ratios of copper redox. Informed by selectivity results, we applied CuHCFe to Li purification and achieved 99.9% purity. Our findings demonstrate an approach to elucidating ion intercalation behavior in order to distinguish and manipulate material properties to optimize separation performance.

Prussian blue analog

Machine learning approach for vibronically renormalized electronic band structures

Here, we present a machine learning (ML) method for efficient computation of vibrational thermal expectation values of physical properties from first principles. Our approach is based on the nonperturbative frozen phonon formulation in which stochastic Monte Carlo algorithm is employed to sample configurations of nuclei in a supercell at finite temperatures based on a first-principles phonon model. A deep-learning neural network is trained to accurately predict physical properties associated with sampled phonon configurations, thus bypassing the time-consuming ab initio calculations. To incorporate the point-group symmetry of the electronic system into the ML model, group-theoretical methods are used to develop a symmetry-invariant descriptor for phonon configurations in the supercell. We apply our ML approach to compute the temperature dependent electronic energy gap of silicon based on density functional theory (DFT). We show that, with less than a hundred DFT calculations for training the neural network model, an order of magnitude larger number of sampling can be achieved for the computation of the vibrational thermal expectation values. Our work highlights the promising potential of ML techniques for finite temperature first-principles electronic structure methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND

Enhanced Catalytic Dechlorination of Polyvinyl Chloride (PVC) and H2 Production Enabled by Synergistic Gaδ+/Ga0 Active Sites in Liquid Metal Particles

Polyvinyl chloride (PVC) is ubiquitous yet challenging to recycle due to its tendency to thermally decompose above 250 °C, releasing toxic, corrosive chlorinated compounds, and its inability to melt. Here, we report a catalytic strategy for PVC upcycling at 160 °C using gallium liquid metal particles (Ga-LMP) featuring a dynamic Ga-GaOOH core–shell architecture. These catalysts enable concurrent dechlorination and hydrogen evolution, yielding up to 7% H2 (based on initial hydrogen atoms in PVC) along with a highly dechlorinated (>95%) carbonaceous solid and aqueous HCl. Mechanistic investigations combining X-ray photoelectron spectroscopy, infrared spectroscopy, solid-state NMR, inelastic neutron scattering, and ab initio molecular dynamics reveal a synergistic interplay between Gaδ+ sites in the GaOOH shell and metallic Ga0 in the core. Cationic Ga initiates C–Cl bond activation and HCl formation, while progressive reduction of the shell exposes Ga0 sites that promote C–H activation and H2 evolution. Control experiments with a Ga salt and bulk Ga liquid metal confirmed that neither oxidation state alone can achieve both transformations efficiently. This work establishes a dynamic dual-site paradigm for liquid metal catalysis, in which the in situ evolution and coexistence of oxidized and metallic species enable sequential and cooperative bond activation pathways. These findings provide a general design principle for novel liquid metal catalysts that target challenging polymer transformations under mild conditions.

Zingg, Benjamin [ORNL] (ORCID:0009000914530153)

Healable Coatings as a Mechanism to Repair Leading Edge Erosion in Wind Energy

Wind turbine blades are highly engineered structures designed to face temperature extremes and high winds. However, erosion of the blade's leading edge and subsequent repair remains a significant and costly challenge for the wind energy industry. Repair of these leading edges can lead to large amounts of downtime for the turbine and significant operational inefficiencies. In this work, the strength of adhesion and healing ability of a commercially available vitrimer (Mallinda's VITRIMAX) was compared to that of a thermoplastic resin, which has previously been demonstrated in wind energy applications (Arkema's Elium) to evaluate their efficacy as surface coatings for wind turbine blades, particularly their leading edges. Vitrimers are a class of inherently reprocessable thermosets, and it was theorized that vitrimer-based leading edge coatings could enable more robust and efficient wind turbine blades with decreased operational downtime and safer maintenance practices. It was found that the VITRIMAX adhered better to the wind blades' surfaces than both the manufacturer's paint and Elium, with increases in pull-off strength of adhesion ranging from 24% to 83% above that of the original paint. Furthermore, the VITRIMAX adhered strongly to the underlying composite of each blade with strength of adhesion values increasing in ranges from 42% to 97% above that of the original paint. Finally, the vitrimer coating showed an 88% decrease in surface roughness compared to end-of-life blade materials, and initial healing demonstrations in which coatings were manually scratched and subsequently healed exhibited an ~84.5% decrease in scratch depths.

Hubbard, Amber [ORNL]

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)