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

All sources NASA NTRS DOE OSTI

NASA NTRS2017

High-Resolution NO2 Observations from the Airborne Compact Atmospheric Mapper: Retrieval and Validation

Nitrogen dioxide (NO2) is a short-lived atmospheric pollutant that serves as an air quality indicator and is itself a health concern. The Airborne Compact Atmospheric Mapper (ACAM) was flown on board the NASA UC-12 aircraft during the Deriving Information on Surface Conditions from Column and Vertically Resolved Observations Relevant to Air Quality Maryland field campaign in July 2011. The instrument collected hyperspectral remote sensing measurements in the 304-910 nm range, allowing daytime observations of several tropospheric pollutants, including nitrogen dioxide (NO2), at an unprecedented spatial resolution of 1.5 × 1.1 sq. km. Retrievals of slant column abundance are based on the differential optical absorption spectroscopy method. For the air mass factor computations needed to convert these retrievals to vertical column abundance, we include high-resolution information for the surface reflectivity by using bidirectional reflectance distribution function data from the Moderate Resolution Imaging Spectroradiometer. We use high-resolution simulated vertical distributions of NO2 from the Community Multiscale Air Quality and Global Modeling Initiative models to account for the temporal variation in atmospheric NO2 to retrieve middle and lower tropospheric NO2 columns (NO2 below the aircraft). We compare NO2 derived from ACAM measurements with in situ observations from NASA's P-3B research aircraft, total column observations from the ground-based Pandora spectrometers, and tropospheric column observations from the space-based Ozone Monitoring Instrument. The high-resolution ACAM measurements not only give new insights into our understanding of atmospheric composition and chemistry through observation of subsampling variability in typical satellite and model resolutions, but they also provide opportunities for testing algorithm improvements for forthcoming geostationary air quality missions.

Lamsal, L. N.↗

DOE OSTI2015

Assessment of current atomic scale modelling methods for the investigation of nuclear fuels under irradiation: Example of uranium dioxide

Multiscale approaches are developed to build more physically based kinetic and mechanical mesoscale models to enhance the predictive capability of fuel performance codes and increase the efficiency of the development of the safer and more innovative nuclear materials needed in the future. Atomic scale methods, and in particular electronic structure and empirical potential methods, form the basis of this multiscale approach. It is therefore essential to know the accuracy of the results computed at this scale if we want to feed them into higher scale models. We focus here on the assessment of the description of interatomic interactions in uranium dioxide using on the one hand electronic structure methods, in particular in the density functional theory (DFT) framework and on the other hand empirical potential methods. These two types of methods are complementary, the former enabling to get results from a minimal amount of input data and further insight into the electronic and magnetic properties, while the latter are irreplaceable for studies where a large number of atoms needs to be considered. We consider basic properties as well as specific ones, which are important for the description of nuclear fuel under irradiation. These are especially energies, which are the main data passed to higher scale models. We limit ourselves to uranium dioxide.

multiscale modeling↗

NASA NTRSDate not supplied

Microscale Constitutive Model Sensitivity on Multiscale Modeling of Fiber Reinforced Composites

Fiber reinforced composites are desirable in applications where low weight and high strength are needed, but are susceptible to microscale variability during manufacturing, making failure predictions difficult. The impact microscale variability has on macroscale mechanical response is difficult to predict due to the computational efficiency needed to simulate many, large, high fidelity, microscale models. In this study, a multiscale approach was taken to model 3-point bend, 4-point bend, and tensile experiments of a unidirectional composite from only having microstructure scans of these samples and constituent properties from literature. These scans were sampled with different sized windows, and statistically equivalent microstructures were generated, then simulated for stiffness, strength, and fracture toughness using an efficient micromechanical. Mesoscale models were created where element sizes equaled microstructure size, and properties were assigned through sampling of microscale simulation results. First, this study showed the effect of using Weibull scaling on constituent matrix strength on macroscale response. Then, a comparison was made between different element sizes and experiments. Finally, model dimensions were fixed, and the effect of randomly distributed local properties alone was examined. Results showed that the scatter of strength and stiffness in the experiments could be predicted well using images of the microscale fiber morphologies and that using stochastic properties produced a 3% coefficient of variation of strength for all experiments.

statistical microstructure↗

DOE OSTI2025

Development and Validation of MALAMUTE model for Electric Field Assisted Sintering of Structural Materials

