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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

TPS Certification by Analysis: Model-Driven Characterization of Properties and Failure in Woven Thermal Protection Systems

Woven, ablative thermal protection system (TPS) materials provide a robust option for aggressive (re)entries and thus have been baselined for the upcoming Mars Sample Return (MSR) mission’s Earth Entry System (EES). The reliability requirements for MSR-EES necessitate understanding of material property variability, which could be significant given the complex structure and anisotropic nature of properties in TPS weaves, as well as the response to potential impact with micrometeoroids or orbital debris during the EES re-entry. The TPS Certification by Analysis effort within the Entry Systems Modeling project seeks to provide computational models and analyses that support the certification against such material-based risks. For the present talk, focus will be given to the characterization of baseline woven TPS material properties and mechanical failure limits, which entails (1) use of computational techniques (e.g., machine learning) to interpret computed tomography images of the weave to generate representative structural models and (2) application of multiscale material modeling approaches to characterize thermomechanical and failure properties.

Justin B. Haskins↗
DOE OSTI2021

A review of microwave-assisted process intensified multiphase reactors

Microwaves provide alternative heating and allow process intensification due to their rapid, volumetric, and selective nature. Recognizing the central role of multiphase reactors in chemical industry, a recent surge in employing microwaves is observed. We review the recent experimental and modeling investigations of microwave heating of multiphase reactors with emphasis on chemical engineering applications. Here we demonstrate that there is accumulated evidence for improved performance via microwave heating and a clear opportunity for further process intensification. In most of the cases, this improved performance stems from a temperature gradient between two phases. We discuss the ongoing debate on the mechanism by which microwaves affect chemical processes exacerbated by the inability of measuring the temperature distribution in a microwave cavity. We outline recent progress in this direction in monolith reactors and needs for future work. We underscore the lack of detailed modeling and simulation tools even for single-phase systems and emphasize the imperative for multiscale predictive modeling to bridge the experimental-modeling gap. Promising results are shown by a few recently published modeling studies that can predict the experimental measurements in complex multiphase reactors. A combination of experimental and modeling tools can provide a comprehensive picture of the microwave multiphase reactors as well as a means toward scale-up and optimization.

42 ENGINEERING↗
DOE OSTI2018

Testing Dynamic Earthquake Rupture Models Generated With Stochastic Stress Drop

In this study we investigate the performance of rupture dynamics modeling using heterogenous stress for a M6.5 crustal earthquake on a strike slip vertical fault. In our simulations we used heterogenous stress drop models with correlated stochastic perturbations following the work of Pitarka et al. (2009). The stress model with multiscale variability is designed to accommodate high-stress drop asperity areas and small-scale spatial stress variations. The stress drop and size of asperities are estimated following the IM2011 (Irikura and Miyake, 2011) recipe. In our spontaneous rupture modeling we use a linear slip weakening friction law, with depth depend strength excess, dynamic stress drop and slip weakening distance Dc. Our strategy in preparing the model parameters with vertical and lateral variability is in agreement with observations of earthquake rupture that have confirmed the existence of a shallow weak rupture zone consisting of ductile crustal material, located in the upper 4 km of the crust, and brittle seismogenic zone, underlying the weak zone. The primary objective of our study is to develop a procedure for preparing dynamic rupture models with multiscale variability, and validate spontaneous rupture modeling of strike-slip faults by comparing the resulting rupture kinematics and near-fault ground motion with empirical models for this type of rupture.

58 GEOSCIENCES↗
DOE OSTI2017

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

Virtual Model Validation of Complex Multiscale Systems: Applications to Nonlinear Elastostatics

