Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “materials modelling”

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

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

At least 1,279 records · Page 71

Oxidation of Reinforced Carbon-Carbon Subjected to Hypervelocity Impact

This paper presents results from arc jet tests conducted at the NASA Johnson Space Center on reinforced carbon-carbon (RCC) samples subjected to hypervelocity impact. The RCC test specimens are representative of RCC components used on the Space Shuttle Orbiter. The arc jet testing established the oxidation characteristics of RCC when hypervelocity projectiles, simulating meteoroid/orbital debris, impact the RCC material. In addition to developing correlations for use in trajectory simulations, we discuss analytical modeling of the increased material oxidation in the impacted area using measured hole growth data. Entry flight simulations are useful in assessing the increased Space Shuttle RCC component degradation as a result of impact damage and the hot gas flow through an enlarging hole into the wing leading-edge cavity.

Curry, Donald M.↗

Simulation Model of A Ferroelectric Field Effect Transistor

An electronic simulation model has been developed of a ferroelectric field effect transistor (FFET). This model can be used in standard electrical circuit simulation programs to simulate the main characteristics of the FFET. The model uses a previously developed algorithm that incorporates partial polarization as a basis for the design. The model has the main characteristics of the FFET, which are the current hysterisis with different gate voltages and decay of the drain current when the gate voltage is off. The drain current has values matching actual FFET's, which were measured experimentally. The input and output resistance in the model is similar to that of the FFET. The model is valid for all frequencies below RF levels. A variety of different ferroelectric material characteristics can be modeled. The model can be used to design circuits using FFET'S with standard electrical simulation packages. The circuit can be used in designing non-volatile memory circuits and logic circuits and is compatible with all SPICE based circuit analysis programs. The model is a drop in library that integrates seamlessly into a SPICE simulation. A comparison is made between the model and experimental data measured from an actual FFET.

MacLeod, Todd C.↗

NASA ACC High Energy Dynamic Impact Methodology and Outcomes

High fidelity analysis methods known as Progressive Damage Analysis (PDA) methods, capable of reliably predicting the onset and progression of damage in composite materials, are being developed for HEDI event simulation. Four state-of-the-art PDA methods are being investigated: LS-DYNA MAT162, LS-DYNA MAT261, Smooth-Particle Galerkin (SPG), and EMU Peridynamics. Validation of these PDA models aims to follow a building block approach, starting with coupon- and component-level development. As part of the overall effort, material characterization testing was performed to bridge gaps between existing experimental data and the material property inputs required to predict ballistic impact behavior at the component-level. The material models were updated with the results of the coupon-level testing and were then used to predict panel behavior and damage in ballistic impact events. To assess the accuracy of the PDA methods, these pre-test predictions were compared to the results of ballistic impact testing.

Hunziker, Kenneth J.↗

Origami-Inspired Folding of Thick, Rigid Panels

To achieve power of 250 kW or greater, a large compression ratio of stowed-to-deployed area is needed. Origami folding patterns were used to inspire the folding of a solar array to achieve synchronous deployment; however, origami models are generally created for near-zero-thickness material. Panel thickness is one of the main challenges of origami-inspired design. Three origami-inspired folding techniques (flasher, square twist, and map fold) were created with rigid panels and hinges. Hinge components are added to the model to enable folding of thick, rigid materials. Origami models are created assuming zero (or near zero) thickness. When a material with finite thickness is used, the panels are required to bend around an increasingly thick fold as they move away from the center of the model. The two approaches for dealing with material thickness are to use membrane hinges to connect the panels, or to add panel hinges, or hinges of the same thickness, at an appropriate width to enable folding.

Trease, Brian P.↗

Modeling the Relationship Between Porosity and Permeability During Oxidation of Ablative Materials

