SEARCH · Search NASA

Results for “Multiscale simulations”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Multiscale and Multifidelity Modeling of a 3D Woven Composite Thermal Protection System

Complex three-dimensional (3D) woven composites have been considered by multiple NASA projects in recent years as a means of offering improved mechanical and thermal performance over traditional laminated composite systems. Parallel efforts have focused on developing simulation capabilities for these systems, which have traditionally and heavily relied on experimental testing to evaluate composite performance. One system is the Heatshield for Extreme Entry Environment Technology (HEEET), which is being considered for the thermal protection system on reentry spacecraft. Optical microscopy was used to characterize the blended carbon and phenolic fiber tows. A section of HEEET insulation layer was imaged with high-resolution micro-computed tomography (microCT) and segmented to separate individual tows, porous matrix, and voids. These data were used to develop multiscale thermomechanical computational models within the NASA Multiscale Analysis Tool (NASMAT). Two NASMAT modeling approaches were considered to capture the details of the 3D woven architecture: a coarse model appropriate for inclusion in multiscale structural analyses and a high-fidelity model created by downsampling the microCT data. Both elastic and thermal properties were computed and compared. The feasibility and challenges associated with modeling complex, hybrid 3D woven composites were also addressed.

NASMAT↗

DOE OSTI2026

Difference in neoclassical edge flows between strongly negative and positive triangularities in the XGC gyrokinetic simulation

The neoclassical baseline study of a strongly negative triangularity (NT) plasma and the corresponding positive triangularity (PT) plasma is performed using the edge-specialized, total-f gyrokinetic code (X-point Gyrokinetic Code, XGC). A DIII-D-like plasma is used, based on the negative triangularity discharge of DIII-D #193 793. An artificial positive triangularity (PT) equilibrium has been constructed to compare the edge rotation physics at the same triangularity strength but with opposite sign, while keeping the same elongation and other geometric parameters. Carbon(6+) ions are added to the deuterium plasma at an experimentally relevant level. By using the experimental profile of carbon toroidal rotation profile as an input, XGC finds that the deuteron rotation is significantly different from the carbon rotation at the inboard and outboard midplanes, mostly caused by the difference in the Pfirsch–Schlüter rotation. More importantly, a significant difference in the X-point orbit loss physics, thus the rotation source, is found between the positive and negative triangularity equilibrium models. However, it is also found that the agreement between the present neoclassical simulation and the experimental NT data is validated only within the middle of the pedestal slope, indicating the importance of edge turbulence. Finally, this study could establish a baseline for the multiphysics, multiscale studies that include turbulence of negative triangularity plasmas.

Deuterium↗

DOE OSTI2019

Multiscale model based design of an energy–intensified novel adsorptive reactor process for the water gas shift reaction

In this work, an adsorptive reactor (AR) process is considered that can energetically intensify the water gas shift reaction (WGSR). To best understand AR process behavior, a multiscale, dynamic, process model is developed. This multiscale model enables the quantification of catalyst and adsorbent effectiveness factors within the reactor environment, obliviating the commonly employed assumption that these factors are constant. Simulations of the AR's alternating adsorption-reaction/desorption operation, using the proposed model, illustrate rapid convergence to a long-term periodic solution. The obtained simulation results quantify the influence of key operating conditions and design parameters (e.g., reactor temperature/pressure, W cat /F CO , W ad /F CO , F H2O /F CO ratios, and pellet size) on the AR's behavior. Finally, they also demonstrate, for pellet diameters used at the industrial scale, significant temporal and axial variation of the catalyst/adsorbent pellet effectiveness factors. Finally, the energetic intensification benefits of the proposed AR process over conventional WGSR packed-bed reactors are quantified.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Enabling AI-Accelerated Multiscale Modeling of Thrombogenesis at Millisecond and Molecular Resolutions on Supercomputers

We report the first congruent integration of HPC, AI, and multiscale modeling (MSM) for solving a mainstream biomechanical problem of thrombogenesis involving 6 million particles at record molecular-scale resolutions in space and at simulation rates of milliseconds per day. The two supercomputers, the IBM Summit-like AiMOS and our University’s SeaWulf, are used for scalability analysis of, and production runs with, the LAMMPS with our customization and AI augmentation and they attained optimal simulation speeds of 3,077 µs/day and 266 µs/day respectively. The long-time and large scales simulations enable the first study of the integrated platelet flowing, flipping, aggregating dynamics in one dynamically-coupled production run. The platelets’ angular and translational speeds, membrane particles’ speeds, and the membrane stress distributions are presented for the analysis of platelets’ aggregations.

