Search NASA⌕ Search

SEARCH · Search NASA

Results for “Heterogeneous interface”

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 685 records · Page 38

High‐Resolution National‐Scale Water Modeling Is Enhanced by Multiscale Differentiable Physics‐Informed Machine Learning

Abstract The National Water Model (NWM) is a key tool for flood forecasting, planning, and water management. Key challenges facing the NWM include calibration and parameter regionalization when confronted with big data. We present two novel versions of high‐resolution (∼37 km 2 ) differentiable models (a type of hybrid model): one with implicit, unit‐hydrograph‐style routing and another with explicit Muskingum‐Cunge routing in the river network. The former predicts streamflow at basin outlets whereas the latter presents a discretized product that seamlessly covers rivers in the conterminous United States (CONUS). Both versions use neural networks to provide a multiscale parameterization and process‐based equations to provide a structural backbone, which were trained simultaneously (“end‐to‐end”) on 2,807 basins across the CONUS and evaluated on 4,997 basins. Both versions show great potential to elevate future NWM performance for extensively calibrated as well as ungauged sites: the median daily Nash‐Sutcliffe efficiency of all 4,997 basins is improved to around 0.68 from 0.48 of NWM3.0. As they resolve spatial heterogeneity, both versions greatly improved simulations in the western CONUS and also in the Prairie Pothole Region, a long‐standing modeling challenge. The Muskingum‐Cunge version further improved performance for basins >10,000 km 2 . Overall, our results show how neural‐network‐based parameterizations can improve NWM performance for providing operational flood predictions while maintaining interpretability and multivariate outputs. The modeling system supports the Basic Model Interface (BMI), which allows seamless integration with the next‐generation NWM. We also provide a CONUS‐scale hydrologic data set for further evaluation and use.

