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 1,315 records · Page 73

The intercrater plains of Mercury and the Moon: Their nature, origin and role in terrestrial planet evolution. Thermal models of Mercury

Recent and more complex thermal models of Mercury and the terrestrial planets are discussed or noted. These models isolate a particular aspect of the planet's thermal history in an attempt to understand that parameter. Among these topics are thermal conductivity, convection, radiogenic sources of heat, other heat sources, and the problem of the molten core and regenerative dynamo.

Leake, M. A.↗

Thermomechanical modeling of the Colorado Plateau-Basin and range transition zone

The Colorado Plateau (CP) basin and range (B & R) boundary is marked by a transition zone on the order of 75 to 150 km in width. As one moves westward across this transition from the CP interior to the B & R there is a variation in the surface topography, surface heat flow, Bouguer gravity, seismicity, and crustal structure. This transition extends eastward into the western CP from the Wastach-Hurricane fault line and is largely coincident with the high plateaus of Utah and the Wasatch Mountains. It has been suggested that this transition zone marks a thermal and tectonic encroachment of the CP by the B & R. A simple two dimensional numerical model of the thermal regime for the transition zone was constructed to test the hypothesis that the observed geophysical signatures across the transition are due to lateral heat conduction from steady state uniform extension within the B & R lithosphere. Surface heat flow, uplift due to flexure from thermal buoyant loading, and regional Bouguer gravity are computed for various extension rates, crustal structures, and compensation depths.

Londe, M. D.↗

Mechanical and thermal properties of planetologically important ices

Two squences of ice composition were proposed for the icy satellites: a dense nebula model and a solar nebula model. Careful modeling of the structure, composition, and thermal history of satellites composed of these various ices requires quantitative information on the density, compressibility, thermal expansion, heat capacity, and thermal conductivity. Equations of state were fitted to the density data of the molecular ices. The unusual thermal and mechanical properties of the molecular and binary ices suggest a larger range of phenomena than previously anticipated, sufficiently complex perhaps to account for many of the unusual geologic phenomena found on the icy satellites.

Croft, Steven K.↗

Polar volatiles on Mercury and the Moon

The possibility that condensed volatiles might be stable in the permanently shadowed polar regions of the Moon and Mercury has been anticipated in a number of theoretical studies. In 1992, VLA-Goldstone and Arecibo observations revealed the presence of bright, depolarized radar features near the poles of Mercury that have been widely interpreted as evidence for the presence of polar ice deposits. Recently acquired high-resolution Arecibo radar images show that the anomalous radar features are concentrated in crater-sized regions whose locations can be made to coincide exactly with the locations of known impact craters in the Mariner-10 images. These new Arecibo images provide an unusual opportunity to learn more about the distribution and history of Mercury's polar ice deposits. We have constructed a thermal model that can predict surface and subsurface temperatures within impact craters on Mercury and the Moon. Included in the calculations are the effects of 1-dimensional subsurface heat conduction, direct sunlight, multiply reflected sunlight within the crater, and re-radiated infrared radiation within the crater. We also use realistic crater topographic profiles for larger flat-floored craters as well as smaller spherical bowl-shaped craters.

Paige, D. A.↗

CFD Simulations of the IHF 13-Inch Nozzle Flow: 55° Sphere-Cone Model, Manufactured Fences and Gaps

This paper reports computational analyses of tests in a high enthalpy arc-jet facility at NASA Ames Research Center. These tests were conducted using 23.6-cm diameter, 55° sphere-cone models placed in a free jet downstream of the 33.0-cm diameter conical nozzle in the 60-MW Interaction Heating Facility. Some of the sphere-cone models include surface features such as manufactured fences and gaps intended to simulate effects of differential recession, all of which disturb the flow, producing augmented heating locally and downstream. Test calibration data were obtained using slug and Gardon gage stagnation calorimeters, and a sphere-cone calorimeter model with six Gardon heat flux gages and five pressure gages. The present analysis comprises computational fluid dynamics simulations of the nonequilibrium flowfield in the facility nozzle and test box, including the models tested, and comparisons with the experimental measurements. Various issues related to testing are considered: diffuser flow capture, flow characterization based on the calorimeter data, effects of the model surface recession, and prediction of surface quantities for the models with surface features. These simulations take into account nonuniform total enthalpy profiles at the nozzle inlet as well as the expansion waves emanating from the nozzle exit and their effects on the model flowfields.