Zhu, Yicong↗

NASA NTRS2001

Designing Adaptive Low-Dissipative High Order Schemes for Long-Time Integrations

A general framework for the design of adaptive low-dissipative high order schemes is presented. It encompasses a rather complete treatment of the numerical approach based on four integrated design criteria: (1) For stability considerations, condition the governing equations before the application of the appropriate numerical scheme whenever it is possible; (2) For consistency, compatible schemes that possess stability properties, including physical and numerical boundary condition treatments, similar to those of the discrete analogue of the continuum are preferred; (3) For the minimization of numerical dissipation contamination, efficient and adaptive numerical dissipation control to further improve nonlinear stability and accuracy should be used; and (4) For practical considerations, the numerical approach should be efficient and applicable to general geometries, and an efficient and reliable dynamic grid adaptation should be used if necessary. These design criteria are, in general, very useful to a wide spectrum of flow simulations. However, the demand on the overall numerical approach for nonlinear stability and accuracy is much more stringent for long-time integration of complex multiscale viscous shock/shear/turbulence/acoustics interactions and numerical combustion. Robust classical numerical methods for less complex flow physics are not suitable or practical for such applications. The present approach is designed expressly to address such flow problems, especially unsteady flows. The minimization of employing very fine grids to overcome the production of spurious numerical solutions and/or instability due to under-resolved grids is also sought. The incremental studies to illustrate the performance of the approach are summarized. Extensive testing and full implementation of the approach is forthcoming. The results shown so far are very encouraging.

Yee, Helen C.↗

NASA NTRS2020

Stochastic Simulation of Mudcrack Damage Formation in an Environmental Barrier Coating

The FEAMAC/CARES program, which integrates finite element analysis (FEA) with the MAC/GMC (Micromechanics Analysis Code with Generalized Method of Cells) and the CARES/Life (Ceramics Analysis and Reliability Evaluation of Structures / Life Prediction) programs, was used to simulate the formation of mudcracks during the cooling of a multilayered environmental barrier coating (EBC) deposited on a silicon carbide substrate. FEAMAC/CARES combines the MAC/GMC multiscale micromechanics analysis capability (primarily developed for composite materials) with the CARES/Life probabilistic multiaxial failure criteria (developed for brittle ceramic materials) and Abaqus (Dassault Systèmes) FEA. In this report, elastic modulus reduction of randomly damaged finite elements was used to represent discrete cracking events. The use of many small-sized low-aspect-ratio elements enabled the formation of crack boundaries, leading to development of mudcrack-patterned damage. Finite element models of a disk-shaped three-dimensional specimen and a twodimensional model of a through-the-thickness cross section subjected to progressive cooling from 1,300 °C to an ambient temperature of 23 °C were made. Mudcrack damage in the coating resulted from the buildup of residual tensile stresses between the individual material constituents because of thermal expansion mismatches between coating layers and the substrate. A two-parameter Weibull distribution characterized the coating layer stochastic strength response and allowed the effect of the Weibull modulus on the formation of damage and crack segmentation lengths to be studied. The spontaneous initiation of cracking and crack coalescence resulted in progressively smaller mudcrack cells as cooling progressed, consistent with a fractal-behaved fracture pattern. Other failure modes such as delamination, and possibly spallation, could also be reproduced. The physical basis assumed and the heuristic approach employed, which involves a simple stochastic cellular automaton methodology to approximate the crack growth process, are described. The results ultimately show that a selforganizing mudcrack formation can derive from a Weibull distribution that is used to describe the stochastic strength response of the bulk brittle ceramic material layers of an EBC.

Nemeth, Noel N.↗

DOE OSTI2020

Computer Simulation of Proton Transport in Fuel Cell Membranes (Final Report)

