Temperature field in multilayer systems with variable thermophysical properties
Mathematical model for heat conduction problems in multilayer systems with variable thermophysical properties
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.
Mathematical model for heat conduction problems in multilayer systems with variable thermophysical properties
Several thermophysical properties of molten silicon measured by the high-temperature electrostatic levitator at JPL are presented. They are density, constant-pressure specific heat capacity, hemispherical total emissivity, and surface tension. Over the temperature range investigated (1350 less than T(sub m) less than 1825 K), the measured liquid density (in g/cc) can be expressed by a quadratic function, pi(P) = p(sub m) - 1.69 x 10(exp -4)(T - T(sub m)) - 1.75 x 10(exp -7)((T - T(sub m))(exp 2) with T(sub m) and p(sub m) being 1687 K and 2.56 g/cc, respectively. The hemispherical total emissivity of molten silicon at the melting temperature was determined to be 0.18, and the constant-pressure specific heat was evaluated as a function of temperature. The surface tension (in 10(exp -3) N/m) of molten silicon over a similar temperature range can be expressed by sigma(T) = 875 - 0.22( T - T(sub m)).
Five thermophysical properties of molten silicon measured by the High Temperature Electrostatic Levitator (HTESL) at JPL are presented. The properties measured are the density, the constant pressure specific heat capacity, the hemispherical total emissivity, the surface tension and the viscosity.
Thermophysical property measurements on ablative materials for Apollo spacecraft propulsion system configuration
Accurate and efficient evaluations of refrigerant thermophysical properties and their partial derivatives are essential for transient simulations of thermofluid systems, where several computations need to be executed at each integration time step. Since the utilization of an Equation of State for retrieving properties based on a pair of independent inputs typically involves numerical iterations in solution procedures, when the input variables differ from the refrigerant state variables employed in dynamic models, a variety of approaches including lookup table interpolation and curve fitting have been developed to explicitly approximate these properties based on the state variables, and consequently eliminate internal iterations. This paper presents an alternative method that exploits derivative-informed neural networks to model refrigerant properties explicitly from inputs of pressure and enthalpy, while ensuring consistent partial derivatives generated by differentiating the neural networks. Computational speed and accuracy of the proposed approach are demonstrated via transient simulations of a discretized heat exchanger model in Modelica, and comparisons against other property evaluation routines. Simulation results indicate that the proposed approach can realize a significant speedup with negligible discrepancies in predicted transients. The method is implemented in an open-source Modelica library.
Density, viscosity, surface tension, electrical resistivity, thermoconductivity, enthalpy, specific heat, and vapor pressure thermophysical properties of potassium
The goal of this project is to gather thermophysical property data – specifically, thermal conductivity, thermal diffusivity, and specific heat – of heat transfer media (HTMs) and containment materials (CMs) used in Topic Area 1 and 2A at high temperatures (>700°C). This goal was met by completing objectives such as qualifying high temperature measurement techniques, defining allowable uncertainty metrics, reporting properties in a public database, and providing ongoing support for prioritized materials. These are described in detail across three budget periods
Reference work provides in tabular and graphical form the thermophysical properties of basic alloys, their constituents and oxides. This is useful for personnel who deal with extreme temperature environments.
Thermophysical properties data for oxygen at pressures below 5000 psia have been extrapolated to higher pressures (5,000-10,000 psia) in the temperature range 100-600 R. The tables include density, entropy, enthalpy, internal energy, speed of sound, specific heat, thermal conductivity, viscosity, thermal diffusivity, Prandtl number, and dielectric constant.
Tabular summary data of the thermophysical properties of fluid helium are given for temperatures from 0.8 to 1500 K, with pressures to 2000 MPa between 75 and 300 K, or to 100 MPa outside of this temperature band. Properties include density, specific heats, enthalpy, entropy, internal energy, sound velocity, expansivity, compressibility, thermal conductivity, and viscosity. The data are calculated from a computer program which is available from the National Institute of Standards and Technology. The computer program is based on carefully fitted state equations for both normal and superfluid helium.
In this paper, we review the current theory of the thermophysical properties of simple liquid metals. The emphasis is on thermodynamic properties, but we also briefly discuss the nonequilibrium properties of liquid metals. We begin by defining a 'simple liquid metal' as one in which the valence electrons interact only weakly with the ionic cores, so that the interaction can be treated by perturbation theory. We then write down the equilibrium Hamiltonian of a liquid metal as a sum of five terms: the bare ion-ion interaction, the electron-electron interaction, the bare electron-ion interaction, and the kinetic energies of electrons and ions. Since the electron-ion interaction can be treated by perturbation, the electronic part contributes in two ways to the Helmholtz free energy: it gives a density-dependent term which is independent of the arrangement of ions, and it acts to screen the ion-ion interaction, giving rise to effective ion-ion pair potentials which are density-dependent, in general. After sketching the form of a typical pair potential, we briefly enumerate some methods for calculating the ionic distribution function and hence the Helmholtz free energy of the liquid: monte Carlo simulations, molecular dynamics simulations, and thermodynamic perturbation theory. The final result is a general expression for the Helmholtz free energy of the liquid metal. It can be used to calculate a wide range of thermodynamic properties of simple metal liquids, which we enumerate. They include not only a range of thermodynamic coefficients of both metals and alloys, but also many aspects of the phase diagram, including freezing curves of pure elements and phase diagrams of liquid alloys (including liquidus and solidus curves). We briefly mention some key discoveries resulting from previous applications of this method, and point out that the same methods work for other materials not normally considered to be liquid metals (such as colloidal suspensions, in which the suspended microspheres behave like ions screened by the salt solution in which they are suspended). We conclude with a brief discussion of some non-equilibrium (i.e., transport) properties which can be treated by an extension of these methods. These include electrical resistivity, thermal conductivity, viscosity, atomic self-diffusion coefficients, concentration diffusion coefficients in alloys, surface tension and thermal emissivity. Finally, we briefly mention two methods by which the theory might be extended to non-simple liquid metals: these are empirical techniques (i.e., empirical two- and three-body potentials), and numerical many-body approaches. Both may be potentially applicable to extremely complex systems, such as nonstoichiometric liquid semiconductor alloys.
The surface tension, viscosity, and density of several bulk metallic glass-forming alloys have been measured using noncontact techniques in the electrostatic levitation facility (ESL) at NASA Marshall Space Flight Center. All three properties show unexpected behavior in the undercooled regime. Similar deviations were previously observed in titanium-based quasicrystal-forming alloys,but the deviations in the properties of the glass-forming alloys are much more pronounced. New results for anomalous thermophysical properties in undercooled glass-forming alloys will be presented and discussed.
Heat transfer within the lunar surface layer depends on several thermophysical properties of the lunar regolith, including the thermal conductivity, the specific heat, the thermal diffusivity, and the thermal parameter. Results of property measurements on simulated lunar materials are presented where appropriate as well as measurements made on the actual samples themselves. The variation of temperature on the moon with depth is considered, taking into account various times of the lunar day. The daily variation in temperature drops to about 1 deg at a depth of only 0.172 meters. The steady temperature on the moon below this depth is 225 K.
In situ determination of effects of ultraviolet radiation, solar wind protons, and temperature on thermophysical properties of solar cell filters and other spacecraft materials
Model can be tested directly, which eliminates costly, time-consuming, and inaccurate method of making test models solely for purpose of determining thermophysical properties. Method is adaptable to numerous modifications and variations.
Here, we present worm-algorithm, path integral Monte Carlo simulations of bulk liquid tritium. The simulations are benchmarked against empirically known thermophysical properties of liquid deuterium and liquid tritium. Results for the pair correlation function, chemical potential, isothermal compressibility, isochoric heat capacity, and single-particle momentum distributions are reported. Given the benchmark comparisons, our predictions of liquid tritium properties are expected to be accurate to within a few percent. Our simulations unambiguously demonstrate the significance of nuclear quantum effects to the properties of liquid tritium. In particular, under saturated vapor pressure, the average molecular kinetic energy of the liquid is found to be more than 60% higher than the value expected from the classical equipartition theorem.
Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.
This two-part paper reports the thermophysical properties of several cold and vacuum plasma sprayed monolithic Cu and Ni-based alloy coatings. Part I presents the electrical and thermal conductivity, thermal diffusivity, and total hemispherical emissivity data while Part II reports the specific heat capacity data for these coatings. Metallic copper alloys, stoichiometric NiAl and NiCrAlY coatings were fabricated by either the cold sprayed or the vacuum plasma spray deposition processes for thermal property measurements between 77 and 1223 K. The temperature dependencies of the thermal conductivities, thermal diffusivities, electrical conductivities and total hemispherical emissivities of these cold and vacuum sprayed monolithic coatings are reported in this paper. The electrical and thermal conductivity data correlate reasonably well for Cu-8%Cr-1%Al, Cu-23%Cr-5%Al and NiAl in accordance with the Wiedemann-Franz (WF) law although a better fit is obtained using the Smith-Palmer relationship. The Lorentz numbers determined from the WF law are close to the theoretical value.