Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multiscale simulation”

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 613 records · Page 34

Solar Activity Across the Scales: From Small-Scale Quiet-Sun Dynamics to Magnetic Activity Cycles

Observations as well as numerical and theoretical models show that solar dynamics is characterized by complicated interactions and energy exchanges among different temporal and spatial scales. It reveals magnetic self-organization processes from the smallest scale magnetized vortex tubes to the global activity variation known as the solar cycle. To understand these multiscale processes and their relationships, we use a two-fold approach: 1) realistic 3D radiative MHD simulations of local dynamics together with high resolution observations by IRIS, Hinode, and SDO; and 2) modeling of solar activity cycles by using simplified MHD dynamo models and mathematical data assimilation techniques. We present recent results of this approach, including the interpretation of observational results from NASA heliophysics missions and predictive capabilities. In particular, we discuss the links between small-scale dynamo processes in the convection zone and atmospheric dynamics, as well as an early prediction of Solar Cycle 25.

Solar↗

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.↗

Stable Machine‐Learning Parameterization of Subgrid Processes in a Comprehensive Atmospheric Model Learned From Embedded Convection‐Permitting Simulations

Modern climate projections often suffer from inadequate spatial and temporal resolution due to computational limitations, resulting in inaccurate representations of sub-grid processes. A promising technique to address this is the multiscale modeling framework (MMF), which embeds a kilometer-resolution cloud-resolving model (CRM) within each atmospheric column of a host climate model to replace traditional convection and cloud parameterizations. Machine learning offers a unique opportunity to make MMF more accessible by emulating the embedded CRM and reducing its substantial computational cost. Although many studies have demonstrated proof-of-concept success of achieving stable hybrid simulations, it remains a challenge to achieve near operational-level success with real geography and comprehensive variable emulation that includes, for example, explicit cloud condensate coupling. In this study, we present a stable hybrid model capable of integrating for at least 5 years with near operational-level complexity, including coarse-grid geography, seasonality, explicit cloud condensate and wind predictions, and land coupling. Our model demonstrates skillful online performance, achieving a 5-year zonal mean tropospheric temperature bias within 2 K, water vapor bias within 1 g/kg, and a precipitation root mean square error of 0.96 mm/day. Key factors contributing to our online performance include an expressive U-Net architecture and physical thermodynamic constraints for microphysics. With microphysical constraints mitigating unrealistic cloud formation, our work is the first to demonstrate realistic multi-year cloud condensate climatology under the MMF framework. Despite these advances, online diagnostics reveal persistent biases in certain regions, highlighting the need for innovative strategies to further optimize online performance.

