Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale problem”

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 19 records

Efficient hybrid explicit-implicit learning for multiscale problems

Splitting method is a powerful method to handle application problems by splitting physics, scales, domain, and so on. Many splitting algorithms have been designed for efficient temporal discretization. Here, in this paper, our goal is to use temporal splitting concepts in designing machine learning algorithms and, at the same time, help splitting algorithms by incorporating data and speeding them up. We propose a machine learning assisted splitting scheme which improves the efficiency of the scheme meanwhile preserves the accuracy. We consider a recently introduced multiscale splitting algorithms, where the multiscale problem is solved on a coarse grid. To approximate the dynamics, only a few degrees of freedom are solved implicitly, while others explicitly. This splitting concept allows identifying degrees of freedom that need implicit treatment. In this paper, we use this splitting concept in machine learning and propose several strategies. First, the implicit part of the solution can be learned as it is more difficult to solve, while the explicit part can be computed. This provides a speed-up and data incorporation for splitting approaches. Secondly, one can design a hybrid neural network architecture because handling explicit parts requires much fewer communications among neurons and can be done efficiently. Thirdly, one can solve the coarse grid component via PDEs or other approximation methods and construct simpler neural networks for the explicit part of the solutions. We discuss these options and implement one of them by interpreting it as a machine translation task. This interpretation of the splitting scheme successfully enables us using the Transformer since it can perform model reduction for multiple time series and learn the connection between them. We also find that the splitting scheme is a great platform to predict the coarse solution with insufficient information of the target model: the target problem is partially given and we need to solve it through a known problem which approximates the target. Our machine learning model can incorporate and encode the given information from two different problems and then solve the target problems. We conduct four numerical examples and the results show that our method is stable and accurate.

97 MATHEMATICS AND COMPUTING↗

Scientific machine learning for closure models in multiscale problems: A review

Here, closure problems are omnipresent when simulating multiscale systems, where some quantities and processes cannot be fully prescribed despite their effects on the simulation's accuracy. Recently, scientific machine learning approaches have been proposed as a way to tackle the closure problem, combining traditional (physics-based) modeling with data-driven (machine-learned) techniques, typically through enriching differential equations with neural networks. This paper reviews the different reduced model forms, distinguished by the degree to which they include known physics, and the different objectives of a priori and a posteriori learning. The importance of adhering to physical laws (such as symmetries and conservation laws) in choosing the reduced model form and choosing the learning method is discussed. The effect of spatial and temporal discretization and recent trends toward discretization-invariant models are reviewed. In addition, we make the connections between closure problems and several other research disciplines: inverse problems, Mori-Zwanzig theory, and multi-fidelity methods. In conclusion, much progress has been made with scientific machine learning approaches for solving closure problems, but many challenges remain. In particular, the generalizability and interpretability of learned models is a major issue that needs to be addressed further.

97 MATHEMATICS AND COMPUTING↗

Scalable Multilevel Uncertainty Quantification Concepts for Extreme-Scale Multiscale Problems

Our work in this project is aimed at making fundamental advances in multilevel and multiscale methods for problems with uncertainties related to flows in porous media. The main thrust of this research is to develop a systematic multiscale, multilevel parallel approaches for problems with uncertainties and use them in the uncertainty quantification. A key emphasis is on problems without an apparent scale separation in space, time, and uncertainty spaces. Multiscale solution methods are currently under active investigation. However, the develop- ment of multiscale methods for problems with uncertainties where the scales and uncertainties are tightly coupled is a research area that is less explored. During this project, we have made some significant contributions to the developments of multiscale methods, the uncertainty quantification for multiscale methods, multiscale model learning, uncertainty quantification in inverse problems, and applications.

97 MATHEMATICS AND COMPUTING↗

NH-PINN: Neural homogenization-based physics-informed neural network for multiscale problems

