Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Multiscale systems”

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 199 records · Page 11

Understanding Solution State Conformation and Aggregate Structure of Conjugated Polymers via Small Angle X-ray Scattering

Donor-acceptor (D-A) conjugated polymers are high-performance organic electronic materials that exhibit complex aggregation behavior. Understanding the solution state conformation and aggregation of conjugated polymers is crucial for controlling morphology during thin-film deposition and the subsequent electronic performance. However, a precise multiscale structure of solution state aggregates is lacking. Here, we present an in-depth small-angle X-ray scattering (SAXS) analysis of the solution state structure of an isoindigo-bithiophene-based D-A polymer (PII-2T) in chlorobenzene and decane as our primary system. Modeling the system as a combination of hierarchical fibrillar aggregates mixed with dispersed polymers, we extract information about conformation and multiscale aggregation and also clarify the physical origin of features often observed but unaddressed or misinterpreted in small-angle scattering patterns of conjugated polymers. The persistence length of the D-A polymer extracted from SAXS agrees well with a theoretical model based on the dihedral potentials. Additionally, we show that the broad high q structure factor peak seen in scattering profiles can be attributed to lamellar stacking occurring within the fibril aggregates and that the low q aggregate scattering is strongly influenced by the polymer molecular weight. Overall, the SAXS profiles of D-A polymers in general exhibit a sensitive dependence on the co-existence of fibrillar aggregate and dispersed polymer chain populations. We corroborate our findings from SAXS with electron microscopy of freeze-dried samples for direct imaging of fibrillar aggregates. Finally, we demonstrate the generality of our approach by fitting the scattering profiles of a variety of D-A polymers based on thieno-isoindigo (PTII-2T), diketopyrrolopyrrole (DPP2T-TT, DPP-BTZ, PDPP2FT-C- 16 ), naphthalenediimide (P(NDI2OD-T2)), and a conjugated block copolymer P3HT-b-DPPT-T. The results presented here establish a picture of the D-A polymer solution state structure and provide a general method of interpreting and analyzing their scattering profiles.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale modeling of metal-hydride interphases—quantification of decoupled chemo-mechanical energies

Abstract The quantification of interphase properties between metals and their corresponding hydrides is crucial for modeling the thermodynamics and kinetics of the hydrogenation processes in solid-state hydrogen storage materials. In particular, interphase boundary energies assume a pivotal role in determining the kinetics of nucleation, growth, and coarsening of hydrides, alongside accompanying morphological evolution during hydrogenation. The total interphase energy arises from both chemical bonding and mechanical strains in these solid-state systems. Since these contributions are usually coupled, it is challenging to distinguish via conventional computational approaches. Here, a comprehensive atomistic modeling methodology is developed to decouple chemical and mechanical energy contributions using first-principles calculations, of which feasibility is demonstrated by quantifying chemical and elastic strain energies of key interfaces within the FeTi metal-hydride system. Derived materials parameters are then employed for mesoscopic micromechanical analysis, predicting crystallographic orientations in line with experimental observations. The multiscale approach outlined verifies the importance of the chemo-mechanical interplay in the morphological evolution of growing hydride phases, and can be generalized to investigate other systems. In addition, it can streamline the design of atomistic models for the quantitative evaluation of interphase properties between dissimilar phases and allow for efficient predictions of their preferred phase boundary orientations.

36 MATERIALS SCIENCE↗

CEGANN: Crystal Edge Graph Attention Neural Network for multiscale classification of materials environment

Abstract We introduce Crystal Edge Graph Attention Neural Network (CEGANN) workflow that uses graph attention-based architecture to learn unique feature representations and perform classification of materials across multiple scales (from atomic to mesoscale) and diverse classes ranging from metals, oxides, non-metals to hierarchical materials such as zeolites and semi-ordered mesophases. CEGANN can classify based on a global, structure-level representation such as space group and dimensionality (e.g., bulk, 2D, clusters, etc.). Using representative materials such as polycrystals and zeolites, we demonstrate its transferability in performing local atom-level classification tasks, such as grain boundary identification and other heterointerfaces. CEGANN classifies in (thermal) noisy dynamical environments as demonstrated for representative zeolite nucleation and growth from an amorphous mixture. Finally, we use CEGANN to classify multicomponent systems with thermal noise and compositional diversity. Overall, our approach is material agnostic and allows for multiscale feature classification ranging from atomic-scale crystals to heterointerfaces to microscale grain boundaries.

