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 919 records · Page 51

Axisymmetric spherical shell models of mantle convection with variable properties and free and rigid lids

Axisymmetric spherical shell numerical simulations of mantle convection were carried out to investigate the influence of two end-member surface stress conditions: stress-free and rigid. These correspond approximately to a subducting or a rigid lithosphere and can be seen as end-member models of the surface of Venus. Our model assumed an effective Rayleigh number of 3 x 10 exp 6, similar to that for earth, and included uniform internal heating and depth-dependent thermal expansivity and thermal conductivity. The simulations utilized a Newtonian viscosity which was constant or varied with depth and/or temperature. We show how the temperature, speed, and vorticity fields change qualitatively and quantitatively with surface temperature, surface stress condition, internal heating and viscosity distribution. We find that a rigid lid and viscosity which increases with depth both promote steady large-scale circulation with smaller-scale circulation in the upper mantle.

Leitch, A. M.↗

Aerothermodynamic Testing and Boundary Layer Trip Sizing of the HIFiRE Flight 1 Vehicle

An experimental wind tunnel test was conducted in the NASA Langley Research Center s 20-Inch Mach 6 Air Tunnel in support of the Hypersonic International Flight Research Experimentation Program. The information in this report is focused on the Flight 1 configuration, the first in a series of flight experiments. This report documents experimental measurements made over a range of Reynolds numbers and angles of attack on several scaled ceramic heat transfer models of the Flight 1 payload. Global heat transfer was measured using phosphor thermography and the resulting images and heat transfer distributions were used to infer the state of the boundary layer on the vehicle windside and leeside surfaces. Boundary layer trips were used to force the boundary layer turbulent, and a brief study was conducted to determine the effectiveness of the trips with various heights. The experimental data highlighted in this test report were used to size and place the boundary layer trip for the flight test vehicle.

Berger, Karen T.↗

Plans for Ice Crystal Icing Tests Using a 3D Heated Test Article at the NASA Icing Research Tunnel

This presentation describes the plans for ice crystal icing and supercooled water icing tests to be conducted in Fall/Winter 2021 at the NASA Icing Research Tunnel. A test article, whose geometry is representative of an inter-compressor duct and strut region of a turbofan engine, has been designed and fabricated for the study of ice crystal icing physics tests. The Simulated Inter-compressor Duct Research Model (SIDRM) will have the ability to heat surfaces to simulate the warm surfaces of the inter-compressor duct within an engine. The test article is instrumented with heaters, heat flux gauges, thermocouples, and pressure taps to measure icing behavior. The aim of these tests is to generate ice shapes on the SIDRM test article under well-characterized conditions. The results of the testing with the SIDRM model will be used to provide icing physics and validation data for development of the engine icing simulation capabilities in GlennICE, the 3D ice prediction tool. This is intended to be a presentation only, and does not have an accompanying paper.

ice-crystal icing↗

Aerothermodynamic Characteristics of Boundary Layer Transition and Trip Effectiveness of the HIFiRE Flight 5 Vehicle

An experimental wind tunnel test was conducted in the NASA Langley Research Center s 20-Inch Mach 6 Tunnel in support of the Hypersonic International Flight Research Experimentation Program. The information in this report is focused on the Flight 5 configuration, one in a series of flight experiments. This report documents experimental measurements made over a range of Reynolds numbers and angles of attack on several scaled ceramic heat transfer models of the Flight 5 vehicle. The heat transfer rate was measured using global phosphor thermography and the resulting images and heat transfer rate distributions were used to infer the state of the boundary layer on the windside, leeside and side surfaces. Boundary layer trips were used to force the boundary layer turbulent, and a study was conducted to determine the effectiveness of the trips with various heights. The experimental data highlighted in this test report were used determine the allowable roughness height for both the windside and side surfaces of the vehicle as well as provide for future tunnel-to-tunnel comparisons.

Berger, Karen T.↗

Lunar Surface Stirling Power Systems Using Isotope Heat Sources