We propose a virtual statistical validation process as an aid to the design of experiments for the validation of phenomenological models of the behavior of material bodies, with focus on those cases in which knowledge of the fabrication process used to manufacture the body can provide information on the micro-molecular-scale properties underlying macroscale behavior. One example is given by models of elastomeric solids fabricated using polymerization processes. We describe a framework for model validation that involves Bayesian updates of parameters in statistical calibration and validation phases. The process enables the quanti cation of uncertainty in quantities of interest (QoIs) and the determination of model consistency using tools of statistical information theory. We assert that microscale information drawn from molecular models of the fabrication of the body provides a valuable source of prior information on parameters as well as a means for estimating model bias and designing virtual validation experiments to provide information gain over calibration posteriors.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Modeling of grain size strengthening in tantalum at high pressures and strain rates

Laser-driven ramp wave compression experiments have been used to investigate the strength (flow stress) of tantalum and other metals at high pressures and high strain rates. Recently this kind of experiment has been used to assess the dependence of the strength on the average grain size of the material, finding no detectable variation with grain size. The insensitivity to grain size has been understood theoretically to result from the dominant effect of the high dislocation density generated at the extremely high strain rates of the experiment. Here we review the experiments and describe in detail the multiscale strength model used to simulate them. The multiscale strength model has been extended to include the effect of geometrically necessary dislocations generated at the grain boundaries during compatible plastic flow in the polycrystalline metal. Lastly, we use the extended model to make predictions of the threshold strain rates and grain sizes below which grain size strengthening would be observed in the laser-driven Rayleigh-Taylor experiments.

36 MATERIALS SCIENCE↗
NASA NTRS2009

Hydrogen Cyanide in the Upper Troposphere: GEM-AQ Simulation and Comparison with ACE-FTS Observations

We investigate the spatial and temporal distribution of hydrogen cyanide (HCN) in the upper troposphere through numerical simulations and comparison with observations from a space-based instrument. To perform the simulations, we used the Global Environmental Multiscale Air Quality model (GEM-AQ), which is based on the threedimensional Gobal multiscale model developed by the Meteorological Service of Canada for operational weather forecasting. The model was run for the period 2004-2006 on a 1.5deg x 1.5deg global grid with 28 hybrid vertical levels from the surface up to 10 hPa. Objective analysis data from the Canadian Meteorological Centre were used to update the meteorological fields every 24 h. Fire emission fluxes of gas species were generated by using year-specific inventories of carbon emissions with 8-day temporal resolution from the Global Fire Emission Database (GFED) version 2. The model output is compared with HCN profiles measured by the Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) instrument onboard the Canadian SCISAT-1 satellite. High values of up to a few ppbv are observed in the tropics in the Southern Hemisphere; the enhancement in HCN volume mixing ratios in the upper troposphere is most prominent in October. Low upper-tropospheric mixing ratios of less than 100 pptv are mostly recorded at middle and high latitudes in the Southern Hemisphere in May-July. Mixing ratios in Northern Hemisphere peak in the boreal summer. The amplitude of the seasonal variation is less pronounced than in the Southern Hemisphere. The comparison with the satellite data shows that in the upper troposphere GEM-AQ perform7s well globally for all seasons, except at northern hi gh and middle latitudes in surnmer, where the model has a large negative bias, and in the tropics in winter and spring, where it exhibits large positive bias. This may reflect inaccurate emissions or possible inaccuracies in the emission profile. The model is able to explain most of the observed variability in the upper troposphere HCN field, includin g the interannual variations in the observed mixing ratio. A complementary comparison with daily total columns of HCN from two middle latitude ground-based stations in Northern Japan for the same simulation period shows that the model captures the observed seasonal variation and also points to an underestimation of model emissions in the Northern Hemisphere in the summer. The estimated average global emission equals 1.3 Tg N/yr. The average atmospheric burden is 0.53 Tg N, and the corresponding lifetime is 4.9 months.

Lupu, A.↗
DOE OSTI2026

Simulations of grain growth in tungsten armor materials under ARC plasma edge operation conditions using an integrated plasma-edge/materials model