Hu, Zeyuan [NVIDIA Corporation, Santa Clara, CA (U↗

Structure- and Physics-Preserving Reductions of Power Grid Models

The large size of multiscale, distribution and transmission, power grids hinders fast systemwide estimation and real-time control and optimization of operations. This paper studies graph reduction methods of power grids that are favorable for fast simulations and follow-up applications. While the classical Kron reduction has been successful in reduced order modeling of power grids with traditional, hierarchical design, the selection of reference nodes for the reduced model in a multiscale, distribution and transmission, network becomes ambiguous. In this work we extend the use of the iterative Kron reduction by utilizing the electric grid's graph topology for the selection of reference nodes, consistent with the design features of multiscale networks. Additionally, we propose further reductions by aggregation of coherent subnetworks of triangular meshes, based on the graph topology and network characteristics, in order to preserve currents and build another power-flow equivalent network. Our reductions are achieved through the use of iterative aggregation of subgraphs that include general tree structures, lines, and triangles. Important features of our reduction algorithms include that (i) the reductions are either equivalent to the Kron reduction or otherwise produce a power-flow equivalent network; (ii) due to the mentioned power-flow equivalence, the reduced network can model the dynamic of the swing equations for a lossless, inductive, steady state network; (iii) the algorithms efficiently utilize hash tables to store the sequential reduction steps. The third feature allows for easy reintroduction of detailed models into the reduced, conceptual network and makes the final reduced order model backward compatible with a sequence of intermediate, partially reduced networks with varying resolution---the ordered sequence of iterative reductions corresponds to a sequence of reduced order models. The performance of our graph reduction algorithms and features of the reduced grids are discussed on a real-word transmission and distribution grid. We produce visualizations of the reduced models through open source libraries and release our reduction algorithms with example code and toy data.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Validation of a High-Fidelity Supersonic Parachute Inflation Dynamics Model and Best Practice

The parachute inflation dynamics (PID) of a Disk-Gap-Band (DGB) parachute system travelling at supersonic velocities in the high-altitude Earth atmosphere is simulated using a high-fidelity, multiphysics, massively-parallel computational framework. The sensitivity of the simulation results with respect to critical modelling and simulation parameters, including the fluid mechanical governing equations (Euler vs Navier-Stokes with RANS turbulence closure vs Navier-Stokes with LES turbulence closer), the interaction between the flow and the suspension lines (accounted for vs not accounted for), the dynamics of the payload (restrained vs free), and the material model of the heterogeneous fabric canopy (Saint-Venant Kirchoff vs mechanics-informed multiscale neural network), is extensively investigated. The parameter investigation is enabled by advancements to the state-of-the-art in the modelling and simulation of the supersonic parachute inflation process, which are presented as part of the computational framework. The predictive quality of the computational framework in the context of supersonic PID is validated by comparison of the simulation results with the experimental results from NASA’s Advanced Supersonic Parachute Research (ASPIRE) project. Through the validation, best practices for the modelling and simulation of supersonic PID are developed and proposed.

Rabinovitch, Jason↗

Information Management Workflow and Tools Enabling Multiscale Modeling Within ICME Paradigm

With the increased emphasis on reducing the cost and time to market of new materials, the need for analytical tools that enable the virtual design and optimization of materials throughout their processing - internal structure - property - performance envelope, along with the capturing and storing of the associated material and model information across its lifecycle, has become critical. This need is also fueled by the demands for higher efficiency in material testing; consistency, quality and traceability of data; product design; engineering analysis; as well as control of access to proprietary or sensitive information. Fortunately, material information management systems and physics-based multiscale modeling methods have kept pace with the growing user demands. Herein, recent efforts to establish workflow for and demonstrate a unique set of web application tools for linking NASA GRC's Integrated Computational Materials Engineering (ICME) Granta MI database schema and NASA GRC's Integrated multiscale Micromechanics Analysis Code (ImMAC) software toolset are presented. The goal is to enable seamless coupling between both test data and simulation data, which is captured and tracked automatically within Granta MI®, with full model pedigree information. These tools, and this type of linkage, are foundational to realizing the full potential of ICME, in which materials processing, microstructure, properties, and performance are coupled to enable application-driven design and optimization of materials and structures.

Materials Engineering↗

A Multiscale Computational Model of the Response of Swine Epidermis After Acute Irradiation

Radiation exposure from Solar Particle Events can lead to very high skin dose for astronauts on exploration missions outside the protection of the Earth s magnetic field [1]. Assessing the detrimental effects to human skin under such adverse conditions could be predicted by conducting territorial experiments on animal models. In this study we apply a computational approach to simulate the experimental data of the radiation response of swine epidermis, which is closely similar to human epidermis [2]. Incorporating experimentally measured histological and cell kinetic parameters into a multiscale tissue modeling framework, we obtain results of population kinetics and proliferation index comparable to unirradiated and acutely irradiated swine experiments [3]. It is noted the basal cell doubling time is 10 to 16 days in the intact population, but drops to 13.6 hr in the regenerating populations surviving irradiation. This complex 30-fold variation is proposed to be attributed to the shortening of the G1 phase duration. We investigate this radiation induced effect by considering at the sub-cellular level the expression and signaling of TGF-beta, as it is recognized as a key regulatory factor of tissue formation and wound healing [4]. This integrated model will allow us to test the validity of various basic biological rules at the cellular level and sub-cellular mechanisms by qualitatively comparing simulation results with published research, and should lead to a fuller understanding of the pathophysiological effects of ionizing radiation on the skin.

FROM↗

Physics-Guided, Physics-Informed, and Physics-Encoded Neural Networks and Operators in Scientific Computing: Fluid and Solid Mechanics

Abstract Advancements in computing power have recently made it possible to utilize machine learning and deep learning to push scientific computing forward in a range of disciplines, such as fluid mechanics, solid mechanics, materials science, etc. The incorporation of neural networks is particularly crucial in this hybridization process. Due to their intrinsic architecture, conventional neural networks cannot be successfully trained and scoped when data are sparse, which is the case in many scientific and engineering domains. Nonetheless, neural networks provide a solid foundation to respect physics-driven or knowledge-based constraints during training. Generally speaking, there are three distinct neural network frameworks to enforce the underlying physics: (i) physics-guided neural networks (PgNNs), (ii) physics-informed neural networks (PiNNs), and (iii) physics-encoded neural networks (PeNNs). These methods provide distinct advantages for accelerating the numerical modeling of complex multiscale multiphysics phenomena. In addition, the recent developments in neural operators (NOs) add another dimension to these new simulation paradigms, especially when the real-time prediction of complex multiphysics systems is required. All these models also come with their own unique drawbacks and limitations that call for further fundamental research. This study aims to present a review of the four neural network frameworks (i.e., PgNNs, PiNNs, PeNNs, and NOs) used in scientific computing research. The state-of-the-art architectures and their applications are reviewed, limitations are discussed, and future research opportunities are presented in terms of improving algorithms, considering causalities, expanding applications, and coupling scientific and deep learning solvers.

Computer Science↗

Tardigrade-examples V0.1.0

Tardigrade-examples is a repository of computational workflows that exercise the Tardigrade software package. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling. These capabilities include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This repository investigates continuum upscaling of various direct numerical simulations (DNS) conducted in Abaqus finite element (FE), Ratel FE, and GEOS material point method (MPM) software. Verification of the upscaling workflow is first investigated by considering DNS of trivial stress states for homogeneous materials, results of which indicate that classical continuum behavior is recovered as expected. DNS of heterogeneous materials are then considered.

Allard, Thomas↗

Quantum-Classical Science for the Plasma-Material Interface in NSTXU (Final Technical Report)

Plasma facing surfaces (PFS) have a profound influence on plasma performance and choice of material for PFSs is a burning issue for tokamak reactors. Using lithium wall conditioning in NSTX has resulted in a beneficial influence on plasma behavior including control of hydrogen recycling and improvements in energy confinement, longer confinement times, edge density control, and the elimination of Edge Localized Modes. A major challenge for the NSTX-U (increased power flux, increased campaign time) as well as for other magnetic fusion devices is to extend high-performance plasmas for very long durations, and to integrate this high performance with plasma facing components (PFCs) that can withstand very high heat and particle fluxes while maintaining structural integrity with minimal retention of fusion fuel. This project directly contributed to filling up the knowledge gaps and data for the dynamics of the plasma-facing interfaces in the NSTX-U at nano and nano-to-micro scales. This was achieved by both deciphering the poorly understood nanoscale phenomenology of plasma-surface interactions and providing the connection and comparison with the experiments which provide understanding at the mesoscopic scales. We modelled the dynamic response and evolution of a mixed material surfaces (D, Li, C, B, O) to bombardment by plasma particles, and investigated microstructure changes, erosion, surface chemistry, deuterium implantation and permeation. The main objectives of the project are i) a comparison of Li and B effects when deposited on carbon, and ii) the role of oxygen and other impurities (e.g. carbon, boron) in the lithium performance. In addition to predicting and understanding the phenomenology of the processes, the project achieved deliverables on the i) plasma induced erosion rates of PFCs (outer wall, divertor, etc.), including chemical and physical sputtering yields at various temperatures (300-700K), ii) deuterium uptake/recycling rates, and iii) accurate probability rates for the processes and reactions in i-ii which were used to understand the experimental results of the J.P. Allain and B. Koel NSTX-U projects. Lithium was a distinctive focus of the project. It easily charges in contact with other fusion material atoms, giving partially or fully its external electron cloud to other atoms in its vicinity and polarizing the mixed material. Our research strategy used experimentally validated multiscale theory applied to the large, multicomponent material surface systems. With the use of US supercomputing facilities we simulated evolution of atomic charges, chemistry, retention of deuterium, sputtering and morphology by introducing the corrections to the classical molecular dynamics, like are electronegativity equalization method (EEM) and quantum-classical molecular dynamics (QCMD). The goal was to understand and quantify the dynamics of the creation and evolution of the plasma-material interface under irradiation by atoms and molecules at boronized carbon, lithiated carbon, and boronized-lithiated carbon for both solid and self-healing liquid-metal divertor.

