Search NASA⌕ Search

SEARCH · Search NASA

Results for “HEAT EQUATION”

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 235 records · Page 13

Numerical simulation of an electrothermal deicer pad

A numerical simulation is developed to investigate the removal of ice from composite aircraft blades by means of electrothermal deicing. The model considers one dimensional, unsteady state heat transfer in the composite blade-ice body. The heat conduction equations are approximated by using the Crank-Nicolson finite difference scheme, and the phase change in the ice layer is handled using the Enthalpy method. To solve the system of equations which result, Gauss-Seidel iteration is used. The simulation computes the temperature profile in the composite blade-ice body, as well as the movement of the ice-water interface, as a function of time. This information can be used to evaluate deicer performance. The simulation can also be used to solve a variety of other heat conduction problems involving composite bodies.

Marano, J. J.↗

Subsurface emissions from Mercury - VLA radio observations at 2 and 6 centimeters

Radio observations of Mercury made with the VLA; once in 1986, and on two dates in February of 1988 are presented. These observations are the first to spatially map both hot regions associated with the theoretical hot poles. These 'hot poles' are separated by 180 deg and are a result of the unusual diurnal heating from Mercury's 3/2 spin-orbit resonance and eccentric orbit. The highest resolution data maps areas of the planet as small as 330 km. Maps of total intensity, brightness temperature, polarized intensity, fractional polarization, depolarization, and spectral index are included. It is found that the subsurface thermal emissions from Mercury are characteristic of blackbody reradiation from the solar insolation over a diurnal cycle. These observations to produce full-disk thermophysical models are used. The one-dimensional, time-dependent heat-diffusion equation for all observed disk elements at each epoch in order to constrain thermophsyical parameters and properties of the subsurface material are solved. Using typical lunar values for several of the parameters, it is possible to reproduce the temperature morphology and most of the observed temperature values. It is found that the best-fit models require a substantial contribution of the heat transport in the subsurface to be radiative in nature. The primary difficulty in the models is in predicting the observed temperature differences as a function of frequency.

Ledlow, Michael J.↗

Moisture and heat budgets of a cirrus cloud from aircraft measurements during FIRE

This study is based on the NCAR King Air aircraft and radiosonde observations on 31 October 1986 during the FIRE in Wisconsin over Oshkosh. Aircraft step-up and spiral descent flights are used to obtain kinematic and thermodynamic data. In the step-up maneuver, six different penetrations were made between 1528 and 1616 UTC. Each penetration was about 30 km long separated in the vertical by about 300 m. The time difference between the two spiral soundings was about 43 min. The aircraft descended at a rate of 1.5 m/s during these spiral soundings. Kinematic, cloud physical, and radiometric observations from various instruments are used to estimate the different terms in the moisture- and heat-budget equations. The results show that the advection terms, estimated using the mean longitudinal wind and vertical velocities, and radiative fluxes are important in forming budgets for the cirrus layers. Ice-crystal growth is significant in the upper layers. The maintenance of cirrus can be attributed to relatively warm and moist air advection, radiative cooling at upper levels, and moisture advection in the vertical. Turbulent heat and moisture fluxes are found to be significant in the low levels of cirrus.

Gultepe, Ismail↗

Navier-Stokes turbine heat transfer predictions using two-equation turbulence closures

Navier-Stokes calculations were carried out in order to predict the heat-transfer rates on turbine blades. The calculations were performed using TRAF2D which is a k-epsilon, explicit, finite volume mass-averaged Navier-Stokes solver. Turbulence was modeled using Coakley's q-omega and Chien's k-epsilon two-equation models and the Baldwin-Lomax algebraic model. The model equations along with the flow equations were solved explicitly on a nonperiodic C grid. Implicit residual smoothing (IRS) or a combination of multigrid technique and IRS was applied to enhance convergence rates. Calculations were performed to predict the Stanton number distributions on the first stage vane and blade row as well as the second stage vane row of the SSME high-pressure fuel turbine. The comparison serves to highlight the weaknesses of the turbulence models for use in turbomachinery heat-transfer calculations.

Ameri, Ali A.↗

Turbulence modeling and experiments

