SEARCH · Engineering Papers
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.
NASA Glenn Electrified Aircraft Propulsion (EAP) Presentation Outline: High-Voltage, High-Power EAP - Standards and Regulations - Safety and Certification - High Power Density, High-Efficiency Electric Machines and Drives - Thermal Management - Fault Management - Superconducting Electric Machines - Materials Research - Modeling Tools - Advanced Controls
Report describes fabrication of graphite-fiber/copper-matrix composite plates, measurements of thermal conductivities of plates at temperatures from ambient to 1,073 K, and measurements of thermal expansions of plates from ambient temperature to 1,050 K. Composites promising lightweight, high-thermal-conductivity materials proposed for use in heat exchangers and other heat-transfer components of power systems in spacecraft and hypersonic aircraft. Graphite/copper also of interest as model composite material.
The Near Earth Asteroid (NEA) Scout is a solar sail mission whose objective is to scout at least one Near Earth Asteroid to help prepare for human missions to Near Earth Asteroids. NEA Scout will launch as a secondary payload on the first SLS-Orion mission. NEA Scout will perform a small trim maneuver shortly after deploy from the spent SLS upper stage using a cold gas propulsion system, but from that point on will depend entirely on the solar sail for thrust. As such, it is important to accurately characterize the thrust of the sail in order to achieve mission success. Additionally, the solar sail creates a relatively large solar disturbance torque that must be mitigated. For early mission design studies a flat plate model of the solar sail with a fixed center of pressure was adequate, but as mission concepts and the sail design matured, greater fidelity was required. Here we discuss the progress to a three-dimensional sail model that includes the effects of tension and thermal deformation that has been derived from a large structural Finite Element Model (FEM) developed by the Langley Research Center. We have found that the deformed sail membrane affects torque relatively much more than thrust; a flat plate model could potentially model thrust well enough to close mission design studies, but a three-dimensional solar sail is essential to control system design. The three-dimensional solar sail model revealed that thermal deformations of unshielded booms would create unacceptably large solar disturbance torques. The original large FEM model was used in control and mission simulations, but was resulted in simulations with prohibitive run times. This led us to adapt the Generalized Sail Model (GSM) of Rios-Reyes. A design reference sail model has been baselined for NEA Scout and has been used to design the mission and control system for the sailcraft. Additionally, since NEA Scout uses reaction wheels for attitude pointing and control, the solar torque model is essentially to successfully design the NEA Scout momentum management control system. We have also updated the estimate of diffusivity used for the aluminized sail material based on optical testing of wrinkled sail material. The model presented here represents the current state of the art of NASA's ability to model solar sail thrust and torque.
Micromechanics-based modeling of composite material behaviors requires an accurate assessment of the constituent properties and behaviors. For the specific case of continuous-fiber-reinforced metal matrix composites (MMC's) manufactured from a foil/fiber/foil process, much emphasis has been placed on characterizing foil-based matrix materials that have been fabricated in the same way as the composite. Such materials are believed to yield mechanical properties and behaviors that are representative of the matrix constituent within the composite (in situ matrix). Therefore, these materials are desired for micromechanics modeling input. Unfortunately, such foils are extremely expensive to fabricate and procure because of the labor-intensive rolling process needed to produce them. As a potential solution to this problem that would maintain appropriately representative in situ properties, the matrix constituent could be characterized with sheet-based materials, which are considerably less expensive to manufacture than foils, are more readily procured, and result in fewer plies to obtain a desired panel thickness. The critical question is, however, does the consolidated sheet material exhibit the same properties and behaviors as do the consolidated foils? Researchers at NASA Lewis Research Center's Life Prediction Branch completed a detailed experimental investigation to answer this question for three titanium alloys commonly used in metal matrix composite form.
The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heatshield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material [1]. PICA is a lightweight carbon fiber/polymeric resin material that offers outstanding performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP) [2]. The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on OpenFOAM (PATO) [4,5,6] using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA [7], is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.
A model was developed to generate quantitative estimates of the risk associated with the release of graphite fibers during fires involving commercial aircraft constructed with graphite fiber composite materials. The model was used to estimate the risk associated with accidents at several U.S. airports. These results were then combined to provide an estimate of the total risk to the nation.
The space environment is detrimental to the exterior of the spacecraft and imposes stringent requirements on spacecraft materials. Space materials are often expected to perform for 15-20 years in the harsh space environment without significant degradation of their optical, electrical, mechanical and thermal properties. Further, as commercial activities increase in the space domain, development of predictive capability for material evolution is critical for agile adaption of novel materials. Therefore, it is very important to understand the effects of the space environments on the spacecraft materials. The Materials on the International Space Station Experiment Flight Facility (MISSE-FF) has flown numerous material samples to investigate the effect of LEO space weather exposure on the performance and durability of materials and devices. We are manifest on the MISSE-16 mission, scheduled to launch in June 2022, to fly fifteen novel and well-characterized materials in the LEO environment for a duration of six months. Changes in spectral reflectivity will be measured in real time throughout the mission with an RGB/IR camera. These time-resolved data will function as a "space-truth" reference for our team's ongoing laboratory-based space weather-material interaction experiments. Correlation of the MISSE-16 data with extensive ground testing of duplicate samples under simulated space weather conditions will enable development of fundamental chemical models for material degradation. This paper discusses preliminary results from the ground test campaign to collect the RGB/IR images for pristine and damaged materials and the development of machine learning algorithms to extract reflectance spectra from the color images.
The space environment is detrimental to the exterior of the spacecraft and imposes stringent requirements on spacecraft materials. Space materials are often expected to perform for 15-20 years in the harsh space environment without significant degradation of their optical, electrical, mechanical and thermal properties. Further, as commercial activities increase in the space domain, development of predictive capability for material evolution is critical for agile adaption of novel materials. Therefore, it is very important to understand the effects of the space environments on the spacecraft materials. The Materials on the International Space Station Experiment Flight Facility (MISSE-FF) has flown numerous material samples to investigate the effect of LEO space weather exposure on the performance and durability of materials and devices. We are manifest on the MISSE-16 mission, scheduled to launch in June 2022, to fly fifteen novel and well-characterized materials in the LEO environment for a duration of six months. Changes in spectral reflectivity will be measured in real time throughout the mission with an RGB/IR camera. These time-resolved data will function as a "space-truth" reference for our team's ongoing laboratory-based space weather-material interaction experiments. Correlation of the MISSE-16 data with extensive ground testing of duplicate samples under simulated space weather conditions will enable development of fundamental chemical models for material degradation. This paper discusses preliminary results from the ground test campaign to collect the RGB/IR images for pristine and damaged materials and the development of machine learning algorithms to extract reflectance spectra from the color images.
RSRM nozzle liner components have been analyzed and tested to explore the occurrence of anomalous material performance known as pocketing erosion. Primary physical factors that contribute to pocketing seem to include the geometric permeability, which governs pore pressure magnitudes and hence load, and carbon fiber high temperature tensile strength, which defines a material limiting capability. The study reports on the results of a coupled thermostructural finite element analysis of Carbon Cloth Phenolic (CCP) material tested at the Laser Hardened Material Evaluation Laboratory (the LHMEL facility). Modeled test configurations will be limited to the special case of where temperature gradients are oriented perpendicular to the composite material ply angle. Analyses were conducted using a transient, one-dimensional flow/thermal finite element code that models pore pressure and temperature distributions and in an explicitly coupled formulation, passes this information to a 2-dimensional finite element structural model for determination of the stress/deformation behavior of the orthotropic fiber/matrix CCP. Pore pressures are generated by thermal decomposition of the phenolic resin which evolve as a multi-component gas phase which is partially trapped in the porous microstructure of the composite. The nature of resultant pressures are described by using the Darcy relationships which have been modified to permit a multi-specie mass and momentum balance including water vapor condensation. Solution to the conjugate flow/thermal equations were performed using the SINDA code. Of particular importance to this problem was the implementation of a char and deformation state dependent (geometric) permeability as describing a first order interaction between the flow/thermal and structural models. Material property models are used to characterize the solid phase mechanical stiffness and failure. Structural calculations were performed using the ABAQUS code. Iterations were made between the two codes involving the dependent variables temperature, pressure and across-ply strain level. Model results comparisons are made for three different surface heat rates and dependent variable sensitivities discussed for the various cases.
Coupled fluid-material response analyses of arc-jet wedge ablation tests conducted in a NASA Ames arc-jet facility are considered. These tests were conducted using blunt wedge models placed in a free jet downstream of the 6-inch diameter conical nozzle in the Ames 60-MW Interaction Heating Facility. The fluid analysis includes computational Navier-Stokes simulations of the nonequilibrium flowfield in the facility nozzle and test box as well as the flowfield over the models. The material response analysis includes simulation of two-dimensional surface ablation and internal heat conduction, thermal decomposition, and pyrolysis gas flow. For ablating test articles undergoing shape change, the material response and fluid analyses are coupled in order to calculate the time dependent surface heating and pressure distributions that result from shape change. The ablating material used in these arc-jet tests was Phenolic Impregnated Carbon Ablator. Effects of the test article shape change on fluid and material response simulations are demonstrated, and computational predictions of surface recession, shape change, and in-depth temperatures are compared with the experimental measurements.
In computing the response of an ablating thermal protection system during atmospheric entry, the aerothermal environment and material response are generally computed separately with a blowing correction term in the material response model to account for the blowing of char and pyrolysis gases. In this work, we apply a coupled approach in which pyrolysis blowing gases, computed in the PATO material response code, are used with a blowing boundary condition in the DPLR hypersonic CFD code. This leads to an iterative method in which blowing products from PATO are input into DPLR to update surface heating estimates. The full iterative method, with the addition of radiative heating estimates using the NEQAIR radiation solver, is shown in Fig. 1. The method is demonstrated on a sphere case with the environment and material properties based on the Mars Science Laboratory entry. Future work includes utilizing this method in computing full 3D material response during the Mars 2020 entry.
The Mars Science Laboratory (MSL) was protected during its atmospheric entry by an instrumented heat-shield using NASA's Phenolic Impregnated Carbon Ablator (PICA) material. PICA is a lightweight carbon fiber/polymeric resin material that offers out-standing performances for protecting probes during planetary entry. The Mars Entry Descent and Landing Instrument (MEDLI) suite on MSL offers unique in-flight validation data for models of material response and atmospheric entry. MEDLI recorded, among other things, time-resolved in-depth temperature data of PICA using thermocouple sensors assembled in the MEDLI Integrated Sensor Plugs (MISP). The objective of this work is to showcase and analyze the coupling between the material response and the aerothermal environment. As shown in Figure 1, the workflow is divided into the following steps. First, the aerothermal properties are computed in the Data Parallel Line Relaxation (DPLR) code [3] and used with the Nonequilibrium air radiation (NEQAIR) program [8] to compute radiative heating. Second, the thermal response inside the material is computed in the Porous material Analysis Toolbox based on Open-FOAM (PATO) using a fixed blowing correction parameter. Third, the pyrolysis gases computed in PATO are used as inputs to a blowing boundary condition within DPLR. Fourth, the new environment properties from DPLR are used in NEQAIR to provide an updated solution, then both the updated aerothermal environment and radiative heating are used in PATO without blowing correction. The third and fourth steps are then repeated until convergence in surface temperature is obtained. Convergence in the radiative heating is generally achieved before surface temperature, at which point the radiative heating is no longer updated. Char mass loss rates are forced to zero to produce a non-receding surface condition. For early time points in the trajectory, where flow around the MSL aeroshell is rarefied, the Direct Simulation Monte Carlo (DSMC) code, SPARTA, is used to compute the aerothermal environment. Iteration between PATO and SPARTA is not performed due to the computational cost of DSMC simulations. Preliminary results of the coupling between PATO and DPLR for the MSL heatshield atmospheric entry model are presented in Figures 2-4 at 65 seconds after entry interface. Figure 2 shows the surface temperature results from an uncoupled simulation in PATO with the blowing correction parameter applied (left) along with the coupled surface temperature after iteration (right). Figure 3 shows the surface temperature along the centerline from windward to leeward for easier comparison. Figure 4 shows the coupled and uncoupled pyrolysis gas blowing rate. Mars 2020 used a similar heatshield consisting of PICA for thermal protection during entry, descent, and landing. In preparation for Mars 2020 post-flight analysis, the predictive material response capability is benchmarked against flight data from MEDLI. This work represents an important milestone toward the development of validated predictive capabilities for designing thermal protection systems for planetary probes.
A summary of the status of this five-year project which is now in its third year of research and development is presented. The goal of the project is the development of several methodologies for probabilistic structural modeling. Probabilistic structural modeling consists of stochastic models of material properties, part geometries, boundary conditions, as well as loading conditions. The current presentation focuses on one methodology - coupling of an advanced finite element structural analysis code with probabilistic modeling strategies. The essential algorithm developments for combining the finite element and probabilistic analysis methods are reported. The validity of the resulting probabilistic structural analysis method is confirmed through a series of test problems with exact results based on Monte Carlo simulations. Additionally, the applicability of the method to a Space Propulsion System (a turbine blade) is demonstrated for static stresses.
A new model is proposed that represents the kinematics of kink-band formation and propagation within the framework of a mesoscale continuum damage mechanics (CDM) model. The model uses the recently proposed deformation gradient decomposition approach to represent a kink band as a displacement jump via a cohesive interface that is embedded in an elastic bulk material. The model is capable of representing the combination of matrix failure in the frame of a misaligned fiber and instability due to shear nonlinearity. In contrast to conventional linear or bilinear strain softening laws used in most mesoscale CDM models for longitudinal compression, the constitutive response of the proposed model includes features predicted by detailed micromechanical models. These features include: 1) the rotational kinematics of the kink band, 2) an instability when the peak load is reached, and 3) a nonzero plateau stress under large strains.
In computing the response of an ablating thermal protection system during atmospheric entry, the aerothermal environment and material response are generally computed separately with a blowing correction term in the material response model to account for the blowing of char and pyrolysis gases [1]. In this work, we apply a coupled approach in which pyrolysis blowing gases, computed in the PATO material response code [2], are used with a blowing boundary condition in the DPLR hypersonic CFD code [3]. This leads to an iterative method in which blowing products from PATO are input into DPLR to update surface heating estimates. The full iterative method, with the addition of radiative heating estimates using the NEQAIR radiation solver [4], is shown in Fig. 1. The method is demonstrated on a sphere case with the environment and material properties based on the Mars Science Laboratory entry. Future work includes utilizing this method in computing full 3D material response during the Mars 2020 entry.
For thermal protection system (heat shield) design for space vehicle entry into earth and other planetary atmospheres, it is essential to know the augmentation of the heat flux due to vehicle surface roughness. At the NASA Ames Hypervelocity Free Flight Aerodynamic Facility (HFFAF) ballistic range, a campaign of heat flux studies on rough models, using infrared camera techniques, has been initiated. Several phenomena can interfere with obtaining good heat flux data when using this measuring technique. These include leakage of the hot drive gas in the gun barrel through joints in the sabot (model carrier) to create spurious thermal imprints on the model forebody, deposition of sabot material on the model forebody, thereby changing the thermal properties of the model surface and unknown in-barrel heating of the model. This report presents developments in launch techniques to greatly reduce or eliminate these problems. The techniques include the use of obturator cups behind the launch package, enclosed versus open front sabot designs and the use of hydrogen gas in the launch tube. Attention also had to be paid to the problem of the obturator drafting behind the model and impacting the model. Of the techniques presented, the obturator cups and hydrogen in the launch tube were successful when properly implemented
Ablative thermal protection system (TPS) materials for spacecraft are composites that are often made out of carbon-based reinforcement and a polymeric matrix. They endure high-temperature oxidation and surface recession when re-entering Earth’s atmosphere. Ablation is the result of many coupled and competing thermal, mechanical, and chemical phenomena, and it is difficult to isolate the role of each on the overall degradation of the TPS. Here we develop an ablation model for material recession coupled explicitly to finite rate carbon oxidation in complex microstructures. In this work, Stochastic PArallel Rarified-gas Time-accurate Analyzer (SPARTA), a direct-simulation Monte Carlo (DSMC) code, is modified to allow oxidation-driven ablation of implicitly defined carbon surfaces. In SPARTA, implicit surfaces are generated from the grid corner point values via a marching cubes algorithm, therefore creating a new set of surface elements every time ablation is performed. The finite-rate oxidation model developed by Gopalan et. al, was adapted to tally surface reactions and other surface data on a per-grid cell basis. The ablation functionality was also adjusted so once the reactions have occurred, the number of reactions leading to CO formation can be converted to corner point reduction values; therefore, carbon removal is directly proportional to surface recession. We also develop robust algorithms which handle the evolution of the flow cells and solid material regions, including split cells (flow cell divided in two by a solid surface). Finally, we demonstrate our implicit chemistry model for 2D and 3D geometries by producing reaction statistics and detailed visualization of oxidation-induced material recession at the microscale.
In this work, a neural network model was developed to predict the constant pressure heat capacity for materials in the rare-earth oxide—silica material space. Several model architectures were trained and tested on heat capacity data generated from first-principles density functional theory calculations. Hyperparameter optimization was performed, and the optimal model was selected for heat capacity predictions. The optimal model architecture was found to have a root-mean-squared error of 5.12 ± 3.37 J/mol-K. The optimal model architecture was then used in a bagging ensemble model trained using the leave-one-group-out method to provide error estimates for model predictions. The out-of-bag score for the ensemble model was 0.997. The predicted heat capacities agree well with the DFT and experimental results and were computed orders of magnitude faster than DFT simulations. Machine learning shows the potential to provide a suitable surrogate model for thermochemical property predictions for candidate environmental barrier coating materials but refining of input material features and model architectures could further improve accuracy for these models.