SEARCH · Engineering Papers

Results for “variational multiscale”

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 145 records · Page 8

All sources NASA NTRS DOE OSTI

DOE OSTI2021

The persistence of large scale structures. Part I. Primordial non-Gaussianity

Abstract We develop an analysis pipeline for characterizing the topology of large scale structure and extracting cosmological constraints based on persistent homology . Persistent homology is a technique from topological data analysis that quantifies the multiscale topology of a data set, in our context unifying the contributions of clusters, filament loops, and cosmic voids to cosmological constraints. We describe how this method captures the imprint of primordial local non-Gaussianity on the late-time distribution of dark matter halos, using a set of N-body simulations as a proxy for real data analysis. For our best single statistic, running the pipeline on several cubic volumes of size 40 (Gpc/h) 3 , we detect f NL loc =10 at 97.5% confidence on ~ 85% of the volumes. Additionally wetest our ability to resolve degeneracies betweenthe topological signature of f NL loc and variation of σ 8 and argue that correctly identifying nonzero f NL loc in this case is possible via an optimal template method. Our method relies on information living at $\mathcal{O}$(10) Mpc/h, a complementary scale with respect to commonly used methods such as the scale-dependent bias in the halo/galaxy power spectrum. Therefore, while still requiring a large volume, our method does not require sampling long-wavelength modes to constrain primordial non-Gaussianity. Moreover, our statistics are interpretable: we are able to reproduce previous results in certain limits and we make new predictions for unexplored observables, such as filament loops formed by dark matter halos in a simulation box.

Astronomy & Astrophysics↗

DOE OSTI2025

Multiphase computational fluid dynamics modeling of reacting flows in absorption columns for carbon capture

First-principles derived computational fluid dynamics (CFD) simulations have been proposed as a fundamental tool for investigating solvent-based CO 2 absorption in packed columns due to their ability to accurately represent the underlying nonlinear, multiscale dynamics. Numerous studies have previously utilized such CFD simulations to investigate hydrodynamics of columns with structured and random packings by assessing the key hydrodynamic metrics such as the interfacial and wetted areas. While mapping such metrics for different conditions is essential to the optimization of absorption columns, it is not sufficient, as the CO 2 capture rate depends also on the coupled, nonlinear dynamics from the underlying chemical reaction kinetics, thermodynamics, and heat-transfer rates. In this work, we present detailed CFD simulation results augmented by incorporating the effects of interfacial physical mass transfer of CO 2 , heat release from chemical reaction kinetics, and thermophysical property variations from resulting temperature gradients. We demonstrate the applicability of the proposed approach in numerically assessing the performance of packed columns by evaluating key hydrodynamic quantities, CO 2 absorption rates, and temperature rise in a reference column with packings that are structurally similar to the Sulzer Mellapak™ 250.Y packing, for different solvent inflow velocities and temperatures. Predictions from simulation results are found to be consistent with the trends in experimental observations from the literature, suggesting that the predictive capabilities of the simulation framework can be leveraged to guide the future development of absorber-column designs and optimized process flowsheets.

Absorption columns↗

DOE OSTI2022

Insight into ideal shear strength of Ni-based dilute alloys using first-principles calculations and correlational analysis

Here the present work examines the effect of alloying elements (denoted X) on the ideal shear strength for 26 dilute Ni-based alloys, Ni11X, as determined by first-principles calculations of pure alias shear deformations. The variations in ideal shear strength are quantitatively explored with correlational analysis techniques, showing the importance of atomic properties such as size and electronegativity. The shear moduli of the alloys are affirmed to show a strong linear relationship with their ideal shear strengths, while the shear moduli of the individual alloying elements were not indicative of alloy shear strength. Through combination with available ideal shear strength data on Mg alloys, a potential application of the Ni alloy data is demonstrated in the search for a set of atomic features suitable for machine learning applications to mechanical properties. As another illustration, the calculated Ni ideal shear strengths play a key role in a predictive multiscale framework for deformation behavior of single crystal alloys at large strains, as shown by simulated stress–strain curves.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Multiscale Modeling Framework for Lithium Nucleation in 3D Porous Carbon Anodes

