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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Thermochemical Ablation Analysis of the Orion Heatshield

The Orion Multi-Purpose Crew Vehicle will one day carry astronauts to the Moon and beyond, and Orion's heatshield is a critical component in ensuring their safe return to Earth. The Orion heatshield is the structural component responsible for absorbing the intense heating environment caused by re-entry to Earth's atmosphere. The heatshield is primarily composed of Avcoat, an ablative material that is consumed during the re-entry process. Ablation is primarily characterized by two processes: pyrolysis and recession. The decomposition of in-depth virgin material is known as pyrolysis. Recession occurs when the exposed surface of the heatshield reacts with the surrounding flow. The Orion heatshield design was changed from an individually filled Avcoat honeycomb to a molded block Avcoat design. The molded block Avcoat heatshield relies on an adhesive bond to keep it attached to the capsule. In some locations on the heatshield, the integrity of the adhesive bond cannot be verified. For these locations, a mechanical retention device was proposed. Avcoat ablation was modelled in CHAR and the in-depth virgin material temperatures were used in a Thermal Desktop model of the mechanical retention device. The retention device was analyzed and shown to cause a large increase in the maximum bondline temperature. In order to study the impact of individual ablation modelling parameters on the heatshield sizing process, a Monte Carlo simulation of the sizing process was proposed. The simulation will give the sensitivity of the ablation model to each of its input parameters. As part of the Monte Carlo simulation, statistical uncertainties on material properties were required for Avcoat. Several properties were difficult to acquire uncertainties for: the pyrolysis gas enthalpy, non-dimensional mass loss rate (B´c), and Arrhenius equation parameters. Variability in the elemental composition of Avcoat was used as the basis for determining the statistical uncertainty in pyrolysis gas enthalpy and B´c. A MATLAB program was developed to allow for faster, more accurate and automated computation of Arrhenius reaction parameters. These parameters are required for a material model to be used in the CHAR ablation analysis program. This MATLAB program, along with thermogravimetric analysis (TGA) data, was used to generate uncertainties on the Arrhenius parameters for Avcoat. In addition, the TGA fitting program was developed to provide Arrhenius parameters for the ablation model of the gap filler material, RTV silicone.

Sixel, William↗
NASA NTRS1993

The importance of properties in modeling

Casting and welding of superalloys, stainless steel and titanium alloys are processes which can be improved through modeling of heat flow, fluid flow, residual stress development, and microstructural evolution. These simulations require inputs of thermophysical data, some of which involves the partially or totally liquid state. In particular, these processes involve melting, flow in the liquid, and solidification. Modeling of such processes can lead to an improved understanding of defects such as shrinkage, inclusions, cracks, incomplete filling (or penetration), macrosegregation, improper grain structure, and deviations from dimensional specifications. Effective modeling can shorten process development time and improve quality. An approach to these problems is to develop efficient models; validate through correlations with thermal, distortion, and microstructural data; run parametric studies; extract knowledge based rules; and apply to adaptive closed loop control systems. With the appropriate pre- and post-processing, such analyses can be made 'user friendly'. This would include graphical user interfaces as well as realistic images and color maps. In such form, these models can be used for sensitivity analyses, which are useful in defining appropriate sensors and in the development of control strategies. Such modeling can be done at several levels, e.g., the MARO level, modeling large scale phenomena such as heat and fluid flow or material deformation; the MICRO level, modeling the development of dendrites, grains or precipitates; or at the NANO level, modeling point defects, dislocations, stacking faults, etc. There are many computational issues associated with these simulations, e.g. computational efficiency and accuracy. In addition, there are many materials issues, not the least of which is the availability of accurate high temperature thermophysical data for complex alloys. This would include latent heat of fusion, temperature dependent heat capacity and thermal conductivity (for liquid and solid), viscosity, surface tension, thermal expansion, mechanical properties, etc. Preliminary data is frequently gathered from the literature; however, this is often not available for modern alloys. If additional data are required, measurements can be used; however, these are costly, time consuming and can be erroneous due to a lack of testing standards or impure materials. Microstructural predictors can be extracted from thermal information, e.g. cooling rate and thermal gradient; the prediction of microstructure is dependent on solidus and liquidus temperature, mushy zone permeability, the solidification curve, volume changes, phase transformations, alloying effects (such as surface tension or viscosity), mold/metal reactions, metal/environment reactions, etc. Defect maps may be needed to predict the onset of shrinkage, hot cracking or 'freckling'. Constants may be needed for stress relaxation, dendrite coarsening, vaporization, etc. Visualization was used as a tool to better comprehend complex data sets associated with the analysis of directional solidification (including crystal growth) and welding. Examples include not only isotherms, but also cooling rate, growth rate and thermal gradient. The latter two are not single valued scalars, but rather time and space dependent vector fields. Efficient models were developed for both casting and welding to predict heat flow and the relationship to dendrite and grain growth. These codes include many of the non-linear effects, e.g. radiation, which dominate these processes. The home-built FDM code(s) were designed to be useful not only to the scientist, but also to the process engineer. Special output can be requested to compare directly to experimental data. Visualization procedures were developed to visualize critical results, e.g. fusion zone width at the surface opposite that where the arc is applied ('penetration'). Both elaborate and simplified distortion analyses were carried out. It is clear that extensive mechanical property data are critical in order to accurately predict residual stress patterns. A scheme is currently being developed to integrate these modeling tools into a set of control algorithms; however, the success of this approach is critically dependent on the availability of accurate high temperature thermophysical data.