36 MATERIALS SCIENCE↗

Results of the Micromorphic Upscaling Workflow for the PSAAP III Year 3 Report

Predicting the mechanical response and failure of heterogeneous materials has proven difficult. Multiscale numerical methods based in higher order continuum theories attempt to bridge the gap between microscale and macroscale structural behavior. Micromorphic continuum theories have shown promise. The Tardigrade software package is an implementation of Eringen’s micromorphic continuum theory with capabilities to support multiscale material modeling workflows. These include homogenization through the Micromorphic Filter, calibration of micromorphic material models, and macroscale simulation in Tardigrade-MOOSE. This work discusses micromorphic upscaling efforts of the University of Colorado Boulder PSAAP III multidisciplinary simulation center (MSC). Verification studies are presented that compare the accuracy of the upscaling workflow with analytical solutions for a trivial stress state and homogeneous material using direct numerical simulations (DNS) conducted in the Ratel finite element method (FEM) and GEOS material point method (MPM) codes. These verification studies consider upscaling using the Micromorphic Filter for a “single filter domain”. Finally, DNS of a heterogeneous composite material is upscaled using a “multiple filter domain” method.

36 MATERIALS SCIENCE↗

GCAM-USA electricity demand results for National Climate Assessment 5

Overview This dataset includes GCAM-USA v 5.3 outputs for the percent change in electricity demand in the U.S. from 2020 to 2050 and from 2020 to 2100 for the thermodynamic global warming scenario "RCP8.5_hotter" and the SSP5 socioeconomic scenario. These results were produced as part of the Integrated Multisector, Multiscale Modeling (IM3) project. Detailed Information The electricity demand is calculated based on the IM3 GCAM-USA simulations. For the purpose of reproducibility, we provide the following data: 1. Raw data: the annual electricity demand for CONUS simulated by IM3 GCAM-USA for the scenario RCP8.5 Hotter - SSP5. 2. R scripts: process raw data, calculate percent change of electricity demand from 2020 to 2050 and from 2020 to 2100, and plot the data over CONUS. 3. Results: figures provided for NCA-5 and the corresponding data table from the R scripts.