Arc-jets↗

Sockeye Validation Support Using the SPHERE Facility

The engineering-scale heat pipe application Sockeye, being developed through the Nuclear En- ergy Advanced Modeling and Simulation program is concurrently undergoing code development and verification and validation activities. The Single Primary Heat Extraction and Removal Emu- lator (SPHERE) facility at Idaho National Laboratory has performed numerous sodium heat pipe experiments, one of which is modeled in this work to support Sockeye validation. Heat pipe data produced by the SPHERE facility was found to be very useful not only for Sockeye validation but also as a useful resource for model development. A Sockeye simulation performed using Sockeye’s conduction model modeled an experimental run at SPHERE and was found to be able to accurately reproduce the steady temperature distribution but revealed discrepancies in transient behavior after shutoff of the power. This effort also revealed a number of challenges of modeling an experimental setup such as SPHERE, for example, the importance of details in the setup such as possible gaps in insulation. Also, continuous collaboration between Sockeye developers and SPHERE experimentalists has been established, which allows future experiments to be refined in response to questions and challenges found by simulating earlier experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CFD modeling of microwave electrothermal thrusters

Microwave-heated plasmas in convergent nozzles are analyzed using a coupled Maxwell and Navier-Stokes solver to examine relevant issues associated with microwave thermal propulsion. Parametric studies are conducted to understand the effect of power, pressure, and plasma location with respect to the nozzle throat. For nozzles in the 0.5 to 3 N range with helium flow, results show that specific impulses up to 550-650 seconds are possible, with further increases being limited by severe wall-heating. Coupling efficiencies of over 90 percent are consistently obtained, with overall efficiencies ranging from 40 percent to 80 percent. Size scale-up studies-done by scaling the frequency from 2.45 GHz to 0.91 GHz-indicate that plasma migration toward the walls occurs more frequently for the lower frequency. Increasing the cavity aspect ratio and detuning the cavity are found to be effective ways of keeping the plasma on axis.

Schwer, Douglas A.↗

An experimental and chemical kinetic modeling study of 4-butoxyheptane combustion