Giamei, A. F.↗
NASA NTRS1981

Survey of the supporting research and technology for the thermal protection of the Galileo Probe

The Galileo Probe, which is scheduled to be launched in 1985 and to enter the hydrogen-helium atmosphere of Jupiter up to 1,475 days later, presents thermal protection problems that are far more difficult than those experienced in previous planetary entry missions. The high entry speed of the Probe will cause forebody heating rates orders of magnitude greater than those encountered in the Apollo and Pioneer Venus missions, severe afterbody heating from base-flow radiation, and thermochemical ablation rates for carbon phenolic that rival the free-stream mass flux. This paper presents a comprehensive survey of the experimental work and computational research that provide technological support for the Probe's heat-shield design effort. The survey includes atmospheric modeling; both approximate and first-principle computations of flow fields and heat-shield material response; base heating; turbulence modelling; new computational techniques; experimental heating and materials studies; code validation efforts; and a set of 'consensus' first-principle flow-field solutions through the entry maneuver, with predictions of the corresponding thermal protection requirements.

Howe, J. T.↗
NASA NTRS2012

Simulating the Response of a Composite Honeycomb Energy Absorber: Full-Scale Impact Testing - Part 2

NASA has sponsored research to evaluate an externally deployable composite honeycomb designed to attenuate loads in the event of a helicopter crash. The concept, designated the Deployable Energy Absorber (DEA), is an expandable Kevlar(Registered TradeMark) honeycomb. The DEA has a flexible hinge that allows the honeycomb to be stowed collapsed until needed during an emergency. Evaluation of the DEA began with material characterization of the Kevlar(Registered TradeMark)-129 fabric/epoxy, and ended with a full-scale crash test of a retrofitted MD-500 helicopter. During each evaluation phase, finite element models of the test articles were developed and simulations were performed using the dynamic finite element code, LS-DYNA(Registered TradeMark). The paper will focus on simulations of two full-scale impact tests involving the DEA, a mass-simulator and a full-scale crash of an instrumented MD-500 helicopter. Isotropic (MAT24) and composite (MAT58) material models, which were assigned to DEA shell elements, were compared. Based on simulations results, the MAT58 model showed better agreement with test.

Fasanella, Edwin L.↗
NASA NTRS1993

Models of classical and recurrent novae

