Search NASA⌕ Search

SEARCH · Search NASA

Results for “gap model”

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 505 records · Page 28

Thermal barrier coating life modeling in aircraft gas turbine engines

Analytical models for predicting ceramic thermal barrier coating (TBC) spalling life in aircraft gas turbine engines are presented. Electron beam-physical vapor deposited (EB-PVD) and plasma sprayed TBC systems are discussed. An overview of the following TBC spalling mechanisms is presented: metal oxidation at the ceramic-metal interface, ceramic-metal interface stress singularities at edges and corners, ceramic-metal interface stresses caused by radius of curvature and interface roughness, material properties and mechanical behavior, temperature gradients, component design features and object impact damage. TBC spalling life analytical models are proposed based on observations of TBC spalling and plausible failure theories. TBC spalling was assumed to occur when the imposed stresses exceed the material strength (at or near the ceramic-metal interface). TBC failure knowledge gaps caused by lack of experimental evidence and analytical understanding are noted. The analytical models are considered initial engineering approaches that capture observed TBC failure trends.

Nissley, David M.↗

Impact of AMS-02 Measurements on Reducing GCR Model Uncertainties

For vehicle design, shield optimization, mission planning, and astronaut risk assessment, the exposure from galactic cosmic rays (GCR) poses a significant and complex problem both in low Earth orbit and in deep space. To address this problem, various computational tools have been developed to quantify the exposure and risk in a wide range of scenarios. Generally, the tool used to describe the ambient GCR environment provides the input into subsequent computational tools and is therefore a critical component of end-to-end procedures. Over the past few years, several researchers have independently and very carefully compared some of the widely used GCR models to more rigorously characterize model differences and quantify uncertainties. All of the GCR models studied rely heavily on calibrating to available near-Earth measurements of GCR particle energy spectra, typically over restricted energy regions and short time periods. In this work, we first review recent sensitivity studies quantifying the ions and energies in the ambient GCR environment of greatest importance to exposure quantities behind shielding. Currently available measurements used to calibrate and validate GCR models are also summarized within this context. It is shown that the AMS-II measurements will fill a critically important gap in the measurement database. The emergence of AMS-II measurements also provides a unique opportunity to validate existing models against measurements that were not used to calibrate free parameters in the empirical descriptions. Discussion is given regarding rigorous approaches to implement the independent validation efforts, followed by recalibration of empirical parameters.

Slaba, T. C.↗

Aera-Mip: Emission Pathways, Remaining Budgets, and Carbon Cycle Dynamics Compatible With 1.5 and 2 °C Global Warming Stabilization

While international climate policies now focus on limiting global warming to well below 2 °C or pursuing a 1.5 °C level of global warming, the climate modelling community has not provided an experimental design in which all Earth system models (ESMs) converge and stabilize at the same prescribed global warming levels. This gap hampers accurate estimations based on comprehensive ESMs of the carbon emission pathways and budgets needed to meet such agreed warming levels and of the associated climate impacts under temperature stabilization. Here, we apply the Adaptive Emission Reduction Approach (AERA) with ESMs to provide such simulations in which all models converge at 1.5 and 2.0 °C warming levels by adjusting their emissions over time. These emission-driven simulations provide a wide range of emission pathways and resulting atmospheric CO 2 projections for a given warming level, uncovering uncertainty ranges that were previously missing in the traditional Coupled Model Intercomparison Project (CMIP) scenarios with prescribed greenhouse gas concentration pathways. Meeting the 1.5 °C warming level requires a 40 % (full model range: 7 % to 76 %) reduction in multi-model mean CO 2 -forcing-equivalent (CO 2 -fe) emissions from 2025 to 2030, a 98 % (57 % to 127 %) reduction from 2025 to 2050, and a stabilization at 1.0 (−1.7 to 2.9) PgC yr−1 from 2100 onward after the 1.5 °C global warming level is reached. Meeting the 2.0 °C warming level requires a 47 % (8 % to 92 %) reduction in multi-model mean CO 2 -fe emissions until 2050 and a stabilization at 1.7 (−1.5 to 2.7) PgC yr−1 from 2100 onward. The on-average positive emissions under stabilized global temperatures are the result of a decreasing transient climate response to cumulative CO 2 -fe emissions over time under stabilized global warming. This evolution is consistent with a slightly negative zero emissions commitment – initially assumed to be zero – and leads to an increase in the post-2025 CO 2 -fe emission budget by a factor of 2.2 (−0.8 to 6.9) by 2150 for the 1.5 °C warming level and a factor of 1.4 (0.9 to 2.4) for the 2.0 °C warming level compared to its first estimate in 2025. The median CO 2 -only carbon budget by 2150, relative to 2020, is 800 GtCO 2 for the 1.5 °C warming level and 2250 GtCO 2 for the 2.0 °C warming level. These median values exceed the median IPCC AR6 estimates by 60 % for the 1.5 °C warming level and 67 % for 2.0 °C. Some of the differences may be explained by the choice of the mitigation scenario for non-CO 2 radiative agents. Our simulations highlight shifts in carbon uptake dynamics under stabilized temperature, such as a cessation of the carbon sinks in the North Atlantic and in tropical forests. On the other hand, the Southern Ocean remains a carbon sink centuries after temperatures stabilize. Overall, this new type of warming-level-based emission-driven simulation offers a more coherent assessment across climate models and opens up a wide range of possibilities for studying both the carbon cycle and climate impacts, such as extreme events, under climate stabilization.