Porous carbon scaffolds offer a promising route for mitigating non-uniform lithium (Li) plating to enhance the safety and longevity of Li metal batteries. However, the influence of microstructural morphology on Li nucleation is not well understood. Here, we present a multiscale modeling framework to investigate how the porous microstructure of carbon materials affects Li nucleation behavior. Ab initio molecular dynamics simulations quantify the nucleation energy barriers of Li on graphene as a function of Li content, surface curvature, and applied potential, providing key parameters for a classical nucleation theory (CNT) model. From macroscale half-cell simulations, we obtained Li concentration and electrical potential profiles to define boundary conditions for mesoscopic simulations. At the mesoscale, three distinct synthetic 3D microstructures with different porosities and characteristic feature sizes are generated to resolve local distributions of Li flux, current density, and mechanical stress. These outputs are integrated into the CNT model to map spatial variation in nucleation rates. Our findings reveal trade-offs between suppressing nucleation rates and achieving spatial uniformity, offering design guidelines for optimizing porous carbon anodes to balance nucleation control and mechanical integrity.

Materials science↗

DOE OSTI2024

Peridynamic Models for Random Media Found by Coarse Graining

Using coarse graining, the upscaled mechanical properties of a solid with small scale heterogeneities are derived. The method maps internal forces at the small scale onto peridynamic bond forces in the coarse grained mesh. These upscaled bond forces are used to calibrate a peridynamic material model with position-dependent parameters. These parameters incorporate mesoscale variations in the statistics of the small scale system. The upscaled peridynamic model can have a much coarser discretization than the original small scale model, allowing larger scale simulations to be performed efficiently. The convergence properties of the method are investigated for representative random microstructures. In conclusion, a bond breakage criterion for the upscaled peridynamic material model is also demonstrated.

36 MATERIALS SCIENCE↗

NASA NTRS2016

Center of Mass Estimation for a Spinning Spacecraft Using Doppler Shift of the GPS Carrier Frequency

A sequential filter is presented for estimating the center of mass (CM) of a spinning spacecraft using Doppler shift data from a set of onboard Global Positioning System (GPS) receivers. The advantage of the proposed method is that it is passive and can be run continuously in the background without using commanded thruster firings to excite spacecraft dynamical motion for observability. The NASA Magnetospheric Multiscale (MMS) mission is used as a test case for the CM estimator. The four MMS spacecraft carry star cameras for accurate attitude and spin rate estimation. The angle between the spacecraft nominal spin axis (for MMS this is the geometric body Z-axis) and the major principal axis of inertia is called the coning angle. The transverse components of the estimated rate provide a direct measure of the coning angle. The coning angle has been seen to shift slightly after every orbit and attitude maneuver. This change is attributed to a small asymmetry in the fuel distribution that changes with each burn. This paper shows a correlation between the apparent mass asymmetry deduced from the variations in the coning angle and the CM estimates made using the GPS Doppler data. The consistency between the changes in the coning angle and the CM provides validation of the proposed GPS Doppler method for estimation of the CM on spinning spacecraft.

Spacecraft↗

NASA NTRS2016

Center of Mass Estimation for a Spinning Spacecraft Using Doppler Shift of the GPS Carrier Frequency

Estimation↗

DOE OSTI2022

Mesoscale structural gradients in human tooth enamel

