Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

An annual model of SSM/I radiobrightness for dry soil

An annual model is presented of the temperature structure within a homogeneous, dry soil halfspace that is subject to both diurnal and annual insolation, radiant heating from the atmosphere, sensible heat exchange with the atmosphere, and radiant cooling. The thermal constitutive properties of the soil are assumed to be constant so that the heat flow equation can be solved analytically. For computational economy, a variable time interval Laplace transform method is developed to predict the temperature.

Liou, Yuei-An↗
NASA NTRS1953

Transient temperature distributions in simple conducting bodies steadily heated through a laminar boundary layer

An analysis is made of the transient heat-conduction effects in three simple semi-infinite bodies: the flat insulated plate, the conical shell, and the slender solid cone. The bodies are assumed to have constant initial temperatures and, at zero time, to begin to move at a constant speed and zero angle of attack through a homogeneous atmosphere. The heat input is taken as that through a laminar boundary layer. Radiation heat transfer and transverse temperature gradients are assumed to be zero. The appropriate heat-conduction equations are solved by an iteration method, the zeroeth-order terms describing the situation in the limit of small time. The method is presented and the solutions are calculated to three orders which are sufficient to give reasonably accurate results when the forward edge has attained one-half the total temperature rise (nose half-rise time). Flight Mach number and air properties occur as parameters in the result. Approximate expressions for the extent of the conduction region and nose half-rise times as functions of the parameters of the problem are presented. (author)

Parker, Hermon M↗
NASA NTRS1957

Theoretical rocket performance of JP-4 fuel with mixtures of liquid ozone and fluorine

Data were estimated by means of a heat correction equation using data for JP-4 fuel with mixtures of oxygen and flourine. The estimated data were checked for several cases by direct calculations. The difference in specific impulse between the estimated and directly calculated values was from 0.2 to 0.8 pound-second per pound. The maximum value of specific impulse was 334.9 pound-seconds per pound for a combustion-chamber pressure of 600 pounds per square inch absolute and an exit pressure of 1 atmosphere.

Huff, Vearl N↗
NASA NTRS1995

Supersonic solutions of the solar-wind equations

We re-examine the inviscid solar-wind equations with heat conduction from below, and establish a fundamentally new approach for finding stellar solutions. Although the problem is fourth order, only two independent integration constants can be assigned, since two boundary conditions that are required to specify well-behaved supersonic solutions determine the values of the other two constants. The solar-wind-model of Noble and Scarf are essentially accurate for practical purposes, but in a fundamental sense they are not self-consistent. At the supersonic level, the ratio of thermal energy to gravitational potential, (kTr/GMm)(sub s), must lie in a narrow range, between 0.4375 and about 0.40, to permit well-behaved solutions for ionized hydrogen. For a supersonic solution this ratio uniquely determines the constants A and epsilon, which in turn govern not only the escape flux and the energy flux per particle but also T'(lambda approaches 0), the temperature gradient at infinity. It is very difficult to find solutions by trial and error, without employing T'(0) as a guide to the iteration of A and epsilon, because not only must (kTr/GMm)(sub s) be chosen within the acceptable range but then either A or epsilon must be guessed with phenomenal accuracy. Otherwise the numerical integration of temperature and expansion velocity do not yield acceptable values at infinity.

Chamberlain, Joseph W.↗
NASA NTRS1997

Numerical Modeling of HgCdTe Solidification: Effects of Phase Diagram, Double-Diffusion Convection and Microgravity Level

Melt convection, along with species diffusion and segregation on the solidification interface are the primary factors responsible for species redistribution during HgCdTe crystal growth from the melt. As no direct information about convection velocity is available, numerical modeling is a logical approach to estimate convection. Furthermore influence of microgravity level, double-diffusion and material properties should be taken into account. In the present study, HgCdTe is considered as a binary alloy with melting temperature available from a phase diagram. The numerical model of convection and solidification of binary alloy is based on the general equations of heat and mass transfer in two-dimensional region. Mathematical modeling of binary alloy solidification is still a challenging numericial problem. A Rigorous mathematical approach to this problem is available only when convection is not considered at all. The proposed numerical model was developed using the finite element code FIDAP. In the present study, the numerical model is used to consider thermal, solutal convection and a double diffusion source of mass transport.

Bune, Andris V.↗
NASA NTRS2000

Radiation Heat Transfer Between Diffuse-Gray Surfaces Using Higher Order Finite Elements

This paper presents recent work on developing methods for analyzing radiation heat transfer between diffuse-gray surfaces using p-version finite elements. The work was motivated by a thermal analysis of a High Speed Civil Transport (HSCT) wing structure which showed the importance of radiation heat transfer throughout the structure. The analysis also showed that refining the finite element mesh to accurately capture the temperature distribution on the internal structure led to very large meshes with unacceptably long execution times. Traditional methods for calculating surface-to-surface radiation are based on assumptions that are not appropriate for p-version finite elements. Two methods for determining internal radiation heat transfer are developed for one and two-dimensional p-version finite elements. In the first method, higher-order elements are divided into a number of sub-elements. Traditional methods are used to determine radiation heat flux along each sub-element and then mapped back to the parent element. In the second method, the radiation heat transfer equations are numerically integrated over the higher-order element. Comparisons with analytical solutions show that the integration scheme is generally more accurate than the sub-element method. Comparison to results from traditional finite elements shows that significant reduction in the number of elements in the mesh is possible using higher-order (p-version) finite elements.