This DOE-supported research grant focused on understanding the nature of proton transport in complex systems such as proton exchange membranes (PEMs). The most unique aspect of the research was the development and implementation of a novel multiscale reactive molecular dynamics (MS-RMD) methodology. In this approach, covalent bonds can dynamically break and form, allowing one to accurately treat the proton hopping process essential to capturing the physics of proton transport. The Voth group applied this method to proton exchange membrane systems, providing insight into their proton transport mechanism. They found that protons can diffuse most rapidly in the water-rich regions, but that protons actually spend so little time in such regions that transport along the hydrophobic – hydrophilic interface controls the membrane performance. In addition, the group worked to increase understanding of acidic solutions, developing methods for simulating and interpreting experimental infrared vibrational spectroscopy for excess protons (acidic solutions). The group also implemented novel tools for developing proton transport reactive MD models using a relative entropy minimization scheme.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2026

Wavelet flow for extragalactic foreground simulations

Extragalactic foregrounds in cosmic microwave background (CMB) observations are both a source of cosmological and astrophysical information and a nuisance to the CMB. Effective field-level modeling that captures their non-Gaussian statistical distributions is increasingly important for optimal information extraction, particularly given the low-noise observations from current and upcoming experiments. Here, we explore the use of Wavelet Flow (WF) models to tackle the novel task of modeling the field-level probability distributions of multi-component CMB secondaries and foregrounds. Specifically, we jointly train correlated CMB lensing convergence (κ) and cosmic infrared background (CIB) maps with a WF model and obtain a network that statistically recovers the input to high accuracy — the trained network generates samples of κ and CIB fields whose average power spectra are within a few percent of the inputs across all scales, and whose Minkowski functionals are similarly accurate compared to the inputs. Leveraging the multiscale architecture of these models, we fine-tune both the model parameters and the priors at each scale independently, optimizing performance across different resolutions. These results demonstrate that WF models can accurately simulate correlated components of CMB secondaries, supporting improved analysis of cosmological data. Our code and trained models can be found on this GitHub repo.

cosmological simulations↗

DOE OSTI2016

Heliospheric current sheet and effects of its interaction with solar cosmic rays

The effects of interaction of solar cosmic rays (SCRs) with the heliospheric current sheet (HCS) in the solar wind are analyzed. A self-consistent kinetic model of the HCS is developed in which ions with quasiadiabatic dynamics can present. The HCS is considered an equilibrium embedded current structure in which two main plasma species with different temperatures (the low-energy background plasma of the solar wind and the higher energy SCR component) contribute to the current. The obtained results are verified by comparing with the results of numerical simulations based on solving equations of motion by the particle tracing method in the given HCS magnetic field with allowance for SCR particles. It is shown that the HCS is a relatively thin multiscale current configuration embedded in a thicker plasma layer. In this case, as a rule, the shear (tangential to the sheet current) component of the magnetic field is present in the HCS. Taking into account high-energy SCR particles in the HCS can lead to a change of its configuration and the formation of a multiscale embedded structure. Parametric family of solutions is considered in which the current balance in the HCS is provided at different SCR temperatures and different densities of the high-energy plasma. The SCR densities are determined at which an appreciable (detectable by satellites) HCS thickening can occur. Possible applications of this modeling to explain experimental observations are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRSDate not supplied

An Integrated Recursive Framework for Arbitrarily Multiscale and Multi-fidelity Modeling

The NASA Multiscale Analysis Tool (NAMSAT) serves as a state-of-the-art, “plug and play,” massively multiscale modeling (M3) platform for hierarchical materials and structures. The development of NASMAT has focused on modularity, upgradability and maintainability, interoperability, and utility. Application Program Interfaces (APIs) have been developed to facilitate the integration NASMAT into other programs (commercial, research, and user-defined) as well as the integration of other codes into NASMAT. NASMAT is an integrated, recursive platform that seamlessly allows for an arbitrary number of scales, and a variety of modeling fidelity, in a single, non-linear analysis. Thus, the macroscopic response of a material can be directly linked to the behavior of the microstructure(s), including subscale defects, making it a suitable tool for integrated computational materials engineering (ICME) because . Examples using NASMAT for simulation of polycrystalline metals at elevated temperatures will be presented.

Multiscale Modeling↗

DOE OSTI2024

Ab initio generalized Langevin equation

