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 343 records · Page 19

gLaSDI: Parametric physics-informed greedy latent space dynamics identification

A parametric adaptive physics-informed greedy Latent Space Dynamics Identification (gLaSDI) method is proposed for accurate, efficient, and robust data-driven reduced-order modeling of high-dimensional nonlinear dynamical systems. In the proposed gLaSDI framework, an autoencoder discovers intrinsic nonlinear latent representations of high-dimensional data, while dynamics identification (DI) models capture local latent-space dynamics. Here, an interactive training algorithm is adopted for the autoencoder and local DI models, which enables identification of simple latent-space dynamics and enhances accuracy and efficiency of data-driven reduced-order modeling. To maximize and accelerate the exploration of the parameter space for the optimal model performance, an adaptive greedy sampling algorithm integrated with a physics-informed residual-based error indicator and random-subset evaluation is introduced to search for the optimal training samples on the fly. Further, to exploit local latent-space dynamics captured by the local DI models for an improved modeling accuracy with a minimum number of local DI models in the parameter space, a -nearest neighbor convex interpolation scheme is employed. The effectiveness of the proposed framework is demonstrated by modeling various nonlinear dynamical problems, including Burgers equations, nonlinear heat conduction, and radial advection. The proposed adaptive greedy sampling outperforms the conventional predefined uniform sampling in terms of accuracy. Compared with the high-fidelity models, gLaSDI achieves 17 to 2,658× speed-up with 1 to 5% relative errors.

97 MATHEMATICS AND COMPUTING↗

Transient thermal response of ablating bodies

A numerical study of transient thermal response of a blunt-nosed axisymmetric body made of Teflon is presented using a two-layer thermal model. It is shown that phase change and transverse heat conduction have a considerable effect on the internal temperature field. Comparison of the numerical results with experimental data shows that the single-layer thermal model does not predict the real feature of the thermal field, whereas the results of the two-layer thermal model agree reasonably well with the experiment.

Arai, N.↗

Quantifying Thermal Output of Energetic Materials (LDRD Final Report)

We present the results of an LDRD project, funded by the Nuclear Deterrence IA, to develop capabilities for quantitative assessment of pyrotechnic thermal output. The thermal battery igniter is used as an exemplar system. Experimental methodologies for thermal output evaluation are demonstrated here, which can help designers and engineers better specify pyrotechnic components , provide thermal output guidelines for new formulations, and generate new metrics for assessing component performance and margin given a known failure condition. A heat-transfer analysis confirms that the dominant mode of energy transfer from the pyrotechnic output plume to the heat pellet is conduction via deposition of hot titanium particles. A simple lumped-parameter model of titanium particle heat transfer and a detailed multi-phase model of deposition heat transfer are discussed. Pyrotechnic function, as defined by "go/no-go" standoff testing of a heat pellet, is correlated with experimentally measured igniter plume temperature, titanium metal particle temperature, and energy deposition. Three high-speed thermal diagnostics were developed for this task. A three-color imaging pyrometer, acquiring 100k images per second on three color channels, is deployed for measurement of titanium particle temperatures. Complimentary measurements of the overall igniter plume emission ("color") temperature were conducted using a transmission-grating spectrograph in line-imaging mode. Heat flux and energy deposition to a cold wall at the heat-pellet location were estimated using an eroding thermocouple probe, with a frequency response of ~5 kHz. Ultimate "go/no-go" function in the igniter/heat-pellet system was correlated with quantitative thermal metrics, in particular surface energy deposition and plume color temperature. Titanium metal-particle and plume color temperatures both experience an upper bound approximated by the 3245-K boiling point of TiO 2 . Average metal-particle temperatures remained nearly constant for all standoff distances at T = 2850 K, ± 300 K, while plume color temperature and heat flux decay with standoff—suggesting that heat-pellet failure results from a drop in metal-particle flux and not particle temperature. At 50% likelihood of heat-pellet failure, peak time-resolved plume color temperatures drop well below TiO 2 boiling to ~2000 - 2200 K, near the TiO 2 melting point. Estimates of peak heat flux decline from up to 1 GW/m 2 for near-field standoffs to below 320 MW/m 2 at 50% failure likelihood.

36 MATERIALS SCIENCE↗

Alternating conduction and convection drying of paper – an experimental analysis with a continuous data acquisition approach

