Search NASA⌕ Search

SEARCH · Search NASA

Results for “Microstructure simulation”

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

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Simulation of Etch Pit Formation in DSMC Through Active Sites in Carbon Fiber Micro-Structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation because of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gases such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various detailed surface reaction mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗
NASA NTRSDate not supplied

Simulation of etch pit formation through active sites in carbon fiber micro-structures

Erosion of carbon due to oxidation does not occur uniformly but through localized etch pit formation as a result of active surface sites. In this work we demonstrate a newly developed capability to capture pitting of carbon fiber microstructures such as FiberForm, which is commonly used as the base material for NASA’s spacecraft ablative thermal protection systems (TPS). The simulations are performed at the meso-scale in order to capture the pit formation and growth using direct simulation Monte Carlo (DSMC), specifically using the Stochastic PArallel Rarefied-gas Time-accurate Analyzer (SPARTA) code. Legacy and latest models both assume uniform reactivity of carbon surface sites with oxygen even at the meso-scale level. However, in reality the carbon surface has locally different reaction rates due to the presence of defects at the atomic scale. The defective nature of these sites enhances their reactivity with atmospheric gases compared to the non-defective sites (2-3 orders of magnitude) and are termed as “active sites”. Thus, these sites tend to be the first to react and eventually get removed through the formation of gasses such as CO, CO2, and CN. Their removal results in all the neighboring atoms becoming defective, thus leading to chain reaction of localized carbon removal and formation of etch pits. Capturing the formation of pits during the ablation simulation of carbon micro-structures is critical to predicting their structural failure. Recently a detailed surface chemistry framework was developed in SPARTA, capable of incorporating various reaction mechanisms such as adsorption, desorption, Eley-Rideal (ER) and Langmuir-Hinshelwood (LH) mechanisms. We have implemented the capability of a single surface having multiple site sets with different reactivities within this framework. We have used this feature to model the presence of active sites on carbon surfaces, whose reactivity is orders of magnitude higher than that of the passive sites due to the presence of defects. The active site fraction is a property of surface elements within SPARTA and is directly proportional to the local reactivity of each surface element. By introducing an initial distribution of the active site fraction across the carbon surface and propagating it in a manner that mimics the evolution of real reacting carbon surfaces, we can capture the formation and growth of etch pits as a result of surface consumption reactions such as oxidation.

DSMC↗
NASA NTRS2011

Modeling and Characterization of Damage Processes in Metallic Materials

This paper describes a broad effort that is aimed at understanding the fundamental mechanisms of crack growth and using that understanding as a basis for designing materials and enabling predictions of fracture in materials and structures that have small characteristic dimensions. This area of research, herein referred to as Damage Science, emphasizes the length scale regimes of the nanoscale and the microscale for which analysis and characterization tools are being developed to predict the formation, propagation, and interaction of fundamental damage mechanisms. Examination of nanoscale processes requires atomistic and discrete dislocation plasticity simulations, while microscale processes can be examined using strain gradient plasticity, crystal plasticity and microstructure modeling methods. Concurrent and sequential multiscale modeling methods are being developed to analytically bridge between these length scales. Experimental methods for characterization and quantification of near-crack tip damage are also being developed. This paper focuses on several new methodologies in these areas and their application to understanding damage processes in polycrystalline metals. On-going and potential applications are also discussed.

Glaessgen, E. H.↗
NASA NTRS2015

Multiscale Modeling of Ceramic Matrix Composites

Results of multiscale modeling simulations of the nonlinear response of SiC/SiC ceramic matrix composites are reported, wherein the microstructure of the ceramic matrix is captured. This micro scale architecture, which contains free Si material as well as the SiC ceramic, is responsible for residual stresses that play an important role in the subsequent thermo-mechanical behavior of the SiC/SiC composite. Using the novel Multiscale Generalized Method of Cells recursive micromechanics theory, the microstructure of the matrix, as well as the microstructure of the composite (fiber and matrix) can be captured.

ceramic matrix composites↗
NASA NTRSDate not supplied

