Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “approach temperature”

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 361 records · Page 20
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Effects of midplane density gradient in the propagation of high-harmonic fast waves considering high temperature anisotropy in the scrape-off layer of NSTX-U

High-Harmonic Fast Wave (HHFW) heating experiments in NSTX have shown that up to 60% of the injected power can be lost in the Scrape-Off Layer (SOL) when the density is above the fast wave cutoff density in front of the antenna, for which the fast wave is able to propagate into the plasma. This work models HHFW propagation in the SOL plasmas of NSTX-U using a 2D divertor SOL profile derived from the pressure assumption and the finite element temperature solution, which accounts for the high anisotropy of heat conduction in a 2D axisymmetric geometry. In this work, the two-dimensional axisymmetric SOL temperature profile is first evaluated by solving the steady-state non-linear heat conduction equation, in which thermal conductivity depends on temperature, using a finite element approach in the Petra-M workbench. A 2D density profile is then obtained from a prescribed density profile at the outer midplane assuming pressure is uniform along a flux tube in the SOL. This approach results in density and temperature profiles in which the strong asymmetric nature of heat conduction between the parallel and perpendicular background magnetic field is successfully captured. Furthermore, this work focuses on investigating the effect of the SOL plasma density profile on parasitic HHFW propagation in the SOL. The simulation results show that the radial gradient of the density profile affects the wavefield propagation in the SOL. As the density profile broadens, the wavefield intensity and its poloidal extent are reduced in the SOL, and the core coupling increases. Similarly, as shown by collisional power deposition, a proxy for power absorption, the fraction of power deposited in the SOL decreases with a broadening profile.

HHFW↗
NASA NTRS1994

Modeling of outgassing and matrix decomposition in carbon-phenolic composites

Work done in the period Jan. - June 1994 is summarized. Two threads of research have been followed. First, the thermodynamics approach was used to model the chemical and mechanical responses of composites exposed to high temperatures. The thermodynamics approach lends itself easily to the usage of variational principles. This thermodynamic-variational approach has been applied to the transpiration cooling problem. The second thread is the development of a better algorithm to solve the governing equations resulting from the modeling. Explicit finite difference method is explored for solving the governing nonlinear, partial differential equations. The method allows detailed material models to be included and solution on massively parallel supercomputers. To demonstrate the feasibility of the explicit scheme in solving nonlinear partial differential equations, a transpiration cooling problem was solved. Some interesting transient behaviors were captured such as stress waves and small spatial oscillations of transient pressure distribution.

Mcmanus, Hugh L.↗
NASA NTRS2006

SIM PlanetQuest: Mili-Kelvin Analysis of the Collector Sub-system Thermal Model

Space Interferometry Mission (SIM) PlanetQuest is part of the NASA Astronomical Search for Origins (ASO) program. It will use a Michelson Interferometer to locate extrasolar planets, among other objectives. The spacecraft's subsystems, temperature requirements, and thermal control approach are discussed. The current SIM thermal design is analyzed using a Milli-Kelvin approach and appears to meet temperature stability requirements on the OOPCC optic. SIM absolute temperature requirements are met through the robustness of the thermal design (i.e., PIC control).

PlanetQuest↗
DOE OSTI2023

Influence of the air-cooled supercritical carbon dioxide recompression cycle design-point on molten salt power tower techno-economic performance

Here, this paper investigates the techno-economic performance of a molten salt power tower integrated with an air-cooled supercritical carbon dioxide recompression cycle. Prior research into the cycle design point showed trade-offs between efficiency, cost, and temperature difference across the heat input, resulting in trade-offs between the design-point generation of the molten salt power tower system and the cost of sensible-heat thermal energy storage. Prior research into cycle off-design showed that colder ambient temperature designs have higher efficiencies on cold days but lower efficiency and capacity degradation on hot days. We present a framework to evaluate how cycle design, cost, and corresponding off-design performance affect the techno-economic performance of the molten salt power tower system. We use a steady-state off-design cycle model assuming fixed shaft speeds, inventory control, and air-cooler fan power control and integrate the cycle model with the molten salt power tower model from the National Renewable Energy Laboratory's System Advisor Model. First, we present a surrogate model that uses a heuristic approach to correct the interpolation when the cycle reaches capacity limitations at ambient temperatures hotter than design. Next, we evaluate how the recompression cycle design-point ambient temperature, air-cooler approach temperature, air-cooler fan power, and recuperator conductance influence cycle and system design metrics. Finally, we complete annual simulations using an hourly "duck curve" pricing schedule and present annual performance metrics to evaluate the cumulative effects of the cycle design parameters and corresponding off-design performance on molten salt power tower techno-economics.

