SEARCH · Search NASA

Results for “Microscale Models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 37 records · Page 2

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Geometric Requirements for Modeling Polycrystalline Representative Volume Elements Using the Generalized Method of Cells

To accurately predict the mechanical response of advanced metallic parts, the grain structure of the metal must be considered and coupled to the macroscale in a multiscale model. The size of the microscale representative volume element must be sufficiently large, and the grid used to discretize the microscale must be adequately refined, in order to ensure an accurate and converged solution. The generalized method of cells is used to predict the mechanical properties of polycrystalline representative volume elements of copper containing randomly oriented grains with the properties of a single face-centered-cubic crystal. Moreover, the results of a parametric study are used to conclude the geometric requirements for an objective microscale model, with respect to the elastic stiffness.

Multiscale Modeling↗

DOE OSTI2021

A Novel Multiphysics Multiscale Multiporosity Shale Gas Transport Model for Geomechanics/Flow Coupling in Steady and Transient States

Summary A novel multiphysics multiscale multiporosity shale gas transport (M3ST) model was developed to investigate shale gas transport in both transient and steady states. The microscale model component contains a kerogen domain and an inorganic matrix domain, and each domain has its own geomechanical and gas transport properties. Permeabilities of various shale cores were measured in the laboratory using a pulse decay permeameter (PDP) with different pore pressure and confining stress combinations. The PDP-measured apparent permeability as a function of pore pressure under two effective stresses was fitted using the microscale M3ST model component based on nonlinear least squares fitting (NLSF), and the fitted model parameters were able to provide accurate model predictions for another effective stress. The parameters and petrophysical properties determined in the steady state were then used in the transient-state, continuum-scale M3ST model component, which performed history matching of the evolutions of the upstream and downstream gas pressures. In addition, a double-exponential empirical model was developed as a powerful alternative to the M3ST model to fit laboratory-measured apparent permeability under various effective stresses and pore pressures. The developed M3ST model and the research findings in this study provided critical insights into the role of the multiphysics mechanisms, including geomechanics, fluid dynamics and transport, and the Klinkenberg effect on shale gas transport across different spatial scales in both steady and transient states.

Engineering↗

DOE OSTI2021

Modeling Electrokinetics of Oxygen Electrodes in Solid Oxide Electrolyzer Cells

A microscale model is presented in this study to simulate electrode kinetics of the oxygen electrode in a solid oxide electrolyzer cell (SOEC). Two mixed ionic/electronic conducting structures are examined for the oxygen producing electrode in this work: single layer porous lanthanum strontium cobalt ferrite (LSCF), and bilayer LSCF/SCT (strontium cobalt tantalum oxide) structures. A yttrium-stabilized zirconia (YSZ) electrolyte separates the hydrogen and oxygen electrodes, as well as a gadolinium doped-ceria (GDC) buffer layer on the oxygen electrode side. Electrochemical reactions occurring at the two-phase boundaries (2PBs) and three-phase boundaries (3PBs) of single-layer LSCF and bilayer LSCF/SCT oxygen electrodes are modeled under various SOEC voltages with lattice oxygen stoichiometry as the key output. The results reveal that there exists a competition in electrode kinetics between 2PBs and 3PBs, but 3PBs are the primary reactive sites for single-layer LSCF oxygen electrode under high voltages. These locations experience the greatest oxygen stoichiometry variations and are therefore the most likely locations for dimensional changes. By applying an active SCT layer over LSCF, the 2PBs become activated to compete with the 3PBs, thus alleviating oxygen stoichiometry variations and reducing the likelihood of dimensional change. This strategy could reduce lattice structural expansion, proving to be valuable for electrode-electrolyte delamination prevention and will be the focus of future work.

30 DIRECT ENERGY CONVERSION↗

NASA NTRSDate not supplied

In-situ Imaging of Pyrolyzing Aerospace Materials