Gould, Dana C.↗
NASA NTRS2000

Combustion Response of AP Composite Propellants

The Cohen & Strand model for AP composite propellants is applied as boundary conditions, one for AP and one for binder, in solving the heat conduction equation in each to compute linear and non-linear combustion response properties for each and for the aggregate propellant. Iterations couple AP and binder through the quasi-steady flame processes. Illustrative results for linear response functions (pressure-coupled and velocity-coupled) are presented for a monomodal AP propellant varying AP size, pressure and crossflow speed, and examples of non-linear responses to arbitrary waveforms are shown. A quantitative comparison with response function data is limited to one well-characterized research formulation. Mechanistic implications are discussed, including recommendations for future modeling work.

Shusser, Michael↗
NASA NTRS1999

Analysis and Testing of High Temperature Fibrous Insulation for Reusable Launch Vehicles

Analytical models were developed to model the heat transfer through high-temperature fibrous insulation used in metallic thermal protection systems on reusable launch vehicles. The optically thick approximation was used to simulate radiation heat transfer through the insulation. Different models for gaseous conduction and solid conduction in the fibers, and for combining the various modes of heat transfer into a local, volume-averaged, thermal conductivity were considered. The governing heat transfer equations were solved numerically, and effective thermal conductivities were calculated from the steady-state results. An experimental apparatus was developed to measure the apparent thermal conductivity of insulation subjected to pressures, temperatures and temperature gradients representative of re-entry conditions for launch vehicles. The apparent thermal conductivity of an alumina fiber insulation was measured at nominal densities of 24, 48 and 96 kg/cu m. Data were obtained at environmental pressures from 10(exp 4) to 760 torr, with the insulation cold side maintained at room temperature and its hot side temperature varying up to 1000 C. The experimental results were used to evaluate the analytical models. The best analytical model resulted in effective thermal conductivity predictions that were within 8% of experimental results.

Daryabeigi, Kamran↗
NASA NTRS2000

Thermographic Imaging of Defects in Anisotropic Composites

Composite materials are of increasing interest to the aerospace industry as a result of their weight versus performance characteristics. One of the disadvantages of composites is the high cost of fabrication and post inspection with conventional ultrasonic scanning systems. The high cost of inspection is driven by the need for scanning systems which can follow large curve surfaces. Additionally, either large water tanks or water squirters are required to couple the ultrasonics into the part. Thermographic techniques offer significant advantages over conventional ultrasonics by not requiring physical coupling between the part and sensor. The thermographic system can easily inspect large curved surface without requiring a surface following scanner. However, implementation of Thermal Nondestructive Evaluations (TNDE) for flaw detection in composite materials and structures requires determining its limit. Advanced algorithms have been developed to enable locating and sizing defects in carbon fiber reinforced plastic (CFRP). Thermal Tomography is a very promising method for visualizing the size and location of defects in materials such as CFRP. However, further investigations are required to determine its capabilities for inspection of thick composites. In present work we have studied influence of the anisotropy on the reconstructed image of a defect generated by an inversion technique. The composite material is considered as homogeneous with macro properties: thermal conductivity K, specific heat c, and density rho. The simulation process involves two sequential steps: solving the three dimensional transient heat diffusion equation for a sample with a defect, then estimating the defect location and size from the surface spatial and temporal thermal distributions (inverse problem), calculated from the simulations.

Plotnikov, Y. A.↗
NASA NTRS2008

An Improved Formulation for Calorimetric Emittance Testing of Spacecraft Thermal Control Coatings

Spacecraft often really heavily on passive thermal control to maintain operating temperature. An important parameter in the spacecraft heat balance equation is the emittance of thermal control coatings as a function of coating temperature. One method for determining the emittance of spacecraft thermal control from elevated temperature to cryogenic temperatures relies on a calorimetric technique. The fundamental equation governing this test method can be found in numerous places in the literature and although it generally provides reasonable results, its formulation is based on a conceptual flaw that only becomes apparent when the sample temperature approaches the wall temperature during testing. This paper investigates the cause for this error and develops the correct formulation for calorimetric emittance testing. Experimental data will also be presented that illustrates the difference between the two formulations and the resulting difference in the calculated emittance.

Kauder, Lonny R.↗
NASA NTRS2008

GFSSP Training Course Lectures

