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 379 records · Page 21

Community Geothermal: Soil Conductivity, Borehole Design, Energy Models, and Load Data for a Residential System Development - Hinesburg, VT

This dataset contains materials from the Coalition for Community-Supported Affordable Geothermal Energy Systems (C2SAGES) project, which evaluated the techno-economic feasibility of a community geothermal system for a residential development in Hinesburg, VT. The dataset includes detailed soil conductivity test reports, energy models, borehole design reports, hourly energy loads for heating, cooling, and hot water, and design layouts. EnergyPlus was used to model building energy loads, and Modelica software was applied for geothermal loop sizing based on these loads and soil conductivity results. Python scripts for network design further refined the models. Key files include PDF reports on borehole design (with projections for 1-year, 15-year, and 30-year systems), soil conductivity test results, EnergyPlus modeling outputs, and 2D/3D design drawings in PDF, DWG, and DXF formats. Python notebooks for network design and OnePipe model files are also provided, with Modelica required for viewing certain files. Outputs and modeling data are in various formats including CSV, JPG, HTML, and IDF, with units and data clearly labeled to support understanding of system design and performance for the proposed geothermal solution.

15 GEOTHERMAL ENERGY↗

Thermal transport in defective and disordered materials

With significant recent advancements in thermal sciences—such as the development of new theoretical and experimental techniques, and the discovery of new transport mechanisms—it is helpful to revisit the fundamentals of vibrational heat conduction to formulate an updated and informed physical understanding. The increasing maturity of simulation and modeling methods sparks the desire to leverage these techniques to rapidly improve and develop technology through digital engineering and multi-scale, electro-thermal models. With that vision in mind, this review attempts to build a holistic understanding of thermal transport by focusing on the often unaddressed relationships between subfields, which can be critical for multi-scale modeling approaches. For example, we outline the relationship between mode-specific (computational) and spectral (analytical) models. We relate thermal boundary resistance models based on perturbation approaches and classic transmissivity based models. We discuss the relationship between lattice dynamics and molecular dynamics approaches along with two-channel transport frameworks that have emerged recently and that connect crystal-like and amorphous-like heat conduction. Throughout, we discuss best practices for modeling experimental data and outline how these models can guide material-level and system-level design.

42 ENGINEERING↗

Heat transfer tests of an 0.006-scale thin-skin space shuttle thermocouple model (41-OTS) in the Langley Research Center unitary plan wind tunnel at M equals 3.7 (IH16)

The results are presented of supersonic heat transfer tests performed on the .006 scale space shuttle vehicle model (41-OTS) in the Langley Research Center Unitary Plan Wind Tunnel. These tests were conducted to parametrically investigate ascent heating of the integrated vehicle and its components. The tests were conducted at a nominal Mach number of 3.7 and Reynolds numbers per foot of 2 and 5 million. The model configurations investigated were the integrated vehicle and each component alone (i.e. orbiter, tank and SRB). All the configurations were run with and without transition strips and through an angle of attack range of 0 deg to minus 5 deg with the exception of the SRB which was tested through an angle of attack range of minus 5 deg to 90 deg. The heat transfer data were obtained from 223 iron constantan thermocouples attached to stainless steel thin-skin areas of the model.

Walstad, D. G.↗

Active Control of Interface Shape During the Crystal Growth of Lead Bromide

A thermal model for predicting and designing the furnace temperature profile was developed and used for the crystal growth of lead bromide. The model gives the ampoule temperature as a function of the furnace temperature, thermal conductivity, heat transfer coefficients, and ampoule dimensions as variable parameters. Crystal interface curvature was derived from the model and it was compared with the predicted curvature for a particular furnace temperature and growth parameters. Large crystals of lead bromide were grown and it was observed that interface shape was in agreement with the shape predicted by this model.

Duval, W. M. B.↗

Ab-initio molecular dynamics study of eutectic chloride salt: MgCl2–NaCl–KCl

Ionic liquid materials are viable candidates as a heat transfer fluid (HTF) in a wide range of applications, notably within concentrated solar power (CSP) technology and molten salt reactors (MSRs). For next-generation CSP and MSR technologies that strive for higher power generation efficiency, a HTF with wide liquid phase range and energy storage capabilities is crucial. Studies have shown that eutectic chloride salts exhibit thermal stability at high temperatures, high heat storage capacity, and are less expensive than nitrate and carbonate salts. However, the experimental data needed to fully evaluate the potential of eutectic chloride salts as a HTF contender are scarce and entail large uncertainties. Considering the high cost and potential hazards associated with the experimental methods used to determine the properties of ionic liquids, molecular modeling can be used as a viable alternative resource. In this study, the eutectic ternary chloride salt MgCl 2 –NaCl–KCl is modeled using ab-initio molecular dynamics simulations (AIMDs) in the liquid phase. Using the simulated data, the thermophysical and transport properties of eutectic chloride salt can be calculated: density, viscosity, heat capacity, diffusion coefficient, and ionic conductivity. For an initial model validation, experimental pair-distribution function data were obtained from X-ray total scattering techniques and compared to the theoretical pair-distribution function. Additionally, theoretical viscosity values are compared to experimental viscosity values for a similar system. The results provide a starting foundation for a MgCl 2 –NaCl–KCl model that can be extended to predict other fundamental properties.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Effective Thermal Conductivity of High Porosity Open Cell Nickel Foam

