Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiscale models”

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 559 records · Page 31

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↗

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↗

Bemnifosbuvir: An HCV NS5B Inhibitor With Multiple Modes of Action

Bemnifosbuvir (BEM) is a potent, pan-genotypic inhibitor targeting the hepatitis C virus (HCV) NS5B polymerase. Its antiviral activity was evaluated in an ascending dose phase I clinical trial involving 30 patients treated once a day for 7 days. After treatment initiation, plasma HCV RNA declined in a biphasic manner with a mean reduction of 2.3 log 10 IU/mL after 24 hours and 4.4 log10 IU/mL by day 7 for the highest dose. Alanine aminotransferase (ALT) also normalized in most patients. We employed a multiscale mathematical model fitted to the HCV RNA and ALT dynamics to quantify the antiviral activity and evaluate the modes of action of BEM. We found that models in which BEM only acted as a typical HCV RNA polymerase inhibitor and reduced the intracellular production of HCV RNA did not fit the data as well as models in which BEM had multiple modes of action, including suppressing viral assembly and secretion and enhancing intracellular HCV RNA degradation. BEM's effectiveness in inhibiting intracellular HCV RNA production increased with dose (150 mg/day: 88.2%, 300 mg/day: 98.8%, 600 mg/day: 99.5%), while inhibition of viral assembly and release was ~95% effective regardless of dose. We observed a dose-dependent enhancement in the degradation of intracellular HCV RNA, with degradation rates 1.5-fold higher in patients receiving 300 mg/day and 2.7-fold higher in those receiving 600 mg/day than in patients receiving 150 mg/day. No significant differences in antiviral activity were detected between HCV genotypes 1b and 3 or between patients with and without compensated cirrhosis.

59 BASIC BIOLOGICAL SCIENCES↗

Data-driven learning of nonlocal physics from high-fidelity synthetic data

A key challenge to nonlocal models is the analytical complexity of deriving them from first principles, and frequently their use is justified a posteriori. Here, we extract nonlocal models from data, circumventing these challenges and providing data-driven justification for the resulting model form. Extracting data-driven surrogates is a major challenge for machine learning (ML) approaches, due to nonlinearities and lack of convexity — it is particularly challenging to extract surrogates which are provably well-posed and numerically stable. Our scheme not only yields a convex optimization problem, but also allows extraction of nonlocal models whose kernels may be partially negative while maintaining well-posedness even in small-data regimes. To achieve this, based on established nonlocal theory, we embed in our algorithm sufficient conditions on the non-positive part of the kernel that guarantee well-posedness of the learnt operator. These conditions are imposed as inequality constraints to meet the requisite conditions of the nonlocal theory. We demonstrate this workflow for a range of applications, including reproduction of manufactured nonlocal kernels; numerical homogenization of Darcy flow associated with a heterogeneous periodic microstructure; nonlocal approximation to high-order local transport phenomena; and approximation of globally supported fractional diffusion operators by truncated kernels.

42 ENGINEERING↗

Multiscale and Machine Learning Modeling for Additive Manufacturing

Additive manufacturing (AM) techniques provide the opportunity to simultaneously design new materials and components with complex structures in less time, enabling faster material developments. Even though compositionally similar, the texture of the materials produced by such techniques is significantly different from conventionally manufactured materials. Additively manufactured materials produces highly heterogeneous microstructure within a single build. Such variations in the microstructure make qualifying AM products challenging for extreme environment applications. Understanding the AM process and its influence on the materials’ microstructures/properties is paramount for evaluating the workability and performance of the manufactured materials. The performance of AM materials for advanced nuclear reactor applications is of interest to the Advanced Materials and Manufacturing Technologies (AMMT) program under the Department of Energy Office of Nuclear Energy. Hence, considering the microstructural variabilities in the AM products and their impact on the performance of the material, it is important to correlate the process conditions to the final product and establish a process-structure-property- performance (PSPP) correlation for AM materials.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A Multiscale Nonlocal Progressive Damage Model for Composite Materials

In this paper, the advantages of a nonlocal progressive damage formulation are described and demonstrated. An approximation of the nonlocal formulation was implemented coupled with the MAT162 composite damage model as a User defined material model in the LS DYNA environment. A comparison of the local model and the nonlocal model is simulated for an 8-ply laminate under tension is carried for increasing mesh densities. The results show the regularization achieved by nonlocal models by providing mesh independent results.

