Search NASA⌕ Search

SEARCH · Search NASA

Results for “Curve fitting”

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 541 records · Page 30

Investigation of Bolt Preload Relaxation for JWST Thermal Heat Strap Assembly Joints with Aluminum-1100 and Indium Gaskets

Accurately predicting fastener preload relaxation in the James Webb Space Telescope (JWST) heat strap assemblies is essential to insure adequate thermal performance during its mission lifecycle. The mechanisms for preload relaxation in the strap joints include Al-1100 material creep, indium gasket flow-out, and embedment of the joint faying surfaces. This report documents the results from a bolted joint relaxation test, including analysis and curve fitting of the test data for predicting preloads five years after initial torque application. The report also includes the derivation of a preload uncertainty factor enveloping both torque/preload application scatter and expected preload relaxation at the end of mission life.

Joint↗

Investigation of Bolt Preload Relaxation for JWST Thermal Heat Strap Assembly Joints with Aluminum-1100 and Indium Gaskets

Accurately predicting fastener preload relaxation in the James Webb Space Telescope (JWST) heat strap assemblies is essential to insure adequate thermal performance during its mission lifecycle. The mechanisms for preload relaxation in the strap joints include Al-1100 material creep, indium gasket flow-out, and embedment of the joint faying surfaces. This report documents the results from a bolted joint relaxation test, including analysis and curve fitting of the test data for predicting preloads five years after initial torque application. The report also includes the derivation of a preload uncertainty factor enveloping both torque/preload application scatter and expected preload relaxation at the end of mission life.

Preload↗

Modal Q Factor and Modal Overlap of Electrically Small Avionics Box

Estimating the electromagnetic field strength in avionics boxes and other small enclosures at the design stage requires an estimate of the Q factor of the cavity modes. When the enclosure is small, it is typically under-moded so that Q measurement techniques which are standard practice in over-moded reverberation chambers may not be a robust measure. Furthermore, practical antenna used to measure Q in a small cavity may have a strong influence on the result obtained, as reported by Tait et al, IEEE Trans. EMC 55 2 2012. This paper reports the results of testing to determine the Q factor of a small aperture enclosure, used in a statistical power balance model to predict the electric field strength. The contributors to the total Q are identified. A novel S(sub 11) curve fitting method to measure modal Q is introduced and compared with the time domain method for measuring Q.

Statistical Electromagnetics↗

Automating Microbial Directed Evolution For Bioengineering Applications

From a micro-biology perspective, directed evolution is a technique that uses controlled environmental pressures to select for a desired phenotype. Directed evolution has the distinct advantage over rational design of not needing extensive knowledge of the genome or pathways associated with a microorganism to induce phenotypes. However, there are currently limitations to the applicability of this technique including being time-consuming, error-prone, and dependent on existing assays that may lack selectivity for the given phenotype. The AADEC (Autonomous Adaptive Directed Evolution Chamber) system is a proof-of-concept instrument to automate and improve the technique such that directed evolution can be used more effectively as a general bioengineering tool. A series of tests using the automated system and comparable by-hand survival assay measurements have been carried out using UV-C radiation and Escherichia coli cultures in order to demonstrate the advantages of the AADEC versus traditional implementations of directed evolution such as random mutagenesis. AADEC uses UV-C exposure as both a source of environmental stress and mutagenesis, so in order to evaluate the UV-C tolerance obtained from the cultures, a manual UV-C exposure survival assay was developed alongside the device to compare the survival fractions at a fixed dosage. This survival assay involves exposing E. coli to UV-C radiation using a custom-designed exposure hood to control the flux and dose. Surviving cells are counted then transferred to the next iteration and so on for several iterations to calculate the survival fractions for each exposure iteration.This survival assay primarily serves as a baseline for the AADEC device, allowing quantification of the differences between the AADEC system over the manual approach. The primary data of comparison is survival fractions; this is obtained by optical density and plate counts in the manual assay and by optical density growth curve fits pre- and post-exposure in the automated case. This data can then be compiled to calculate trends over the iterations to characterize increasing UV-C resistance of the E.coli strains. The observed trends are statistically indistinguishable through several iterations from both sources.