Physics-informed neural network (PINN) is a data-driven approach to solving equations. It is successful in many applications; however, the accuracy of the PINN is not satisfactory when it is used to solve multiscale equations. Homogenization approximates a multiscale equation by a homogenized equation without multiscale property; it includes solving cell problems and the homogenized equation. The cell problems are periodic, and we propose an oversampling strategy that significantly improves the PINN accuracy on periodic problems. The homogenized equation has a constant or slow dependency coefficient and can also be solved by PINN accurately. We hence proposed a 3-step method, neural homogenization based PINN (NH-PINN), to improve the PINN accuracy for solving multiscale problems with the help of homogenization.

97 MATHEMATICS AND COMPUTING↗

Finite Dimensional Approximations for Continuum Multiscale Problems

The completed research project concerns the development of novel computational techniques for modeling nonlinear multiscale physical and biological phenomena. Specifically, it addresses the theoretical development and applications of the homogenization theory (coarse graining) approach to calculation of the effective properties of highly heterogenous biological and bio-inspired materials with many spatial scales and nonlinear behavior. This theory studies properties of strongly heterogeneous media in problems arising in materials science, geoscience, biology, etc. Modeling of such media raises fundamental mathematical questions, primarily in partial differential equations (PDEs) and calculus of variations, the subject of the PI’s research. The focus of completed research was on mathematical models of biological and bio-inspired materials with the common theme of multiscale analysis and coarse grain computational techniques. Biological and bio-inspired materials offer the unique ability to create environmentally clean functional materials used for energy conversion and storage. These materials are intrinsically complex, with hierarchical organization occurring on many nested length and time scales. The potential to rationally design and tailor the properties of these materials for broad energy applications has been hampered by the lack of computational techniques, which are able to bridge from the molecular to the macroscopic scale. The project addressed the challenge of computational treatments of such complex materials by the development of a synergistic approach that combines innovative multiscale modeling/analysis techniques with high performance computing.

97 MATHEMATICS AND COMPUTING↗

Mathematical Foundations for Nonlocal Interface Problems: Multiscale Simulations of Heterogeneous Materials (Final LDRD Report)

Nonlocal models provide a much-needed predictive capability for important Sandia mission applications, ranging from fracture mechanics for nuclear components to subsurface flow for nuclear waste disposal, where traditional partial differential equations (PDEs) models fail to capture effects due to long-range forces at the microscale and mesoscale. However, utilization of this capability is seriously compromised by the lack of a rigorous nonlocal interface theory, required for both application and efficient solution of nonlocal models. To unlock the full potential of nonlocal modeling we developed a mathematically rigorous and physically consistent interface theory and demonstrate its scope in mission-relevant exemplar problems.

97 MATHEMATICS AND COMPUTING↗

Mixed Arlequin method for multiscale poromechanics problems: Mixed Arlequin method for multiscale poromechanics problems

An Arlequin poromechanics model is introduced to simulate the hydro-mechanical coupling effects of fluid-infiltrated porous media across different spatial scales within a concurrent computational framework. A two-field poromechanics problem is first recast as the twofold saddle point of an incremental energy functional. We then introduce Lagrange multipliers and compatibility energy functionals to enforce the weak compatibility of hydro-mechanical responses in the overlapped domain. Here, to examine the numerical stability of this hydro-mechanical Arlequin model, we derive a necessary condition for stability, the twofold inf–sup condition for multi-field problems, and establish a modified inf–sup test formulated in the product space of the solution field. We verify the implementation of the Arlequin poromechanics model through benchmark problems covering the entire range of drainage conditions. Finally, through these numerical examples, we demonstrate the performance, robustness, and numerical stability of the Arlequin poromechanics model.

42 ENGINEERING↗

Nuclear Energy Advanced Modeling and Simulation (NEAMS) Accident Tolerant Fuels High Impact Problem: Coordinate Multiscale U 3 Si 2 Modeling