Kodagali, Karan↗

Multi-scale fission product release model with comparison to AGR data

TRistructural ISOtropic (TRISO) particle fuel is central to several advanced, high-temperature reactor designs. Each particle consists of a fuel kernel encapsulated by three layers of carbon and ceramics that prevent the release of fission products and ensure physical integrity. Despite outstanding retention properties, fission product release has been observed from intact particles. To better understand and quantify fission product release from TRISO particles, a multiscale, mechanistic model of fission product transport is being developed by the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. Previous work focused on silver (Ag) transport and improved Ag release predictions. The work described in this report builds on this experience to better understand cesium (Cs) transport in silicon carbide (SiC), the main barrier to the release of fission products. Atomistic simulations provide bulk and grain boundary (GB) Cs diffusivities in SiC, which are used by phase field simulations in the mesoscale code Marmot to determine the temperature, microstructure, and irradiation-dependent Cs diffusivity at the mesoscale in SiC. This approach attributes the different temperature regimes experimentally observed for Cs diffusivities in SiC to a transition from bulk-dominated diffusivity at high temperatures to a GB-dominated regime at low temperatures, providing new insight. The multiscale, mechanistic effective diffusivity is then implemented in the fuel performance code BISON and further validated by comparing Cs release predictions from Advanced Gas Reactor (AGR)-1 and AGR-2 post-irradiation measurements. The new model improves BISON’s predictability. This document also reports improvements made on Ag transport modeling by accounting for different GB types having different diffusivities. Moreover, this report details preliminary efforts to model palladium (Pd) attack of the SiC at the mesoscale using a phase field approach. Pd attack and its impact on accelerated Ag transport remains a misunderstood phenomenon, and we use the model to demonstrate that the formation of lamellae that has been observed in experiments can be explained by the reaction of Pd with SiC to form alternating layers of graphite and Pd 2 Si. This effort aims to improve our understanding of the reaction and eventually provide a model for BISON to account for Pd penetration and its effects on fission product release.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Physics-Informed Machine Learning Model for Ceramic Matrix Composite Creep

A physics-informed recurrent neural network (RNN) based surrogate model is developed to emulate the nonlinear, time-dependent constitutive behavior of ceramic matrix composites (CMCs) driven by matrix damage and constituent creep at the microscale. Physics-informed constraints are introduced into the surrogate model through regularization to ground the prediction in physics and improve its predictive capabilities. Training data is generated using the high-fidelity generalized method of cells (HFGMC) approach which calls appropriate creep and damage models for each of the constituents. This coupling permits simulating the nonlinear behavior of CMCs based on constituent response at the microscale along with microstructural features such as fiber and porosity volume fraction and fiber radius. The microscale repeating unit cell is loaded under creep fatigue conditions to replicate the material loading experienced in a turbine engine. Therefore, the RNN-based surrogate model is tasked with predicting, as a function of variable input stress sequence, temperature, and microstructural features, the resulting strain history response while satisfying physical constraints related to creep rate, isochoric inelastic deformation, and strain energy density. The trained surrogate model is shown to effectively match the strain history over quantified distributions of microstructural features and relevant loading regimes and temperatures. Neural network based surrogate models can offer efficient alternatives to running computationally intensive multiscale material models to simulate the nonlinear response of large structural models. Therefore, the presented work provides evidence towards the feasibility of developing, training, and running such models for CMCs with complex microstructures, nonlinear time-dependent material response, and under non-monotonic loading conditions.

ceramic matrix composites↗

Description and Evaluation of the Multiscale Online Nonhydrostatic AtmospheRe CHemistry Model (NMMB-MONARCH) Version 1.0: Gas-Phase Chemistry at Global Scale

