Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Microscale modeling”

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 19 records
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Bridging the Gap Between Microscale Modeling and Additive Manufacturing for TPS

An overview of the now open-source NASA ARC software PuMA will be provided. On top of the well-documented ability to import and compute material properties from micro-CT images, PuMA has the ability to design microstructures and compute/predict material properties such as: porosity, permeability, tortuosity, thermal and electrical conductivity, tensile strength, etc. Therefore, a whole new range of capabilities is available, and provides users with the ability to build TPS materials (both fibrous and woven types) and optimize their properties based on missions and requirements. This talk will aim to provide a seed to start bridging the gap between microscale modeling of TPS materials and additive manufacturing.

Conductivity↗
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↗
NASA NTRS1977

The multiple nesting of mesoscale, submesoscale, and nonhydrostatic microscale numerical models for a case of mesocyclogenesis and severe storm development

A fundamental approach involving the integration of two different numerical models over several different spatial mesh lengths is considered. Model I is a hydrostatic, incompressible, and adiabatic primitive (Newtonian) equation set with no boundary layer forcing, i.e., surface heat or moisture flux or external frictional stresses. Typically this model is initialized from conventional radiosonde data and run for a 12-15 hour forecast period over a 40 km horizontal mesh. If phase I of the tornado developmental processes develops, the data is 'frozen' in space and time, interpolated to a 20 km mesh, and integrated for a shorter time period (3 hours) with a shorter time step. Model II is a nonhydrostatic compressible 18-level complement of model I. If model II indicates the development of the intensification of phase V dynamical processes, the data is frozen and interpolated to finer mesh lengths until the tornadic circulation system is approximated

Kaplan, M. L.↗
NASA NTRS2015

Multiscale and Multiphysics Modeling of Additive Manufacturing of Advanced Materials

The objective of this proposed project is to research and develop a prediction tool for advanced additive manufacturing (AAM) processes for advanced materials and develop experimental methods to provide fundamental properties and establish validation data. Aircraft structures and engines demand materials that are stronger, useable at much higher temperatures, provide less acoustic transmission, and enable more aeroelastic tailoring than those currently used. Significant improvements in properties can only be achieved by processing the materials under nonequilibrium conditions, such as AAM processes. AAM processes encompass a class of processes that use a focused heat source to create a melt pool on a substrate. Examples include Electron Beam Freeform Fabrication and Direct Metal Deposition. These types of additive processes enable fabrication of parts directly from CAD drawings. To achieve the desired material properties and geometries of the final structure, assessing the impact of process parameters and predicting optimized conditions with numerical modeling as an effective prediction tool is necessary. The targets for the processing are multiple and at different spatial scales, and the physical phenomena associated occur in multiphysics and multiscale. In this project, the research work has been developed to model AAM processes in a multiscale and multiphysics approach. A macroscale model was developed to investigate the residual stresses and distortion in AAM processes. A sequentially coupled, thermomechanical, finite element model was developed and validated experimentally. The results showed the temperature distribution, residual stress, and deformation within the formed deposits and substrates. A mesoscale model was developed to include heat transfer, phase change with mushy zone, incompressible free surface flow, solute redistribution, and surface tension. Because of excessive computing time needed, a parallel computing approach was also tested. In addition, after investigating various methods, a Smoothed Particle Hydrodynamics Model (SPH Model) was developed to model wire feeding process. Its computational efficiency and simple architecture makes it more robust and flexible than other models. More research on material properties may be needed to realistically model the AAM processes. A microscale model was developed to investigate heterogeneous nucleation, dendritic grain growth, epitaxial growth of columnar grains, columnar-to-equiaxed transition, grain transport in melt, and other properties. The orientations of the columnar grains were almost perpendicular to the laser motion's direction. Compared to the similar studies in the literature, the multiple grain morphology modeling result is in the same order of magnitude as optical morphologies in the experiment. Experimental work was conducted to validate different models. An infrared camera was incorporated as a process monitoring and validating tool to identify the solidus and mushy zones during deposition. The images were successfully processed to identify these regions. This research project has investigated multiscale and multiphysics of the complex AAM processes thus leading to advanced understanding of these processes. The project has also developed several modeling tools and experimental validation tools that will be very critical in the future of AAM process qualification and certification.

Liou, Frank↗
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↗
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↗
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

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

The Effect of Roughness Model on Scattering Properties of Ice Crystals.

We compare stochastic models of microscale surface roughness assuming uniform and Weibull distributions of crystal facet tilt angles to calculate scattering by roughened hexagonal ice crystals using the geometric optics (GO) approximation. Both distributions are determined by similar roughness parameters, while the Weibull model depends on the additional shape parameter. Calculations were performed for two visible wavelengths (864 nm and 410 nm) for roughness values between 0.2 and 0.7 and Weibull shape parameters between 0 and 1.0 for crystals with aspect ratios of 0.21, 1 and 4.8. For this range of parameters we find that, for a given roughness level, varying the Weibull shape parameter can change the asymmetry parameter by up to about 0.05. The largest effect of the shape parameter variation on the phase function is found in the backscattering region, while the degree of linear polarization is most affected at the side-scattering angles. For high roughness, scattering properties calculated using the uniform and Weibull models are in relatively close agreement for a given roughness parameter, especially when a Weibull shape parameter of 0.75 is used. For smaller roughness values, a shape parameter close to unity provides a better agreement. Notable differences are observed in the phase function over the scattering angle range from 5deg to 20deg, where the uniform roughness model produces a plateau while the Weibull model does not.

Weibull density functions↗