Fusion power plant designs feature extreme material performance requirements for structural material candidates. In addition to conventional alloys, more advanced composites and oxide dispersion strengthened (ODS) alloys are being explored, however, achieving the desired microstructures to maximize performance using traditional manufacturing methods can be challenging. The advanced manufacturing (AM) electric field-assisted sintering (EFAS) technique offers improved control over the final microstructure through higher heating and cooling rates and moderate pressures. Modeling and simulation tools show promise in elucidating the process-structure-property-performance (PSPP) correlation for AM-produced parts, including the EFAS process. An inherently multiscale process, the EFAS technique aligns well with the multiscale modeling capability of the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE)[cite]. We present here an electro-thermo-mechanical approach to modeling the EFAS process using the MOOSE Application Library for Advanced Manufacturing UTilitiEs (MALAMUTE) code. Prediction of the field and gradient distributions across the EFAS tooling is required to accurately describe the conditions for the lower-scale microstructural evolution models. In this work we present the MALAMUTE model developed to predict the electrical potential, temperature, and mechanical stress distribution across the EFAS graphite tooling and part at the larger engineering-scale. Validation of the MALAMUTE engineering-scale model is completed using data from experimental densification and pre-densified runs of iron powder via EFAS at 1000oC. These runs were conducted using a Thermal Technology DCS-5 EFAS system. Data collected during the experiment runs include the direct current (DC) supplied to the graphite tooling, the temperature of the graphite tooling as measured with a pyrometer, and the force applied to the top of the graphite tooling stack, and the data were recorded every 10 seconds. Our validation approach used the current and force data from the EFAS run as boundary condition inputs to the MAMALUTE simulation; the temperature data were used to evaluate the MALAMUTE EFAS model prediction. Results of the MALAMUTE simulations are employed to connect the external pyrometer temperature measurement to the temperature profile across the part undergoing consolidation. We investigate the impact of material property variation and mesh deformation on the temperature profile as predicted by MALAMUTE. We conclude by highlighting projects where the MALAMUTE EFAS modeling and simulation capabilities will be used to assist experimental design.

36 - MATERIALS SCIENCE↗

DOE OSTI2022

Variational multiscale reinforcement learning for discovering reduced order closure models of nonlinear spatiotemporal transport systems

Abstract A central challenge in the computational modeling and simulation of a multitude of science applications is to achieve robust and accurate closures for their coarse-grained representations due to underlying highly nonlinear multiscale interactions. These closure models are common in many nonlinear spatiotemporal systems to account for losses due to reduced order representations, including many transport phenomena in fluids. Previous data-driven closure modeling efforts have mostly focused on supervised learning approaches using high fidelity simulation data. On the other hand, reinforcement learning (RL) is a powerful yet relatively uncharted method in spatiotemporally extended systems. In this study, we put forth a modular dynamic closure modeling and discovery framework to stabilize the Galerkin projection based reduced order models that may arise in many nonlinear spatiotemporal dynamical systems with quadratic nonlinearity. However, a key element in creating a robust RL agent is to introduce a feasible reward function, which can be constituted of any difference metrics between the RL model and high fidelity simulation data. First, we introduce a multi-modal RL to discover mode-dependant closure policies that utilize the high fidelity data in rewarding our RL agent. We then formulate a variational multiscale RL (VMRL) approach to discover closure models without requiring access to the high fidelity data in designing the reward function. Specifically, our chief innovation is to leverage variational multiscale formalism to quantify the difference between modal interactions in Galerkin systems. Our results in simulating the viscous Burgers equation indicate that the proposed VMRL method leads to robust and accurate closure parameterizations, and it may potentially be used to discover scale-aware closure models for complex dynamical systems.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Variational quantum simulation of the critical Ising model with symmetry averaging

Here we investigate the use of deep multiscale entanglement renormalization ansatz (DMERA) circuits as a variational ansatz. We use the exactly solvable one-dimensional critical transverse-field Ising model as a test bed. Numerically exact simulation of the quantum circuit ansatz can in this case be carried out to hundreds of qubits by exploiting efficient classical algorithms for simulating matchgate circuits. We find that, for this system, the DMERA strongly outperforms a standard quantum approximate optimization algorithm (QAOA)–style ansatz, and that a major source of systematic error in correlation functions approximated using the DMERA is the breaking of the translational and Kramers-Wannier symmetries of the transverse-field Ising model. We are able to reduce this error by up to four orders of magnitude by symmetry averaging, without incurring additional cost in qubits or circuit depth. Here, we propose that this technique for mitigating systematic error could be applied to noisy intermediate-scale quantum (NISQ) simulations of physical systems with other symmetries.