14 SOLAR ENERGY↗
DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗
DOE OSTI2026

An International Laboratory Comparison Study on Approximating the Enthalpy of Adsorption via the Clausius‐Clapeyron Approach

Materials-based gas capture and storage is an increasingly important area of research. Robust and accurate determination of material properties is required for judicial selection of materials for specific applications and for engineering materials–based systems at scale. One key property is the strength of the adsorbate–adsorbent interaction often quantified via the isosteric enthalpy of adsorption. The heat of adsorption can be measured directly through calorimetry; however, a more widely used approach is to apply the Clausius-Clapeyron (CC) equation to adsorption isotherms collected at different temperatures. While this approach appears to be straightforward, there exist multiple variants in the application of the methodologies employed. This raises the question on how these variations may or may not affect the determined results. Presented here is a discussion of the most common methodologies and a comparison of indirect determinations (via CC) of the isosteric enthalpy of adsorption by different laboratories on identical material. Included in that comparison are discussions on the measurement and analysis reproducibility. Importantly, details of the methodologies are shown to be critical when comparing enthalpies among laboratories, and different methodologies contribute to significant discrepancies and artifacts in the results. Recommendations are provided to promote robust determination and the reporting thereof.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

The Role of Land‐Atmosphere Feedbacks in Midlatitude Wintertime Surface Temperature Variability

Accurately representing synoptic near-surface temperature variability is crucial to predict weather extremes, yet models remain biased. Existing studies primarily attribute wintertime midlatitude near-surface temperature variability to tropospheric large-scale advection, assuming minimal land influence. However, nudging the model's circulation toward observations yields little improvement in wintertime temperature variance over Northern Hemisphere land, suggesting that land-atmosphere interactions also warrant attention. We introduce a new scaling framework for temperature variance that incorporates local land-atmosphere feedbacks. Comparing our framework to the mixing length approach—which links temperature variance to the meridional temperature gradient and air parcel displacement (mixing length)—shows that land-atmosphere feedbacks are inherently embedded in the mixing length, a connection previously overlooked. Roles of land–atmosphere feedbacks are evaluated via model experiments with perturbed circulation, land, or both. We find that longwave radiative damping dominates temperature variance responses over meridional temperature gradient when both land and circulation are perturbed.

atmosheric science↗
NASA NTRS1980

Real-gas effects 1: Simulation of ideal gas flow by cryogenic nitrogen and other selected gases

The thermodynamic properties of nitrogen gas do not thermodynamically approximate an ideal, diatomic gas at cryogenic temperatures. Choice of a suitable equation of state to model its behavior is discussed and the equation of Beattie and Bridgeman is selected as best meeting the needs for cryogenic wind tunnel use. The real gas behavior of nitrogen gas is compared to an ideal, diatomic gas for the following flow processes: isentropic expansion; normal shocks; boundary layers; and shock wave boundary layer interactions. The only differences in predicted pressure ratio between nitrogen and an ideal gas that may limit the minimum operating temperatures of transonic cryogenic wind tunnels seem to occur at total pressures approaching 9atmospheres and total temperatures 10 K below the corresponding saturation temperature, where the differences approach 1 percent for both isentropic expansions and normal shocks. Several alternative cryogenic test gases - air, helium, and hydrogen - are also analyzed. Differences in air from an ideal, diatomic gas are similar in magnitude to those of nitrogen. Differences for helium and hydrogen are over an order of magnitude greater than those for nitrogen or air. Helium and hydrogen do not approximate the compressible flow of an ideal, diatomic gas.