An integrated model of grain growth deuterium-exposed tungsten polycrystals, consisting of a two-dimensional vertex dynamics model fitted to atomistic data, has been developed to assess the grain growth kinetics of deuterium-exposed polycrystalline tungsten (W). The model tracks the motion of grain boundaries under the effect of driving forces stemming from grain boundary curvature and differential deuterium concentration accumulation. Here, we apply the model to experimentally synthesized tungsten polycrystals under deuterium-saturated conditions relevant to the ARC concept design. The results indicate rapid grain growth kinetics in the near-surface region adjacent to the plasma, where the temperature reaches 1400 K, whereas the microstructure remains stable deeper in the material with the lower temperature of 1000 K. The combined modeling and analysis further reveal that monolithic tungsten produced via conventional fabrication routes is highly susceptible to grain coarsening at temperatures exceeding 1000 K, largely independent of the magnitude of the applied driving force. Moreover, the accumulation of deuterium near grain boundaries has a pronounced inhibitory effect on grain boundary migration. High-angle grain boundaries ( > 50°) contribute more significantly to the overall grain growth process.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Seamlessly joining length scales: From atomistic thermal graphs to anisotropic continuum conductivity

Thermal transport in complex solids is governed by local structure, defects, and anisotropy, yet most continuum models still rely on oversimplified and homogenized conductivities. Here, we bridge atomistic and continuum descriptions by building finite element (FE) models directly from the site-projected thermal conductivity (SPTC), an atomic-level decomposition of the Green–Kubo thermal conductivity. We introduce a toolkit, the “Simulator Collection for Atomic-to-Continuum Scales (SCACS)”, which uses a graph neural network to predict SPTC on large atomic structures, coarse-grains these fields into anisotropic conductivity tensors, and embeds them into the heat-flow FE equation with a customized, anisotropy-aware adaptive mesh refinement scheme. Applied to silicon nanostructures, the resulting FE models act as representative volume elements, reproduce bulk conductivities, and capture interfacial and defect-driven anisotropy while maintaining thermodynamic consistency. Additionally, SCACS predicts experimental conductance trends and fields. This work demonstrates a general route for transferring atomistic transport information into device-scale thermal simulations with physics-based approximations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

A Stochastic Reduced-Order Model for Statistical Microstructure Descriptors Evolution

Integrated computational materials engineering (ICME) models have been a crucial building block for modern materials development, relieving heavy reliance on experiments and significantly accelerating the materials design process. However, ICME models are also computationally expensive, particularly with respect to time integration for dynamics, which hinders the ability to study statistical ensembles and thermodynamic properties of large systems for long time scales. To alleviate the computational bottleneck, we propose to model the evolution of statistical microstructure descriptors as a continuous-time stochastic process using a non-linear Langevin equation, where the probability density function (PDF) of the statistical microstructure descriptors, which are also the quantities of interests (QoIs), is modeled by the Fokker–Planck equation. In this work, we discuss how to calibrate the drift and diffusion terms of the Fokker–Planck equation from the theoretical and computational perspectives. The calibrated Fokker–Planck equation can be used as a stochastic reduced-order model to simulate the microstructure evolution of statistical microstructure descriptors PDF. Considering statistical microstructure descriptors in the microstructure evolution as QoIs, we demonstrate our proposed methodology in three integrated computational materials engineering (ICME) models: kinetic Monte Carlo, phase field, and molecular dynamics simulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Model America - Arizona extract from ORNL's AutoBEM v1.1

