Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale modeling”

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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Overview of the MOSAiC expedition—Atmosphere

With the Arctic rapidly changing, the needs to observe, understand, and model the changes are essential. To support these needs, an annual cycle of observations of atmospheric properties, processes, and interactions were made while drifting with the sea ice across the central Arctic during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition from October 2019 to September 2020. An international team designed and implemented the comprehensive program to document and characterize all aspects of the Arctic atmospheric system in unprecedented detail, using a variety of approaches, and across multiple scales. These measurements were coordinated with other observational teams to explore cross-cutting and coupled interactions with the Arctic Ocean, sea ice, and ecosystem through a variety of physical and biogeochemical processes. This overview outlines the breadth and complexity of the atmospheric research program, which was organized into 4 subgroups: atmospheric state, clouds and precipitation, gases and aerosols, and energy budgets. Atmospheric variability over the annual cycle revealed important influences from a persistent large-scale winter circulation pattern, leading to some storms with pressure and winds that were outside the interquartile range of past conditions suggested by long-term reanalysis. Similarly, the MOSAiC location was warmer and wetter in summer than the reanalysis climatology, in part due to its close proximity to the sea ice edge. The comprehensiveness of the observational program for characterizing and analyzing atmospheric phenomena is demonstrated via a winter case study examining air mass transitions and a summer case study examining vertical atmospheric evolution. Overall, the MOSAiC atmospheric program successfully met its objectives and was the most comprehensive atmospheric measurement program to date conducted over the Arctic sea ice. The obtained data will support a broad range of coupled-system scientific research and provide an important foundation for advancing multiscale modeling capabilities in the Arctic.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2023

NASA GRC ICME Schema for Materials Data Management: An Executive Summary

Integrated Computational Materials Engineering (ICME) has received a growing emphasis in attention due its potential impact on rapid material design, reduction in cost and time to market for new applications, and the promise of ‘fit-for-purpose’ materials coupled with recent advances in high performance computing and material characterization tools. However, for an organization to implement ICME practices for material discovery and design, a series of both technical and cultural challenges must be overcome to foster an environment that enables efficient, traceable, and predictive multiscale simulations of material behavior to enable virtual design of materials. In 2016, NASA sponsored a 2040 Vision study to define the potential 25-year future state required for integrated multiscale modeling of materials and systems to improve both the associated time and cost for aerospace and aeronautical innovation. The study envisions a cyber-physical-social ecosystem of experimentally validated computational models, tools, and techniques, along with the associated digital tapestry, that can enable rapid, optimized, ‘fit-for-purpose’ design of materials, components, and systems. A key requirement for such an ecosystem is the development of a robust information management system for materials across their full lifecycle, including material pedigree, experimental (real) and virtual (simulation) data, developed material models, and the implementation of models in engineering applications, such that process-structure-property-performance relationships can be established, thereby enabling the virtual design and optimization of materials. Such an information management system must be able to effectively capture: i) material information at each length scale; ii) test data and analysis; iii) associated material models; and iv) material and model deployment in engineering applications. These systems must also provide traceability between experimental and virtual representations of the material to ensure, when appropriate, the material digital twin is maintained. Additionally, this robust material information management system must be able to seamlessly connect with both commercial and an organization’s in-house software tools, be they analysis tools, other material databases, product lifecycle management (PLM) or simulation data management (SDM) tools, etc., such that automation of the design and analysis of a material across multiple length scales is possible. In this paper, an executive summary of the NASA GRC ICME Schema for materials information management is presented. The database best practices and schema design philosophy specifically for ICME materials data management and an overview description of each element in the schema is given, along with its associated role in an ICME workflow. Additionally, auxiliary tools that interact with the database and provide judicious automation with regards to importing, exporting, and analyzing materials data are presented. Such tools are critical to an ICME ecosystem, not only for their role in enabling optimization, but also in relieving users of tedious manual tasks, thus helping to promote adoption and combat the cultural challenges organizations face in enabling ICME.

Materials↗
DOE OSTI2024

Limitations of Hydrogen Detection After 150 Years of Research on Hydrogen Embrittlement