For many years, NASA has used the decay of plutonium-238 (Pu-238) (in the form of the General Purpose Heat Source (GPHS)) as a heat source for Radioisotope Thermoelectric Generators (RTGs), which have provided electrical power for many NASA missions. While RTGs have an impressive reliability record for the missions in which they have been used, their relatively low thermal to electric conversion efficiency and the scarcity of plutonium-238 (Pu-238) has led NASA to consider other power conversion technologies. NASA is considering returning both robotic and human missions to the lunar surface and, because of the long lunar nights (14.75 Earth days), isotope power systems are an attractive candidate to generate electrical power. NASA is currently developing the Advanced Stirling Radioisotope Generator (ASRG) as a candidate higher efficiency power system that produces greater than 160 W with two GPHS modules at the beginning of life (BOL) (~32% efficiency). The ASRG uses the same Pu-238 GPHS modules, which are used in RTG, but by coupling them to a Stirling convertor provides a four-fold reduction in the number of GPHS modules. This study considers the use of americium-241 (Am-241) as a substitute for the Pu-238 in Stirling- convertor-based Radioisotope Power Systems (RPS) for power levels from tens of watts to 5 kWe. The Am-241 is used as a substitute for the Pu-238 in GPHS modules. Depending on power level, different Stirling heat input and removal systems are modeled. It was found that substituting Am-241 GPHS modules into the ASRG reduces power output by about one-fifth while maintaining approximately the same system mass. In order to obtain the nominal 160 W of electrical output of the Pu-238 ASRG requires 10 Am-241 GPHS modules. Higher power systems require changing from conductive coupling heat input and removal from the Stirling convertor to either pumped loops or heat pipes. Liquid metal pumped loops are considered as the primary heat transportation on the hot end and water pumped loop/heat pipe radiator is considered for the heat rejection side for power levels above 1 kWe.

Schmitz, Paul C.↗

Material Characterization and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

TPS↗

Material Properties and Modeling of Room Temperature Vulcanizing Silicone

Room Temperature Vulcanizing silicone (RTV) is a high-temperature adhesive that has successfully been used as a gap-filler between Thermal Protection System (TPS) tiles for heatshields on numerous missions. It is also used to bond instrumentation plugs such as temperature and pressure sensors into the heatshields. While RTV has been traditionally assumed to be a non-porous and non-ablating material, numerous experiments have shown that RTV pyrolyzes and becomes highly porous as it is heated. Heating RTV has also shown swelling, or intumescence, which can pose unique problems that lead to roughness induced boundary-layer transition, surface oxide formation and contamination of heat shield sensors. Therefore, it is crucial to understand and model the intumescence phenomenon of RTV. As data for RTV material properties is limited, the first step in modeling RTV is to collect material properties such as pyrolysis mass-loss, microstructure change, virgin and char porosity, etc. which was performed in our initial study. Additionally, thermomechanical properties such as Young’s modulus and Poisson ratio are required for modeling the intumescence of RTV, which were taken from literature and the coefficient of thermal expansion was collected using in-situ heating and Micro Computed Tomography (µ-CT) in previous studies. Finally, numerous other properties such as pyrolysis gas properties, virgin and char thermal conductivity and specific heat were compiled from previous experiments and literature into a material database that can be used for simulations. In Porous Material Analysis Toolbox based on OpenFOAM (PATO) [4], structural mechanics coupled with material response was used for simulating the intumescence of RTV as it is heated. However, since the permeability of the material is very low, the pyrolysis gas creates an internal pressure build-up as the material is being heated, significantly contributing to the deformation of the material. To correctly characterize this phenomenon, additional physics models were implemented into PATO's stress analysis solver, and results were compared with RTV dilatometry test data as a preliminary verification case. Future work will include experiments of RTV at the Plasmatron X facility and the in-situ heating cell with µ-CT, and improvement of simulation tools to more accurately model RTV intumescence.

PATO↗

Local Aerodynamic Heat Transfer and Boundary-Layer Transition on Roughened Sphere-Ellipsoid Bodies at Mach Number 3.0

A wind-tunnel investigation was made to determine heat-transfer distributions on three steel sphere-ellipsoid bodies with surface roughnesses of 5, 100, and 200 microinches. Tests were conducted in the Langley 9- by 6-foot thermal structures tunnel at a Mach number of 3.0, free-stream Reynolds numbers (based on model spherical diameter) of 4.25 x 10(exp 6) and 2.76 x l0(exp 6), and at a stagnation temperature of 650 F. Pressure distributions were obtained also on a fourth model. The results indicated that the combination of surface roughness and boundary-layer cooling tended to promote early transition and nullify the advantages attributable to the blunt shape of the model for reducing local temperatures. Good correlation between experimental heating rates and those calculated from laminar theory was achieved up to the start of boundary-layer transition. The correlation also was good with the values predicted by turbulent theory for surface stations downstream from the 45 deg. station.

Deveikis, William D.↗

Icing Physics Studies Using the 3D SIDRM Test Article: Aerodynamic and Supercooled Liquid Icing Analysis

