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 1,099 records · Page 61
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation and contact interfaces, and example simulations are included. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗
NASA NTRS2016

Flow Characterization Studies of the 10-MW TP3 Arc-Jet Facility: Probe Sweeps

This paper reports computational simulations and analysis in support of calibration and flow characterization tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted in the NASA Ames 10-MW TP3 facility using flat-faced stagnation calorimeters at six conditions corresponding to the steps of a simulated flight heating profile. Data were obtained using a conical nozzle test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot pressure and heat flux probes were also performed at these arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. Two different sets of pitot pressure and heat probes were used: 9.1-mm sphere-cone probes (nose radius of 4.57 mm or 0.18 in) with null-point heat flux gages, and 15.9-mm (0.625 in) diameter hemisphere probes with Gardon gages. The probe survey data clearly show that the test flow in the TP3 facility is not uniform at most conditions (not even axisymmetric at some conditions), and the extent of non-uniformity is highly dependent on various arc-jet parameters such as arc current, mass flow rate, and the amount of cold-gas injection at the arc-heater plenum. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested. Comparisons of computations with the experimental measurements show reasonably good agreement except at the extreme low pressure conditions of the facility envelope.

Goekcen, Tahir↗
NASA NTRS2016

Flow Characterization Studies of the 10-MW TP3 Arc-Jet Facility: Probe Sweeps

This paper reports computational simulations and analysis in support of calibration and flow characterization tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted in the NASA Ames 10-MW TP3 facility using flat-faced stagnation calorimeters at six conditions corresponding to the steps of a simulated flight heating profile. Data were obtained using a conical nozzle test configuration in which the models were placed in a free jet downstream of the nozzle. Experimental surveys of arc-jet test flow with pitot pressure and heat flux probes were also performed at these arc-heater conditions, providing assessment of the flow uniformity and valuable data for the flow characterization. Two different sets of pitot pressure and heat probes were used: 9.1-mm sphere-cone probes (nose radius of 4.57 mm or 0.18 in) with null-point heat flux gages, and 15.9-mm (0.625 in) diameter hemisphere probes with Gardon gages. The probe survey data clearly show that the test flow in the TP3 facility is not uniform at most conditions (not even axisymmetric at some conditions), and the extent of non-uniformity is highly dependent on various arc-jet parameters such as arc current, mass flow rate, and the amount of cold-gas injection at the arc-heater plenum. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested. Comparisons of computations with the experimental measurements show reasonably good agreement except at the extreme low pressure conditions of the facility envelope.

CFD↗
NASA NTRSDate not supplied

Modeling and Analysis of Stirling Power Convertors

