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,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS1966

Protuberance heating test program

Results are presented of the protuberance heating test program. Four general protuberance shapes on a flat plate were tested. Presentation and evaluation of the data both on the protuberance and in the wake regions are made. The test program is an extension of the general protuberance heat transfer test. The additional series of tests was conducted to define the extent of wake heating and to assess the effects of Reynolds number variation on heating both on and around the protuberances. The protuberance models were mounted near the forward end of a six-foot instrumented test plate with stringers that simulated interstage and skirt structure of the Saturn S-4B stage. The tests were performed at Mach numbers of 2.49, 3.51, and 4.44. Reynolds numbers per foot of 3 million and 1.5 million were used for the two lower Mach numbers and 3 million for a Mach number of 4.44. The test Mach numbers simulated the Saturn S-4B flight conditions during the most severe aerodynamic heating period. The test Reynolds numbers were somewhat higher than the flight values, but lower values could not be used because of tunnel and instrumentation limitations. Oil flow runs were made on two representative models at various combinations of Mach number and Reynolds number to help define the extent of wake heating.

Sieker, W. D.↗
NASA NTRS1975

Heating of solid targets with laser pulses

Analytical and numerical solutions to the heat-conduction equation are obtained for the heating of absorbing media with pulsed lasers. The spatial and temporal form of the temperature is determined using several different models of the laser irradiance. Both surface and volume generation of heat are discussed. It is found that if the depth of thermal diffusion for the laser-pulse duration is large compared to the optical-attenuation depth, the surface- and volume-generation models give nearly identical results. However, if the thermal-diffusion depth for the laser-pulse duration is comparable to or less than the optical-attenuation depth, the surface-generation model can give significantly different results compared to the volume-generation model. Specific numerical results are given for a tungsten target irradiated by pulses of different temporal durations and the implications of the results are discussed with respect to the heating of metals by picosecond laser pulses.

Bechtel, J. H.↗
NASA NTRS1978

Electron heating within interaction zones of simple high-speed solar wind streams

In the present paper, electron heating within the high-speed portions of three simple stream-stream interaction zones is studied to further our understanding of the physics of heat flux regulation in interplanetary space. To this end, the thermal signals present in the compressions at the leading edges of the simple high-speed streams are analyzed, showing that the data are inconsistent with the Spitzer conductivity. Instead, a polynomial law is found to apply. Its implication concerning the mechanism of interplanetary heat conduction is discussed, and the results of applying this conductivity law to high-speed flows inside of 1 AU are studied. A self-consistent model of the radial evolution of electrons in the high-speed solar wind is proposed.

Feldman, W. C.↗
NASA NTRS1984

Estimation of electric fields and current from ground-based magnetometer data

Recent advances in numerical algorithms for estimating ionospheric electric fields and currents from groundbased magnetometer data are reviewed and evaluated. Tests of the adequacy of one such algorithm in reproducing large-scale patterns of electrodynamic parameters in the high-latitude ionosphere have yielded generally positive results, at least for some simple cases. Some encouraging advances in producing realistic conductivity models, which are a critical input, are pointed out. When the algorithms are applied to extensive data sets, such as the ones from meridian chain magnetometer networks during the IMS, together with refined conductivity models, unique information on instantaneous electric field and current patterns can be obtained. Examples of electric potentials, ionospheric currents, field-aligned currents, and Joule heating distributions derived from ground magnetic data are presented. Possible directions for future improvements are also pointed out.

Kamide, Y.↗
NASA NTRS1985

Analysis of ocean surface wind fields using SEASAT-A scatterometer data

In addition to the SEASAT forecast impact studies previously conducted with the GLA model, a study was also made of the global analyses of surface wind stress, wind stress curl, and surface sensible and latent heat fluxes, which are generated as a by-product of the assimilation process. Both instantaneous analyses of these quantities and 15-day averages of each of these fields were found to be in good qualitative agreement with published fields.

Kalnay, E.↗
NASA NTRS1991

Studies of shock/shock interaction on smooth and transpiration-cooled hemispherical nosetips in hypersonic flow

Experimental studies are conducted to examine the utilization of transpiration cooling to reduce the peak-heating loads in areas of shock/shock interaction. Smooth and transpiration-cooled nosetip models, 12 inches in diameter, were employed in these studies, which focused on defining the pressure distributions and heat transfer in type III and IV interaction areas. Transpiration cooling was determined to significantly increase the size of the shock layer and to move the peak-heating point around the body. A transpiration-cooling rate of more than 30 percent of the freestream maximum flux did not lower the peak-heating level more than 10 percent, but the integrated heating loads were reduced.

Holden, M. S.↗
NASA NTRS1992

Turbine disk cavity aerodynamics and heat transfer

