Search NASA⌕ Search

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.

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

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

BISON Development and Validation for Priority LWR-ATF concepts

Over the years, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) (2015-2018, 2020) and Consortium for Advanced Simulation of Light Water Reactors (CASL) (2019) programs have provided support for development of Accident Tolerant Fuel (ATF) material models in the BISON fuel performance code. Since the beginning, the goal has been to utilize a multiscale modeling approach to gain a physical understanding of the fuel concepts of interest and to develop mechanistic models in the absence of a large amount of experimental data. This work builds upon that of previous years. In particular we present newly updated fission gas release models for both gas behavior in Cr 2 O 3 -doped UO 2 and U 3 Si 2 fuels, and a new creep model for U 3 Si 2 . The validation exercises completed last year are revisited with the latest models and the results updated. A brief summary of recent modeling activities for FeCrAl cladding is also provided.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2024

Filament Wound Composite Analysis Using the NASA Multiscale Analysis Tool (NASMAT) and Finite Element Analysis

Fiber reinforced composite materials, owing to their tailorable thermomechanical and functional properties, allow one to produce a structure that is stronger, stiffer, and lighter than its metal counterpart while performing the same function, yielding a more efficient structure. This not only allows for the improvement of current technologies like aircraft structures, but also enables new technologies like gaseous hydrogen storage for mobility applications, which are otherwise impractical when manufactured using traditional metals due to weight and space restrictions or material embrittlement. However, the use of composites imposes greater design and manufacturing challenges on an engineer, since they are heterogenous, having a distinct structure across multiple length scale, behave generally anisotropically at the structural level and require complex manufacturing and processing methods. Capturing this complex behavior requires detailed numerical simulations, including the modeling of microstructural features like undulations, voids, and fiber alignment. In this paper, multiple repeating unit cells (RUCs), representing filament wound composites, are developed (via a script provided in the Appendix) and analyzed. The refinement of these RUCs is varied, and the analyses are performed using both the Abaqus finite element software and the NASA Multiscale Analysis Tool (NASMAT). A study is undertaken to compare the predicted effective elastic properties of the wound RUC to a laminate representation of the wound RUC, which neglects the undulations. Additionally, two different sets of periodic boundary conditions (PBCs) have been examined. One approximates the real boundary conditions using a standard approach and the other represents the PBCs exactly through the use of an offset. Lastly, a comparison of the local elastic stress fields is made among the models and approaches. Since wound structures are often approximated as laminated structures, it is important to understand the degree to which this assumption is valid, namely by first comparing the elastic constants and local elastic fields. This will provide, on the one hand, information concerning the bulk mechanical behavior and, on the other hand, insights concerning local load distributions and likely damage initiation sites.

Composites↗
DOE OSTI2025

Observation of a Knotted Electron Diffusion Region in Earth's Magnetotail Reconnection

Magnetic reconnection is a fundamental plasma process that alters the magnetic field topology and releases magnetic energy. Most numerical simulations and spacecraft observations assume a two‐dimensional diffusion region, with the electron diffusion region (EDR) embedded in the same plane as the ion diffusion region (IDR) and a uniform guide field throughout. Using observations from Magnetospheric Multiscale mission, we report a non‐coplanar, knotted EDR in Earth's magnetotail current sheet. The reconnection plane of the knotted EDR deviates by approximately 38° from that of the IDR, with the guide field exhibiting both a 38° directional shift and a twofold increase in amplitude. Moreover, the Hall magnetic field is bipolar in the EDR but quadrupolar in the IDR, indicating different Hall current structures at electron and ion scales. These observations highlight the importance of three‐dimensional effects and illustrate the complexity of multiscale coupling between the EDR and IDR during reconnection studies.

diffusion region↗
NASA NTRS2004

Complexity and Intermittent Turbulence in Space Plasmas

Sporadic and localized interactions of coherent structures arising from plasma resonances can be the origin of "complexity" of the coexistence of non- propagating spatiotemporal fluctuations and propagating modes in space plasmas. Numerical simulation results are presented to demonstrate the intermittent character of the non-propagating fluctuations. The technique of the dynamic renormalization-group is introduced and applied to the study of scale invariance of such type of multiscale fluctuations. We also demonstrate that the particle interactions with the intermittent turbulence can lead to the efficient energization of the plasma populations. An example related to the ion acceleration processes in the auroral zone is provided.