The behavior of novae may be divided roughly into two separate stages: quiescence and outburst. However, at closer inspection, both stages cannot be separated. It should be attempted to explain features in both stages with a similar model. Various simple models to explain the observed light and spectral observations during post optical maximum activity are conceivable. In instantaneous ejection models, all or nearly all material is ejected in a time that is short compared with the duration of post optical maximum activity. Instantaneous ejection type 1 models are those where the ejected material is in a fairly thin shell, the thickness of which remains small. In the instantaneous ejection type 2 model ('Hubble Flow'), a thick envelope is ejected instantaneously. This envelope remains thick as different parts have different velocities. Continued ejection models emphasize the importance of winds from the nova after optical maximum. Ejection is supposed to occur from one of the components of the central binary, and one can imagine a general swelling of one of the components, so that something resembling a normal, almost stationary, stellar photosphere is observed after optical maximum. The observed characteristics of recurrent novae in general are rather different from those of classical novae, thus, models for these stars need not be the same.

Friedjung, Michael↗
NASA NTRS1982

The composition of ultra-heavy cosmic rays

Preliminary results on the composition of ultraheavy nuclei in the cosmic rays are now available from experiments on HEAO-3 and Ariel-VI satellites. These results are compared with the composition predicted for reasonable propagation models using source abundances similar to those of solar system material. The effects of varying model parameters and of preferential acceleration are discussed. Discrepancies between predicted and observed abundances are larger than can be accounted for solely by first ionization potential effects and are interpreted in terms of a nucleosynthesis process.

Protheroe, R. J.↗
NASA NTRS1987

Stability and rheology of dispersions of silicon nitride and silicon carbide

The relationship between the surface and colloid chemistry of commercial ultra-fine silicon carbide and silicon nitride powders was examined by a variety of standard characterization techniques and by methodologies especially developed for ceramic dispersions. These include electrokinetic measurement, surface titration, and surface spectroscopies. The effects of powder pretreatment and modification strategies, which can be utilized to augment control of processing characteristics, were monitored with these technologies. Both silicon carbide and nitride were found to exhibit silica-like surface chemistries, but silicon nitride powders possess an additional amine surface functionality. Colloidal characteristics of the various nitride powders in aqueous suspension is believed to be highly dependent on the relative amounts of the two types of surface groups, which in turn is determined by the powder synthesis route. The differences in the apparent colloidal characteristics for silicon nitride powders cannot be attributed to the specific absorption of ammonium ions. Development of a model for the prediction of double-layer characteristics of materials with a hybrid site interface facilitated understanding and prediction of the behavior of both surface charge and surface potential for these materials. The utility of the model in application to silicon nitride powders was demonstrated.

Feke, Donald L.↗
NASA NTRS1988

A constitutive model with damage for high temperature superalloys

A unified constitutive model is searched for that is applicable for high temperature superalloys used in modern gas turbines. Two unified inelastic state variable constitutive models were evaluated for use with the damage parameter proposed by Kachanov. The first is a model (Bodner, Partom) in which hardening is modeled through the use of a single state variable that is similar to drag stress. The other (Ramaswamy) employs both a drag stress and back stress. The extension was successful for predicting the tensile, creep, fatigue, torsional and nonproportional response of Rene' 80 at several temperatures. In both formulations, a cumulative damage parameter is introduced to model the changes in material properties due to the formation of microcracks and microvoids that ultimately produce a macroscopic crack. A back stress/drag stress/damage model was evaluated for Rene' 95 at 1200 F and is shown to predict the tensile, creep, and cyclic loading responses reasonably well.

Sherwood, J. A.↗
NASA NTRS2005

Lifing of Engine Components