Adaptive Emission Reduction Approach↗

Thermal barrier coating life modeling in aircraft gas turbine engines

Analytical models useful for predicting ceramic thermal barrier coating (TBC) spalling life in aircraft gas turbine engines are presented. Electron beam-physical vapor deposited (EB-PVD) and plasma sprayed TBC systems are discussed. TBC spalling was attributed to a combination of mechanisms such as metal oxidation at the ceramic-metal interface, ceramic-metal interface stress concentrations at free surfaces due to dissimilar materials, ceramic-metal interface stresses caused by local radius of curvature and interface roughness, material properties and mechanical behavior, transient temperature gradients across the ceramic layer and component design features. TBC spalling life analytical models were developed based on observations of TBC failure modes and plausible failure theories. TBC failure was assumed to occur when the imposed stresses exceeded the material strength (at or near the ceramic-metal interface). TBC failure knowledge gaps caused by lack of experimental evidence and analytical understanding are noted. The analytical models are considered initial engineering approaches that capture observed TBC failure trends.

Nissley, D. M.↗

Active Turbine Tip Clearance Control Trade Space Analysis of an Advanced Geared Turbofan Engine

Tip clearance within the high pressure turbine of a gas turbine engine is a significant factor in engine performance and efficiency. In the pursuit of higher efficiency, aero-engine designs are migrating toward compact gas turbine (CGT) technology that seeks to increase the bypass ratio of the gas turbine engine without increasing the size of the fan, which is constrained by its underwing location. The reduced size of CGTs invoke concern over increased sensitivity of engine performance due to turbine tip clearance gap that makes an argument for advanced tip clearance mitigation and control techniques to be employed. This paper evaluates the tip clearance trade space for a conceptual geared turbofan engine with a CGT core. This is accomplished through a modeling and simulation approach that includes a sensitivity analysis of engine performance in response to high pressure turbine tip clearance as well as an evaluation of the sensitivity of tip clearance to various design parameters, including material properties and component cooling characteristics. Also included is a parametric study of actuators that provides preliminary requirements for implementation of active turbine tip clearance control actuation systems. The results produced from these studies are meant to be informative, with special emphasis on the demonstration of a systematic approach. The modeling approach appears to capture expected trends. The studies suggest that the tip clearance gap will have a greater impact on the new CGT engines and that a relatively slow, actively controlled actuation system may be sufficient as long as it has control authority to both open and close the tip clearance gap.

tip clearance control↗

The Effects of a 1998 Observing System Change on MERRA-2-Based Ozone Profile Simulations