Since the events at the Fukushima-Daiichi nuclear power plant in March 2011 significant research has unfolded at national laboratories, universities and other institutions into alternative materials that have potential enhanced accident tolerance when compared to traditional \uo~fuel zircaloy clad fuel rods. One of the potential replacement fuels is uranium silicide (\usi) for its higher thermal conductivity and uranium density. The lower melting temperature is of potential concern during postulated accident conditions. Another disadvantage for \usi~ is the lack of experimental data under power reactor conditions. Due to the aggressive development schedule for inserting some of the potential materials into lead test assemblies or rods by 2022~\cite{bragg-sitton_2014} multiscale multiphysics modeling approaches have been used to provide insight into these materials. \\ \noindent The purpose of this letter report is to highlight the multiscale modeling effort for \usi~fuel as part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Accident Tolerant Fuel (ATF) High Impact Problem (HIP). This fiscal year two new models for \usi~fuel have been incorporated into the BISON fuel performance code~\cite{Williamson2012} based upon lower length scale simulations using the GRASS-SST rate theory code. The models are briefly described in the following sections.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Nuclear Energy Advanced Modeling and Simulation (NEAMS) Accident Tolerant Fuels High Impact Problem: Coordinate Multiscale FeCrAl Modeling

Since the events at the Fukushima-Daiichi nuclear power plant in March 2011 significant research has unfolded at national laboratories, universities and other institutions into alternative materials that have potential enhanced ac- cident tolerance when compared to traditional UO2 fuel zircaloy clad fuel rods. One of the potential replacement claddings are iron-chromium-alunimum (FeCrAl) alloys due to their increased oxidation resistance [1–4] and higher strength [1, 2]. While the oxidation characteristics of FeCrAl are a benefit for accident tolerance, the thermal neu- tron absorption cross section of FeCrAl is about ten times that of Zircaloy. This neutronic penalty necessitates thinner cladding. This allows for slightly larger pellets to give the same cold gap width in the rod. However, the slight increase in pellet diameter is not sufficient to compensate for the neutronic penalty and enriching the fuel beyond the current 5% limit appears to be necessary [5]. Current estimates indicate that this neutronic penalty will impose an increase in fuel cost of 15-35% [1, 2]. In addition to the neutronic disadvantage, it is anticipated that tritium release to the coolant will be larger because the permeability of hydrogen in FeCrAl is about 100 times higher than in Zircaloy [6]. Also, radiation-induced hardening and embrittlement of FeCrAl need to be fully characterized experimentally [7]. Due to the aggressive development schedule for inserting some of the potential materials into lead test assemblies or rods by 2022 [8] multiscale multiphysics modeling approaches have been used to provide insight into these the use of FeCrAl as a cladding material. The purpose of this letter report is to highlight the multiscale modeling effort for iron-chromium-alunimum (FeCrAl) cladding alloys as part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program through its Accident Tolerant Fuel (ATF) High Impact Problem (HIP). The approach taken throughout the HIP is to utilize lower length scale approaches (e.g., density functional theory, cluster dynamics, rate theory, phase field, and Visco-Plastic- Self-Consistent (VPSC)) to develop more physically informed models at the engineering scale for use in the BISON [9] fuel performance code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Multiscale Universal Interface: A concurrent framework for coupling heterogeneous solvers

Graphical abstract: - Abstract: Concurrently coupled numerical simulations using heterogeneous solvers are powerful tools for modeling multiscale phenomena. However, major modifications to existing codes are often required to enable such simulations, posing significant difficulties in practice. In this paper we present a C++ library, i.e. the Multiscale Universal Interface (MUI), which is capable of facilitating the coupling effort for a wide range of multiscale simulations. The library adopts a header-only form with minimal external dependency and hence can be easily dropped into existing codes. A data sampler concept is introduced, combined with a hybrid dynamic/static typing mechanism, to create an easily customizable framework for solver-independent data interpretation. The library integrates MPI MPMD support and an asynchronous communication protocol to handle inter-solver information exchange irrespective of the solvers' own MPI awareness. Template metaprogramming is heavily employed to simultaneously improve runtime performance and code flexibility. We validated the library by solving three different multiscale problems, which also serve to demonstrate the flexibility of the framework in handling heterogeneous models and solvers. In the first example, a Couette flow was simulated using two concurrently coupled Smoothed Particle Hydrodynamics (SPH) simulations of different spatial resolutions. In the second example, we coupled the deterministic SPH method with the stochastic Dissipative Particle Dynamics (DPD) method to study the effect of surface grafting on the hydrodynamics properties on the surface. In the third example, we consider conjugate heat transfer between a solid domain and a fluid domain by coupling the particle-based energy-conserving DPD (eDPD) method with the Finite Element Method (FEM)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multiscale high-order/low-order (HOLO) algorithms and applications

