Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat conduction modeling”

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 595 records · Page 33

Dual Heat Pulse, Dual Layer Thermal Protection System Sizing Analysis and Trade Studies for Human Mars Entry Descent and Landing

NASA has been recently updating design reference missions for the human exploration of Mars and evaluating the technology investments required to do so. The first of these started in January 2007 and developed the Mars Design Reference Architecture 5.0 (DRA5). As part of DRA5, Thermal Protection System (TPS) sizing analysis was performed on a mid L/D rigid aeroshell undergoing a dual heat pulse (aerocapture and atmospheric entry) trajectory. The DRA5 TPS subteam determined that using traditional monolithic ablator systems would be mass expensive. They proposed a new dual-layer TPS concept utilizing an ablator atop a low thermal conductivity insulative substrate to address the issue. Using existing thermal response models for an ablator and insulative tile, preliminary hand analysis of the dual layer concept at a few key heating points indicated that the concept showed potential to reduce TPS masses and warranted further study. In FY09, the followon Entry, Descent and Landing Systems Analysis (EDL-SA) project continued by focusing on Exploration-class cargo or crewed missions requiring 10 to 50 metric tons of landed payload. The TPS subteam advanced the preliminary dual-layer TPS analysis by developing a new process and updated TPS sizing code to rapidly evaluate mass-optimized, full body sizing for a dual layer TPS that is capable of dual heat pulse performance. This paper describes the process and presents the results of the EDL-SA FY09 dual-layer TPS analyses on the rigid mid L/D aeroshell. Additionally, several trade studies were conducted with the sizing code to evaluate the impact of various design factors, assumptions and margins.

McGuire, Mary Kathleen↗

Cosmic-ray heating of cooling flows - A critical analysis

It is shown that a combination of MHD wave-mediated cosmic ray heating and thermal conduction could balance cooling in intracluster media and substantially reduce the rate of inflow. The appropriate system of steady state equations is solved, including a new self-consistent formulation for the cosmic-ray diffusivity. Models which can produce substantial positive temperature gradients in static configurations are found when conduction is reduced by a factor of 10 or more. These models have too-flat thermal pressure profiles compared with observations. It is found that cosmic-ray heating is unlikely either to stabilize positive density perturbations against condensation or to contribute appreciably to the powering of the optical filaments.

Loewenstein, Michael↗

Estimation of surface heat and moisture fluxes over a prairie grassland. IV - Impact of satellite remote sensing of slow canopy variables on performance of a hybrid biosphere model

Numerical experiments are conducted using the Ex-BATS model of Crosson and Weng (1993), which is an adaptation the Dickinson (1983, 1984) and Dickinson et al. (1986) biosphere model BATS. The purpose of these experiments is the assessment of the Ex-BATS performance when using remotely sensed data for the estimation of three key canopy variables retrieved from NOAA-AVHRR measurements: the total surface albedo, the leaf area index (LAI), and the nondiurnally varying component of stomatal resistance, r(s). The results of the simulations, which cover the entire FIFE 1987 time period, show that the satellite retrievals of r(s) are only 20 to 30 percent less accurate than the idealized results of the control experiment. The performance of the model which used satellite retrieval of the surface albedo and LAI was essentially equivalent to the hypothetical version.

Crosson, William L.↗

Evaluation of a Thermodynamic Ice Crystal Icing Model Using Experimental Ice Accretion Data