Bioengineering↗

Blunt-Body Paradox and Improved Application of Transient-Growth Framework

The “Reshotko-Tumin transition criterion" based on optimal transient growth successfully correlates laboratory measurements of roughness induced transition over blunt body configurations. Even though transient growth has not been conclusively linked to the measured onset of transition, the above correlation denotes the only available physics-based model for subcritical transition in blunt body flows, since the latter do not support any modal instabilities at typical experimental conditions. Unlike other established models based on empirical curve fits that are valid for a specific subclass of datasets, the optimal-growth-based transition criterion appears to provide a reasonable correlation with measurements in various wind tunnel and ballistic range facilities and for a broad range of surface temperature ratios. This paper is focused on optimal growth calculations that improve upon significant shortcomings of the computations underlying the Reshotko-Tumin correlation. The improved framework is applied to leeward transition over a spherical section forebody that was tested in the Mach 6 Adjustable Contour Expansion wind tunnel at Texas A&M University. The computed results highlight the significance of nonparallel basic state evolution, curvature terms, and an optimization procedure that varies both inflow and outflow locations of the transient growth interval. More important, the results indicate that the modified correlation is very close to its original form, and hence, that the accuracy of the transient-growth-based transition criterion is not compromised by using a more thorough theoretical framework. Yet the results also show that the optimal energy gain up to the predicted transition onset location can be rather small, highlighting the need to further investigate the optimal growth criterion for additional experimental configurations and to also uncover the in-depth physics underlying blunt body transition.

Freestream conditions↗

Effect of Reactant Pressure on Proton Exchange Membrane Fuel Cell Performance

It is an objective of the National Aeronautics and Space Administration (NASA) to return to the moon and to create sustained lunar operations. Since many applicable lander and upper stage vehicles utilize cryogenic hydrogen and oxygen propellants, it is appealing to incorporate fuel cells into potential mission concepts. A fuel cell could generate electricity, heat, and water from the residual hydrogen and oxygen propellants. Such a concept depends on the capability of the fuel cell to utilize dry, propellant-grade reactants at reduced pressures compared to typical fuel cell operating conditions. This report describes a NASA evaluation of two passive water removal, non-flow-through proton exchange membrane (PEM) fuel cells: a twelve-cell 50 sq. cm active area stack and a seven-cell 150 sq. cm active area stack. Both stacks operated on research grade hydrogen and oxygen reactants supplied at pressures ranging from 20 psia up to the 45-48 psia design operating point, reliably supporting current densities of at least 450 mA/sq. cm at reduced pressures. Polarization curves, fitted with a semi-empirical equation for overpotential, show an increase in effective ohmic resistance and decrease in reversible cell potential, hydrogen crossover current, and Tafel slope when operating at lower pressures, resulting in a reduction in cell voltage at all current densities. Fitted parameters are compared to published values and rationalized with first principles explanations.

fuel cells↗

Characterization of Collection Efficiency of the Common Research Model Midspan Wing Section in the IRT

