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.
The results are presented of aerodynamic heating tests conducted on a 0.0175-scale model of the space shuttle orbiter. The tests were conducted at a nominal Mach number of 8 and model angles of attack varying between 20 deg and 45 deg.
Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.
Closed cell polyurethane foam for cryogenic insulation, determining thermal conductivity and net heat flow by analytical model and test
The structure and energy balance of the chromosphere-corona transition region is investigated by means of a static, planar model which is compared with the results of XUV-resonance-line observations. In this model, the transition region is heated by thermal conduction from the corona and cooled by radiative losses. Comparison of the model with observational results implies that this is the dominant process in the energy balance of the transition region, and that the base of the transition region is inherently nonstatic and/or nonplanar. The model explains the observational finding of Noyes et al. (1970) that the number density and the downward heat flux both increase by the same factor from quiet regions to active regions.
Spacecraft thermal margins depend on the temperature penalty of installed interfaces, yet catalog conductivity omits bondline thickness, mating surfaces, preload, and environment. The Testbed for Advanced Interface Materials in Vacuum (TAIMV) addresses this engineering-data gap by converting ambient/vacuum temperature fields from a four-coupon stack into quantities usable in spacecraft thermal models. Its 17-node steady-state model resolves axial conduction, parasitic fixture paths, grease-filled contact, radiation, and ambient convection. Four monolithic checks gave axial heat-rate ratios of 0.90–1.02, and Braycote calibration self-recovery gave 2.20% casewise mean absolute percentage error (MAPE) and 2.85% maximum difference. With transferable fixture terms frozen, Krytox LVP tested cross-material transfer against a manufacturer-derived external relation. The six Krytox cases gave 2.28% MAPE, 5.10% maximum difference, and full-field residual root-mean-square values of 0.553 °C ambient and 1.306 °C vacuum. Installed-joint resistance spanned 6.39–7.15 × 10⁻⁴ m²·K/W for Braycote and 4.09–4.49 × 10⁻⁴ m²·K/W for Krytox, corresponding to 2.33–6.66 K per modeled interface over the tested heat-flux range. TAIMV therefore supplies directly usable installed-joint resistance and conductance, plus apparent installed conductivity and explicitly model-conditioned grease conductivity/contact terms.
Heat transfer in high-temperature, high-porosity, flexible refractory ceramic fibrous insulation felts is investigated. Heat transfer in these insulation materials consists of combined gas conduction, solid conduction, and radiation modes, with the precise theoretical modeling of the latter two modes being formidable. A semi-empirical model that requires inverse methods and steady-state thermal test data to infer some of the required model parameters is further developed in this study, and applied to five insulation materials for temperatures between 300 K and 1900 K. The steady-state thermal test setup at NASA Langley Research Center with recent modifications to increase its testing capability to 1900 K is discussed. Design considerations to ensure one-dimensional heat transfer in the test setup are described. Test data and corresponding thermal models for alumina and zirconia-based fibrous insulation felts are presented. Furthermore, test data and thermal models on two fibrous insulation samples containing additives to further suppress either radiation or gas conduction modes of heat transfer are presented. Previously published alumina-based insulation data are also re-processed using the updated modeling methodology. The significance of various heat transfer modes in typical insulation samples is discussed and used to provide general guidance on optimum insulation layups.
Energy-efficient scientific applications require insight into how high performance computing system features impact the applications' power and performance. This insight can result from the development of performance and power models. Here, in this article, we use the modeling and prediction tool MuMMI (Multiple Metrics Modeling Infrastructure) and 10 machine learning methods to model and predict performance and power consumption and compare their prediction error rates. We use an algorithm-based fault-tolerant linear algebra code and a multilevel checkpointing fault-tolerant heat distribution code to conduct our modeling and prediction study on the Cray XC40 Theta and IBM BG/Q Mira at Argonne National Laboratory and the Intel Haswell cluster Shepard at Sandia National Laboratories. Our experimental results show that the prediction error rates in performance and power using MuMMI are less than 10% for most cases. By utilizing the models for runtime, node power, CPU power, and memory power, we identify the most significant performance counters for potential application optimizations, and we predict theoretical outcomes of the optimizations. Based on two collected datasets, we analyze and compare the prediction accuracy in performance and power consumption using MuMMI and 10 machine learning methods.
Membrane distillation (MD) is a thermal desalination process that is gaining attention for treating hypersaline brines. One recent approach uses composite membranes with a thermally conductive layer that delivers heat to the membrane–feed interface, where it drives evaporation. Though this increases single-pass recovery, it comes with the challenge of promoting lateral heat conduction through the thin membrane. Herein, we develop a 3D, computational fluid dynamics (CFD) model that simulates conjugate heat, mass, and momentum transport in the feed channel and composite membrane. We then use the CFD to verify a simpler numerical model that approximates the feed velocity field analytically. We validate the numerical model experimentally, and use it to investigate the impacts of the conductive layer thickness, feed channel geometry, and operating conditions on temperature and concentration polarization, permeate production, and specific energy consumption. Overall, we find that lateral heating increases permeate production at the expense of increased concentration polarization. In extreme cases, the concentration increases more than four-fold along the membrane surface. Furthermore, we show however, that the feed flow rate and conductive layer thickness can be tailored to mitigate concentration polarization, for only a small reduction in permeate production.
Not Available
Design and development of spacecraft radiator using forced convection, conduction, and radiation for heat transfer under steady state conditions
This paper presents a steady-state thermal model of a hot-wire instrument applicable to atmospheric measurement of water content in clouds. In this application, the power required to maintain the wire at a given temperature is used to deduce the water content of the cloud. The model considers electrical resistive heating, axial conduction, convection to the flow, radiation to the surroundings, as well as energy loss due to the heating, melting, and evaporation of impinging liquid and or ice. All of these parameters can be varied axially along the wire. The model further introduces a parameter called the evaporation potential which locally gauges the maximum fraction of incoming water that evaporates. The primary outputs of the model are the steady-state power required to maintain a spatially-average constant temperature as well as the variation of that temperature and other parameters along the wire. The model is used to understand the sensitivity of the hot-wire performance to various flow and boundary conditions including a detailed comparison of dry air and wet (i.e. cloud-on) conditions. The steady-state power values are compared to experimental results from a Science Engineering Associates (SEA) Multi-Element probe, a commonly used water-content measurement instrument. The model results show good agreement with experiment for both dry and cloud-on conditions with liquid water content. For ice, the experimental measurements under read the actual water content due to incomplete evaporation and splashing. Model results, which account for incomplete evaporation, are still higher than experimental results where the discrepancy is attributed to splashing mass-loss which is not accounted in the model.
This paper presents the correlation between the through-thickness thermal diffusivity and the porosity of composites. Finite element analysis (FEA) was used to determine the transient thermal response of composites that were subjected to laser heating. A series of finite element models were built and thermal responses for isotropic and orthographic materials with various thermal diffusivities subjected to different heating conditions were investigated. Experiments were conducted to verify the models and to estimate the unknown parameters such as the amount of heat flux. The analysis and experimental results show good correlation between thermal diffusivity and porosity in the composite materials. They also show that both laser and flash heating can be used effectively to obtain thermal diffusivity. The current infrared thermography system is developed for use with flash heating. The laser heating models and the FEA results can provide useful tools to develop practical thermal diffusivity measurement scheme using laser heat.
A two-fluid model of the solar wind was introduced by Sturrock and Hartle (1966) and Hartle and Sturrock (1968). In these studies the proton energy equation was integrated neglecting the heat conductive term. Later several authors solved the equations for the two-fluid solar wind model keeping the proton heat conductive term. Methods where the equations are integrated simultaneously outward and inward from the critical point were used. The equations were also integrated inward from a large heliocentric distance. These methods have been applied to cases with low coronal base electron densities and high base temperatures. In this paper we present a method of integrating the two-fluid solar wind equations using an iteration procedure where the equations are integrated separately and the proton flux is kept constant during the integrations. The technique is applicable for a wide range of coronal base densities and temperatures. The method is used to carry out a parameter study of the two-fluid solar wind.
The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is an electrically heated, helium cooled, experimental facility. The HTTF was modeled using RELAP-7 in both 1-D and 3-D for validation. Due to core symmetry, 1/6th of the HTTF core was modeled in 3-D using the coupled heat conduction and forced convection capabilities of RELAP-7. This served as a validation case for RELAP-7 and its capabilities to model advanced nuclear technologies such as high temperature gas-cooled reactors.
A computer program implements a thermal model of an insulated wire carrying electric current and surrounded by a vacuum. The model includes the effects of Joule heating, conduction of heat along the wire, and radiation of heat from the outer surface of the insulation on the wire. The model takes account of the temperature dependences of the thermal and electrical properties of the wire, the emissivity of the insulation, and the possibility that not only can temperature vary along the wire but, in addition, the ends of the wire can be thermally grounded at different temperatures. The resulting second-order differential equation for the steady-state temperature as a function of position along the wire is highly nonlinear. The wire is discretized along its length, and the equation is solved numerically by use of an iterative algorithm that utilizes a multidimensional version of the Newton-Raphson method.
The High Temperature Test Facility (HTTF) is an electrically heated, helium cooled, experimental facility located at Oregon State University (OSU). Modeling of the HTTF PG-26 transient using a multi-app approach utilizing MOOSE and RELAP-7 was conducted. A 3-D model of the HTTF core, reflector, core barrel, and RPV was generated, and 3-D heat conduction throughout the structure was coupled to 1-D fluid flow results from RELAP-7. The model was able to accurately predict peak ceramic core temperatures experienced during the transient which helps validate the capability of RELAP-7 to model advanced gas-cooled nuclear reactors.
Convective and radiative heating environments have been computed for a three-dimensional ellipsled configuration which would perform an aerocapture maneuver at Neptune. This work was performed as part of a one-year Neptune aerocapture spacecraft systems study that also included analyses of trajectories, atmospheric modeling, aerodynamics, structural design, and other disciplines. Complementary heating analyses were conducted by separate teams using independent sets of aerothermodynamic modeling tools (i.e. Navier-Stokes and radiation transport codes). Environments were generated for a large 5.50 m length ellipsled and a small 2.88 m length ellipsled. Radiative heating was found to contribute up to 80% of the total heating rate at the ellipsled nose depending on the trajectory point. Good agreement between convective heating predictions from the two Navier-Stokes solvers was obtained. However, the radiation analysis revealed several uncertainties in the computational models employed in both sets of codes, as well as large differences between the predicted radiative heating rates.
Abstract Salt offers an optimal medium for the permanent isolation of heat‐producing radioactive waste due to its impermeability, high thermal conductivity, and ability to close fractures through creep. A thorough understanding of the thermal‐hydrological‐mechanical (THM) processes, encompassing brine migration, is fundamental for secure radioactive waste disposal within salt formations. At the Waste Isolation Pilot Plant (WIPP), we conducted joint in situ geophysical monitoring experiments during active heating to investigate brine migration near excavations. This experiment incorporated electrical resistivity tomography (ERT) alongside high‐resolution fiber‐optic‐based distributed temperature sensing within a controlled heating experiment. Additionally, discrete element model (DEM) based numerical simulations were conducted to simulate THM processes during heating, providing a more mechanistic understanding of the coupled processes leading to the observed changes in geophysical measurements. During heating, resistivity shifts near the heater were reasonably explained by temperature effects. However, in more distant, cooler regions, the resistivity decrease exceeded predictions based solely on temperature. DEM simulations highlighted brine migration, propelled by pore pressure gradients, as the likely primary factor contributing to the additional resistivity decline beyond temperature effects. The comparison between the predicted ERT responses and observations was much improved when considering the effects of brine migration based on the DEM simulations. These geophysical and simulation findings shed light on brine migration in response to salt heating, enhancing our understanding of the coupled THM processes in salt for safe radioactive waste disposal.