In conventional multi-cylinder drying of paper and board, both conductive drying from steam-heated dryer cylinders and convective drying by flowing air over the paper surface in the pockets are used. Conductive drying from steam-heated drying cylinders is a critical component in providing the necessary thermal energy to paper and board as they dry. Steam temperature and internal and external resistances at the contacting surface are critical process parameters influencing the conductive drying process. An experimental setup was developed to study the alternating conductive and convective drying of paper and board. Paper sheet moisture, temperature, and temperature distribution within the heated platen and the instantaneous heat flux as the sheet was being dried were measured. The instantaneous heat flux, contact heat transfer coefficient, and drying rates were determined as drying proceeds. Experimental results, as well as comparisons to literature and commercial data, are presented. The conductive heat transfer coefficients determined were compared to traditional correlations normally used in the modeling of paper drying. Similarly, the convective heat and mass transfer coefficients are also determined and compared to literature data. In addition to the evaluation of alternating conductive and convective drying characteristics of paper and board, the potential inclusion of auxiliary energy components will also be included. Experimental results from the conduction and convection drying system are presented. Furthermore, this data will be useful in process development, intensification of manufacturing processes, and modeling and simulation of paper drying processes.

42 ENGINEERING↗

Resin Film Infusion (RFI) Process Modeling for Large Transport Aircraft Wing Structures

This investigation completed the verification of a three-dimensional resin transfer molding/resin film infusion (RTM/RFI) process simulation model. The model incorporates resin flow through an anisotropic carbon fiber preform, cure kinetics of the resin, and heat transfer within the preform/tool assembly. The computer model can predict the flow front location, resin pressure distribution, and thermal profiles in the modeled part. The formulation for the flow model is given using the finite element/control volume (FE/CV) technique based on Darcy's Law of creeping flow through a porous media. The FE/CV technique is a numerically efficient method for finding the flow front location and the fluid pressure. The heat transfer model is based on the three-dimensional, transient heat conduction equation, including heat generation. Boundary conditions include specified temperature and convection. The code was designed with a modular approach so the flow and/or the thermal module may be turned on or off as desired. Both models are solved sequentially in a quasi-steady state fashion. A mesh refinement study was completed on a one-element thick model to determine the recommended size of elements that would result in a converged model for a typical RFI analysis. Guidelines are established for checking the convergence of a model, and the recommended element sizes are listed. Several experiments were conducted and computer simulations of the experiments were run to verify the simulation model. Isothermal, non-reacting flow in a T-stiffened section was simulated to verify the flow module. Predicted infiltration times were within 12-20% of measured times. The predicted pressures were approximately 50% of the measured pressures. A study was performed to attempt to explain the difference in pressures. Non-isothermal experiments with a reactive resin were modeled to verify the thermal module and the resin model. Two panels were manufactured using the RFI process. One was a stepped panel and the other was a panel with two 'T' stiffeners. The difference between the predicted infiltration times and the experimental times was 4% to 23%.

Loos, Alfred C.↗

Empirical estimation of densities in NaCl-KCl-UCl 3 and NaCl-KCl-YCl 3 molten salts using Redlich-Kister expansion

Densities of molten KCl-NaCl-UCl 3 and KCl-NaCl-YCl 3 ternary systems have been estimated using a multidimensional Redlich-Kister model. Temperature and composition dependent Redlich-Kister functions have been used to generate binary interaction parameters in the outlined ternary salt systems. These binary interactions have been used in the extrapolation to ternary system densities. The results of the density extrapolations by Muggianu interpolation scheme provide agreement within 2–3% for the NaCl-KCl-YCl 3 liquids and 11% in NaCl-KCl-UCl 3 liquids compared to the available experimental data. Modeling NaCl-KCl-UCl 3 molten phase density with a ternary interaction parameter improved the agreement within 4%. Thermophysical modeling used in this study has shown promising results for use in other material properties, such as viscosity, thermal conductivity, and heat capacity of the molten salts. Lastly, the outlined modeling method applied in these specific molten salt ternaries can be used for quaternary or higher multicomponent molten salt systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The effects of Reynolds number, rotor incidence angle, and surface roughness on the heat transfer distribution in a large-scale turbine rotor passage