This paper presents a preliminary study for the characterization of collection efficiency from icing tests conducted in the Icing Research Tunnel at NASA Glenn Research Center. A test method previously developed for measuring the attachment line maximum collection efficiency of a swept NACA 0012 airfoil model at zero angle of attack was applied to the leading-edge region of a 65%-scale version of the Common Research Model midspan wing section. A correlation for the stagnation line maximum collection efficiency as a function of the modified inertia parameter was obtained with LEWICE3D simulations utilizing a discrete number of drop diameters. It was then compared with the collection efficiency measurement data obtained in the IRT. For the experimental collection efficiency, two ice shape digitization procedures were utilized to extract 2-D chord-wise ice shape profiles, i.e., the Maximum Combined Cross Section or MCCS, and the Minimum Combined Cross Section or Min CCS, at selected span-wise locations from the 3-D scanned ice shapes. The preliminary result showed that the rime ice thickness method can be used to characterize the collection efficiency distribution, for conditions free of ice erosion, in the main ice shape region. However, a high-order statistical processing of the ice scan is needed for estimating the mean ice thickness in areas where feathers are prevalent. From the limited comparison of the experimental and LEWICE3D collection efficiency data in the CRM65 MS model leading edge area, it was shown that the attachment line maximum collection efficiency is reasonably estimated by the best curve-fit correlation in the range of modified inertia parameter tested. A tighter alignment control of the iced and cold clean model scans is needed to improve the comparison. Further evaluation of this correlation is recommended to assess its applicability for Common Research Model type swept wing icing scaling analysis.

Ice Scailing↗

Characterization of Collection Efficiency of the Common Research Model Midspan Wing Section in the IRT

This paper presents a preliminary study for the characterization of collection efficiency from icing tests conducted in the Icing Research Tunnel at NASA Glenn Research Center. A test method previously developed for measuring the attachment line maximum collection efficiency of a swept NACA 0012 airfoil model at zero angle of attack was applied to the leading-edge region of a 65%-scale version of the Common Research Model midspan wing section. A correlation for the stagnation line maximum collection efficiency as a function of the modified inertia parameter was obtained with LEWICE3D simulations utilizing a discrete number of drop diameters. It was then compared with the collection efficiency measurement data obtained in the IRT. For the experimental collection efficiency, two ice shape digitization procedures were utilized to extract 2-D chord-wise ice shape profiles, i.e., the Maximum Combined Cross Section or MCCS, and the Minimum Combined Cross Section or MinCCS, at selected span-wise locations from the 3-D scanned ice shapes. The preliminary result showed that the rime ice thickness method can be used to characterize the collection efficiency distribution, for conditions free of ice erosion, in the main ice shape region. However, a high-order statistical processing of the ice scan is needed for estimating the mean ice thickness in areas where feathers are prevalent. From the limited comparison of the experimental and LEWICE3D collection efficiency data in the CRM65 MS model leading edge area, it was shown that the attachment line maximum collection efficiency is reasonably estimated by the best curve-fit correlation in the range of modified inertia parameter tested. A tighter alignment control of the iced and cold clean model scans as well as the model angle of attack are needed to improve the comparison. Further evaluation of this correlation is recommended to assess its applicability for Common Research Model type swept wing icing scaling analysis.

Icing Scaling↗

Electric Aircraft Propulsion Power System Impedance Modelling Methodology

This paper describes techniques to be used to model megawatt Electrified Aircraft Propulsion (MW EAP) components, specifically motor/generators and inverters, DC-DC converters, and long power leads; resulting in impedance-based models which require no information from the component manufacturers. A custom curve fitting method is described, and an example power supply model is fit using measured impedance data. A method called subscale measurement, employing measurement equipment which has much lower power capability than the MW EAP system under test, is described, simulated, and used to generate an impedance transfer function model of a load. Finally, stability analysis using a Nyquist based approach is performed on the combined supply and load system, using the generated models.

Megawatt electrified aircraft propulsion↗

Data Processing of Miniaturized Laser Heterodyne Radiometer (mini-LHR) Ground Instrument Retrievals