Tracking morphological changes of materials during heating is crucial to understand its response in fire protection, biofuel production, thermal protection systems (TPS) for hypersonic flight. As materials are heated, they undergo physical and chemical changes due to water loss, stretching or shrinking, pyrolysis and chemical reactions in the ambient environment. The effects of these changes can have a profound impact on the material’s performance, indicated by changes in on the porosity and volume. While materials such as wood shrink as they pyrolyze and lose mass, others swell due to their inherent characteristics when exposed to heat [1]. This study focuses on experiments conducted at the Advanced Light Source (ALS) beamline 8.3.2, where in situ micro-computed tomography (µ-CT) is performed on materials as they are being pyrolyzed. Through in situ µ-CT, the change in total volume and porosity can be obtained in real-time, allowing for better understanding of the underlying thermophysical and chemical processes as a function of temperature. This study also focuses on the implementation of the Porous Microstructure Analysis software (PuMA) [2] to obtain thermal conductivity, permeability, and other properties of the material from the 3D tomographies. The information gained from these tomographies will supplement microscale model development of material morphological change and will aid macroscale modeling for high-temperature applications. For this study, Room Temperature Vulcanizing silicone (RTV) [3-5] is heated from room temperature to 1000°C using an infrared lamp heating system, and tomographies are continuously collected as the sample is heated. The tomographies are then segmented to obtain solid and void phases, from which estimates of pore size, porosity and total volume are extracted as a function of temperature. PuMA is deployed on the segmented tomographies to obtain thermal conductivity, permeability, and other properties as a function of temperature. Preliminary results show that RTV first intumesces (swells) as pyrolysis begins, due to build-up of pyrolysis gases in closed pores, and then shrinks significantly as more open pores are formed and the pyrolysis gases outgas. Pore network visualization of the tomographies using OpenPNM [6] showed the increase in pore connectivity with increase in temperature. Future work will focus on using PuMA to obtain macroscopic properties of RTV as a function of temperature.

Tomography↗

DOE OSTI2024

Modeling and observations of North Atlantic cyclones: Implications for U.S. Offshore wind energy

To meet the Biden-Harris administration's goal of deploying 30 GW of offshore wind power by 2030 and 110 GW by 2050, expansion of wind energy into U.S. territorial waters prone to tropical cyclones (TCs) and extratropical cyclones (ETCs) is essential. This requires a deeper understanding of cyclone-related risks and the development of robust, resilient offshore wind energy systems. Here, this paper provides a comprehensive review of state-of-the-science measurement and modeling capabilities for studying TCs and ETCs, and their impacts across various spatial and temporal scales. We explore measurement capabilities for environments influenced by TCs and ETCs, including near-surface and vertical profiles of critical variables that characterize these cyclones. The capabilities and limitations of Earth system and mesoscale models are assessed for their effectiveness in capturing atmosphere–ocean–wave interactions that influence TC/ETC-induced risks under a changing climate. Additionally, we discuss microscale modeling capabilities designed to bridge scale gaps from the weather scale (a few kilometers) to the turbine scale (dozens to a few meters). We also review machine learning (ML)-based, data-driven models for simulating TC/ETC events at both weather and wind turbine scales. Special attention is given to extreme metocean conditions like extreme wind gusts, rapid wind direction changes, and high waves, which pose threats to offshore wind energy infrastructure. Finally, the paper outlines the research challenges and future directions needed to enhance the resilience and design of next-generation offshore wind turbines against extreme weather conditions.

17 WIND ENERGY↗

NASA NTRSDate not supplied

Predictive Modeling of Carbon Ablators