Modeling and Analysis of Stirling Power Convertors Luis A. Rodriguez1 Steven M. Geng, Terry V. Reid, Scott D. Wilson NASA Glenn Research Center, Cleveland, OH, 44135, USA NASA Glenn’s Thermal Energy Conversion Branch is supporting the development of the next generation free-piston Stirling power convertors. American Superconductor (AMSC) and Sunpower Inc. are the two firms under contract to develop the Flexure Isotope Stirling Convertor (FISC) and the Sunpower Robust Stirling Convertor (SRSC), respectively. To comprehend and forecast convertor performance, Sage, ANSYS® Maxwell, and ANSYS® Fluent were used to model the Stirling thermodynamic cycle, alternator electromagnetics, and piston and displacer dynamics. I. Introduction Stirling convertors are being developed by NASA as a potential steady source of electrical power for NASA’s future scientific space missions. Currently, NASA Glenn Research Center has two corporations under contract, American Superconductor (AMSC) and Sunpower Inc., for the development of the next generation of free-piston Stirling convertors for dynamic radioisotope power systems. AMSC is developing the Flexure Isotope Stirling Convertor (FISC), which uses flexures to prevent side motion and rubbing of the piston. Similarly, Sunpower Inc, is developing the Sunpower Robust Stirling Convertor (SRSC). The SRSC uses gas bearings to prevent radial contact of the moving piston. As convertor development continues, it is increasingly important to understand and predict the interactions of components in the system, how they respond to one another, and how they perform as a response to changes in operating conditions. A suitable and enlightening way to demonstrate and foresee these interactions is with the use of accurate modeling software. Sage, ANSYS® Maxwell, and ANSYS® Fluent are the current modeling tools used by NASA to analytically determine convertor performance. Sage is a one-dimensional object-oriented commercial software package used for modeling and optimizing Stirling convertors for Dynamic Radioisotope Power Systems (DRPS) and it is one of the most accurate Stirling convertor codes in use by NASA. This code is the successor to GLIMPS (Globally-Implicit Stirling Cycle Simulation) and GLOP (GLIMPS Optimization) software created by Gedeon Associates [1]. Model input parameters are typically material/gas type, component physical dimensions, temperatures, frequency, charge pressure, and number of time/space nodes. Sage is used to model both the FISC’s and SRSC’s Stirling cycle thermodynamics and piston/displacer dynamics. Performance maps were created and analyzed for both power systems to better understand the relationship between the following conditions: cold-end temperature, hot-end temperature, piston/displacer amplitudes, pressure drop, and thermal input power. The synergy between these conditions will help determine parameter sensitivity. ANSYS® Maxwell was used to create a three-dimensional (3-D) axisymmetric model for both FISC and SRSC alternators. The significant physical components included in each model are the magnets, magnet carrier, outer/inner laminations, and the coil. Inputs to the model are piston amplitude, piston frequency, alternator load, coil resistance, tuning capacitance, and specific material properties. The alternator models calculate terminal voltage, current, piston/current phase, voltage/current phase, coil inductance, terminal power and efficiency. The RI2 losses, core (hysteresis and eddy) losses, and magnet/can eddy losses are also a part of the final results. ANSYS® Fluent is used to build 3-D computational fluid dynamic (CFD) models to examine the Stirling cycle thermodynamics for both the FISC and SRSC systems. Three-dimensional Computer Aided Design (CAD) models were used to create the physical components of each convertor. Steady-state simulations were conducted for hardware testing, prediction of environmental losses during testing, and generation of radiation look-up tables. The model inputs to the aforementioned analysis are the material properties and boundary thermal conditions. The steady-state model calculates temperature and heat flow distributions. Transient 3-D calculations were also part of the CFD analysis. In this study a physically reduced version of the FISC is used to obtain a prediction of available engine power. For the gas bearing SRSC, the transient effort is used to obtain a prediction of bearing pad performance and its sensitivity to micro-channel geometric variation. The model inputs to the transient simulations are the piston amplitude, displacer amplitude, frequency, displacer/piston phase angle, dynamic deforming CFD grid, temperature boundary conditions, and user defined files describing motion profile of piston/displacer. The results of the model are temperature distributions, heat distributions, and PV power produced at pre-determined conditions.

Luis A Rodriguez↗
NASA NTRSDate not supplied

Plume-Surface Interaction: Preliminary Observations from a Physics Focused Ground Test

Near surface operations conducted by spacecraft using rocket propulsion, such as during landing or the initial portion of ascent, may induce surface interactions that pose a risk to the spacecraft itself or nearby assets. NASA’s Space Technology Mission Directorate is conducting a multi-year project to mature the capability to predict plume-surface interactions (PSI) and reduce uncertainty through modeling, simulation, and ground testing. The Physics Focused Ground Test (PFGT), conducted in summer 2021, aimed to collect PSI data for plume, erosion, and ejecta physics to characterize PSI behaviors across a range of parameters relevant to the validation of computational modeling and with consideration to flight-relevant, though not flight-scale, environments. PFGT is a sub-scale, intrusive half-plane, inert-gas test conducted in a 15 foot-diameter vacuum chamber using a supersonic, heated, gaseous nitrogen plume. Tests were conducted with six regolith simulants, varying in complexity from spherical glass beads to BP-1 lunar soil simulant, and varied vacuum chamber ambient pressures to simulate Martian and lunar conditions. Nozzle height and mass flow rate were also varied to observe PSI behaviors and transitions of interest. Three high speed cameras captured crater formation and ejecta behavior during each test. An overview of this experiment is presented along with preliminary observations and analysis.

Wesley A Chambers↗
NASA NTRS1981

Application of Monte Carlo techniques to transient thermal modeling of cavity radiometers having diffuse-specular surfaces