Chang, Tom↗
DOE OSTI2026

Full Waveform Inversion Tomography for Improved Seismic Waveform Simulations and Lithospheric Structure of Northwest China and Surrounding Regions

We report a new model of three-dimensional (3D) anisotropic seismic wavespeed structure of the lithospheric of northwestern China, western Mongolia and east }|ern Kazakhstan from full waveform inversion tomography using openly available broadband data. To investigate the efficacy of a conservative multiscale inversion methodology, we used a simple isotropic starting model.

58 GEOSCIENCES↗
DOE OSTI2017

Multiscale Modeling of Deformation, Fracture and Failure of Fusion Materials and Structures Final Report

The main objective of the research performed at UCLA under grant number DE-FG02-03ER54708 was to develop a physical understanding of the underlying mechanisms that control high-temperature deformation and fracture of fusion structural materials. The key aspect is the synergistic interaction between thermo-mechanical damage phenomena (e.g. creep, fatigue, and fracture), and irradiation damage phenomena (helium and hydrogen bubbles, and irradiation-induced defects). Another thrust area of the research was on the development of constitutive equations for the mechanical behavior of irradiated structures in fusion components (mainly the First Wall and Plasma-Facing components), and in utilizing them in descriptions in large-scale Finite Element Modeling of fusion structures. A parallel experimental effort on thermo-mechanical and plasma-induced damage is designed to validate and confront the developed theoretical and computational predictions. The models have been experimentally validated through experimental tests in simulated plasma transients on specially-designed W samples. This grant started at UCLA on November 14, 2003 and has had successive continuations till November 14, 2017. Research performed under the grant was to develop a range of hierarchical models for the post-elastic deformation, fracture and failure of fusion structural materials. The developed multiscale modeling approach has been based on rigorous mathematical, physical, and computational methods at the forefront of computational materials science. At the fundamental and smallest length scale (nm-μm), we developed advanced rate theory and Monte Carlo approaches to model microstructure evolution, non-equilibrium phase transformations, and dislocation-defect interactions. Microscopic and mesoscopic models of radiation hardening, ductile-to-brittle-transition, post-elastic deformation, plastic instabilities, and fracture processes are based on Dislocation Dynamics (DD) and Grain Growth Dynamics (GGD). At the continuum level, we focused on modeling fracture and failure mechanisms for Virtual Integrated Testing (VISTA), and for understanding the limits of both materials-by-design versus structural component design through advancing the new approach of multiscale-multiphysics. The overall research has been strongly coupled with the U.S. national experimental program to ascertain and verify the range of investigated phenomena. In the following, we give a summary of research progress supported by this grant since its inception till November 14, 2017. Detailed research results are available in our publications, listed in the publications section of this report.

36 MATERIALS SCIENCE↗
DOE OSTI2021

First Principles Simulations of Electrified Interfaces in Electrochemistry

This chapter discusses some of the recent advances made in the first principles modeling of electrochemical catalysts. It also discusses the key development, namely the ability to explicitly treat the effects of surface electrification due to electrochemical processes and applied voltages in a computationally efficient manner. The chapter introduces the thermodynamics and statistical mechanics of electrified metal-solution interfaces. It then discusses the structure of the electrode-electrolyte interface and the effects of applied voltages, followed by an overview of the first-principles model and a motivating example. The chapter also provides a brief summary of classical thermodynamics and describes the basics of thermodynamic detour. It considers the thermodynamics of macroscopic systems and introduces several useful fundamental relations. Here, the chapter also considers a system to exist in a certain macrostate, which is a particular thermodynamic state specified by a set of fixed properties such as constant particle number, constant volume, and constant temperature.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Thermoviscoelastic modelling of high strain thin-ply composites by means of multiscale plate and beam model

High strain thin-ply (HS-TPC) technology is being increasingly adopted for high-performanceaerospace applications. Albeit many of these structures such as deployable booms can bemodeled as one-dimensional beam problems, there is a lack of themoviscoelastic beam mod-els to efficiently and accurately simulate HS-TPC. This work will use mechanics of struc-ture genome (MSG) to construct linear thermoviscoelastic beam models that can homogenizethree-dimensional heterogeneous materials made of constituents with time- and temperature-dependent behavior. The formulation derives the transient strain energy based on integralformulation for thermorheologically simple materials subject to finite temperature changeswith the restriction that the strain is small. A lenticular boom is used as a numerical exampleto verify the MSG-based linear thermoviscoelastic beam model against MSG-based shell/platemodel, which has already been validated against experimental data provided by NASA, anddirect numerical simulations performed in a finite element commercial package.

