Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale methods”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A multiresolution adaptive wavelet method for nonlinear partial differential equations

We report the multiscale complexity of modern problems in computational science and engineering can prohibit the use of traditional numerical methods in multi-dimensional simulations. Therefore, novel algorithms are required in these situations to solve partial differential equations (PDEs) with features evolving on a wide range of spatial and temporal scales. To meet these challenges, we present a multiresolution wavelet algorithm to solve PDEs with significant data compression and explicit error control. We discretize in space by projecting fields and spatial derivative operators onto wavelet basis functions. We provide error estimates for the wavelet representation of fields and their derivatives. Then, our estimates are used to construct a sparse multiresolution discretization which guarantees the prescribed accuracy. Additionally, we embed a predictor-corrector procedure within the temporal integration to dynamically adapt the computational grid and maintain the accuracy of the solution of the PDE as it evolves. We present examples to highlight the accuracy and adaptivity of our approach.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Water Imbibition and Oil Recovery in Shale: Dynamics and Mechanisms Using Integrated Centimeter-to-Nanometer-Scale Imaging

Water imbibition, and the associated oil displacement, is an important process in shale oil reservoirs after hydraulic fracturing and in water-based enhanced oil recovery (EOR). Current techniques for water imbibition measurement are mostly “black-box”-type methods. A more explicit understanding of the water imbibition/oil recovery dynamics and geological controls is in demand. In this paper, a multiscale imaging technique that covers centimeter to nanometer scale (i.e., core to pore scale), integrating neutron radiography, microcomputed tomography (micro-CT), and scanning electron microscope (SEM) is applied to investigate the water imbibition depth and rate and the cause of heterogeneity of imbibition in shale samples. The dynamic processes of water imbibition in the 1-in. (25.4-mm) core sample were explicitly demonstrated, and the imbibition along the matrix and imbibition through microfractures are distinguished through neutron radiography image analysis. The causes of observed imbibition heterogeneity were further investigated through micro-CT and SEM image analysis for 1.5-mm diameter miniplug samples from different laminas of the 1-in. core samples. Imbibition depth and rate were calculated on the basis of image analysis as well. Estimation of oil recovery through water imbibition in shale matrix was performed for an example shale field. Here this innovative and integrated multiscale imaging technique provides a “white/gray-box” method to understand water imbibition and water-oil displacement in shale. The wide span of the length scale (from centimeter to nanometer) of this technique enables a more comprehensive, accurate, and specific understanding of both the core-scale dynamics and pore-scale mechanisms of water imbibition, oil recovery, and matrix-fracture interaction.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2024

Molecular Dynamics Simulation of Complex Reactivity with the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) Software Package

Simulating chemically reactive phenomena such as proton transport on nanosecond to microsecond and beyond time scales is a challenging task. Ab initio methods are unable to currently access these time scales routinely, and traditional molecular dynamics methods feature fixed bonding arrangements that cannot account for changes in the system’s bonding topology. The Multiscale Reactive Molecular Dynamics (MS-RMD) method, as implemented in the Rapid Approach for Proton Transport and Other Reactions (RAPTOR) software package for the LAMMPS molecular dynamics code, offers a method to routinely sample longer time scale reactive simulation data with statistical precision. RAPTOR may also be interfaced with enhanced sampling methods to drive simulations toward the analysis of reactive rare events, and a number of collective variables (CVs) have been developed to facilitate this. Key advances to this methodology, including GPU acceleration efforts and novel CVs to model water wire formation are reviewed, along with recent applications of the method which demonstrate its versatility and robustness.

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

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

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

Local-nonlocal coupling in Emu/PDMS

A technique called the splice method for coupling local to peridynamic subregions of a body is described. The method relies on ghost nodes, whose values of displacement are interpolated from nearby physical nodes, to make each subregion visible to the other. In each time step, the nodes in each subregion treat the nodes in the other subregion as boundary conditions. Adaptively changing the subregions is possible through the creation and deletion of ghost nodes. Example problems in 2D and 3D illustrate how the method is used to perform multiscale modeling of fracture and impact events within a larger structure.