Modeling the Effective Elasticity of Anisotropic Porous Materials

The development and optimization of composite materials designed for thermal protection of NASA’s spacecraft require understanding their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the microscale on realistic representations of these composites. The digital microstructures are generated either synthetically or through X-ray micro-computed tomography reconstructions. One of the main challenges in the prediction of the structural response of heatshield materials is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the microscale. These are driven by the fiber layout within the microstructure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled using a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA) that was presented at the 10th Ablation Workshop. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases, and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Elasticity↗
NASA NTRSDate not supplied

Thermochemical Interactions of Yttria-Stabilized Zirconia and Molten Lunar Regolith Simulants

Oxygen produced from lunar resources through in-situ resource utilization (ISRU) is critical to maintaining a permanent human presence on the lunar surface. Molten regolith electrolysis and carbothermal reduction are two promising ISRU techniques for generating oxygen directly from lunar regolith, which is primarily a mixture of oxide minerals; however, both processes require operating temperatures of 1600C to melt lunar regolith and dissociate the molten oxides. These conditions limit the use of many oxide refractory materials, such as Al2O3 and MgO, due to rapid degradation resulting from reactions between the refractory materials and molten lunar regolith. Yttria-stabilized zirconia (YSZ) is a promising refractory oxide to provide containment of molten regolith while demonstrating limited reactivity. This work focuses on corrosion studies of YSZ powders and dense YSZ crucibles in contact with molten lunar mare and highlands regolith simulants at 1600C. The interactions between YSZ and molten regolith are characterized using SEM/EDS, XRD, and EBSD with an emphasis on elemental and microstructural analysis to assess reactivity and degradation of YSZ. Due to lunar regolith’s similar composition to calcium-magnesium-aluminosilicates (CMAS) and YSZ’s usage as a thermal barrier coating, these interactions can serve to inform YSZ/CMAS behavior by simulating cases of elevated CMAS/YSZ ratios and for higher than intended gas turbine temperatures.

Kevin Yu↗
NASA NTRSDate not supplied

Investigation of the Effect of Etch Pits on the Material Properties of Carbon Fiber Structures

It is well known that etch pits are formed on carbon surfaces during their consumption in reactive gas environments such as during ablation of carbon-based thermal protection system (TPS) materials. These pits start out from atomic defects on the carbon surface, grow in size and then ultimately lead to the complete removal of the solid carbon material. However, the effect of such etch pits on the material properties of carbon ablators such as FiberForm are poorly understood, as the maximum size of pits are on the order of microns in size, which makes them hard to observe experimentally. Hence, we have developed a module within the Porous Microstructure Analysis code PuMA that can generate etch pits on an arbitrary material surface, with a prescribed distribution of pit size and density. We use this newly developed module to calculate material properties of FiberForm, which is the base material for one of the most commonly used TPS materials - Phenolic Impregnated Carbon Ablator (PICA). The detailed microstructure of FiberForm in PuMA will either be generated synthetically or obtained from X-ray microtomography. Our simulations show that material properties such as the thermal conductivity and tortuosity decrease with increasing degree of pitting. This will be helpful to more accurately predict the degradation of carbon-based TPS during ablation, and furthermore may allow us to identify previously unknown TPS failure modes due to pitting.

PuMA↗
NASA NTRSDate not supplied

Unraveling Solar Wind Space Weathering of Carbon-Rich Asterods: Low-Flux VS. High-Flux Ion Irradiation of Murchison