The effective thermal conductivity of high-porosity open cell nickel foam samples was measured over a wide range of temperatures and pressures using a standard steady-state technique. The samples, measuring 23.8 mm, 18.7 mm, and 13.6 mm in thickness, were constructed with layers of 1.7 mm thick foam with a porosity of 0.968. Tests were conducted with the specimens subjected to temperature differences of 100 to 1000 K across the thickness and at environmental pressures of 10(exp -4) to 750 mm Hg. All test were conducted in a gaseous nitrogen environment. A one-dimensional finite volume numerical model was developed to model combined radiation/conduction heat transfer in the foam. The radiation heat transfer was modeled using the two-flux approximation. Solid and gas conduction were modeled using standard techniques for high porosity media. A parameter estimation technique was used in conjunction with the measured and predicted thermal conductivities at pressures of 10(exp -4) and 750 mm Hg to determine the extinction coefficient, albedo of scattering, and weighting factors for modeling the conduction thermal conductivity. The measured and predicted conductivities over the intermediate pressure values differed by 13%.

Sullins, Alan D.↗

Heat transfer model for predicting squib ignition times

A squib ignition model based on transient heat condition from the hot bridgewire to the pyrotechnic is described. No Arrhenius-type chemical reaction is included. Instead, a thermal contact resistance is postulated to exist between the hot bridgewire and the pyrotechnic. Ignition is assumed to occur when a 2.5 micron layer of pyrotechnic next to the bridgewire reaches a characteristic ignition temperature for that pyrotechnic. This model was applied to the JPL squib, which uses a 50 micron (0.002-in.) diameter Tophet A bridgewire to ignite a boron, potassium perchlorate mix. A computer program was utilized that solves the transient heat condition problem with the boundary conditions stipulated by the model. The thermal contact conductance at the interface was determined by trial and error so that the experimentally determined ignition time for one firing condition would be properly predicted by the model. The agreement was quite good for tests run between -129 C and +93.3 C at current levels of 3.5 and 5 A. Axial heat conduction along the bridgewire is shown to be negligible.

Sernas, V.↗

Inhibition of conductive heat flow by magnetic constriction in the corona and transition region - Dependence on the shape of the constriction

Attention is given to the process by which the shape of the constriction of field lines inhibiting the conduction of heat down from the corona acts in conjunction with the degree of constriction to inhibit the heat flow. The analytical model for the heat flow involves tapered flux tubes in which the plasma properties are constant on cross sections, while the plasma is static and the only energy transfer is by thermal conduction. It is determined that only two specific solutions to the model are applicable to the solar atmosphere: the steady state, which is appropriate for quiet regions and for active ones that are not flaring, and the time-dependent case, in which no heat enters the hot end, which is appropriate for conductive cooling of flare loops.

Moore, R. L.↗

Model of hot-film sensor with substrate

A detailed model is constructed to investigate parameters influencing a hot-film anemometer system used to study cross-flow phenomena on swept wings. This anemometer is designed to detect vortices set up by cross-flow. The nature of the vortices is not well known, and a well-verified instrumentation model is considered essential to allow accurate interpretation of experimental results. Previous investigators have indicated that heat conduction to the film's substrate is significant. Two variations on the inclusion of substrate effects in the model are discussed. In the first, the substrate is assumed to have a constant temperature throughput: in the second more elaborate model, the temperature of the substrate is allowed to vary. This variation is related to the convective heat removal at the substrate/gas interface.

Mcree, Griffith J.↗

Finite Element A Posteriori Error Estimation for Heat Conduction