42 ENGINEERING↗
DOE OSTI2020

Multimodal, Multidimensional, and Multiscale X-ray Imaging at the National Synchrotron Light Source II

Over the last couple of decades, the synchrotron radiation research community has witnessed tremendous advancement in the field of x-ray imaging and microscopy. Continuing enhancement of the light sources’ brightness, advances in x-ray focusing optics, incorporation of precision instruments, and development of innovative imaging techniques are some of the leading contributors for the rapid progress. New imaging and microscopy beamlines, such as the ones at the National Synchrotron Light Source II (NSLS-II), a Department of Energy (DOE) Office of Science user facility located at DOE’s Brookhaven National Laboratory, are capable of performing more sophisticated and complicated measurements than ever before, either on their own or used together as a suite of tools. These sophisticated and complex measurements exhibit attributes for multimodal, multidimensional, and multiscale imaging. Recent popularity of these methods is strongly driven by the current trends in materials characterizations, where researchers desire to map out hierarchical materials structure over a large range of length scales, to understand structure-property correlation, to quantify materials structures in 3D, and/or to perform operando or in situ experiments. It is important to emphasize that “materials” under investigation are not only synthesized materials but also the natural materials. In this article, we describe multimodal, multidimensional and multiscale x-ray imaging capabilities of the NSLS-II beamlines and how these methods are used to tackle complex scientific problems.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Beyond quantum cluster theories: multiscale approaches for strongly correlated systems

The degrees of freedom that confer to strongly correlated systems their many intriguing properties also render them fairly intractable through typical perturbative treatments. For this reason, the mechanisms responsible for their technologically promising properties remain mostly elusive. Computational approaches have played a major role in efforts to fill this void. In particular, dynamical mean field theory and its cluster extension, the dynamical cluster approximation have allowed significant progress. However, despite all the insightful results of these embedding schemes, computational constraints, such as the minus sign problem in quantum Monte Carlo (QMC), and the exponential growth of the Hilbert space in exact diagonalization (ED) methods, still limit the length scale within which correlations can be treated exactly in the formalism. A recent advance aiming to overcome these difficulties is the development of multiscale many body approaches whereby this challenge is addressed by introducing an intermediate length scale between the short length scale where correlations are treated exactly using a cluster solver such QMC or ED, and the long length scale where correlations are treated in a mean field manner. At this intermediate length scale correlations can be treated perturbatively. This is the essence of multiscale many-body methods. Furthermore, we will review various implementations of these multiscale many-body approaches, the results they have produced, and the outstanding challenges that should be addressed for further advances.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Automated and Accelerated Continuum Model Development for Electrochemical Systems (Abbreviated Report)

Despite the availability of computational resources and advancements in numerical computing capabilities, the multiscale models core to understanding, predicting the behaviors of, and designing energy and environmental systems involving porous media are still 1.) developed through by-hand derivations and 2.) limited by many methodological assumptions employed during model derivation. As a result, the advancement of effective media models for engineering DOE mission-critical systems (e.g., batteries, flow batteries, electrolyzers, geothermal systems, subsurface chemical storage systems, etc.) is slow (i.e., it takes years for models to traverse from stages of “development” to “practical utilization”), hindering our ability to effectively optimize such systems and stay at the cutting-edge of the energy frontier. In this work, we aimed to address these limitations by 1.) automating and accelerating multiscale model derivation via symbolic computing and 2.) develop a novel multiscale modeling methodology for flow and transport through porous media that avoids the typical assumptions hindering previous models. As a result of our efforts, we 1.) developed a hybrid symbolic-numeric code called Fouriera for fully-automating the implementation of multiphysical and phase-field models via the Fourier spectral method for materials science research, and 2.) advanced a multiscale modeling methodology called The Method of Finite Averages that rigorously predicts the behaviors of flow and transport through heterogeneous porous media under the influence of non-local effects and strong advection. Ultimately, these deliverables provide strong foundations from which further efforts can advance multiscale modeling tools and capabilities that do not intrinsically rely on 1.) the speed and mathematical capabilities of humans, nor 2.) the methodological assumptions limiting current models.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Toward Multiscale, Multimaterial 3D Printing