Oak Ridge National Laboratory (ORNL) has developed the Automatic Building Energy Modeling (AutoBEM) software suite to process multiple types of data, extract building-specific descriptors, generate building energy models, and simulate them on High Performance Computing (HPC) resources. For more information, see AutoBEM-related publications (bit.ly/AutoBEM).Two sets of sample data are provided for 2,555,152 buildings located within the boundary of Arizona in the United States:Data (846.3MB *.csv) - minimalist list of each building (rows) for the following fields (columns) • ID - unique building ID • Centroid - building center location in latitude/longitude (from Footprint2D) • Footprint2D - building polygon of 2D footprint (lat1/lon1_lat2/lon2_...) • State_abbr - state name • Area - estimate of total conditioned floor area (ft2) • Area2D - footprint area (ft2) • Height - building height (ft) • NumFloors - number of floors (above-grade) • WWR_surfaces - percent of each facade (pair of points from Footprint2D) covered by fenestration/windows (average 14.5% for residential, 40% for commercial buildings) • CZ - ASHRAE Climate Zone designation • BuildingType - DOE prototype building designation (IECC=residential) as implemented by OpenStudio-standards • Standard - building vintage • Sample Models (114GB*.zip by county) - OpenStudio and EnergyPlus building energy models named according to IDThis data is made free and openly available in hopes of stimulating any simulation-informed use case. Data is provided as-is with no warranties, express or implied, regarding fitness for a particular purpose. We wish to thank our sponsors which include Oak Ridge National Laboratory (ORNL), U.S. Dept. of Energy’s (DOE) Building Technologies Office (BTO), Office of Electricity (OE), and Biological and Environmental Research (BER).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Use of Remote Sensing and In-situ Observations to Develop and Evaluate Improved Representations of Convection and Clouds for the Accelerated Climate Model for Energy. Final report

Representation of unresolved cloud processes has been one of the most challenging problems in numerical modeling of the atmosphere. During the last decades, considerable efforts have been spent to better represent the collective effects of those processes in terms of resolved scales (see Arakawa 2004 for a review). Still, the progress on this problem has been unacceptably slow as pointed out by Randall et al. (2003). The large part of the difficulty is inherent in the artificial separation of atmospheric processes into the explicitly simulated resolved scales and the highly parameterized subgrid scales. Arakawa (2004) and Arakawa et al. (2011) suggested two plausible approaches to eliminate the problem of artificial separation of atmospheric processes: “parameterize everything” and “resolve everything”.

58 GEOSCIENCES↗
NASA NTRS2021

Mineral Dust Cycle in the Multiscale Online Nonhydrostatic AtmospheRe CHemistry Model (Monarch) Version 2.0

We present the dust module in the Multiscale Online Non-hydrostatic AtmospheRe CHemistry model (MONARCH) version 2.0, a chemical weather prediction system that can be used for regional and global modeling at a range of resolutions. The representations of dust processes in MONARCH were upgraded with a focus on dust emission (emission parameterizations, entrainment thresholds, considerations of soil moisture and surface cover), lower boundary conditions (roughness, potential dust sources), and dust–radiation interactions. MONARCH now allows modeling of global and regional mineral dust cycles using fundamentally different paradigms, ranging from strongly simplified to physics-based parameterizations. We present a detailed description of these updates along with four global benchmark simulations, which use conceptually different dust emission parameterizations, and we evaluate the simulations against observations of dust optical depth. We determine key dust parameters, such as global annual emission/deposition flux, dust loading, dust optical depth, mass-extinction efficiency, single-scattering albedo, and direct radiative effects. For dust-particle diameters up to 20 µm, the total annual dust emission and deposition fluxes obtained with our four experiments range between about 3500 and 6000 Tg, which largely depend upon differences in the emitted size distribution. Considering ellipsoidal particle shapes and dust refractive indices that account for size-resolved mineralogy, we estimate the global total (longwave and shortwave) dust direct radiative effect (DRE) at the surface to range between about −0.90 and −0.63 W m−2 and at the top of the atmosphere between −0.20 and −0.28 W m−2. Our evaluation demonstrates that MONARCH is able to reproduce key features of the spatiotemporal variability of the global dust cycle with important and insightful differences between the different configurations.

Dust cycle↗
DOE OSTI2025

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

A Novel Multiphysics Multiscale Multiporosity Shale Gas Transport Model for Geomechanics/Flow Coupling in Steady and Transient States