In-flight icing is an important safety issue and is a factor that affects aircraft design and performance. Newer regulations are driving a need for improvements in airframe and engine icing simulation capability. Experimental data is required for development of icing physics models and simulation validation. To that end, this paper presents the analysis of the supercooled liquid icing data subset from tests conducted in 2022 at the NASA Icing Research Tunnel that studied both supercooled water and ice-crystal icing. The test article that was utilized replicated 3D geometrical features of an inter-compressor duct and strut region of a turbofan engine. The surfaces of the Simulated Inter-compressor Duct Research Model (SIDRM) can be heated to simulate the warm surfaces of the turbofan inter-compressor duct. The test article is instrumented with pressure taps, heaters, heat flux gauges, and thermocouples, while a 3D laser scanner, cameras, and a scale to measure ice mass were utilized to characterize the icing behavior. The aim of these tests was to generate ice accretions on the SIDRM test article under well-characterized supercooled liquid icing and ice crystal icing conditions. This paper discusses measurements related to aerodynamic testing and supercooled liquid icing tests that were conducted. Aerodynamic measurements were analyzed and compared to computational simulations and were found to be in good agreement for the range of airspeeds (50 to 230 knots) and angles of attack (0 to 4°) tested. Various parametric sweeps were conducted during the supercooled liquid icing portion of the test entry (cloud median volumetric diameter ranged from 15 to 90 µm, total air temperature from -3 to -17 °C, angle of attack from 0 to 4°, and accretion time from 5 to 20 min). These sweeps were performed to measure that parameter’s impact on ice accretion size, location (icing extent), characteristics (such as glaze/rime ice and shedding behavior), and test article surface temperature. Analysis of the test data showed that clouds composed of larger drops, colder air temperatures, smaller angles of attack, and longer spray times were the primary parameters that resulted in accretions with greater ice mass. Test article angle of attack and cloud droplet size influenced the location of ice accretion as these two parameters directly impact collection efficiency. With respect to icing characteristics, total air temperature dictated icing type, and smaller cloud drop size along with warmer air temperatures resulted in greater amounts of ice shedding. Surface temperature increased during ice accretion from the release of latent (fusion) heat, where total air temperature and cloud drop size impacted the amount of surface temperature change. The icing measurements collected during the SIDRM tests will be used to develop and validate 3D computational engine icing tools, such as GlennICE, that predictively assesses the onset and growth of ice. One of the goals of the sponsoring NASA project is to develop simulation models and tools that can assist in the design and certification of engines for flight in icing conditions in a cost effective way.

supercooled liquid icing↗

Zero-gravity transient thermal mixing simulation

The experimental program described is an outgrowth of independent investigations into alternate redesign concepts for the Apollo SM cryogenic oxygen storage system. The experiments were continued, after the redesign was established, to provide physical insight into transient thermal mixing in zero-gravity and to aid in the characterization of the system performance in flight. Zero-gravity heat transfer and fluid mixing were simulated experimentally through an analogy between unsteady heat conduction and species diffusion. To further support numerical analyses of the cryogenic oxygen storage system, the experimental investigation was extended to include a cubical tank geometry, representative of existing numerical models. In general, the transient flow patterns in the cubical tank are far more complex than those of the spherical tank and the extent of fluid mixing is significantly greater but less repeatable.

Lands, J. F., Jr.↗

The polarization and directivity of solar-flare hard X-ray bremsstrahlung from a thermal source

The polarization and directivity of hard X-ray bremsstrahlung from a thermal source consisting of a region in which a thermal flux drives a pair of steel collisionless conduction fronts were evaluated. The conduction fronts are symmetrically driven from a central region, heated by the flare energy dissipation process. By comparing results with similar calculations based on a nonthermal thick-target electron beam model of the source, it was aimed to determine the degree to which the observed polarization and directivity of solar flare hard X-rays favor either model. Results which exhibit significant polarization and directivity of the hard X-ray radiation emitted by the source are produced by using Maxwellian electron-phase-space distribution functions modified to take into account a directional heat flux and a steady direct current in the X-ray source, and a fully relativistic treatment of the bremsstrahlung emission process. The results are consistent with solar hard X-ray anisotropy and polarization observations to date, although these observations are too crude to be conclusive.

Emslie, A. G.↗

Carbon-fiber technology

The state of the art of PAN based carbon fiber manufacture and the science of fiber behavior is surveyed. A review is given of the stabilization by oxidation and the subsequent carbonization of fibers, of the apparent structure of fibers deduced from scanning electron microscopy, from X-ray scattering, and from similarities with soft carbons, and of the known relations between fiber properties and heat treatment temperature. A simplified model is invoked to explain the electrical properties of fibers and recent quantum chemical calculations on atomic clusters are used to elucidate some aspects of fiber conductivity. Some effects of intercalation and oxidative modification of finished fibers are summarized.

Hansen, C. F.↗

An experimental study of massive blowing from a nosetip during Jovian entry

The fluid dynamic structure and stability of the boundary layer and shock layer ahead of a highly blowing heat shield in an environment simulating entry into the Jovian atmosphere were investigated. The experimental studies were conducted in the Calspan 96 in shock tunnel with a freestream of a H2/He mixture at Mach numbers above 10. Measurements of heat transfer and pressure were made on the porous model surface, and high-speed schlieren photography was used to examine the structure and stability of the shock layer. Results are presented which show how the heat transfer and shock layer stability vary for gases of different molecular weights for a range of Reynolds numbers and mass-addition rates.