Climate Change↗

Bridging scales in multiphysics VCSEL modeling

On our way to develop a comprehensive simulator of vertical-cavity surface-emitting lasers including carrier transport, optical modes, light-matter interaction and heat conduction, we address some critical multiscale aspects of the adopted coupling strategy. The quantum corrections to the semiclassical carrier transport framework are the bridge from nanometer to micrometer scales. In this paper, such corrections are shown to be fundamental in view of predicting the mode competition ruling VCSEL operation. Nevertheless, they can generate unforeseen features such as swirling electron flows in the active region, which are discussed here in detail. The simulation approach has been finally tested through a successful comparison with a large set of experimental results.

42 ENGINEERING↗

Multiscale Development and Validation of the Stainless Steel Alloy Corrosion (SStAC) Tool for High Temperature Engine Materials

The goal of this project was to create a simulation tool to assist in part design, reducing costs by eliminating the need for conservative material selection, and enabling alloy optimization to improve corrosion resistance. Thus, it helps decrease costs for existing engines and assists in the material and engine design for the engines of the future. To meet our stated goal, we developed the Stainless-Steel Alloy Corrosion (SStAC) tool, which: (1) Models corrosion of engine valves fabricated from 21-2N valve steel in an engine environment at temperatures up to 800 °C. (2) Runs in 1D for a fast estimate of the corrosion rate, and in 2D or 3D for more detailed simulations that represent any part geometry and predict the precise location and rate of corrosion and its impact on the mechanical and thermal behavior. (3) Was implemented using the open-source Multiphysics Object Oriented Simulation Environment (MOOSE). (4) Couples the corrosion model with mechanics and thermal transport models. The SStAC tool can predict the corrosion rate of 21-2N engine valves with less than 10% error.

36 MATERIALS SCIENCE↗

Effects of surface vacancy-adatom pair formation on PFC tungsten surface morphological response

Abstract We report a simulation study of the effect of He-irradiation-induced surface vacancy-adatom pair formation on the surface morphological evolution of plasma-facing component (PFC) tungsten and examine a number of factors that impact such evolution. Our analysis is based on self-consistent dynamical simulations according to an atomistically-informed, continuum-scale surface evolution model that has been developed following a hierarchical multiscale modeling strategy and can access the spatiotemporal scales of relevance to fuzz formation. The model accounts for the flux of surface adatoms generated as a result of the surface vacancy-adatom pair formation effect upon He implantation, which contributes to the anisotropic growth of surface nanostructural features due to the different rates of adatom diffusion along and across step edges of islands on the tungsten surface. We have carried out atomic-scale computations of optimal diffusion pathways along and across island step edges on the W(110) surface and calculated Ehrlich–Schwoebel (ES) barriers in adatom diffusion along and across such step edges. This aspect of surface adatom diffusion contributes to anisotropic surface atomic fluxes, terrace and step diffusive currents, and has been incorporated into our PFC surface evolution model, which predicts the formation of preferentially aligned nanoridge stripe patterns on the PFC surface. We establish that these anisotropic diffusive currents accelerate nanotendril growth on the PFC surface and the onset of surface nanostructure pattern formation. We also explore systematically the dependence of the PFC surface morphological response on the surface temperature and He ion incident flux, characterize in detail the resulting surface topographies and growth kinetics, and compare the predicted surface morphologies with experimental observations. Our simulation predictions for the emerging surface nanostructure patterns under certain plasma exposure conditions are consistent with experimental findings in the literature.

