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 595 records · Page 33

ATF material model development and validation for priority cladding concepts

Development, implementation, and validation of material and behavior models for accident tolerant fuel (ATF) concepts in the Bison fuel performance code began in 2014 in response to the events that occurred at the Fukushima Daichii nuclear power plant in March 2011. Early on the focus was on U 3 Si 2 fuel and FeCrAl cladding as part of a high impact problem through the Nuclear Energy and Advanced Modeling Simulation (NEAMS) program. Then, developments for Cr 2 O 3 -doped UO 2 fuel, and SiC-SiC and Cr-coated zirconium-based claddings began based upon industry interests. In late fiscal year 2018 the Consortium for Advanced Simulation of Light Water Reactors (CASL) took over further ATF work in Bison in support of the Nuclear Regulatory Commission (NRC) engagement. Discussions with the NRC identified their list of priority fuel and cladding concepts, which included Cr 2 O 3 -doped UO 2 and U 3 Si 2 fuels, and Cr-coated zirconium-based and FeCrAl claddings. In particular, the NRC suggested that reports similar in form to NUREG/CR-7024 that was developed for traditional LWR materials UO 2 and zirconium-based claddings (i.e., Zircaloy-4, M5 ® , ZIRLO TM ) be created for the priority ATF concepts. The approach to ATF capability development in Bison since the beginning has been two-fold: (1) empirical correlations and (2) multiscale model development. Both approaches have uncertainty inherent to them. Uncertainty in empirical correlations is bounded by the experimental data upon which with the correlation was developed. Models developed through a multiscale approach have uncertainty associated with the lower length scale calculations and input parameters that must be propagated to the engineering scale model in Bison. In this report, the recommended models, their range of applicability (e.g., temperature, fluence), and associated uncertainty for the NRC priority clad concepts Cr-coated zirconium based clad and FeCrAl are presented in a manner similar to the aforementioned NUREG. In addition, in the absence of experimental data for validation, demonstration cases are presented that investigate the ballooning and burst behavior of Cr-coated and FeCrAl clads as well as tritium permeation in FeCrAl. A portion of the study by Gamble et al. on FeCrAl failure is also revisited with the latest updated models. First, the nominal results are shown using the as implemented models. Then a sensitivity analysis and uncertainty quantification is performed to assess the impact of the uncertainties in the models on the parameters of interest. For the Cr-clad ballooning case the parameter of interest is the time to burst. For the revisited FeCrAl study, the parameter of interest is the pressure at burst.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Viral kinetic modeling: state of the art

Viral kinetic modeling has led to increased understanding of the within host dynamics of viral infections and the effects of therapy. Here we review recent developments in the modeling of viral infection kinetics with emphasis on two infectious diseases: hepatitis C and influenza. We review how viral kinetic modeling has evolved from simple models of viral infections treated with a drug or drug cocktail with an assumed constant effectiveness to models that incorporate drug pharmacokinetics and pharmacodynamics, as well as phenomenological models that simply assume drugs have time varying-effectiveness. We also discuss multiscale models that include intracellular events in viral replication, models of drug-resistance, models that include innate and adaptive immune responses and models that incorporate cell-to-cell spread of infection. Overall, viral kinetic modeling has provided new insights into the understanding of the disease progression and the modes of action of several drugs. In conclusion, we expect that viral kinetic modeling will be increasingly used in the coming years to optimize drug regimens in order to improve therapeutic outcomes and treatment tolerability for infectious diseases.

59 BASIC BIOLOGICAL SCIENCES↗

tell: a Python package to model future total electricity loads in the United States