The successful development of advanced aerospace engines depends greatly on the capabilities of high performance materials and structures. Advanced materials, such as nickel based single crystal alloys, metal foam, advanced copper alloys, and ceramics matrix composites, have been engineered to provide higher engine temperature and stress capabilities. Thermal barrier coatings have been developed to improve component durability and fuel efficiency, by reducing the substrate hot wall metal temperature and protecting against oxidation and blanching. However, these coatings are prone to oxidation and delamination failures. In order to implement the use of these materials in advanced engines, it is necessary to understand and model the evolution of damage of the metal substrate as well as the coating under actual engine conditions. The models and the understanding of material behavior are utilized in the development of a life prediction methodology for hot section components. The research activities were focused on determining the stress and strain fields in an engine environment under combined thermo-mechanical loads to develop life prediction methodologies consistent with the observed damage formation of the coating and the substrates.

Source record↗
NASA NTRS2005

Dynamics Impact Tolerance of Shuttle RCC Leading Edge Panels Using LS-DYNA

This paper describes a research program conducted to enable accurate prediction of the impact tolerance of the shuttle Orbiter leading-edge wing panels using physics-based codes such as LS-DYNA, a nonlinear, explicit transient dynamic finite element code. The shuttle leading-edge panels are constructed of Reinforced-Carbon-Carbon (RCC) composite material, which is used because of its thermal properties to protect the shuttle during reentry into the Earth's atmosphere. Accurate predictions of impact damage from insulating foam and other debris strikes that occur during launch required materials characterization of expected debris, including strain-rate effects. First, analytical models of individual foam and RCC materials were validated. Next, analytical models of foam cylinders impacting 6- in. x 6-in. RCC flat plates were developed and validated. LS-DYNA pre-test models of the RCC flat plate specimens established the impact velocity of the test for three damage levels: no-detectable damage, non-destructive evaluation (NDE) detectable damage, or visible damage such as a through crack or hole. Finally, the threshold of impact damage for RCC on representative Orbiter wing panels was predicted for both a small through crack and for NDE-detectable damage.

Fasanella, Edwin L.↗
NASA NTRS2005

Cost Modeling for low-cost planetary missions

This presentation will provide an overview of the JPL parametric cost models used to estimate flight science spacecrafts and instruments. This material will emphasize the cost model approaches to estimate low-cost flight hardware, sensors, and instrumentation, and to perform cost-risk assessments. This presentation will also discuss JPL approaches to perform cost modeling and the methodologies and analyses used to capture low-cost vs. key cost drivers.

cost model↗
NASA NTRS2008

Dynamic Impact Tolerance of Shuttle RCC Leading Edge Panels using LS-DYNA

This paper describes a research program conducted to enable accurate prediction of the impact tolerance of the shuttle Orbiter leading-edge wing panels using 'physics-based- codes such as LS-DYNA, a nonlinear, explicit transient dynamic finite element code. The shuttle leading-edge panels are constructed of Reinforced-Carbon-Carbon (RCC) composite material, which issued because of its thermal properties to protect the shuttle during re-entry into the Earth's atmosphere. Accurate predictions of impact damage from insulating foam and other debris strikes that occur during launch required materials characterization of expected debris, including strain-rate effects. First, analytical models of individual foam and RCC materials were validated. Next, analytical models of individual foam cylinders impacting 6-in. x 6-in. RCC flat plates were developed and validated. LS-DYNA pre-test models of the RCC flat plate specimens established the impact velocity of the test for three damage levels: no-detectable damage, non-destructive evaluation (NDE) detectable damage, or visible damage such as a through crack or hole. Finally, the threshold of impact damage for RCC on representative Orbiter wing panels was predicted for both a small through crack and for NDE-detectable damage.

Fasanella, Edwin↗
NASA NTRSDate not supplied

Assessment of Pioneer Venus Entry Heating with Coupled Radiation and Ablation