GFSSP has been extended to model conjugate heat transfer Fluid Solid Network Elements include: a) Fluid nodes and Flow Branches; b) Solid Nodes and Ambient Nodes; c) Conductors connecting Fluid-Solid, Solid-Solid and Solid-Ambient Nodes. Heat Conduction Equations are solved simultaneously with Fluid Conservation Equations for Mass, Momentum, Energy and Equation of State. The extended code was verified by comparing with analytical solution for simple conduction-convection problem The code was applied to model: a) Pressurization of Cryogenic Tank; b) Freezing and Thawing of Metal; c) Chilldown of Cryogenic Transfer Line; d) Boil-off from Cryogenic Tank.

Majumdar, Alok K.↗
NASA NTRS2012

Analysis and Sizing for Transient Thermal Heating of Insulated Aerospace Vehicle Structures

An analytical solution was derived for the transient response of an insulated structure subjected to a simplified heat pulse. The solution is solely a function of two nondimensional parameters. Simpler functions of these two parameters were developed to approximate the maximum structural temperature over a wide range of parameter values. Techniques were developed to choose constant, effective thermal properties to represent the relevant temperature and pressure-dependent properties for the insulator and structure. A technique was also developed to map a time-varying surface temperature history to an equivalent square heat pulse. Equations were also developed for the minimum mass required to maintain the inner, unheated surface below a specified temperature. In the course of the derivation, two figures of merit were identified. Required insulation masses calculated using the approximate equation were shown to typically agree with finite element results within 10%-20% over the relevant range of parameters studied.

Blosser, Max L.↗
NASA NTRS2017

Pressurization of a Flightweight, Liquid Hydrogen Tank: Evaporation & Condensation at a Liquid/Vapor Interface

This paper presents an analysis and simulation of evaporation and condensation at a motionless liquid/vapor interface. A 1-D model equation, emphasizing heat and mass transfer at the interface, is solved in two ways, and incorporated into a subgrid interface model within a CFD simulation. Simulation predictions are compared with experimental data from the CPST Engineering Design Unit tank, a cryogenic fluid management test tank in 1-g. The numerical challenge here is the physics of the liquid/vapor interface; pressurizing the ullage heats it by several degrees, and sets up an interfacial temperature gradient that transfers heat to the liquid phase-the rate limiting step of condensation is heat conducted through the liquid and vapor. This physics occurs in thin thermal layers O(1 mm) on either side of the interface which is resolved by the subgrid interface model. An accommodation coefficient of 1.0 is used in the simulations which is consistent with theory and measurements. This model is predictive of evaporation/condensation rates, that is, there is no parameter tuning.

Stewart, Mark E. M.↗
NASA NTRS2023

CFD Integrated Transition Modeling for High-Speed Flows via Coupled OVERFLOW-LASTRAC Analysis

This work details ongoing efforts at the NASA Langley Research Center to develop and validate a general-use CFD tool that includes built-in predictions of boundary-layer transition in high-speed flows. Existing tools for the coupling of the NASA OVERFLOW structured overset RANS solver and LASTRAC stability analysis code have been extended to capture boundary-layer transition in high-speed flows driven by either Mack’s first-mode or secondmode instabilities and/or crossflow instabilities. The efficacy of this coupled approach has been demonstrated by examining a variety of supersonic and hypersonic test cases for which experimental validation data is available. The intermittency prescription parameters inherited from prior low-speed applications are found to be suitable for some, but not all high-speed flow scenarios. The method is found to be robust for straight cone configurations and improvements necessary to accurately capture transition on geometries with spatially complex amplification factor envelopes with intermediate regions of slow N-factor variation are examined. The automated, coupled analysis is demonstrated for a 3d supersonic test case and found to perform well within the limits of the linear stability analysis on which it relies. Finally, a preliminary investigation of the method’s robustness to non-ideal CFD meshes is conducted.

High-speed flow↗
NASA NTRSDate not supplied

Composite Bond Line Measurements Based on A Bayesian Analysis of Flash Thermography Data

For bonded composite materials, an accurate characterization of the adhesive bond line is needed to predict failure modes and fracture toughness. In this paper, bond line thickness was estimated from data obtained using through transmission flash thermography. The forward model that predicts back surface temperature is based on a three layer heat diffusion equation with varying diffusivity and flux boundary conditions. The corresponding inverse problem of estimating bond line thickness from measurement data was solved using a Bayesian approach that assumed Gaussian priors for the bond line thickness and thermal diffusivity of the adherends. Finally, the outputs of the thermography based method were compared to measurements that were collected using a micrometer and ultrasound testing.

Flash thermography↗
NASA NTRSDate not supplied

Composite Bond Line Measurements Based on a Bayesian Analysis of Flash Thermography Data

For bonded composite materials, an accurate characterization of the adhesive bond line is needed to predict failure modes and fracture toughness. In this paper, bond line thickness was estimated from data obtained using through transmission flash thermography. The forward model that predicts back surface temperature is based on a three layer heat diffusion equation with varying diffusivity and flux boundary conditions. The corresponding inverse problem of estimating bond line thickness from measurement data was solved using a Bayesian approach that assumed Gaussian priors for the bond line thickness and thermal diffusivity of the adherends. Finally, the outputs of the thermography based method were compared to measurements that were collected using a micrometer and ultrasound testing.

Flash thermography↗