The purpose of the Total ELectricity Load (tell) model is to generate 21st century profiles of hourly electricity load (demand) across the Conterminous United States (CONUS). tell loads reflect the impact of climate and socioeconomic change at a spatial and temporal resolution adequate for input to an electricity grid operations model. tell uses machine learning to develop profiles that are driven by projections of climate/meteorology and population. tell also harmonizes its results with United States (U.S.) state-level, annual projections from a national- to global-scale energy-economy model. This model accounts for a wide range of other factors affecting electricity demand, including technology change in the building sector, energy prices, and demand elasticities, which stems from model coupling with the U.S. version of the Global Change Analysis Model (GCAM-USA). tell was developed as part of the Integrated Multisector Multiscale Modeling (IM3) project. IM3 explores the vulnerability and resilience of interacting energy, water, land, and urban systems in response to compound stressors, such as climate trends, extreme events, population, urbanization, energy system transitions, and technology change

24 POWER TRANSMISSION AND DISTRIBUTION↗

Asynchronous Reciprocal Coupling of Martini 2.2 Coarse-Grained and CHARMM36 All-Atom Simulations in an Automated Multiscale Framework

The appeal of multiscale modeling approaches is predicated on the promise of combinatorial synergy. However, this promise can only be realized when distinct scales are combined with reciprocal consistency. Here, we consider multiscale molecular dynamics (MD) simulations that combine the accuracy and macromolecular flexibility accessible to fixed-charge all-atom (AA) representations with the sampling speed accessible to reductive, coarse-grained (CG) representations. AA-to-CG conversions are relatively straightforward because deterministic routines with unique outcomes are achievable. Conversely, CG-to-AA conversions have many solutions due to a surge in the number of degrees of freedom. While automated tools for biomolecular CG-to-AA transformation exist, we find that one popular option, called Backward, is prone to stochastic failure and the AA models that it does generate frequently have compromised protein structure and incorrect stereochemistry. Although these shortcomings can likely be circumvented by human intervention in isolated instances, automated multiscale coupling requires reliable and robust scale conversion. In this work, we detail an extension to Multiscale Machine-learned Modeling Infrastructure (MuMMI), including an improved CG-to-AA conversion tool called sinceCG. This tool is reliable (~98% weakly correlated repeat success rate), automatable (no unrecoverable hangs), and yields AA models that generally preserve protein secondary structure and maintain correct stereochemistry. We describe how the MuMMI framework identifies CG system configurations of interest, converts them to AA representations, and simulates them at the AA scale while on-the-fly analyses provide feedback to update CG parameters. Application to systems containing the peripheral membrane protein RAS and proximal components of RAF kinase on complex eight-component lipid bilayers with ~1.5 million atoms is discussed in the context of MuMMI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Continuum- and particle-based modeling of shapes and dynamics of red blood cells in health and disease

Here, we review recent advances in multiscale modeling of the mechanics of healthy and diseased red blood cells (RBCs), and blood flow in the microcirculation. We cover the traditional continuum-based methods but also particle-based methods used to model both the RBCs and the blood plasma. We highlight examples of successful simulations of blood flow including malaria and sickle cell anemia.

59 BASIC BIOLOGICAL SCIENCES↗

Modeling the self-assembly of ordered nanoporous materials

This report describes progress on a collaborative project on the multiscale modeling of the assembly processes in the synthesis of nanoporous materials. Such materials are of enormous importance in modern technology with application in the chemical process industries, biomedicine and biotechnology as well as microelectronics. The project focuses on two important classes of materials: i) microporous crystalline materials, such as zeolites, and ii) ordered mesoporous materials. In the first case the pores are part of the crystalline structure, while in the second the structures are amorphous on the atomistic length scale but where surfactant templating gives rise to order on the length scale of 2 - 20 nm. We have developed a modeling framework that encompasses both these kinds of materials. Our models focus on the assembly of corner sharing silica tetrahedra in the presence of structure directing agents. We emphasize a balance between sufficient realism in the models and computational tractibility given the complex many-body phenomena. We use both on-lattice and off-lattice models and the primary computational tools are Monte Carlo simulations with sampling techniques and ensembles appropriate to specific situations. Our modeling approach is the first to capture silica polymerization, nanopore crystallization, and mesopore formation through computer-simulated self assembly.