Model simulations of ozone (O3) driven by meteorological reanalyses are useful for filling observational gaps and interpreting observed O3 variability and trends. However, the transport circulation of reanalysis products is impacted by changes to the observing system (the data assimilated into the reanalyses). We examine the impacts of these changes on simulated O3 from two models, Global Modeling Initiative (GMI) Chemistry Transport Model (GMI CTM) and Modern‐Era Retrospective Analysis for Research Applications version 2 (MERRA‐2) GMI Replay (M2 GMI Replay) simulation, using observations from global ozonesondes (>50,000 profiles) and satellites from 1980 to 2016. Both models are constrained by meteorology from the NASA MERRA‐2 reanalysis, and both use versions of NASA's GMI chemical mechanism. We focus on an observing system change affecting simulated O3 after 1998, associated with the assimilation of temperature and humidity data from new microwave profiling satellites. A large post‐1998 O3 increase, mainly confined to 15–20‐km altitude, of ~10 Dobson units (DU) in midlatitudes occurs in the GMI CTM, worsening the bias compared to observations. In contrast, an increase in M2 GMI Replay simulation O3 of ~10 DU is observed only near −60° latitude, reducing the bias compared to observations. The GMI CTM O3 high biases display a Quasi‐Biennial Oscillation (QBO)‐like periodicity that result from excessive transport from the tropical stratosphere to the midlatitude lower stratosphere during the QBO westerly phase. We quantify O3 discontinuities caused by MERRA‐2 observing system changes and demonstrate how the MERRA‐2 Global Modeling Initiative Replay simulation dampens the effects of these changes and QBO‐driven artifacts on simulated lower stratospheric and total O3. We caution against using simulations driven by a reanalysis to derive multidecadal O3 trends, especially prior to 1998.

Replay↗

Evidence of a Discontinuous Disk Structure Around the Herbig Ae Star HD 139614

The formation and evolution of a planetary system are intrinsically linked to the evolution of the primordial accretion disk and its dust and gas content. A new class of pre-main sequence objects has been recently identified as pre-transitional disks. They present near-infrared excess coupled to a flux deficit at about 10 microns and a rising mid-infrared and far-infrared spectrum. These features suggest a disk structure with inner and outer dust components, separated by a dust-depleted region (or gap). This could be the result of particular planet formation mechanisms that occur during the disk evolution. We here report on the first interferometric observations of the disk around the Herbig Ae star HD 139614. Its infrared spectrum suggests a flared disk, and presents pre-transitional features, namely a substantial near-infrared excess accompanied by a dip around 6 microns and a rising mid-infrared part. In this framework, we performed a study of the spectral energy distribution (SED) and the mid-infrared VLTI/MIDI interferometric data to constrain the spatial structure of the inner dust disk region and assess its possibly multi-component structure. We based our work on a temperature-gradient disk model that includes dust opacity.While we could not reproduce the SED and interferometric visibilities with a one-component disk, a better agreement was obtained with a two-component disk model composed of an optically thin inner disk extending from 0.22 to 2.3 AU, a gap, and an outer temperature-gradient disk starting at 5.6 AU. Therefore, our modeling favors an extended and optically thin inner dust component and in principle rules out the possibility that the near-infrared excess originates only from a spatially confined region. Moreover, the outer disk is characterized by a very steep temperature profile and a temperature higher than 300 K at its inner edge. This suggests the existence of a warm component corresponding to a scenario where the inner edge of the outer disk is directly illuminated by the central star. This is an expected consequence of the presence of a gap, thus indicative of a "pre-transitional" structure.

primordial accretion disk↗

Large-scale spectral structure with a gap in the stably stratified atmosphere

The analysis of the large-scale structure of turbulence in a stratified atmosphere is important for the prediction and modelling of turbulence. The Monin-Obukhov similarity consideration is well suited for the modelling of small-scale turbulence. But the large-scale turbulence possesses new features and parameters, as connected with the appearance of a spectral gap in the stably stratified atmosphere, for the understanding of which an analytical investigation is needed. By a group-kinetic method, the anisotropic form of the transfer function and its governing eddy viscosity are derived. On this basis, the spectral laws k to the -5/3rd, k to the -1st, k to the -3rd are obtained for the subranges of inertia, shear, and buoyancy, respectively. A gap connects the last two subranges by the variable anisotropy, as a result of the fall of the spectrum of the vertical velocity fluctuations. The variable anisotropy causes the deviations from the spectral laws k to the -1st and k to the -3rd. An experiment supports the analytical prediction.