Introduction: Space weathering processesalter the microstructural, compositional, and optical properties of airless planetary surfaces.Solar wind space weatheringforms ion-damaged rims within the outermost ~100 nm of regolith grains through atomic-displacements and implantation of energetic H+and He+ions[1]. Transmission electron microscope(TEM)studies of space weathered Itokawa and lunar olivine particles reveal predominantly nanocrystalline solar wind-damaged rims with only minor regions of amorphization[2-5]. This finding differs from laboratoryion irradiation experiments in which samples become amorphous at fluences 2-3 orders of magnitude lower than the independentlyconstrained fluences for still-crystalline natural samples [6-7].Ion flux is generally considered to have a secondaryeffectcompared to fluence in these irradiation experiments. However, solar wind irradiation of natural samples occurs at a flux that is ~4-5 orders of magnitude lower than whatisused inlaboratory simulations (e.g., [6,8-9]). Constraining how the very low ion flux of solar wind influences the microstructural, chemical, and spectral alteration of silicate minerals will improve our understanding of solar wind space weathering and aid in the analysis of returned samples from the Hayabusa2 and OSIRIS-REx missions. Here we report results from low-and high-flux ion irradiation experiments on the Murchison (CM2)meteoritewithcomparable total fluences. Methods: We performed low-flux and high-flux He+(4 keV) and H+(1 keV) irradiation experiments on separate, dry-cut Murchison slabs. Low-flux He+irradiation reached a total fluence of 2.1 x 1016ions/cm2(~400 years of exposure at Bennu) using a flux of 3.6 x 1011ions/cm2/s; low-flux H+irradiation reached a total fluence of 3.9 x 1016ions/cm2(~30 years) using a flux of 6.6 x 1011ions/cm2/s. High-flux He+irradiation reached a total fluence of 2.0 x 1016ions/cm2(~400 years)using a flux of 9.1 x 1012ions/cm2/s; high-flux H+irradiation reached a total fluence of 5.8 x 1016ions/cm2(~50 years)using flux of 8.4 x 1012ions/cm2/s. By irradiating to similar total fluences in both the low-flux and high-flux experiments, we can better isolate and evaluatethe effects of ion flux. Differences between unirradiatedand irradiated surfaces are determined using threeanalytical techniques:(1) in situX-ray photoelectron spectroscopy(XPS) reveals changes in surface chemistry,(2) visible to near-infrared spectroscopy (VNIR; 0.35-2.5 μm) highlights changes in spectraltrends, such as albedo and slope, and (3) TEM and energy-dispersive X-ray spectroscopy (EDS) identifieschanges in nanoscale structure and elemental composition within olivine and matrix electron-transparent thin sections prepared using focused ion beam scanning electron microscopy (FIB-SEM). Results and Discussion: XPS analyses of the low-fluxand high-fluxHe+-and H+-irradiated samples showa minorreduction in surface carbon contentand chemicalreduction ofFe3+to Fe2+. The latter result is consistent with the higher Fe2+/Fe3+ratios observed in phyllosilicate-rich space weathered rims of Ryugu grains [10].VNIR analyses suggest that the low-flux H+-irradiated spectrum has a slightly bluer slope (over 0.65-2.5 μm) compared to the unirradiated VNIR spectrum from [9] while the low-flux He+-irradiated spectral slope is comparable to the unirradiated surface. Both the high-flux He+-and H+-irradiated VNIR spectra exhibitslightlyhigheralbedos at shorter wavelengths relativeto their unirradiated counterparts, with He+irradiation yielding a greater change. Thus far, the low-flux He+-irradiated matrix and olivine FIBsections as well as the low-flux H+-irradiated olivine FIBsection have been examined with TEM. The surfaces of all three samplesshowless vesiculation than high-flux/high-fluenceirradiated matrix and olivine samples from [9]. The depth of phyllosilicate amorphization in thelow-fluxHe+-irradiated matrix FIBsection varieslaterally, with crystalline domains occurring as shallow as ~25 nm in some regions. Thelow-flux He+-irradiated olivine FIBsection has acontinuous,~40-90 nm thick ion-affected layer with polycrystalline and amorphous microstructures. The low-flux H+-irradiated olivine sampleexhibits a distinct, ~10-40 nm ion-affected layer whose uppermost ~20 nmin thicker regions is predominantly amorphous.Continuing to compare results from this work to those from [9] will provide a comprehensive understanding of how ion flux impacts the microstructural, chemical, and optical modification of various mineral phases and, in turn,help identify space weathering features from parent body alteration in Hayabusa2 and OSIRIS-REx returned samples.

D L Laczniak↗
NASA NTRS2007

