Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Grain boundary segregation”

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.

192 records · Page 5
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

Thermal Cycling and Isothermal Deformation Response of Polycrystalline NiTi: Simulations vs. Experiment

A recent microstructure-based FEM model that couples crystal-based plasticity, the B2<-> MB190 phase transformation and anisotropic elasticity at the grain scale is calibrated to recent data for polycrystalline NiTi (49.9 at.% Ni). Inputs include anisotropic elastic properties, texture and differential scanning calorimetry data, as well as a subset of recent isothermal deformation and load-biased thermal cycling data. The model is assessed against additional experimental data. Several experimental trends are captured - in particular, the transformation strain during thermal cycling monotonically increases and reaches a peak with increasing bias stress. This is achieved, in part, by modifying the martensite hardening matrix proposed by Patoor et al. [Patoor E, Eberhardt A, Berveiller M. J Phys IV 1996;6:277]. Some experimental trends are underestimated - in particular, the ratcheting of macrostrain during thermal cycling. This may reflect a model limitation that transformation-plasticity coupling is captured on a coarse (grain) scale but not on a fine (martensitic plate) scale.

Phase Transformations↗
DOE OSTI2023

Predicting microstructurally sensitive fatigue‐crack path in WE43 magnesium using high‐fidelity numerical modeling and three‐dimensional experimental characterization

Abstract Microstructurally small fatigue‐crack growth in polycrystalline materials is highly three‐dimensional due to sensitivity to local microstructural features (e.g., grains). One requirement for modeling microstructurally sensitive crack propagation is establishing the criteria that govern crack evolution, including crack deflection. Here, a high‐fidelity finite‐element modeling framework is used to assess the performance and validity of various crack‐growth criteria, including slip‐based metrics (e.g., fatigue‐indicator parameters), as potential criteria for predicting three‐dimensional crack paths in polycrystalline materials. The modeling framework represents cracks as geometrically explicit discontinuities and involves voxel‐based remeshing, mesh‐gradation control, and a crystal‐plasticity constitutive model. The predictions are compared to experimental measurements of WE43 magnesium samples subject to fatigue loading, for which three‐dimensional grain structures and fatigue‐crack surfaces were measured post‐mortem using near‐field high‐energy x‐ray diffraction microscopy and x‐ray computed tomography. Findings from this work are expected to improve the predictive capabilities of simulations involving microstructurally small fatigue‐crack growth in polycrystalline materials.

Engineering↗
NASA NTRSDate not supplied

Improvements to RANS Modeling for Aeroheating Predictions on Blunt Bodies

Accurate predictions of aeroheating are critical for designing thermal protection systems for planetary entry vehicles. For larger vehicles, turbulence in the boundary layer can substantially increase convective heating. This turbulence must be accurately modeled to ensure the thermal protection system is sufficient. The majority of hypersonic turbulence model development and validation focuses on boundary layers developing over flat-plates or sharp cones; these cases are substantially different than the boundary layer that develops over the heatshield of a blunt body traveling at hypersonic speeds. Planetary missions often use blunt body geometries, such as the 70-degree sphere-cone favored by Mars missions or the 45-degree sphere-cone planned for the upcoming DAVINCI mission. Due to smaller vehicle size and the lower velocities in the stagnation region, planetary entry vehicles have relatively low Reynolds numbers. Surface curvature and high enthalpy gradients create additional challenges. These difficulties must be addressed to obtain high accuracy needed for the ambitious planetary missions in the upcoming decade. This work focuses on both assessing and improving Reynolds-averaged Navier-Stokes (RANS) turbulence models for blunt-body geometries typical of planetary entry vehicles, with a focus on one-equation and two-equation formulations.

Mars2020↗
NASA NTRS2026

Remote Sensing-Driven Hydrodynamic Modeling in Data-Scarce Regions: Integrating ICESat-2, Sentinel-2, SWOT and Re-analysis Models for Coastal Monitoring

