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 631 records · Page 35

Afterbody Heating Characteristics of a Proposed Mars Sample Return Orbiter

Aeroheating wind-tunnel tests were conducted on a 0.028 scale model of an orbiter concept considered for a possible Mars sample return mission. The primary experimental objectives were to characterize hypersonic near wake closure and determine if shear layer impingement would occur on the proposed orbiter afterbody at incidence angles necessary for a Martian aerocapture maneuver. Global heat transfer mappings, surface streamline patterns, and shock shapes were obtained in the NASA Langley 20-inch Mach 6 Air and CF4 Tunnels for post-normal shock Reynolds numbers (based on forebody diameter) ranging from 1,400 to 415,000, angles of attack ranging from -5 to 10 degrees at 0, 3, and 6 deg sideslip, and normal-shock density, ratios of 5 and 12. Laminar, transitional, and turbulent shear layer impingement on the cylindrical afterbody was inferred from the measurements and resulted in a localized heating maximum that ranged from 40 to 75% of the reference forebody stagnation point heating. Comparison of laminar heating prediction to experimental measurement along the orbiter afterbody highlight grid alignment challenges associated with numerical simulation of three-dimensional separated wake flows.

Horvath, Thomas J.↗

Consequences of a conductive boundary on the local cloud. I - No dust

Using boundary conditions appropriate to the Local Cloud embedded in the Local Bubble, the evaporation of the cloud is modeled assuming steady flow. Included in the model are the effects of radiative cooling and nonequilibrium ionization structure in the interface, and pressure support and heat conduction inhibition due to the magnetic field. The soft X-ray and EUV emission from the interface are calculated and compared with observations. It is found that the emission in the softest X-ray band, the Be band, is harder than that calculated for an equilibrium, normal abundance, 10 to the 6th K plasma. Column densities of C IV, N V, O VI, Si IV, and Si III in the interface are calculated. The column density for O VI is found to be very close to observed values for nearby stars, while all others are too small to contribute significantly to observed values. The ionization of the Local Cloud due to radiation from the interface is calculated. It is found that a significant amount of ionizing radiation is produced, dominating all other known sources. The hardness of the emission, however, results in ionization of helium as well as hydrogen. Thus, the type of boundary modeled cannot explain the ionization of the very local interstellar medium inferred from solar backscatter measurements.

Slavin, Jonathan D.↗

MOCVD of GaAs in a horizontal reactor - Modeling and growth

A two-dimensional model for metalorganic chemical vapor deposition of GaAs in a horizontal reactor is presented. The model is characterized by the following parameters: reactor geometry and operating pressure, thermal boundary conditions, ratio of reactants, chemical reactions, total inlet gas flow rate, as well as molecular weights, thermal conductivities, heat capacities, viscosities, and binary diffusion coefficients of the gas-phase species. Film thickness profiles predicted by the model are compared with those of GaAs thin films grown in the modeled reactor. Results obtained show a good agreement between the predictions and data over the entire length of the deposition region for the low pressure and high flow rate run. Attention is also given to the reactor design and growth conditions.

Clark, Ivan O.↗

A Model for Hydrogen Thermal Conductivity and Viscosity Including the Critical Point

In order to conduct a thermal analysis of heat transfer to liquid hydrogen near the critical point, an accurate understanding of the thermal transport properties is required. A review of the available literature on hydrogen transport properties identified a lack of useful equations to predict the thermal conductivity and viscosity of liquid hydrogen. The tables published by the National Bureau of Standards were used to perform a series of curve fits to generate the needed correlation equations. These equations give the thermal conductivity and viscosity of hydrogen below 100 K. They agree with the published NBS tables, with less than a 1.5 percent error for temperatures below 100 K and pressures from the triple point to 1000 KPa. These equations also capture the divergence in the thermal conductivity at the critical point

Wagner, Howard A.↗

Matrix Graphite Material Models In Pebbles and Compacts For Bison