Orzuri Rique Garaizar↗
DOE OSTI2023

Mesoscale shock structure in particulate composites

Multiscale experiments in heterogeneous materials and the knowledge of their physics under shock compression are limited. Here, this study examines the multiscale shock response of particulate composites comprised of soda-lime glass particles in a PMMA matrix using full-field high speed digital image correlation (DIC) for the first time. Normal plate impact experiments, and complementary numerical simulations, are conducted at stresses ranging from 1.1 - 3.1 GPa to elucidate the mesoscale mechanisms responsible for the distinct shock structure observed in particulate composites. The particle velocity from the macroscopic measurement at continuum scale shows a relatively smooth velocity profile, with shock thickness decreasing with an increase in shock stress, and the composite exhibits strain rate scaling as the second power of the shock stress. In contrast, the mesoscopic response was highly heterogeneous, which led to a rough shock front and the formation of a train of weak shocks traveling at different velocities. Additionally, the normal shock was seen to diffuse the momentum in the transverse direction, affecting the shock rise and the rounding-off observed at the continuum scale measurements. The numerical simulations indicate that the reflections at the interfaces, wave scattering, and interference of these reflected waves are the primary mechanisms for the observed rough shock fronts.

36 MATERIALS SCIENCE↗
DOE OSTI2018

First-Principles Integrated Adsorption Modeling for Selective Capture of Uranium from Seawater by Polyamidoxime Sorbent Materials

Nuclear power is a relatively carbon-free energy source that has the capacity to be utilized today in an effort to stem the tides of global warming. The growing demand for nuclear energy, however, could put significant strain on our uranium ore resources, and the mining activities utilized to extract that ore can leave behind long-term environmental damage. A potential solution to enhance the supply of uranium fuel is to recover uranium from seawater using amidoximated adsorbent fibers. This technology has been studied for decades but is currently plagued by the material’s relatively poor selectivity of uranium over its main competitor vanadium. In this work, we investigate the binding schemes between uranium, vanadium, and the amidoxime functional groups on the adsorbent surface. Using quantum chemical methods, binding strengths are approximated for a set of complexation reactions between uranium and vanadium with amidoxime functionalities. Those approximations are then coupled with a comprehensive aqueous adsorption model developed in this work to simulate the adsorption of uranium and vanadium under laboratory conditions. Experimental adsorption studies with uranium and vanadium over a wide pH range are performed, and the data collected are compared against simulation results to validate the model. It was found that coupling ab initio calculations with process level adsorption modeling provides accurate predictions of the adsorption capacity and selectivity of the sorbent materials. Furthermore, this work demonstrates that this multiscale modeling paradigm could be utilized to aid in the selection of superior ligands or ligand compositions for the selective capture of metal ions. Furthermore, this first-principles integrated modeling approach opens the door to the in silico design of next-generation adsorbents with potentially superior efficiency and selectivity for uranium over vanadium in seawater.

36 MATERIALS SCIENCE↗
DOE OSTI2020

NRC Multiphysics Analysis Capability Deployment FY2020: Part 3

This report details progress and activities of Idaho National Laboratory (INL) on the Nuclear Regulatory Commission (NRC) project “Development and Modeling Sup- port for Advanced Non-Light Water Reactors.” The tasks completed for this report are: Task2c: Explicit modeling of pebble transient temperature response. In this simulation, the 400 MWth Pebble-Bed Modular Reactor (PBMR) design, PBMR- 400, experiences a 20-second power ramp from 100% to 150% power. This is followed by a similar reduction in the power back to 100%. Several multiscale pebble coupling approaches are tested with one pebble per mesh element in the active core region. The results show good conservation behavior and the stability of the coupling.; Extended scope part 1: An assessment of the computational efficiency of the Discontinuous Finite Element Method (DFEM) heat transfer solver shows good scalability. The DFEM solver is a factor of 4 more expensive in solution time than the Finite Element Method (FEM) solver for heat transfer problems due to the increased number of degrees of freedom. Nonetheless, the DFEM approach provides the user with the flexibility to model gap heat transfer problems.; Extended scope part 2: The GapHeatTransferInterfaceMaterial was improved to give the user increased flexibility with the modeling of heat transfer through gaps with the DFEM solver. A number of gap parameters can now be coupled both through functions and variables.; Extended scope part 3: Demonstration of how the gap width between hexagonal fuel cells can be calculated during a heat-up transient and used in the GapHeatTransferInterface model. A full-domain DFEM model with gap expansion is coupled to a SubApp that models the thermal expansion of the base plate. The results show the expected physical behavior, although have not been fully bench-marked at this point in time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Predictive multiscale modeling for Unidirectional Carbon Fiber Reinforced Polymers