Hydrodynamic models in coastal and estuarine systems are typically constrained by sparse bathymetry, boundary, and validation data, especially in regions where field campaigns are costly or impractical. Here we develop and test a fully satellite-driven framework for hydrodynamic modeling in South Africa’s Langebaan Lagoon without using any local in situ measurements. Bathymetry is derived by training multispectral Sentinel-2 reflectance against ICESat-2 ATL24 photon-derived depths using an XGBoost model optimized with Bayesian search. The final satellite-derived bathymetry reproduces independent ATL24 points with RMSE = 0.45 m and R 2 = 0.97. This bathymetry was used in a depth-averaged Delft3D Flexible Mesh model driven at the open boundary by TPXO tidal harmonics and by ERA5 winds. We validate modeled water surface elevation against 16 SWOT low-rate (250 m, unsmoothed) passes in 2023. SWOT–model comparisons yield an overall RMSE of 0.11 m and R 2 = 0.61, with typical point differences <0.10 m (∼7% of the 1.5 m tidal range), and showed consistent spatial gradients in water level from the offshore boundary, through Saldanha Bay, and into the lagoon. At the offshore boundary, TPXO and SWOT sea surface heights agree closely (R 2 = 0.86). A simple phase adjustment of ∼26,min between TPXO and SWOT lowers the RMSE from 0.18,m to 0.11,m, showing that phase offset accounts for some of the discrepancy, with additional errors likely linked to non-tidal signals. Our results demonstrate that combining passive optical, photon-counting LiDAR, radar interferometry, and global tidal/atmospheric models enables robust, transferrable hydrodynamic modeling in data-scarce coastal systems, offering a cost-effective pathway for monitoring.

ICESat-2↗
NASA NTRS2026

LAI Assimilation Schedules to Constrain Uncertain Cultivars and Soils in CERES-Maize

In anticipation of large-domain crop model applications where precise local configuration and calibration is not possible, we describe benefits and potential drawbacks of employing a crop pest module to achieve leaf area index (LAI) assimilation into a high performing CERES-Maize crop model configuration at a field experiment site in Perry, Iowa. Simulation experiments explore the use of MODIS satellite-derived LAI to constrain and adjust LAI to counter imprecise cultivar and soil configurations often occurring in the absence of high-quality local information. Simulations using single-day, window, and continuous LAI replacement across 5 cultivars and 2 soil calibration approaches for nine corn rotation years from 2004 to 2020 led to different yield outcomes and reverberations throughout the field environment. Evaluating variance and mean bias, results indicate minimal interventions in early vegetative and grain-filling stages were more beneficial than use of continuous LAI adjustments, as they minimized disruptions to the internal resource balances governing plant stresses and grain production. LAI adjustment was particularly helpful in constraining growth related to uncertain thermal unit requirements and leaf tip appearance rates (P1 and PHINT cultivar parameters, respectively). Findings underscore the need to assimilate additional state variables to ensure internal biophysical coherence. This approach shows promise for applications spanning wider domains with prediction time pressure where detailed configuration, more complex assimilation methods, or recalibration of crop model parameters may not be practical.

Phenology↗
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↗
NASA NTRS2014

Analytical Model for Steady Flow through a Finite Channel with One Porous Wall with Arbitrary Variable Suction or Injection

This paper presents an exact solution of two-dimensional laminar flow through a finite length channel with one porous wall. It improves upon previous solutions by (1) satisfying the no-slip boundary condition at the channel dead end, (2) adding a turbulent term to the porous wall boundary condition, (3) allowing for arbitrary variable suction or injection across the porous wall, and (4) model validation against new cryogenic liquid hydrogen and oxygen experimental data. Of particular interest in the current work is the modeling of cryogenic propellant flow through a porous liquid acquisition device (LAD) screen and channel inside a propellant tank. First, a detailed review of the literature is presented for previously attempted solutions to channel flow with one porous wall. Next, the governing equations, boundary conditions, and model assumptions are used to derive the analytical flow solution and present general model results for pressure and velocity fields within the channel. Then, the model solution is compared with horizontal LAD channel flow data in liquid oxygen as well as vertical LAD channel flow data in an inverted outflow configuration in liquid hydrogen. Model results are used to update the static cryogenic bubble point pressure model with a dynamic bubble point term which factors in enhanced convection and cooling at the screen during propellant outflow. Convective heat transfer at the LAD screen during outflow is also quantified by comparing model and data. The new analytical flow solution with the dynamic bubble point model is shown to compare well with available cryogenic experimental data