A combined experimental and computational program was conducted to examine the heat transfer distribution in a turbine rotor passage geometrically similiar to the Space Shuttle Main Engine (SSME) High Pressure Fuel Turbopump (HPFTP). Heat transfer was measured and computed for both the full-span suction and pressure surfaces of the rotor airfoil as well as for the hub endwall surface. The primary objective of the program was to provide a benchmark-quality data base for the assessment of rotor passage heat transfer computational procedures. The experimental portion of the study was conducted in a large-scale, ambient temperature, rotating turbine model. Heat transfer data were obtained using thermocouple and liquid-crystal techniques to measure temperature distributions on the thin, electrically-heated skin of the rotor passage model. Test data were obtained for various combinations of Reynolds number, rotor incidence angle and model surface roughness. The data are reported in the form of contour maps of Stanton number. These heat distribution maps revealed numerous local effects produced by the three-dimensional flows within the rotor passage. Of particular importance were regions of local enhancement produced on the airfoil suction surface by the main-passage and tip-leakage vortices and on the hub endwall by the leading-edge horseshoe vortex system. The computational portion consisted of the application of a well-posed parabolized Navier-Stokes analysis to the calculation of the three-dimensional viscous flow through ducts simulating the a gas turbine passage. These cases include a 90 deg turning duct, a gas turbine cascade simulating a stator passage, and a gas turbine rotor passage including Coriolis forces. The calculated results were evaluated using experimental data of the three-dimensional velocity fields, wall static pressures, and wall heat transfer on the suction surface of the turbine airfoil and on the end wall. Particular attention was paid to an accurate modeling of the passage vortex and to the development of the wall boundary layers including crossflow.

Blair, Michael F.↗

Experimental Investigation of Roughness Effects on Transition Onset and Turbulent Heating Augmentation on a Hemisphere at Mach 6 and Mach 10

An experimental investigation of the effects of distributed surface roughness on boundary-layer transition and turbulent heating has been conducted. Hypersonic wind tunnel testing was performed using hemispherical models with surface roughness patterns simulating those produced by heat shield ablation. Global aeroheating and transition onset data were obtained using phosphor thermography at Mach 6 and Mach 10 over a range of roughness heights and free stream Reynolds numbers sufficient to produce laminar, transitional and turbulent flow. Upstream movement of the transition onset location and increasing heating augmentation over predicted smooth-wall levels were observed with both increasing roughness heights and increasing free stream Reynolds numbers. The experimental heating data are presented herein, as are comparisons to smooth-wall heat transfer distributions from computational flow-field simulations. The transition onset data are also tabulated, and correlations of these data are presented.

Hollis, Brian R.↗

Mechanical energy flow models of rods and beams

It has been proposed that the flow of mechanical energy through a structural/acoustic system may be modeled in a manner similar to that of flow of thermal energy/in a heat conduction problem. If this hypothesis is true, it would result in relatively efficient numerical models of structure-borne energy in large built-up structures. Fewer parameters are required to approximate the energy solution than are required to model the characteristic wave behavior of structural vibration by using traditional displacement formulations. The energy flow hypothesis is tested in this investigation for both longitudinal vibration in rods and transverse flexural vibrations of beams. The rod is shown to behave approximately according to the thermal energy flow analogy. However, the beam solutions behave significantly differently than predicted by the thermal analogy unless locally-space-averaged energy and power are considered. Several techniques for coupling dissimilar rods and beams are also discussed. Illustrations of the solution accuracy of the methods are included.

Wohlever, J. C.↗

Multispectral Emission of the Sun during the First Whole Sun Month: Magnetohydrodynamic Simulations

We demonstrate that a three-dimensional magnetohydrodynamic (MHD) simulation of the corona can model its global plasma density and temperature structure with sufficient accuracy to reproduce many of the multispectral properties of the corona observed in extreme ultraviolet (EW) and X-ray emission. The key ingredient to this new type of global MHD model is the inclusion of energy transport processes (coronal heating, anisotropic thermal conduction, and radiative losses) in the energy equation. The calculation of these processes has previously been confined to one-dimensional loop models, idealized two-dimensional computations, and three-dimensional active region models. We refer to this as the thermodynamic MHD model, and we apply it to the time period of Carrington rotation 1913 (1996 August 22 to September 18). The form of the coronal heating term strongly affects the plasma density and temperature of the solutions. We perform our calculation for three different empirical heating models: (1) a heating function exponentially decreasing in radius; (2) the model of Schrijver et al.; and (3) a model reproducing the heating properties of the quiet Sun and active regions. We produce synthetic emission images from the density and temperature calculated with these three heating functions and quantitatively compare them with observations from E W Imaging Telescope on the Solar and Heliospheric Observatory and the soft X-ray telescope on Yohkoh. Although none of the heating models provide a perfect match, heating models 2 and 3 provide a reasonable match to the observations.

Lionello, Roberto↗

Validation of PICA Ablation and Thermal-Response Model at Low Heat Flux