The outstanding mechanical and chemical properties of dental enamel emerge from its complex hierarchical architecture. An accurate, detailed multiscale model of the structure and composition of enamel is important for understanding lesion formation in tooth decay (dental caries), enamel development (amelogenesis) and associated pathologies (e.g., amelogenesis imperfecta or molar hypomineralization), and minimally invasive dentistry. Although features at length scales smaller than 100 nm (individual crystallites) and greater than 50 µm (multiple rods) are well understood, competing field of view and sampling considerations have hindered exploration of mesoscale features, i.e., at the level of single enamel rods and the interrod enamel (1 to 10 µm). Here, we combine synchrotron X-ray diffraction at submicrometer resolution, analysis of crystallite orientation distribution, and unsupervised machine learning to show that crystallographic parameters differ between rod head and rod tail/interrod enamel. This variation strongly suggests that crystallites in different microarchitectural domains also differ in their composition. Thus, we use a dilute linear model to predict the concentrations of minority ions in hydroxylapatite (Mg 2+ and CO 3 2- /Na + ) that plausibly explain the observed lattice parameter variations. While differences within samples are highly significant and of similar magnitude, absolute values and the sign of the effect for some crystallographic parameters show interindividual variation that warrants further investigation. In conclusion, by revealing additional complexity at the rod/interrod level of human enamel and leaving open the possibility of modulation across larger length scales, these results inform future investigations into mechanisms governing amelogenesis and introduce another feature to consider when modeling the mechanical and chemical performance of enamel.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2025

Quantitative approaches for multiscale structural analysis with atomic resolution electron microscopy

Atomic-resolution imaging with scanning transmission electron microscopy is a powerful tool for characterizing the nanoscale structure of materials, in particular features such as defects, local strains, and symmetry-breaking distortions. In addition to advanced instrumentation, the effectiveness of the technique depends on computational image analysis to extract meaningful features from complex datasets recorded in experiments, which can be complicated by the presence of noise and artifacts, small or overlapping features, and the need to scale analysis over large representative areas. Here, we present image analysis approaches which synergize real and reciprocal space information to efficiently and reliably obtain meaningful structural information with picometer scale precision across hundreds of nanometers of material from atomic-resolution electron microscope images. Damping superstructure peaks in reciprocal space allows symmetry-breaking structural distortions to be disentangled from other sources of inhomogeneity and measured with high precision. Real-space fitting of the wavelike signals resulting from Fourier filtering enables absolute quantification of lattice parameter variations and strain, as well as the uncertainty associated with these measurements. Implementations of these algorithms are made available as an open source python package.

36 MATERIALS SCIENCE↗

DOE OSTI2022

Multiscale characterization and comparison of historical and modern nuclear graphite grades

Beginning with Chicago Pile I, graphite has been used as a moderator material in nuclear power stations and is considered a potential material for use in future Generation IV advanced reactors. The microstructure of graphite is responsible for much of its mechanical and thermo-physical properties, and how it responds to irradiation. To understand graphite microstructure, it is necessary to understand its porosity at the macro- and micro-scales; and to understand its porosity, it is necessary to characterize the morphological connectivity of the void content and the two main phases of graphite: filler and binder. Here, using several microscopy and analytical techniques, a detailed examination of the heterogeneity, microstructure and pore structure of different graphite grades and their binder and filler phases is presented. Significant differences were found between coarser and finer nuclear grades. Coarse grades have a more diverse range of filler particles, pores and thermal cracks. Finer grades have a more well-defined pore size distribution, fewer variations of filler particles sizes and do not contain as many large thermal cracks. Fine grades tend to have a well-connected network of pores whereas coarser grades contain a larger content of closed porosity. The framework developed within this work can be applied and used to assess the various graphite grades that would down-select materials for specific use in graphite moderated reactor designs.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

NASA NTRSDate not supplied

Multiscale Modeling of Fracture Strength in Fibrous Thermal Protection System Materials

This work presents a multiscale modeling approach to predict the fracture strength of fibrous Thermal Protection System (TPS) materials. The model assumes that system failure is initiated at the joints between individual fibers. We investigated three distinct TPS compositions: amorphous silica, alumina and aluminosilicate fibers. Molecular dynamics (MD) simulations were employed to determine the fracture strength values of these fiber joints for both material systems. These fracture strength values were then integrated into simulations of 3D randomly populated fiber structures, where tensile load transfer occurs through the fiber joints. These microscale properties are upscaled through a renormalization approach [1] to predict macroscale tensile strength of 3D random fiber networks, accounting for joint-dominated failure and effective load-bearing area. The study concludes by demonstrating the resulting strength variation as a function of material composition, fiber density, and morphology. We also show validation of results by comparing them against explicit fiber finite element (FE) modeling [2] where fiber joint fracture is represented by cohesive elements.