Micromechanics of Spray-On Foam Insulation

Understanding the thermo-mechanical response of the Space Shuttle External Tank spray-on foam insulation (SOFI) material is critical, to NASA's Return to Flight effort. This closed-cell rigid polymeric foam is used to insulate the metallic Space Shuttle External Tank, which is at cryogenic temperatures immediately prior to and during lift off. The shedding of the SOFI during ascent led to the loss of the Columbia, and eliminating/minimizing foam lass from the tank has become a priority for NASA as it seeks to resume scheduled space shuttle missions. Determining the nature of the SOFI material behavior in response to both thermal and mechanical loading plays an important role as any structural modeling of the shedding phenomenon k predicated on knowledge of the constitutive behavior of the foam. In this paper, the SOFI material has been analyzed using the High-Fidelity Generalized Method of Cells (HFGMC) micromechanics model, which has recently been extended to admit a triply-periodic 3-D repeating unit cell (RUC). Additional theoretical extensions that mere made in order to enable modeling of the closed-cell-foam material include the ability to represent internal boundaries within the RUC (to simulated internal pores) and the ability to impose an internal pressure within the simulated pores. This latter extension is crucial as two sources contribute to significant internal pressure changes within the SOFI pores. First, gas trapped in the pores during the spray process will expand or contract due to temperature changes. Second, the pore pressure will increase due to outgassing of water and other species present in the foam skeleton polymer material. With HFGMC's new pore pressure modeling capabilities, a nonlinear pressure change within the simulated pore can be imposed that accounts for both of these sources, in addition to stmdar&-thermal and mechanical loading; The triply-periodic HFGMC micromechanics model described above was implemented within NASA GRC's MAC/GMC software package, giving the model access to a range of nonlinear constitutive models for the polymeric foam skeleton material. A repeating unit cell architecture was constructed that, while relatively simple, still accounts for the geometric anisotropy of the porous foam microstructure and its thin walls and thicker edges. With the lack of reliable polymeric foam skeleton materia1 properties, many simulations were executed aimed at backing out these material properties. Then, using these properties, predictions of the thermo-mechanical behavior of the foam, including calculated internal applied pressure profiles, were performed and compared with appropriate experimental data.

Bednarcyk, Brett A.↗
NASA NTRS2026

Controlled Parametric Forcing During Directional Solidification of a Bulk Organic Alloy Under Microgravity

The response of dendritic microstructures to step-like pulling velocity conditions is investigated using microgravity directional solidification experiments conducted on DECLIC-DSI combined with phase-field simulations. Under a constant pulling velocity of 1.5 µm/s, the evolution toward steady-state growth is characterized in terms of primary spacing, dendrite drift, and tip dynamics. For the first time, side-view observations enabled direct measurement of tip radius and sidebranching frequency. When step-like oscillations of the pulling velocity are imposed, the dendritic array exhibits a strongly period-dependent response: short periods lead to rapid tip adaptation, whereas longer periods induce a phase lag between tip position and morphology, resulting in progressive tip flattening and, above a critical period, interface destabilization and dendrite splitting. Quantitative phase-field simulations, including a realistic thermal field and stochastic noise, reproduce the experimental observations and provide insight into the governing mechanisms, highlighting the role of characteristic relaxation times, sequence-dependent effects, and the irreversible reorganization of the microstructure following splitting.

Microgravity↗
NASA NTRSDate not supplied

Investigation of the effect of etch pits on the material properties of carbon fiber structures

