SEARCH · Engineering Papers
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.
Thermal conductivity of closed cell polyurethane foam insulation, analyzing steady state heat flow, gas diffusion, condensing, mole fraction and temperature
The outflow of coronal plasma into interplanetary space is a consequence of the coronal heating process. Therefore the formation of the corona and the acceleration of the solar wind should be treated as a single problem. Traditionally the mass or particle flux emanating from the extended corona has been thought of as being determined by the coronal temperature or scale height and the coronal (base) density. This argument follows from considerations of the momentum balance of the corona-wind system from which one obtains models of a close to hydrostatic corona out to the critical point where the flow becomes supersonic. With this approach to the acceleration of the wind is has been difficult to reconcile the relatively small variation observed in the proton flux at 1 AU with the predicted exponential dependence of the proton flux on the coronal temperature. In this talk we would like to emphasize another approach in which coronal energetics play the primary role. The deposition of energy into the corona through some 'mechanical' energy flux is balanced by the various energy sinks available to the corona and the sum of these processes determine the coronal structure, i.e. its temperature and density. The corona loses energy through heat conduction into the transition region, through radiative losses, and through the gravitational potential energy and kinetic energy put into the solar wind itself. We will show from a series of models of the chromosphere transition region-corona-solar wind system that most of the energy deposited in a magnetically open region will go into the solar wind, with roughly half going into kinetic energy and half into lifting the plasma out of the solar gravity field. The coronal base density will adjust itself in such a way that the heat conductive flux flowing into the transition region is radiated away in the upper chromosphere. The coronal temperature is set by the requirements that most of the deposited energy goes into accelerating the solar wind; the coronal scale height will adjust itself so that the solar wind energy losses conform to the amplitude of the input energy. These processes are modified by the 'mode' of energy deposition, and we will show the effects on coronal structure of changing the parameters describing coronal heating as well as the effects of including a helium fluid in the models. However, the location, scale height and/or form of the energy deposition (i.e. heating or direct acceleration) are not too important for the solar wind, the coronal density and temperature structure will vary with the 'mode' of energy deposition, but the solar wind mass flux depends mainly on the amplitude of the energy flux.
A new method to model vibrational power flow and predict the resulting energy density levels in uniform rods and beams is investigated. This method models the flow of vibrational power in a manner analogous to the flow of thermal power in a heat conduction problem. The classical displacement solutions for harmonically excited, hysteretically damped rods and beams are used to derive expressions for the vibrational power flow and energy density in the rod and beam. Under certain conditions, the power flow in these two structural elements will be shown to be proportional to the energy density gradient. Using the relationship between power flow and energy density, an energy balance on differential control volumes in the rod and beam leads to a Poisson's equation which models the energy density distribution in the rod and beam. Coupling the energy density and power flow solutions for rods and beams is also discussed. It is shown that the resonant behavior of finite structures complicates the coupling of solutions, especially when the excitations are single frequency inputs. Two coupling formulations are discussed, the first based on the receptance method, and the second on the travelling wave approach used in Statistical Energy Analysis. The receptance method is the more computationally intensive but is capable of analyzing single frequency excitation cases. The traveling wave approach gives a good approximation of the frequency average of energy density and power flow in coupled systems, and thus, is an efficient technique for use with broadband frequency excitation.
Traditional Air Plasma Sprayed (APS) ZrO2-Y2O3 Thermal Barrier Coatings (TBC's) and Low Pressure Plasma Sprayed (LPPS) ZrO2-Y2O3/Ni-Cr-Al-Y cermet coatings were tested in a H2/O2 rocked engine. The traditional ZrO2-Y2O3 (TBC's) showed considerable metal temperature reductions during testing in the hydrogen-rich environment. A thermal model was developed to predict the thermal response of the tubes with the various coatings. Good agreement was observed between predicted temperatures and measured temperatures at the inner wall of the tube and in the metal near the coating/metal interface. The thermal model was also used to examine the effect of the differences in the reported values of the thermal conductivity of plasma sprayed ZrO2-Y2O3 ceramic coatings, the effect of 100 micron (0.004 in.) thick metallic bond coat, the effect of tangential heat transfer around the tube, and the effect or radiation from the surface of the ceramic coating. It was shown that for the short duration testing in the rocket engine, the most important of these considerations was the effect of the uncertainty in the thermal conductivity of temperatures (greater than 100 C) predicted in the tube. The thermal model was also used to predict the thermal response of the coated rod in order to quantify the difference in the metal temperatures between the two substrate geometries and to explain the previously-observed increased life of coatings on rods over that on tubes. A thermal model was also developed to predict heat transfer to the leading edge of High Pressure Fuel Turbopump (HPFTP) blades during start-up of the space shuttle main engines. The ability of various TBC's to reduce metal temperatures during the two thermal excursions occurring on start-up was predicted. Temperature reductions of 150 to 470 C were predicted for 165 micron (0.0065 in.) coatings for the greater of the two thermal excursions.
Previous data and high and low dispersion IUE observations of Chi 1 Ori, Alpha Cen A, Xi Boo A, Alpha Cen B, and Epsilon Eri are used to model the atmospheres of these stars from the high chromosphere to the corona. The energy lost by radiation and transferred by thermal conduction is investigated to establish the heating requirements, and results similar to those found for the solar atmosphere are found. The observed emission measurement distribution below 100,000 K can be described by a model in which Alfven wave energy input, observed through nonthermal line broadening, is balanced by radiation losses. Comparison of the coronal, transition region, and chromospheric pressures shows scaling relations which are compatible with previous flux correlations.
Dr. Albert's current research is centered on transfer processes in porous media, including air-snow exchange in the Polar Regions and in soils in temperate areas. Her research includes field measurements, laboratory experiments, and theoretical modeling. Mary conducts field and laboratory measurements of the physical properties of natural terrain surfaces, including permeability, microstructure, and thermal conductivity. Mary uses the measurements to examine the processes of diffusion and advection of heat, mass, and chemical transport through snow and other porous media. She has developed numerical models for investigation of a variety of problems, from interstitial transport to freezing of flowing liquids. These models include a two-dimensional finite element code for air flow with heat, water vapor, and chemical transport in porous media, several multidimensional codes for diffusive transfer, as well as a computational fluid dynamics code for analysis of turbulent water flow in moving-boundary phase change problems.
The lowest surface temperature possible for the existance of spray evaporative cooling is determined experimentally to be a linear function of the impinging spray mass flux. A conduction-controlled analytical model of droplet evaporation gives fairly good agreement with experimental measurements at atmospheric pressure. At reduced pressures droplet evaporation rates are decreased significantly such that an optimum operating pressure exists for each desired surface heat flux. The initiation of the 'Leidenfrost state' provides the upper surface temperature bound for spray evaporative cooling.
A laboratory heating test simulating hypersonic heating was conducted on a heat-sink type structure to provide basic thermal stress measurements. Six NASTRAN models utilizing various combinations of bar, shear panel, membrane, and plate elements were used to develop calculated thermal stresses. Thermal stresses were also calculated using a beam model. For a given temperature distribution there was very little variation in NASTRAN calculated thermal stresses when element types were interchanged for a given grid system. Thermal stresses calculated for the beam model compared similarly to the values obtained for the NASTRAN models. Calculated thermal stresses compared generally well to laboratory measured thermal stresses. A discrepancy of signifiance occurred between the measured and predicted thermal stresses in the skin areas. A minor anomaly in the laboratory skin heating uniformity resulted in inadequate temperature input data for the structural models.
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.
Knowledge of the thermal conductivity of boron carbides is necessary to evaluate its potential for high temperature thermoelectric energy conversion applications. The thermal diffusivity of hot pressed boron carbide B/sub 1-x/C/sub x/ samples as a function of composition, temperature and temperature cycling was measured. These data in concert with density and specific heat data yield the thermal conductivities of these materials. The results in terms of a structural model to explain the electrical transport data and novel mechanisms for thermal conduction are discussed.
Presumed end-of-life atmospheric reentry of the GA SP-100 system was studied to assess dispersal feasibility and associated hazards. Reentry was studied by sequential use of an orbital trajectory and a heat analysis computer program. Two heating models were used. The first model assumed a thermal equilibrium condition between the stagnation point aerodynamic heating and the radiative cooling of the skin material surface. The second model allowed for infinite conductivity of the skin material. Four reentering configurations were studied representing stages of increased SP-100 breakup: (1) radiator, shield and reactor, (2) shield and reactor, (3) reactor with control drums, and (4) reactor without control drums. Each reentering configuration was started from a circular orbit at 116 km having an inertial velocity near Mach 25. The assumed failing criterion was the attainment of melting temperature of a critical system component. The reentry analysis revealed breakup of the vessel in the neighborhood of 61 km altitude and scattering of the fuel elements. Subsequent breakup of the fuel elements was not predicted. Oxidation of the niobium skin material was calculated to cause an increase in surface temperature of less than ten percent. The concept of thermite analogs for enhancing reactor reentry dispersal was assessed and found to be feasible in principle. A conservative worst-case hazards analysis was performed for radioactive and nonradioactive toxic SP-100 materials assumed to be dispersed during end-of-life reentry. The hazards associated with this phase of the SP-100 mission were calculated to be insignificant.
The dynamic formation of prominences in the initial magnetothermal equilibrium and their stability to sideward displacements are investigated focusing on the structure of the 2D solar atmosphere in the presence of coronal arcades or loops. A model based on 2D magnetohydrodynamic equations takes into account gravity, compressible flows, heating, radiation, anisotropic thermal conduction, and coupling to a deep chromosphere. It is found that prominences in simple arcades characterized by magnetic field with significant curvature at the apex are unstable to a lateral displacement.
A number of wire mesh and metal felt test samples, with a range of porosities, yield generic correlations for friction factor, Nusselt number, enhanced axial conduction ratio, and overall heat flux ratio. This information is directed primarily toward stirling cycle regenerator modelers, but will be of use to anyone seeking to better model fluid flow through these porous materials. Behind these results lies an oscillating-flow test rig, which measures pumping dissipation and thermal energy transport in sample matrices, and several stages of data-reduction software, which correlate instantaneous values for the above dimensionless groups. Within the software, theoretical model reduces instantaneous quantifies from cycle-averaged measurables using standard parameter estimation techniques.
A flight investigation has been conducted to study the heat transfer to swept-wing leading edges. A rocket-powered model was used for the investigation and provided data for Mach number ranges of 1.78 to 2.99 and 2.50 to 4.05 with corresponding free-stream Reynolds number per foot ranges of 13.32 x 10(exp 6) to 19.90 x 10(exp 6) and 2.85 x 10(exp 6) to 4.55 x 10(exp 6). The leading edges employed were cylindrically blunted wedges ', three of which were swept 450 with leading-edge diameters of 1/4, 1/2, and 3/4 inch and one swept 36-750 with a leading-edge diameter of 1/2 inch. In the high Reynolds number range, measured values of heat transfer were found to be much higher than those predicted by laminar theory and at the larger values of leading-edge diameter were approaching the values predicted by turbulent theory. For the low Reynolds number range a comparison between measured and theoretical heat transfer showed that increasing the leading-edge diameter resulted in turbulent flow on the cylindrical portion of the leading edge.
We report on our progress towards the development of arrays of X-ray microcalorimeters as candidates for the high resolution x-ray spectrometer on the Constellation-X mission. The microcalorimeter arrays (30 x 30) with appropriate pixel sizes (0.25 mm. x 0.25 mm) and high packing fractions (greater than 96%) are being developed. Each individual pixel has a 10 micron thick Bi X-ray absorber that is shaped like a mushroom to increase the packing fraction, and a Mo/Au proximity effect superconducting transition edge sensor (TES). These are deposited on a 0.25 or 0.5 micron thick silicon nitride membrane with slits to provide a controllable weak thermal link to the sink temperature. Studies are underway to model, test and optimize the TES pixel uniformity, critical current, heat capacity and the membrane thermal conductance in the array structure. Fabrication issues and procedures, and results of our efforts based on these optimizations will be provided.
A unique multilayer insulation concept for orbital cryogenic storage was experimentally evaluated at NASA Marshall Space Flight Center (MSFC) using the Multipurpose Hydrogen Test Bed (MHTB). A combination of foam/Multi layer Insulation (MLI) was used. The MLI (45 layers of Double Aluminized Mylar (DAM) with Dacron net spacers) was designed for an on-orbit storage period of 45 days and included several unique features such as: a variable layer density and larger but fewer DAM perforations for venting during ascent to orbit. The focus of this paper is on analytical modeling of the variable density MLI performance during orbital coast periods. The foam/MLI combination model is considered to have five segments. The first segment represents the foam layer. The second, third, and fourth segments represent the three layers of MLI with different layer densities and number of shields. Finally, the last segment is considered to be a shroud that surrounds the last MLI layer. The hot boundary temperature is allowed to vary from 164 K to 305 K. To simulate MLI performance, two approaches are considered. In the first approach, the variable density MLI is modeled layer by layer while in the second approach, a semi-empirical model is applied. Both models account for thermal radiation between shields, gas conduction, and solid conduction through the separator materials. The heat flux values predicted by each approach are compared for different boundary temperatures and MLI systems with 30, 45, 60, and 75 layers.
An experimental program was generated by the Icing Technology Branch at NASA Glenn Research Center to validate two ice protection simulation codes: (1) LEWICE/Thermal for transient electrothermal de-icing and anti-icing simulations, and (2) ANTICE for steady state hot gas and electrothermal anti-icing simulations. An electrothermal ice protection system was designed and constructed integral to a 36 inch chord NACA0012 airfoil. The model was fully instrumented with thermo-couples, RTD'S, and heat flux gages. Tests were conducted at several icing environmental conditions during a two week period at the NASA Glenn Icing Research Tunnel. Experimental results of running-wet and evaporative cases were compared to the ANTICE computer code predictions and are presented in this paper.
We introduce a new framework for Machine-Learning (ML) based Auroral Ionosphere Model (ML-AIM). ML-AIM solves a current continuity equation by utilizing the ML model of Field Aligned Currents (FACs) of Kunduri et al., 2020 (https://doi.org/10.1029/2020JA027908), the FAC-derived aurora conductance model of Robinson et al., 2020 (https://doi.org/10.1029/2020JA028008), and the solar irradiance conductance model of Moen & Brekke (1993). The ML-AIM inputs are 60min time histories of solar wind plasma, interplanetary magnetic fields (IMF), and geomagnetic indices, and its outputs are ionospheric electric potential, electric fields, Pederson/Hall currents, and Joule Heating. We conduct two ML-AIM simulations for a weak geomagnetic activity on 14 May 2013 and a geomagnetic storm on 7-8 September 2017. ML-AIM produces reasonable ionospheric potential patterns such as two cell convection patterns and the enhancement of electric potentials during active times. The cross polar cap potential drop from ML-AIM is also comparable to the ones from the Weimer 2005 model, Super Dual Auroral Radar Network (SuperDARN), and Defense Meteorological Satellite Program (DMSP) F17 observations. ML-AIM is unique in a sense that it predicts ionospheric responses to the time-varying solar wind and geomagnetic conditions, while other traditional empirical model like Weimer 2005 is designed to provide static ionospheric conditions under steady solar wind/IMF conditions. In future, ML-AIM will include ML-based models of aurora precipitation and ionospheric conductance, improving its performance during active times.