Larsen, S. E.↗

Towards Co-evolution of Membrane Transport and Metabolism

Protocellular boundaries were inextricably connected to the metabolism they encapsulated: to be inheritable, early metabolism must have led to an increased rate of growth and division of vesicles and, similarly, transport through vesicle boundaries must have supported the evolution of metabolism. Even though explaining how this coupling emerged and evolved in the absence of the complex machinery of modern cells is one of the key issues in studies on the origin of life, little is known about the biochemical and biophysical processes that might have been involved. This gap in our knowledge is a major impediment in efforts to construct scenarios for the origin of life and laboratory models of protocells. A combination of experimental and computational studies carried out by us and our collaborators is aimed at helping to close this gap. Properties of membranes might have contributed to the selection of RNA as an early biopolymer. A kinetic mechanism was proposed (Sacerdote & Szostak, 2005) in which ribose was supplied more quickly than other aldopentoses to primordial cells for preferential incorporation of ribonucleotides into nucleic acids. This proposal is based on a finding that ribose permeates membranes an order of magnitude faster than its diastereomers, arabinose and xylose. Our computer simulations, which yield permeation rates in excellent agreement with experiment, and kinetic modeling explain this phenomenon in terms of inter- and intramolecular interactions involving exocyclic hydroxyl groups attached to carbon atoms of the pyranose ring (Wei and Pohorille, 2009). They also constrain scenarios for the formation of the earliest nucleic acids (Wei and Pohorille, 2013). In one scenario, sugars permeate protocellular walls and subsequently are used to synthesize nucleic acids inside protocells. As long as this process proceeds at the rate faster than 6x10(exp -3)/s, ribose derivatives will be available for synthesis easier than their diastereomers. If nucleosides or their activated derivatives are synthesized outside protocells and subsequently transported across protocellular membranes the kinetic mechanism does not apply because all diastereomers, which have their sugars in the furanose rather than pyranose form, permeate the membrane at approximately the same rate. Properties of membranes might have been also coupled to metabolism involving peptides. Recently, Adamala and Szostak (2013) have shown that a dipeptide inside fatty-acid vesicles catalyzes the formation of another dipeptide that binds to vesicle walls and, by doing so, promotes their growth at the expense of other vesicles. This coupling of metabolism, permeability of vesicles and their growth is the first demonstration of evolutionary advantage imparted by small, membrane-bound peptides. Building on this work we have calculated the rate at which different blocked amino acids are delivered to a protocell for synthesis of dipeptides. We have further shown that the dipeptides are located at the water-membrane interface rather than in the center of the bilayer. On these basis it is anticipated that other dipeptides containing aromatic, but not necessarily hydrophobic amino acids (e.g. tyrosine) could have the same catalytic effects. Insight from these studies allows for estimating the rate of vesicle growth and the rates of dipeptide synthesis required to keep the system in balance. These results, in combination with our earlier studies, lead to a general scenario for evolution from membrane-bound dipeptides to ion channels in the origin of life.

Metabolism↗

Bulk Growth of Wide Band Gap II-VI Compound Semiconductors by Physical Vapor Transport

The mechanism of physical vapor transport of II-VI semiconducting compounds was studied both theoretically, using a one-dimensional diffusion model, as well as experimentally. It was found that the vapor phase stoichiometry is critical in determining the vapor transport rate. The experimental heat treatment methods to control the vapor composition over the starting materials were investigated and the effectiveness of the heat treatments was confirmed by partial pressure measurements using an optical absorption technique. The effect of residual (foreign) gas on the transport rate was also studies theoretically by the diffusion model and confirmed experimentally by the measurements of total pressure and compositions of the residual gas. An in-situ dynamic technique for the transport rate measurements and a further extension of the technique that simultaneously measured the partial pressures and transport rates were performed and, for the first time, the experimentally determined mass fluxes were compared with those calculated, without any adjustable parameters, from the diffusion model. Using the information obtained from the experimental transport rate measurements as guideline high quality bulk crystal of wide band gap II-VI semiconductor were grown from the source materials which undergone the same heat treatment methods. The grown crystals were then extensively characterized with emphasis on the analysis of the crystalline structural defects.