1-dimensional spin chains↗

DOE OSTI2017

Turbulent mass transfer caused by vortex induced reconnection in collisionless magnetospheric plasmas

Magnetic reconnection is believed to be the main driver to transport solar wind into the Earth’s magnetosphere when the magnetopause features a large magnetic shear. However, even when the magnetic shear is too small for spontaneous reconnection, the Kelvin–Helmholtz instability driven by a super-Alfvénic velocity shear is expected to facilitate the transport. Although previous kinetic simulations have demonstrated that the non-linear vortex flows from the Kelvin–Helmholtz instability gives rise to vortex-induced reconnection and resulting plasma transport, the system sizes of these simulations were too small to allow the reconnection to evolve much beyond the electron scale as recently observed by the Magnetospheric Multiscale (MMS) spacecraft. Here in this paper, based on a large-scale kinetic simulation and its comparison with MMS observations, we show for the first time that ion-scale jets from vortex-induced reconnection rapidly decay through self-generated turbulence, leading to a mass transfer rate nearly one order higher than previous expectations for the Kelvin–Helmholtz instability.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Optimal renormalization of multiscale systems

While model order reduction is a promising approach in dealing with multi-scale time-dependent systems that are too large or too expensive to simulate for long times, the resulting reduced order models can suffer from instabilities. We have recently developed a time-dependent renormalization approach to stabilize such reduced models. In the current work, we extend this framework by introducing a parameter that controls the time-decay of the memory of such models and optimally selecting this parameter based on limited fully resolved simulations. First, we demonstrate our framework on the inviscid Burgers equation whose solution develops a finite-time singularity. Our renormalized reduced order models are stable and accurate for long times while using for their calibration only data from a full order simulation before the occurrence of the singularity. Furthermore, we apply this framework to the 3D Euler equations of incompressible fluid flow, where the problem of finite-time singularity formation is still open and where brute force simulation is only feasible for short times. Our approach allows us to obtain for the first time a perturbatively renormalizable model which is stable for long times and includes all the complex effects present in the 3D Euler dynamics. We find that, in each application, the renormalization coefficients display algebraic decay with increasing resolution, and that the parameter which controls the time-decay of the memory is problem-dependent.

Price, Jacob↗

DOE OSTI2023

Autonomous MultiScale Library

AMSLib provides infrastructure to tightly couple multi-scale physics simulation code with ML surrogate model inference. It provides a wholistic runtime execution paradigm to supports uncertainty quantification, surrogate model inference, persistent data storing throughout the execution of a simulation.

Bhatia, Harsh↗

DOE OSTI2018

Integrated Particle- and Reactor-Scale Simulation of Pine Pyrolysis in a Fluidized Bed

We report results from a multiscale computational modeling study of biomass fast pyrolysis in an experimental laboratory reactor that combined the hydrodynamics predicted by a two-fluid model (TFM) with predictions from a finite element model (FEM) of heat and mass transfer and chemical reactions within biomass particles. The experimental pyrolyzer consisted of a 2-inch (5.1 cm)-diameter bubbling fluidized bed reactor (FBR) fed with milled pine pellets. The predicted FBR hydrodynamics included estimates of the residence times that the gas and biomass particles spend in the reactor before they exit. A single-particle FEM model was constructed based on geometry and heat transfer properties determined from optical and X-ray Computed Tomography (XCT) measurements of wood and char particles collected from the experimental FBR, along with previously proposed pyrolysis reaction kinetics. Taken together, the combined TFM and FEM simulation results predicted net bio-oil yields at the reactor exit that agree well with experimental observations, without any arbitrary fitting parameters. The combined computational models also provided practical information about the most important reactor and feedstock parameters.

09 BIOMASS FUELS↗

DOE OSTI2025

Mechanistic nuclear fuel performance modeling of uranium nitride