It is well known that etch pits are formed on carbon surfaces during their consumption in reactive gas environments such as during ablation of carbon-based thermal protection system (TPS) materials. These pits start out from atomic defects on the carbon surface, grow in size and then ultimately lead to the complete removal of the solid carbon material [1, 2]. However, the effect of such etch pits on the material properties of carbon ablators such as FiberForm are poorly understood, as the maximum size of pits are on the order of microns in size, which makes them hard to observe experimentally. Hence, we have developed a module within the Porous Microstructure Analysis code PuMA [3] that can generate etch pits on an arbitrary material surface, with a prescribed distribution of pit size and density. Figure 1 shows a FiberForm sample with artificially generated pits of varying radii and pit density. We use this newly developed module to calculate material properties of FiberForm, which is the base material for one of the most commonly used TPS materials - Phenolic Impregnated Carbon Ablator (PICA). The detailed microstructure of FiberForm in PuMA will either be generated synthetically or obtained from X-ray microtomography [4]. Our simulations show that material properties such as the thermal conductivity and tortuosity decrease with increasing degree of pitting. This will be helpful to more accurately predict the degradation of carbon-based TPS during ablation, and furthermore may allow us to identify previously unknown TPS failure modes due to pitting.

PuMA↗
NASA NTRS2017

Tortuosity Computations of Porous Materials using the Direct Simulation Monte Carlo

Low-density carbon fiber preforms, used as thermal protection systems (TPS) materials for planetary entry systems, have permeable, highly porous microstructures consisting of interlaced fibers. Internal gas transport in TPS is important in modeling the penetration of hot boundary-layer gases and the in-depth transport of pyrolysis and ablation products. The gas effective diffusion coefficient of a porous material must be known before the gas transport can be modeled in material response solvers; however, there are very little available data for rigid fibrous insulators used in heritage TPS.The tortuosity factor, which reflects the efficiency of the percolation paths, can be computed from the effective diffusion coefficient of a gas inside a porous material and is based on the micro-structure of the material. It is well known, that the tortuosity factor is a strong function of the Knudsen number. Due to the small characteristic scales of porous media used in TPS applications (typical pore size of the order of 50 micron), the transport of gases can occur in the rarefied and transitional regimes, at Knudsen numbers above 1. A proper way to model the gas dynamics at these conditions consists in solving the Boltzmann equation using particle-based methods that account for movement and collisions of atoms and molecules.In this work we adopt, for the first time, the Direct Simulation Monte Carlo (DSMC) method to compute the tortuosity factor of fibrous media in the rarefied regime. To enable realistic simulations of the actual transport of gases in the porous medium, digitized computational grids are obtained from X-ray micro-tomography imaging of real TPS materials. The SPARTA DSMC solver is used for simulations. Effective diffusion coefficients and tortuosity factors are obtained by computing the mean-square displacement of diffusing particles.We first apply the method to compute the tortuosity factors as a function of the Knudsen number for computationally designed materials such as random cylindrical fibers and packed bed of spheres with prescribed porosity. Results are compared to literature values obtained using random walk methods in the rarefied and transitional regime and a finite-volume method for the continuum regime. We then compute tortuosity factors for a real carbon fiber material with a transverse isotropic structure (FiberForm), quantifying differences between through-thickness and in-plain tortuosities at various Knudsen regimes.

Tortuosity↗
NASA NTRS1993

A study of subsurface crack initiation produced by rolling contact fatigue

Results of subsurface crack initiation studies produced by pure rolling contact fatigue in 7075-T6 aluminum alloy are presented in this article. Microstructural changes and subsequent crack initiation below the contacting surface in cylindrical test specimens subjected to repeated rolling contact are illustrated. The rolling conditions are simulated in a three-dimensional elastic-plastic finite element model in order to estimate the plastic strains and residual stresses in the test material. The numerically estimated distribution of plastic strains in the model correlate well with the extent of microstructural changes observed in the test specimen. Results also indicate that a combination of plastic strains and low values of residual stresses is conducive to subsurface crack initiation and growth.

Kumar, Arun M.↗
NASA NTRS1991

Effect of a prior stretch on the aging response of an Al-Cu-Li-Ag-Mg-Zr alloy