Experiments were conducted to define the nature of the aerodynamics and heat transfer for the flow within the disk cavities and blade attachments of a large-scale model, simulating the Space Shuttle Main Engine (SSME) turbopump drive turbines. These experiments of the aerodynamic driving mechanisms explored the following: (1) flow between the main gas path and the disk cavities; (2) coolant flow injected into the disk cavities; (3) coolant density; (4) leakage flows through the seal between blades; and (5) the role that each of these various flows has in determining the adiabatic recovery temperature at all of the critical locations within the cavities. The model and the test apparatus provide close geometrical and aerodynamic simulation of all the two-stage cavity flow regions for the SSME High Pressure Fuel Turbopump and the ability to simulate the sources and sinks for each cavity flow.

Johnson, B. V.↗
NASA NTRS1993

Thermal Conductivity And Expansion Of Graphite/Copper

Report describes fabrication of graphite-fiber/copper-matrix composite plates, measurements of thermal conductivities of plates at temperatures from ambient to 1,073 K, and measurements of thermal expansions of plates from ambient temperature to 1,050 K. Composites promising lightweight, high-thermal-conductivity materials proposed for use in heat exchangers and other heat-transfer components of power systems in spacecraft and hypersonic aircraft. Graphite/copper also of interest as model composite material.

Mcdanels, David L.↗
NASA NTRS1998

Observation and Modeling of the Solar Transition Region: A Quasi-Static Loops Model with Implications for Heating the Lower Transition Region - 1

We report on observations of the solar atmosphere in several extreme ultraviolet and far-ultraviolet bandpasses obtained by the Multi-Spectral Solar Telescope Array, a rocket borne spectroheliograph, on flights in 1987, 1991, and 1994, spanning the last solar maximum. Quiet sun emission observed in the 171 A - 175 A bandpass, which includes lines of 0Ov, O vi, Fe ix, and Fe x, has been analyzed to test models of the temperatures and geometries of the structures responsible for this emission. Analyses of intensity variations above the solar limb reveal scale heights consistent with a quiet sun plasma temperature of 500 000 K less than or equal to T(sub e) less than or equal to 800 000 K. Intensity modulations in the quiet sun are observed to occur on a scale comparable to the supergranular scale. The structures responsible for the quiet sun EUV emission are modeled as small quasi-static loops. We find that the emission predicted by loop models with maximum temperatures between 700 000 K and 900 000 K are consistent with our observations. We also present a preliminary comparison of the predictions of our models with observations of diagnostic spectral line ratios obtained from previous observers. We discuss the implications a distribution of loops of the type we model here would have for heating the lower transition region. Finally, in fight of the models we calculate here, we briefly review the current state of knowledge concerning the contributions thermal conduction from coronal (T(sub e) greater than or equal to 10(exp 6) K) and upper transition region (10(exp 5) K less than T(sub e) less than 10(esp 6) K) structures make to lower transition region emission. We argue that the evidence which has lead many authors to conclude that the interface of hotter and cooler plasmas makes a negligible contribution to lower transition region emission is much less compelling in light of recent observations and analyses. We further argue that it is the interface of chromospheric material with structures such as loops that have sub-coronal peak temperatures (i.e. less than 900 000 K) that makes the dominant contribution to lower transition region emission in the quiet sun.

Oluseyi, Hakeem M.↗
NASA NTRS2013

Implicit Coupling Approach for Simulation of Charring Carbon Ablators

This study demonstrates that coupling of a material thermal response code and a flow solver with nonequilibrium gas/surface interaction for simulation of charring carbon ablators can be performed using an implicit approach. The material thermal response code used in this study is the three-dimensional version of Fully Implicit Ablation and Thermal response program, which predicts charring material thermal response and shape change on hypersonic space vehicles. The flow code solves the reacting Navier-Stokes equations using Data Parallel Line Relaxation method. Coupling between the material response and flow codes is performed by solving the surface mass balance in flow solver and the surface energy balance in material response code. Thus, the material surface recession is predicted in flow code, and the surface temperature and pyrolysis gas injection rate are computed in material response code. It is demonstrated that the time-lagged explicit approach is sufficient for simulations at low surface heating conditions, in which the surface ablation rate is not a strong function of the surface temperature. At elevated surface heating conditions, the implicit approach has to be taken, because the carbon ablation rate becomes a stiff function of the surface temperature, and thus the explicit approach appears to be inappropriate resulting in severe numerical oscillations of predicted surface temperature. Implicit coupling for simulation of arc-jet models is performed, and the predictions are compared with measured data. Implicit coupling for trajectory based simulation of Stardust fore-body heat shield is also conducted. The predicted stagnation point total recession is compared with that predicted using the chemical equilibrium surface assumption

Ablation↗
NASA NTRS2009