Efforts to build a Predictive Material Modeling (PMM) framework from the micro-scale to the macro-scale are presented in this abstract. To reduce the need for extensive testing, accelerate the design cycle process, and reduce uncertainty margins applied to final designs, NASA is developing simulation and modeling tools that enable characterization of material properties and response to high-enthalpy environments. The Porous Microstructure Analysis (PuMA) code has been developed for computing macroscale (volume averaged) properties of porous materials using microscale images from micro-computed tomography (micro-CT). Microscale modeling requires a realistic representation of a material microstructure; these are obtained either synthetically during the design of the material or through X-ray micro-CT. Volume averaged properties are then used to inform macroscale material response models, such as those implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO) software, also actively developed by NASA. The computational model in PATO is a generic heat and mass transfer model for porous reactive materials containing several solid phases and a single gas phase. The detailed chemical interactions occurring between the solid phases and the gas phase are modeled at the pore scale assuming local thermal equilibrium. These tools were developed to efficiently interface with other pre-existing codes such as SPARTA (direct simulation Monte Carlo), DPLR (hypersonic CFD), NEQAIR (radiative transport) and DAKOTA (uncertainty quantification and optimization). Detailed flight data (Mars Science Laboratory [MSL] Entry Descent and Landing Instrument [MEDLI]) is critical for validating these computational tools for NASA applications. Examples of modeling ablative material response using these codes will be presented including 3D simulations of the full-scale heatshield of the MSL capsule. The simulations demonstrate the ability of the modern material response code, PATO, to handle the material response of geometrically complex and large domains, through the use of massively parallel computations.

Thermal Protection Systems↗

NASA NTRSDate not supplied

Predictive Modeling of Carbon Ablators Using Micro and Macro-Scale Modeling

Thermal Protection Systems↗

NASA NTRS1983

A nested-grid mesoscale numerical weather prediction model modified for Space Shuttle operational requirements

A nested-grid mesoscale atmospheric simulation system (MASS) is tested over Florida for the case of intense seabreeze-induced convection. The goal of this modeling system is to provide real-time aviation weather support which is designed to fit local terminal operations such as those supporting NASA's STS. Results from a 58 km and a 14.5 km nested-grid simulation show that this version of the MASS is capable of simulating many of the basic characteristics of convective complexes during periods of relatively weak synoptic scale flow regimes. However, it is noted that extensive development work is required with nested-grid cumulus and planetary boundary layer parameterization schemes before many of the meso-beta scale features such as thunderstorm downdraft-produced bubble high pressure centers can be accurately simulated. After these schemes are properly tuned, MASS can be utilized to initialize microscale modeling systems.

Kaplan, M. L.↗

DOE OSTI2025

Numerical Investigation of Susceptor-Catalyst Design for Ethylene Generation in Radio Frequency Based Reactors

Ethylene is the most widely produced petrochemical component in the world. Whether reactors are heated directly or indirectly via steam, manufacturers use economies of scale to overcome inherent thermodynamic inefficiencies when burning fossil fuels. While allowing large-scale operations to use alternative sources of energy and raw materials, new methods of supplying energy to chemical reactor systems can reduce the energy waste produced by conventional processes. One viable method for effectively supplying energy to reactor systems is electromagnetic (EM) induction heating. Using the properties of radio frequency (RF) waves, heterogeneous catalysts can be precisely targeted for heating inside reactors. Site-selective heating can greatly lower the energy requirements of the process by supplying heat to reaction sites while reducing needless heat transfer elsewhere. In this study, a microscale model was used to help create guidelines for susceptors and catalysts to improve ethylene production. The oxidative dehydrogenation of ethane is investigated by utilizing several catalysts and potential catalyst/susceptor combinations, with heat provided by an EM susceptor. Having the susceptors and catalyst function separately results in higher gradients in both heat and mass transfer, which drives transport through the catalyst region, while still providing adequate heating for endothermic reactions. However, using a susceptive core with a catalyst covering produces the maximum ethylene concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

A Multifunctional Isostructural Bilayer Oxygen Evolution Electrode for Durable Intermediate-Temperature Electrochemical Water Splitting