This paper presents a predictive multiscale modeling scheme for Unidirectional (UD) Carbon Fiber Reinforced Polymers (CFRP). A bottom-up modeling procedure is discussed for predicting the performance of UD structures. UD material responses are computed from high fidelity Representative Volume Elements (RVEs). A data-driven Reduced Order Modeling approach compresses RVEs into Reduced Order Models so material responses can be computed in a concurrent fashion along with the structural level simulation. The approach presented in this paper is validated against experimental data and is believed to provide design guidance for future fiber reinforced polymers development.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Role of compression metallization in UO 2 fission-product energy cascade track: Multiscale electron-phonon analyses

While the electronic stoppage of charged fission fragments is relatively well understood, the subsequent energy cascade is not. Recent efforts to investigate this cascade and predict the resulting damage have used a two-temperature model (TTM) of the electronic and phononic systems coupled with a classical molecular dynamics (MD) simulation of the crystal lattice. In order to accurately predict the track radius produced by a fission fragment in UO 2 , this model (TTM + MD) requires that UO 2 , an insulator, have metallic properties, e.g., a substantial electron thermal conductivity and heat capacity. However, it has been predicted that UO 2 becomes metallic under large pressures, and we perform ab initio (DFT-HSE) simulations to support this prediction. We show that the average U-U bond length decreases within and near the ion track during TTM + MD simulations, supporting the use of volume contraction to model the pressurized UO 2 cell. Additionally, we evaluate the electron, phonon, and electron-phonon coupling properties of UO 2 for variations in the pressure. In particular, we calculate the electronic heat capacity and thermal conductivity, and the electron-phonon energy coupling for use in subsequent TTM + MD simulations. The ab initio parameterized TTM + MD simulations provide a set of the track radii predictions which bracket and include the experimentally observed radii. The accuracy of the ab initio parameterized TTM + MD simulations depends on the pressure and degree of electron-phonon non-equilibrium assumed during the ab initio calculations. We suggest improvements to the current TTM + MD methodology in light of these results. Here we show that the pressure-induced transition of UO 2 from insulator to metal and subsequent energy transfer from the electronic to phononic systems can accurately explain radiation damage during swift, heavy ion stoppage in UO 2 . We make some additional observations regarding the accumulation and recombination of damage along the ion track and make comparison to the common SRIM model of ion stoppage and damage accumulation.

36 MATERIALS SCIENCE↗
DOE OSTI2015

On the electron dynamics during island coalescence in asymmetric magnetic reconnection

We present an analysis of the electron dynamics during rapid island merging in asymmetric magnetic reconnection. We consider a doubly periodic system with two asymmetric transitions. The upper layer is an asymmetric Harris sheet of finite width perturbed initially to promote a single reconnection site. The lower layer is a tangential discontinuity that promotes the formation of many X-points, separated by rapidly merging islands. Across both layers, the magnetic field and the density have a strong jump, but the pressure is held constant. Our analysis focuses on the consequences of electron energization during island coalescence. We focus first on the parallel and perpendicular components of the electron temperature to establish the presence of possible anisotropies and non-gyrotropies. Thanks to the direct comparison between the two different layers simulated, we can distinguish three main types of behavior characteristic of three different regions of interest. The first type represents the regions where traditional asymmetric reconnections take place without involving island merging. The second type of regions instead shows reconnection events between two merging islands. Finally, the third regions identify the regions between two diverging island and where typical signature of reconnection is not observed. Electrons in these latter regions additionally show a flat-top distribution resulting from the saturation of a two-stream instability generated by the two interacting electron beams from the two nearest reconnection points. Finally, the analysis of agyrotropy shows the presence of a distinct double structure laying all over the lower side facing the higher magnetic field region. This structure becomes quadrupolar in the proximity of the regions of the third type. The distinguishing features found for the three types of regions investigated provide clear indicators to the recently launched Magnetospheric Multiscale NASA mission for investigating magnetopause reconnection involving multiple islands.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Multiresolution molecular mechanics: Implementation and efficiency

