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 757 records · Page 42

Computational and Experimental Study of Film-Cooling Effectiveness With and Without Downstream Vortex Generators

Downstream vortex generators that involve a pair of rectangular plates arranged in an open V-shape placed just downstream of each film-cooling hole were shown to create flow and vortical structures that entrain lifted film-cooling flow back to the surface and increase its lateral spreading on the surface (GT2020–14317). In this study, computations and measurements were performed to examine the flow mechanisms this vortex generator induces to improve film-cooling effectiveness of a flat plate with the cooling flow emanating from one row of inclined holes. Parameters studied include blowing ratio (BR = 0.75 and 1.0) and temperature ratio (TR = 1.07 and 1.9). Here, the computational study is based on steady Reynolds-averaged Navier–Stokes (RANS) closed by the shear-stress transport (SST) turbulence model with and without conjugate analysis. The experimental study was conducted by using a conjugate heat transfer test rig with a plenum, where cooling flow is introduced. Measurements made include velocity and temperature profiles upstream and downstream of the film-cooling holes as well as the temperature at several locations on the hot and cold sides of the film-cooled flat plate. The computational study was validated by comparing computed results with those from measurements at BR = 0.75 and 1.0 and TR = 1.9. Computational and experimental results are presented to show the effects of BR and TR on the flow structures and how those structures improve the effectiveness of film cooling with and without the downstream vortex generators and with and without conjugate heat transfer.

42 ENGINEERING↗

Oakridge Geothermal Resources Assessment

Through the Communities Local Energy Action Program (LEAP), NREL provided technical assistance to a coalition of stakeholders from Oakridge, Oregon. The coalition included community-based non-profit organizations, the city government, the local utility, and others. Technical assistance provides analysis and information to support Oakridge stakeholders with their goals to increase energy reliability and resilience in the community, while promoting economic development. The Geothermal Resource Assessment provides technical analysis related to the geothermal resource in the Oakridge area, possible technology applications, and information related to further geothermal exploration. Technical analysis includes a review of a previous geothermal study conducted for Oakridge, thermal modeling to estimate reservoir temperatures, and technoeconomic simulation of a combined heat and power application in the Oakridge area.

15 GEOTHERMAL ENERGY↗

Apparatus for determining thermophysical properties of test specimens

Apparatus is described for directly measuring the quantity square root of pck of a test specimen such as a wind tunnel model where p is density, c is the specific heat and k is the thermal conductivity of the specimen. The test specimen and a reference specimen are simultaneously subjected to the heat from a heat source. A thermocouple is attached to the reference specimen for producing a first electrical analog signal proportional to the heat rate Q that the test specimen is subjected to and an infrared radiometer that is aimed at the test specimen produces a second electrical analog signal proportional to the surface temperature T of the test specimen. An analog-to-digital converter converts the first and second electrical analog signals to digital signals. These digital signals are applied to a computer for determining the quantity.

Creel, T. R., Jr.↗

Theoretical aspects of tidal and planetary wave propagation at thermospheric heights