The cores and reflectors in high-temperature gas-cooled reactors (HTGRs) are made of graphite materials, with the graphite acting as a moderator, a fuel host matrix, or the foundation for various structural components. This study aims to survey the models in the literature for graphite materials being used as host matrices in pebble/fuel compacts and to implement those surveyed models into Bison to conduct an early assessment of graphite's thermo-mechanical response under various reactor conditions. In this study, thermal (e.g., thermal conductivity, and specific heat capacity) and mechanical (e.g., elastic properties, thermal expansion, irradiation-induced dimensional changes, and irradiation-induced creep) material models for various graphite grades (e.g., H-451, IG-110, G-348, 2020, A3-3, and A3-27) are incorporated into Bison. Two benchmark problems are then exercised utilizing these new graphite-related capabilities: (1) modeling an Advanced Gas Reactor (AGR)-2 fuel compact, and (2) modeling the debonding of a particle-matrix interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Transient Analysis of Thermal Protection System for X-33 Aircraft using MSC/NASTRAN

X-33 is an advanced technology demonstrator vehicle for the Reusable Launch Vehicle (RLV) program. The thermal protection system (TPS) for the X-33 is composed of complex layers of materials to protect internal components, while withstanding severe external temperatures induced by aerodynamic heating during high speed flight. It also serves as the vehicle aeroshell in some regions using a stand-off design. MSC/NASTRAN thermal analysis capability was used to predict transient temperature distribution (within the TPS) throughout a mission, from launch through the cool-off period after landing. In this paper, a typical analysis model, representing a point on the vehicle where the liquid oxygen tank is closest to the outer mold line, is described. The maximum temperature difference between the outer mold line and the internal surface of the liquid oxygen tank can exceed 1500 F. One dimensional thermal models are used to select the materials and determine the thickness of each layer for minimum weight while insuring that all materials remain within the allowable temperature range. The purpose of working with three dimensional (3D) comprehensive models using MSC/NASTRAN is to assess the 3D radiation effects and the thermal conduction heat shorts of the support fixtures.

Miura, Hirokazu↗

Heating of the nighttime D region by very low frequency transmitters

VLF signals propagating in the Earth-ionosphere waveguide are used to probe the heated nighttime D region over three U.S. Navy very low frequency (VLF,3-30 kHz) transmitters. Ionospheric cooling and heating are observed when a transmitter turns off and on in the course of normal operations. Heating by the 24.0-kHz NAA transmitter in Cutler, Maine, (1000 kW radiated power) was observed by this method in 41 of 52 off/on episodes during December 1992, increasing the amplitude and retarding the phase of the 21.4-kHz NSS probe wave propagating from Annapolis, Maryland, to Gander, Newfoundland, by as much as 0.84 dB and 5.3 deg, respectively. In 6 of these 41 episodes, the amplitude of the 28.5-kHz NAU probe wave propagating from Puerto Rico to Gander was also perturbed by as much as 0.29 dB. The latter observations were unexpected due to the greater than 770 km distance between NAA and the NAU-Gander great circle path. Heating by the NSS (21.4 kHz, 265 kW) and NLK (24.8 kHz, 850 kW) transmitters was observed serendipitously in data from earlier measurements of the amplitudes of VLF signals propagating in the Earth-ionosphere waveguide. A three-dimensional model of wave absorption and electron heating in a magnetized, weakly ionized plasma is used to calculate the extent nad shape of the collision frequency (i.e., electron temperature) enhancement above a VLF transmitter. The enhancements are annular, with a geomagnetic north-south asymmetry and a radius at the outer half-maximum of the collision frequency enhancement of about 150 km. Heating by the NAA transmitter is predicted to increase the nighttime D region electron temperature by as much as a factor of 3. The calculated changes in the D region conductivity are used in a three-dimensional model of propagation in the Earth-ionosphere wavelength to predict the effect of the heated patch on a subionospheric VLF probe wave. The range of predicted scattered field amplitudes is in general consistent with the observed signal perturbations. Discrepanices in the predictions are attributed to lack of knowledge of the D region electron density profile along the probe wave great circle paths.

Rodriguez, Juan V.↗

Optimizing Conductivity and Hardness in Cu–Nb Nanolamellar Composites Fabricated Through Accumulative Roll Bonding Without Intermittent Heat Treatments

