SEARCH · Search NASA

Results for “surface model”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

On the role of α-alumina in the origin of life: Surface-driven assembly of amino acids

We investigate the hypothesis that mineral/water interfaces played a crucial catalytic role in peptide formation by promoting the self-assembly of amino acids. Using classical molecular dynamics simulations, we demonstrate that the α-alumina(0001) surface exhibits an affinity of 4k B Tfor individual glycine or GG dipeptide molecules due to hydrogen bonds. In simulations with multiple glycine molecules, surface-bound glycine enhances further adsorption, leading to the formation of long chains connected by hydrogen bonds between the carboxyl and amine groups of glycine molecules. We find that the likelihood of observing chains longer than 10 glycine units increases by at least five orders of magnitude at the surface compared to the bulk. This surface-driven assembly is primarily due to local high density and alignment with the alumina surface pattern. Together, these results propose a model for how mineral surfaces can induce configuration-specific assembly of amino acids, thereby promoting condensation reactions.

Science & Technology - Other Topics↗

DOE OSTI2025

Dataset, Code, and Models for Training Deep Learning Potentials for Low Temperature Plasma-Surface Interactions

This repository contains datasets, training scripts, and finished models, and test simulations used in the development of DeepREBO— a machine-learned interatomic potential trained to emulate the REBO2 empirical potential. The data was generated to study deep potential development for simulations of plasma-surface interactions. It uses an active learning framework, starting from a minimal dataset and iteratively expanding it. Included are those generated datasets, the trained models, and simulations used to evaluate the performance of the training process. This resource supports reproducibility and provides a reference framework for training deep potentials in plasma-surface interaction studies.

active learning↗

DOE OSTI2025

HDG-2 & VIC Design Updates

HDG-2 and VIC Design updates. Including background, HDG-1 challenges, HDG-2 Mitigation efforts, Effect of updated ring model, Pressure boundary surface finish, and conclusions. Vendor Irradiation Capsule (VIC) to include, background, design overview, testing and neutronics results, and models.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Advancing process-based flood frequency analysis for assessing flood hazard and population flood exposure

Recent studies have showcased the use of process-based hydrological models with Stochastic Storm Transposition (SST) techniques to conduct Flood Frequency Analysis (FFA). This framework, referred hereby FFA-SST, has proved to be a robust strategy to estimate peak flows of specific annual exceedance probability (e.g., 100-year peak flow) that can reflect natural and anthropogenic disturbances, including changes in land use and meteorological patterns. With the objective of advancing the FFA-SST framework, this study presents for the first time the use of an Integrated Surface-Subsurface Hydrological Model (ISSHM) to conduct FFA-SST by extending the analysis from peak flow responses to flood extent, enabling a unique view and analysis of flood hazard and population flood exposure at the basin scale. As a proof-of-concept, we used the ISSHM, Advanced Terrestrial Simulator (Amanzi-ATS) model, and the SST model, RainyDay, to conduct FFA-SST by simulating the flood response to 5,000 annual synthetic storm events in a 2,227 $km^2$ Southeast Texas watershed. We demonstrate that ATS, without site-specific calibration, provides a robust process-based representation of peak flows, flood extent, streamflow, evapotranspiration, soil moisture content, and water storage changes. Our results and analyses, covering frequency curves up to a 500-year return period for peak flows, basin inundation fractions, and the number of people exposed to flooding, offer a unique perspective to analyze flood impacts across spatial scales. Overall, this study provides critical insights for flood risk management by extending the FFA-SST framework to include both flood hazard and population flood exposure analyses at the basin scale. Such an approach will empower stakeholders and disaster emergency agencies with a more comprehensive understanding of flood impacts across the entire basin domain, facilitating informed decision-making for flood risk assessment and management.

58 GEOSCIENCES↗

DOE OSTI2025

Leptodermic corrections to the TOV equations and nuclear astrophysics within the effective surface approximation

The macroscopic model for a neutron star (NS) as a liquid drop at the equilibrium is used to extend the Tolman-Oppenheimer-Volkoff (TOV) equations taking into account the gradient terms responsible for the system surface. The parameters of the Schwarzschild metric in the spherical case are found with these surface corrections to the known leading (zero) order of the leptodermic approximation a/R << 1, where a is the NS effective-surface (ES) thickness, and R is the effective NS radius. The energy density $\mathscr{E}$ is considered in a general form including the functions of the particle number density and of its gradient terms. The macroscopic gravitational component $Φ$(ρ) of the energy density is taken into account in the simplest form as expansion in powers of $ρ$ – $\overline{ρ}$, where $\overline{ρ}$ is the saturation density, up to second order, in terms of its contributions to the separation particle energy and incompressibility. Density distributions ρ across the NS ES in the normal direction to the ES, which are derived in the simple analytical form at the same leading approximation, was used for the derivation of the modified TOV (MTOV) equations by accounting for their NS surface corrections. As a result, the MTOV equations are analytically solved at first order and the results are compared with the standard TOV approach of the zero order.