Biological materials and organisms possess the fundamental ability to self-organize, through which different components are assembled from the molecular level up to hierarchical structures with superior mechanical properties and multifunctionalities. These complex composites inspire material scientists to design new engineered materials by integrating multiple ingredients and structures over a wide range. Additive manufacturing, also known as 3D printing, has advantages with respect to fabricating multiscale and multi-material structures. The need for multifunctional materials is driving 3D printing techniques toward arbitrary 3D architectures with the next level of complexity. In this paper, the aim is to highlight key features of those 3D printing techniques that can produce either multiscale or multimaterial structures, including innovations in printing methods, materials processing approaches, and hardware improvements. Several issues and challenges related to current methods are discussed. Ultimately, the authors also provide their perspective on how to realize the combination of multiscale and multimaterial capabilities in 3D printing processes and future directions based on emerging research.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Controlled patterning of crystalline domains by frontal polymerization

Materials with hierarchical architectures that combine soft and hard material domains with coalesced interfaces possess superior properties compared with their homogeneous counterparts. These architectures in synthetic materials have been achieved through deterministic manufacturing strategies such as 3D printing, which require an a priori design and active intervention throughout the process to achieve architectures spanning multiple length scales. Here we harness frontal polymerization spin mode dynamics to autonomously fabricate patterned crystalline domains in poly(cyclooctadiene) with multiscale organization. This rapid, dissipative processing method leads to the formation of amorphous and semi-crystalline domains emerging from the internal interfaces generated between the solid polymer and the propagating cure front. The size, spacing and arrangement of the domains are controlled by the interplay between the reaction kinetics, thermochemistry and boundary conditions. Small perturbations in the fabrication conditions reproducibly lead to remarkable changes in the patterned microstructure and the resulting strength, elastic modulus and toughness of the polymer. Furthermore, this ability to control mechanical properties and performance solely through the initial conditions and the mode of front propagation represents a marked advancement in the design and manufacturing of advanced multiscale materials. Drawing inspiration from biological systems in which structural complexity develops through dissipative reaction–diffusion processes, this study explores a transformative synthetic manufacturing strategy aimed at harnessing the principles underpinning morphogenic growth, unlocking new avenues for advanced materials design and fabrication. Synthetic coupled reaction-transport processes offer a versatile yet relatively underexplored method to manipulate the spatial attributes of synthetic materials10. Here we introduce an innovative manufacturing approach based on frontal ring-opening metathesis polymerization (FROMP) that draws parallels with morphogenic growth and development, enabling the formation of patterned microstructures within polymeric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Micrometer: Micromechanics transformer for predicting full field mechanical responses of heterogeneous materials

Predicting mechanical responses of heterogeneous materials across scales remains a significant challenge. Traditional computational methods often struggle with complex and multiscale nature of these materials, limiting their effectiveness in real-world applications. Here, in this paper, we introduce Micrometer, a vision transformer based deep learning model designed to predict full field mechanical responses of heterogeneous materials, bridging the gap between computer vision and solid mechanics problems. We show that Micrometer, trained on a large-scale high-resolution dataset of 2D fiber-reinforced composites, can achieve state-of-the-art performance in predicting microscale strain fields across a wide range of material properties and loading conditions. Our model demonstrates accuracy and computational efficiency in applications such as computational homogenization and multiscale modeling, reducing computational time by up to two orders of magnitude compared to conventional numerical solvers while maintaining less than 1 % errors in predicting macroscale stress fields. Furthermore, we showcase Micrometer’s adaptability through transfer learning experiments on new materials with limited data, highlighting its potential to tackle diverse scenarios in computational solid mechanics. These results represent a significant step towards AI-driven innovation in materials science, addressing the limitations of traditional numerical methods and paving the way for more efficient simulations of heterogeneous materials across various industrial applications.

Composite materials↗
DOE OSTI2025

Multiscale Explanation of the Missing Gallium Vacancy in Gallium Arsenide