A viable alternative to the net exchange method of radiative analysis which is equally applicable to diffuse and diffuse-specular enclosures is presented. It is particularly more advantageous to use than the net exchange method in the case of a transient thermal analysis involving conduction and storage of energy as well as radiative exchange. A new quantity, called the distribution factor is defined which replaces the angle factor and the configuration factor. Once obtained, the array of distribution factors for an ensemble of surface elements which define an enclosure permits the instantaneous net radiative heat fluxes to all of the surfaces to be computed directly in terms of the known surface temperatures at that instant. The formulation of the thermal model is described, as is the determination of distribution factors by application of a Monte Carlo analysis. The results show that when fewer than 10,000 packets are emitted, an unsatisfactory approximation for the distribution factors is obtained, but that 10,000 packets is sufficient.

Mahan, J. R.↗
NASA NTRS1993

Turbulence modeling of free shear layers for high performance aircraft

In many flowfield computations, accuracy of the turbulence model employed is frequently a limiting factor in the overall accuracy of the computation. This is particularly true for complex flowfields such as those around full aircraft configurations. Free shear layers such as wakes, impinging jets (in V/STOL applications), and mixing layers over cavities are often part of these flowfields. Although flowfields have been computed for full aircraft, the memory and CPU requirements for these computations are often excessive. Additional computer power is required for multidisciplinary computations such as coupled fluid dynamics and conduction heat transfer analysis. Massively parallel computers show promise in alleviating this situation, and the purpose of this effort was to adapt and optimize CFD codes to these new machines. The objective of this research effort was to compute the flowfield and heat transfer for a two-dimensional jet impinging normally on a cool plate. The results of this research effort were summarized in an AIAA paper titled 'Parallel Implementation of the k-epsilon Turbulence Model'. Appendix A contains the full paper.

Sondak, Douglas↗
NASA NTRS1990

Characteristics of the Martian atmosphere surface layer

Elements of various terrestrial boundary layer models are extended to Mars in order to estimate sensible heat, latent heat, and momentum fluxes within the Martian atmospheric surface ('constant flux') layer. The atmospheric surface layer consists of an interfacial sublayer immediately adjacent to the ground and an overlying fully turbulent surface sublayer where wind-shear production of turbulence dominates buoyancy production. Within the interfacial sublayer, sensible and latent heat are transported by non-steady molecular diffusion into small-scale eddies which intermittently burst through this zone. Both the thickness of the interfacial sublayer and the characteristics of the turbulent eddies penetrating through it depend on whether airflow is aerodynamically smooth or aerodynamically rough, as determined by the Roughness Reynold's number. Within the overlying surface sublayer, similarity theory can be used to express the mean vertical windspeed, temperature, and water vapor profiles in terms of a single parameter, the Monin-Obukhov stability parameter. To estimate the molecular viscosity and thermal conductivity of a CO2-H2O gas mixture under Martian conditions, parameterizations were developed using data from the TPRC Data Series and the first-order Chapman-Cowling expressions; the required collision integrals were approximated using the Lenard-Jones potential. Parameterizations for specific heat and binary diffusivity were also determined. The Brutsart 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 definition of the Monin-Obukhov length was modified to properly account for the buoyancy forces arising from water vapor gradients in the Martian atmospheric boundary layer. 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.↗
NASA NTRS1994

Al III, Si IV, and C IV absorption toward zeta Ophiuchi: Evidence for photionized and collisionally ionized gas

We present Goddard High-Resolution Spectrograph observations at 3.5 km/s resolution and signal-to-noise ratios of 30 to 60 for the Al III, Si IV, and N V absorption lines in the far-ultraviolet spectrum of the O9.5 V star zeat Ophiuchi. The measurement reveal three types of highly ionized gas along the 140 pc line of sight. (1) Narrow components of Al III (b = 4.3 km/s, the mean value of (v(helio)) = -7.8 km/s; b = 3.2 km/s, the mean value of (v(sub helio)) = -14.4 km/s) and Si IV (b = 5.3 km/s, the mean value of (v(sub helio)) = -15.0 km/s) trace photionized gas in the expanding H II region surrounding zeta Oph. The observed magnitude and direction of the velocity offset between the Al III and Si IV profiles can be explained by models of H II regions that incorporate expansion. Narrow C IV absorption associated with the H II region is not detected. Predictions of the expected amounts of Si IV and C IV overestimate the column densities of these ions by factors of 30 and more than 10, respectively. The discrepancy may be due to the effects of elemental depletions in the gas and/or to the interaction of the stellar wind with surrounding matter. (2) Broad (b = 15 to 18 km/s) and weak Si IV and C IV absorption components are detected near the mean value of (v(sub helio)) = -26 km/s. The high-ionization species associated with these absorption components are probably produced by electron collisional ionization in a heated gas. This absorption may be physically related to the zeta Oph bow shock ot to a cloud complex situated within the local interstellar medium at d less than 60 pc. The C IV to Si IV column density ratio in this gas is 8, a factor of 6 less than conductive interface models predict, but this discrepancy may be removed by considering the effects of self-photoionization within the cooling gas in the model calculations. (3) A broad (b = 13 km/s) and weak C IV absorption feature detected at the mean value of (v(sub helio)) = -61 km/s is not seen in other species. We tentatively ascribe this absorption to gas in a postshock ragion of an optically thin shock in the zeta Oph stellar wind.