Su, Ching-Hua↗

Bulk Growth of Wide Band Gap II-VI Compound Semiconductors by Physical Vapor Transport

The mechanism of physical vapor transport of II-VI semiconducting compounds was studied both theoretically, using a one-dimensional diffusion model, as well as experimentally. It was found that the vapor phase stoichiometry is critical in determining the vapor transport rate. The experimental heat treatment methods to control the vapor composition over the starting materials were investigated and the effectiveness of the heat treatments was confirmed by partial pressure measurements using an optical absorption technique. The effect of residual (foreign) gas on the transport rate was also studies theoretically by the diffusion model and confirmed experimentally by the measurements of total pressure and compositions of the residual gas. An in-situ dynamic technique for the transport rate measurements and a further extension of the technique that simultaneously measured the partial pressures and transport rates were performed and, for the first time, the experimentally determined mass fluxes were compared with those calculated, without any adjustable parameters, from the diffusion model. Using the information obtained from the experimental transport rate measurements as guideline high quality bulk crystal of wide band gap II-VI semiconductor were grown from the source materials which undergone the same heat treatment methods. The grown crystals were then extensively characterized with emphasis on the analysis of the crystalline structural defects.

Su, Ching-Hua↗

The Land Surface in Earth System Reanalysis

Current global atmospheric reanalysis products such as ERA-5, MERRA-2, and JRA-55 provide estimates of land surface states and fluxes, including soil moisture, soil temperature, snow mass, latent and sensible heat fluxes, and runoff. These land surface estimates are based on land surface process models and, depending on the reanalysis product, on precipitation observations or the assimilation of land surface observations of soil moisture and snow conditions from satellites and in situ networks. The land models used in the current global atmospheric reanalysis systems are generally focused on water and energy balance processes, with prescribed static land cover and prescribed climatological vegetation conditions. That is, current global atmospheric reanalysis products do not reflect changes in land use or land cover, nor do they provide dynamic estimates of the land carbon cycle. These and other gaps, including the omission of anthropogenic processes such as irrigation in the land model, need to be addressed in the development of the next generation Earth system reanalysis systems. This presentation provides an overview of the state-of-the-art and future developments of land surface estimates in global reanalysis systems.

Rolf Reichle↗

Spectral Power-law Formation by Sequential Particle Acceleration in Multiple Flare Magnetic Islands

We present a first-principles model of pitch-angle and energy distribution function evolution as particles are sequentially accelerated by multiple flare magnetic islands. Data from magnetohydrodynamic (MHD) simulations of an eruptive flare/coronal mass ejection provide ambient conditions for the evolving particle distributions. Magnetic islands, which are created by sporadic reconnection at the self-consistently formed flare current sheet, contract and accelerate the particles. The particle distributions are evolved using rules derived in our previous work. In this investigation, we assume that a prescribed fraction of particles sequentially "hops" to another accelerator and receives an additional boost in energy and anisotropy. This sequential process generates particle number spectra that obey an approximate power law at mid-range energies and presents low- and high-energy breaks. We analyze these spectral regions as functions of the model parameters. We also present a fully analytic method for forming and interpreting such spectra, independent of the sequential acceleration model. The method requires only a few constrained physical parameters, such as the percentage of particles transferred between accelerators, the energy gain in each accelerator, and the number of accelerators visited. Our investigation seeks to bridge the gap between MHD and kinetic regimes by combining global simulations and analytic kinetic theory. The model reproduces and explains key characteristics of observed flare hard X-ray spectra as well as the underlying properties of the accelerated particles. Our analytic model provides tools to interpret high-energy observations for missions and telescopes, such as RHESSI, FOXSI, NuSTAR, Solar Orbiter, EOVSA, and future high-energy missions.

magnetic reconnection↗

A Conceptual Framework to Integrate Biodiversity, Ecosystem Function, and Ecosystem Service Models