Jaehyun Cho↗

DOE OSTI2021

Towards a Unified theory of Fractional and Nonlocal Vector Calculus

Nonlocal and fractional-order models capture effects that classical partial differential equations cannot describe; for this reason, they are suitable for a broad class of engineering and scientific applications that feature multiscale or anomalous behavior. This has driven a desire for a vector calculus that includes nonlocal and fractional gradient, divergence and Laplacian type operators, as well as tools such as Green’s identities, to model subsurface transport, turbulence, and conservation laws. In the literature, several independent definitions and theories of nonlocal and fractional vector calculus have been put forward. Some have been studied rigorously and in depth, while others have been introduced ad-hoc for specific applications. The goal of this work is to provide foundations for a unified vector calculus by (1) consolidating fractional vector calculus as a special case of nonlocal vector calculus, (2) relating unweighted and weighted Laplacian operators by introducing an equivalence kernel, and (3) proving a form of Green’s identity to unify the corresponding variational frameworks for the resulting nonlocal volume-constrained problems. Here, the proposed framework goes beyond the analysis of nonlocal equations by supporting new model discovery, establishing theory and interpretation for a broad class of operators, and providing useful analogues of standard tools from the classical vector calculus.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

A Unified Theory of Fractional Nonlocal and Weighted Nonlocal Vector Calculus

Nonlocal and fractional-order models capture effects that classical partial differential equations cannot describe; for this reason, they are suitable for a broad class of engineering and scientific applications that feature multiscale or anomalous behavior. This has driven a desire for a vector calculus that includes nonlocal and fractional gradient, divergence and Laplacian type operators, as well as tools such as Green's identities, to model subsurface transport, turbulence, and conservation laws. In the literature, several independent definitions and theories of nonlocal and fractional vector calculus have been put forward. Some have been studied rigorously and in depth, while others have been introduced ad-hoc for specific applications. The goal of this work is to provide foundations for a unified vector calculus by (1) consolidating fractional vector calculus as a special case of nonlocal vector calculus, (2) relating unweighted and weighted Laplacian operators by introducing an equivalence kernel, and (3) proving a form of Green's identity to unify the corresponding variational frameworks for the resulting nonlocal volume-constrained problems. The proposed framework goes beyond the analysis of nonlocal equations by supporting new model discovery, establishing theory and interpretation for a broad class of operators, and providing useful analogues of standard tools from the classical vector calculus.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Deep learning-based spatio-temporal estimate of greenhouse gas emissions using satellite data

Accurate estimation of greenhouse gases (GHGs) emissions is very important for developing mitigation strategies to climate change by controlling and reducing GHG emissions. This project aims to develop multiple deep learning approaches to estimate anthropogenic greenhouse gas emissions using multiple types of satellite data. NO2 concentration is chosen as an example of GHGs to evaluate the proposed approach. Two sentinel satellites (sentinel-2 and sentinel-5P) provide multiscale observations of GHGs from 10-60m resolution (sentinel-2) to ~kilometer scale resolution (sentinel-5P). Among multiple deep learning (DL) architectures evaluated, two best DL models demonstrate that key features of spatio-temporal satellite data and additional information (e.g., observation times and/or coordinates of ground stations) can be extracted using convolutional neural networks and feed forward neural networks, respectively. In particular, irregular time series data from different NO 2 observation stations limit the flexibility of long short-term memory architecture, requiring zero-padding to fill in missing data. However, deep neural operator (DNO) architecture can stack time-series data as input, providing the flexibility of input structure without zero-padding. As a result, the DNO outperformed other deep learning architectures to account for time-varying features. Overall, temporal patterns with smooth seasonal variations were predicted very well, while frequent fluctuation patterns were not predicted well. In addition, uncertainty quantification using conformal inference method is performed to account for prediction ranges. Overall, this research will lead to a new groundwork for estimating greenhouse gas concentrations using multiple satellite data to enhance our capability of tracking the cause of climate change and developing mitigation strategies.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2026