We introduce a machine learning–based approach called ab initio generalized Langevin equation (AIGLE) to model the dynamics of slow collective variables (CVs) in materials and molecules. In this scheme, the parameters are learned from atomistic simulations based on ab initio quantum mechanical models. Force field, memory kernel, and noise generator are constructed in the context of the Mori–Zwanzig formalism, under the constraint of the fluctuation–dissipation theorem. Combined with deep potential molecular dynamics and electronic density functional theory, this approach opens the way to multiscale modeling in a variety of situations. Here, we demonstrate this capability with a study of two mesoscale processes in crystalline lead titanate, namely the field-driven dynamics of a planar ferroelectric domain wall, and the dynamics of an extensive lattice of coarse-grained electric dipoles. In the first case, AIGLE extends the reach of ab initio simulations to a regime of noise-driven motions not accessible to molecular dynamics. In the second case, AIGLE deals with an extensive set of CVs by adopting a local approximation for the memory kernel and retaining only short-range noise correlations. The scheme is computationally more efficient than molecular dynamics by several orders of magnitude and mimics the microscopic dynamics at low frequencies where it reproduces accurately the dominant far-infrared absorption frequency.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Multiscale Modeling in the Clinic: Drug Design and Development

A wide range of length and time scales are relevant to pharmacology, especially in drug development, drug design and drug delivery. Therefore, multi-scale computational modeling and simulation methods and paradigms that advance the linkage of phenomena occurring at these multiple scales have become increasingly important. Multi-scale approaches present in silico opportunities to advance laboratory research to bedside clinical applications in pharmaceuticals research. This is achievable through the capability of modeling to reveal phenomena occurring across multiple spatial and temporal scales, which are not otherwise readily accessible to experimentation. The resultant models, when validated, are capable of making testable predictions to guide drug design and delivery. In this review we describe the goals, methods, and opportunities of multi-scale modeling in drug design and development. We demonstrate the impact of multiple scales of modeling in this field. We indicate the common mathematical techniques employed for multi-scale modeling approaches used in pharmacology and present several examples illustrating the current state-of-the-art regarding drug development for: Excitable Systems (Heart); Cancer (Metastasis and Differentiation); Cancer (Angiogenesis and Drug Targeting); Metabolic Disorders; and Inflammation and Sepsis. We conclude with a focus on barriers to successful clinical translation of drug development, drug design and drug delivery multi-scale models.

Clancy, Colleen E.↗

DOE OSTI2018

Computational Homogenization of Polycrystals

This paper reviews the current state of the art in the simulation of the mechanical behavior of polycrystalline materials by means of computational homogenization. The key ingredients of this modeling strategy are presented in detail starting with the parameters needed to describe polycrystalline microstructures and the digital representation of such microstructures in a suitable format to perform computational homogenization. The different crystal plasticity frameworks that can describe the physical mechanisms of deformation in single crystals (dislocation slip and twinning) at the microscopic level are presented next. This is followed by the description of computational homogenization methods based on mean-field approximations by means of the viscoplastic self-consistent approach, or on the full-field simulation of the mechanical response of a representative polycrystalline volume element by means of the finite element method or the fast Fourier transform-based method. Multiscale frameworks based on the combination of mean-field homogenization and the finite element method are presented next to model the plastic deformation of polycrystalline specimens of arbitrary geometry under complex mechanical loading. Examples of application to predict the strength, fatigue life, damage, and texture evolution under different conditions are presented to illustrate the capabilities of the different models. Lastly, current challenges and future research directions in this field are summarized.

36 MATERIALS SCIENCE↗

DOE OSTI2015

A systematic multiscale modeling and experimental approach to protect grain boundaries in magnesium alloys from corrosion

A multiscale modeling Internal State Variable (ISV) constitutive model was developed that captures the fundamental structure-property relationships. The macroscale ISV model used lower length scale simulations (Butler-Volmer and Electronics Structures results) in order to inform the ISVs at the macroscale. The chemomechanical ISV model was calibrated and validated from experiments with magnesium (Mg) alloys that were investigated under corrosive environments coupled with experimental electrochemical studies. Because the ISV chemomechanical model is physically based, it can be used for other material systems to predict corrosion behavior. As such, others can use the chemomechanical model for analyzing corrosion effects on their designs.

36 MATERIALS SCIENCE↗

DOE OSTI2016

Energy Storage Annual Progress Report for FY15