36 MATERIALS SCIENCE↗

Hyporheic‐Zone Processes and Stream Oxygen Dynamics: Insights From a Multiscale Reactive Transport Model

Aquatic ecosystem metabolism encapsulates the daily fixation (gross primary production, GPP d ) and mineralization (ecosystem respiration, ER d ) of organic carbon. In fluvial systems, these are commonly estimated by inverse solutions to field observations using a model that describes oxygen concentrations varying in the water column in response to metabolic fluxes and air‐water gas exchange controlled by a rate coefficient (K 600 ). The most common conceptual model is the single‐station metabolism (SSM) model. The simplicity and flexibility of this conceptualization make it attractive; however, it implicitly assumes that all the processes that consume oxygen in fluvial systems can be lumped into a bulk estimate of respiration with poorly understood consequences for estimates of GPP d , ER d , and K 600 . Here, we focus on the implications of using SSM conceptualization when estimating metabolic fluxes from oxygen dynamics in channels where hyporheic exchange occurs. We use a new multiscale numerical model for reactive transport in streams that represents hyporheic exchange and streambed heterotrophic respiration. Nondimensionalization of this model reveals dimensionless groups that collectively control oxygen dynamics. Numerical experiments offer a mechanistic understanding of the impacts of hyporheic exchange on diel oxygen dynamics revealing that potential biases arise from neglecting mass transfer limitations. Specifically, we found that hyporheic exchange significantly affects diel oxygen dynamics, even for nonreactive streambed sediments. Moreover, while the SSM performs well in many situations, we find conditions where significant bias is produced by hyporheic exchange, even when oxygen data are well‐fitted. These situations pose a major challenge in the interpretation of metabolism assessment estimates.