Uranium mononitride (UN) is a nuclear fuel candidate for advanced reactor designs and an alternative being considered for light water reactors due to its higher thermal conductivity and uranium density than UO 2 . As with any nuclear fuel, swelling and fission gas release are important factors for safety, while also being some of the hardest phenomena to predict with a high degree of confidence. Getting a grasp on the gas swelling behavior and release is crucial to lower the barrier for UN utilization. An accelerated swelling rate at high temperatures observed experimentally, sometimes referred to as “breakaway swelling,” further complicates the prediction of fuel performance of UN. A mechanistic model has been developed using a multiscale approach to describe the intragranular and intergranular fission gas behavior. Lower-length-scale calculations have been employed to inform models of the gas and self-diffusion behavior, resolution rate, and bubble shape. Leveraging previous work on high burnup UO 2 , two populations of intragranular bubbles are considered; small bulk bubbles and larger bubbles located along dislocations. The dislocation bubbles were found to be crucial to the overall swelling behavior, and the breakaway swelling transition was associated with the transition in the gas atom diffusion mechanism from an irradiation-induced athermal diffusion regime at lower temperatures to an intrinsic thermal equilibrium regime at higher temperatures, accelerating the growth of the dislocation bubbles. Similarly, the threshold for fission gas release was associated with the grain boundary vacancy diffusivity surpassing the gas atom diffusivity at sufficiently high temperatures, allowing the over-pressurized grain boundary bubble to grow in size and interconnect. Using thermo-mechanical models with the fission gas model, two integral fuel pin assessment cases were simulated. Finally, this work demonstrates the ability of a multiscale approach to accelerate the understanding of advanced fuel forms when experimental data is limited.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2020

Simulator for Hydrologic Unstructured Domains (SHUD v1.0): numerical modeling of watershed hydrology with the finite volume method

Abstract. Hydrologic modeling is an essential strategy for understanding and predicting natural flows, particularly where observations are lacking in either space or time or where complex terrain leads to a disconnect in the characteristic time and space scales of overland and groundwater flow. However, significant difficulties remain for the development of efficient and extensible modeling systems that operate robustly across complex regions. This paper introduces the Simulator for Hydrologic Unstructured Domains (SHUD), an integrated, multiprocess, multiscale, flexible-time-step model, in which hydrologic processes are fully coupled using the finite volume method. SHUD integrates overland flow, snow accumulation/melt, evapotranspiration, subsurface flow, groundwater flow, and river routing, thus allowing physical processes in general watersheds to be realistically captured. SHUD incorporates one-dimensional unsaturated flow, two-dimensional groundwater flow, and a fully connected river channel network with hillslopes supporting overland flow and baseflow. The paper introduces the design of SHUD, from the conceptual and mathematical description of hydrologic processes in a watershed to the model's computational structures. To demonstrate and validate the model performance, we employ three hydrologic experiments: the V-catchment experiment, Vauclin's experiment, and a model study of the Cache Creek Watershed in northern California. Ongoing applications of the SHUD model include hydrologic analyses of hillslope to regional scales (1 m2 to 106 km2), water resource and stormwater management, and interdisciplinary research for questions in limnology, agriculture, geochemistry, geomorphology, water quality, ecology, climate and land-use change. The strength of SHUD is its flexibility as a scientific and resource evaluation tool where modeling and simulation are required.

58 GEOSCIENCES↗

DOE OSTI2019

M3MS-19IN0101201 Modeling Intra-granular Fission Gas Bubble Evolution in Uranium Dioxide during Transients

The SciDAC milestone M3MS-19IN0101201 on Modeling Intra-granular Fission Gas Bubble Evolution in Uranium Dioxide during Transients was completed with the development, testing and implementation in the Bison code of an engineering model for intra-granular fission gas bubble evolution that is applicable to transient conditions. The model extends the previous work on modeling bubble evolution during normal operating conditions to include additional complexity that arises during transient situations, in particular, the so-called intra-granular bubble coarsening. This phenomenon refers to the appearance of a second population of coarsened bubbles with diameters of tens to hundreds of nm, which co-exist with nanometric bubbles present also during normal operating conditions and are responsible for a large portion of transient fuel swelling. The developed model considers growth of bubbles along dislocations, associated with the absorption of vacancies available in the dislocation core region, and favored by pipe diffusion of fission gas atoms. The model also utilizes the formulation for the gas atom resolution rate developed in the present SciDAC project through molecular dynamics simulations by other institutions. The present work therefore demonstrates an operational multiscale approach, based on collaboration within the project. The model is validated to an extensive experimental dataset of ramp tested uranium dioxide fuel samples. Additionally, we present new applications of the base model for bubble evolution during normal operating conditions to Bison engineering simulations of Cr2O3-doped fuel rods irradiated in the Halden reactor. Simulation results in terms of integral fission gas release in the fuel rods are compared to the available experimental data, pointing out significant improvements in the engineering predictions compared to application of the empirical model that was available in Bison before the SciDAC developments.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2021