Combined Modeling of Acceleration, Transport, and Hydrodynamic Response in Solar Flares: The Numerical Model - 1

Acceleration and transport of high-energy particles and fluid dynamics of atmospheric plasma are interrelated aspects of solar flares, but for convenience and simplicity they were artificially separated in the past. We present here self consistently combined Fokker-Planck modeling of particles and hydrodynamic simulation of flare plasma. Energetic electrons are modeled with the Stanford unified code of acceleration, transport, and radiation, while plasma is modeled with the Naval Research Laboratory flux tube code. We calculated the collisional heating rate directly from the particle transport code, which is more accurate than those in previous studies based on approximate analytical solutions. We repeated the simulation of Mariska et al. with an injection of power law, downward-beamed electrons using the new heating rate. For this case, a -10% difference was found from their old result. We also used a more realistic spectrum of injected electrons provided by the stochastic acceleration model, which has a smooth transition from a quasi-thermal background at low energies to a non thermal tail at high energies. The inclusion of low-energy electrons results in relatively more heating in the corona (versus chromosphere) and thus a larger downward heat conduction flux. The interplay of electron heating, conduction, and radiative loss leads to stronger chromospheric evaporation than obtained in previous studies, which had a deficit in low-energy electrons due to an arbitrarily assumed low-energy cutoff. The energy and spatial distributions of energetic electrons and bremsstrahlung photons bear signatures of the changing density distribution caused by chromospheric evaporation. In particular, the density jump at the evaporation front gives rise to enhanced emission, which, in principle, can be imaged by X-ray telescopes. This model can be applied to investigate a variety of high-energy processes in solar, space, and astrophysical plasmas.

STOCHASTIC ACCELERATION MODEL↗
DOE OSTI2021

CFD modeling of pre-spark heat release in a boosted direct-injection spark-ignition engine

Accurate predictions of low-temperature heat release (LTHR) are critical for modeling auto-ignition processes in internal combustion engines. While LTHR is typically obscured by deflagration, extremely late ignition phasing can lead to LTHR prior to the spark, a behavior known as pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection spark-ignition engine was studied using 3-D computational fluid dynamics (CFD) and detailed chemical kinetics. The turbulent combustion was modeled via a hybrid approach that incorporates the G-equation model for tracking the turbulent flame front, and the well-stirred reactor model with detailed chemistry for assessing the low-temperature reactions in unburnt gas. Simulations were conducted using Co-Optima alkylate and E30 fuels at operating conditions characterized by different PSHR intensities. The predicted in-cylinder pressure and heat release rate were found to agree well with experiments. It was found the estimate of previous-cycle trapped residuals is of utmost importance for capturing PSHR correctly. A simulation best practice was developed which keeps the detailed chemistry solver active throughout the entire simulation, allowing to track the evolution of intermediate species from one cycle to the next. Following the validation, the dynamics of PSHR were analyzed in detail employing the pressure-temperature (P-T) trajectory framework. It was shown that PSHR correlated with the first-stage ignition delay of the fuel, hence showing close relation to the in-cylinder P-T trajectory and the chemical kinetics. Besides, it was indicated that LTHR is a self-limiting process that has the effect of attenuating the thermal stratification in the combustion chamber. Furthermore, it was observed the occurrence of PSHR caused the P-T trajectory of end-gas to overlap with the negative temperature coefficient region of the fuel’s ignition-delay maps. This effect was more significant in the fuel-rich regions where engine knock tendency would be generally higher, with potential implications on knock control and mitigation.

42 ENGINEERING↗
NASA NTRS1987

Effects of thermospheric motions on the polar wind - A time-dependent numerical study

This paper presents the results of a numerical study investigating the effects of low-altitude short-duration plasma heating episodes (such as horizontal frictional heating) on transient heavy ion outflows from the polar ionosphere using a time-dependent model. In a previous set of calculations, where the effects of high-altitude transient heat sources were investigated, the low-altitude neutral atmosphere (z less than 800 km) acted as a very efficient heat sink, absorbing most of the energy conducted to this region and thus preventing large differences between the ion and neutral temperatures. One purpose of the present study was to investigate whether realistic, specified low-altitude frictional heating rates, based on published experimental data, can result in elevated ion temperatures and lift a part of the heavy ion population over the gravitational barrier before the extra energy is lost to the neutral atmosphere. In these model calculations the specified heating generated an upflowing O(+) disturbance (with a normalized peak flux of about 100 million/sq cm per sec) moving upward along the open magnetic field line with a velocity of about 2 km/s. The spatial extent of the disturbance was approximately equal to the disturbance velocity times the duration of the heat source. The model results provide a possible mechanism to explain recent European Incoherent Scatter observations of high-latitude O(+) upwellings.

Gombosi, T. I.↗
NASA NTRS1992