A simple semiquantitative model is presented which allows analytic solutions of tidal and planetary wave propagation at thermospheric heights. This model is based on perturbation approximation and mode separation. The effects of viscosity and heat conduction are parameterized by Rayleigh friction and Newtonian cooling. Because of this simplicity, one gains a clear physical insight into basic features of atmospheric wave propagation. In particular, we discuss the meridional structures of pressure and horizontal wind (the solutions of Laplace's equation) and their modification due to dissipative effects at thermospheric heights. Furthermore, we solve the equations governing the height structure of the wave modes and arrive at a very simple asymptotic solution valid in the upper part of the thermosphere. That 'system transfer function' of the thermosphere allows one to estimate immediately the reaction of the thermospheric wave mode parameters such as pressure, temperature, and winds to an external heat source of arbitrary temporal and spatial distribution. Finally, the diffusion effects of the minor constituents due to the global wind circulation are discussed, and some results of numerical calculations are presented.

Volland, H.↗

A Model of the Soft X-ray Background as a Blast Wave Viewed from Inside

The suggestion that the soft X-ray background arises in part from the Sun which is inside a large supernova blastwave was examined by models of spherical blastwaves. The models can produce quantitative fits to both surface brightnesses and energy band ratios when t = 10 to the 5th power E sub o = 5 x 10 to the 50th power ergs, and n sub approx. 0.004 cm to the -3 power. The models are generalized by varying the relative importance of factors such as thermal conduction, Coulomb heating of electrons, and external pressure; and to allow the explosions to occur in preexisting cavities with steep density gradients, or by examination of the effects of large obstructions or other anisotrophies in the ambient medium.

Edgar, R. J.↗

The loss and depth of CO2 ice in comet nuclei

An analytical model was developed to simulate the material differentiation of a cometary nucleus composed of water ice, putative unclathrated CO2 ice and silicate dust in specified proportions. Selective sublimation of any free CO2 ice present in a comet would produce a surface layer of water ice and dust overlying the original CO2 rich material. This surface layer reduces the temperature of buried CO2 ice and restricts the outflow of gaseous CO2. On each orbit, water sublimation at smaller heliocentric distances temporarily reduces the thickness of the water ice and dust layer and liberates dust. The model includes the effects of nucleus rotation, arbitrary orientation of the rotation axis, latitude, heat conduction into the deep interior of the nucleus and restriction of CO2 gas outflow by the water ice and dust layer. The effects of the permeability of the surface water ice layer, the nucleus rotation rate, and the latitude were investigated. Comparison of these and similar results with observations could yield information regarding the permeability and chemical composition of cometary material and suggest sampling strategies to minimize fractionation effects.

Fanale, F. P.↗

Rate limits in silicon sheet growth - The connections between vertical and horizontal methods

Meniscus-defined techniques for the growth of thin silicon sheets fall into two categories: vertical and horizontal growth. The interactions of the temperature field and the crystal shape are analyzed for both methods using two-dimensional finite-element models which include heat transfer and capillarity. Heat transfer in vertical growth systems is dominated by conduction in the melt and the crystal, with almost flat melt/crystal interfaces that are perpendicular to the direction of growth. The high axial temperature gradients characteristic of vertical growth lead to high thermal stresses. The maximum growth rate is also limited by capillarity which can restrict the conduction of heat from the melt into the crystal. In horizontal growth the melt/crystal interface stretches across the surface of the melt pool many times the crystal thickness, and low growth rates are achievable with careful temperature control. With a moderate axial temperature gradient in the sheet a substantial portion of the latent heat conducts along the sheet and the surface of the melt pool becomes supercooled, leading to dendritic growth. The thermal supercooling is surpressed by lowering the axial gradient in the crystal; this configuration is the most desirable for the growth of high quality crystals. An expression derived from scaling analysis relating the growth rate and the crucible temperature is shown to be reliable for horizontal growth.

Thomas, Paul D.↗

Microgravity ignition experiment

The purpose of this project is to develop a flight ready apparatus of the microgravity ignition experiment for the GASCan 2 program. This involved redesigning, testing, and making final modifications to the existing apparatus. The microgravity ignition experiment is intended to test the effect of microgravity on the time to ignition of a sample of alpha-cellulose paper. An infrared heat lamp is used to heat the paper sample within a sealed canister. The interior of the canister was redesigned to increase stability and minimize conductive heat transfer to the sample. This design was fabricated and tested and a heat transfer model of the paper sample was developed.

Motevalli, Vahid↗

Quantitative Global Heat Transfer in a Mach-6 Quiet Tunnel

This project developed quantitative methods for obtaining heat transfer from temperature sensitive paint (TSP) measurements in the Mach-6 quiet tunnel at Purdue, which is a Ludwieg tube with a downstream valve, moderately-short flow duration and low levels of heat transfer. Previous difficulties with inferring heat transfer from TSP in the Mach-6 quiet tunnel were traced to (1) the large transient heat transfer that occurs during the unusually long tunnel startup and shutdown, (2) the non-uniform thickness of the insulating coating, (3) inconsistencies and imperfections in the painting process and (4) the low levels of heat transfer observed on slender models at typical stagnation temperatures near 430K. Repeated measurements were conducted on 7 degree-half-angle sharp circular cones at zero angle of attack in order to evaluate the techniques, isolate the problems and identify solutions. An attempt at developing a two-color TSP method is also summarized.

Sullivan, John P.↗

MHD Modeling of Coronal Loops: the Transition Region Throat

Context. The expansion of coronal loops in the transition region may considerably influence the diagnostics of the plasma emission measure. The cross-sectional area of the loops is expected to depend on the temperature and pressure, and might be sensitive to the heating rate. Aims. The approach here is to study the area response to slow changes in the coronal heating rate, and check the current interpretation in terms of steady heating models. Methods. We study the area response with a time-dependent 2D magnetohydrodynamic (MHD) loop model, including the description of the expanding magnetic field, coronal heating and losses by thermal conduction, and radiation from optically thin plasma. We run a simulation for a loop 50 Mm long and quasi-statically heated to about 4 millikelvin. Results. We find that the area can change substantially with the quasi-steady heating rate, e.g., by approx. 40% at 0.5 millikelvin as the loop temperature varies between 1 millikelvin and 4 millikelvin, and, therefore, affects the interpretation of the differential emission measure vs. temperature (DEM(T)) curves.

magnetohydrodynamics (MHD)↗

Analysis of cure in composites processing

Finite element analysis is a general numerical tool for solving the field equations of engineering practice, and this paper demonstrates its use in modeling the nonisothermal cure of pultruded composite material. A very simple grid is used in this case to model a narrow strip of material, and this grid is then solved using a time-stepping transient algorithm to simulate the passage of the strip along the pultruder die. As time proceeds, heat is conducted into the strip from the heated boundaries at the die walls, and cure proceeds at a rate dependent on the local temperature. The computer model can be used to minimize the time needed for sufficient cure, and helps avoid such processing errors as undercure or thermal degradation.

Aylward, L.↗

A model to separate conduction and radiation in high temperature laser flash measurements for semi-transparent materials

Laser flash analysis (LFA) is a widely used technique in the scientific and industrial communities for measuring the thermal diffusivity and conductivity of materials. However, traditional fitting methods based solely on heat conduction do not account for radiation heat transfer in transparent and semi-transparent materials, which becomes significant at high temperatures. In this study, we have developed a coupled radiation-conduction model that considers the optical properties of materials to accurately extract thermal diffusivity from LFA measurements. We conducted LFA measurements on fused silica glass samples of varying thicknesses and within the temperature range of 25 to 1000 °C and fitted the data with the model. Since the coupled model requires optical properties as input parameters, which are often unavailable for new materials, we have also proposed a modified transparent model. Compared to the coupled model, results from the modified transparent model generally overestimate thermal conductivity, but the errors are acceptable in silica glass with small thickness. Here, we also investigated the applicability of the modified transparent model with materials of different optical properties.

14 SOLAR ENERGY↗

A Modular Conjugate Heat Transfer Optimization Framework for Thermal Management of Electric Aircraft

Conjugate heat transfer (CHT) analysis and optimization is a powerful method for improving thermal management, as it simultaneously resolves the temperature distribution in both fluid and solid domains. This paper presents a modular, discrete adjoint-based CHT optimization capability integrated within the OpenMDAO/MPhys framework. A unique feature of the proposed framework is its flexibility to extend to multidisciplinary optimization, including aero-structural-thermal applications. The fluid domain is modeled using a finite-volume Computational Fluid Dynamics (CFD) solver, and the solid domain with a conduction heat transfer solver. A mixed Neumann-Dirichlet boundary condition is developed to enable full submersion of the solid geometry within the fluid domain, while ensuring consistent temperature and heat flux coupling at the CHT interface. Gradient-based optimization is performed; the gradients are efficiently computed using the discrete adjoint solvers implemented in DAFoam. To demonstrate the method, this paper considers two cases related to electric aircraft thermal management: a U-bend heat exchanger and an actively cooled battery pack. The U-bend case aims to minimize pressure loss while maximizing heat flux by changing the pipe geometry. The optimized design reduces pressure loss by 52.7% and increases total heat flux by 2.3%. In the battery pack case, a 3-by-3 cell configuration is cooled by ambient airflow, with constant heat generation prescribed in the cells. The battery casing shape serves as the design variable, and the objective function is a weighted sum of pressure loss and pack weight, subject to a maximum temperature constraint. The optimized design achieves a 44.6% reduction in pressure loss and a 1.5% reduction in weight, while satisfying the thermal constraint. To ensure the reliability of the optimized designs, this study validates coarse-mesh, steady-state predictions against fine-mesh unsteady simulations, demonstrating consistency within acceptable errors. This work demonstrates the potential of the developed framework to enable rapid, high-fidelity design of thermal management systems for electric aircraft.

heat transfer↗

Environmental Loss Characterization of an Advanced Stirling Convertor (ASC-E2) Insulation Package Using a Mock Heater Head

The U.S. Department of Energy (DOE) and Lockheed Martin Space Systems Company (LMSSC) have been developing the Advanced Stirling Radioisotope Generator (ASRG) for use as a power system for space science missions. This generator would use two highefficiency Advanced Stirling Convertors (ASCs), developed by Sunpower Inc. and NASA Glenn Research Center (GRC). As part of ground testing of these ASCs, different operating conditions are used to simulate expected mission conditions. These conditions require achieving a specified electrical power output for a given net heat input. While electrical power output can be precisely quantified, thermal power input to the Stirling cycle cannot be directly measured. In an effort to improve net heat input predictions, the Mock Heater Head was developed with the same relative thermal paths as a convertor using a conducting rod to represent the Stirling cycle and tested to provide a direct comparison to numerical and empirical models used to predict convertor net heat input. The Mock Heater Head also served as the pathfinder for a higher fidelity version of validation test hardware, known as the Thermal Standard. This paper describes how the Mock Heater Head was tested and utilized to validate a process for the Thermal Standard.

Schifer, Nicholas A.↗

CFD modeling of crystallization fouling with CO 2 desorption incorporated for a falling-film evaporator in thermal desalination

Seawater flowing around horizontal tubes in falling-film evaporators is a common configuration for thermal desalination. Heat is transferred from in-tube steam condensation to the shell-side seawater by conduction through the tube wall and scaling layer, and conduction and convection to the evaporating liquid film. CO 2 is simultaneously released from seawater and mixed with the produced vapor. Fouling coupled with CO 2 desorption has deleterious impacts on both evaporation and condensation. The evaluation of spatiotemporal dependent crystallization fouling and CO 2 desorption is vital to selecting optimal operating conditions. Here, in this work, a CO 2 desorption model is integrated into CFD modeling for predictions of coupled heat and mass transfer, to understand and predict scale formation and CO 2 desorption with local and transient profiles of temperature, carbonate species concentrations, pH, and total alkalinity. The porosity of the scale layer is experimentally determined as 67.7% and invoked in the calculation of effective thermal conductivity. The simulation results reveal that high steam temperatures increase seawater total alkalinity and accelerate scale formation. The scale thickness on the bottom tube reaches 0.30 μm and 5.24 μm for steam temperatures of 60 °C and 80 °C, respectively. High salinity leads to a large CO 2 desorption rates. The CO 2 desorption rate increases 43% when the seawater salinity increases from 35 to 55 g/kg. The effects of operating conditions on carbonate speciation and pH have been compared and analyzed. This model can serve as a comprehensive tool for the design of thermal desalination systems and optimal operation.

42 ENGINEERING↗

Boundary-Layer Transition on Hollow Cylinders in Supersonic Free Flight as Affected by Mach Number and a Screwthread Type of Surface Roughness

The effects of Mach number and surface-roughness variation on boundary-layer transition were studied using fin-stabilized hollow-tube models in free flight. The tests were conducted over the Mach number range from 2.8 to 7 at a nominally constant unit Reynolds number of 3 million per inch, and with heat transfer to the model surface. A screwthread type of distributed two-dimensional roughness was used. Nominal thread heights varied from 100 microinches to 2100 microinches. Transition Reynolds number was found to increase with increasing Mach number at a rate depending simultaneously on Mach number and roughness height. The laminar boundary layer was found to tolerate increasing amounts of roughness as Mach number increased. For a given Mach number an optimum roughness height was found which gave a maximum laminar run greater than was obtained with a smooth surface.

James, Carlton S.↗