Here, the combustion kinetics of a novel oxygenated bioblendstock for diesel, 4-butoxyheptane (4-BH), was investigated experimentally using a flow reactor and a heated, high-pressure shock tube. The flow reactor experiments employed oxygen as the oxidizer and helium as the diluent with oxidation conducted at atmospheric pressure and 10 bar for temperatures from 400 to 1000 K at 20-K intervals. The fuel, oxidizer, and diluent flow rates were varied at different temperatures to maintain a constant initial fuel mole fraction of 1000 ppm, with stoichiometric equivalence ratio, and a residence time of 2.0 s. The reacted gas was fed to two separate GC systems that could qualitatively and quantitatively detect product species. Additionally, real fuel-air ignition delay time (IDT) data were collected using a heated, high-pressure shock-tube facility. Fuel lean (φ = 0.5) and stoichiometric (φ = 1.0) mixtures were investigated at 10 atm as well as at 30 atm for the fuel lean case for temperatures between 847 and 1259 K. A detailed chemical kinetics mechanism was developed to model the product distribution from the flow reactor and IDTs from the shock tube. The proposed model was able to predict the double NTC behavior in flow reactor experiments reasonably well. Model predictions at low temperatures were observed to be highly sensitive to the rate constants of ketohydroperoxide (KHP) decomposition in the case of the OOH group in α position which were modeled based on existing literature studies on ethers. It was noted that in the absence of theoretical or experimental studies, the rate constants for KHP decomposition used in the literature were empirically set. Additional studies are required to address the gap in model prediction obtained in this study and to reduce the uncertainty in kinetics models for ether oxidation. Predicted product concentrations and IDTs showed some quantitative agreement with experimental data, but the overall reactivity of the IDTs is underpredicted. Additionally, significant deviation is observed for the IDT results at 10 atm for the stoichiometric case with minor deviations for the other cases. The reaction pathways to the missing products were then further analyzed theoretically through quantum-mechanical calculations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Simplified Model of VIPER Thermal Management System. Part II: Integrated Vehicle

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) thermal management system (TMS) relies on four loop heat pipes (LHPs) to transport electronic waste heat to the vehicle cooling radiative surface and avoid overheating. The TMS has also ten constant conductance heat pipes (CCHPs) dedicated to balance the thermal load within the internal environment where the avionics boxes are mounted, also called warm electronic box (WEB), and to transport the heat from two of the science payload instruments. The TMS also uses two thermal straps to thermally link the batteries to the WEB. These thermal components, in addition to heaters, thermostat, multi-layer insulation (MLIs), and isolators forms the core of the VIPER TMS. The complex heat transport balance managed by the TMS is challenging to characterize and model. The more fidelity and granularity of a model, the more costly the computational resources needed and the longer the simulation and modeling time. When the priority is to provide quick but reliable assessments of the thermal performance or real time thermal feedback for training of console operators, simplified modeling tools are needed. To satisfy that need, this paper describes the effort to develop and correlate a model of VIPER TMS based on control volume approach. The correlation effort in particular focuses on hibernation, cold thermal balance, and hot thermal balance data from the integrated vehicle thermal vacuum (TVAC) test. Thus, the correlated model captures the heat leaks during hibernations and the performance at two extremes, bounding, operating scenarios.

Loop Heat Pipe↗

A Simplified Model of VIPER Thermal Management System. Part II: Integrated Vehicle

NASA’s Volatiles Investigating Polar Exploration Rover (VIPER) thermal management system (TMS) relies on four loop heat pipes (LHPs) to transport electronic waste heat to the vehicle cooling radiative surface and avoid overheating. The TMS has also ten constant conductance heat pipes (CCHPs) dedicated to balance the thermal load within the internal environment where the avionics boxes are mounted, also called warm electronic box (WEB), and to transport the heat from two of the science payload instruments. The TMS also uses two thermal straps to thermally link the batteries to the WEB. These thermal components, in addition to heaters, thermostat, multi-layer insulation (MLIs), and isolators forms the core of the VIPER TMS. The complex heat transport balance managed by the TMS is challenging to characterize and model. The more fidelity and granularity of a model, the more costly the computational resources needed and the longer the simulation and modeling time. When the priority is to provide quick but reliable assessments of the thermal performance or real time thermal feedback for training of console operators, simplified modeling tools are needed. To satisfy that need, this paper describes the effort to develop and correlate a model of VIPER TMS based on control volume approach. The correlation effort in particular focuses on hibernation, cold thermal balance, and hot thermal balance data from the integrated vehicle thermal vacuum (TVAC) test. Thus, the correlated model captures the heat leaks during hibernations and the performance at two extremes, bounding, operating scenarios.

Thermal Modeling↗

Icing Physics Studies Using the 3D SIDRM Test Article: Ice Crystal Icing Analysis