The Miniaturized Laser Heterodyne Radiometer (mini-LHR) is a passive ground instrument that observes the mole fraction of carbon dioxide (CO) and methane (CH)in the atmospheric column by measuring their absorption of sunlight at 1.6 microns. A laser heterodyne radiometer is similar in design to the super heterodyne radio receiver that is well known by ham radio enthusiasts. While not previously a commercial technique, laser heterodyne radiometers have a history of measuring atmospheric trace gases that started in the 1960s. With the commercial availability of inexpensive, low-power, thumbnail-sized lasers developed for the telecommunications industry, it was possible to miniaturize this technique and ultimately commercialize it. The mini-LHR has been under development at NASA GSFC since 2009. During that time, in addition to signicant technical improvements, processing has also evolved and been streamlined. Here we present details of the processing approach for raw mini-LHR data to produce 30- and 60-minute data products of CH and CO column mole fractions. Processing occurs in two general stages: a python-based pre-processing of raw data, followed by ingestion into a Planetary Spectrum Generator (PSG) retrieval algorithm. Raw data processing involves removal of outliers, correcting for changes in air mass throughout the day, averaging scans, and ultimately converting averaged scans into transmittance vs. wavelength. The PSG retrieval simulates a spectra for the time/day/location of the scan with meteorological inputs from Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) data set and then perturbs the concentrations of CO and CH to obtain a t based on an iterative least-squares curve fitting procedure.

Giancarlo Roberto Zambrano↗

Trend Detection of Atmospheric Time Series: Incorporating Appropriate Uncertainty Estimates and Handling Extreme Events

This paper is aimed at atmospheric scientists without formal training in statistical theory. Its goal is to, 1) provide a critical review of the rationale for trend analysis of the time series typically encountered in the field of atmospheric chemistry; 2) describe a range of trend-detection methods; and 3) demonstrate effective means of conveying the results to a general audience. Trend detections in atmospheric chemical composition data are often challenged by a variety of sources of uncertainty, which often behave differently to other environmental phenomena such as temperature, precipitation rate, or stream flow, and may require specific methods depending on the science questions to be addressed. Some sources of uncertainty can be explicitly included in the model specification, such as autocorrelation and seasonality, but some inherent uncertainties are difficult to quantify, such as data heterogeneity and measurement uncertainty due to the combined effect of short- and long-term natural variability, instrumental stability, and aggregation of data from sparse sampling frequency. Failure to account for these uncertainties might result in an inappropriate inference of the trends and their estimation errors. On the other hand, the variation in extreme events might be interesting for different scientific questions, for example, the frequency of extremely high surface ozone events and their relevance to human health. In this study we aim to, 1) review trend detection methods for addressing different levels of data complexity in different chemical species; 2) demonstrate that the incorporation of scientifically interpretable covariates can outperform pure numerical curve fitting techniques in terms of uncertainty reduction and improved predictability; 3) illustrate the study of trends based on extreme quantiles that can provide insight beyond standard mean or median based trend estimates; and 4) present an advanced method of quantifying regional trends based on the inter-site correlations of multi-site data. All demonstrations are based on time series of observed trace gases relevant to atmospheric chemistry, but the methods can be applied to other environmental data sets.

Trace gas↗

Specific Power and Efficiency Projections of Electric Machines and Circuit Protection Exploration for Aircraft Applications

The purpose of this paper is to generate specific power and efficiency projections through the year 2050 for electric machines for aircraft applications. A general literature review was performed to identify the types of electric machines that are commonly used and which types have the biggest potential for future aircraft applications due to their high specific power and efficiency. A database with historical data was built to include parameters such as weight [kg], rated power [kW], specific power [kW/kg], RPM, efficiency, year, motor cooling type, application type and motor type to allow for trend identification and accurate projections. Once the data was gathered, multiple curve fits on the historical data were generated and extrapolated to produce the projections for specific power according to conservative, nominal and aggressive projection scenarios. A different process was followed for the efficiency projections due to the scattered nature of the data. A state of the art (SoA) value for efficiency was identified through literature review and was used to create the conservative, nominal and aggressive projections for the time frames of 2030, 2040, and 2050. The efficiency and the specific power projections of EMs for 2050 are 0.989 and 50kW/kg respectively. This paper will also be examining circuit protection as it is an additional component of electric powertrains.

electric machines↗

Statistical Process Control and Analysis on the Water Content Measurements in NASA Glenn’s Icing Research Tunnel