Phenolic Impregnated Carbon Ablator (PICA) was the forebody heatshield material on the Stardust sample-return capsule and is also a primary candidate material for the Mars Science Lander (MSL), the Orion Crew Module, and the SpaceX Dragon vehicle. As part of the heatshield qualification for Orion, physical and thermal properties of virgin and charred PICA were measured, and an ablation and thermal response model was developed. We validated the model by comparing it with recession and temperature data from stagnation arcjet tests conducted over a wide range of stagnation heat flux of 107 to 1102 W/sq cm. The effect of orthotropic thermal conductivity was evident in the thermal response of the arcjet models. In general, model predictions compared well with the data; however, the uncertainty of the recession prediction was greatest for heat fluxes below 200 W/sq cm. More recent MSL testing focused on the low heat flux regime of 45 to 250 W/sq cm. The new results confirm the recession uncertainty, especially for pressures below 6 kPa. In this work we focus on improving the model predictions for MSL and Orion tests below 250 W/sq cm.

Milos, Frank S.↗

Using Thermal Desktop to Model the Effects of Plume Heating on MLI

Gateway draws on experience from International Space Station (ISS) and geostationary satellites for many aspects of design, including Multilayer Insulation (MLI) material selection. MLI is used on many external surfaces, some of which can be exposed to high plume heating from visiting vehicle during docking and undocking. Due to Gateway’s smaller size and closer proximity of sensitive hardware to plumes from visiting vehicles, the expected and contingency plume heating loads are higher than those experienced by the ISS. Because of MLI’s low mass and thermal conductivity between layers, high plume heating loads can potentially heat MLI past its thermal limits. MLI is typically modeled as a single node of insulation with a given effective emissivity. This method is not adequate for modeling plume heating impacts because the temperature of each individual layer cannot be determined. Modeling MLI is complicated because there are several layers and the conduction between each layer is undefined. It is further complicated by the variety of material options for MLI and their respective mass, optics, transmissivity, and temperature limits. This paper shows modeling of a MLI blanket exposed to solar and plume heating loads using Thermal Desktop. Several sensitivity studies were run including examining distance between layers, varying effective emissivity, effects of different MLI materials, and plume heating loads and durations.

MLI↗

Using Thermal Desktop to Model the Effects of Plume Heating on MLI

Gateway draws on experience from International Space Station (ISS) and geostationary satellites for many aspects of design, including Multilayer Insulation (MLI) material selection. MLI is used on many external surfaces, some of which can be exposed to high plume heating from visiting vehicle during docking and undocking. Due to Gateway’s smaller size and closer proximity of sensitive hardware to plumes from visiting vehicles, the expected and contingency plume heating loads are higher than those experienced by the ISS. Because of MLI’s low mass and thermal conductivity between layers, high plume heating loads can potentially heat MLI past its thermal limits. MLI is typically modeled as a single node of insulation with a given effective emissivity. This method is not adequate for modeling plume heating impacts because the temperature of each individual layer cannot be determined. Modeling MLI is complicated because there are several layers and the conduction between each layer is undefined. It is further complicated by the variety of material options for MLI and their respective mass, optics, transmissivity, and temperature limits. This paper shows modeling of a MLI blanket exposed to solar and plume heating loads using Thermal Desktop. Several sensitivity studies were run including examining distance between layers, varying effective emissivity, effects of different MLI materials, and plume heating loads and durations.

MLI↗

An Optothermal Raman Model Accounting for Convection, Radiation, and Substrate Effects for Evaluation of the Thermal Conductivity of Two‐Dimensional Coatings

Thermal conductivity of two‐dimensional (2D) materials is important to characterize as this will govern thermal transport physics in the various thermal and energy applications benefitting from the unique properties of these materials. Here, analytical and numerical techniques are presented to determine the thermal conductivity and interfacial conductance of a 2D material coating using confocal micro‐Raman spectroscopy as a noncontact diagnostic. While several methods have been proposed for calculating the values of thermal conductivity, interfacial conductance, and spatial temperature profile, these models often do not consider effects of radiation, convection, and substrate resistance on the temperature and heat flux profiles. We present a model to calculate the thermal conductivity and interfacial conductance which accounts for convection, radiation, and substrate effects to characterize a variety of 2D material coatings, which is demonstrated using large area graphene transferred onto copper and nickel substrates. Convection and radiation effects are found to have a negligible effect on the temperature profile of supported coatings, while the substrate effects have a considerable effect and therefore impact the quantification of thermal conductivity and interfacial conductance in addition to limiting the spatial resolution of the technique. These findings will allow for more accurate and reproducible extraction of the thermal conductivity of graphene and other 2D material coatings when employing noncontact optothermal Raman spectroscopy methods.