Navier Stokes Equations↗
NASA NTRS2026

Surprisingly Enriched CO and CH4 in Halley-type Comet 13P/Olbers: Clues to Its Interstellar Heritage

The native ice composition of Halley-type Comet 13P/Olbers (hereafter 13P) was measured post-perihelion, on UT 2024 August 16 using iSHELL, the 1–5 μm cross-dispersed, high-resolution facility spectrograph at the NASA InfraRed Telescope Facility. At that time, the heliocentric distance (Rh) of 13P was 1.37 au and the geocentric distance (Δ) was 1.99 au. Using a slit that delivered spectral resolving power RP λ/Δλ 4.5 × 104, two instrument settings targeted rovibrational emissions of seven parent molecules released into the coma through sublimation of ices in the nucleus. These settings provided broad and contiguous or nearly contiguous spectral coverage, one simultaneously measuring H2O, CO, and OCS, and the other simultaneously measuring CH4, C2H6, CH3OH, and H2CO. Meaningful abundances were obtained for all seven species. Our study revealed substantial enrichments of CO and CH4, the two most volatile species systematically targeted at IR wavelengths in comets. Relative to their respective mean abundance ratios among measured comets from the Oort cloud reservoir, CO in 13P was enriched by a factor of about six and CH4 was enriched by a factor of about three. Together with near-average abundances measured for CH3OH and C2H6, this suggests large endowments of CO and CH4 on interstellar grains, yet it could indicate inhibited surface chemistry on the pre-cometary grains subsequently contained in the nucleus or 13P. Assuming that this represents the primordial composition of 13P, it could reflect very low temperatures in its natal environment (∼10 K or lower).

Michael A Disanti↗
DOE OSTI2026

Ionic Precursors Transformed Into Vinyl Acetate Synthesis Catalyst via Reaction-Driven Restructuring

Conventional preparation of supported bimetallic catalysts relies on solution-mediated metal salt immobilization and pre-formation of alloy nanoparticles before reaction. Here, we report a fundamentally different synthesis strategy of using a physical mixture of salt precursors to generate an active catalyst during reaction. The catalytic structure is generated in situ from Pd3(OAc)6, Au(OH)3, and KOAc through H2 treatment and reaction-driven restructuring under vinyl acetate monomer (VAM) synthesis conditions. Ascertained from in situ X-ray diffraction and operando infrared spectroscopy analyses, reduction treatment produces segregated Pd and Au domains, and subsequent exposure to a VAM reaction mixture triggers dynamic extraction of Pd from the metal surface. This latter process, mediated by acetate-assisted redox cycles, facilitates Pd migration toward Au domains to form a near-surface localized Pd50Au50 alloy phase. Monometallic Pd domains serve as a reservoir of Pd to the alloy phase, leading to and sustaining a more Pd-enriched active surface and a higher population of accessible Pd sites, compared to a conventionally prepared K-PdAu/SiO2 catalyst. Consequently, this leads to a twofold increase in the VAM formation rate, demonstrating highly active bimetallic catalysts can be generated through the gas-phase treatment of physically mixed ionic precursors.

Cha, Byeong Jun [Rice University]↗
DOE OSTI2026

Equilibrium and non-equilibrium effects in high pressure phase transformations of carbon