The best way of verifying turbulence is to do a direct comparison between the various terms and their models. The success of this approach depends upon the availability of the data for the exact correlations (both experimental and DNS). The other approach involves numerically solving the differential equations and then comparing the results with the data. The results of such a computation will depend upon the accuracy of all the modeled terms and constants. Because of this it is sometimes difficult to find the cause of a poor performance by a model. However, such a calculation is still meaningful in other ways as it shows how a complete Reynolds stress model performs. Thirteen homogeneous flows are numerically computed using the second order closure models. We concentrate only on those models which use a linear (or quasi-linear) model for the rapid term. This, therefore, includes the Launder, Reece and Rodi (LRR) model; the isotropization of production (IP) model; and the Speziale, Sarkar, and Gatski (SSG) model. Which of the three models performs better is examined along with what are their weaknesses, if any. The other work reported deal with the experimental balances of the second moment equations for a buoyant plume. Despite the tremendous amount of activity toward the second order closure modeling of turbulence, very little experimental information is available about the budgets of the second moment equations. Part of the problem stems from our inability to measure the pressure correlations. However, if everything else appearing in these equations is known from the experiment, pressure correlations can be obtained as the closing terms. This is the closest we can come to in obtaining these terms from experiment, and despite the measurement errors which might be present in such balances, the resulting information will be extremely useful for the turbulence modelers. The purpose of this part of the work was to provide such balances of the Reynolds stress and heat flux equations for the buoyant plume.

Shabbir, Aamir↗

A finite element analysis of the freeze/thaw behavior of external artery heat pipes

A two-dimensional finite element model was used to determine the freeze/thaw characteristics of an external artery heat pipe. During startup, the working fluid, which was located in the liquid channel and the circumferential wall grooves, experienced a phase transformation from a solid to a liquid state. The transient heat conduction equations with moving interfacial conditions were solved using the appropriate initial boundary conditions. The modelling results include the cross-sectional temperature distribution and the interfacial or melt front position as a function of time. A fixed grid approach was adopted in the model for the phase-change process during thawing of frozen working fluid. The interfacial position between the liquid and solid regions was found by balancing the latent heat caused by interfacial movement with the heat addition or extraction at the related grid points.

Lu, X. J.↗

Boundary Layer Instabilities Over a Cone-Cylinder-Flare Model at Mach 6

Computations are performed to investigate the boundary-layer instabilities over a sharp cone-cylinder-flare model at zero degrees angle of attack. The model geometry and the flow conditions are selected to match the experiments conducted in the Boeing/AFOSR Mach 6 Quiet Tunnel (BAM6QT) at Purdue University. The geometry consists of a nominally sharp 5-degree half-angle cone, followed by a cylindrical segment and then a 10-degree flare. An axisymmetric separation bubble is generated as a result of the laminar shock/boundary-layer interaction in the cylinder-flare region. The comparison of the laminar flow solution and the schlieren images shows a remarkable agreement between the respective locations of both the boundarylayer edge and the reattachment shock. The predicted heat flux distribution is also in agreement with the measured values downstream of the reattachment location. The analysis of convective and global instabilities is performed for flare half angles equal to 8, 10, and 12 degrees and nosetip radii equal to 0.1, 1, and 5 mm. The linear amplification of first and second Mack mode instabilities that begin to amplify in the cone region are computed with a combination of the parabolized stability equations (PSE) and the harmonic linearized Navier-Stokes equations (HLNSE). The predicted frequency spectra of the surface pressure fluctuations associated with both planar and oblique instability waves are compared with the measured spectra at the various locations of the PCB and Kulite sensors. The comparison shows that the computational analysis captures the distinct lobes within the disturbance amplification spectra measured in the experiments, but some differences in amplification characteristics are noted at low frequencies. Overall, the oblique disturbances are found to be more amplified than the planar disturbances. To our knowledge, this represents the first successful comparison between convective instability analysis and measured surface pressure fluctuations for a hypersonic configuration with a separation bubble. Finally, the global instability analysis shows that the laminar flow becomes supercritical for flare half angles larger than 8 degrees. The unstable global mode for the experimental configuration of a 10 degrees flare and a sharp nosetip cone corresponds to a stationary three-dimensional disturbance that is concentrated in the recirculation region and achieves its maximum growth rate for an azimuthal wavenumber of 5.

Boundary layer transition↗

An approximation theory for the identification of linear thermoelastic systems