Global biodiversity and ecosystem service models typically operate independently. Ecosystem service projections may therefore be overly optimistic because they do not always account for the role of biodiversity in maintaining ecological functions. We review models used in recent global model intercomparison projects and develop a novel model integration framework to more fully account for the role of biodiversity in ecosystem function, a key gap for linking biodiversity changes to ecosystem services. We propose two integration pathways. The first uses empirical data on biodiversity–ecosystem function relationships to bridge biodiversity and ecosystem function models and could currently be implemented globally for systems and taxa with sufficient data. We also propose a trait-based approach involving greater incorporation of biodiversity into ecosystem function models. Pursuing both approaches will provide greater insight into biodiversity and ecosystem services projections. Integrating biodiversity, ecosystem function, and ecosystem service modeling will enhance policy development to meet global sustainability goals.

biodiversity models↗

Model and full-scale large transport airframe noise

Results are presented for an airframe noise experiment intended to obtain airframe noise data of a 0.03-scale model of a large transport aircraft in a flow facility for comparison with similar airframe noise data from full-scale flyover tests of the same aircraft. The test model simulated the landing and cruise configuration, including individual flap systems. Considerable noise is found to be associated with the extension of the model's leading edge flap system and to be reduced across the noise spectra when the gaps between the leading-edge flaps and the leading edge of the wing are sealed. Amplitude and frequency scaling relationships are developed which correlate respective data for a variety of aircraft configurations. Normalized model landing gear noise failed to simulate full-scale data.

Shearin, J. G.↗

Extremely low nonalloyed and alloyed contact resistance using an InAs cap layer on InGaAs by molecular-beam epitaxy

Extremely low alloyed and nonalloyed ohmic contact resistances have been formed on n-type InAs/In(0.53)Ga(0.47)As/In(0.52)Al(0.48)As structures grown on InP(Fe) by molecular-beam epitaxy. To insure the accuracy of the small contact resistances measured, an extended transmission line model was used to extrapolate contact resistances from test patterns with multiple gap spacings varying from 1 to 20 microns. For a 150-A-thick InAs layer doped to 2 x 10 to the 18th/cu cm and a 0.1-micron-thick InGaAs layer doped to 1 x 10 to the 18th/cu cm, a specific contact resistance of 2.6 x 10 to the -8th ohm-asterisk sq cm was measured for the nonalloyed contact, while a resistance less than 1.7 x 10 to the -8th ohm-asterisk sq cm is reported for the alloyed contact. Conventional Au-Ge/Ni/Au was used for the ohmic metal contact and alloying was performed at 500 C for 50 s in flowing H2. Using a thermionic field emission model, the barrier height at the InAs/InGaAs interface was calculated to be 20 meV.

Peng, C. K.↗

Climate Change and a Global City: An Assessment of the Metropolitan East Coast Region

The objective of the research is to derive an assessment of the potential climate change impacts on a global city - in this case the 31 county region that comprises the New York City metropolitan area. This study comprises one of the regional components that contribute to the ongoing U.S. National Assessment: The Potential Consequences of Climate Variability and Change and is an application of state-of-the-art climate change science to a set of linked sectoral assessment analyses for the Metro East Coast (MEC) region. We illustrate how three interacting elements of global cities react and respond to climate variability and change with a broad conceptual model. These elements include: people (e.g., socio- demographic conditions), place (e.g., physical systems), and pulse (e.g., decision-making and economic activities). The model assumes that a comprehensive assessment of potential climate change can be derived from examining the impacts within each of these elements and at their intersections. Thus, the assessment attempts to determine the within-element and the inter-element effects. Five interacting sector studies representing the three intersecting elements are evaluated. They include the Coastal Zone, Infrastructure, Water Supply, Public Health, and Institutional Decision-making. Each study assesses potential climate change impacts on the sector and on the intersecting elements, through the analysis of the following parts: 1. Current conditions of sector in the region; 2. Lessons and evidence derived from past climate variability; 3. Scenario predictions affecting sector; potential impacts of scenario predictions; 4. Knowledge/information gaps and critical issues including identification of additional research questions, effectiveness of modeling efforts, equity of impacts, potential non-local interactions, and policy recommendations; and 5. Identification of coping strategies - i.e., resilience building, mitigation strategies, new technologies, education that affects decision-making, and better preparedness for contingencies.

Rosenzweig, Cynthia↗

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↗