We review the state of the art in the formulation, implementation, and performance of so-called high-order/low-order (HOLO) algorithms for challenging multiscale problems. HOLO algorithms attempt to couple one or several high-complexity physical models (the high-order model, HO) with low-complexity ones (the low-order model, LO). The primary goal of HOLO algorithms is to achieve nonlinear convergence between HO and LO components while minimizing memory footprint and managing the computational complexity in a practical manner. Key to the HOLO approach is the use of the LO representations to address temporal stiffness, effectively accelerating the convergence of the HO/LO coupled system. The HOLO approach is broadly underpinned by the concept of nonlinear elimination, which enables segregation of the HO and LO components in ways that can effectively use heterogeneous architectures. The accuracy and efficiency benefits of HOLO algorithms are demonstrated with specific applications to radiation transport, gas dynamics, plasmas (both Eulerian and Lagrangian formulations), and ocean modeling. Across this broad application spectrum, HOLO algorithms achieve significant accuracy improvements at a fraction of the cost compared to conventional approaches. It follows that HOLO algorithms hold significant potential for high-fidelity system scale multiscale simulations leveraging exascale computing.

97 MATHEMATICS AND COMPUTING↗

Mapping reactive flow patterns in monolithic nanoporous catalysts

The development of high-efficiency porous catalyst membranes critically depends on our understanding of where the majority of the chemical conversions occur within the porous structure. This then requires mapping of chemical reactions and mass transport inside the complex nanoscale architecture of porous catalyst membranes which is a multiscale problem in both the temporal and spatial domains. In order to address this problem, we developed a multiscale mass transport computational framework based on the lattice Boltzmann method that allows us to account for catalytic reactions at the gas–solid interface by introducing a new boundary condition. In good agreement with experiments, the simulations reveal that most catalytic reactions occur near the gas-flow facing side of the catalyst membrane if chemical reactions are fast compared to mass transport within the porous catalyst membrane.

36 MATERIALS SCIENCE↗

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↗

Operator learning for predicting multiscale bubble growth dynamics

We report simulating and predicting multiscale problems that couple multiple physics and dynamics across many orders of spatiotemporal scales is a great challenge that has not been investigated systematically by deep neural networks (DNNs). Herein, we develop a framework based on operator regression, the so-called deep operator network (DeepONet), with the long-term objective to simplify multiscale modeling by avoiding the fragile and time-consuming “hand-shaking” interface algorithms for stitching together heterogeneous descriptions of multiscale phenomena. To this end, as a first step, we investigate if a DeepONet can learn the dynamics of different scale regimes, one at the deterministic macroscale and the other at the stochastic microscale regime with inherent thermal fluctuations. Specifically, we test the effectiveness and accuracy of the DeepONet in predicting multirate bubble growth dynamics, which is described by a Rayleigh–Plesset (R–P) equation at the macroscale and modeled as a stochastic nucleation and cavitation process at the microscale by dissipative particle dynamics (DPD). First, we generate data using the R–P equation for multirate bubble growth dynamics caused by randomly time-varying liquid pressures drawn from Gaussian random fields (GRFs). Our results show that properly trained DeepONets can accurately predict the macroscale bubble growth dynamics and can outperform long short-term memory networks. We also demonstrate that the DeepONet can extrapolate accurately outside the input distribution using only very few new measurements. Subsequently, we train the DeepONet with DPD data corresponding to stochastic bubble growth dynamics. Although the DPD data are noisy and we only collect sparse data points on the trajectories, the trained DeepONet model is able to predict accurately the mean bubble dynamics for time-varying GRF pressures. Taken together, our findings demonstrate that DeepONets can be employed to unify the macroscale and microscale models of the multirate bubble growth problem, hence providing new insight into the role of operator regression via DNNs in tackling realistic multiscale problems and in simplifying modeling with heterogeneous descriptions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