Gomez‐Velez, Jesus D. [Oak Ridge National Laborato↗

Revealing EDL-driven reduction mechanisms in binary, ternary, and quaternary fluorinated electrolytes via an integrated MD–DFT–ML framework

Accurately predicting solid electrolyte interphase (SEI) formation requires explicitly resolving the electric double layer (EDL) structure, which deviates significantly from that of the bulk electrolyte. Although an established molecular dynamics (MD) and Density Functional Theory (DFT) framework can model SEI formation by evaluating reduction reactions of local clusters in the EDL, it suffers from a combinatorial computational bottleneck. To overcome this limitation, we introduce a machine-learning-accelerated simulation workflow (MD–DFT–ML), integrating a gradient-boosted regression model trained on EDL composition data to efficiently predict reduction potentials. We apply this framework to seven fluorinated electrolytes comprising fluorinated anions, a fluorinated ester solvent, two types of diluent (ion-solvating ester vs. non-solvating ether), and an FEC additive. The analysis shows that the EDL selectively accumulates cation-binding species; consequently, the non–cation-binding ether diluent rarely enters the EDL and makes minimal contributions to SEI formation. DFT calculations on statistically representative EDL clusters provide reduction potentials and fluorine-release pathways, while the ML model, which substantially reduces the DFT workload, predicts cluster reduction energies with a mean absolute error of 0.1 eV. The combined MD–DFT–ML approach also quantifies contributions from different sources to LiF formation in the SEI. This methodology establishes a generalizable route for multiscale modeling electrolyte and interphase design for next-generation electrochemical energy-storage systems.

DFT-MD-ML workflow↗

C3F: Collaborative Container-based Model Coupling Framework

Solving complex real-world grand challenge problems requires in-depth collaboration of researchers from multiple disciplines. Such collaboration often involves harnessing multiscale and multi-dimensional data and combining models from different fields to simulate systems. However, the progress on this front has been limited mainly due to significant gaps in domain knowledge and tools that are typically employed in silos of the domains. Researchers from different fields face considerable barriers to understanding and reusing each other’s data/models in order to collaborate effectively. For example, in solving the global sustainability problems, researchers from hydrology, climate science, agriculture, and economics need to run their respective models to study different components of the global and local food, energy and water systems while, at the same time, need to interact with other researchers and integrate the results of one model with another. Developing this kind of model coupling workflow calls for (1) a large amount of data being processed and exchanged across domains and organizations, (2) identifying and processing the output of one model to make it ready for integration into another model, (3) controlling the workflow dynamically so that it runs until a certain convergence condition or other criteria is met, and (4) close collaboration among the modelers to explore, tune, and test the configuration and data transformation needed to link the models. We have developed C3F, a flexible collaborative model coupling framework to help researchers accelerate their model integration and linking efforts by leveraging advanced cyberinfrastructure such as high-performance computing and virtual containers. In this paper, we describe our experience and lessons learned in developing this cyberinfrastructure solution to support the linking of Water Balance Model (WBM) and SIMPLE-G agricultural economic model in an NSF funded INFEWS project and a DOE-funded Program on Coupled Human and Earth Systems (PCHES) to study the implications of groundwater scarcity for food-energy-water systems. The C3F model coupling framework can be extended to facilitate other model linkages as well.

containerization↗

Cooperative Research and Development Agreement among National Energy Technology Laboratory, Advanced Manufacturing LLC, and Connecticut Center for Advanced Technology [Abstract]

The National Energy Technology Laboratory (NETL), Advanced Manufacturing LLC (AMLLC) and Connecticut Center for Advanced Technology (CCAT) will collaborate to develop and manufacture cost-effective, oxide dispersion-strengthened, NiCrFeCo-rich high entropy alloys that are superior to Nickel-based superalloys (e.g. IN740) for repair or replacement service in extreme environments. The new alloys can result in higher efficiency and reduced cost of advanced power systems such as advanced ultra-supercritical (AUSC) boilers and gas turbines. High throughput, multiscale computational modeling will be performed to accelerate alloy discovery by interrogating the intrinsic properties of the alloys including thermodynamic, kinetic, and mechanical properties. Validation and evaluation will be performed, including additive manufacturing, microstructure characterization, and mechanical properties tests.

36 MATERIALS SCIENCE↗

Understanding Ionizing Radiation-Induced Speciation, Chemistry, and Transport in Nuclear Materials

The research to be presented will cover an overview of three projects supported by the United States Department of Energy Office of Science Solar Photochemistry Program: (i) Radiation-Induced Late Actinide Redox Chemistry (2019–2022); (ii) Understanding Radiation-Induced Iodine Speciation, Chemistry, and Transport in High Temperature Molten Salts (2022–2025); and (iii) Radiation-Induced Chemistry of Nuclear Materials (2023–2026). The overarching aim of these distinct projects is to develop a deep fundamental mechanistic understanding of ionizing-radiation-induced processes occurring in solutions over multiple time and length scales. Investigated solutions range from aqueous to organic to molten salt liquids in the presence and absence of ligands and various redox-active ions and molecules, including the actinides. The ultimate goal of these projects is to utilize the new foundational knowledge gained to construct quantitative multiscale computer models with the capacity to predict radiation-induced behavior in complex systems of relevance to nuclear energy technologies and beyond.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Multiscale modeling high-order methods and data-driven modeling

Projection-based reduced-order models (ROMs) comprise a promising set of data-driven approaches for accelerating the simulation of high-fidelity numerical simulations. Standard projection-based ROM approaches, however, suffer from several drawbacks when applied to the complex nonlinear dynamical systems commonly encountered in science and engineering. These limitations include a lack of stability, accuracy, and sharp a posteriori error estimators. This work addresses these limitations by leveraging multiscale modeling, least-squares principles, and machine learning to develop novel reduced-order modeling approaches, along with data-driven a posteriori error estimators, for dynamical systems. Theoretical and numerical results demonstrate that the two ROM approaches developed in this work - namely the windowed least-squares method and the Adjoint Petrov - Galerkin method - yield substantial improvements over state-of-the-art approaches. Additionally, numerical results demonstrate the capability of the a posteriori error models developed in this work.

97 MATHEMATICS AND COMPUTING↗

Convergence and Quantum Advantage of Trotterized MERA for Strongly-Correlated Systems

Strongly-correlated quantum many-body systems are difficult to study and simulate classically. We recently proposed a variational quantum eigensolver (VQE) based on the multiscale entanglement renormalization ansatz (MERA) with tensors constrained to certain Trotter circuits. Here, we determine the scaling of computation costs for various critical spin chains which substantiates a polynomial quantum advantage in comparison to classical MERA simulations based on exact energy gradients or variational Monte Carlo. Algorithmic phase diagrams suggest an even greater separation for higher-dimensional systems. Hence, the Trotterized MERA VQE is a promising route for the efficient investigation of strongly-correlated quantum many-body systems on quantum computers. Furthermore, we show how the convergence can be substantially improved by building up the MERA layer by layer in the initialization stage and by scanning through the phase diagram during optimization. For the Trotter circuits being composed of single-qubit and two-qubit rotations, it is experimentally advantageous to have small rotation angles. We find that the average angle amplitude can be reduced considerably with negligible effect on the energy accuracy. Benchmark simulations suggest that the structure of the Trotter circuits for the TMERA tensors is not decisive; in particular, brick-wall circuits and parallel random-pair circuits yield very similar energy accuracies.

Miao, Qiang [Duke Quantum Center, Duke University,↗

Multiscale Patterning from Competing Interactions and Length Scales

We live in an era of complexity marked by impressive new tools powering the scientific method to accelerate discovery, prediction, and control of increasingly complex systems. In common with many disciplines and societal challenges and opportunities, materials and condensed matter sciences are beneficiaries. The volume and fidelity of experimental, computational, and visualization data available, and tools to rapidly interpret them, are remarkable. Conceptual frameworks, including multiscale, multiphysics modeling of this complexity, are fueled by the data and, in turn, guide directions for future experimental and computational strategies. In this spirit, we discuss the importance of competing interactions, length scales, and constraints as pervasive sources of spatio-temporal complexity. We use representative examples drawn from materials and condensed matter, including the important role of elasticity in some technologically important quantum materials.

36 MATERIALS SCIENCE↗

Solvent selection for biphasic extraction of 5-hydroxymethylfurfural via multiscale modeling and experiments

We introduce a comprehensive conceptual framework for selecting solvents for reactive extraction in biphasic organic-water systems and demonstrate it for the separation of HMF (5-hydroxylmethylfurfural), a platform chemical produced in the acid-catalyzed dehydration of hexoses. We first perform in silico screening of ~2500 solvents, from the ADFCRS-2018 database using the ADF COSMO-RS implementation, and classification, based on the solvent partition coefficient. We then determine experimentally the partition coefficients for HMF, fructose, and products of HMF rehydration (levulinic acid (LA), and formic acid (FA)), the mutual water-organic solvent solubilities, and the separation factors in >50 select solvents spanning multiple homologous series at room temperature and a typical reaction temperature with in situ sampling. We find that COSMO-RS is excellent for screening purposes (typical error in most cases within a factor of ~2). Increased temperatures lead to significant reduction in partitioning, and room temperature measurements are clearly inadequate for solvent selection. Upon down selecting classes of solvents based on separation performance, we perform experimental thermal stability and reaction compatibility studies of a small set of solvents at relevant reactive-extraction temperatures. We discover that many substituted phenols exhibit an order-of-magnitude increase in partitioning compared to conventional solvents due chiefly to hydrogen bond interactions and show the necessary stability but retain a significant fraction of water and LA, factors that need to be considered in technoeconomic analysis. In contrast, anilines, aldehydes, and acids are good to excellent regarding separation but incompatible with this specific reaction medium. Furthermore, this multifaceted framework can be extended to other biomass-derived products and processes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multiscale Shear Properties and Flow Performance of Milled Woody Biomass

One dominant challenge facing the development of biorefineries is achieving consistent system throughput with highly variant biomass feedstock quality and handling performance. Current handling unit operations are adapted from other sectors (primarily agriculture), where some simplifying assumptions about granular mechanics and flow performance do not translate well to a highly compressible and anisotropic material with nonlinear time- and stress-dependent properties. This work explores the shear and frictional properties of loblolly pine at multiple experimental test apparatus and particle scales to elucidate a property window that defines the shear behavior over a range of material attributes (particle size, size distribution, moisture content, etc.). In general, it was observed that the bulk internal friction and apparent cohesion depend strongly on both the stress state of the sample in granular shear testers and the overall particle size and distribution span. For equipment designed to characterize the quasi-static shear stress failure of bulk materials ranging from 50 to 1,000 ml in test volume, similar test results were observed for finely milled particles (50% passing size of 1.4 mm) with a narrow size distribution (span between 10 and 90% passing size of 0.9 mm), while stress chaining and over-torque issues persisted for the bench-scale test apparatus for larger particle sizes or widely dispersed sample sizes. Measurement of the anisotropic particle–particle friction ranged from coefficients of approximately 0.20 to 0.45 and resulted in significantly higher and more variable friction measurements for larger particle sizes and in perpendicular alignment orientations. To supplement these laboratory-scale properties, this work explores the flow of loblolly pine and Douglas fir through a pilot-scale wedge-shaped hopper and a screw feeder. For the gravity-driven hopper flow, the critical arching distance and mass discharge rate ranged from approximately 10 to 30 mm and 2 to 16 tons/hour, respectively, for both materials, where the arching distance depends strongly on the overall particle size and depends less on the hopper inclination angle. Comparatively, the auger feeder was found to be much more impacted by the size of the particles, where smaller particles had a more consistent and stable flow while consuming less power.

09 BIOMASS FUELS↗

Computer Simulation of Proton Transport in Fuel Cell Membranes (Final Report)

This DOE-supported research grant focused on understanding the nature of proton transport in complex systems such as proton exchange membranes (PEMs). The most unique aspect of the research was the development and implementation of a novel multiscale reactive molecular dynamics (MS-RMD) methodology. In this approach, covalent bonds can dynamically break and form, allowing one to accurately treat the proton hopping process essential to capturing the physics of proton transport. The Voth group applied this method to proton exchange membrane systems, providing insight into their proton transport mechanism. They found that protons can diffuse most rapidly in the water-rich regions, but that protons actually spend so little time in such regions that transport along the hydrophobic – hydrophilic interface controls the membrane performance. In addition, the group worked to increase understanding of acidic solutions, developing methods for simulating and interpreting experimental infrared vibrational spectroscopy for excess protons (acidic solutions). The group also implemented novel tools for developing proton transport reactive MD models using a relative entropy minimization scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Integrating Ultra-Coarse-Grained Protein Models into Accessible Workflows for Multiscale Molecular Dynamics

To capture protein conformational transitions using molecular dynamics (MD), several simulation resolutions covering different spatial and temporal scales are typically needed. All-atom (AA) simulations provide fine resolution, but are computationally infeasible for large systems over longer durations. Coarse-grained (CG) and ultra-coarse-grained (UCG) models have a lower resolution and computational cost while still being able to conserve essential protein features. Prior work on a Multiscale Machinelearned Modeling Infrastructure (MuMMI) combined both AA and CG simulations to study RAS-RAF protein interactions, leveraging CG models for longer time scales and using AA to investigate unusual conformations in greater detail. However, MuMMI is still resource-intensive, and this study aims to maximize exploration of the protein conformational space while reducing computational cost. In this paper, we build on prior work that integrates UCG models based on heterogeneous elastic network modeling (hENM) into the MuMMI workflow. We demonstrate that UCG models enable accurate sampling of protein conformations, focusing on simulating RAS-RAF protein interactions. Using higher-resolution CG Martini simulation data, we can automatically refine intramolecular interactions in UCG models. We present a scalable Python package that uses fluctuations observed in higher-resolution CG Martini simulations to estimate bond coefficients of the UCG model. We built novel machine learning-based backmapping methods to recover more detailed CG Martini structures from UCG structures, using diffusion models to learn the mapping between scales. Finally, we present UCG-mini-MuMMI, an accessible and less compute-intensive version of MuMMI as a resource for the scientific community. Incorporating UCG models into MD studies is applicable to a broad range of systems and proteins, and our study offers insights into the advantages and limitations of these methods.

Chemical structure↗

Subseasonal Tropical Convection Characteristics in the Energy Exascale Earth System Model Version 2

Accurate simulation of subseasonal tropical moist convection remains a key challenge for Earth system models. The difficulties stem from the reliance of cumulus cloud processes on model parameterizations and the need to represent the multiscale nature of interactions among clouds, radiation, moisture, circulation, and surface energy fluxes. Equatorial convection drives circulation anomalies that can affect weather patterns and extremes globally, motivating efforts to better understand and simulate these tropical disturbances. Here, a detailed review of subseasonal tropical convective behavior as simulated in the Energy Exascale Earth System Model version 2 (E3SMv2) is presented, with comparison to its predecessor version 1 (E3SMv1) and reference data sets. Model structural changes to the deep convective trigger function and surface fluxes, along with parametric tuning of the cloud and microphysics schemes, together result in an improved depiction of organized tropical convection across scales. In particular, E3SMv2 exhibits a more realistic Madden‐Julian oscillation (MJO) and low‐frequency Kelvin waves—owing to a sharper time mean equatorial meridional moisture gradient and improved convection‐circulation coupling —as well as a better depiction of MJO Northern Hemisphere teleconnections. Despite these improvements, subseasonal precipitation variance continues to be strongly underestimated in E3SMv2. Use of a cloud plume model also reveals that the coupling between daily averaged tropical precipitation and lower tropospheric instability in E3SM is inconsistent with observations, a bias that could potentially impact the simulation of intraseasonal disturbances.

54 ENVIRONMENTAL SCIENCES↗

Multiscale Modeling of the Mechanical Response of Silicon Carbide Composite Within the Accelerated Fuel Qualification Framework

The accelerated fuel qualification (AFQ) framework has been used for the initial development of multiscale modeling of silicon carbide (SiC) fiber reinforced composite (SiC-SiC). The AFQ framework provides a methodology to leverage physics-informed multiscale modeling along with a reduced set of empirical test data to reduce the time and cost of licensing and qualification of new nuclear fuel systems while maintaining the overall nuclear power plant safety case. SiC-SiC is being proposed for in-core applications, most notably fuel cladding, for current and next-generation nuclear reactors because of its high temperature stability, irradiation tolerance, and ability to withstand many accident conditions. As these composites exhibit multiscale architectures and complex microstructure-based fracture mechanics, it is an appealing use case for the AFQ methodology. While the end goal of this work is a single multiscale model that can be used for predictive in-core performance, current focus is on the individual various length scale models. Four individual models have been initially developed from microscale to engineering system level to capture key physics-based effects across different length scales. These models include a microscale homogenized tow model, a mesoscale fast Fourier transform–based weave model that integrates the homogenized tow model, a mesoscale finite element–based weave model, and a system-level BISON fuel performance model. Results of these models have undergone an initial comparison with separate-effects test data showing a good match to experimental results. By using the AFQ framework during model development, several near-term benefits have been secured including a reduction in development time for the SiC-SiC cladding, more targeted irradiation testing, and a better understanding of uncertainty.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Unraveling kinking: A plasticity enhancing failure mode in high strength nano metallic laminates

Kinking is an important and plasticity-enhancing deformation/failure mode in numerous mechanically anisotropic materials including high-strength nano metallic laminates (NMLs). However, our current limited understanding of the mechanics of kinking and its dependence on microstructural attributes is insufficient for thoroughly comprehending and eventually being able to control failure behaviors of materials. In this study, we investigate kinking dependencies on microstructural attributes in NMLs via in situ micropillar compression, multiscale microstructure characterization, dislocation dynamic simulations, and crystal plasticity modeling. Additionally, by examining several NML systems (Cu/Fe, Ag/Fe, Al-4Mg/Fe), we demonstrate that the development of internal stresses during loading activates local layer-parallel glide triggering kinking in NMLs. Furthermore, this work reveals the effect of key microstructural features including layer thickness, layer waviness, interface barrier strength, and work hardening capacity on kink band formation in NMLs. More broadly, our efforts represent a generically applicable approach for probing large-strain deformation behavior of complex materials via synergetic modeling and experimental efforts.

36 MATERIALS SCIENCE↗