This research investigates residual-based a posteriori error estimates for finite element approximations of heat conduction in single-layer and multi-layered materials. The finite element approximation, based upon hierarchical modelling combined with p-version finite elements, is described with specific application to a two-dimensional, steady state, heat-conduction problem. Element error indicators are determined by solving an element equation for the error with the element residual as a source, and a global error estimate in the energy norm is computed by collecting the element contributions. Numerical results of the performance of the error estimate are presented by comparisons to the actual error. Two methods are discussed and compared for approximating the element boundary flux. The equilibrated flux method provides more accurate results for estimating the error than the average flux method. The error estimation is applied to multi-layered materials with a modification to the equilibrated flux method to approximate the discontinuous flux along a boundary at the material interfaces. A directional error indicator is developed which distinguishes between the hierarchical modeling error and the finite element error. Numerical results are presented for single-layered materials which show that the directional indicators accurately determine which contribution to the total error dominates.

Lang, Christapher G.↗

Assessment of a One-Dimensional Finite Element Charring Ablation Material Response Model for Phenolic-Impregnated Carbon Ablator

In this study, mathematical formulations to model charring ablation problems were numerically implemented using finite element analysis (FEA) with ABAQUS, which account for the material decomposition and progressive surface removal in the heat conduction and the surface energy balance equations. FEA was performed for a one-dimensional model to predict the temperature and ablation histories of a phenolic-impregnated carbon ablator sample (i.e., a common heat shield material for hypersonic vehicles and spacecraft) subjected to oxy-acetylene torch flame (i.e., 0.8 SLPM acetylene gas to 2.7 SLPM oxygen gas). The recovery enthalpy and convective heat transfer coefficient for the ablation model were calculated based on gas compositions and two assumed surface conditions (i.e., equilibrium and frozen). Simulations using the calculated recovery enthalpy and convective heat transfer coefficient resulted in a recession rate of 6.38 times (equilibrium) and 14.08 times (frozen) higher than the experimental data, despite fair agreement of the surface temperature. In addition, the effect of the heat transfer coefficient was investigated through a steady-state ablation analysis. The results of the analysis indicate that there is not one single value for the heat transfer coefficient that would allow the prediction to match both measured recession rate and surface temperature. Possible reasons for such an inconsistency are provided and discussed.

Thermal Protection System↗

Phase change paint tests to investigate effects of TPS tiles on heating rates of the Rockwell space shuttle orbiter (test OH4C, model 21-0)

Information and data from wind tunnel tests conducted on 0.0175-scale models of the space shuttle orbiter are presented. The primary objective of the tests was to evaluate aerodynamic heating effects of the tiles in the thermal protection system (TPS). Tile gap depth and flow orientation effects on the TPS were investigated. Tile patterns were cut into the undersides of the orbiter models to simulate the gaps. One model was left smooth for comparison.

Quan, M.↗

Thermal stress tectonics on the satellites of Saturn and Uranus

Thermal stress histories of the Saturnian and Uranian satellites are investigated. To this end, the thermal evolution of an icy satellite subjected to accretional and radiogenic heating, thermal conduction, and solid-state convection is modeled, and changes in the internal stress that occur during satellite evolution are examined. Results show that internal temperature changes that occur during normal evolution of many of the satellites of Saturn and Uranus can be expected to generate large extensional stresses in the satellites' outer regions. These stresses arise from three sources: (1) radiogenic warming, causing thermal expansion of materials in the satellite's deep interior; (2) radiogenic warming in larger satellites that can induce a phase transition from ice II to ice I and to produce a volume increase in the deep interior; and (3) accretional heating depositing heat in the satellite'e outer regions.

Hillier, John↗

The nature of the residual Martian polar caps

A model of the behavior of the Martian polar caps is described which incorporates the heating effects of the atmosphere, as well as insolation and conduction. This model is used to try to match the observed regression curves of the polar caps, and it predicts that all the seasonally condensed CO2 will be lost by around the summer solstice. The implication is that the residual caps are composed of water ice which, it is found by further modeling, should be stable during the Martian summers. However, it is also argued that this model may be too simplistic, and that the effects of wind in redistributing the seasonal condensate may lead to sufficient thicknesses of CO2 in the central polar region to allow the year-long existence of CO2 without significantly changing the retreat characteristics of the cap, and it is, therefore, concluded that at the present, the nature of the residual caps cannot be reliably determined.

Briggs, G. A.↗

Prediction of Ablation Rates from Solid Surfaces Exposed to High Temperature Gas Flow