The Energy Storage research and development (R&D) subprogram within the DOE Vehicle Technologies Office (VTO) provides support and guidance for projects focusing on batteries for plug-in electric vehicles (PEVs) in support of the EV Everywhere Grand Challenge. PEVs could have a significant impact on the nation's goal of reducing dependence on imported oil and gaseous pollutant emissions. The Energy Storage program targets overcoming technical barriers to enable market success, including: (1) significantly reducing battery cost; (2) increasing battery performance (power, energy, durability); (3) reducing battery weight and volume; and (4) increasing battery tolerance to abusive conditions such as short circuit, overcharge, and crush. The National Renewable Energy Laboratory (NREL) supports the VTO's Energy Storage program by evaluating the thermal performance of cells and packs, developing electrochemical-thermal models to accelerate the design cycle for developing batteries, investigating the behavior of lithium-ion batteries under abuse conditions such as crush, enhancing the durability of electrodes by coatings such as atomic layer deposition, synthesis of materials for higher energy density batteries, and conducting techno-economic analysis of batteries in various electric-drive vehicles. This report describes the progress made by NREL on the research and development projects funded by the DOE VTO Energy Storage subprogram in FY15.

25 ENERGY STORAGE↗

NASA NTRS2009

Prediction of Material Properties of Nanostructured Polymer Composites Using Atomistic Simulations

Atomistic models of epoxy polymers were built in order to assess the effect of structure at the nanometer scale on the resulting bulk properties such as elastic modulus and thermal conductivity. Atomistic models of both bulk polymer and carbon nanotube polymer composites were built. For the bulk models, the effect of moisture content and temperature on the resulting elastic constants was calculated. A relatively consistent decrease in modulus was seen with increasing temperature. The dependence of modulus on moisture content was less consistent. This behavior was seen for two different epoxy systems, one containing a difunctional epoxy molecule and the other a tetrafunctional epoxy molecule. Both epoxy structures were crosslinked with diamine curing agents. Multifunctional properties were calculated with the nanocomposite models. Molecular dynamics simulation was used to estimate the interfacial thermal (Kapitza) resistance between the carbon nanotube and the surrounding epoxy matrix. These estimated values were used in a multiscale model in order to predict the thermal conductivity of a nanocomposite as a function of the nanometer scaled molecular structure.

Hinkley, J.A.↗

NASA NTRS2019

GPS Based Autonomous Navigation Study for the Lunar Gateway

This paper describes and predicts the performance of a conceptual autonomous GPS-based navigation system for NASA's planned lunar Gateway. This system is based on the flight-proven Magnetospheric Multiscale (MMS) GPS navigation system, augmented with an earth-pointed high-gain antenna and with an option for an atomic clock. High-fidelity simulations, calibrated against MMS flight data and making use of GPS transmitter patterns from the GPS Antenna Characterization Experiment (ACE) project are developed for operation of the system in the Gateway Near-Rectilinear Halo Orbit (NRHO). The results indicate that GPS can provide an autonomous, realtime navigation capability with comparable, or superior, performance to traditional Deep Space Network approach with eight hours of tracking per day.

Winternitz, Luke B.↗

DOE OSTI2017

Multiscale Investigation on Biofilm Distribution and Its Impact on Macroscopic Biogeochemical Reaction Rates

Biofilms are critical locations for biogeochemical reactions in the subsurface environment. The occurrence and distribution of biofilms at microscale as well as their impacts on macroscopic biogeochemical reaction rates are still poorly understood. This paper investigated the formation and distributions of biofilms in heterogeneous sediments using multiscale models, and evaluated the effects of biofilm heterogeneity on local and macroscopic biogeochemical reaction rates. Sediment pore structures derived from X-ray computed tomography were used to simulate the microscale flow dynamics and biofilm distribution in the sediment column. The response of biofilm formation and distribution to the variations in hydraulic and chemical properties was first examined. One representative biofilm distribution was then utilized to evaluate its effects on macroscopic reaction rates using nitrate reduction as an example. The results revealed that microorganisms primarily grew on the surfaces of grains and aggregates near preferential flow paths where both electron donor and acceptor were readily accessible, leading to the heterogeneous distribution of biofilms in the sediments. Furthermore, the heterogeneous biofilm distribution decreased the macroscopic rate of biogeochemical reactions as compared with those in homogeneous cases. Operationally considering the heterogeneous biofilm distribution in macroscopic reactive transport models such as using dual porosity domain concept can significantly improve the prediction of biogeochemical reaction rates. Overall, this study provided important insights into the biofilm formation and distribution in soils and sediments as well as their impacts on the macroscopic manifestation of reaction rates.

42 ENGINEERING↗