Hall, R. M.↗
DOE OSTI2025

Interactions Enhance Ramp Reversal Memory in Locally Phase Separated Materials

The ramp-reversal memory (RRM) effect in metal–insulator transition metal oxides (TMOs), a non-volatile resistance change induced by repeated temperature cycling, has attracted considerable interest in neuromorphic computing and non-volatile memory devices. Our previous defect motion model successfully explained RRM in vanadium dioxide (VO 2 ), capturing observed critical temperature shifts and memory accumulation throughout the sample. However, this approach lacked interactions between metallic and insulating domains. Here, we extend our model by combining a correlated Random Field Ising Model with defect diffusion-segregation, enabling accurate hysteresis modeling while predicting the relationship between RRM and domain interactions. Our simulations demonstrate that the maximum RRM occurs when the turnaround temperature approaches the inflection point. This peak in RRM vs. turnaround temperature is consistent with prior transport measurements, as well as our own optical measurements reported here. Significantly, we find that increasing nearest-neighbor interactions enhances the maximum memory effect, thus providing a clear mechanism for optimizing RRM performance. Since our model employs minimal assumptions, we predict that RRM should be a widespread phenomenon in materials exhibiting patterned phase coexistence of electronic domains. This work not only advances fundamental understanding of memory behavior in TMOs but also establishes a much-needed theoretical framework for optimizing device applications.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dielectric polymers for high-temperature capacitive energy storage

Polymers are the preferred materials for dielectrics in high-energy-density capacitors. The electrification of transport and growing demand for advanced electronics require polymer dielectrics capable of operating efficiently at high temperatures. In this review, we critically analyze the most recent development in the dielectric polymers for high-temperature capacitive energy storage applications. While general design considerations are discussed, emphasis is placed on the elucidation of the structural dependence of the high-field dielectric and electrical properties and the capacitive performance, including discharged energy density, charge–discharge efficiency and cyclability, of dielectric polymers at high temperatures. Here, advantages and limitations of current approaches to high-temperature dielectric polymers are summarized. Challenges along with future research opportunities are highlighted at the end of this article.

25 ENERGY STORAGE↗
NASA NTRS1999

"Ultra"-Fast Fracture Strength of Advanced Structural Ceramic Materials Studied at Elevated Temperatures

The accurate determination of inert strength is important in reliable life prediction of structural ceramic components. At ambient temperature, the inert strength of a brittle material is typically regarded as free of the effects of slow crack growth due to stress corrosion. Therefore, the inert strength can be determined either by eliminating active species, especially moisture, with an appropriate inert medium, or by using a very high test rate. However, at elevated temperatures, the concept or definition of the inert strength of brittle ceramic materials is not clear, since temperature itself is a degrading environment, resulting in strength degradation through slow crack growth and/or creep. Since the mechanism to control strength is rate-dependent viscous flow, the only conceivable way to determine the inert strength at elevated temperatures is to utilize a very fast test rate that either minimizes the time for or eliminates slow crack growth. Few experimental studies have measured the elevated-temperature, inert (or "ultra"-fast fracture) strength of advanced ceramics. At the NASA Lewis Research Center, an experimental study was initiated to better understand the "ultra"-fast fracture strength behavior of advanced ceramics at elevated temperatures. Fourteen advanced ceramics - one alumina, eleven silicon nitrides, and two silicon carbides - have been tested using constant stress-rate (dynamic fatigue) testing in flexure with a series of stress rates including the "ultra"-fast stress rate of 33 000 MPa/sec with digitally controlled test frames. The results for these 14 advanced ceramics indicate that, notwithstanding possible changes in flaw populations as well as flaw configurations because of elevated temperatures, the strength at 33 000 MPa/sec approached the room-temperature strength or reached a higher value than that determined at the conventional test rate of 30 MPa/sec. On the basis of the experimental data, it can be stated that the elevated-temperature, inert strength of an advanced ceramic material can be defined as the strength where no slow crack growth takes place at the temperature.