Hydrogen's significance in contemporary society lies in its remarkable energy density, yet its integration into the worldwide energy grid presents a substantial challenge. Exposing materials to hydrogen environments leads to degradation of mechanical properties, damage, and failure. While the current approach for assessing hydrogen's impact on materials involves mainly multiscale modeling and mechanical testing, there exists a significant deficiency in detecting the intricate interactions between hydrogen and materials at the nanoatomic scales and under in situ conditions. This perspective review highlights the experimental endeavors aimed at bridging this gap, pointing toward the imminent need for new experimental techniques that can detect and map hydrogen in materials’ microstructures and their site‐specific dependencies.

Tunes, Matheus A.↗
DOE OSTI2025

Direct Imaging of Hydrogen‐Driven Dislocation and Strain Field Evolution in a Stainless Steel Grain

Hydrogen embrittlement (HE) poses a significant challenge to the durability of materials used in hydrogen production and utilization. Disentangling the competing nanoscale mechanisms driving HE often relies on simulations and electron-transparent sample techniques, limiting experimental insights into hydrogen-induced dislocation behavior in bulk materials. This study employs in situ Bragg coherent X-ray diffraction imaging to track three-dimensional (3D) dislocation and strain field evolution during hydrogen charging in a bulk grain of austenitic 316 stainless steel. Tracking a single dislocation reveals hydrogen-enhanced mobility and relaxation, consistent with dislocation dynamics simulations. Subsequent observations reveal dislocation unpinning and climb processes, likely driven by osmotic forces. Additionally, nanoscale strain analysis around the dislocation core directly measures hydrogen-induced elastic shielding. These findings experimentally validate theoretical predictions and offer mechanistic insights into hydrogen-driven dislocation behavior. The quantified nanoscale phenomena serve as critical inputs for multiscale modeling frameworks to predict bulk material responses and accelerate the development of HE-resistant alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Theoretical Investigation of the Oxygen Reduction Reaction over Platinum Catalysts Supported by Multi-Edged Vertically Aligned Carbon Nanofiber for Electrocatalyst Preparation

Vertically aligned carbon nanofibers (VACNFs) are promising supports for oxygen reduction reaction (ORR) electrocatalysts in fuel cells. Although experimentally these catalytic systems have shown great potential, there is lack of molecular understanding of the catalytic sites and reaction mechanisms. This work investigated the origin of the ORR reactivities of the platinum catalysts on multi-edged VACNFs (Pt/VACNF) using a multiscale modeling approach combining Density Functional Theory (DFT) and classical Molecular Dynamics (MD) simulations. Based on the ReaxFF potential, here all nanoscale Pt particles (Pt 55 , P 100 , and Pt 147 ) are stabilized by the open edges located axially along the VACNF walls. The calculated first-shell coordination numbers, CN, of surface Pt atoms are 6.63, 7.27, and 7.85, respectively, suggesting that the percentage of low coordination sites increases as the particle size decreases. The adsorption energies of OOH, O, and OH on Pt-55 were systematically probed using DFT calculations. These adsorption energies retain a linear correlation against the generalized coordination numbers ($\overline{CN}$). For Pt nanoparticles supported on VACNF, we found that the OOH and OH bind stronger than on Pt (111) by 0.14 and 0.17 eV, respectively, which can hinder the ORR activity with lower limiting potential than Pt (111). Our theoretical prediction is in good agreement with the linear sweeping voltammetry that revealed a left shift of the half-wave potential.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Solar Energy Research Institute for India and the United States (SERIIUS): Lessons and Results from a Binational Consortium