Sembach, Kenneth R.↗
NASA NTRS1973

Thermal protection system gap heating rates of the Rockwell International flat plate heat transfer model (OH2A/OH2B)

Heat transfer data for the Rockwell International Flat Plate Thermocouple Model are presented. The model simulated the Space Shuttle Vehicle Thermal Protection System. Data were recorded for locations in and around various size gaps for various gap orientation configurations. The test was conducted at Mach 5.1 for free-stream Reynolds number per foot values from 500,000 to 1,500,000.

Foster, T. F.↗
NASA NTRS1982

Thermoelectric refrigeration for temperatures below 100 K: A study of titanium sesquioxide

Previous measurements of the specific heat of V-doped Ti2O3 at low temperatures were explained by a model which also suggested the material would have a high thermoelectric figure-of-merit. The sample preparation, experimental apparatus, and the results of measurements on the thermal conductivity, thermoelectric power, and electrical resistivity of a single crystal Ti2O3 - 4% V sample are described. The results are used to derive the thermoelectric figure-of-merit between 5 and 300 K. The figure-of-merit is much smaller than expected and of little practical value because of the very high phonon thermal conductivity.

Redebaugh, R.↗
NASA NTRS1988

Overview of numerical codes developed for predicted electrothermal deicing of aircraft blades

An overview of the deicing computer codes that have been developed at the University of Toledo under sponsorship of the NASA-Lewis Research Center is presented. These codes simulate the transient heat conduction and phase change occurring in an electrothermal deicier pad that has an arbitrary accreted ice shape on its surface. The codes are one-dimensional rectangular, two-dimensional rectangular, and two-dimensional with a coordinate transformation to model the true blade geometry. All modifications relating to the thermal physics of the deicing problem that have been incorporated into the codes will be discussed. Recent results of reformulating the codes using different numerical methods to increase program efficiency are described. In particular, this reformulation has enabled a more comprehensive two-dimensional code to run in much less CPU time than the original version. The code predictions are compared with experimental data obtained in the NASA-Lewis Icing Research Tunnel with a UH1H blade fitted with a B. F. Goodrich electrothermal deicer pad. Both continuous and cyclic heater firing cases are considered. The major objective in this comparison is to illustrate which codes give acceptable results in different regions of the airfoil for different heater firing sequences.

Keith, Theo G.↗
NASA NTRS1991

Experimental investigation of passive infrared ice detection for helicopter applications

A technique is proposed to remotely detect rotor icing on helicopters. Using passive infrared (IR) thermometry it is possible to detect the warming caused by latent heat released as supercooled water freezes. During icing, the ice accretion region on the blade leading edge will be warmer than the uniced trailing edge resulting in a chordwise temperature profile characteristic of icing. Preliminary tests were conducted on a static model in the NASA Icing Research Tunnel for a variety of wet (glaze) and dry (rime) ice conditions. The characteristic chordwise temperature profiles were observed with an IR thermal video system and confirmed with thermocouple measurements. A prototype detector system was built consisting of a single point IR pyrometer, and experiments were run on a small scale rotor model. Again the characteristic chordwise temperature profiles were observed during icing, and the IR system was able to remotely detect icing. Based on the static and subscale rotor tests the passive IR technique is promising for rotor ice detection.

Dershowitz, Adam↗
NASA NTRS1991

Experimental investigation of passive infrared ice detection for helicopter applications