This paper presents the evaluation of a thermodynamic ice crystal icing model, previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. It has been proposed that there are two types of distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. In the former, ice accretion occurs where a freeze fraction (0 to 1) of melted ice crystals freezes on a surface, along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized by having strong adhesion to the surface. In the latter, icing occurs from accumulated unmelted ice on a surface, where a melt fraction (0 to 1) dictates the amount of unmelted impinged ice. This melt-dominated icing is characterized by weakly bonded surface adhesion. The experimentally observed ice growth rates suggest that only a small fraction of the impinging ice remains on the surface, implying a mass loss mechanism such as splash, runback, bounce, or erosion. This mass loss parameter must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface. This loss parameter, however, along with the freeze and melt fraction, are the only experimental parameters that are currently not measured directly. Using reported icing growth rates from published ice crystal icing experiments, a methodology is proposed to determine these unknown parameters. This work takes reported ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) in 2016 and tests NASA collaborated on with the National Research Council (NRC) of Canada in 2012 that examined the fundamental physics of ice crystal icing. Those research efforts sought to generate icing conditions representative of those that occur inside a jet engine when ingesting ice crystals. This paper presents the fundamental equations of the thermodynamic model, the methodology used to determine the aforementioned unknown icing parameters, and results from model evaluation using experimental data. In addition, this paper builds on the previously proposed model by adding a transient conduction term to explain ice growth behavior at the onset of experimental tests that was observed to be different from steady-state ice growth that occurred later in the test run.With the addition of this energy term, this becomes a quasi-steady model. A key finding from this work suggests that mass loss fractions can exceed 0.90 for steady ice growth periods. In addition, due to conductive heat fluxes when using a warmer-than-freezing airfoil, lower mass loss fraction values were calculated during the initial transient period.

Ice-Crystal Ice Accretion↗

Heat Flow, Thermal Conductivity, and the Plausibility of the White Mars Hypothesis

Due to the low thermal conductivity of CO2 ice and clathrate vs. water ice, we find that liquid water reservoirs would not be confined to the deep subsurface as predicted by the controversial White Mars model, even assuming low global heat flow. Additional information is contained in the original extended abstract.

Urquhart, M. L.↗

Calculating the High-Latitude Ionospheric Electrodynamics Using A Machine Learning-Based Field-Aligned Current Model