36 MATERIALS SCIENCE↗

Evaluation of kriging based surrogate models constructed from mesoscale computations of shock interaction with particles

Macro-scale computations of shocked particulate flows require closure laws that model the exchange of momentum/energy between the fluid and particle phases. Closure laws are constructed in this work in the form of surrogate models derived from highly resolved mesoscale computations of shock-particle interactions. The mesoscale computations are performed to calculate the drag force on a cluster of particles for different values of Mach Number and particle volume fraction. Two Kriging-based methods, viz. the Dynamic Kriging Method (DKG) and the Modified Bayesian Kriging Method (MBKG) are evaluated for their ability to construct surrogate models with sparse data; i.e. using the least number of mesoscale simulations. It is shown that if the input data is noise-free, the DKG method converges monotonically; convergence is less robust in the presence of noise. The MBKG method converges monotonically even with noisy input data and is therefore more suitable for surrogate model construction from numerical experiments. This work is the first step towards a full multiscale modeling of interaction of shocked particle laden flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Concrete in engineered barriers for radioactive waste disposal facilities - phenomenological study and assessment of long term performance

The paper aims to highlight recent developments at the Belgian Nuclear Research Centre SCK-CEN in experimental and numerical study of the coupled physical-chemical behaviour of concrete subject to chemical degradation. The discussion mainly focusses on three inter-linked research projects covering novel experimental methods to study the alteration of hydraulic and transport properties during carbonation and calcium leaching, a pore scale numerical model to capture microstructural changes due to the above degradation processes and a generic multiscale model aimed at determining evolution of the properties of a macrostructure over the long term. The paper also describes supplementary continuum scale numerical studies concerning concrete-clay interactions and geochemical impact on the physical structure of concrete. Preliminary findings from these studies show encouraging results such as the development of novel leaching, water permeability and diffusion apparatus, a robust pore scale model based on Lattice-Boltzmann method and a mesoscale study focused on the importance of interfacial transition zones on the effective diffusivity for linear and nonlinear diffusion problems. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Multi-Scale Modeling Framework for Mercury Biogeochemistry

Multi-Scale modeling of mercury (Hg) geochemical speciation and reactions has been performed by integrating atomistic quantum chemical calculations with continuum scale speciation models. Major progress has been made in the improvement of quantum chemical models to calculate critical thermodynamic data for Hg complexes in aquatic environments. Rapid and reliable quantum chemical approaches have been developed for calculating acid dissociation constants (pK a ) and stability constants (log K), with calculated mean unsigned errors of 0.5 and 1.5 log units, respectively for ligand molecules and Hg complexes. At the continuum scale, systematic analysis of uncertainty propagation in mercury (Hg) speciation modeling has been conducted and was used to identify environmental conditions under which thermodynamic constant uncertainties are significant and recommended to be accounted for. The integrated framework for multi-scale modeling of mercury geochemistry is open to the research community through the web-based multiscale modeling aqueous speciation resource, AQUA-MER. The improved quantum chemical approaches for thermodynamic constant calculations are accessible through AQUA-MER and can be used to provide the missing constants in the continuum scale speciation calculations. In addition to low molecular mass Hg complex speciation, modeling natural aquatic environments also involve the transport of high molecular weight dissolved organic matter (DOM) in reactive flows simultaneously with equilibrium and kinetic reactions. To this end, atomistic MD simulations were performed to capture the details of aggregation, mechanisms and distribution of functional groups in DOM at the molecular level. The elemental composition and calculated bulk properties of the DOM models are in close agreement with experimental measurements. A travel-time based reactive transport model in the hyporheic zone of stream corridors was established for the multicomponent Hg-DOM-S system and implemented through PFLOTRAN.

54 ENVIRONMENTAL SCIENCES↗

Leveraging MARVEL and SPHERE to Demonstrate NEAMS Thermal Hydraulics Codes