High strength/high conductivity Cu–XNb nanolamellar bimetallic composites (where X < 50 pct) were successfully processed via accumulative roll bonding. No intermittent heat treatments were used which increases industrial applicability. Estimates of conductivity from published models were validated at length scales of ~ 150 nm but were shown to be four times too low for composites with layer thicknesses < 40 nm. Nanoscale composites with conductivity of ~ 47 pct IACS and 99 Rockwell B hardness were fabricated.

36 MATERIALS SCIENCE↗

Design of a prototypical natural circulation water-based reactor cavity cooling system (RCCS) for a pebble-bed generic FHR

A prototypical natural circulation water-based reactor cavity cooling system (RCCS) for the UC Berkeley Mark-1 pebble-bed gFHR is designed based on one-dimensional thermal hydraulics modeling and optimization implemented in an in-house MATLAB code. The model employs a lumped core consisting of fuel pebble and graphite pebble regions, while the graphite reflector, vessel, and RCCS are represented as separate but energy-coupled regions. The model is derived based on steady state energy balance equations accounting for conductive heat transfer from the lumped core to the vessel and subsequent radiative heat transfer to the RCCS and convective heat transfer in the water. Further, mass flow rate is calculated based on momentum and thermal energy balance in the RCCS. Starting with ANL’s water NSTF as baseline, effects of RCCS design parameters such as size of the plate between pipes, pipe diameter, source-to-sink distance, surface emissivity, and pitch from the core are examined in an integral effects framework which accounts for subsequent variation of the mass flow rate, number of pipes, convection coefficient, mean region and surface temperatures, radiative view factors, and water outlet temperature. The importance of including a prototypical reactor design in RCCS design calculations is emphasized as the calculations show that temperature drops from the core to the RCCS fluid to transfer the heat tend to be the limitation not the heat removal capacity of the fluid. A closed-loop design is obtained with capability to safely remove up to 0.72% of the nominal reactor power at maximum estimated peak conditions and 0.40% at shutdown. The design is based on physics calculations and does not account for economic optimization.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Contribution of the Hydroxyl Group on Interfacial Heat Conduction of Monohydric Alcohols: A Molecular Dynamics Study

Monohydric alcohols have been used as promising phase change materials (PCMs) for low-temperature latent heat storage. However, the heat storage/retrieval rates are limited due to the low thermal conductivity of such alcohols. In this work, nonequilibrium molecular dynamics (NEMD) simulations were performed to study the microscopic heat conduction in example monohydric alcohols, i.e., 1-dodecanol (C12H26O), 1-tetradecanol (C14H30O), and 1-hexadecanol (C16H34O). A simplified ideal crystal model was proposed to exploit the potential for improving the thermal conductivity of monohydric alcohols. The effect of ideal crystalline structures, especially the contribution of the hydroxyl group, on the microscopic heat conduction process was analyzed. The thermal conductivity of the ideal crystals of the various monohydric alcohols was predicted to be more than twice as compared to that of their respective solids. The major thermal resistance in the ideal crystals was found around the molecular interfaces, as a result of the excellent heat conduction performance along the linear molecular chains. The calculated vibrational density of states (VDOS) and interfacial heat transfer were then investigated. When the interfaces are surrounded by hydroxyl groups as walls, strong hydrogen bond (HB) interactions were observed. The interfacial heat transfer coefficient of the ideal crystalline structures of 1-tetradecanol was found to reach up to ~735.6 MW/m2 W. It was elucidated that the high interfacial heat transfer rate is clearly related to the stronger intermolecular interactions.

Engineering↗

Convection sensitivity and thermal analyses for indium and indium-lead mixing experiment (74-18)

Sounding rocket Experiment 74-18 was designed to demonstrate the effects of the Black Brandt rocket acceleration levels (during the low-g coast phase of its flight) on the motion of a liquid metal system to assist in preflight design. Some post flight analyses were also conducted. Preflight studies consisted of heat transfer analysis and convection sensitivity and convection modeling analyses which aided in the: (1) final selection of fluid materials (indium-lead melts rather than paraffins); (2) design and timing of heater and quench system; and (3) preflight predictions of the degree of lead penetration into the pure indium segment of the fluid. Postflight studies involved: (1) updating the convection sensitivity calculations by utilizing actual flight gravity levels; and (2) modeling the mixing in the flight samples.