2D coatings↗

The influence of thermal inertia on Mars' seasonal pressure variation and the effect of the weather component

Using a Leighton-Murray type diurnal and seasonal Mars thermal model, we found that it is possible to reproduce the seasonal variation in daily-averaged pressures (approximately 680-890 Pa) measured by Viking Lander 1 (VL1), during years without global dust storms, with a standard deviation of less than 5 Pa. In this simple model, surface CO2, frost condensation, and sublimation rates at each latitude are determined by the net effects of radiation, latent heat, and heat conduction in subsurface soil layers. An inherent assumption of our model is that the seasonal pressure variation is due entirely to the exchange of mass between the atmosphere and polar caps. However, the results of recent Mars GCM modeling have made it clear that there is a significant dynamical contribution to the seasonal pressure variation. This 'weather' component is primarily due to large-scale changes in atmospheric circulation, and its magnitude depends somewhat on the dust content of the atmosphere. The overall form of the theoretical weather component at the location of VL1, as calculated by the AMES GCM, remains the same over the typical range of Mars dust opacities.

Wood, S. E.↗

Verification of MOOSE/Bison's Heat Conduction Solver Using Combined Spatiotemporal Convergence Analysis

Bison is a computational physics code that uses the finite element method to model the thermo-mechanical response of nuclear fuel. Since Bison is used to inform high-consequence decisions, it is important that its computational results are reliable and predictive. One important step in assessing the reliability and predictive capabilities of a simulation tool is the verification process, which quantifies numerical errors in a discrete solution relative to the exact solution of the mathematical model. One step in the verification process—called code verification—ensures that the implemented numerical algorithm is a faithful representation of the underlying mathematical model, including partial differential or integral equations, initial and boundary conditions, and auxiliary relationships. In this paper, the code verification process is applied to spatiotemporal heat conduction problems in Bison. Simultaneous refinement of the discretization in space and time is employed to reveal any potential mistakes in the numerical algorithms for the interactions between the spatial and temporal components of the solution. For each verification problem, the correct spatial and temporal order of accuracy is demonstrated for both first- and second-order accurate finite elements and a variety of time-integration schemes. Furthermore, these results provide strong evidence that the Bison numerical algorithm for solving spatiotemporal problems reliably represents the underlying mathematical model in MOOSE. The selected test problems can also be used in other simulation tools that numerically solve for conduction or diffusion.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

JWST ISIM Harness Thermal Evaluation

The James Webb Space Telescope (JWST) will be a large infrared telescope with a 6.5-meter primary mirror. Launch is planned for 2013. JWST wl1 be the premier observatory of the next decade serving thousands of astronomers worldwide. The Integrated Science Instrument Module (ISIM) is the unit that will house thc four main JWST instruments. The ISIM enclosure passively cooled to 37 Kelvin and has a tightly managed thermal budget. A significant portion of the ISIM heat load is due to parasitic heat gains from the instrument harnesses. These harnesses provide a thermal path from the Instrument Electronics Control (IEC) to the ISIM. Because of the impact of this load to the ISIM thermal design, understanding the harness parasitic heat gains is critical. To this effect, a thermal test program has been conducted in order to characterize these parasitic loads and verify harness thermal models. Recent parasitic heat loads tests resulted in the addition of a dedicated multiple stage harness radiator. In order for the radiator to efficiently reject heat from the harness, effective thermal contact conductance values for multiple harnesses had to be determined. This presentation will describe the details and the results of this test program.

Kobel, Mark↗

On the preferential acceleration and heating of solar wind heavy ions

It is pointed out that the properties of solar wind ions heavier than H(+) have presented a puzzle to theorists for some time. The most recent theoretical efforts to explain the observations are considered, taking into account studies conducted by Dusenbery and Hollweg (1981) and Marsch et al. (1982). The first model of a wave-driven, three-fluid, supersonic solar wind is considered. Details of the resonant wave-particle interaction are discussed and the relevant equations for the model are derived. An investigation is conducted of the effects of resonant cyclotron acceleration and heating by a spectrum of nondispersive ion cyclotron waves on the model. It is found that there are serious difficulties with the hypothesis that the resonant cyclotron interaction is responsible for the preferential heating and acceleration of heavy ions in the solar wind. The model presented here is unable to produce the required energization by this mechanism. However, this model should be viewed as an approximate one, subject to possible correction by more advanced calculations.

Isenberg, P. A.↗