Atomistic/continuum coupling methods combine accurate atomistic methods and efficient continuum methods to simulate the behavior of highly ordered crystalline systems. Coupled methods utilize the advantages of both approaches to simulate systems at a lower computational cost, while retaining the accuracy associated with atomistic methods. Many concurrent atomistic/continuum coupling methods have been proposed in the past; however, their true computational efficiency has not been demonstrated. The present work presents an efficient implementation of a concurrent coupling method called the Multiresolution Molecular Mechanics (MMM) for serial, parallel, and adaptive analysis. First, we present the features of the software implemented along with the associated technologies. The scalability of the software implementation is demonstrated, and the competing effects of multiscale modeling and parallelization are discussed. Then, the algorithms contributing to the efficiency of the software are presented. These include algorithms for eliminating latent ghost atoms from calculations and measurement-based dynamic balancing of parallel workload. The efficiency improvements made by these algorithms are demonstrated by benchmark tests. The efficiency of the software is found to be on par with LAMMPS, a state-of-the-art Molecular Dynamics (MD) simulation code, when performing full atomistic simulations. Speed-up of the MMM method is shown to be directly proportional to the reduction of the number of the atoms visited in force computation. Finally, an adaptive MMM analysis on a nanoindentation problem, containing over a million atoms, is performed, yielding an improvement of 6.3–8.5 times in efficiency, over the full atomistic MD method. For the first time, the efficiency of a concurrent atomistic/continuum coupling method is comprehensively investigated and demonstrated.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Assessing the Time Dependence of Reconnection With Poynting's Theorem: MMS Observations

We investigate the time dependence of electromagnetic–field–to–plasma energy conversion in the electron diffusion region of asymmetric magnetic reconnection. To do so, we consider the terms in Poynting's theorem. In a steady state there is a perfect balance between the divergence of the electromagnetic energy flux and the conversion between electromagnetic field and particle energy. This energy balance is demonstrated with a particle–in–cell simulation of reconnection.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Minimization of Cathode|Solid-Electrolyte Interfacial Delamination through the Application of Interphase Layers

Next-generation lithium-ion batteries are expected to use solid electrolytes (SEs) to enable higher energy density and extreme fast-charge capabilities. One major mode of degradation at the cathode|SE interface is delamination between the cathode active materials and SEs, which leads to performance decay. Experimental observations indicate that implementation of interphase layers can minimize the cathode|SE delamination induced capacity fade. A multiscale computational methodology is developed here to investigate the applicability of boron substituted lithium carbonate (Li 2+x B x C 1–x O 3 , x = 0.5, or LBCO) to minimize the delamination at the cathode|SE interface. Atomistic simulations indicate that the fracture energies at both the cathode|LBCO and LBCO|SE interfaces are higher than those at the cathode|SE interface, which reduces the extent of delamination. Mesoscale simulations indicate that, apart from increasing the fracture energy, decreasing the evolution of strain energy by lowering the elastic modulus of the interphase layer can also minimize the extent of delamination at the cathode|SE interface. However, the adoption of an interphase layer with high ionic conductivity is necessary to minimize the ohmic losses during operation at higher current densities. This study provides guidance on selecting interphase layers with specific properties and thicknesses to minimize both interfacial delamination and impedance growth.

LBCO↗
NASA NTRSDate not supplied

Modeling of Multiscale Solar Dynamics for Understanding Drivers of Space Weather

Understanding the solar dynamics is critical for improving our capabilities to forecast the evolution of space weather conditions. We take advantage of currently available computational capabilities to model solar dynamics from the deep interior to the corona and investigate mechanisms that may drive space weather conditions. The simulations are performed using the 3D radiative MHD code StellarBox. Comparison of synthetic spectroscopic observables obtained from numerical simulations and actual observations allows us to uncover physical processes associated with observed phenomena. To facilitate a transition from modeling short-term physical phenomena to developing a reliable forecast-oriented model, we suggest using the data assimilation approach. It allows us to cross-analyze dynamo model solutions and observations and to consider possible uncertainties and errors. In this presentation, we briefly summarize current multi-scale modeling capabilities and results and discuss ongoing developments to build a reliable physics-based forecast-oriented model of solar activity.

Irina N Kitiashvili↗