The Icing Research Tunnel at NASA Glenn follows the recommended practice for calibration outlined in SAE’s ARP5905. The calibration team has followed the schedule of a full calibration every five years with a check calibration done every six months following. The liquid water content of the IRT has maintained stability within the stated specifications of variation within +/- 10% of the curve fit equation generated from calibration data. Using past measurements and data trends, IRT characterization engineers wanted to develop methods for the ability to know when data were not within variation. Trends can be observed in the liquid water content measurement process by constructing statistical process control charts. This paper describes data processing procedures for the Multi-Element Sensor in the IRT, including collision efficiency corrections, canonical correlation analysis, process for rejection of data, and construction of control charts. Data are presented to display the control capability to meet defined liquid water content specifications of the IRT with the Multi-Element Sensor mounted in the center of the test section.

Statistical Process Control↗

Statistical Process Control and Analysis on the Water Content Measurements in NASA Glenn’s Icing Research Tunnel

The Icing Research Tunnel at NASA Glenn follows the recommended practice for calibration outlined in SAE’s ARP5905. The calibration team has followed the schedule of a full calibration every five years with a check calibration done every six months following. The liquid water content of the IRT has maintained stability within the stated specifications of variation within +/- 10% of the curve fit equation generated from calibration data. Using past measurements and data trends, IRT characterization engineers wanted to develop methods for the ability to know when data were not within variation. Trends can be observed in the liquid water content measurement process by constructing statistical process control charts. This paper describes data processing procedures for the Multi-Element Sensor in the IRT, including collision efficiency corrections, canonical correlation analysis, process for rejection of data, and construction of control charts. Data are presented to display the control capability to meet defined liquid water content specifications of the IRT with the Multi-Element Sensor mounted in the center of the test section.

Statistical Process Control↗

Characterization of Large Drop Velocity in the NASA Icing Research Tunnel

This paper presents experimental work conducted in the Icing Research Tunnel at NASA Glenn Research Center to characterize the velocity of large drops in the test section. Some icing spray clouds with large drops were generated with Mod1 nozzles at low nozzle air pressure of 2 to 4 psig for various tunnel air speeds. Drop diameters and drop velocities were measured via high-resolution imaging with a Particle Imaging Particle Tracking Velocimetry probe developed by Artium Technologies. The probe was mounted at four different locations aligned with the centerline of the test section from near the end of the contraction to the constant height test section part of the tunnel. CFD analyses were performed. It showed that the probe head geometry affects the local air flow in the measurement area of the probe between the prongs and in front of the probe. Initial analysis of the air velocity data during the test also indicated that the probe mounting stand has blockage effect on the local tunnel air velocity measurement by a pitot-static probe affixed on the mounting plate next to the stand. Those findings were later verified in the Icing Research Tunnel using a new pitot-static probe design with a linear motion system to measure the local tunnel air velocity with and without the probe. As a result, additional drop trajectory simulations were performed with the airflow moving towards the probe head with a Langmuir-D 7-bin distribution cloud containing large drops. The simulation results helped identify a critical drop-size threshold of 300 μm above which the velocities of larger drops are negligibly affected by the adverse pressure gradient generated by the probe head due to their large drop inertia. From the dimensional analysis of the drop velocity measurement data obtained, it showed that at the tunnel test section reference location a generalized empirical correlation was developed for the non-dimensional drop velocity as a function of the non-dimensional corrected drop diameter independent of the actual spraybar pressure settings and the tunnel air speeds. The generalized curve-fit correlation showed that the drop velocity was universally asymptotic to about 86 percent of the corresponding tunnel air speed at the test section reference location for the largest drop diameter captured by the probe. Further evaluation of this correlation is recommended to assess its applicability for Supercooled Large Drop icing scaling applications in the Icing Research Tunnel.

SLD↗

Comparison and Application of Battery Modeling Methods for Conceptual Electric Aircraft