Choi, Sung R.↗
DOE OSTI2021

Inherent issues regarding the use of in situ x-ray diffraction measurements to determine temperature in shock-compressed metals

Temperature determination in shock compressed solids constitutes an important and long-standing scientific need. Since the reduction of Bragg diffracted peaks due to temperature increase, using the Debye-Waller factor, is well established, we examined the use of this approach to determine temperatures in shock compressed gold and platinum by representing the shocked state as a superposition of density and temperature changes. Comparison of the calculated and measured diffraction peaks did not show good agreement, because x-ray diffraction (XRD) profiles in the shock compressed state are not governed solely by density and temperature changes. XRD results are also influenced measurably by shock wave induced microstructural changes. Our results demonstrate that contributions from microstructural changes need to be incorporated and modeled in the theoretical analysis to use XRD measurements for reliable temperature determination in shock compressed solids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

The interactive indoor-outdoor building energy modeling for enhancing the predictions of urban microclimates and building energy demands

There is a lack of an urban building energy modeling framework that considers the influence of surrounding buildings and local urban climate on building thermal performance. This can lead to inaccurate results since the thermal performance of individual buildings is heavily influenced by their surrounding built and climatic environment. This study establishes an interactive indoor-outdoor building energy modeling method to enhance the predictions of urban microclimates and building energy demands by coupling an urban physics model with a physics-based building energy model. Validation of the interactive coupling scheme uses field measurement datasets. Parametric simulation and analysis are conducted to understand the influence of the roof-to-canyon width ratio, canyon orientation, and ground vegetation fraction on canyon temperature, building energy consumption, and energy demand. Furthermore, the impacts of building energy model complexity (e.g., detailed vs. simplified building models) and coupling approaches on canyon temperature and building energy profiles are demonstrated using two case study buildings. In comparison with the one-way coupling approach, cooling energy consumption predicted with the dynamic two-way coupling approach varies by 3.5% and 0.5% for the detailed medium office building model and high-rise building model, respectively, and peak cooling demand varies by 8.4% and 7.0% for the detailed medium office building model and high-rise building model, respectively. Here this study also suggests that adopting a complex two-way coupling approach with environmental data exchange at various elevations is necessary for modeling tall buildings at the urban scale.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Finite-temperature electron-capture rates for neutron-rich nuclei near N = 50 and effects on core-collapse supernova simulations

In this paper the temperature dependence of stellar electron-capture (EC) rates is investigated, with a focus on nuclei near N = 50, just above Z = 28, which play an important role during the collapse phase of core-collapse supernovae (CCSN). Two new microscopic calculations of stellar EC rates are obtained from relativistic and nonrelativistic finite-temperature quasiparticle random-phase approximation approaches, for a conventional grid of temperatures and densities. In both approaches, EC rates due to Gamow-Teller transitions are included. In the relativistic calculation, contributions from first-forbidden transitions are also included and add strongly to the EC rates. The new EC rates are compared with large-scale shell-model calculations for the specific case of 86 Kr, providing insight into the finite-temperature effects on the EC rates. At relevant thermodynamic conditions for core collapse, the discrepancies between the different calculations of this paper are within about one order of magnitude. Numerical simulations of CCSN are performed with the spherically symmetric gr1d simulation code to quantify the impact of such differences on the dynamics of the collapse. These simulations also include EC rates based on two parametrized approximations. A comparison of the neutrino luminosities and enclosed mass at core bounce shows that differences between simulations with different sets of EC rates are relatively small (≈ 5 %), suggesting that the EC rates used as inputs for these simulations have become well constrained.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1983

Evaluation of Three Constitutive Models for the Prediction of Hastelloy X Elevated Temperature Cyclic Response

The approach for the method development assumes that, for a thermally loaded structure, the overall strain history is defined by linear elastic analysis. The local stress history at a fatigue critical location is then determined from a one dimensional material behavior model and the local strain and temperature conditions. Three material models are currently being evaluated to assess their ability to predict relevant high temperature cyclic material response characteristics. They are: a time independent classical plasticity and creep representation, a time dependent viscoplastic model capable of predicting combined creep and plasticity effects, and an approximate elastic analysis approach that uses a series of stress-strain curves and a cyclic hardening model to determine reverse plasticity.