An abstract approximation framework and convergence theory for the identification of thermoelastic systems is developed. Starting from an abstract operator formulation consisting of a coupled second order hyperbolic equation of elasticity and first order parabolic equation for heat conduction, well-posedness is established using linear semigroup theory in Hilbert space, and a class of parameter estimation problems is then defined involving mild solutions. The approximation framework is based upon generic Galerkin approximation of the mild solutions, and convergence of solutions of the resulting sequence of approximating finite dimensional parameter identification problems to a solution of the original infinite dimensional inverse problem is established using approximation results for operator semigroups. An example involving the basic equations of one dimensional linear thermoelasticity and a linear spline based scheme are discussed. Numerical results indicate how the approach might be used in a study of damping mechanisms in flexible structures.

Rosen, I. G.↗

Influences of CO2 sublimation/condensation processes on the long-term evolution of the Martian atmosphere

The accumulation process of frozen CO2 within the Martian regolith is analyzed. The boundary conditions are taken to be the periodically variable temperature at the Mars surface and the constant heat flux at the bottom of the regolith layers, some hundred meters thick. The mean value of temperature as well as its amplitude are assumed to be dependent on latitude and could be variable with the geological time scale. The flux of gaseous CO2 within the regolith is controlled by the mechanism of gas transport through the porous medium. An appropriate initial geometry of distribution of pores is assumed. The porosity and flux change due to condensation/sublimation processes. The equations of heat and mass transfer are solved numerically for a given latitude. Next, the local (over the parallel of latitude belt) annual balance of CO2 is calculated. The procedure is repeated over the whole Martian surface and next the total annual net flux is calculated.

Kossacki, K.↗

Extension of boundary-layer heat-transfer theory to cooled turbine blades

An equation for average heat transfer of a surface was derived when the boundary layer changed from laminar to turbulent. Influences on the heat transfer through a laminar boundary layer of Mach number, temperature ratio (gas temperature divided by wall temperature), and exponents of gas-property temperature relations were shown to be relatively small for air with Mach numbers less than 2 and temperature ratios between 1 and 4. Good agreement was obtained with experimental results from cylinders, an airfoil, and turbine blades.

HEAT TRANSFER, AERODYNAMIC↗

Thermal Modelling of the in-Situ Consolidation of Automated Fiber Placement of Thermoplastic Composites

NASA is developing carbon fiber reinforced thermoplastic composites processing methods under the Hi-rate Composites Aircraft Manufacturing (HiCAM) project. The in-situ consolidation automated fiber placement (AFP) of thermoplastics (ICAT) process has high potential to increase manufacturing throughput and lower costs as it combines the repeatability and fast lay down rates of AFP with the thermoformability of thermoplastic materials to achieve part consolidation out of the autoclave. Physics-based process models are employed in the ICAT process development to understand, fundamentally, the thermal response of the carbon fiber (CF)/polyaryletherketone (PAEK) material during the rapid heating and cooling associated with laser-assisted AFP. A one-dimensional (1-D) closed-form model was developed capable of predicting the temperature profile through-the-thickness of the tape material. In addition, a two-dimensional (2-D) model was developed to predict the heat transfer through-the-thickness and in the laydown x-direction of the AFP head motion. The solution of the Fourier heat-transfer equations in the 2-D model is approximated using the explicit finite difference numerical method. The temperature profile during placement of PAEK slit-tape materials at various laser target temperatures and placement speeds were measured during ICAT process development trials at Electroimpact®, Inc. †† and the resulting experimental data is compared with the predictions of the 1-D and 2-D models.

Thermoplastics↗

Turbine Vane External Heat Transfer. Volume 2. Numerical Solutions of the Navier-stokes Equations for Two- and Three-dimensional Turbine Cascades with Heat Transfer

The application of the time-dependent ensemble-averaged Navier-Stokes equations to transonic turbine cascade flow fields was examined. In particular, efforts focused on an assessment of the procedure in conjunction with a suitable turbulence model to calculate steady turbine flow fields using an O-type coordinate system. Three cascade configurations were considered. Comparisons were made between the predicted and measured surface pressures and heat transfer distributions wherever available. In general, the pressure predictions were in good agreement with the data. Heat transfer calculations also showed good agreement when an empirical transition model was used. However, further work in the development of laminar-turbulent transitional models is indicated. The calculations showed most of the known features associated with turbine cascade flow fields. These results indicate the ability of the Navier-Stokes analysis to predict, in reasonable amounts of computation time, the surface pressure distribution, heat transfer rates, and viscous flow development for turbine cascades operating at realistic conditions.

Yang, R. J.↗

HEATING OF THE UPPER ATMOSPHERE

Theoretical time-dependent model of the upper atmosphere through solutions to heat conduction equation for an equilibrium gas in a gravitational field

UPPER ATMOSPHERE↗