Growing interest in design and optimization of electrified aircraft propulsion concepts prompts the need for accurate, flexible, and efficient methods to model battery systems. Presented in this paper are three battery modeling methods that have been used at NASA’s Glenn Research Center, each representing different mathematical or electrical approaches. Thévenin equivalent circuit, normalization, and curve-fitting methods are compared against battery cell test data for the X-57 Maxwell electric aircraft technology demonstrator. The methods are then applied in a simple multidisciplinary optimization context using NASA’s Six-Passenger Electric Quadrotor concept to determine their applicability and performance. The normalization method achieves the highest accuracy for steady and unsteady discharge rates with a voltage mean error percentage of 0.423% and 1.186%, respectively. Optimal quadrotor mission range between the models varies up to 0.5 nmi, identifying current battery modeling methods as a potentially significant contributor to mission analysis error. A set of relevant tools and techniques for conceptual battery modeling are identified in this paper, with conclusions made on the utility of each modeling approach for various design challenges.

power↗

Characterization of Large Drop Velocity in the NASA Icing Research Tunnel

This presentation presents experimental work conducted in the Icing Research Tunnel at NASA Glenn Research Center to characterize the velocity of large drops in the tunnel test section. Some icing spray clouds with large-sized drops were generated with Mod1 nozzles at low nozzle air pressure of 2 to 4 psig for various tunnel air speeds. Drop diameters and drop velocities were measured via high-resolution imaging with a Particle Imaging Particle Tracking Velocimetry probe developed by Artium Technologies. The probe was mounted at four different locations aligned with the centerline of the test section from near the end of the contraction to the constant height test section part of the tunnel. CFD analyses were performed. It showed that the probe head geometry affects the local air flow between the prongs and in front of the probe. Initial analysis of the air velocity data during the test also indicated that the probe mounting stand has blockage effect on the local tunnel air velocity measurement by a Pitot static probe affixed on the mounting plate. Those findings were later verified in the Icing Research Tunnel using a new Pitot probe design with a linear motion system to measure the local tunnel air velocity with and without the probe. As a result, additional drop trajectory simulations as airflow moving towards the probe head were run for a Langmuir-D 7-bin drop size distribution of a spray cloud with a nominally large value of medium volumetric diameter. The simulation results helped identify a critical drop-size threshold of 300 µm above which the velocities of larger drops are not affected by the adverse pressure gradient generated by the probe head due to their large drop inertia. From the dimensional analysis of the drop velocity measurement data obtained, it showed that at the tunnel test section reference location a generalized empirical correlation was developed for the non-dimensional drop velocity as a function of the non-dimensional corrected drop diameter independent of the actual spraybar pressure settings and the tunnel air speeds. The generalized curve-fit correlation showed that the drop velocity was universally asymptotic to about 86 percent of the corresponding tunnel air speed at the test section reference location for the largest drop diameter captured by the probe. Further evaluation of this correlation is recommended to assess its applicability for Supercooled Large Drop icing scaling applications in the Icing Research Tunnel.

SLD↗

Aerothermal Analysis of the Rocket Lab Venus Probe Heatshield