Song, Yalan [Civil and Environmental Engineering T↗

Dispatch Manager for NEML2 Constitutive Model Calculations Embedded in MOOSE

This report describes the extended capabilities of the NEML2 constitutive modeling library, including a flexible and efficient work dispatching system designed to leverage both CPU and GPU resources. This enhancement addresses one of the primary computational challenges in large-scale simulations: the ability to distribute and execute batches of material model evaluations across heterogeneous computing devices. The new dispatch system introduces a modular set of dispatcher and scheduler classes that coordinate the flow of data and execution between devices. The dispatcher is responsible for efficiently packaging work, managing device-specific memory operations, and synchronizing results. This modularity allows for extensibility, making it straightforward to integrate additional computing backends in the future. From an implementation standpoint, the dispatcher system interfaces seamlessly with NEML2's existing models. They handle device-aware tensor operations, optimize memory transfers, and support asynchronous execution when applicable. This design ensures that batches of material points can be evaluated concurrently, substantially improving throughput compared to previous single-device or serial implementations. These improvements not only enhance the raw performance of NEML2 but also improve its usability in multiscale and high-fidelity simulations, where the simultaneous evaluation of large material point batches is critical. Benchmarks included in the report demonstrate the system’s scalability, highlighting its effectiveness when leveraging modern GPU architectures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Earth-Science Data Co-Locating Tool

This software is used to locate Earth-science satellite data and climate-model analysis outputs in space and time. This enables the direct comparison of any set of data with different spatial and temporal resolutions. It is written in three separate modules that are clearly separated for their functionality and interface with other modules. This enables a fast development of supporting any new data set. In this updated version of the tool, several new front ends are developed for new products. This software finds co-locatable data pairs for given sets of data products and creates new data products that share the same spatial and temporal coordinates. This facilitates the direct comparison between the two heterogeneous datasets and the comprehensive and synergistic use of the datasets.

Lee, Seungwon↗

Synergistic Solvent-Surface Interactions Enable Alkyne Semihydrogenation at Palladium

Enabling higher yield and better selectivity for fine-chemical synthesis through heterogeneous catalysis is intricately linked to the interplay of active sites, reaction conditions, and mass transfer influence provided by the catalyst. Alkyne semihydrogenation is ubiquitous in the production of bulk chemicals in the pharmaceutical, polymer, or fine-chemical industries, but product selectivity remains a major challenge. Here, in this study, we demonstrate that the design of catalysts encompassing nickel (Ni) foams as contiguous monolith supports, decorated with ultralow loading of Pd/PdO x nanoparticles on a carbonized polydopamine interface and tuned with a thin layer of Al 2 O 3 , in conjunction with an optimized reaction environment leads to highly selective alkyne semihydrogenation. The reactions demonstrate good functional group tolerance and applicability to flow reactor systems. Combined computational and experimental studies are presented to describe the synergistic effect between the solvent-surface interaction and the degree of Pd surface reduction that are necessary to promote this selectivity. The system highlights the opportunity for catalyst-solvent codesign as a benign alternative to more complex reactants featuring extrinsic poisons or less-favored dopants.

atomic layer deposition↗

Linear Reciprocating Tribometer for In Situ Neutron Reflectometry of Soft Matter

Neutron reflectometry is a technique for measuring structure near planar interfaces that has been previously used to non-destructively characterize the polymer density of hydrated, dilute, and soft materials. Previous investigations have conducted neutron reflectometry measurements of liquids, gels, emulsion, and polymer solutions at rest, in compression, and subject to shear stress. However, correlating structure with tribological properties of soft materials presents significant experimental challenges for prior instruments due to wall slip, sample thickness, and structural heterogeneity (e.g., depth-wise gradients). A linear reciprocating tribometer offers several advantages for in situ neutron reflectometry studies, including uniform velocity profiles, constant shear stress over large regions of interest, and independent control of normal force and sliding velocity during measurements. This work outlines basic considerations for the design of a custom linear reciprocating tribometer that operates in a neutron beamline and includes commissioning measurements. The tribometer is designed to compress soft and hydrated materials against linearly reciprocating silicon disks. The three key design considerations for this tribometer are (1) safety, (2) neutron transmission, and (3) sample positioning. This instrument design will enable in situ studies of soft matter and illuminate the role of interfacial structure on tribological phenomena.

42 ENGINEERING↗

Evolution of Local Microstructures (ELMS): Spatial Instabilities of Coarsening

This work examines the diffusional growth of discrete phase particles dispersed within a matrix. Engineering materials are microstructurally heterogeneous, and the details of the microstructure determine how well that material performs in a given application. Critical to the development of designing multiphase microstructures with long-term stability is the process of Ostwald ripening. Ripening, or phase coarsening, is a diffusion-limited process which arises in polydisperse multiphase materials. Growth and dissolution occur because fluxes of solute, driven by chemical potential gradients at the interfaces of the dispersed phase material, depend on particle size. The kinetics of these processes are "competitive," dictating that larger particles grow at the expense of smaller ones, overall leading to an increase of the average particle size. The classical treatment of phase coarsening was done by Todes, Lifshitz, and Slyozov, (TLS) in the limit of zero volume fraction, V(sub v), of the dispersed phase. Since the publication of TLS theory there have been numerous investigations, many of which sought to describe the kinetic scaling behavior over a range of volume fractions. Some studies in the literature report that the relative increase in coarsening rate at low (but not zero) volume fractions compared to that / 2 1/ 3 predicted by TLS is proportional to V(sub v)(exp 1/2), whereas others suggest V(sub v)(exp 1/3). This issue has been resolved recently by simulation studies at low volume fractions in three dimensions by members of the Rensselaer/MSFC team.

Glicksman, Martin E.↗

Evolution of Local Microstructures: Spatial Instabilities of Coarsening Clusters

This work examines the diffusional growth of discrete phase particles dispersed within a matrix. Engineering materials are microstructurally heterogeneous, and the details of the microstructure determine how well that material performs in a given application. Critical to the development of designing multiphase microstructures with long-term stability is the process of Ostwald ripening. Ripening, or phase coarsening, is a diffusion-limited process which arises in polydisperse multiphase materials. Growth and dissolution occur because fluxes of solute, driven by chemical potential gradients at the interfaces of the dispersed phase material, depend on particle size. The kinetics of these processes are "competitive," dictating that larger particles grow at the expense of smaller ones, overall leading to an increase of the average particle size. The classical treatment of phase coarsening was done by Todes, Lifshitz, and Slyozov, (TLS) in the limit of zero volume fraction, V(sub v), of the dispersed phase. Since the publication of TLS theory there have been numerous investigations, many of which sought to describe the kinetic scaling behavior over a 0 range of volume fractions. Some studies in the literature report that the relative increase in coarsening rate at low (but not zero) volume fractions compared to that predicted by TLS is proportional to v(sub v)(exp 1/2), whereas others suggcest V(sub v)(exp 1/3). This issue has been resolved recently by simulation studies at low volume fractions in three dimensions by members of the Rensselaer/MSFC team. Our studies of ripening behavior using large-scale numerical simulations suggest that although there are different circumstances which can lead to either scaling law, the most important length scale at low volume fractions is the diffusional analog of the Debye screening length. The numerical simulations we employed exploit the use of a recently developed "snapshot" technique, and identifies the nature of the coarsening dynamics at various volume fractions. Preliminary results of numerical and experimental investigations, focused on the growth of finite particle clusters, provide important insight into the nature of the transition between the two scaling regimes. The companion microgravity experiment centers on the growth within finite particle clusters, and follows the temporal dynamics driving microstructural evolution, using holography.

Frazier, Donald O.↗

Physics-Based Methods of Failure Analysis and Diagnostics in Human Space Flight

The Integrated Health Management (IHM) for the future aerospace systems requires to interface models of multiple subsystems in an efficient and accurate information environment at the earlier stages of system design. The complexity of modern aeronautic and aircraft systems (including e.g. the power distribution, flight control, solid and liquid motors) dictates employment of hybrid models and high-level reasoners for analysing mixed continuous and discrete information flow involving multiple modes of operation in uncertain environments, unknown state variables, heterogeneous software and hardware components. To provide the information link between key design/performance parameters and high-level reasoners we rely on development of multi-physics performance models, distributed sensors networks, and fault diagnostic and prognostic (FD&P) technologies in close collaboration with system designers. The main challenges of our research are related to the in-flight assessment of the structural stability, engine performance, and trajectory control. The main goal is to develop an intelligent IHM that not only enhances components and system reliability, but also provides a post-flight feedback helping to optimize design of the next generation of aerospace systems. Our efforts are concentrated on several directions of the research. One of the key components of our strategy is an innovative approach to the diagnostics/prognostics based on the real time dynamical inference (DI) technologies extended to encompass hybrid systems with hidden state trajectories. The major investments are into the multiphysics performance modelling that provides an access of the FD&P technologies to the main performance parameters of e.g. solid and liquid rocket motors and composite materials of the nozzle and case. Some of the recent results of our research are discussed in this chapter. We begin by introducing the problem of dynamical inference of stochastic nonlinear models and reviewing earlier results. Next, we present our analytical approach to the solution of this problem based on the path integral formulation. The resulting algorithm does not require an extensive global search for the model parameters, provides optimal compensation for the effects of dynamical noise, and is robust for a broad range of dynamical models. In the following Section the strengths of the algorithm are illustrated illustrated by inferring the parameters of the stochastic Lorenz system and comparing the results with those of earlier research. Next, we discuss a number of recent results in application to the development of the IHM for aerospace system. Firstly, we apply dynamical inference approach to a solution of classical three tank problems with mixed unknown continuous and binary parameters. The problem is considered in the context of ground support system for filling fuel tanks of liquid rocket motors. It is shown that the DI algorithm is well suited for successful solution of a hybrid version of this benchmark problem even in the presence of additional periodic and stochastic perturbation of unknown strength. Secondly, we illustrate our approach by its application to an analysis of the nozzle fault in a solid rocket motor (SRM). The internal ballistics of the SRM is modelled as a set of one-dimensional partial differential equations coupled to the dynamics of the propellant regression. In this example we are specifically focussed on the inference of discrete and continuous parameters of the nozzle blocking fault and on the possibility of an application of the DI algorithm to reducing the probability of "misses" of an on-board FD&P for SRM. In the next section re-contact problem caused by first stage/upper stage separation failure is discussed. The reaction forces imposed on the nozzle of the upper stage during the re-contact and their connection to the nozzle damage and to the thrust vector control (TVC) signal are obtained. It is shown that transient impact induced torquean be modelled as a response of an effective damped oscillator. A possible application of the DI algorithm to the inference of damage parameters and predicting fault dynamics ahead of time using the actuator signal is discussed. Finally, we formulate Bayesian inferential framework for development of the IHM system for in-flight structural health monitoring (SHM) of composite materials. We consider the signal generated by piezoelectric actuator mounted on composite structure generating elastic waves in it. The signal received by the sensor is than compared with the baseline signal. The possibility of damage inference is discussed in the context of development of the SHM.

Smelyanskiy, Vadim N.↗

Microstructural heterogeneity of the buffer layer of TRISO nuclear fuel particles

Tristructural isotropic (TRISO) nuclear fuel particles contain a layered spherical shell designed to retain fission products; however, failure occurs in rare cases—commonly initiated in the porous pyrocarbon buffer layer. Achieving a comprehensive understanding of the buffer-initiated failure mechanisms requires detailed characterization of the buffer porosity and its heterogeneous distribution across multiple length scales. Here we performed FIB-SEM tomography across the buffer layer thickness (~100 µm) to produce 3D reconstructions of the buffer microstructure with 50 nm spatial resolution. We found an average overall porosity of ~14%, which does not solely account for the low density of the buffer (50% of the theoretical density). Additionally, the local porosity and its fluctuation increase from the kernel interface towards the inner pyrocarbon (IPyC) layer, which we attribute to the chemical vapor deposition process conditions during the TRISO particle fabrication. Detailed characterization of the porous microstructure—including analysis of the pore size, distribution, shape, and orientation—provides insight into the process-structure-property-performance relations of TRISO nuclear fuel particles and will inform multiscale models designed to predict the failure of TRISO particles under irradiation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Comprehensive Stress Drop Map From Trench to Depth in the Northern Chilean Subduction Zone

We compute stress drops for earthquakes in Northern Chile recorded between 2007 and 2021. By applying two analysis techniques, (a) the spectral ratio (SR) method and (b) the spectral decomposition (SDC) method, a stress drop map for the subduction zone consisting of 51,510 stress drop values is produced. We build an extended set of empirical Green’s functions (EGF) for the SR method by systematic template matching. Outputs are used to compare with results from the SDC approach, where we apply cell-wise obtained global EGF's to compensate for the structural heterogeneity of the subduction zone. We find a good consistency of results of the two methods. The increased spatial coverage and quantity of stress drop estimates from the SDC method facilitate a consistent stress drop mapping of the different seismotectonic domains. Albeit only small differences of median stress drop, strike-perpendicular depth sections clearly reveal systematic variations, with earthquakes at different seismotectonic locations exhibiting distinct values. In particular, interface seismicity is characterized by the lowest observed median value, whereas upper plate earthquakes show noticeably higher stress drop values. Intermediate depth earthquakes show comparatively high average stress drop and a rather strong depth-dependent increase of median stress drop. Additionally, we observe spatio-temporal variability of stress drops related to the occurrence of the two megathrust earthquakes in the study region. The presented study is the first coherent large scale 3D stress drop mapping of the Northern Chilean subduction zone. It provides an important component for further detailed analysis of the physics of earthquake ruptures.

58 GEOSCIENCES↗

Probing charge density in materials with atomic resolution in real space

The charge distribution in materials at the nanoscale can often explain the origin of macroscopic properties such as localized conductivity or the plasmonic response and illuminate more fundamental changes in the microscopic structure such as changes in chemical bonding characteristics. Previously, direct visualization of the charge density with high spatial resolution was often a missing link in the formation of structure–property relationships, especially in heterogeneous materials systems. Furthermore, recent advancements in microscopy technology have enabled researchers to visualize the charge distribution in materials down to subatomic length scales. In this Technical Review, we discuss the developments in high-resolution real-space charge distribution imaging using diffraction techniques and electron microscopy, with a focus on the recent advancement of four-dimensional scanning transmission electron microscopy, electron holography, and applications to materials interfaces.

Electronic properties and materials↗

The surface and through crack problems in layered orthotropic plates

An analytical method is developed for a relatively accurate calculation of Stress Intensity Factors in a laminated orthotropic plate containing a through or part-through crack. The laminated plate is assumed to be under bending or membrane loading and the mode 1 problem is considered. First three transverse shear deformation plate theories (Mindlin's displacement based first-order theory, Reissner's stress-based first-order theory, and a simple-higher order theory due to Reddy) are reviewed and examined for homogeneous, laminated and heterogeneous orthotropic plates. Based on a general linear laminated plate theory, a method by which the stress intensity factors can be obtained in orthotropic laminated and heterogeneous plates with a through crack is developed. Examples are given for both symmetrically and unsymmetrically laminated plates and the effects of various material properties on the stress intensity factors are studied. In order to implement the line-spring model which is used later to study the surface crack problem, the corresponding plane elasticity problem of a two-bonded orthotropic plated containing a crack perpendicular to the interface is also considered. Three different crack profiles: an internal crack, an edge crack, and a crack terminating at the interface are considered. The effect of the different material combinations, geometries, and material orthotropy on the stress intensity factors and on the power of stress singularity for a crack terminating at the interface is fully examined. The Line Spring model of Rice and Levy is used for the part-through crack problem. The surface crack is assumed to lie in one of the two-layered laminated orthotropic plates due to the limitation of the available plane strain results. All problems considered are of the mixed boundary value type and are reduced to Cauchy type of singular integral equations which are then solved numerically.

Erdogan, Fazil↗

A carbon isotope mass balance for an anoxic marine sediment: Isotopic signatures of diagenesis

A carbon isotope mass balance was determined for the sediments of Cape Lookout Bight, NC to constrain the carbon budgets published previously. The diffusive, ebullitive and burial fluxes of sigma CO2 and CH4, as well as the carbon isotope signatures of these fluxes, were measured. The flux-weighted isotopic signature of the remineralized carbon (-18.9 plus or minus 2.7 per mil) agreed with the isotopic composition of the remineralized organic carbon determined from the particulate organic carbon (POC) delta(C-13) profiles (-19.2 plus or minus 0.2), verifying the flux and isotopic signature estimates. The measured delta(C-13) values of the sigma CO2 and CH4 diffusive fluxes were significantly different from those calculated from porewater gradients. The differences appear to be influenced by methane oxidation at the sediment-water interface, although other potential processes cannot be excluded. The isotope mass balance provides important information concerning the locations of potential diagenetic isotope effects. Specifically, the absence of downcore change in the delta(C-13) value of the POC fraction and the identical isotopic composition of the POC and the products of remineralization indicate that no isotopic fractionation is expressed during the initial breakdown of the POC, despite its isotopically heterogeneous composition.

Boehme, Susan E.↗

TAMM: Tensor algebra for many-body methods

Tensor algebra operations such as contractions in computational chemistry consume a significant fraction of the computing time on large-scale computing platforms. The widespread use of tensor contractions between large multi-dimensional tensors in describing electronic structure theory has motivated the development of multiple tensor algebra frameworks targeting heterogeneous computing platforms. In this paper, we present Tensor Algebra for Many-body Methods (TAMM), a framework for productive and performance-portable development of scalable computational chemistry methods. TAMM decouples the specification of the computation from the execution of these operations on available high-performance computing systems. With this design choice, the scientific application developers (domain scientists) can focus on the algorithmic requirements using the tensor algebra interface provided by TAMM, whereas high-performance computing developers can direct their attention to various optimizations on the underlying constructs, such as efficient data distribution, optimized scheduling algorithms, and efficient use of intra-node resources (e.g., graphics processing units). The modular structure of TAMM allows it to support different hardware architectures and incorporate new algorithmic advances. We describe the TAMM framework and our approach to the sustainable development of scalable ground- and excited-state electronic structure methods. We present case studies highlighting the ease of use, including the performance and productivity gains compared to other frameworks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Shock state distributions in porous tantalum and characterization with multipoint velocimetry

Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, in this work, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SPADES (Scalable Parallel Discrete Events Simulation) [SWR-24-99]

SPADES (Solver for PArallel Discrete Event Simulation) is an open-source parallel discrete event simulation (PDES) package built on the AMReX library. Targeted at solving discrete event systems in parallel, this software package aims to be performance portable and scalable on heterogeneous computing architectures, e.g., graphic processing units (GPU). SPADES implements optimistic synchronization with rollback through an implementation of the Time Warp algorithm. An alternative conservative synchronization approach is also implemented using the Lower Bound on Incoming Time Stamp. In our implementation, logical processes are represented as cells in a grid and event messages are represented as particles. SPADES supports various parallel decomposition strategies, including the use of the Message Passing Interface (MPI) and OpenMP threading. All major GPU architectures (e.g., Intel, AMD, NVIDIA) are supported through the use of performance portability functionalities implemented in AMReX. The SPADES software is released in NREL Software Record SWR-24-99 “SPADES (Scalable Parallel Discrete Events Simulation)”.

Henry de Frahan, Marc [National Renewable Energy L↗

Cryogenic Carbon Monoxide Oxidation on Cuprous Oxide

Performing oxidation reactions at low temperatures using earth-abundant materials is crucial for advancing solutions for sustainable chemistry. CO oxidation serves as a benchmark reaction to characterize oxidation and to advance fundamental concepts in surface chemistry. While there are several examples of CO oxidation occurring on metal oxides at low temperatures, from 300 K to ∼200 K, reactivity in the cryogenic temperature regime typically requires a metal nanoparticle on a metal oxide. Here, we show oxygen atoms on the (111) facet of Cu 2 O react with CO to form CO 2 at temperatures below 100 K. Combining spectroscopic experimental evidence with calculations, we propose a low barrier path for CO oxidation at reconstructed surface sites on Cu 2 O(111). This finding is a rare example of an earth-abundant metal oxide, in this case copper, that can provide highly reactive multifunctional sites, enabling both adsorption and reaction fundamental steps toward the efficient heterogeneous oxidation of chemicals.

36 MATERIALS SCIENCE↗

High-rate electroreduction of carbon monoxide to multi-carbon products

Carbon monoxide electrolysis has previously been reported to yield enhanced multi-carbon (C 2+ ) Faradaic efficiencies of up to ~55%, but only at low reaction rates. This is due to the low solubility of CO in aqueous electrolytes and operation in batch-type reactors. Here, we present a high-performance CO flow electrolyser with a well controlled electrode–electrolyte interface that can reach total current densities of up to 1 A cm –2 , together with improved C 2+ selectivities. Computational transport modelling and isotopic C 18 O reduction experiments suggest that the enhanced activity is due to a higher surface pH under CO reduction conditions, which facilitates the production of acetate. Here, at optimal operating conditions, we achieve a C 2+ Faradaic efficiency of ~91% with a C 2+ partial current density over 630 mA cm –2 . Further investigations show that maintaining an efficient triple-phase boundary at the electrode–electrolyte interface is the most critical challenge in achieving a stable CO/CO 2 electrolysis process at high rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