Chemical differentiation of a convecting planetary interior: Consequences for a one-plate planet such as Venus

Chemically depleted mantle forming a buoyant, refractory layer at the top of the mantle can have important implications for the evolution of the interior and surface. On Venus, the large apparent depths of compensation for surface topographic features might be explained if surface topography were supported by variations in the thickness of a 100-200 km thick chemically buoyant mantle layer or by partial melting in the mantle at the base of such a layer. Long volcanic flows seen on the surface may be explained by deep melting that generates low-viscosity MgO-rich magmas. The presence of a shallow refractory mantle layer may also explain the lack of volcanism associated with rifting. As the depleted layer thickens and cools, it becomes denser than the convecting interior and the portion of it that is hot enough to flow can mix with the convecting mantle. Time dependence of the thickness of a depleted layer may create episodic resurfacing events as needed to explain the observed distribution of impact craters on the venusian surface. We consider a planetary structure consisting of a crust, depleted mantle layer, and a thermally and chemically well-mixed convecting mantle. The thermal evolution of the convecting spherical planetary interior is calculated using energy conservation: the time rate of change of thermal energy in the interior is equated to the difference in the rate of radioactive heat production and the rate of heat transfer across the thermal boundary layer. Heat transfer across the thermal boundary layer is parameterized using a standard Nusselt number-Rayleigh number relationship. The radioactive heat production decreases with time corresponding to decay times for the U, Th, and K. The planetary interior cools by the advection of hot mantle at temperature T interior into the thermal boundary layer where it cools conductively. The crust and depleted mantle layers do not convect in our model so that a linear conductive equilibrium temperature distribution is assumed. The rate of melt production is calculated as the product of the volume flux of mantle into the thermal boundary layer and the degree of melting that this mantle undergoes. The volume flux of mantle into the thermal boundary layer is simply the heat flux divided by amount of heat lost in cooling mantle to the average temperature in the thermal boundary layer. The degree of melting is calculated as the temperature difference above the solidus, divided by the latent heat of melting. A maximum degree of melting is prescribed corresponding to the maximum amount of basaltic melt that the mantle can initially generate. As the crust thickens, the pressure at the base of the crust becomes high enough and the temperature remains low enough for basalt to transform to dense eclogite.

Parmentier, E. M.↗
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Surface Heater Fabrication Using Micro-Lithography for Transpiration Cooling Heat Transfer Coefficient Measurements

Abstract Transpiration cooling is able to provide more uniform coolant coverage than film cooling to effectively protect the component surface from contacting the hot gas. Due to numerous coolant ejection outlets within a small area at the target surface, the experimental thermo-fluid investigation on transpiration cooling becomes a significant challenge. Two classic methods to investigate film cooling, the steady-state foil heater method and the transient thermography technique, fail for transpiration cooling because the foil heater would block numerous coolant outlets, and the semi-infinite solid conduction model no longer holds for porous plates. In this study, a micro-lithography method to fabricate a silver coil pattern on top of the additively manufactured polymer porous media as the surface heater was proposed. The circuit was deliberately designed to cover the solid surface in a combination of series connection, and parallel connection to ensure the power in each unit cell area at the target surface was identical. With uniform heat flux generation, the steady-state tests were conducted to obtain distributions of a pair of parameters, adiabatic cooling effectiveness, and heat transfer coefficient (HTC). The results showed that the adiabatic cooling effectiveness could reach 0.65 with a blowing ratio lower than 0.5. Meanwhile, the heat transfer coefficient ratio (hf/h0) of transpiration cooling was close to 1 with a small blowing ratio at 0.125. A higher HTC ratio was observed for smaller pitch-to-diameter cases due to more turbulence intensity generated at the target surface.

Engineering↗
NASA NTRS2017

Radiative Heat Transfer Modeling in Fibrous Porous Media

Phenolic-Impregnated Carbon Ablator (PICA) was developed at NASA Ames Research Center as a lightweight thermal protection system material for successful atmospheric entries. The objective of the current work is to compute the effective radiative conductivity of fibrous porous media, such as preforms used to make PICA, to enable the efficient design of materials that can meet the thermal performance goals of forthcoming space exploration missions.

Radiative↗
NASA NTRS1985

The perturbed neutral circulation in the vicinity of a symmetric stable auroral arc

In the present investigation, a sophisticated high-resolution time-dependent numerical model is used to study the neutral response to changes in the E region electrodynamical forcing and particle heating associated with the appearance of a stable auroral arc. The initial condition is a nighttime ionosphere in the presence of a background precipitation of electrons associated with diffuse auroras. A simulation study is conducted of the perturbation in the neutral circulation. The simulations were performed on the basis of a model which describes nonlinear, nonhydrostatic viscous flow on a rotating atmosphere.

Walterscheid, R. L.↗