In-flight icing is an important safety issue and is a factor that affects aircraft design and performance. Newer regulations are driving a need for improvements in airframe and engine icing simulation capability. Experimental data is required for the development of icing physics models and simulation validation. This paper presents the analysis of the ice crystal icing data subset from tests conducted in 2022 at the NASA Icing Research Tunnel that studied both supercooled liquid and ice-crystal icing. The test article that was utilized replicated 3D geometrical features of an inter-compressor duct and strut region of a turbofan engine. The surfaces of the Simulated Inter-compressor Duct Research Model (SIDRM) can be heated to simulate the warm surfaces of the turbofan inter-compressor duct. The test article is instrumented with pressure taps, heaters, heat flux gauges, and thermocouples, while a 3D laser scanner, cameras, and a weighing scale to measure ice mass were utilized to characterize the icing behavior. The aim of these tests was to generate ice accretions on the SIDRM test article under well-characterized icing conditions. To that end, ice crystal cloud characterization tests were conducted that provided cloud property values at the test section. The ice crystal accretion tests investigated how different test condition parameters impacted ice accretion size, location, and characteristics (such as physical attributes). A key finding was that a sufficiently high surface heat flux at the leading edge can provide a continuous supply of liquid melt to critical accretion zones downstream, resulting in elevated ice mass accretions. Another key finding was that warmer heater settings in critical accretion zones can suppress ice accretion. These findings are consistent with icing behavior observed in full-scale engine icing tests conducted previously at the NASA Glenn Research Center. Test runs investigating the contribution of liquid melt generated at the leading edge found that the melt can be transported downstream to critical accretion regions in two forms, via runback and splash. The form of melt supplied downstream affected accretion behavior and features. Sharkteeth-shaped ice accretions were generated at certain test conditions, a physical attribute that has been observed in full-scale engine icing tests. The ice accretion data will be used to develop and validate 3D computational engine icing tools, such as GlennICE, that predictively assesses the onset and growth of ice. One of the goals of the sponsoring NASA project is to develop simulation models and tools that can assist in the design and certification of engines for flight in icing conditions in a cost‑effective way.

ice crystal icing↗

Icing Physics Studies Using the 3D SIDRM Test Article: Ice Crystal Icing Analysis

In-flight icing is an important safety issue and is a factor that affects aircraft design and performance. Newer regulations are driving a need for improvements in airframe and engine icing simulation capability. Experimental data is required for the development of icing physics models and simulation validation. This paper presents the analysis of the ice crystal icing data subset from tests conducted in 2022 at the NASA Icing Research Tunnel that studied both supercooled liquid and ice-crystal icing. The test article that was utilized replicated 3D geometrical features of an inter-compressor duct and strut region of a turbofan engine. The surfaces of the Simulated Inter-compressor Duct Research Model (SIDRM) can be heated to simulate the warm surfaces of the turbofan inter-compressor duct. The test article is instrumented with pressure taps, heaters, heat flux gauges, and thermocouples, while a 3D laser scanner, cameras, and a weighing scale to measure ice mass were utilized to characterize the icing behavior. The aim of these tests was to generate ice accretions on the SIDRM test article under well-characterized icing conditions. To that end, ice crystal cloud characterization tests were conducted that provided cloud property values at the test section. The ice crystal accretion tests investigated how different test condition parameters impacted ice accretion size, location, and characteristics (such as physical attributes). A key finding was that a sufficiently high surface heat flux at the leading edge can provide a continuous supply of liquid melt to critical accretion zones downstream, resulting in elevated ice mass accretions. Another key finding was that warmer heater settings in critical accretion zones can suppress ice accretion. These findings are consistent with icing behavior observed in full-scale engine icing tests conducted previously at the NASA Glenn Research Center. Test runs investigating the contribution of liquid melt generated at the leading edge found that the melt can be transported downstream to critical accretion regions in two forms, via runback and splash. The form of melt supplied downstream affected accretion behavior and features. Sharkteeth-shaped ice accretions were generated at certain test conditions, a physical attribute that has been observed in full-scale engine icing tests. The ice accretion data will be used to develop and validate 3D computational engine icing tools, such as GlennICE, that predictively assesses the onset and growth of ice. One of the goals of the sponsoring NASA project is to develop simulation models and tools that can assist in the design and certification of engines for flight in icing conditions in a cost‑effective way.

ice crystal icing↗

Prediction of orbiter RSI tile gap heating ratios from NASA/Ames double wedge model test