We introduce a new framework called Machine Learning (ML) based Auroral Ionospheric electrodynamics Model (ML-AIM). ML-AIM solves a current continuity equation by utilizing the ML model of Field Aligned Currents of Kunduri et al. (2020, https://doi.org/10.1029/2020JA027908), the FAC-derived auroral conductance model of Robinson et al. (2020, https://doi.org/10.1029/2020JA028008), and the solar irradiance conductance model of Moen and Brekke (1993, https://doi.org/10.1029/92gl02109). The ML-AIM inputs are 60-min time histories of solar wind plasma, interplanetary magnetic fields (IMF), and geomagnetic indices, and its outputs are ionospheric electric potential, electric fields, Pedersen/Hall currents, and Joule Heating. We conduct two ML-AIM simulations for a weak geomagnetic activity interval on 14 May 2013 and a geomagnetic storm on 7–8 September 2017. ML-AIM produces physically accurate ionospheric potential patterns such as the two-cell convection pattern and the enhancement of electric potentials during active times. The cross polar cap potentials (ΦPC) from ML-AIM, the Weimer (2005, https://doi.org/10.1029/2004ja010884) model, and the Super Dual Auroral Radar Network (SuperDARN) data-assimilated potentials, are compared to the ones from 3204 polar crossings of the Defense Meteorological Satellite Program F17 satellite, showing better performance of ML-AIM than others. ML-AIM is unique and innovative because it predicts ionospheric responses to the time-varying solar wind and geomagnetic conditions, while the other traditional empirical models like Weimer (2005, https://doi.org/10.1029/2004ja010884) designed to provide a quasi-static ionospheric condition under quasi-steady solar wind/IMF conditions. Plans are underway to improve ML-AIM performance by including a fully ML network of models of aurora precipitation and ionospheric conductance, targeting its characterization of geomagnetically active times.

auroral electrodynamics↗

Rocket thrust chamber thermal barrier coatings

A research program was conducted to generate data and develop analytical techniques to predict the performance and reliability of ceramic thermal barrier coatings in high heat flux environments. A finite element model was used to analyze the thermomechanical behavior of coating systems in rocket thrust chambers. Candidate coating systems (using a copper substrate, NiCrAlY bond coat and ZrO2.8Y2O3 ceramic overcoat) were selected for detailed study based on photomicrographic evaluations of experimental test specimens. The effects of plasma spray application parameters on the material properties of these coatings were measured and the effects on coating performance evaluated using the finite element model. Coating design curves which define acceptable operating envelopes for seleted coating systems were constructed based on temperature and strain limitations. Spray gun power levels was found to have the most significant effect on coating structure. Three coating systems were selected for study using different power levels. Thermal conductivity, strain tolerance, density, and residual stress were measured for these coatings. Analyses indicated that extremely thin coatings ( 0.02 mm) are required to accommodate the high heat flux of a rocket thrust chamber and ensure structural integrity.

Batakis, A. P.↗

Combustion Response of AP Composite Propellants

The Cohen & Strand model for AP composite propellants is applied as boundary conditions, one for AP and one for binder, in solving the heat conduction equation in each to compute linear and non-linear combustion response properties for each and for the aggregate propellant. Iterations couple AP and binder through the quasi-steady flame processes. Illustrative results for linear response functions (pressure-coupled and velocity-coupled) are presented for a monomodal AP propellant varying AP size, pressure and crossflow speed, and examples of non-linear responses to arbitrary waveforms are shown. A quantitative comparison with response function data is limited to one well-characterized research formulation. Mechanistic implications are discussed, including recommendations for future modeling work.

Shusser, Michael↗

Characteristics of the Martian atmosphere surface layer

Researchers extend elements of various terrestrial boundary layer models to Mars in order to estimate sensible heat, latent heat, and momentum fluxes within the Martian atmospheric surface layer. To estimate the molecular viscosity and thermal conductivity of a CO2-H2O gas mixture under Martian conditions, parameterizations were developed. Parameterizations for specific heat and and binary diffusivity were also determined. The Prandtl and Schmidt numbers derived from these thermophysical properties were found to range from 0.78 - 1.0 and 0.47 - 0.70, respectively, for Mars. Brutsaert's model for sensible and latent heat transport within the interfacial sublayer for both aerodynamically smooth and rough airflow was experimentally tested under similar conditions, validating its application to Martian conditions. For the surface sublayer, the researchers modified the definition of the Monin-Obukhov length to properly account for the buoyancy forces arising from water vapor gradients in the Martian atmospheric boundary layer. This length scale was then utilized with similarity theory turbulent flux profiles with the same form as those used by Businger et al. and others. It was found that under most Martian conditions, the interfacial and surface sublayers offer roughly comparable resistance to sensible heat and water vapor transport and are thus both important in determining the associated fluxes.

Clow, G. D.↗

Thermal Analysis of a Carbon Fiber Rope Barrier for Use in the Reusable Solid Rocket Motor Nozzle Joint-2

This study provides for development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the "open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media "Darcy" relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations, Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Clayton, J. Louie↗

X-38 Experimental Aeroheating at Mach 10

This report provides an update of the hypersonic aerothermodynamic wind tunnel test program conducted at the NASA Langley Research Center in support of the X-38 program. Global surface heat transfer distributions were measured on 0.0177 and 0.0236 scale models of the proposed X-38 configuration at Mach 10 in air. The parametrics that were investigated primarily include freestream unit Reynolds numbers of 0.6 to 2.2 million per foot and body flap deflections of 15, 20, and 25 deg for an angle-of-attack of 40 deg. The model-scale variance was tested to obtain laminar, transitional, and turbulent heating levels on the defected bodyflaps. In addition, a limited investigation of forced boundary layer transition through the use of discrete roughness elements was performed. Comparisons of the present experimental results to computational predictions and previous experimental data were conducted Laminar, transitional, and turbulent heating levels were observed on the deflected body flap, which compared favorably to the computational results and to the predicted heating based on the flight aerothermodynamic database.

Berry, Scott A.↗

Thermal Analysis of a Carbon Fiber Rope Barrier for Use in the Reusable Solid Rocket Motor Nozzle Joint-2

This study provides development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the 'open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media 'Darcy' relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations. Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Clayton, J. Louie↗

Are Solar Activity Variations Amplified by the QBO: A Modeling Study

Solar cycle activity effects (SCAE) in the lower and middle atmosphere, reported in several studies, are difficult to explain on the basis of the small changes in solar radiation that accompany the 11-year cycle. It is therefore natural to speculate that dynamical processes may come into play to produce a leverage. Such a leverage may be provided by the Quasi-Biennial Oscillation (QBO) in the zonal circulation of the stratosphere, which has been linked to solar activity variations. Driven primarily by wave mean flow interaction, the QBO period and its amplitude are variable but are also strongly influenced by the seasonal cycle in the solar radiation. This influence extends to low altitudes and is referred to as 'downward control'. Small changes in the solar radiative forcing may produce small changes in the period and phase of the QBO, but these in turn may produce measurable differences in the wind field. Thus, the QBO may be an amplifier of solar activity variations and a natural conduit of these variations to lower altitudes. To test this hypothesis, we conducted experiments with a 2D version of our Numerical Spectral Model that incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). Solar cycle radiance variations (SCRV) are accounted for by changing the radiative heating rate on a logarithmic scale from 0.1% at the surface to 1% at 50 km to 10% at 100 km. With and without SCRV, but with the same GW flux, we then conduct numerical experiments to evaluate the magnitude of the SCAE in the zonal circulation. The numerical results indicate that, under certain conditions, the SCAE is significant and can extend to lower altitudes where the SCRV is small. For a modeled QBO period of 30 months, we find that the seasonal cycle in the solar forcing acts as a strong pacemaker to lock up the phase and period of the QBO. The SCAE then shows up primarily as a distinct but relatively weak amplitude modulation. But with a different QBO period between 30 and 34 (or less than 30, presumably) months, the seasonal phase lock is weak. Solar flux variations in the seasonal cycle then cause variations in the QBO period and phase. These amplify the SCAE to produce relatively large variations in the wind field. The SCAE in this case extends to mid-latitudes.