The ablative materials used in thermal protection systems (TPS) undergo oxidation during atmospheric entry which leads to an in-depth change in both permeability and porosity. These properties have a significant affect on heat transfer in a TPS during entry. X-ray micro-tomography has provided 3D images capturing the micro-structure of TPS materials. In this study, we use micro-tomography based simulations to create high-fidelity models relating permeability to porosity during oxidation of FiberForm, the carbon fiber preform of the Phenolic Impregnated Carbon Ablator (PICA) often used as a TPS material. The goal of this study is to inform full-scale models and reduce uncertainty in TPS modeling.

Ablativ↗

Computational Simulation of the Formation and Material Behavior of Ice

Computational methods are described for simulating the formation and the material behavior of ice in prevailing transient environments. The methodology developed at the NASA Lewis Research Center was adopted. A three dimensional finite-element heat transfer analyzer was used to predict the thickness of ice formed under prevailing environmental conditions. A multi-factor interaction model for simulating the material behavior of time-variant ice layers is presented. The model, used in conjunction with laminated composite mechanics, updates the material properties of an ice block as its thickness increases with time. A sample case of ice formation in a body of water was used to demonstrate the methodology. The results showed that the formation and the material behavior of ice can be computationally simulated using the available composites technology.

Tong, Michael T.↗

Organic Solid Matter as a Coloring Agent in Outer Solar System Bodies

Small bodies in the outer Solar System OSS, exhibit a range of color, or slope of the reflectance in the photovisual spectral region, ranging from neutral to very red, sometimes with and sometimes without distinct absorption bands. These objects range in geometric albedo from 0.03 to 1.0, with the higher albedo objects typically showing clear evidence of water ice. Water ice has also been found in a few objects with albedo 0. 1 or less. We explore here the identification of the material or materials that color these icy and non-icy surfaces through scattering models that incorporate minerals, meteoritic material, and organic solids (tholins) produced ID the laboratory by energy deposition in ices and gases. These models must match not only the color in the photovisual region, but the spectral reflectance properties throughout the near-infrared. Among some classes of objects, such as Kuiper Belt objects, the coloring agent may be a single material that is present in greater or lesser abundance, thus accounting for the range in color from neutral to very red. This may also apply to the Centaur objects, the Jovian Trojans, and the outer-main belt asteroids, each taken as a separate class. If so, each class may be colored to varying degrees by a different material, or they all might be colored by a common material that is widespread throughout the OSS, from 3 to 50 AU, and beyond. In this paper, we model the reflectances of "Kuiper Belt objects, Centaurs, Trojans, outer ARAB asteroids, and planetary satellites. Our models show that the reddest surfaces cannot be colored by minerals or meteoritic materials, but can be matched throughout the photovisual and near-infrared by organic solids, specifically certain tholins.

Cruikshank, D. P.↗

An Update to the NASA Reference Solar Sail Thrust Model

An optical model of solar sail material originally derived at JPL in 1978 has since served as the de facto standard for NASA and other solar sail researchers. The optical model includes terms for specular and diffuse reflection, thermal emission, and non-Lambertian diffuse reflection. The standard coefficients for these terms are based on tests of 2.5 micrometer Kapton sail material coated with 100 nm of aluminum on the front side and chromium on the back side. The original derivation of these coefficients was documented in an internal JPL technical memorandum that is no longer available. Additionally more recent optical testing has taken place and different materials have been used or are under consideration by various researchers for solar sails. Here, where possible, we re-derive the optical coefficients from the 1978 model and update them to accommodate newer test results and sail material. The source of the commonly used value for the front side non-Lambertian coefficient is not clear, so we investigate that coefficient in detail. Although this research is primarily designed to support the upcoming NASA NEA Scout and Lunar Flashlight solar sail missions, the results are also of interest to the wider solar sail community.

Heaton, Andrew F.↗

Effects of elevated temperature on the viscoplastic modeling of graphite/polymeric composites