Multimodal Nanoscale Mapping of Local Structure and CO 2 Adsorption in Metal–Organic Frameworks

Diamine functionalization of the metal−organic framework Mg 2 (dobpdc) (dobpdc 4− = 4,4′-dioxidobiphenyl-3,3′-dicarboxylate) significantly enhances its selectivity for CO 2 capture from flue gases and air. The structure and CO 2 capacity of such materials are typically assessed using bulk techniques that rely on averaging signal over large ensembles of unit cells, obscuring local heterogeneities, such as variations in CO 2 occupancy across individual nanocrystals. To resolve this limitation, we demonstrate a multimodal, nanoscale characterization of Mg 2 (dobpdc) appended with 1,3-diaminopropane. By employing recently developed characterization techniques at progressively smaller length scales, we uncover insights from correspondingly smaller populations of unit cells. First, we use parallel-beam 3D electron diffraction (3D ED) to identify a prominent expansion in lattice parameters upon desorption of CO 2 , as observed at the level of single nanocrystals. Second, we use convergent-probe 4D scanning transmission electron microscopy (4D-STEM) to quantify associated differences in lattice strain as a function of gas loading and diamine appending. These measurements sample small subvolumes within individual nanocrystals. Finally, we apply infrared scattering scanning near-field optical microscopy (IR s- SNOM) to confirm variable CO 2 chemisorption across adsorption sites at the surface of single nanocrystals. This multimodal, multiscale approach allows us to map heterogeneity within individual nanocrystals. Collectively, these findings emphasize the importance of local, nanoscale characterization of metal−organic frameworks in revealing previously unresolvable features that impact their performance.