Mengel, J. G.↗

Metallurgical and Thermal Properties of 57Ni–40Ti–3Hf (wt%)

Selected properties of 57Ni-40Ti-3Hf (wt%) were measured to further develop this material for use in aerospace rolling-element bearings. The average grain size and microindentation hardness of the material at room temperature were 90 μm and 549 HV, respectively. From room temperature to 500 °C, electrical resistivity increased from 8.9×10 –7 to 11.3×10 –7 Ω⋅m, thermal expansion coefficients varied from 1.1×10 –5 to 1.2×10 –5 mm/mm⋅°C, specific heat increased from 0.44 to 0.53 W⋅sec/g⋅K, thermal diffusivity rose from approximately 0.003 to 0.005 cm 2 /sec and thermal conductivity increased from approximately 0.09 to 0.19 W/cm⋅K. The behavior of the material can be explained by well-known metallurgical theory. These results are expected to provide a basis for design and modeling of future aerospace components.

intermetallics↗

Thermal Management Tools for Propulsion System Trade Studies and Analysis

Energy-related subsystems in modern aircraft are more tightly coupled with less design margin. These subsystems include thermal management subsystems, vehicle electric power generation and distribution, aircraft engines, and flight control. Tighter coupling, lower design margins, and higher system complexity all make preliminary trade studies difficult. A suite of thermal management analysis tools has been developed to facilitate trade studies during preliminary design of air-vehicle propulsion systems. Simulink blocksets (from MathWorks) for developing quasi-steady-state and transient system models of aircraft thermal management systems and related energy systems have been developed. These blocksets extend the Simulink modeling environment in the thermal sciences and aircraft systems disciplines. The blocksets include blocks for modeling aircraft system heat loads, heat exchangers, pumps, reservoirs, fuel tanks, and other components at varying levels of model fidelity. The blocksets have been applied in a first-principles, physics-based modeling and simulation architecture for rapid prototyping of aircraft thermal management and related systems. They have been applied in representative modern aircraft thermal management system studies. The modeling and simulation architecture has also been used to conduct trade studies in a vehicle level model that incorporates coupling effects among the aircraft mission, engine cycle, fuel, and multi-phase heat-transfer materials.