The Single Primary Heat Extraction and Removal Emulator (SPHERE) and Microreactor Applications Research Validation and Evaluation (MARVEL) were utilized to demonstrate thermal hydraulics codes of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. A recent long-duration test of a sodium heat pipe conducted at the SPHERE facility was utilized for the validation of the heat pipe code Sockeye. The first 10 hours of operation featured a startup of the heat pipe from room temperature; this portion was modeled with Sockeye and compared to the experimental data. Measurements included several thermocouple measurements across the length of the heat pipe surface. Two of Sockeye’s heat pipe models were applied to experiment, as well as several models of the cooling at the condenser end of the heat pipe. Sockeye results gave a decent match to the experimental data, with notable discrepancy in the condenser section, where experimental data indicated a significant inactive length of the heat pipe. The MARVEL microreactor was modeled through a coupled simulation of System Analysis Module (SAM) and SubChannel Module (SCM). This multiscale coupled simulation was performed as an exercise to demonstrate the compatibility and applicability of these codes for thermal-hydraulic analysis of microreactors. The coupling between the SCM and SAM codes for multiscale modeling of MARVEL, was achieved through a domain overlapping approach. In this example a transient SAM simulation provided boundary conditions to SCM, and SCM provided pressure drop to SAM.

97 - MATHEMATICS AND COMPUTING↗

MultiSector Dynamics: 2023 Inaugural Workshop Report

Preface The MultiSector Dynamics (MSD) Community of Practice (CoP) hosted an inaugural workshop on October 3-5, 2023 at the University of California, Davis, to bring together members of the MSD community of practice to advance understanding of the co-evolution of human and natural systems, and to build the next generation of tools that bridge sectors, scales, and systems to realize a more resilient and equitable future. The theme of the workshop was "Advancing Complex Adaptive Human-Earth Systems Science in a World of Interconnected Risks". This document outlines the motivation for the workshop, its goals and objectives, the application process, the agenda, overviews of the training sessions offered to the workshop participants and a summary of each breakout session. The MSD workshop report further discusses the feedback from workshop participants and presents some reflections and next steps. The MSD Workshop organizers thank the DOE Office of Science, Earth and Environmental System Modeling, MultiSector Dynamics program area for financial support of its activities through the Integrated Multisector Multiscale Modeling (IM3) project. For more information related to the broader DOE MultiSector Dynamics Program please see https://climatemodeling.science.energy.gov/program-area/multisector-dynamics. D.L.M. and C.M.B. acknowledge support from the Laboratory Directed Research and Development Program of Oak Ridge National Laboratory (ORNL), managed by UT-Battelle, LLC, for the US Department of Energy (DOE). Disclaimer This report was prepared as an account of work sponsored by an agency of the United States Government. Neither theUnited States Government nor any agency thereof, nor Battelle Memorial Institute, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, complete- ness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial products, process, or service by trade name,trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof, or Battelle Memorial Institute. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof. Pacific Northwest National Laboratory operated by Battelle for the United States Department of Energy Available from:Office of Scientific and Technical Information http://www.OSTI.gov multisectordynamics.org This work is made available under the terms of the Creative Commons Attribution- NonCommercial 4.0 International (CC BY-NC 4.0) https://creativecommons.org/licenses/by-nc/4 Suggested citation: Monier, E., Reed, P.M., Vernon, C.R., Hadjimichael, A., Brelsford, C.M., Burleyson, C.B., Dyreson, A.R., Fletcher, S.M., Giang, A., Gupta, R.S., Jackson, N.D., Jones, A.D., Lamontagne, J.R., McCollum, D.L., Morris, J.F., Moss, R.H., Peng, W., Saari, R.K., Srikrishnan, V., Szinai, J.K., Yoon, J. (2024) MultiSector Dynamics: 2023 Inaugural Workshop Report. MSD-LIVE Data Repository. doi:10.57931/2371710.

Monier, Erwan↗

Prediction of local concentration fields in porous media with chemical reaction using a multi scale convolutional neural network