This book describes the development, functioning, and results of a successful binational program to promote significant scientific advances in Earth-abundant photovoltaics (PV) and concentrated solar power (CSP), advanced process/manufacturing technologies, multiscale modeling and reliability testing, and analysis of integrated solar energy systems. SERIIUS is a consortium between India and the United States dedicated to developing new solar technologies and assessing their potential impact in the two countries. The consortium consists of nearly 50 institutions including academia, national laboratories, and industry, with the goal of developing significant new technologies in all areas of solar deployment. In addition, the program focused on workforce development through graduate students, post-doctoral students, and an international exchange program. Particular emphasis was placed on the following efforts: creating disruptive technologies in PV and CSP through high-impact fundamental and applied research and development (R&D); identifying and quantifying the critical technical, economic, and policy issues for solar energy development and deployment in India; overcoming barriers to technology transfer by teaming research institutions and industry in an effective project structure; building a new platform for binational collaboration using a formalized R&D project structure, along with effective management, coordination, and decision processes; creating a sustainable network and workforce development program from which to build large collaborations and fostering a collaborative culture and outreach programs. This includes using existing and new methodologies for collaboration based on advanced electronic and web-based communication to facilitate functional international teams. The book summarizes the general lessons learned from these experiences.

14 SOLAR ENERGY↗
DOE OSTI2021

Split representation of adaptively compressed polarizability operator

The polarizability operator plays a central role in density functional perturbation theory and other perturbative treatment of first principle electronic structure theories. The cost of computing the polarizability operator generally scales as O(Ne4) where Ne is the number of electrons in the system. The recently developed adaptively compressed polarizability operator (ACP) formulation [L. Lin, Z. Xu and L. Ying, Multiscale Model. Simul. 2017] reduces such complexity to O(Ne 3 ) in the context of phonon calculations with a large basis set for the first time, and demonstrates its effectiveness for model problems. In this paper, we improve the performance of the ACP formulation by splitting the polarizability into a near singular component that is statically compressed, and a smooth component that is adaptively compressed. The new split representation maintains the O(Ne 3 ) complexity, and accelerates nearly all components of the ACP formulation, including Chebyshev interpolation of energy levels, iterative solution of Sternheimer equations, and convergence of the Dyson equations. For simulation of real materials, we discuss how to incorporate nonlocal pseudopotentials and finite temperature effects. In this work, we demonstrate the effectiveness of our method using one-dimensional model problem in insulating and metallic regimes, as well as its accuracy for real molecules and solids.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Enhanced liquid metal wetting on oxide surfaces via patterned particles

Liquid metals typically do not intrinsically wet oxide surfaces due to chemical dissimilarity, preventing the formation of well-bonded interfaces in many applications. This paper investigates how the wetting and spreading of liquid silver (Ag) droplets on oxide surfaces can be enhanced via patterned nickel (Ni) particles using multiscale modeling. First, a force field (FF) for an Ag/yttria-stabilized-zirconia system was developed by fitting the interfacial binding energy to density-functional theory (DFT) calculations. Then, molecular dynamics (MD) simulations using these FF parameters revealed that the good wetting of Ag on Ni caused Ag droplets to be attracted to the Ni particles, promoting spreading. Thus, the optimal hexagonal Ni pattern was shown to enhance the wetting area by 224% compared to a bare oxide surface. Second, an analytical model predicting the maximum possible enhancement in wetting area as a function of the intrinsic wetting angle, size of the Ag droplet, and size of the Ni particles was derived and validated with MD simulations, to generalize this method to other substrates. Here, both the analytical model and MD results suggested that the ability of Ni patterns to enhance the Ag wetting area becomes more significant as the wettability of substrate becomes worse.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mobility of twinning dislocations in copper up to supersonic speeds

Understanding the mobility of twinning dislocations is important for multiscale modeling of crystal plasticity, especially at high strain rates, where such dislocations may reach transonic or supersonic speeds. Here, we used molecular dynamics simulations to investigate the relationship between dislocation velocity and the applied resolved shear stress of an edge twinning dislocation in copper up to supersonic speeds. The twinning dislocation mobility relation is composed of two branches separated by a band of forbidden velocities. The lower velocity branch is limited by the first transverse sound speed ~2000 m/s while the upper branch stretches from ~3500 m/s in the transonic regime to supersonic velocities. Twinning dislocations cannot undergo uniform steady-state motion at velocities within the forbidden band. Our simulation results also reveal that edge twinning dislocation motion in copper is kink-mediated. We discuss the implications of our findings for the motion of twins, twinning dislocations, and twinning dislocation kinks in copper.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evaluating grid stress and reliability in future electricity grids across a range of demand, generation mix, and weather trends