McCarthy, Kevin↗

Porosity and the ecology of icy satellites

The case for a significant role for porosity in the structure and evolution of icy bodies in the Solar System has been difficult to establish. We present a relevant new data set and a series of structure models including a mechanical compression, not thermal creep, model for porosity that accounts satisfactorily for observed densities, moments of inertia, geologic activity, and sizes of tectonic features on icy satellites. Several types of observational data sets have been used to infer significant porosity, but until recently, alternative explanations have been preferred. Our first area of concern is the occurrence of cryovolcanism as a function of satellite radius; simple radiogenic heating models of icy satellites suggest minimum radii for melting and surface cryovolcanism to be 400 to 500 km, yet inferred melt deposits are seen on satellites half that size. One possible explanation is a deep, low conductivity regolith which lowers conductivity and raises internal temperatures, but other possibilities include tidal heating or crustal compositions of low conductivity. Our second area of concern is the occurrence and magnitude of tectonic strain; tectonic structures have been seen on icy satellites as small as Mimas and Proteus. The structures are almost exclusively extensional, with only a few possible compression Al features, and inferred global strains are on the order of 1 percent expansion. Expansions of this order in small bodies like Mimas and prevention of late compressional tectonics due to formation of ice mantles in larger bodies like Rhea are attained only in structure models including low-conductivity, and thus possibly high porosity, crusts. Thirdly, inferred moments of inertia less than 0.4 in Mimas and Tethys can be explained by high-porosity crusts, but also by differentiation of a high density core. Finally, the relatively low densities of smaller satellites like Mimas and Miranda relative to larger neighbors can be explained by deep porosity, but also by bulk compositional differences. Recent work has strengthened the case for significant porosity. Halley's nucleus was found to have a density near 0.6 g/cu cm, Janus and Epimethus were proposed to have densities near 0.7 g/cu cm, densities almost certainly due to high porosity. The irregular-spherical shape transition of icy satellites was quantitatively explained by low conductivity regoliths. A creative structure/thermal history model for Mimas simultaneously accounts quantitatively for Mimas' low density and moment of inertia by invoking initial high-porosity and subsequent compaction in the deep interior by thermal creep. The main problem with this promising model is that approximately 7 percent predicts a reduction in Mimas' radius, implying significant compressional failure and prevention of extensional tectonics, in contradiction to the observed extensional features and inferred 1 percent expansion in radius.

Croft, Steven K.↗

The Prediction of Broadband Shock-Associated Noise from Dualstream and Rectangular Jets Using RANS CFD

Supersonic jets operating off-design produce broadband shock-associated noise. Broadband shock-associated noise is characterized by multiple broadband peaks in the far-field and is often the dominant source of noise towards the sideline and upstream direction relative to the jet axis. It is due to large scale coherent turbulence structures in the jet shear layers interacting with the shock cell structure. A broadband shock-associated noise model recently developed by the authors predicts this noise component from solutions to the Reynolds averaged Navier-Stokes equations using a two-equation turbulence model. The broadband shock-associated noise model is applied to dualstream and rectangular nozzles operating supersonically, heated, and off-design. The dualstream jet broadband shock-associated noise predictions are conducted for cases when the core jet is supersonic and the fan jet is subsonic, the core jet is subsonic and the fan jet is supersonic, and when both jet streams operate supersonically. Rectangular jet predictions are shown for a convergent-divergent nozzle operating both over- and under-expanded for cold and heated conditions. The original model implementation has been heavily modified to make accurate predictions for the dualstream jets. It is also argued that for over-expanded jets the oblique shock wave attached to the nozzle lip contributes little to broadband shock-associated noise. All predictions are compared with experiments.

Miller, Steven A.↗