Magner, A. G. [Institute for Nuclear Research, Kyi↗

DOE OSTI2024

The Shortwave Cloud-SST Feedback Amplifies Multi-Decadal Pacific Sea Surface Temperature Trends: Implications for Observed Cooling

Climate models struggle to produce sea surface temperature (SST) gradient trends in the tropical Pacific comparable to those seen recently in nature. Here, we find that the magnitude of the cloud-SST feedback in the subtropical Southeast Pacific is correlated across models with the magnitude of Eastern Pacific multi-decadal SST variability. A heat-budget analysis reveals coupling between cloud-radiative effects, circulation, and SST gradients in driving multi-decadal variability in the Eastern Pacific. Using this relationship and observed feedback estimates, we find that internal Eastern Pacific SST variability is underestimated in most models. Adjusting for model bias increases the likelihood of generating a cooling trend at least as large as observations in preindustrial control simulations by ~56% on average. If models underestimate climate “noise,” as our results suggest, this bias should be accounted for when attributing the relative importance of forced versus unforced changes in the climate.

58 GEOSCIENCES↗

DOE OSTI2025

Investigation of Electrocatalytic CO 2 Reduction on MXene Materials via First-Principles Simulations

Computational studies of CO 2 reduction to yield various products were carried out on the basal plane and edges of different MXene materials. The impact of vacancies upon Mo 2 TiC 2 T x , W 2 TiC 2 T x , and Ti 3 C 2 T x (T x = O and OH) was also examined for both edge and basal sites. Initial calibrations were carried out to generate surfaces with optimal oxide/hydroxide ratios and proper termination sites upon which various vacancy sites were explored to ensure an accurate model of the surfaces’ resting states. From this work, Mo 2 TiC 2 O was determined to exhibit the lowest theoretical overpotential for methane as determined by volcano plot analyses. At a large enough vacancy concentration, the CO 2 reduction reaction (CO 2 RR) is predicted to outcompete the hydrogen evolution reaction (HER) as the predominant reaction on the surface. In conclusion, when examining the edge of Mo 2 TiC 2 O, stronger CO 2 binding was exhibited to such an extent that the reaction was predicted to terminate after the generation of formate on the edge.

Lund, Colton [Argonne National Laboratory (ANL), A↗

DOE OSTI2025

Cross-Scale Catalyst Modeling Applied to H 2 Storage and Release via Formic Acid

Here, we propose the Systems-to-Atoms (S2A) modeling framework that integrates the kinetics of reaction chemistry and structural configurations across various length scales with the aim of establishing a versatile template for multiscale modeling of reactive flow problems and to predict the operando activity of catalyst materials. The approach encompasses a microkinetic model to analyze surface reactions on individual facets of catalyst nanoparticles coupled with the computation of average surface reaction rates for catalyst nanoparticles of specific size distributions. Macro-homogeneous surface reaction kinetics are derived as a function of catalyst loading and used as input parameters for the continuum-scale reactor model. The cross-scale framework enables the optimization of catalyst utilization through reactor design and operating strategy. To demonstrate the framework, we studied the storage and release of hydrogen from formic acid, a promising liquid organic hydrogen carrier (LOHC), over Pd, Pt, and Cu catalysts. The framework predicts observed trends in formic acid dehydrogenation activity for catalysts with comparable weight loadings and metal particle diameters, demonstrating satisfactory quantitative alignment. Finally, the seamless transmission of parameter uncertainties between scales is also discussed.

08 HYDROGEN↗

DOE OSTI2026

Community-Engaged Modeling of Urban Flood Adaptation Pathways

Climate change is intensifying the hydrologic cycle, leading to more frequent and severe rainfall-driven (pluvial) flooding in urban areas. In the mid-Atlantic US cities, aging and under-designed stormwater infrastructure is increasingly strained by these events, resulting in recurring damage to property and disruptions to transportation networks. In this study, we combine community engagement with hydrologic modeling to develop and evaluate potential urban flood adaptation strategies. Over a three-year period, local technical experts and community representatives met regularly to discuss flooding concerns, identify priorities, and co-develop adaptation strategies. These discussions informed the development of an urban flooding model (EPA Storm Water Management Model) for the Baltimore Harbor watershed, the focus location of this study. The flooding model integrates complex surface and subsurface stormwater infrastructure data, local expert knowledge, and community insights. We simulate stakeholder-prioritized adaptations, such as green and gray infrastructure strategies. Model results demonstrate that enhanced infrastructure maintenance is the most effective adaptation for reducing flood depths, but has varied effects across the watershed, and can increase flooding in some locations. Spatially concentrated greening provides limited benefit to the watershed as a whole, but moderate benefit in community priority areas. Together, these adaptations have the potential to reduce flood depths by as much as 58% in some locations, greatly reducing property damage and mobility impacts, primary concerns of stakeholders. Future work will implement robust optimization tools to search for adaptations which meet stakeholder objectives and perform highly under varied future climate conditions. This work contributes to the expanding literature on collaborative modeling and demonstrates that community-engaged approaches can enhance model credibility and generate more actionable insights for communities seeking to strengthen climate resilience.

Spangler, Ava [Pennsylvania State University] (ORC↗

DOE OSTI2026

Community-Informed Urban Flood Modeling for Impact Mitigation

Baltimore MD↗

DOE OSTI2026

A bespoke model of Arctic river basins based on hillslope delineation: Model Archive

This dataset is a model archive of the paper A bespoke model of Arctic river basins based on hillslope delineation (in prep), which introduces a watershed decomposition and parameterization method for large scale permafrost hydrology simulation. With this dataset, this study aims to address the research question: whether a computationally efficient hillslope-based modeling framework can reliably simulate discharge at Arctic river-basin scales. This dataset contains model input and output data for five modeling scenarios at a study site located in the Sagavanirktok River basin. The five modeling scenarios include three modeling cases under temperate conditions using full 3D, decomposed 3D, and decomposed 2D modeling strategies; and two modeling cases under actual Arctic conditions with permafrost using full 3D and decomposed 2D modeling strategies. Simulations were performed using the Advanced Terrestrial Simulator (ATS, v1.6 for three temperate scenarios and v1.5 for two Arctic scenarios), a physics-rich integrated surface–subsurface hydrologic model with cryo-hydrology features. For the three temperate models, simulations were conducted for the period of 10/01/1993 - 09/30/2002; and for the two Arctic models, simulations were conducted for the period of 01/01/1994 - 12/31/2002. To facilitate reproducibility of simulations, all datasets are organized hierarchically. The dataset contains: (1) Mesh files (.exo) for full 3D model, decomposed 3D models, and decomposed 2D models, located in huc/190604020802_gauge15906000/mesh/. Mesh files can be visualized through Paraview or read by Python. (2) Climate forcings (.h5) for full 3D model and decomposed 3D/2D models are located in huc/190604020802_gauge15906000/daymet_onePiece/, and huc/190604020802_gauge15906000/vp_pr_revised_daymet_1980_2006_with_wind/ separately. Accessible by Python. (3) Raw measured gage discharge (.csv) from USGS, located in huc/190604020802_gauge15906000/gaged_basin15906000_discharge_usgs/. Accessible by Python. (4) Delineated subdomain raster (.tif) and shape files (.shp), and the final parameterized results (.npy) for decomposed models, located in huc/190604020802_gauge15906000/data_preprocessed-meshing. Accessible by Python. (5) Temperate models are located in nonpermaf_huc190604020802_gauge15906000/, which includes three cases: decomposed 2D models (inside model_0*-hillslope_*), decomposed 3D models (inside model_1*-subcatchment_*), and full 3D model (inside model_2*-onepiece_*). Two step spin-up results (checkpoint_final.h5) are located in model_*1-*_spinup_steadystate and model_*2-*_spinup_cycle, separately, which are used to initialize real transient models. The input files (.xml) and output results (.dat) of the real transient models are located in model_*3-*_transient/. Especially, for two example hillslope models (ID=-11 and 11), additional h5py files are included in model_03-hillslope_transient/hillslope-11/, model_03-hillslope_transient/hillslope11, model_13-subcatchment_transient/subcatchment-11/, model_13-subcatchment_transient/subcatchment/11, respectively, which are used to plot the saturation figure (Figure 5) in the manuscript. Accessible by Python. (6) Arctic models are located in huc190604020802_gauge15906000/, which includes two cases: decomposed 2D models (inside model_04-hillslope_transient), and full 3D model (inside model_05-onepiece_transient_mannp1_ra). Three step spin-up results (checkpoint_final.h5) are located in model_01-column_freezeup/, model_02-column_spinup/, model_03-hillslope_spinup/, respectively, which are used to initialize real 2D transient hillslope models. The input files (.xml) and output results (.dat) of transient 2D hillslope models are located in model_04-hillslope_transient/. The input files (.xml) and output results (.dat) of the full 3D transient model is located in model_05-onepiece_transient_mannp1_ra/. The full 3D transient model is initialized by model_02-column_spinup/. Accessible by Python. (7) The MOSART routed discharge results (.csv) under Arctic conditions is located in huc190604020802_gauge15906000/MOSART/. Accessible by Python. (8) All Python codes (.py) used to parameterize full 3D model to decomposed 2D models are located in script/. These codes fit with watershed workflow (a watershed delineation tool) v1.4 under the branch gaob/v1.4 from https://github.com/gaobhub/watershed-workflow.git.

EARTH SCIENCE > CRYOSPHERE↗

DOE OSTI2025

The Role of Snowmelt and Subsurface Heterogeneity in Headwater Hydrology of a Mountainous Catchment in Colorado: A Model‐Data Integration Approach

Mountainous headwater streams are sustained by both snowmelt‐driven streamflow and groundwater discharge in the Upper Colorado River Basin. However, predicting headwater stream discharge magnitude and peak flow timing is challenging in mountainous terrains, where snowmelt rates vary with vegetation type and elevation, and heterogeneous subsurface physical properties influence groundwater storage and its release. We used a model‐data integration approach to investigate the roles of snowmelt and subsurface structure in stream discharge and groundwater level. We ran an ensemble of 100 integrated surface‐subsurface hydrologic models for a mountainous headwater catchment near Crested Butte, Colorado, USA. We also evaluated and calibrated these models against observed data sets, including snow depth measurements using distributed temperature probes, stream discharge, and groundwater levels. Calibration with multiple data sources using neural density estimators has further constrained uncertainty in subsurface properties and snowmelt rates. Results indicated that observed slower snowmelt rates in evergreen forests delayed the peak flow and baseflow onset. In upstream areas with lower subsurface permeability, water was stored within the subsurface but was not released as interflow or shallow groundwater flow, and thereby not contributing to downstream streamflow during recession limb periods. Double peaks in groundwater occurred in areas with spatial subsurface heterogeneity, in our case due to the contrast between granodiorite and Mancos shale. These process‐based insights into groundwater and snowmelt dynamics in mountainous headwaters will help improve predictions of headwater hydrology.

Wang, Lijing [University of Connecticut, Storrs, C↗

DOE OSTI2025

Parametric Finite Element Analysis of Naturally Corroded Steel Specimens Using 3D Surface Laser Scans

Corrosion is considered a uniform thickness reduction design guideline of the maritime industry. However, additionally, the corroded and irregular morphology of the surface affects the steel's load-bearing capacity and its impact on the strength and elongation behaviour of the steel is not yet fully understood. These effects on the local behaviour of steel structures under tensile loading were investigated with tensile tests on naturally corroded steel specimens and nonlinear finite element simulations including the corroded surface morphology with a uniform surface idealation. The models also include the deformed specimen shape. The developed approach led to highly accurate parametric finite element models predicting the ultimate tensile strength and longitudinal position of fracture. The results show that all included aspects are essential for accurate simulations, while solely the maximum available surface resolution was not as decisive.

corrosion↗

DOE OSTI2025

Microstructure Scale Lithium-Ion Battery Modeling: Part I. On Through-Plane Heterogeneity, Impact of Mesh Representation, and Differences between Macro- and Microscale Models

Li-ion battery performance and degradation are strongly correlated with the electrode microstructures and can be modeled at different scales, each with their own limitations. Herein, we compare predictions achieved with a macro- and a micro-scale model, that is, respectively, neglecting or considering the microstructural heterogeneity of the composite electrodes, on virtual numerically generated and real microstructures. While both models are in relative agreement at the low charge rates, differences arise for fast charging scenarios and especially for the real, highly heterogenous, microstructures. The microscale model predicts that electrolyte concentration saturation and depletion, respectively, at the back of the cathode and of the anode are exacerbated, and that lithium plating occurs earlier for real microstructures. The present work also indicates that the mesh representation significantly impacts the microscale model predictions, and consequently that microscale models should add surface area as a parameter to consider explicitly surface roughness. This article is the first of a series, with subsequent entries further investigating in-plane heterogeneities, lithium plating, and the impact of microstructure representativity on model predictions.

25 ENERGY STORAGE↗

DOE OSTI2024

The HHV-6B U20 glycoprotein binds ULBP1, masking it from recognition by NKG2D and interfering with natural killer cell activation

Human Herpesvirus 6B (HHV-6B) impedes host immune responses by downregulating class I MHC molecules (MHC-I), hindering antigen presentation to CD8+ T cells. Downregulation of MHC-I disengages inhibitory receptors on natural killer (NK) cells, resulting in activation and killing of the target cell if NK cell activating receptors such as NKG2D have engaged stress ligands upregulated on the target cells. Previous work has shown that HHV-6B downregulates three MHC-like stress ligands MICB, ULBP1, and ULBP3, which are recognized by NKG2D. The U20 glycoprotein of the related virus HHV-6A has been implicated in the downregulation of ULBP1, but the precise mechanism remains undetermined. We set out to investigate the role of HHV-6B U20 in modulating NK cell activity. We used HHV-6B U20 expressed as a recombinant protein or transduced into target cells, as well as HHV-6B infection, to investigate binding interactions with NK cell ligands and receptors and to assess effects on NK cell activation. Small-angle X-ray scattering was used to align molecular models derived from machine-learning approaches. We demonstrate that U20 binds directly to ULBP1 with sub-micromolar affinity. Transduction of U20 decreases NKG2D binding to ULBP1 at the cell surface but does not decrease ULBP1 protein levels, either at the cell surface or in toto. HHV-6B infection and soluble U20 have the same effect. Transduction of U20 blocks NK cell activation in response to cell-surface ULBP1. Structural modeling of the U20 – ULBP1 complex indicates some similarities to the m152-RAE1γ complex.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2026

Robustness of the Ensemble Score Filter to the Type of Assimilated Observation Networks

Recent advances in data assimilation (DA) have focused on developing more flexible approaches that can better accommodate nonlinearities in models and observations. However, it remains unclear how the performance of these advanced methods depends on the observation network characteristics. In this study, we present initial experiments with the surface quasi‐geostrophic model, in which we compare a recently developed ensemble filter using score‐based diffusion models with the standard Local Ensemble Transform Kalman Filter (LETKF). Our results show that the analysis solutions respond differently to the number, spatial distribution, and nonlinear fraction of assimilated observations. We also find notable changes in the multiscale characteristics of the analysis errors. Given that standard DA techniques will eventually be replaced by more advanced methods, we hope this study sets the ground for future efforts to reassess the value of Earth observing systems in the context of newly emerging algorithms.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

The Atlantic Meridional Overturning Circulation’s Response to CO2 Increase: Assessing the Roles of Surface Flux and Oceanic Advection Feedbacks

Abstract The Atlantic meridional overturning circulation (AMOC) is projected to slow down in climate models due to greenhouse gas emissions. However, the physical mechanisms determining the rate of the projected AMOC slowdown remain unclear. Accordingly, this study isolates the roles of oceanic advection and surface flux feedbacks that might accelerate or decelerate the AMOC’s weakening using carbon dioxide (CO 2 ) quadrupling simulations in the CESM1.2 model. Surface flux feedbacks are isolated in partially coupled experiments in which either all surface flux components or the momentum flux responses to AMOC’s weakening that might provide feedback are suppressed, while a tracer decomposition of ocean density anomalies isolates the advection feedbacks. Comparing the ocean density components in the experiments shows that the AMOC’s response is initially determined by CO 2 -induced anomalous surface heat fluxes, and afterward, feedbacks determine its response. In the fully coupled case, surface heat flux feedback strongly promotes AMOC slowdown and causes its near shutdown, while a weaker but active AMOC is maintained when the surface flux feedback is inhibited in the partially coupled case. The positive surface heat flux feedback works by canceling out the negative oceanic heat advection feedback on deep-water formation in the subpolar North Atlantic (SPNA). With the heat advection feedback thus reduced, the positive salinity advection feedback becomes the dominant contributor to SPNA density changes and deep-water formation. In the partially coupled case, negative ocean heat advection feedback and the CO 2 -induced subtropical Atlantic saline anomalies imported into the SPNA play a stabilizing role. The results highlight the importance of SPNA salinity gradients and gyre circulation strength in determining the AMOC’s response rate or recovery.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Understanding surfaces and interfaces in nanocomposites of silicone and barium titanate through experiments and modeling

Barium titanate (BTO) is a ferroelectric perovskite used in electronics and energy storage systems because of its high dielectric constant. Decreasing the BTO particle size was shown to increase the dielectric constant of the perovskite, which is an intriguing but contested result. We investigated this result by fabricating silicone-matrix nanocomposite specimens containing BTO particles of decreasing diameter. Furthermore, density functional theory modeling was used to understand the interactions at the BTO particle surface. Combining results from experiments and modeling indicated that polymer type, particle surface interactions, and particle surface structure can influence the dielectric properties of polymer-matrix nanocomposites containing BTO.

composite↗