In-depth gap heating ratios for Orbiter RSI tile sidewalls were predicted based on near steady state temperature measurements obtained from double wedge model tests. An analysis was performed to derive gap heating ratios which would result in the best fit of test data; provide an assessment of open gap response, and supply the definition of gap filler requirements on the Orbiter. A comparison was made of these heating ratios with previously derived ratios in order to verify the extrapolation of the wing glove data to Orbiter flight conditions. The analysis was performed with the Rockwell TPS Multidimensional Heat Conduction Program for a 3-D, 2.0-inch thick flat RSI tile with 255 nodal points. The data from 14 tests was used to correlate with the analysis. The results show that the best-fit heating ratios at the station farthest upstream on the model for most gap depths were less than the extrapolated values of the wing glove model heating ratios. For the station farthest downstream on the model, the baseline heating ratios adequately predicted or over-predicted the test data.

Source record↗

Fuel Cell Thermal Management Through Conductive Cooling Plates

An analysis was performed to evaluate the concept of utilizing conductive cooling plates to remove heat from a fuel cell stack, as opposed to a conventional internal cooling loop. The potential advantages of this type of cooling system are reduced stack complexity and weight and increased reliability through the reduction of the number of internal fluid seals. The conductive cooling plates would extract heat from the stack transferring it to an external coolant loop. The analysis was performed to determine the required thickness of these plates. The analysis was based on an energy balance between the thermal energy produced within the stack and the heat removal from the cooling plates. To accomplish the energy balance, the heat flow into and along the plates to the cooling fluid was modeled. Results were generated for various numbers of cells being cooled by a single cooling plate. The results provided cooling plate thickness, mass, and operating temperature of the plates. It was determined that utilizing high-conductivity pyrolitic graphite cooling plates can provide a specific cooling capacity (W/kg) equivalent to or potentially greater than a conventional internal cooling loop system.

Colozza, Anthony J.↗

Heat current fluctuations and anomalous transport in low-dimensional carbon lattices

Molecular dynamics simulations and methods of importance sampling are used to study the heat transport of low-dimensional carbon lattices. For both carbon nanotubes and graphene sheets, heat transport is found to be anomalous, violating Fourier’s law of conduction with a system size dependent thermal conductivity and concomitant nonlinear temperature profiles. For carbon nanotubes, the thermal conductivity is found to increase as the square root of the length of the nanotube, while for graphene sheets the thermal conductivity is found to increase as the logarithm of the length of the sheet over the system sizes considered. Furthermore, the particular length dependence and nonlinear temperature profiles place carbon lattices into a universality class with nonlinear lattice models, and suggest that heat transport through carbon nanostructures is better described by a Levy walk rather than simple diffusion.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Determination of Complex Microcalorimeter Parameters with Impedance Measurements

The proper understanding and modeling of a microcalorimeter s response requires the accurate knowledge of a handful of parameters, such as C, G, alpha, . . . . While a few of these, such 8s the normal state resistance and the total thermal conductance to the heat bath (G) are directly determined from the DC IV characteristics, some others, notoriously the heat capacity (C) and alpha, appear in degenerate combinations in most measurable quantities. The case of a complex microcalorimeter, i.e. one in which the absorber s heat capacity is connected by a finite thermal impedance to the sensor, and subsequently by another thermal impedance to the heat bath, results in an added ambiguity in the determination of the individual C's and G's. In general, the dependence of the microcalorimeter s complex impedance on these parameters varies with frequency. This variation allows us to determine the individual parameters by fitting the prediction of the microcalorimeter model to the impedance data. We describe in this paper our efforts at characterizing the Goddard X-ray microcalorimeters. Using the parameters determined with this method we them compare the pulse shape and noise spectra predicted by the microcalorimeter model to data taken with the same devices.

Saab, T.↗

Modulation of the Seasonal Cycle of the Antarctic Sea Ice Extent By Sea Ice Processes and Feedbacks With the Ocean and the Atmosphere