Holden, M. S.↗

A frost formation model and its validation under various experimental conditions

A numerical model that was used to calculate the frost properties for all regimes of frost growth is described. In the first regime of frost growth, the initial frost density and thickness was modeled from the theories of crystal growth. The 'frost point' temperature was modeled as a linear interpolation between the dew point temperature and the fog point temperature, based upon the nucleating capability of the particular condensing surfaces. For a second regime of frost growth, the diffusion model was adopted with the following enhancements: the generalized correlation of the water frost thermal conductivity was applied to practically all water frost layers being careful to ensure that the calculated heat and mass transfer coefficients agreed with experimental measurements of the same coefficients.

Dietenberger, M. A.↗

Parametric Weight Comparison of Advanced Metallic, Ceramic Tile, and Ceramic Blanket Thermal Protection Systems

A parametric weight assessment of advanced metallic panel, ceramic blanket, and ceramic tile thermal protection systems (TPS) was conducted using an implicit, one-dimensional (I-D) finite element sizing code. This sizing code contained models to account for coatings fasteners, adhesives, and strain isolation pads. Atmospheric entry heating profiles for two vehicles, the Access to Space (ATS) vehicle and a proposed Reusable Launch Vehicle (RLV), were used to ensure that the trends were not unique to a certain trajectory. Ten TPS concepts were compared for a range of applied heat loads and substructural heat capacities to identify general trends. This study found the blanket TPS concepts have the lightest weights over the majority of their applicable ranges, and current technology ceramic tiles and metallic TPS concepts have similar weights. A proposed, state-of-the-art metallic system which uses a higher temperature alloy and efficient multilayer insulation was predicted to be significantly lighter than the ceramic tile stems and approaches blanket TPS weights for higher integrated heat loads.

Myers, David E.↗

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work is limited to very simple models. In this project, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. The proposed work is based on previously funded NASA work that allowed us to develop a two-dimensional numerical method for boiling flows and to demonstrate the ability of the method to simulate film boiling. While numerical simulations of multi-fluid flows have been advanced in a major way during the last five years, or so, similar capability for flows with phase change are still in their infancy. Although the feasibility of the proposed approach has been demonstrated, it has yet to be extended and applied to fully three-dimensional simulations. Here, a fully three-dimensional, parallel, grid adaptive code will be developed. The numerical method will be used to study nucleate boiling in microgravity, with particular emphasis on two aspects of the problem: 1) Examination of the growth of bubbles at a wall nucleation site and the instabilities of rapidly growing bubbles. Particular emphasis will be put on accurately capturing the thin wall layer left behind as a bubble expands along a wall, on computing instabilities on bubble surfaces as bubbles grow, and on quantifying the effects of both these phenomena on heat transfer; and 2) Examination of the effect of shear flow on bubble growth and heat transfer.

Tryggvason, Gretar↗

The Material Response of PICA-NuSil to a Hypersonic Flow Environment

Phenolic Impregnated Carbon Ablator (PICA) is NASA’s baseline thermal protection system for missions to Mars. PICA is friable; therefore, NuSil ® CV-1144-0 is applied to the surface of flight hardware to mitigate the contamination of spacecraft sensors. NuSil is composed of a polysiloxane resin system that can transform into an oxidation-resistant coating upon heating. Therefore, a pathfinder campaign was conducted at the Hypersonic Materials Environmental Test System (HyMETS) at the NASA Langley Research Center to investigate the material response of instrumented sphere-cone models made from PICA coated with NuSil ® . Analysis of the post-test data is presented and suggests that NuSil ® impacts the material response of test articles subjected to oxidizing atmospheres, as evidenced by a reduction of the measured surface temperature and in-depth temperature response. Finally, a four-step ablation mechanism is proposed based on observations from high-speed video, optical emission spectroscopy, thermal analysis, and post-test surface characterization.

PICA PICA-NuSil Siloxane Silicone↗

Correlation and analysis of elevon seal test results

The heating rate and pressure distributions in and around the elevon coves were correlated. An analytical approach incorporating such a correlation into a thermal structure computer code to examine the elevon seal and cover structure in flight was developed. Emphasis was placed on the design of a thermal model that predicts gas and wall temperatures in the cove area. The elevon cove area was modelled analytically by a parallel plate flow channel, incorporating effects of axial wall conduction, axial and environmental radiant energy exchanges, and effects of laminar and turbulent convection entrance regions. Comparisons were made between computed and measured values of (1) cove wall temperatures; (2) local gas temperature within the cove; and (3) local heating rates.

Keshock, E. G.↗