The objective of this work is to reconstruct thermocouple temperature measurements taken at the stagnation point of the Pioneer Venus large probe using uncoupled and coupled ablation calculations. The time history of the large probe forebody entry heating was obtained from a trajectory-based analysis that employed a material response solver to model equilibrium and nonequilibrium ablation chemistry at the surface. In the uncoupled approach, equilibrium ablation chemistry was assumed in the material response calculations, but the material response outputs were not returned to the flowfield solver. Therefore, the effects of char and pyrolysis gas injection into the boundary layer due to ablation (e.g., convective blockage and radiative energy absorption) did not influence the flowfield solutions. In the coupled approach, nonequilibrium ablation chemistry was performed and the material response outputs were returned to the flowfield solver. Therefore, the effects of char and pyrolysis gas injection into the boundary layer were accounted for in the flowfield solutions. Both uncoupled and coupled approaches produced reconstructed temperatures that were within 10% of the flight measurements through the entry heat pulse, with the coupled results comparing more favorably with the measurements. Overall, the reconstructed thermocouple temperatures in this work significantly improved on the results of previous reconstructions for the large probe.

Jannuel V V Cabrera↗
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 NTRS2016

Development of a Subcell Based Modeling Approach for Modeling the Architecturally Dependent Impact Response of Triaxially Braided Polymer Matrix Composites

Understanding the high velocity impact response of polymer matrix composites with complex architectures is critical to many aerospace applications, including engine fan blade containment systems where the structure must be able to completely contain fan blades in the event of a blade-out. Despite the benefits offered by these materials, the complex nature of textile composites presents a significant challenge for the prediction of deformation and damage under both quasi-static and impact loading conditions. The relatively large mesoscale repeating unit cell (in comparison to the size of structural components) causes the material to behave like a structure rather than a homogeneous material. Impact experiments conducted at NASA Glenn Research Center have shown the damage patterns to be a function of the underlying material architecture. Traditional computational techniques that involve modeling these materials using smeared homogeneous, orthotropic material properties at the macroscale result in simulated damage patterns that are a function of the structural geometry, but not the material architecture. In order to preserve heterogeneity at the highest length scale in a robust yet computationally efficient manner, and capture the architecturally dependent damage patterns, a previously-developed subcell modeling approach where the braided composite unit cell is approximated as a series of four adjacent laminated composites is utilized. This work discusses the implementation of the subcell methodology into the commercial transient dynamic finite element code LS-DYNA (Livermore Software Technology Corp.). Verification and validation studies are also presented, including simulation of the tensile response of straight-sided and notched quasi-static coupons composed of a T700/PR520 triaxially braided [0deg/60deg/-60deg] composite. Based on the results of the verification and validation studies, advantages and limitations of the methodology as well as plans for future work are discussed.

finite element method↗
NASA NTRS1999

Simulating the Formation of Lunar Floor-Fracture Craters Using Elastoviscoplastic Relaxation