The behavior of carbon in the range 1–100 GPa and 1–10 kK is central to problems in planetary interiors, inertial confinement fusion targets, and high-pressure synthesis of carbon-based materials, but experiments in this regime are difficult and often provide only indirect constraints on phase behavior. As a result, phase boundary loci, structure, and limits of metastability at high pressure remain uncertain. In this work, machine-learning enhanced atomistic simulations are used to address this knowledge gap. We determine the melt line up to 100 GPa, the graphite-diamond phase boundary up to the melt line, and analyze structure of the coexisting phases. We show that the coexisting liquid evolves smoothly with pressure without evidence for a first-order liquid–liquid transition. Orientation-resolved graphite melting simulations indicate that basal-plane interfaces develop a dewetting layer and undergo layer-by-layer melting, producing kinetic hysteresis and an apparent orientation dependence of the melt line. Non-equilibrium quenches from the melt are used to construct a kinetically limiting graphite–diamond phase boundary for rapid quenches from above the melt line, and show that graphite is metastable at pressures of up to ≈ 25 GPa. These results provide bounds on equilibrium and metastable behavior in carbon relevant for interpreting high-pressure experiments and for designing synthesis pathways to specific carbon microstructures.

Lyu, Yanjun [Department of Materials Science and E↗
DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗
DOE OSTI2023

Fast and Non‐equilibrium Uptake of Hydrogen by Pd Icosahedral Nanocrystals

Abstract We report for the first time that Pd nanocrystals can absorb H via a “single‐phase pathway” when particles with a proper combination of shape and size are used. Specifically, when Pd icosahedral nanocrystals of 7‐ and 12‐nm in size are exposed to H atoms, the H‐saturated twin boundaries can divide each particle into 20 smaller single‐crystal units in which the formation of phase boundaries is no longer favored. As such, absorption of H atoms is dominated by the single‐phase pathway and one can readily obtain PdH x with anyx in the range of 0–0.7. When switched to Pd octahedral nanocrystals, the single‐phase pathway is only observed for particles of 7 nm in size. We also establish that the H‐absorption kinetics will be accelerated if there is a tensile strain in the nanocrystals due to the increase in lattice spacing. Besides the unique H‐absorption behaviors, the PdH x (x=0–0.7) icosahedral nanocrystals show remarkable thermal and catalytic stability toward the formic acid oxidation due tothe decrease in chemical potential for H atoms in a Pd lattice under tensile strain.

Chemistry↗
DOE OSTI2026

Interfacial magnetization enhancement and the perpendicularly magnetized component in Fe16N2 thin films

A detailed study of depth-resolved magnetization and microstructure of Fe16N2 thin films on MgO (001) and MgAl2O4 (001) substrates and Fe and Cr seed layers is presented. Two aspects of the magnetic properties of Fe16N2 thin films are discussed. First, magnetization enhancement at the interface is observed. Strain and nitrogen deficiency are discussed as possible interfacial mechanisms contributing to this enhancement. Second, the perpendicularly magnetized component (PMC) is identified in Fe16N2 thin films. Correlation with microstructural observations suggests that the PMC is associated with V-shaped grains that are not fully confined within the continuous Fe16N2 layer.

Hang, Xudong [University of Minnesota]↗
NASA NTRS2026

Coma Physics of an Interstellar Object: JWST Spatial-Spectral Mapping of 3I/ATLAS

We report a survey of molecular emission from cometary volatiles using the James Webb Space Telescope (JWST) toward interstellar object 3I/ATLAS carried out on UT 2025 December 22 and 23 at a heliocentric distance (\rh{}) of $2.37-2.41$ au. These measurements of CO, \ce{CO2}, \ce{H2O}, \ce{CH3OH}, and \ce{CH4} sampled molecular chemistry in 3I/ATLAS as it receded from its encounter with our Sun and entered the vicinity of the \ce{H2O} ice line --- the region between \rh{} = $2-3$ au where the temperature becomes too low for H$_2$O to vigorously sublime and CO and \ce{CO2} begin to control the overall activity. CO was the most abundant molecule, followed by \ce{H2O} and \ce{CO2}, whose molecular abundances with respect to CO were $(40.5\pm3.1)\%$ and ($41.6\pm0.3)\%$, respectively. This work presents spatial-spectral maps of column density and rotational temperature as a function of distance from the nucleus for all detected species. The spatial distributions of both quantities were highly anisotropic for the apolar species in the coma of 3I/ATLAS, yet were more nearly symmetric for the polar molecules. These results demonstrate how volatiles were segregated in the nucleus ices of 3I/ATLAS and reveal heating and cooling mechanisms in its coma. Derived maps of the ortho-to-para ratio (OPR) for \ce{H2O} were flat with increasing distance from the nucleus and consistent with a coma-averaged value $\mathrm{OPR}=2.7\pm0.2$, slightly less than the expected equilibrium value of three.

Nathan X Roth↗
NASA NTRSDate not supplied

Evaluations of Additively Manufactured Superalloy Lattice Blocks

Cast lattice block structures made up of high-temperature superalloys were previously shown to offer advantages in density, efficiency, and tailored properties. However, inconsistencies in casting die-fill, composition, and microstructure were identified as persistent issues with the casting process path. Three candidate lattice designs were screened through additive manufacturing of a superalloy. Laser powder bed fusion was applied with superalloy powder to produce lattice blocks for three-dimensional octahedral unit-cell designs of three heights. Properties of the structural designs versus additively manufactured lattices were compared, including structural dimensions, densities, defects, grain structures, crystallographic texture, and surface residual stresses. Findings indicate this concept could produce useful superalloy structures for high-temperature-capable static components requiring tailored density, and thermophysical and mechanical properties.

superalloy↗
NASA NTRS2026

Developing a Pyrolysis Gas Thermal Blocking Model for Reentry Demise

In NASA’s Object Reentry Survival Analysis Tool (ORSAT), aerodynamic drag and aerothermal heating coefficients are computed for each of the free-molecular, continuum, and transitional flow regimes using analytical and semi-analytical methods. These heating coefficients were derived for typical metallic materials that melt and do not have a strong gas-phase contribution to the flow in the boundary layer. Modern satellites typically feature fiber-reinforced polymer (FRP) components, such as solar array booms, facesheets of sandwich panels, or overwraps for composite-overwrapped pressure vessels (COPV). These FRP materials do not behave the same as metals in the reentry environment, but instead will pyrolyze and develop significant volumes of gas into the boundary layer. Accurately predicting the reentry demise of FRP components is critical to assessing the reentry casualty risk for modern spacecraft. Research in recent years has shown that this demisability can depend heavily on how the expulsion of gaseous pyrolysis products through the outer surface of the material affects the heat flux at the surface. The ODPO has been developing a reduced-order model of the effect of pyrolysis gas blowing on the heat flux based on correlations between a blowing factor and a non-dimensional heat flux to be incorporated in the upcoming version 7.3 of the Object Reentry Survivability Analysis Tool (ORSAT). This presentation discusses the progress of this development project and the challenges remaining for generalizing the model across families of FRP materials.

Benton Greene↗
NASA NTRS2026

Unraveling the Complex History of Aqueous Alteration and Metamorphism in CV Chondrite, Camel Donga 040

Camel Donga (CD) 040 is a CV breccia collected from Australia in 1988. Initial studies of CD 040 revealed two lithologies, both containing chondrules up to 1 mm in diameter (some with fine-grained rims) and abundant troilite, but each with remarkably different olivine compositions. Lithology 1 was characterized as petrologic type 3 material, while lithology 2 was described as being “heated to ~1100°C”. CD 040 was initially classified as an ungrouped C3; however, subsequent analyses of oxygen isotopes and bulk rock compositions confirmed that both lithologies were CV chondrite material. Though CD 040 was one of the earliest examples of a CV chondrite containing metamorphosed material, it was never extensively studied. Here we report the results of a detailed study of CD 040, focusing primarily on the petrography and mineral chemistry of metamorphosed lithology.

T L Dunn↗