To support the development of new materials for the design of next generation supersonic transports, a research program is underway at NASA to assess the long term durability of advanced polymer matrix composites (PMC's). One of main objectives of the program was to explore the effects of elevated temperature (23 to 200 C) on the constitutive model's material parameters. To achieve this goal, test data on the observed nonlinear, stress-strain behavior of IM7/5260 and IM7/8320 composites under tension and compression loading were collected and correlated against temperature. These tests, conducted under isothermal conditions using variable strain rates, included such phenomena as stress relaxation and short term creep. The second major goal was the verification of the model by comparison of analytical predictions and test results for off axis and angle ply laminates. Correlation between test and predicted behavior was performed for specimens of both material systems over a range of temperatures. Results indicated that the model provided reasonable predictions of material behavior in load or strain controlled tests. Periods of loading, unloading, stress relaxation, and creep were accounted for.

Gates, Thomas S.↗

Prediction of Stiffness and Fatigue Lives of Polymer Matrix Composite Laminates Using Artificial Neural Networks

Machine learning (ML) models are increasingly being used in many engineering fields due to the advancements in ML algorithms and availability of high-speed computing power. One of the most popular ML class of models is artificial neural networks (ANN). ML is increasingly being used in the design and analysis of composite materials and structures, specifically in the constitutive modeling of composite materials with the focus on greatly accelerating multiscale analyses of composite materials and structures through development of surrogate models. Towards that end, both Python and MATLAB-based neural nets have been developed to predict initial stiffness and fatigue life of an eight-ply symmetric polymer matrix composite laminate. Two types of neural networks, a Multilayer Perceptron (MLP) and a Recurrent Neural Network (RNN), have been developed for both platforms. Results show that the both neural net types can provide an excellent estimate of initial stiffness as well as fatigue life of eight-ply symmetric polymer matrix composite laminate. RNNs are better able to capture the shape of the fatigue curve of a laminate. This tool can be very useful for system level studies to obtain an estimate of desired properties and life of PMC composite laminates. The associated surrogate models could also be used in composite multiscale analyses to replace the actual physics-based calculations at lower scales and thereby significantly increase the computational efficiency of such analyses and thus make multiscale analyses a viable industrial tool for large scale structural problems.

Composite↗

Stiffness and Fatigue Life Estimator for Polymer Composite Laminates Using Machine Learning

Machine learning (ML) models are increasingly being used in many engineering fields due to the advancements in ML algorithms and availability of high-speed computing power. One of the most popular ML class of models is artificial neural networks (ANN). ML is increasingly being used in the design and analysis of composite materials and structures, specifically in the constitutive modeling of composite materials with the focus on greatly accelerating multiscale analyses of composite materials and structures through development of surrogate models. Towards that end, Python-based neural nets have been developed to predict initial stiffness and fatigue life of an eight-ply symmetric polymer matrix composite laminate. Two types of neural networks, a Multilayer Perceptron (MLP) and a Recurrent Neural Network (RNN), have been established. Results show that both neural net type algorithms can provide an excellent estimate of initial laminate stiffness as well as fatigue life of eight-ply symmetric polymer matrix composite laminates (PMCs). RNNs are better able to capture the shape of the fatigue curve of a laminate. The resulting tool and GUI can be very useful for system level studies to obtain an estimate of desired properties and life of PMC composite laminates. Further, the associated surrogate models can also be used in composite multiscale analyses to replace the actual physics-based calculations at lower scales and thereby significantly increase the computational efficiency of such analyses and thus make micromechanics-based multiscale analyses a viable industrial tool for large scale structural problems.

multiscale analysis↗

Stiffness and Fatigue Life Estimator for Polymer Composite Laminates Using Machine Learning

Machine learning (ML) models are increasingly being used in many engineering fields due to the advancements in ML algorithms and availability of high-speed computing power. One of the most popular ML class of models is artificial neural networks (ANN). ML is increasingly being used in the design and analysis of composite materials and structures, specifically in the constitutive modeling of composite materials with the focus on greatly accelerating multiscale analyses of composite materials and structures through development of surrogate models. Towards that end, Python-based neural nets have been developed to predict initial stiffness and fatigue life of an eight-ply symmetric polymer matrix composite laminate. Two types of neural networks, a Multilayer Perceptron (MLP) and a Recurrent Neural Network (RNN), have been established. Results show that both neural net type algorithms can provide an excellent estimate of initial laminate stiffness as well as fatigue life of eight-ply symmetric polymer matrix composite laminates (PMCs). RNNs are better able to capture the shape of the fatigue curve of a laminate. The resulting tool and GUI can be very useful for system level studies to obtain an estimate of desired properties and life of PMC composite laminates. Further, the associated surrogate models can also be used in composite multiscale analyses to replace the actual physics-based calculations at lower scales and thereby significantly increase the computational efficiency of such analyses and thus make micromechanics-based multiscale analyses a viable industrial tool for large scale structural problems.

multiscale analysis↗

Uncertainty Quantification of Classical Theories of Dendritic Growth Kinetics Applied to Nickel-Based Alloys

The solidification velocity in a model nickel-alloy single crystal during laser spot melting was recently characterized using synchrotron X-ray imaging. The measured solidification velocity was found to exceed the absolute stability threshold predicted by the Kurz-Giovanola-Trivedi (KGT) model. The discrepancies between the model and experiments motivate the further assessment of accurate material properties. This work quantifies the impact material property uncertainty has on model predictions of the absolute stability threshold velocities. Properties from the literature are reviewed and compared to those calculated using computational thermodynamics to provide uncertainty estimates on input properties to the KGT model. Global sensitivity analysis is used to quantify the influence of each uncertain input property on the predicted threshold velocity. This work supports the understanding of the nickel-alloy solidification during powder bed fusion additive manufacturing and identifies the solidification material properties that are the most important to assess from first-principles computations and experiments.

Computational thermodynamics↗

Thermal performance of 625-kg/cu m elastomeric ablative materials on spherically blunted 0.44-radian cones

Spherically blunted 0.44-radian (25 deg) half-angle conical models coated with elastomeric ablative materials were tested in supersonic arc-heated wind tunnels to evaluate performance of the ablators over a range of conditions typical of lifting entry. Four test conditions were combinations of stagnation point-heat transfer rates of 2.3 and 4.5 MW/m2 and stagnation pressures of 20 and 2kN/m2. Afterbody values of heat transfer rate and pressure were 0.05 to 0.20 of stagnation point values. Stagnation enthalpy varied from 4.4 to 25 MJ/kg (1900 to 11000 Btu/lbm) and free-stream Mach number was in a range from 3.5 to 4. Ablative materials retained the spherical nose shape throughout tests at the lower heat transfer level, but receded, assuming a flattened nose shape, during tests at the high heat transfer level. The residue layer that formed on the conical after-body was weak, friable, and extensively cracked. The reference ablative material, which contained phenolic microspheres, generally retained the conical shape on the model afterbody. However, a modified ablator, in which phenolic microspheres were replaced with silica microspheres, deformed and separated from the undegraded material, and thereby produced a very uneven surface. Substrate temperatures and ablator recession were in good agreement with values computed by a numerical analysis.

Champman, A. J.↗

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

silicone intumescence↗

Detection of Subsurface Material Separation in Shuttle Orbiter Slip-Side Joggle Region of the Wing Leading Edge using Infrared Imaging Data from Arc Jet Tests

The objective of the present study was to determine whether infrared imaging (IR) surface temperature data obtained during arc-jet tests of Space Shuttle Orbiter s reinforced carbon-carbon (RCC) wing leading edge panel slip-side joggle region could be used to detect presence of subsurface material separation, and if so, to determine when separation occurs during the simulated entry profile. Recent thermostructural studies have indicated thermally induced interlaminar normal stress concentrations at the substrate/coating interface in the curved joggle region can result in local subsurface material separation, with the separation predicted to occur during approach to peak heating during reentry. The present study was an attempt to determine experimentally when subsurface material separations occur. A simplified thermal model of a flat RCC panel with subsurface material separation was developed and used to infer general surface temperature trends due to the presence of subsurface material separation. IR data from previously conducted arc-jet tests on three test specimens were analyzed: one without subsurface material separation either pre or post test, one with pre test separation, and one with separation developing during test. The simplified thermal model trend predictions along with comparison of experimental IR data of the three test specimens were used to successfully infer material separation from the arc-jet test data. Furthermore, for the test specimen that had developed subsurface material separation during the arc-jet tests, the initiation of separation appeared to occur during the ramp up to the peak heating condition, where test specimen temperature went from 2500 to 2800 F.

Daryabeigi, Kamran↗

Testing the Predictive Capability of the High-Fidelity Generalized Method of Cells Using an Efficient Reformulation

The High-Fidelity Generalized Method of Cells is a new micromechanics model for unidirectionally reinforced periodic multiphase materials that was developed to overcome the original model's shortcomings. The high-fidelity version predicts the local stress and strain fields with dramatically greater accuracy relative to the original model through the use of a better displacement field representation. Herein, we test the high-fidelity model's predictive capability in estimating the elastic moduli of periodic composites characterized by repeating unit cells obtained by rotation of an infinite square fiber array through an angle about the fiber axis. Such repeating unit cells may contain a few or many fibers, depending on the rotation angle. In order to analyze such multi-inclusion repeating unit cells efficiently, the high-fidelity micromechanics model's framework is reformulated using the local/global stiffness matrix approach. The excellent agreement with the corresponding results obtained from the standard transformation equations confirms the new model's predictive capability for periodic composites characterized by multi-inclusion repeating unit cells lacking planes of material symmetry. Comparison of the effective moduli and local stress fields with the corresponding results obtained from the original Generalized Method of Cells dramatically highlights the original model's shortcomings for certain classes of unidirectional composites.

Arnold, Steven M.↗

Thermo-Poro-Mechanical Modeling of RTV Intumescence

Room temperature vulcanizing (RTV) silicone is a high-temperature adhesive used as a gap-filler between heatshield tiles in numerous entry missions. Its propensity to intumesce, or swell upon exposure to heat, is a well-known effect that needs to be carefully quantified during design. At tile interfaces of charring ablators, intumescence, combined with differential recession, could cause the gap filler to protrude past the ablator outer mold line, forming a “fence”. Fencing can in turn cause transition to turbulence of the flow wetting the heat shield, leading to augmented surface heating. Recent experiments conducted at the Plasmatron X facility, the high enthalpy wind tunnel of the Center for Hypersonics and Entry Systems Studies, have shown prominent fencing of RTV gap fillers in PICA, under both nitrogen and air plasmas. Similar observations are well known in the arcjet literature. Further experiments under controlled environment, performed using in situ X-ray micro-computed tomography (micro-CT) at the Advanced Light Source (ALS), have shown heating rate-dependent swelling and shrinkage of RTV during pyrolysis. To simulate RTV intumescence, a novel model was introduced in the Porous Materials Analysis Toolbox based on OpenFOAM, PATO, to account for pore-pressure buildup within both closed- and open-pores. The governing equation for the thermo-poro-mechanical response were developed, assuming linear elasticity for the charring silicone. A new multi-pyrolysis model that tracks non-monotonic advancement of material properties with pyrolysis was proposed. This model addresses the limitations of state-of-the-art ablator models to capture the different stages of thermal degradation and coupled thermomechanics. Swelling of RTV was simulated using the new thermo-poro-mechanical model and compared against in situ micro-CT data. Results showed good agreement in intumescence height and temperature profiles at all heating rates, indicating that the key factor contributing to RTV swelling is the internal pressure build-up within closed- and open-pores. As RTV is cured into a soft (rubbery) compound with low-porosity and permeability, initial temperature increase and pyrolysis gas production cause a significant increase of internal pressure, causing a pronounced volume growth. As thermal degradation progresses, rigidization of the silicone occurs due to char hardening which counteract volume shrinkage after gas pressure relief. Overall, our model shows that accounting for changes in properties such porosity, permeability and key thermomechanical coefficients is crucial for capturing the RTV volume change during ablation and enable a predictive capability for heatshield tile interface response. A plan for future calibration of thermomechanical properties evolution during degradation is discussed, as a key next step to close the new model.

RTV↗