The study of solute transport in porous media is of interest in many chemical engineering systems. Some example applications include packed bed catalytic reactors, filtration devices, and batteries. The pore scale modeling of these systems is time consuming and may require large computing resources, for this reason computational fluid dynamics (CFD) simulations are not practical if a large number of simulations is required, like in multiscale modeling, where a model at a large scale calls for pore scale simulations. It has been shown that neural networks can be trained with a dataset of flow simulations and then predict fields orders of magnitude faster, and with less computational resources, in new domains. However, it is crucial to provide the neural network with an effective description of the domain and the undergoing operating conditions to be able to train models that generalize accurately in unseen samples. Therefore, research is needed to employ neural networks in new complex systems. The appropriate training of a network for predicting coupled flow and solute transport processes is an outstanding problem due to the complex interplay between geometry and operating conditions. In this work, we train a multi scale convolutional neural network (MSNet) with a diverse dataset of simulations of transport and chemical reaction in porous media to predict the local concentration fields in images of porous media. Our dataset contains a wide diversity of sphere pack arrangements under different operating conditions (Péclet and Reynolds numbers). Further, we train a robust model by employing different input descriptors that represent the medium and the different operating conditions of each system. Our trained model is able to provide nearly instantaneous predictions, compared to around twenty hours of the CFD workflow, with less than 3.5% error on new geometries and transport conditions. Thus the model could be easily integrated in a multiscale workflow where fast response is needed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

2040 Vision: A Roadmap for Integrated, Multiscale Materials and System Modeling and Simulation

Over the last few decades, advances in high-performance computing, new materials characterization methods, and, more recently, an emphasis on integrated computational materials engineering (ICME) and additive manufacturing have been a catalyst for multiscale modeling and simulation-based design of materials and structures in the aerospace industry. While these advances have driven significant progress in the development of aerospace components and systems, that progress has been limited by persistent technology and infrastructure challenges that must be overcome to realize the full potential of integrated materials and systems design and simulation modeling throughout the supply chain. As a result, NASA's Transformational Tools and Technology (TTT) Project sponsored a study (performed by a team led by Pratt & Whitney) to define the potential 25-year future state required for integrated multiscale modeling of materials and systems (e.g., load-bearing structures) to accelerate the pace and reduce the expense of innovation in future aerospace and aeronautical systems. This talk will describe the findings of this 2040 Vision study (e.g., the 2040 vision state; the required interdependent core technical work areas, Key Element (KE); associated critical gaps and actions to close those gaps; and major recommendations). This roadmap is a community consensus document and is a result of over 450 professionals input obtain via: i) four society workshops (AIAA, NAFEMS, and two TMS), ii) community-wide survey, and iii) the establishment of 9 expert panels (one per key element) consisting on average 10 non-team members from academia, government and industry to review, update content, and prioritize gaps and actions. It is NASA's hope and desire that this roadmap provides the strategic guidance to both public and private research and development decision makers to make the proposed 2040 vision state a reality and thereby provide a significant advancement in the United States global competitiveness.

Multiscale Modeling↗

A Method for Combining Experimentation and Molecular Dynamics Simulation to Improve Cohesive Zone Models for Metallic Microstructures

Fracture processes within a material begin at the nanometer length scale at which the formation, propagation, and interaction of fundamental damage mechanisms occur. Physics-based modeling of these atomic processes quickly becomes computationally intractable as the system size increases. Thus, a multiscale modeling method, based on the aggregation of fundamental damage processes occurring at the nanoscale within a cohesive zone model, is under development and will enable computationally feasible and physically meaningful microscale fracture simulation in polycrystalline metals. This method employs atomistic simulation to provide an optimization loop with an initial prediction of a cohesive zone model (CZM). This initial CZM is then applied at the crack front region within a finite element model. The optimization procedure iterates upon the CZM until the finite element model acceptably reproduces the near-crack-front displacement fields obtained from experimental observation. With this approach, a comparison can be made between the original CZM predicted by atomistic simulation and the converged CZM that is based on experimental observation. Comparison of the two CZMs gives insight into how atomistic simulation scales.