Bourgeois, S. V.↗

Shuttle antenna radome technology test program. Volume 2: Development of S-band antenna interface design

The effects of the Thermal Protection Subsystem (TPS) contamination on the space shuttle orbiter S band quad antenna due to multiple mission buildup are discussed. A test fixture was designed, fabricated and exposed to ten cycles of simulated ground and flight environments. Radiation pattern and impedance tests were performed to measure the effects of the contaminates. The degradation in antenna performance was attributed to the silicone waterproofing in the TPS tiles rather than exposure to the contaminating sources used in the test program. Validation of the accuracy of an analytical thermal model is discussed. Thermal vacuum tests with a test fixture and a representative S band quad antenna were conducted to evaluate the predictions of the analytical thermal model for two orbital heating conditions and entry from each orbit. The results show that the accuracy of predicting the test fixture thermal responses is largely dependent on the ability to define the boundary and ambient conditions. When the test conditions were accurately included in the analytical model, the predictions were in excellent agreement with measurements.

Kuhlman, E. A.↗

Convective thinning of the lithosphere - A mechanism for the initiation of continental rifting

A model of lithospheric thinning, in which heat is convected to the base and conducted within the lithosphere, is presented. An analytical equation for determinining the amount of thinning attainable on increasing the heat flux from the asthenosphere is derived, and a formula for lithosphere thickness approximations as a function of time is given. Initial and final equilibrium thicknesses, thermal diffusivity, transition temperature profile, and plume temperature profile are all factors considered for performing rate of thinning determinations. In addition, between initial and final equilibrium states, lithospheric thinning occurs at a rate which is inversely proportional to the square root of the time. Finally, uplift resulting from thermal expansion upon lithospheric thinning is on the order of 10 to the 2nd to 10 to the 3rd m.

Spohn, T.↗

Analytical modeling of helium II in forced flow conditions

A numerical model was developed to examine the flow behavior of He II under forced flow conditions. The model basically solves the governing equations in the two-fluid model of He II, and can be used to assess the pump requirement for transferring He II in spaceborne systems. Predictions of liquid temperature and pressure drop agree well with empirical data obtained from small scale tests. At high transfer rates, the model predicts that the forced convection dominates the Gorter-Mellink conduction in heat transfer, and the classical turbulent flow correlation is sufficient for calculating the pressure drop across the transfer line.

Ng, Y. S.↗

Diurnal variations on a plasmaspheric flux tube - Light ion flows and F region temperature enhancements

The paper concentrates on the diurnal variations on a plasmaspheric flux tube modeled using a time-dependent multispecies one-stream interhemispheric model for plasma flows. The model takes into account the effects of ionization, charge exchange, recombination, collisions, heat conduction, and allows for external heat sources. The simulation is done for June solstice conditions during solar minimum. Focus is placed on the presence of large downward H(+) velocities at about 320-km altitude in the winter (southern) hemisphere, in early morning when the summer hemisphere is sunlit but the winter hemisphere is dark. In addition, an upward H(+) flux is seen in the Southern Hemisphere at altitudes above 2000 km when the sun rises in the northern end.

Guiter, S. M.↗

Heat transfer in rotating serpentine passages with smooth walls

Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multi-pass, smooth-wall heat transfer model with both radially inward and outward flow. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages (coolant-to-wall temperature ratio, Rossby number, Reynolds number and radius-to-passage hydraulic diameter ratio). These four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. It was found that both Coriolis and buoyancy effects must be considered in turbine blade cooling designs and that the effect of rotation on the heat transfer coefficients was markedly different depending on the flow direction. Local heat transfer coefficients were found to decrease by as much as 60 percent and increase by 250 percent from no rotation levels. Comparisons with a pioneering stationary vertical tube buoyancy experiment showed reasonably good agreement. Correlation of the data is achieved employing dimensionless parameters derived from the governing flow equations.

Wagner, J. H.↗

One-Dimensional Hybrid Satellite Track Model for the Dynamics Explorer 2 (DE 2) Satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. satellite track code but also includes elements of Roble et al. global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the MSIS-90 model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule beating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