This document provides an overview of the aerothermodynamic analyses performed by the Aerothermodynamics Branch at NASA Langley Research Center for the Rocket Lab Venus Probe (RLVP). In addition to defining the baseline heating environment to the heatshield, this document pursues the experimental validation of key physical models at RLVP-relevant conditions. This experimental validation analysis, which captures the model form uncertainty, is used as one of two primary components of the margin assessment, where the other component is the parametric uncertainty. These model form (experimental) and parametric uncertainty components are used to construct a spatial and time varying margin for the heating to the RLVP heatshield. The margin is evaluated as the sum of the parametric and model form uncertainty components. The model form uncertainty is defined as the difference between the RLVP-relevant measurements and their simulations, using the upper limit uncertainty bounds for both the measurements and simulations in the comparisons. The differences in the dominant physics in the stagnation region and flank lead to the separate RLVP-relevant measurements for assessing the model form uncertainty in these two regions. These regions are addressed as follows: Stagnation Region Heating Environment: For the high-temperature stagnation-region, both the radiative heating and impact of blowing on convective heating are significant, while the impacts of turbulence and roughness are negligible. Coupled radiation and ablation LAURA/HARA solutions with ray-tracing provide the radiative heating over the entire vehicle, including the contributions from the Venus atmosphere and ablation species. Non-ablating LAURA simulations provide the convective heating. During the material-response computation typically used for TPS sizing, this non-ablating convective heating is corrected for the impact of ablation using the blowing correction. Coupled ablation LAURA simulations that capture finite-rate sur-face processes show that this blowing correction may be non-conservative over most of the heatshield. This non-conservatism is due to hydrogen recombination in the finite-rate surface model, which tends to increase the coupled ablation convective heating to near the non-ablating values, therefore making any reduction in the non-ablating value through the blowing correction non-conservative. This non-conservatism due to H catalysis is captured in the parametric component of the margin. The best available ground-test measurements that capture the impact of blowing on stagnation region convective heating, at RLVP-relevant conditions, indicate that the current blowing reduction model is non-conservative by up to 20% at RLVP-relevant blowing rates (the coupled blowing simulations were also non-conservative). Because of the relatively low velocity of the ground tests and the non-Venus atmospheric chemistry, these measurements do not capture the chemistry and therefore do not inform the uncertainty due to H catalysis. However, they do capture the fluid mechanics of blowing. The non-conservatism of the blowing correction implied by these measurements is covered by the model form component of the margin, which leads to total margin values over 50%. For the radiative heating, the shock-tube informed bias approach suggests a model form uncertainty of roughly 20%, while the parametric uncertainty analysis suggests values over 100%. The combined stagnation-point radiation margin of over 100% leads to peak margined radiative heating values of over300 W/cm2, which remains small relative to the peak margined convective heating of nearly 2000 W/cm2. Based on this analysis, at the stagnation point, the peak margined heat rate is 2203 W/cm2 and the margined total heat load is 31.5 kJ/cm2 for the current nominal trajectory. Flank Heating Environment: The forebody flank (and near-shoulder) heating environment is dominated by the impact of turbulence, roughness augmentation, and ablation on the convective heating. An extensive collection of ground test measurements with RLVP-relevant turbulence and roughness is studied to show that the maximum difference between the simulated and measured convective heating is 5%. However, with the exception of the Holden measurements from the 1980s, these measurements do not include roughness elements extending into the supersonic region of the boundary layer, which is likely to occur for RLVP (due to the 45 degree sphere-cone geometry). The interaction between the supersonic flow and roughness could cause convective heating augmentation fundamentally different than for locally subsonic flow. Although these Holden measurements are consistent with the other measurements considered, another path was pursued to assure that the rough-ness height extending into supersonic flow does not fundamentally change the roughness augmentation. This additional path was a computational effort to resolve the roughness elements in the CFD grid, so that the interaction be-tween the roughness elements and locally supersonic flow may be simulated in detail. This roughness-resolved CFD simulation is feasible because of the RLVP forebody TPS’s patterned roughness, which may be approximated analytically, and because of the axisymmetric nominal flow field, which allows a narrow surface region to be simulated and therefore make the computational expense feasible. These grid-resolved roughness simulations, which are performed at actual RLVP flight conditions, result in heating augmentation values that are below the design approach for roughness augmentation. This provides evidence that the design approach for RLVP roughness augmentation is sufficient. Based on this analysis, at this flank or near-shoulder location, the peak margined total heat rate is 2088 W/cm2and the margined total heat load is 26.0 kJ/cm2for the current nominal trajectory. Heat flux, shear, pressure and heat transfer coefficient at the RLVP stagnation point and near shoulder location are evaluated for the entire trajectory, and curved fit to a functional form of F=AρB∞UC∞. These simplified relationships for the nominal and margined aerothermal environments are referred to as aerothermal indicators, and presented at the end of this document.

Christopher O Johnston↗