Multiscale Models for Fibril Formation: Rare Events Methods, Microkinetic Models, and Population Balances

Amyloid fibrils are thought to grow by a two-step dock-lock mechanism. However, previous simulations of fibril formation (i) overlook the bi-molecular nature of the docking step and obtain rates with first-order units, or (ii) superimpose the docked and locked states when computing the potential of mean force for association and thereby muddle the docking and locking steps. Here, we developed a simple microkinetic model with separate locking and docking steps and with the appropriate concentration dependences for each step. We constructed a simple model comprised of chiral dumbbells that retains qualitative aspects of fibril formation. We used rare events methods to predict separate docking and locking rate constants for the model. The rate constants were embedded in the microkinetic model, with the microkinetic model embedded in a population balance model for “bottom-up” multiscale fibril growth rate predictions. These were compared to “top-down” results using simulation data with the same model and multiscale framework to obtain maximum likelihood estimates of the separate lock and dock rate constants. We used the same procedures to extract separate docking and locking rate constants from experimental fibril growth data. Our multiscale strategy, embedding rate theories, and kinetic models in conservation laws should help to extract docking and locking rate constants from experimental data or long molecular simulations with correct units and without compromising the molecular description.

Shayesteh Zadeh, Armin (ORCID:0000000299896710)↗

NASA NTRS2012

Multiscale Multifunctional Progressive Fracture of Composite Structures

A new approach is described for evaluating fracture in composite structures. This approach is independent of classical fracture mechanics parameters like fracture toughness. It relies on computational simulation and is programmed in a stand-alone integrated computer code. It is multiscale, multifunctional because it includes composite mechanics for the composite behavior and finite element analysis for predicting the structural response. It contains seven modules; layered composite mechanics (micro, macro, laminate), finite element, updating scheme, local fracture, global fracture, stress based failure modes, and fracture progression. The computer code is called CODSTRAN (Composite Durability Structural ANalysis). It is used in the present paper to evaluate the global fracture of four composite shell problems and one composite built-up structure. Results show that the composite shells. Global fracture is enhanced when internal pressure is combined with shear loads. The old reference denotes that nothing has been added to this comprehensive report since then.

Chamis, C. C.↗

DOE OSTI2023

Quantifying the Agyrotropy of Proton and Electron Heating in Turbulent Plasmas

An important aspect of energy dissipation in weakly collisional plasmas is that of energy partitioning between different species (e.g., protons and electrons) and between different energy channels. Here we analyse pressure–strain interaction to quantify the fractions of isotropic compressive, gyrotropic, and nongyrotropic heating for each species. An analysis of kinetic turbulence simulations is compared and contrasted with corresponding observational results from Magnetospheric Multiscale Mission data in the magnetosheath. In assessing how protons and electrons respond to different ingredients of the pressure–strain interaction, we find that compressive heating is stronger than incompressive heating in the magnetosheath for both electrons and protons, while incompressive heating is stronger in kinetic plasma turbulence simulations. Concerning incompressive heating, the gyrotropic contribution for electrons is dominant over the nongyrotropic contribution, while for protons nongyrotropic heating is enhanced in both simulations and observations. Variations with plasma β are also discussed, and protons tend to gain more heating with increasing β.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2023

Quantifying the Agyrotropy of Proton and Electron Heating in Turbulent Plasmas

Yan Yang↗

DOE OSTI2015

A Multiscale Model for the Quasi-Static Thermo-Plastic Behavior of Highly Cross-Linked Glassy Polymers

In this paper, we present experimentally validated molecular dynamics predictions of the quasi- static yield and post-yield behavior for a highly cross-linked epoxy polymer under gen- eral stress states and for different temperatures. In addition, a hierarchical multiscale model is presented where the nano-scale simulations obtained from molecular dynamics were homogenized to a continuum thermoplastic constitutive model for the epoxy that can be used to describe the macroscopic behavior of the material. Three major conclusions were achieved: (1) the yield surfaces generated from the nano-scale model for different temperatures agree well with the paraboloid yield crite- rion, supporting previous macroscopic experimental observations; (2) rescaling of the entire yield surfaces to the quasi-static case is possible by considering Argon’s theoretical predictions for pure compression of the polymer at absolute zero temperature; (3) nano- scale simulations can be used for an experimentally-free calibration of macroscopic con- tinuum models, opening new avenues for the design of materials and structures through multi-scale simulations that provide structure-property-performance relationships.

36 MATERIALS SCIENCE↗