Karstens, Sarah L. [University of California, Berk↗

DOE OSTI2025

Laser Powder Directed Energy Deposition of Steels for Nuclear Applications

This comprehensive investigation examines the structure–property relationships in two nuclear alloy systems—Alloy 709 (A709) austenitic stainless steel and Grade 92 (G-92) ferritic/martensitic (F/M) steel—manufactured via directed energy deposition (DED) for sodium-cooled fast reactor applications. This study establishes the fundamental mechanisms for controlling microstructures for optimizing the mechanical performance of additively manufactured nuclear materials through systematic heat treatment optimization and multiscale characterization. As-deposited A709 steel develops a complex multiscale strengthening architecture consisting of a fine cellular solidification structure with diameter of 2-3 µm within10–50 µm grains, elevated dislocation densities from rapid thermal cycling, and grain boundary precipitates that activate concurrent Hall–Petch, dislocation, and precipitation hardening mechanisms to achieve exceptional properties [yield strength (YS): 603 MPa, ultimate tensile strength (UTS): 844 MPa, Vickers hardness: 220 HV] that achieve a 44% superior strength compared to that of the wrought material. Heat treatments produce different results. Solution annealing (SA) dissolves the cellular structure and reduces the hardness to 190 HV. Precipitation treatment (PT) keeps the cellular structure but adds carbides, allowing the hardness to reach 205 HV. The best approach combines both treatments (SA+PT) and creates uniform precipitate distributions with M 23 C 6 carbides at the grain boundaries and MX carbonitrides in the matrix, achieving a hardness of 195 HV. However, directional differences persist, with a 12%–15% strength variation between orientations due to the inherited layered microstructural architecture that survives aggressive heat treatment. While tensile testing at 550°C demonstrates 40%–50% thermal softening with dynamic strain aging, DED A709 steel still maintains a 71% higher YS than that of the wrought material. Ion irradiation studies (100–400 dpa) of DED A709 steel reveal progressive radiation damage with increasing void density and radiation-induced segregation causing nickel enrichment and chromium depletion, which will ultimately compromise mechanical properties. As-deposited G-92 exhibits exceptional strength (UTS: 1650–1700 MPa, 430 HV) through a complex microstructure containing both ferrite and martensite phases, a high geometrically necessary dislocation (GND) density (17.04×10 14 /m 2 ), and fine carbides. Heat treatments create distinct changes. Normalizing produces fresh martensite with the highest hardness (460 HV) and an increased GND density (20.23×10 14 /m 2 ). Tempering develops dual precipitation systems and reduces the hardness to 290 HV. The optimal approach uses sequential normalizing plus tempering, achieving balanced properties with the lowest hardness (250 HV) and a reduced GND density (11.01×10 14 /m 2 ). A processing-dependent anisotropy is observed: horizontal specimens achieve superior ductile behavior, while vertical specimens exhibit brittle failure. A tempering heat treatment successfully mitigates this anisotropic behavior by transforming the hard martensitic as-deposited structure into tempered martensite enabling both horizontal and vertical specimens to exhibit similar stress–strain characteristics with visible necking behavior. Remarkably, testing at 550°C reveals a reversal in the anisotropy, where as-deposited specimens achieve near isotropy with superior thermal stability (a 15%–20% strength reduction), while tempered specimens develop an orientation dependence with a 25%–30% strength reduction. Both alloy systems demonstrate that DED processing creates specimens with a superior strength through refined microstructural features, though with distinct strengthening mechanisms—austenitic through cellular structures and precipitates versus F/M through phase transformations and precipitates. Heat treatment optimization requires alloy-specific approaches, with A709 benefiting from controlled precipitation while G-92 requires careful phase transformation control. The results show that DED manufacturing can produce nuclear materials with exceptional performance, but directional effects and temperature-dependent behavior must be carefully considered for reactor component design and qualification.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS2009

The Magnetospheric Multiscale Missions Fast Plasma Investigations Dual Electron Spectrometer Development

The Magnetospheric Multiscale mission (MMS) is designed to examine magnetic reconnection that occurs on both the Earths dayside magnetopause and in the magnetotail region on Earths night side. In order to resolve fine structures of the three dimensional electron distributions in both regions, the Fast Plasma Investigation's (FPI) Dual Electron Spectrometer (DES) is designed to measure electron distributions with a time resolution of 30 ms. In order to achieve this unprecedented sampling rate, the DES will have eight individual spectrometers each sampling 180 x 45 degree sections of the sky. Because of the field of view limitations of top hat analyzers, each spectrometer will use electro-static deflectors to change its look direction. The engineering model of the DES has been fabricated and tested. We will present the results of measurements for fields of view, angular FVVHM responses, dE/E, analyzer constant, and geometric factors for all deflection states. We will compare these results to simulation results and discuss causes of the response variations.

Shappirio, M.↗

DOE OSTI2023

Effect of hydrogen on plasticity of $α$-Fe: A multi-scale assessment

A multi-scale study was carried out to quantify the effect of interstitial hydrogen concentration on plasticity in $α$-Fe. Here, in this work, the influence of hydrogen on the screw dislocation glide behavior was examined across several length-scales. The insights obtained were integrated to provide an accurate continuum description for the effect of hydrogen on the dislocation based plasticity in polycrystalline $α$-Fe. At the outset of this work, a new Fe—H interatomic potential was formulated that enhanced the atomistic estimation of the variation in dislocation glide behavior in presence of hydrogen. Next, the dislocation core reconstruction observed due to the addition of hydrogen using atomistic simulations was validated with the help of large-scale DFT calculations based on the DFT-FE framework. Several atomistic simulations were carried out to comprehensively quantify the effect of hydrogen on the non-Schmid behavior exhibited during the dislocation glide in $α$-Fe. Finally, crystal plasticity simulations were carried out to understand the effect of hydrogen on the meso-scale deformation behavior of polycrystalline $α$-Fe.

36 MATERIALS SCIENCE↗