Irradiation of gallium arsenide (GaAs) produces immobile vacancies and mobile interstitials. Yet, after decades of experimental investigation, the immobile Ga vacancy continues to evade detection, raising the question: where is the Ga vacancy? Static first-principles calculations predict a Ga vacancy should be readily observed. We find that short-time dynamical evolution of primary defects is the key to explaining this conundrum. Using a dynamical multiscale atomistically informed device engineering (AIDE) method, we discover that during the initial displacement damage, the Ga vacancy (3-/2-) defect level pins the Fermi level near the midgap, producing oppositely charged vacancies and interstitials. Driven by Coulomb attraction, fast As interstitials preferentially annihilate Ga vacancies. The Ga vacancy population plummets below detectable limits—and the now unpinned Fermi level recovers—before being experimentally observed. This dynamical model solves the mystery of the missing Ga vacancy and reveals the importance of a multiscale approach to explore the dynamical chemical behavior in experimentally inaccessible short-time regimes.

Diaz, Leopoldo [Sandia National Laboratories (SNL-↗
DOE OSTI2021

Scale-bridging with the extended/generalized finite element method for linear elastodynamics

This paper presents an extended/generalized finite element method for bridging scales in linear elastodynamics in the absence of scale separation. More precisely, the GFEMgl framework is expanded to enable the numerical solution of multiscale problems through the automated construction of specially-tailored shape functions, thereby enabling high-fidelity finite element modeling on simple, fixed finite element meshes. Furthermore, this introduces time-dependencies in the shape functions in that they are subject to continuous adaptation with time. The temporal aspects of the formulation are investigated by considering the Newmark-β time integration scheme, and the efficacy of mass lumping strategies is explored in an explicit time-stepping scheme. This method is demonstrated on representative wave propagation examples as well as a dynamic fracture problem to assess its accuracy and flexibility.

36 MATERIALS SCIENCE↗
DOE OSTI2023

QMMM 2023: A program for combined quantum mechanical and molecular mechanical modeling and simulations

Combined quantum mechanical and molecular mechanical (QM/MM) methods play an important role in multiscale modeling and simulations. QMMM 2023 is a general-purpose program for single-point calculations, geometry optimizations, transition-state optimizations, and molecular dynamics (MD) at the QM/MM level. It calls a QM package and an MM package to perform the required single-level calculations and combines them into a QM/MM energy by a variety of schemes. QMMM 2023 supports GAMESS-US, Gaussian, and ORCA as QM packages and Tinker as the MM package. Four types of treatments are available for embedding the QM subsystem in the MM environment: mechanical embedding with gas-phase calculations of the QM region, electronic embedding that allows polarization of the QM region by the MM environment, polarizable embedding for mutual polarization of the QM and MM regions, and flexible embedding for both mutual polarization and partial charge transfer between the QM and MM regions. Boundaries between QM and MM regions that pass through covalent bonds can be treated by several methods, including the redistributed charge (RC) scheme, redistributed charge and dipole (RCD) scheme, balanced-RC scheme, balanced-RCD scheme, screened charge scheme that takes account of charge penetration effects, and smeared charge scheme that delocalizes the MM charges near the QM–MM boundary. Geometry optimization can be done using the optimizer implemented in QMMM 2023 or the Berny optimizer in Gaussian through external calls to Gaussian. Molecular dynamics simulations can be performed at the pure-MM level, pure-QM level, fixed-partitioning QM/MM level, and adaptive-partitioning QM/MM level. As a result, the adaptive-partitioning treatments permit on-the-fly relocation of the QM–MM boundary by dynamically reclassifying atoms or groups into the QM or MM subsystems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Hybrid Domain Overlapping Method for Coupling System Thermal Hydraulics and CFD Codes

A hybrid multiscale coupling methodology based on a domain overlapping approach has been developed for coupling System Thermal Hydraulics (STH) and Computational Fluid Dynamics (CFD) codes. The method has been implemented between the modern STH code SAM and the CFD code NekRS, using the coupling tool Cardinal. The coupling aims to extend the STH code's applicability to scenarios where local momentum and energy transfers are important yet difficult for STH codes to capture, such as three-dimensional mixing. Two coupling strategies are implemented and compared: a hybrid domain overlapping method and the conventional domain decomposition method. Here, the strategies are applied to two closed-loop applications, and the present method shows superior stability behavior when compared to the domain decomposition method. Then, the present coupling method is validated against experimental data from a double T-junction experiment. The present STH/CFD coupling shows improved agreement with experimental data when compared to STH standalone simulations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Accelerating manufacturing for biomass conversion via integrated process and bench digitalization: a perspective

We present a perspective for accelerating biomass manufacturing via digitalization. We summarize the challenges for manufacturing and identify areas where digitalization can help. A profound potential in using lignocellulosic biomass and renewable feedstocks, in general, is to produce new molecules and products with unmatched properties that have no analog in traditional refineries. Discovering such performance-advantaged molecules and the paths and processes to make them rapidly and systematically can transform manufacturing practices. Furthermore, we discuss retrosynthetic approaches, text mining, natural language processing, and modern machine learning methods to enable digitalization. Laboratory and multiscale computation automation via active learning are crucial to complement existing literature and expedite discovery and valuable data collection without a human in the loop. Such data can help process simulation and optimization select the most promising processes and molecules according to economic, environmental, and societal metrics. We propose the close integration between bench and process scale models and data to exploit the low dimensionality of the data and transform the manufacturing for renewable feedstocks.

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

The Schwarz Alternating Method for the Seamless Coupling of Nonlinear Reduced Order Models and Full Order Models

Projection-based model order reduction allows for the parsimonious representation of full order models (FOMs), typically obtained through the discretization of a set of partial differential equations (PDEs) using conventional techniques (e.g., finite element, finite volume, finite difference methods) where the discretization may contain a very large number of degrees of freedom. As a result of this more compact representation, the resulting projection-based reduced order models (ROMs) can achieve considerable computational speedups, which are especially useful in real-time or multi-query analyses. One known deficiency of projection-based ROMs is that they can suffer from a lack of robustness, stability and accuracy, especially in the predictive regime, which ultimately limits their useful application. Another research gap that has prevented the widespread adoption of ROMs within the modeling and simulation community is the lack of theoretical and algorithmic foundations necessary for the “plug-and-play” integration of these models into existing multi-scale and multi-physics frameworks. This paper describes a new methodology that has the potential to address both of the aforementioned deficiencies by coupling projection-based ROMs with each other as well as with conventional FOMs by means of the Schwarz alternating method [41]. Leveraging recent work that adapted the Schwarz alternating method to enable consistent and concurrent multiscale coupling of finite element FOMs in solid mechanics [35, 36], we present a new extension of the Schwarz framework that enables FOM-ROM and ROM-ROM coupling, following a domain decomposition of the physical geometry on which a PDE is posed. In order to maintain efficiency and achieve computation speed-ups, we employ hyper-reduction via the Energy-Conserving Sampling and Weighting (ECSW) approach [13]. We evaluate the proposed coupling approach in the reproductive as well as in the predictive regime on a canonical test case that involves the dynamic propagation of a traveling wave in a nonlinear hyper-elastic material.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Machine-learning-based dynamic-importance sampling for adaptive multiscale simulations

Multiscale simulations are a well-accepted way to bridge the length and time scales required for scientific studies with the solution accuracy achievable through available computational resources. Traditional approaches either solve a coarse model with selective refinement or coerce a detailed model into faster sampling, both of which have limitations. Here, we present a paradigm of adaptive, multiscale simulations that couple different scales using a dynamic-importance sampling approach. Our method uses machine learning to dynamically and exhaustively sample the phase space explored by a macro model using microscale simulations and enables an automatic feedback from the micro to the macro scale, leading to a self-healing multiscale simulation. As a result, our approach delivers macro length and time scales, but with the effective precision of the micro scale. Our approach is arbitrarily scalable as well as transferable to many different types of simulations. Overall, our method made possible a multiscale scientific campaign of unprecedented scale to understand the interactions of RAS proteins with a plasma membrane in the context of cancer research running over several days on Sierra, which is currently the second-most-powerful supercomputer in the world.

59 BASIC BIOLOGICAL SCIENCES↗