Summary A novel multiphysics multiscale multiporosity shale gas transport (M3ST) model was developed to investigate shale gas transport in both transient and steady states. The microscale model component contains a kerogen domain and an inorganic matrix domain, and each domain has its own geomechanical and gas transport properties. Permeabilities of various shale cores were measured in the laboratory using a pulse decay permeameter (PDP) with different pore pressure and confining stress combinations. The PDP-measured apparent permeability as a function of pore pressure under two effective stresses was fitted using the microscale M3ST model component based on nonlinear least squares fitting (NLSF), and the fitted model parameters were able to provide accurate model predictions for another effective stress. The parameters and petrophysical properties determined in the steady state were then used in the transient-state, continuum-scale M3ST model component, which performed history matching of the evolutions of the upstream and downstream gas pressures. In addition, a double-exponential empirical model was developed as a powerful alternative to the M3ST model to fit laboratory-measured apparent permeability under various effective stresses and pore pressures. The developed M3ST model and the research findings in this study provided critical insights into the role of the multiphysics mechanisms, including geomechanics, fluid dynamics and transport, and the Klinkenberg effect on shale gas transport across different spatial scales in both steady and transient states.

Engineering↗
DOE OSTI2015

Building Virtual Ecosystems: Computational Challenges for Mechanistic Modeling of Terrestrial Environments Workshop, Germantown, Maryland, March 26–27, 2014

The mission of the Climate and Environmental Sciences Division of the Office of Biological and Environmental Research (BER) within the U.S. Department of Energy’s (DOE) Office of Science is “to advance a robust, predictive understanding of Earth’s climate and environmental systems and to inform the development of sustainable solutions to the nation’s energy and environmental challenges.” This formidable challenge requires quantification of stocks and controls on states, fluxes, and residence times of water, carbon, and other key elements through all components of the terrestrial system. These components include vegetation, soils, the deep vadose zone, groundwater, and surface water. To achieve this level of predictive understanding, a new generation of multiscale, multiphysics models is needed for terrestrial systems—models that incorporate process couplings and feedbacks between various “pools” (i.e., vegetation, soils, subsurface aquifers, and surface waters) across wide ranges of spatial and temporal scales. To explore the potential of a new generation of multiscale, multiphysics models for revolutionizing the understanding of terrestrial ecosystem dynamics, BER held the Computational Challenges for Mechanistic Modeling of Terrestrial Environments workshop on March 26–27, 2014, in Germantown, Maryland. The workshop brought together 29 researchers with diverse expertise, including hydrologists, environmental scientists, ecologists, microbiologists, plant scientists, and computational scientists. Through a combination of invited talks, breakout sessions, and report-back discussions, workshop attendees identified challenges and research opportunities needed to develop a more seamless and continuous framework for mechanistic modeling of terrestrial environments extending from the bedrock to the atmospheric boundary layer and from single-plant systems to fields of crops to watersheds and river basins.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

Parameterizing the Morse potential for coarse-grained modeling of blood plasma

Multiscale simulations of fluids such as blood represent a major computational challenge of coupling the disparate spatiotemporal scales between molecular and macroscopic transport phenomena characterizing such complex fluids. In this paper, a coarse-grained (CG) particle model is developed for simulating blood flow by modifying the Morse potential, traditionally used in Molecular Dynamics for modeling vibrating structures. The modified Morse potential is parameterized with effective mass scales for reproducing blood viscous flow properties, including density, pressure, viscosity, compressibility and characteristic flow dynamics of human blood plasma fluid. The parameterization follows a standard inverse-problem approach in which the optimal micro parameters are systematically searched, by gradually decoupling loosely correlated parameter spaces, to match the macro physical quantities of viscous blood flow. The predictions of this particle based multiscale model compare favorably to classic viscous flow solutions such as Counter-Poiseuille and Couette flows. It demonstrates that such coarse grained particle model can be applied to replicate the dynamics of viscous blood flow, with the advantage of bridging the gap between macroscopic flow scales and the cellular scales characterizing blood flow that continuum based models fail to handle adequately.

60 APPLIED LIFE SCIENCES↗