The overarching goal of the proposed research is to address SOEC’s degradation problem by advancing a new isostructural highly electrocatalytically active bilayer oxygen evolution reaction (OER) electrode, consisting of a LSCF (La 1-x Sr x Co 1-y Fe y O 3-δ ) core and a SCT (SrCo 0.9 Ta 0.1 O 3-δ ) shell, to achieve high and sustainable rate of oxygen evolution matching operating current densities without encountering delamination. To realize this goal, the project has adopted a combined experimental and theoretical approach to conduct research in the following six areas closely associated with SOPO tasks: 1) Development of electrocatalytically active bilayer oxygen electrodes (SOPO task-1) 2) Development of new symmetric three electrode cell (STEC) methodology to extract electrokinetic data of oxygen electrodes (SOPO task-2) 3) Quantification of electrokinetics of bilayer oxygen electrodes and correlation with degradation and delamination (SOPO task-2) 4) Performances of bilayer oxygen electrodes under fuel cells and electrolyzers modes (SOPO task-3) 5) Microscale modeling of oxygen electrode/electrolyte interface in solid oxide electrolysis cells (SOPO task-4) 6) Prediction of crack growth rate at oxygen electrode/electrolyte interface in solid oxide electrolysis cells (SOPO task-4)

08 HYDROGEN↗

NASA NTRSDate not supplied

Statistical Descriptors of Composite Fiber Aggregation

This study introduces a method of characterizing fiber aggregation and resin rich regions in composite microstructures. Microscale models of representative elements (RVE) need to be indicative of the extend of clustering (i.e. close fiber-to-fiber interaction) and resin rich “pools” which may impact the overall strength and performance of a composite structure. This algorithm was used to evaluate different unidirectional 2-D microstructure scans, which will be compared to their manufacturing method or any special treatment processes. These cluster and pool scan statistics can be used as criteria to judge statistical equivalency of artificially constructed microstructures.

Statistics↗

NASA NTRSDate not supplied

Addressing Critical Knowledge Gaps on Wildland Fires with UAS Technology

To better predict and respond to extreme fire behavior, wind shear, and superheated gases, there is a need for enhanced tactical microclimate wind forecasting. Collecting real-time or near-real-time three-dimensional atmospheric data during wildland fire suppression is vital for both ground firefighters and aviation safety. The use of balloons for soundings is not allowed due to aircraft operations, so new technology must be used. To address this need, the NASA FireSense Project has invested in co-developing and transitioning advanced technologies for atmospheric data collection to operational platforms to support decisions for wildland fire management. One of these technology investments has been with Uninhabited Aerial Systems (UAS) for atmospheric soundings. In collaboration with MITRE Corporation, a technology demonstration of Uninhabited Aerial Systems (UAS) for atmospheric soundings was held in Missoula, MT. This demonstration featured a NASA-designed payload on a Freefly Alta-X UAS, consistent with USFS UAS operations, balloon-borne soundings for data validation, and microscale modeling efforts. Results of this effort will be presented including comparison of forecasted conditions to data collection values as well as the impact on forecasts given real-time mixing heights and dispersion levels.

Jennifer Fowler↗

NASA NTRS2018

Modeling Fiber Kinking at the Microscale and Mesoscale

A computational micromechanics (CMM) model is employed to interrogate the assumptions of a recently developed mesoscale continuum damage mechanics (CDM) model for fiber kinking. The CMM model considers an individually discretized three dimensional fiber and surrounding matrix accounting for nonlinearity in the fiber, matrix plasticity, fiber/matrix interface debonding, and geometric nonlinearity. Key parameters of the CMM model were measured through experiments. In particular, a novel experimental technique to characterize the in situ longitudinal compressive strength of carbon fibers through indentation of micropillars is presented. The CDM model is formulated on the basis of Budiansky's fiber kinking theory (FKT) with a constitutive deformation-decomposition approach to alleviate mesh size sensitivity. In contrast to conventional mesoscale CDM models that prescribe a constitutive response directly, the response of the proposed model is an outcome of material nonlinearity and large rotations of the fiber direction following FKT. Comparison of the predictions from the CMM and CDM models shows remarkable correlation in strength, post-peak residual stress, and fiber rotation, with less than 10% difference between the two models in most cases. Additional comparisons are made with several fiber kinking models proposed in the literature to highlight the efficacy of the two models. Finally, the CMM model is exercised in parametric studies to explore opportunities to improve the longitudinal compression strength of a ply through the use of nonconventional microstructures.

Herraez, Miguel↗

DOE OSTI2021