Lunar floor-fracture craters formed during the height of mare basalt emplacement. Due to a general temporal and spatial relation with the maria, these craters, numbering some 200, may be diagnostic of the thermal structure of the crust during this time. As the name suggests, these craters exhibit brittle failure, generally limited to the central floor region. That, and a shallower depth than fresh lunar craters, has led to two main theories as to their formation: laccolith emplacement under the crater and viscous relaxation. The implications of each model for the state of the Moon's crust during this time are quite different, so the viability of each model must be checked. Laccolith emplacement has been treated elsewhere. However, previous attempts to study the relaxation of the craters have assumed only a uniform, Newtonian viscous response of the near surface to the topographic driving forces, and simply postulated that the fractures resulted from tensile stresses associated with floor uplift. Here, we use a more sophisticated rheological model that includes not only non-Newtonian viscous behavior (i.e., the viscosity is stress-dependent), but also incorporates elastic behavior and a plastic component to the rheology to directly simulate the formation of the floor fractures. The results of our simulations show that while elastoviscoplastic relaxation is potentially viable for larger floor-fracture craters, it is not viable for craters with diameters < or = 60 km, the size of the majority of floor-fracture craters. We employ the finite element method, a numerical technique well suited for boundary-value problems, via the commercially available MARC software package. To test the viability of topographic relaxation, our goal is to prepare the simulations as to maximize the amount of relaxation. We take advantage of the natural axisymmetry of craters, simulating one radial plane. Initial shapes are based on data for fresh craters from Pike. To simplify implementation, a fourth order polynomial is used for the basin, while a third order inverse function is used for the rim. This form closely approximates the long-wavelength behavior of complex craters, while ignoring higher-frequency topography, save the rim. This approximation is appropriate because crater relaxation is strongly controlled by long-wavelength topography. Loading is accomplished assuming a uniform gravity field (1.62 m/s-square) and a uniform density of 2900 kg/cubic m. The initial stress state is set to be hydrostatic, with an additional pressure term to account for any overlying topography. This additional pressure term is tapered exponentially with depth. While the simulations quickly settle on a preferred stress state, and while the final solution is fairly insensitive to the choice of the e-folding depth of the taper, selecting an e-folding depth close to the diameter of the crater sets the initial stress state near the preferred state. We assume a diuranally averaged surface temperature of -20 C, and allow temperature to increase with depth at a rate of 50 K//km. Assuming a thermal conductivity of 2 W/in/K, this gradient translates to a heat flow of 100mW/square m, an extremely high value for the Moon. Temperature, of course, will not increase without bound. To maximize relaxation, we allow our temperature profile to increase linearly until it reaches the solidus (assumed to be 1200C) at a depth of 24.4 km, at which point it is kept constant. The presence of melt will drop the bulk viscosity; however, we have no rheological control for partial melts. Therefore, we make no attempt to simulate this situation. Elastoviscoplastic rheological model. In general, geologic materials can behave in three main ways: elastically, viscously (via solid-state creep), and brittly (plasticity is a continuum approach to simulate this phenomenon). We combine these three deformation mechanisms in an extended Maxwell solid, where the total strain can be broken down into a simple summation of the elastic, creep, and plastic strains. In relaxation phenomena in general, the system takes advantage of any means possible to eliminate deviatoric stresses by relaxing away the topography. Previous analyses have only modeled the viscous response. Comparatively, the elastic response in our model can augment the relaxation, to a point. This effect decreases as the elastic response becomes stiffer; indeed, in the limit of infinite elastic Young's modulus (and with no plasticity), the solution converges on the purely viscous solution. Igneous rocks common to the lunar near-surface have Young's modulii in the range of 10-100 GPa. To maximize relaxation, we use a Young's modulus of 10 GPa. (There is negligible sensitivity to the other elastic modulus, the Poisson's ratio; we use 0.25.) For the viscous response, we use a flow law for steady-state creep in thoroughly dried Columbia diabase, because the high plagioclase (about 70 vol%) and orthopyroxene (about 17 vol%) content is similar to the composition of the lunar highland crust as described by remote sensing and sample studies: noritic anorthosite. This flow law is highly non-Newtonian, i.e., the viscosity is highly stress dependent. That, and the variability with temperature, stands in strong contrast to previous examinations of lunar floor-fracture crater relaxation. To model discrete, brittle faulting, we assume "Byerlee's rule," a standard geodynamical technique. We implement this "rule" with an-angle of internal friction of about 40 deg, and a higher-than-normal cohesion of about 3.2 MPa (to approximate the breaking of unfractured rock). The actual behavior of geologic materials is more complex than in our rheological model, so the uncertainties in the plasticity do not represent the state-of-the-art error. Additional information is contained in the original.

Dombard, A. J.↗
NASA NTRSDate not supplied

Investigating Material Behavior in Atmospheric Entry Conditions: Arc-Jet Testing Insights from Meteorite Ablation to High-Temperature Coatings.

Arc-jet testing is an indispensable tool in elucidating the complex interactions materials undergo when subjected to the extreme thermal and mechanical stresses of atmospheric entry. By replicating these harsh conditions, arc-jets enable researchers to explore critical phenomena such as ablation, melting, and spallation, generating invaluable data that informs the development of sophisticated, physics-based models. These models are pivotal not only in predicting material performance for spacecraft re-entry systems but also in assessing the threat posed by celestial bodies as they encounter Earth’s atmosphere. One of the leading facilities for such investigations is the Hypersonic Materials Environmental Test System (HyMETS) at NASA’s Langley Research Center. A pioneering test campaign conducted at HyMETS focused on unraveling the ablation mechanisms of an ordinary chondrite meteorite (Tamdakht H5) and a terrestrial analog (basalt). These studies unveiled distinct material behaviors: Tamdakht demonstrated a remarkably stable melt flow, with mass loss primarily driven by the volatilization of elements such as iron and sodium, whereas basalt exhibited more aggressive surface degradation due to the rapid decomposition of hydrated minerals, culminating in significant spallation and surface material ejection. Beyond meteorite analysis, HyMETS has also been integral in the exploration of advanced Thermal Protection Systems (TPS) such as Phenolic Impregnated Carbon Ablator (PICA), a material with extensive flight heritage in missions like Stardust, Mars Science Laboratory, Mars 2020, and Osiris-REx. PICA, when coated with NuSil CV-1144-0—a polysiloxane resin designed to prevent particle shedding—undergoes a remarkable transformation upon heating. The resin pyrolyzes to form a thin, oxidation-resistant silicon oxycarbide layer that profoundly influences the material’s thermal response through a sophisticated four-stage process. Initially, the silicon oxycarbide acts as a formidable barrier, impeding reactive interactions with the PICA char and effectively suppressing surface temperature rise and material recession. However, as the heating persists, the protective layer decomposes via carbothermal reduction, exposing the underlying char, which in turn leads to a dramatic spike in surface temperature and accelerated material erosion. In the final stage, equilibrium is reached, with recession rates aligning with those of the virgin material. This intricate understanding of PICA-NuSil behavior under extreme conditions offers valuable insights, enhancing material response models and advancing the development of next-generation TPS for space exploration.

Chemistry↗
NASA NTRS1996

Buoyancy-Driven Heat Transfer During Application of a Thermal Gradient for the Study of Vapor Deposition at Low Pressure Using and Ideal Gas

A mathematical model has been developed to determine heat transfer during vapor deposition of source materials under a variety of orientations relative to gravitational accelerations. The model demonstrates that convection can occur at total pressures as low as 10-2 mm Hg. Through numerical computation, using physical material parameters of air, a series of time steps demonstrates the development of flow and temperature profiles during the course of vapor deposition. These computations show that in unit gravity vapor deposition occurs by transport through a fairly complicated circulating flow pattern when applying heat to the bottom of the vessel with parallel orientation with respect to the gravity vector. The model material parameters for air predict the effect of kinematic viscosity to be of the same order as thermal diffusivity, which is the case for Prandtl number approx. 1 fluids. Qualitative agreement between experiment and the model indicates that 6-(2-methyl-4-nitroanilino)-2,4-hexadiyn-l-ol (DAMNA) at these pressures indeed approximates an ideal gas at the experiment temperatures, and may validate the use of air physical constants. It is apparent that complicated nonuniform temperature distribution in the vapor could dramatically affect the homogeneity, orientation, and quality of deposited films. The experimental test i's a qualitative comparison of film thickness using ultraviolet-visible spectroscopy on films generated in appropriately oriented vapor deposition cells. In the case where heating of the reaction vessel occurs from the top, deposition of vapor does not normally occur by convection due to a stable stratified medium. When vapor deposition occurs in vessels heated at the bottom, but oriented relative to the gravity vector between these two extremes, horizontal thermal gradients induce a complex flow pattern. In the plane parallel to the tilt axis, the flow pattern is symmetrical and opposite in direction from that where the vessel is positioned vertically. The ground-based experiments are sufficient preliminary tests of theory and should be of significant interest regarding vapor deposited films in microgravity.

Frazier, D. O.↗