The reliability of power grids in the future will depend on how system planners account for the integration of new technologies, extreme weather events, and uncertainties in demand growth from increased electrification and data centers. This study introduces an open-source, multisectoral, multiscale modeling framework that projects grid stress and reliability trends between 2020 and 2055 in the Western Interconnection of the United States. The framework integrates global to national energy-water-land dynamics with power plant siting and hourly grid operations modeling. We analyze future wholesale electricity price shocks and unserved energy events across eight scenarios spanning a range of population growth and economic change, generation mixes, and weather conditions. Our results show future grids with high percentage of non-renewable generation and strong economic growth are characterized by higher reliability and lower wholesale electricity prices than lower growth scenarios because of larger reliance on dispatchable generators and lower fossil fuel extraction costs. Scenarios with high percentage of renewable resources have lower median but more volatile wholesale electricity prices as well as more frequent and severe unserved energy events compared to scenarios relying more on dispatchable generators. These events occur because higher proportion of solar and wind energy causes net demand curves to deepen during midday (duck curves get progressively severe), exacerbating the challenge of meeting demand during summer evening peaks. This study suggests that robust and co-optimized transmission and energy storage planning could help maintain low wholesale electricity prices and high reliability levels in future electricity grids across uncertainties in generation mixes.

Electric grid reliability↗
DOE OSTI2022

Protein-based models offer mechanistic insight into complex nickel metalloenzymes

There are ten nickel enzymes found across biological systems, each with a distinct active site and reactivity that spans reductive, oxidative, and redox–neutral processes. We focus on the reductive enzymes, which catalyze reactions that are highly germane to the modern-day climate crisis: [NiFe] hydrogenase, carbon monoxide dehydrogenase, acetyl coenzyme A synthase, and methyl coenzyme M reductase. The current mechanistic understanding of each enzyme system is reviewed along with existing knowledge gaps, which are addressed through the development of protein-derived models, as described here. Finally, this opinion is intended to highlight the advantages of using robust protein scaffolds for modeling multiscale contributions to reactivity and inspire the development of novel artificial metalloenzymes for other small molecule transformations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Mesh-based super-resolution of fluid flows with multiscale graph neural networks

A graph neural network (GNN) approach is introduced in this work which enables mesh-based three-dimensional super-resolution of fluid flows. In this framework, the GNN is designed to operate not on the full mesh-based field at once, but on localized meshes of elements (or cells) directly. To facilitate mesh-based GNN representations in a manner similar to spectral (or finite) element discretizations, a baseline GNN layer (termed a message passing layer, which updates local node properties) is modified to account for synchronization of coincident graph nodes, rendering compatibility with commonly used element-based mesh connectivities. Furthermore, the architecture is multiscale in nature, and is comprised of a combination of coarse-scale and fine-scale message passing layer sequences (termed processors) separated by a graph unpooling layer. The coarse-scale processor embeds a query element (alongside a set number of neighboring coarse elements) into a single latent graph representation using coarse-scale synchronized message passing over the element neighborhood, and the fine-scale processor leverages additional message passing operations on this latent graph to correct for interpolation errors. Demonstration studies are performed using hexahedral mesh-based data from Taylor–Green Vortex and backward-facing step flow simulations at Reynolds numbers of 1600 and 3200. Through analysis of both global and local errors, the results ultimately show how the GNN is able to produce accurate super-resolved fields compared to targets in both coarse-scale and multiscale model configurations. Reconstruction errors for fixed architectures were found to increase in proportion to the Reynolds number. Geometry extrapolation studies on a separate cavity flow configuration show promising cross-mesh capabilities of the super-resolution strategy.

Backward-facing step↗
DOE OSTI2022

hPIC2: A hardware-accelerated, hybrid particle-in-cell code for dynamic plasma-material interactions