A mathematical model and a solution algorithm is developed to study the physics of high temperature heat transfer and material ablation and identify the problems associated with the flow of hydrogen gas at very high temperatures and velocities through pipes and various components of Nuclear Thermal Rocket (NTR) motors. Ablation and melting can be experienced when the inner solid surface of the cooling channels and the diverging-converging nozzle of a Nuclear Thermal Rocket (NTR) motor is exposed to hydrogen gas flow at temperatures around 2500 degrees Kelvin and pressures around 3.4 MPa. In the experiments conducted on typical NTR motors developed in 1960s, degradation of the cooling channel material (cracking in the nuclear fuel element cladding) and in some instances melting of the core was observed. This paper presents the results of a preliminary study based on two types of physics based mathematical models that were developed to simulate the thermal-hydrodynamic conditions that lead to ablation of the solid surface of a stainless steel pipe exposed to high temperature hydrogen gas near sonic velocities. One of the proposed models is one-dimensional and assumes the gas flow to be unsteady, compressible and viscous. An in-house computer code was developed to solve the conservations equations of this model using a second-order accurate finite-difference technique. The second model assumes the flow to be three-dimensional, unsteady, compressible and viscous. A commercial CFD code (Fluent) was used to solve the later model equations. Both models assume the thermodynamic and transport properties of the hydrogen gas to be temperature dependent. In the solution algorithm developed for this study, the unsteady temperature of the pipe is determined from the heat equation for the solid. The solid-gas interface temperature is determined from an energy balance at the interface which includes heat transfer from or to the interface by conduction, convection, radiation, and ablation. Two different ablation models are proposed to determine the heat loss from the solid surface due to the ablation of the solid material. Both of them are physics based. Various numerical simulations were carried out using both models to predict the temperature distribution in the solid and in the gas flow, and then predict the ablation rates at a typical NTR motor hydrogen gas temperature and pressure. Solid mass loss rate per foot of a pipe was also calculated from these predictions. The results are presented for fully developed turbulent flow conditions in a sample SS pipe with a 6 inch diameter.

Akyuzlu, Kazim M.↗

Aerodynamic heating in hypersonic flows

Aerodynamic heating in hypersonic space vehicles is an important factor to be considered in their design. Therefore the designers of such vehicles need reliable heat transfer data in this respect for a successful design. Such data is usually produced by testing the models of hypersonic surfaces in wind tunnels. Most of the hypersonic test facilities at present are conventional blow-down tunnels whose run times are of the order of several seconds. The surface temperatures on such models are obtained using standard techniques such as thin-film resistance gages, thin-skin transient calorimeter gages and coaxial thermocouple or video acquisition systems such as phosphor thermography and infrared thermography. The data are usually reduced assuming that the model behaves like a semi-infinite solid (SIS) with constant properties and that heat transfer is by one-dimensional conduction only. This simplifying assumption may be valid in cases where models are thick, run-times short, and thermal diffusivities small. In many instances, however, when these conditions are not met, the assumption may lead to significant errors in the heat transfer results. The purpose of the present paper is to investigate this aspect. Specifically, the objectives are as follows: (1) to determine the limiting conditions under which a model can be considered a semi-infinite body; (2) to estimate the extent of errors involved in the reduction of the data if the models violate the assumption; and (3) to come up with correlation factors which when multiplied by the results obtained under the SIS assumption will provide the results under the actual conditions.

Reddy, C. Subba↗

Temporal behaviour of the thermal model of hard X-ray bursts

A simple, analytic model is presented of a hot, thermal hard X-ray source, continuously heated, bounded by ion-acoustic conduction fronts, and expanding in a loop. The model is used to investigate the assumption that the 'rise time' of the X-ray emission is approximately given by the loop length divided by the ion-sound speed appropriate to the peak temperature. It is found that a freely-expanding source does not behave in this way; instead, the rise time is symptomatic of the timescale for primary energy release. If the energy release rate does not fall significantly before the source fills the loop, however, then this assumption may be approximately satisfied, if a condition on the temporal behavior of the energy release is satisfied. Finally, some remarks on the relative timing of temperature and emission measure peaks are made, and possible further applications mentioned of the results presented herein.

Mackinnon, A. L.↗

Decoupled Method for Reconstruction of Surface Conditions From Internal Temperatures On Ablative Materials With Uncertain Recession Model

Obtaining measurements of flight environments on ablative heat shields is both critical for spacecraft development and extremely challenging due to the harsh heating environment and surface recession. Thermocouples installed several millimeters below the surface are commonly used to measure the heat shield temperature response, but an ill-posed inverse heat conduction problem must be solved to reconstruct the surface heating environment from these measurements. Ablation can contribute substantially to the measurement response making solutions to the inverse problem strongly dependent on the recession model, which is often poorly characterized. To enable efficient surface reconstruction for recession model sensitivity analysis, a method for decoupling the surface recession evaluation from the inverse heat conduction problem is presented. The decoupled method is shown to provide reconstructions of equivalent accuracy to the traditional coupled method but with substantially reduced computational effort. These methods are applied to reconstruct the environments on the Mars Science Laboratory heat shield using diffusion limit and kinetically limited recession models.

Oliver, A. Brandon↗