Transport Modeling of Kinetically Limited Microscale Extraction Systems: Droplet and Supported Liquid Membrane Separations

Microfluidic droplet extraction devices can be used to efficiently measure solute phase transfer kinetics in liquid–liquid extraction systems. In this work, a droplet extraction apparatus was constructed using off-the-shelf and custom 3D printed parts for less than 110 USD. The interfacial kinetics of the transfer of Pu(IV) between an aqueous nitric acid phase and a 5 v/v% tridodecylamine (TDA) in toluene phase were characterized by determining the extent of Pu(IV) transfer between phases at different flow rates. Forward interfacial rate constants calculated from this data using an analytical model are 5.6 × 10 –3 mm s –1 under extraction (aqueous-to-organic solute transfer) conditions and 3.3 × 10 –3 mm s –1 under stripping (organic-to-aqueous solute transfer) conditions. The analytical model is valid when the rate of solute transfer across the interface between the two liquid phases has pseudo-first-order reversible reaction kinetics and is slow relative to the rate of solute transport to the interface. The change in the outlet concentration with residence time determined using the analytical model was compared with the results of a numerical model accounting for advection/diffusion to the interface in addition to kinetically limited interfacial solute transfer. These calculations revealed that for systems with forward interfacial rate constants greater than approximately 10 –2 mm s –1 , the analytical model must be supplemented with a numerical transport model accounting for an additional advection/diffusion mass transfer resistance. This threshold corresponds to a Damköhler number of one. The measured interfacial kinetic rate constants for Pu(IV)/TDA were used in a supported liquid membrane (SLM) numerical model, which showed that the Pu(IV) phase transfer kinetics are fast enough for use in a practical separation.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

NASA NTRSDate not supplied

In-Situ Scanning Electron Microscope Experiments for Microscale Mechanical Testing and Validated Modeling of Fiber Reinforced Thermoplastics

A novel, in-situ, scanning electron microscope (SEM) mechanical testing capability for materials at the microscale which provides experimental validation to a machine learning (ML) toolset for full-field validation of physics-based micromechanics models is being developed by researchers at NASA Glenn Research Center. These are enabling technologies for the integration of multiscale digital twins for materials into system level models which will result in the improved performance, material discovery, reduced production cost and time, rapid characterization, and prognostic structural health monitoring (SHM) for materials and structures for extreme environments in support of NASA space exploration missions. In order to bridge the material structure-to-system gap for digital twins, physics-based models must be experimentally validated at multiple length scales. Seminal microscale experiments, conducted at the Air Force Research Laboratory (AFRL), were limited to transverse compression of single-layer, unidirectional thermoset polymer matrix composite (PMC) micropillar specimens [1]. The early phases of the current project followed those initial results and setup to reproduce the compression testing of PMC material on the custom-built piezoelectric actuated micromechanical testing rig built by MicroTesting Solutions LLC. In this work, samples of thermoplastic PMC material were first machined into 3 mm cubes, and then further machining and final milling was done using a Focused Ion Beam (FIB). The initial experiment was done on a pillar roughly 20 µm x 20 µm x 40 µm tall. Additional pillars were milled with final sizes ranging from 20 µm x 20 µm x 40 µm tall to 40 µm x 40 µm x 65 µm tall. A speckle pattern for in-situ full-field measurements using Digital Image Correlation (DIC) was applied with platinum, which was coated on the surface, and then the FIB was used to mill away some of the coating to produce an irregular pattern of Pt on the pillar surface. The samples were loaded into the custom testing rig and placed into the SEM and loaded under compression until failure. Images were collected in the SEM during testing. Post-processing of the images was conducted using DIC to obtain full-field displacement and strain measurements elucidating the role of the matrix as well as fiber-fiber interaction at the microscale within the composite subjected to compression loading well into the non-linear regime of the material. Moreover, the evolution of fiber-matrix debonding and matrix cracking is observed in-situ at the microscale. This data, along with images segmented with a newly developed ML toolset [2], was used to create and validate physics-based micromechanics models. An image of the failed micropillar is shown in Figure 1. The techniques developed in the initial compression experiment was tailored to the validation needs of the models and expanded to include different sized samples as well as possibly tension and fatigue.