This paper presents a comprehensive description and benchmark evaluation of the tropospheric gas-phase chemistry component of the Multiscale Online Nonhydrostatic AtmospheRe CHemistry model (NMMBMONARCH), formerly known as NMMB/BSC-CTM, that can be run on both regional and global domains. Here, we provide an extensive evaluation of a global annual cycle simulation using a variety of background surface stations (EMEP, WDCGG and CASTNET), ozonesondes (WOUDC, CMD and SHADOZ), aircraft data (MOZAIC and several campaigns), and satellite observations (SCIAMACHY and MOPITT).We also include an extensive discussion of our results in comparison to other state-of-the-art models. We note that in this study, we omitted aerosol processes and some natural emissions (lightning and volcano emissions). The model shows a realistic oxidative capacity across the globe. The seasonal cycle for CO is fairly well represented at different locations (correlations around 0.3-0.7 in surface concentrations), although concentrations are underestimated in spring and winter in the Northern Hemisphere, and are overestimated throughout the year at 800 and 500 hPa in the Southern Hemisphere. Nitrogen species are well represented in almost all locations, particularly NO2 in Europe (root mean square error - RMSE - below 5 ppb). The modeled vertical distributions of NOx and HNO3 are in excellent agreement with the observed values and the spatial and seasonal trends of tropospheric NO2 columns correspond well to observations from SCIAMACHY, capturing the highly polluted areas and the biomass burning cycle throughout the year. Over Asia, the model underestimates NOx from March to August, probably due to an underestimation of NOx emissions in the region. Overall, the comparison of the modeled CO and NO2 with MOPITT and SCIAMACHY observations emphasizes the need for more accurate emission rates from anthropogenic and biomass burning sources (i.e., specification of temporal variability).

atmospheric chemistry↗

Regional-Scale Modeling Parameterizations for Secondary Organic Aerosol Formation from Isoprene Epoxydiols: Experimentally Based Evaluation and Optimization

Isoprene is an abundant volatile organic compound emitted from broadleaf forests. Under low nitric oxide concentrations, isoprene is photochemically oxidized to form gas-phase isoprene epoxydiols (IEPOX). In the presence of acidified sulfate aerosols, IEPOX enhances the secondary organic aerosol (SOA) formation. Predictions of IEPOX-SOA in regional-scale models, e.g., the Community Multiscale Air Quality Model (CMAQ), are uncertain due to homogeneous aerosol assumptions, underpredictions of water uptake (hygroscopicity), and aerosol surface area. Here, we used experimental measurements of IEPOX-SOA tracers, 2-methyltetrols (2-MT) and 2-methyltetrol sulfates (2-MTS), formed at initial IEPOX-to-inorganic sulfate ratios ranging from 1–10.5, at ∼50% relative humidity to constrain key IEPOX-SOA parameters: phase separation, organic shell diffusivity (D org ), acidity, hygroscopic growth, mass accommodation, and kinetics. The base CMAQ parametrization overpredicted experimental IEPOX-SOA with an average normalized mean bias (NMB average ) of 1.63. CMAQ with phase separation underpredicted IEPOX-SOA (NMB average = −0.71). Using the phase-separated model, CMAQ model performance was optimized (NMB average = 0.077) with an increased D org = 2 × 10 –16 m 2 s –1 and increased rate constants (k 2-MT = 1 × 10 –3 M 2 s –1 , k 2-MTS = 8.83 × 10 –3 M 2 s –1 ). The optimized model explicitly accounted for hygroscopic growth by utilizing experimentally derived growth rates, improving aerosol surface area predictions. Our model highlights the importance of the aerosol mixing state (homogeneous versus phase-separated), aerosol size dynamics, and hygroscopic growth in modeling heterogeneous reactive uptake of IEPOX.

aerosols↗

Arrested Phase Separation in Double-Exchange Models: Large-Scale Simulation Enabled by Machine Learning

In this work, we present large-scale dynamical simulations of electronic phase separation in the single-band double-exchange model based on deep-learning neural-network potentials trained from small-size exact diagonalization solutions. We uncover an intriguing correlation-induced freezing behavior as doped holes are segregated from half filled insulating background during equilibration. While the aggregation of holes is stabilized by the formation of ferromagnetic clusters through Hund’s coupling between charge carriers and local magnetic moments, this stabilization also creates confining potentials for holes when antiferromagnetic spin-spin correlation is well developed in the background. The dramatically reduced mobility of the self-trapped holes prematurely disrupts further growth of the ferromagnetic clusters, leading to an arrested phase separation. Implications of our findings for phase separation dynamics in materials that exhibit colossal magnetoresistance effect are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Modeling Images of Natural 3D Surfaces: Overview and Potential Applications