Recently, a family of Al-Cu-Li alloys containing minor amounts of Ag, Mg, and Zr and having desirable combinations of strength and toughness were developed. The Weldalite (trademark) alloys exhibit a unique characteristic in that with or without a prior stretch, they obtain significant strength-ductility combinations upon natural and artificial aging. The ultra-high strength (approximately 690 MPa yield strength) in the peak-aged tempers (T6 and T8) were primarily attributed to the extremely fine T(sub 1) (Al2CuLi) or T(sub 1)-type precipitates that occur in these alloys during artificial aging, whereas the significant natural aging response observed is attributed to strengthening from delta prime (Al3Li) and GP zones. In recent work, the aging behavior of an Al-Cu-Li-Ag-Mg alloy without a prior stretch was followed microstructurally from the T4 to the T6 condition. Commercial extrusions, rolled plates, and sheets of Al-Cu-Li alloys are typically subjected to a stretching operation before artificial aging to straighten the extrusions and, more importantly, introduce dislocations to simulate precipitation of strengthening phases such as T(sub 1) by providing relatively low-energy nucleation sites. The goals of this study are to examine the microstructure that evolves during aging of an alloy that was stretch after solution treatment and to compare the observations with those for the unstretched alloy.

Kumar, K. S.↗
NASA NTRS2026

Modeling Study of Hatch Spacing’s Effect on Grain Morphology in Repairing Damaged Metal Parts with Welding and Hatch Spacing’s Potential Use in Lunar Exploration

The welding process is a potential way of repairing a damaged metal component, especially cavity damage caused by a harsh environment like the lunar environment, which is characterized by large temperature differences and reduced gravity. The adjustment of welding parameter (e.g., hatch spacing) can improve production efficiency in the repair process. Seen from the microstructural level, hatch spacing sensitivity affects the metallic grain evolution and morphology in the welding process, which can further influence a repaired part’s mechanical properties; however, the study of hatch spacing’s effect on microstructure is challenging. Traditional experimental procedures are costly and time-consuming, and any change in hatch spacing value needs roll-back of experimental procedure. A modeling study can address the above challenges in experimental observation. In this research, a modeling approach based on the Kinetic Monte Carlo (KMC) Potts theory was used to simulate grain evolution and morphology with three hatch spacings. Through quantifying and analyzing the predicted grain morphologies, the effect of hatch spacing on microstructure in a welding-fabricated part was investigated. The predicted grain morphologies were validated with an EBSD image of welding microstructure, which has been published before. The primary grain morphologies were columnar grains with a small amount of fine equiaxed grains formed in the scanning path centerline. When increasing the hatch spacing, the columnar grains become larger and more lengthy, while the effect of hatch spacing on the equiaxed grains is not obvious.

Welding for repairing↗
NASA NTRS2015

Microstructural Influence on Deformation and Fatigue Life of Composites Using the Generalized Method of Cells

A fully coupled deformation and damage approach to modeling the response of composite materials and composite laminates is presented. It is based on the semi-­‐analytical generalized method of cells (GMC) micromechanics model as well as its higher fidelity counterpart, HFGMC, both of which provide closed-form constitutive equations for composite materials as well as the micro scale stress and strain fields in the composite phases. The provided constitutive equations allow GMC and HFGMC to function within a higher scale structural analysis (e.g., finite element analysis or lamination theory) to represent a composite material point, while the availability of the micro fields allow the incorporation of lower scale sub­‐models to represent local phenomena in the fiber and matrix. Further, GMC's formulation performs averaging when applying certain governing equations such that some degree of microscale field accuracy is surrendered in favor of extreme computational efficiency, rendering the method quite attractive as the centerpiece in a integrated computational material engineering (ICME) structural analysis; whereas HFGMC retains this microscale field accuracy, but at the price of significantly slower computational speed. Herein, the sensitivity of deformation and the fatigue life of graphite/epoxy PMC composites, with both ordered and disordered microstructures, has been investigated using this coupled deformation and damage micromechanics based approach. The local effects of fiber breakage and fatigue damage are included as sub‐models that operate on the microscale for the individual composite phases. For analysis of laminates, classical lamination theory is employed as the global or structural scale model, while GMC/HFGMC is embedded to operate on the microscale to simulate the behavior of the composite material within each laminate layer. A key outcome of this study is the statistical influence of microstructure and micromechanics idealization (GMC or HFGMC) on the overall accuracy of unidirectional and laminated composite deformation and fatigue response.

Microscale Models↗
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↗