Hochhalter, J. D.↗

Rapid freezing of water under dynamic compression

Understanding the behavior of materials at extreme pressures is a central issue in fields like aerodynamics, astronomy, and geology, as well as for advancing technological grand challenges such as inertial confinement fusion. Dynamic compression experiments to probe high-pressure states often encounter rapid phase transitions that may cause the materials to behave in unexpected ways, and understanding the kinetics of these phase transitions remains an area of great interest. We examine experimental and theoretical/computational efforts to study the freezing kinetics of water to a high-pressure solid phase known as ice VII. We first present a detailed analysis of dynamic compression experiments in which water has been observed to freeze on sub-microsecond time scales to ice VII. This is followed by a discussion of the limitations of currently available molecular and continuum simulation methods in modeling these experiments. We then describe how our phase transition kinetics models, which are based on classical nucleation theory, provide a more physics-based framework that overcomes some of these limitations. We give suggestions on future experimental and modeling work on the liquid–ice VII transition, including an outline of the development of a predictive multiscale model in which molecular and continuum simulations are intimately coupled.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Evaluation of present-day extreme precipitation over the United States: an inter-comparison of convection and dynamic permitting configurations of E3SMv1

Abstract Accurate simulation of the present-day characteristics of mean and extreme precipitation at regional scales remains a challenge for Earth system models, which is due in part to deficiencies in model physics such as convective parameterization (CP), and coarse resolution. High horizontal resolution (HR, ∼25 km) and multiscale modeling framework (MMF, i.e. replacing conventional CP with embedded km-scale cloud-resolving models) are two promising directions that could help improve the interaction between subgrid-scale physical processes and large-scale climate. Here, we evaluate simulated extreme precipitation over the United States (US) across three configurations (i.e. low-resolution [LR], HR, and MMF) of the Energy Exascale Earth System Model (E3SMv1) and intercompare them against two gridded observation datasets (climate prediction center daily US precipitation and integrated multi-satellite retrievals for global precipitation measurement). We assess the model’s ability to simulate very heavy seasonal precipitation (illustrated by the difference between the 99th and 90th percentile values) as well as the spatial distributions of several extreme precipitation indices defined by the expert team on climate change detection and indices. Our results show that both the dry (i.e. consecutive dry days (CDD)) and wet (i.e. consecutive wet days, maximum 5 day precipitation, and very wet days) extremes evaluated herein show some improvement as well as degradation with MMF and HR relative to LR. These results vary across seasons and US subregions. For instance, only the very heavy precipitation of winter is improved with MMF and HR. Both configurations alleviate the well-known drizzling bias evident in LR across both winter and summer in many parts of the US, largely due to the overall improvement in intensity and frequency of precipitation. Additionally, our results suggest that while E3SMv1-MMF has higher intensity rates when it does rain, it has too many CDD during the summer, contributing to a low mean precipitation bias.

54 ENVIRONMENTAL SCIENCES↗

A Hybrid Multi-Scale Model of Crystal Plasticity for Handling Stress Concentrations

Microstructural effects become important at regions of stress concentrators such as notches, cracks and contact surfaces. A multiscale model is presented that efficiently captures microstructural details at such critical regions. The approach is based on a multiresolution mesh that includes an explicit microstructure representation at critical regions where stresses are localized. At regions farther away from the stress concentration, a reduced order model that statistically captures the effect of the microstructure is employed. The statistical model is based on a finite element representation of the orientation distribution function (ODF). As an illustrative example, we have applied the multiscaling method to compute the stress intensity factor K I around the crack tip in a wedge-opening load specimen. The approach is verified with an analytical solution within linear elasticity approximation and is then extended to allow modeling of microstructural effects on crack tip plasticity.

36 MATERIALS SCIENCE↗