The seasonal cycle of the Antarctic sea ice extent is strongly asymmetric, with a relatively slow increase after the summer minimum followed by a more rapid decrease after the winter maximum. This cycle is intimately linked to the seasonal cycle of the insolation received at the top of the atmosphere, but sea ice processes as well as the exchanges with the atmosphere and ocean may also play a role. To quantify these contributions, a series of idealized sensitivity experiments have been performed with an eddy-permitting (∘) NEMO-LIM3 (Nucleus for European Modelling of the Ocean–Louvain-la-Neuve sea ice model version 3) Southern Ocean configuration, including a representation of ice shelf cavities, in which the model was either driven by an atmospheric reanalysis or coupled to the COSMO-CLM2 regional atmospheric model. In those experiments, sea ice thermodynamics and dynamics as well as the exchanges with the ocean and atmosphere are strongly perturbed. This perturbation is achieved by modifying snow and ice thermal conductivities, the vertical mixing in the ocean top layers, the effect of freshwater uptake and release upon sea ice growth and melt, ice dynamics, and surface albedo. We find that the evolution of sea ice extent during the ice advance season is largely independent of the direct effect of the perturbation and appears thus mainly controlled by initial state in summer and subsequent insolation changes. In contrast, the melting rate varies strongly between the experiments during the retreat, in particular if the surface albedo or sea ice transport are modified, demonstrating a strong contribution of those elements to the evolution of ice coverage through spring and summer. As with the advance phase, the retreat is also influenced by conditions at the beginning of the melt season in September. Atmospheric feedbacks enhance the model winter ice extent response to any of the perturbed processes, and the enhancement is strongest when the albedo is modified. The response of sea ice volume and extent to changes in entrainment of subsurface warm waters to the ocean surface is also greatly amplified by the coupling with the atmosphere.

regional atmosphere ocean sea ice model↗

An "Inefficient Fin" Non-Dimensional Parameter to Measure Gas Temperatures Efficiently

A gas containment vessel that is not in thermal equilibrium with the bulk gas can affect its temperature measurement. The physical nature of many gas dynamics experiments often makes the accurate measurement of temperature a challenge. The environment itself typically requires that the thermocouple be sheathed, both to protect the wires and hot junction of the instrument from their environment, and to provide a smooth, rigid surface for pressure sealing of the enclosure. However, that enclosure may also be much colder than the gas to be sensed, or vice-versa. Either way, the effect of such gradients is to potentially skew the temperature measurements themselves, since heat may then be conducted by the instrument. Thermocouple designers traditionally address this problem by insulating the sheath from the thermocouple leads and hot junction as much as possible. The thermocouple leads are typically packed in a ceramic powder inside the sheath, protecting them somewhat from temperature gradients along the sheath, but there is no effective mechanism to shield the sheath from the enclosure body itself. Standard practice dictates that thermocouples be used in installations that do not present large thermal gradients along the probe. If this conduction dominates heat transfer near the tip of the probe, then temperature measurements may be expected to be skewed. While the same problem may be experienced in the measurement of temperature at various points within a solid in a gradient, it tends to be aggravated in the measurements of gas temperature, since heat transfer dependent on convection is often less efficient than conduction along the thermocouple. The proposed solution is an inefficient fin thermocouple probe. Conventional wisdom suggests that in many experiments where gas flows through an enclosure (e.g., flow in pipe, manifold, nozzle, etc.), the thermocouple be introduced flush to the surface, so as not to interfere with the flow. In practice, however, many such experiments take place where the flow is already turbulent, so that a protruding thermocouple probe has a negligible effect on the flow characteristics. The key question then becomes just how far into the flow should a thermocouple protrude in order to properly sense the gas temperature at that point. Modeling the thermocouple as an "inefficient fin" directly addresses this question. The appropriate assumptions in this case are: one-dimensional conduction along the fin; steady-state, constant, and homogeneous thermal conductivity; negligible radiation; and a uniform, constant heat transfer coefficient over the probe surface. It is noted that the nature of the ceramic-filled probe makes the key assumption of homogeneous thermal conductivity that much more conservative.

Lemieux, Patrick↗