Generative models of natural images have long been used in computer vision. However, since they only describe the of 2D scenes, they fail to capture all the properties of the underlying 3D world. Even though such models are sufficient for many vision tasks a 3D scene model is when it comes to inferring a 3D object or its characteristics. In this paper, we present such a generative model, incorporating both a multiscale surface prior model for surface geometry and reflectance, and an image formation process model based on realistic rendering, the computation of the posterior model parameter densities, and on the critical aspects of the rendering. We also how to efficiently invert the model within a Bayesian framework. We present a few potential applications, such as asteroid modeling and Planetary topography recovery, illustrated by promising results on real images.

Jalobeanu, Andre↗

Materials Informatics at NASA GRC: Machine Learning Surrogate Modeling, Data Management, and Integrated Toolsets for Establishing/Maintaining the Digital Thread

Integrated Computational Materials Engineering (ICME) has recently received widespread attention due to its promises in reducing dependence on physical testing for engineering design by relying on simulation, reducing both time and cost to market for various applications. ICME however requires validated multiscale material models, which heavily depend on available test data with full material and test pedigree, including material processing, test and measurement equipment, raw data collection, and analysis methodology and results that is findable and usable, along with integrated, efficient toolsets for effectively passing information across various length and time scales across such models. At the NASA Glenn Research Center under the Transformational Tools and Technologies Project, significant recent efforts have been directed towards establishing the required cyberinfrastructure to enable optimized ICME processes and the design of “fit-for-purpose” materials to achieve the goals outlined in the NASA Vision 2040 report. Such efforts include development of multiscale physics-based material models, which can be used to train highly efficient surrogate machine learning models, development of best practices and infrastructure for effective, traceable materials information management, and development of toolsets that integrate with physics-based codes, machine learning models, and an information management system to enable high throughput of materials data collection and analysis, establishment of digital twins and the digital thread, and automation of the ICME design process for material optimization.

Machine Learning↗

Numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

For this work, we numerically evaluated photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively) on stream corridor transport. Correctly separating HSZ and SSZ effects is critical to estimating the hydro-biogeochemical function of a stream because HSZs and SSZs expose solutes to significantly different biogeochemical conditions, like sunlight exposure, microbial processes, and oxygen availability. Our numerical experiments used a multiscale river-corridor transport model implemented in the ATS code, which accommodates multiple storage zones with distinct travel time distributions and biogeochemical reactions. For parameter inferences, we used Bayesian inverse modeling. We found that breakthrough curves for photo-decaying tracers from day and night injection can delineate surface and hyporheic transient storage contributions, but only when interpreted jointly through a two-storage zone model. Numerical experiments that used only daytime injection or interpreted breakthrough curves with a single storage zone model yielded good fit to breakthrough curves, but parameter estimates were biased and controlling processes misattributed, examples of good model fits for the wrong reasons. Using those biased parameter estimates in reactive transport simulations resulted in significantly different projections of denitrification, which underscores the potential for stream function to be mischaracterized if tracer tests are interpreted through an inappropriately simplified model for transient storage. More generally, this study highlights the role of modeling in evaluating the experimental design and identifying the potential of system mischaracterization and its implications.

54 ENVIRONMENTAL SCIENCES↗

Modeling The Effects of Loading Scenario and Thermal Expansion Coefficient on Potential Failure of Cryo-compressed Hydrogen Vessels

A multiscale thermomechanical model for a simplified Type-3 cryogenic compressed hydrogen (H2) storage vessel is developed in this paper. The model accounts for the temperature-dependent elastic-plastic behavior of the vessel carbon/epoxy composite overwrap and aluminum alloy liner. The homogenized thermo-elastic-plastic behavior for the individual laminas of the vessel layup is obtained by an incremental Eshelby-Mori-Tanka approach associated with a micromechanical failure criterion to predict lamina failure while a standard elastic-plastic constitutive model is used to describe the behavior of a typical aluminum alloy assumed for the liner. The vessel response to external loadings is achieved by a finite element method. Four loading scenarios representing four thermomechanical cycles applied to the vessel are analyzed to evaluate constituent and lamina stresses as well as the associate failure criterion during the cycle according to these scenarios. The model can provide helpful guidance to mitigate thermal stresses by an adequate selection of loading scenario, optimizing the layup and by tailoring thermomechanical properties of the resin matrix.

36 MATERIALS SCIENCE↗