Moreno, V.↗
NASA NTRS2004

Ice Particle Impact on Cloud Water Content Instrumentation

Determining the total amount of water contained in an icing cloud necessitates the measurement of both the liquid droplets and ice particles. One commonly accepted method for measuring cloud water content utilizes a hot wire sensing element, which is maintained at a constant temperature. In this approach, the cloud water content is equated with the power required to keep the sense element at a constant temperature. This method inherently assumes that impinging cloud particles remain on the sensing element surface long enough to be evaporated. In the case of ice particles, this assumption requires that the particles do not bounce off the surface after impact. Recent tests aimed at characterizing ice particle impact on a thermally heated wing section, have raised questions about the validity of this assumption. Ice particles were observed to bounce off the heated wing section a very high percentage of the time. This result could have implications for Total Water Content sensors which are designed to capture ice particles, and thus do not account for bouncing or breakup of ice particles. Based on these results, a test was conducted to investigate ice particle impact on the sensing elements of the following hot-wire cloud water content probes: (1) Nevzorov Total Water Content (TWC)/Liquid Water Content (LWC) probe, (2) Science Engineering Associates TWC probe, and (3) Particle Measuring Systems King probe. Close-up video imaging was used to study ice particle impact on the sensing element of each probe. The measured water content from each probe was also determined for each cloud condition. This paper will present results from this investigation and attempt to evaluate the significance of ice particle impact on hot-wire cloud water content measurements.

Emery, Edward F.↗
NASA NTRS2004

Processing Techniques Developed to Fabricate Lanthanum Titanate Piezoceramic Material for High-Temperature Smart Structures

Piezoelectric ceramic materials are potential candidates for use as actuators and sensors in intelligent gas turbine engines. For piezoceramics to be applied in gas turbine engines, they will have to be able to function in temperatures ranging from 1000 to 2500 F. However, the maximum use temperature for state-of-the-art piezoceramic materials is on the order of 300 to 400 F. Research activities have been initiated to develop high-temperature piezoceramic materials for gas turbine engine applications. Lanthanum titanate has been shown to have high-temperature piezoelectric properties with Curie temperatures of T(sub c) = 1500 C and use temperatures greater than 1000 C. However, the fabrication of lanthanum titanate poses serious challenges because of the very high sintering temperatures required for densification. Two different techniques have been developed at the NASA Glenn Research Center to fabricate dense lanthanum titanate piezoceramic material. In one approach, lower sintering temperatures were achieved by adding yttrium oxide to commercially available lanthanum titanate powder. Addition of only 0.1 mol% yttrium oxide lowered the sintering temperature by as much as 300 C, to just 1100 C, and dense lanthanum titanate was produced by pressure-assisted sintering. The second approach utilized the same commercially available powders but used an innovative sintering approach called differential sintering, which did not require any additive.

Goldsby, Jon C.↗
NASA NTRSDate not supplied

Direct Reconstruction of Ablative Thermal Protection System Aeroheating Using A Green's Function Approach

A Green’s function inverse heat transfer (IHT) approach is used to reconstruct the surface heating conditions on ablative thermal protection system (TPS) materials from embedded heat flux sensor and temperature probe measurements. The approach models the temperature time-history at the measurement location as a discrete linear system, allowing for the heat flux boundary condition to be recovered directly without the need for time-marching schemes. The effects of material decomposition and pyrolysis gas transport are modeled using an energy source/sink analogue. The performance of the reconstruction approach is analyzed on a 1D test case representative of an atmospheric entry heating scenario. The approach can recover the TPS surface heat flux to within 4% of the input heating condition with a computation time of 2-3 seconds (>3 orders of magnitude faster than current time-marching IHT methods). As a byproduct of the surface heating reconstruction, the algorithm also captures the surface pyrolysis gas mass flux and solid decomposition at multiple through-thickness locations within the TPS.

Kenneth McAfee↗