The exascale era of high performance computing promises to bring the field of computational plasma physics ever closer to the goal of accurate multiscale modeling. Such computers will rely on hardware acceleration to offload work to dedicated components, notably general-purpose graphics processing units (GPUs). However, devices from different manufacturers require software to be written with different parallel programming models, greatly increasing the code maintenance burden of applications designed to perform on more than one such device. hPIC2 is a hybrid plasma simulation code developed with the Kokkos performance portability framework to target the architectures that will drive exascale computing for the foreseeable future. As a hybrid simulation code, hPIC2 investigates the simultaneous use of various plasma models on the same domain, at the same time. hPIC2 also optionally couples to RustBCA, which accurately models ion-material interactions using the binary collision approximation (BCA) method. In conclusion, hPIC2 therefore achieves scalable performance on a variety of computing architectures when simulating complex and diverse plasmas, particularly near plasma-material interfaces.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Parallel transport dynamics for mixed quantum states with applications to time-dependent density functional theory

Direct simulation of the von Neumann dynamics for a general (pure or mixed) quantum state can often be expensive. One prominent example is the real-time time-dependent density functional theory (rt-TDDFT), a widely used framework for the first principle description of many-electron dynamics in chemical and materials systems. Practical rt-TDDFT calculations often avoid the direct simulation of the von Neumann equation, and solve instead a set of Schrödinger equations, of which the dynamics is equivalent to that of the von Neumann equation. However, the time step size employed by the Schrödinger dynamics is often much smaller. Here, in order to improve the time step size and the overall efficiency of the simulation, we generalize a recent work of the parallel transport (PT) dynamics for simulating pure states [An, Lin, Multiscale Model. Simul. 18, 612, 2020] to general quantum states. The PT dynamics provides the optimal gauge choice, and can employ a time step size comparable to that of the von Neumann dynamics. Going beyond the linear and near adiabatic regime in previous studies, we find that the error of the PT dynamics can be bounded by certain commutators between Hamiltonians, density matrices, and their derived quantities. Such a commutator structure is not present in the Schrödinger dynamics. We demonstrate that the parallel transport-implicit midpoint (PT-IM) method is a suitable method for simulating the PT dynamics, especially when the spectral radius of the Hamiltonian is large. The commutator structure of the error bound, and numerical results for model rt-TDDFT calculations in both linear and nonlinear regimes, confirm the advantage of the PT dynamics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

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↗
DOE OSTI2022

Additively manufactured strain sensors for in-pile applications

Accurate, real-time monitoring of strain in fuel, cladding, and structural components of nuclear reactors is critical to better understand radiation induced phenomena during reactor tests and operations. The data provided is crucial to verify physics-based, multiscale modeling and simulation efforts, which aim to shorten the timeline for the development of new nuclear materials. Resistive strain gauges have limited performance during in-pile experiments due to the harsh operating conditions and limited physical space between fuel and cladding components. In this work, aerosol jet printing using silver nanoparticle inks was used to fabricate interdigitated electrode capacitive strain gauges on aluminum alloy 6061 tensile specimens. Here to simulate the temperatures of a traditional light water reactor, the capacitive strain gauges were tested with a mechanical test frame up to 300 °C and compared to commercially available bondable resistive strain gauges. The printed capacitive strain gauges exhibited a gauge factor of 1.0 and showed higher reproducibility and predictability of strain sensing performance than the resistive strain gauge. The results demonstrate the potential of aerosol jet printing to fabricate strain sensors with predictable performance and reduced invasiveness for high-temperature applications with confined spacing.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Grain2mesh: A Python and cubit mesh generator from unprocessed mesoscale images

Predicting bulk behavior from microscale features constitutes a key objective in multiscale modeling research, often involving numerical models composed of finite elements that capture the diversity of constituent phases, shapes, and orientations within the material. The Grain2mesh toolbox allows the user to input unprocessed mesoscopic images for automatic segmentation, pre-processing, quality control, and numerical mesh generation. The numerical mesh generation incorporates Cubit routines to generate robust multi-phase mesh structure for use in computational mechanics solvers. The python classes developed contain detailed documentation and examples to support standard usage and case-specific alternative options.

58 GEOSCIENCES↗