A technique is proposed to remotely detect rotor icing on helicopters. Using passive infrared (IR) thermometry, it is possible to detect the warming caused by latent heat released as supercooled water freezes. During icing, the ice accretion region on the blade leading edge will be warmer than the uniced trailing edge, resulting in a chordwise temperature profile characteristic of icing. Preliminary tests were conducted on a static model in the NASA Icing Research Tunnel for a variety of wet (glaze) and dry (rime) ice conditions. The characteristic chordwise temperature profiles were observed with an IR thermal video system and confirmed with thermocouple measurements. A prototype detector system was built consisting of a single point IR pyrometer. Experiments were run on a small scale rotor model. Again, the characteristic chordwise temperature profiles were observed during icing, and the IR system was able to remotely detect icing. Based on the static and subscale rotor tests, the passive IR technique is promising for rotor ice detection.

Dershowitz, Adam↗
NASA NTRS1992

Analysis of surface roughness generation in aircraft ice accretion

Patterns of roughness evolution have been studied analysis of high magnification video observations of accreting ice surfaces provided by the NASA Lewis Research Center. Three distinct patterns of surface roughness generation have been identified within the parametric regions studied. They include: Rime, Multi-Zone Glaze, and Uniform Glaze. Under most icing conditions, a brief period of transient rime ice growth was observed caused by heat conduction into the body. The resulting thin rime layer explains previously observed insensitivity of some ice accretions to substrate insensitivity of some ice accretions to substrate surface chemistry and may provide justification for simplifying assumptions in ice accretion sailing and modeling effects.

Hansman, R. J., Jr.↗
NASA NTRS2003

Response Modeling of Lightweight Charring Ablators and Thermal Radiation Testing Results

Under NASA's In-Space Propulsion/Aerocapture Program, ARA conducted arc-jet and thermal-radiation ablation test series in 2003 for advanced development, characterization, and response modeling of SRAM-20, SRAM-17, SRAM-14, and PhenCarb-20 ablators. Testing was focused on the future Titan Explorer mission. Convective heating rates (CW) were as high as 153 W/sq cm in the IHF and radiation rates were 100 W/sq cm in the Solar Tower Facility. The ablators showed good performance in the radiation environment without spallation, which was initially a concern, but they also showed higher in-depth temperatures when compared to analytical predictions based on arc-jet thermal-ablation response models. More testing in 2003 is planned in both of these facility to generate a sufficient data base for Titan TPS engineering.

Congdon, William M.↗
NASA NTRS1996

Comparison of Experimental and Computational Aerothermodynamics of a 70-deg Sphere-Cone

Numerical solutions for hypersonic flows of carbon-dioxide and air around a 70-deg sphere-cone have been computed using an axisymmetric non-equilibrium Navier-Stokes solver. Freestream flow conditions for these computations were equivalent to those obtained in an experimental blunt-body heat-transfer study conducted in a high-enthalpy, hypervelocity expansion tube. Comparisons have been made between the computed and measured surface heat-transfer rates on the forebody and afterbody of the sphere-cone and on the sting which supported the test model. Computed forebody heating rates were within the estimated experimental uncertainties of 10% on the forebody and 15% in the wake except for within the recirculating flow region of the wake.

Hollis, Brian R.↗
NASA NTRS2007

Testing of the Crew Exploration Vehicle in NASA Langley's Unitary Plan Wind Tunnel

As part of a strategic, multi-facility test program, subscale testing of NASA s Crew Exploration Vehicle was conducted in both legs of NASA Langley s Unitary Plan Wind Tunnel. The objectives of these tests were to generate aerodynamic and surface pressure data over a range of supersonic Mach numbers and reentry angles of attack for experimental and computational validation and aerodynamic database development. To provide initial information on boundary layer transition at supersonic test conditions, transition studies were conducted using temperature sensitive paint and infrared thermography optical techniques. To support implementation of these optical diagnostics in the Unitary Wind Tunnel, the experiment was first modeled using the Virtual Diagnostics Interface software. For reentry orientations of 140 to 170 degrees (heat shield forward), windward surface flow was entirely laminar for freestream unit Reynolds numbers equal to or less than 3 million per foot. Optical techniques showed qualitative evidence of forced transition on the windward heat shield with application of both distributed grit and discreet trip dots. Longitudinal static force and moment data showed the largest differences with Mach number and angle of attack variations. Differences associated with Reynolds number variation and/or laminar versus turbulent flow on the heat shield were very small. Static surface pressure data supported the aforementioned trends with Mach number, Reynolds number, and angle of attack.

Murphy, Kelly J.↗