Laura Wilson↗

NASA NTRS1976

The excess fission xenon problem in lunar samples

This paper presents noble-gas (neon, krypton, and xenon) data obtained in stepwise heating of five lunar breccias: 14063, 14082, 15205, 15405, and 15445. These data are discussed in terms of spallation gas content, cosmic-ray exposure ages, and content of fission xenon. None of these samples is rich in solar-wind gases, and none contains more fission xenon than would be expected for in situ fission of U-238 and Pu-244. The essential features of the excess-fission-xenon phenomenon seen in several Apollo 14 samples are reviewed, as are models which might account for this phenomenon and criteria by which these models may be evaluated. The lack of excess fission xenon in 14063 and 14082 eliminates models relying on microscale xenon redistribution. A model is proposed in which the excess fission xenon originates in local material and percolates through and is adsorbed on regolith material.

Drozd, R. J.↗

NASA NTRS2006

Molecular-dynamics Simulation-based Cohesive Zone Representation of Intergranular Fracture Processes in Aluminum

A traction-displacement relationship that may be embedded into a cohesive zone model for microscale problems of intergranular fracture is extracted from atomistic molecular-dynamics simulations. A molecular-dynamics model for crack propagation under steady-state conditions is developed to analyze intergranular fracture along a flat 99 [1 1 0] symmetric tilt grain boundary in aluminum. Under hydrostatic tensile load, the simulation reveals asymmetric crack propagation in the two opposite directions along the grain boundary. In one direction, the crack propagates in a brittle manner by cleavage with very little or no dislocation emission, and in the other direction, the propagation is ductile through the mechanism of deformation twinning. This behavior is consistent with the Rice criterion for cleavage vs. dislocation blunting transition at the crack tip. The preference for twinning to dislocation slip is in agreement with the predictions of the Tadmor and Hai criterion. A comparison with finite element calculations shows that while the stress field around the brittle crack tip follows the expected elastic solution for the given boundary conditions of the model, the stress field around the twinning crack tip has a strong plastic contribution. Through the definition of a Cohesive-Zone-Volume-Element an atomistic analog to a continuum cohesive zone model element - the results from the molecular-dynamics simulation are recast to obtain an average continuum traction-displacement relationship to represent cohesive zone interaction along a characteristic length of the grain boundary interface for the cases of ductile and brittle decohesion. Keywords: Crack-tip plasticity; Cohesive zone model; Grain boundary decohesion; Intergranular fracture; Molecular-dynamics simulation

Yamakov, Vesselin I.↗

DOE OSTI2024

Microscale mechanical property variations of Al-substituted LLZO: insights from compression testing and molecular dynamics simulations

Ceramic solid electrolytes based on LLZO (Li 7 La 3 Zr 2 O 12 ) are promising candidates for all-solid-state batteries due to their high ionic conductivity and good apparent stability vs. lithium metal, however they are prone to mechanical failure. Lithium metal intrusions, alongside cell stack pressure, transition polycrystalline solid electrolyte grains into a compressed state that promotes crack propagation and fracture. Here this work examines the mechanical response of Al-substituted LLZO to compressive forces by measuring ultimate strength under pillar compression with a flat punch tip. Failure modes characterized by in situ scanning electron microscopy show diverse splitting patterns arising from internal porosity, grain boundaries, and slip planes. Large correlated variations in compressive strength (0.93–2.63 GPa) and Young's modulus (72.1–150.97 GPa) are observed across microscale regions of the solid electrolyte. Molecular dynamics simulations of LLZO with different porosities describe the variation of compressive strength and Young's modulus, and enable a microscale porosity model to be fit accounting for Young's modulus reduction across the solid electrolyte. Overall, the results indicate the importance of microscale mechanical testing of ceramic solid electrolytes to identify preferential sites for mechanical degradation and Li intrusion, and ensure the robust design of all-solid-state lithium metal batteries.

25 ENERGY STORAGE↗