36 MATERIALS SCIENCE↗

Micromechanics Modeling of Textiles for Re-Entry Parachute Applications

Recent flight test projects and NASA missions have highlighted the challenges associated with accurately and efficiently modeling the behavior of parachute deployment systems needed for parachute design. Moreover, parachute deployment has been identified as one of the higher risk components for such missions. The analysis of textile fabrics used for atmospheric entry is inherently complex due to the multiple scales present in the fabric structure, including individual fiber filaments at the microscale, yarn bundles of fibers at the mesoscale, and the overall woven fabric at the macroscale. Computational tools for simulating fabric behavior must be able to account for the different mechanisms present at each scale without sacrificing computational efficiency. This work examines the generalized multiscale method of cells micromechanics theory, which has previously been used for the analysis of reinforced composite structures, to unreinforced textile fabrics. Modifications to the existing composite multiscale framework, implemented in NASA’s Multiscale Analysis Tool (NASMAT), include the specific mechanics unique to unreinforced textile fabrics, and overcoming the assumptions of a fixed fiber angle. It looks to assess the feasibility of using the NASMAT tool for efficient prediction of the response of unreinforced fabrics to loading such that it can ultimately be applied to fluid structure interaction tools for the prediction of parachute deployment systems. In this work, fabric behavior is simulated in NASMAT through homogenization of a triply periodic repeating unit cell, where the geometry of the subcells can change as a function of loading to represent the relative rotation and uncrimping that can occur in fabric tows. Predictions from the amended NASMAT code are compared to experimental data for uniaxial and off-axis tension to verify the ability of the code to incorporate lower-scale mechanics in prediction of unreinforced fabrics under loading.

Micromechanics↗

Multiscale Modeling of Carbon/Phenolic Composite Thermal Protection Materials: Atomistic to Effective Properties

Next generation ablative thermal protection systems are expected to consist of 3D woven composite architectures. It is well known that composites can be tailored to achieve desired mechanical and thermal properties in various directions and thus can be made fit-for-purpose if the proper combination of constituent materials and microstructures can be realized. In the present work, the first, multiscale, atomistically-informed, computational analysis of mechanical and thermal properties of a present day - Carbon/Phenolic composite Thermal Protection System (TPS) material is conducted. Model results are compared to measured in-plane and out-of-plane mechanical and thermal properties to validate the computational approach. Results indicate that given sufficient microstructural fidelity, along with lowerscale, constituent properties derived from molecular dynamics simulations, accurate composite level (effective) thermo-elastic properties can be obtained. This suggests that next generation TPS properties can be accurately estimated via atomistically informed multiscale analysis.

Materials Science↗

The Strong Impact of Weak Horizontal Convergence on Continental Shallow Convection

Idealized large-eddy simulation (LES) is a basic tool for studying three-dimensional turbulence in the planetary boundary layer. LES is capable of providing benchmark solutions for parameterization development efforts. However, real small-scale atmospheric flows develop in heterogeneous and transient environments with locally varying vertical motions inherent to open multiscale interactive dynamical systems. These variations are often too subtle to detect them by state-of-the-art remote and in situ measurements, and are typically excluded from idealized simulations. The present study addresses the impact of weak [i.e., O(10-6) s-1] short-lived low-level large-scale convergence/divergence perturbations on continental shallow convection. The results show a strong response of shallow nonprecipitating convection to the applied weak large-scale dynamical forcing. Evolutions of CAPE, mean liquid water path, and cloud-top heights are significantly affected by the imposed convergence/divergence. In contrast, evolving cloud-base properties, such as the area coverage and mass flux, are only weakly affected. To contrast those impacts with microphysical sensitivity, the baseline simulations are perturbed assuming different observationally based cloud droplet number concentrations and thus different rainfall. For the tested range of microphysical perturbations, the imposed convergence/divergence provides significantly larger impact than changes in the cloud microphysics. Simulation results presented here provide a stringent test for convection parameterizations, especially important for large-scale models progressing toward resolving some nonhydrostatic effects.

Meteorology & Atmospheric Sciences↗