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,135 records · Page 63

Analytical solution and parameter estimation for heat of wetting and vapor adsorption during spontaneous imbibition in tuff

An analytical expression is derived for the thermal response observed during spontaneous imbibition of water into a dry core of zeolitic tuff. Sample tortuosity, thermal conductivity, and thermal source strength are estimated from fitting an analytical solution to temperature observations during a single laboratory test. Here, the closed-form analytical solution is derived using Green's functions for heat conduction in the limit of “slow” water movement; that is, when advection of thermal energy with the wetting front is negligible. The solution has four free fitting parameters and is efficient for parameter estimation. Laboratory imbibition data used to constrain the model include a time series of the mass of water imbibed, visual location of the wetting front through time, and temperature time series at six locations. The thermal front reached the end of the core hours before the visible wetting front. Thus, the predominant form of heating during imbibition in this zeolitic tuff is due to vapor adsorption in dry zeolitic rock ahead of the wetting front. The separation of the wetting front and thermal front in this zeolitic tuff is significant, compared to wetting front behavior of most materials reported in the literature. This work is the first interpretation of a thermal imbibition response to estimate transport (tortuosity) and thermal properties (including thermal conductivity) from a single laboratory test.

36 MATERIALS SCIENCE↗

Improved multiphysics model of the High Temperature Engineering Test Reactor for the simulation of loss-of-forced-cooling experiments

We present a multiphysics model of the High Temperature Engineering Test Reactor for comparison with past and predict future loss-of-forced-cooling (LOFC) experiments. The approach selected combines (1) 3-D full-core superhomogenization-corrected neutronics, (2) 3-D full-core homogenized or semi-heterogeneous heat transfer (macroscale), (3) 2-D axisymmetric fuel rod heat transfer (mesoscale), and (4) 1-D thermal-hydraulics channels. Although large sensitivities remain and the available experimental data are insufficient for use in a rigorous validation exercise, the overall transient behavior is successfully reproduced. The novelty of our work includes (1) a new macroscale-mesoscale heat transfer coupling approach relying on gap conductance to drastically speed up numerical convergence by two orders of magnitude, (2) determination of a radial effective thermal conductivity that yields reasonably close results between the homogenized and semi-heterogeneous macroscale thermal models, and (3) a preliminary study of the reactor's early behavior following a LOFC event, thus enabling further assessment of numerical models against fission power measurements.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Two-dimensional thermal modeling of power monolithic microwave integrated circuits (MMIC's)

Numerical simulations of the two-dimensional temperature distributions for a typical GaAs MMIC circuit are conducted, aiming at understanding the heat conduction process of the circuit chip and providing temperature information for device reliability analysis. The method used is to solve the two-dimensional heat conduction equation with a control-volume-based finite difference scheme. In particular, the effects of the power dissipation and the ambient temperature are examined, and the criterion for the worst operating environment is discussed in terms of the allowed highest device junction temperature.

Fan, Mark S.↗

Modeling laser interactions with aluminum and tantalum targets using a hybrid atomistic-continuum model

A hybrid atomistic-continuum method can model the microstructure evolution of metals subjected to laser irradiation. This method combines classical molecular dynamics (MD) simulations with the two-temperature model (TTM) to account for the laser energy absorption and heat diffusion behavior. Accurate prediction of the temperature evolution in the combined MD-TTM method requires reliable accuracy in electron heat capacity, electron thermal conductivity, and electron–phonon coupling factor across the temperatures generated. This study uses the electronic density of states (DOS) obtained from first-principle calculations. The calculated electron temperature-dependent parameters are used in MD-TTM simulations to study the laser metal interactions in FCC and BCC metals and the phenomenon of laser shock loading and melting. This study uses FCC Al and BCC Ta as model systems to demonstrate this capability. When subjected to short pulsed laser shocks, the dynamic failure behavior predicted using temperature-dependent parameters is compared with the experimentally reported single-crystal and nanocrystalline Al and Ta systems. The MD-TTM simulations also investigate laser ablation and melting behavior of Ta to compare with the ablation threshold reported experimentally. In conclusion, this manuscript demonstrates that integrating the temperature-dependent parameters into MD-TTM simulations leads to the accurate modeling of the laser–metal interaction and allows the prediction of the kinetics of the solid–liquid interface.

36 MATERIALS SCIENCE↗

Sensitivity and vulnerability to summer heat extremes in major cities of the United States

Many cities are experiencing increases in extreme heat because of global temperature rise combined with the urban heat island effect. The heterogeneity of urban morphology also leads to fine-scale variability in potential for heat exposure. Yet, how this rise in temperature and local variability together impacts urban residents differently at exposure-relevant scales is still not clear. Here we map the Universal Thermal Climate Index, a more complete indicator of human heat stress at an unprecedentedly fine spatial resolution (1 m), for 14 major cities in the United States using urban microclimate modeling. We examined the different heat exposure levels across different socioeconomic and racial/ethnic groups in these cities, finding that income level is most consistently associated with heat stress. We further conducted scenario simulations for a hypothetical 1 °C increase of air temperature in all cities. Results show that a 1 °C increase would have a substantial impact on human heat stress, with impacts that differ across cities. The results of this study can help us better evaluate the impact of extreme heat on urban residents at decision-relevant scales.

54 ENVIRONMENTAL SCIENCES↗

Investigations of Mechanisms Associated with Nucleate Boiling Under Microgravity Conditions

In this work a building block type of approach is used so that a basic understanding of the processes that contribute to nucleate boiling heat fluxes under microgravity conditions can be developed. This understanding will lead to development of a mechanistic model for nucleate boiling heat transfer which could eventually be used as a design tool in space applications. Task Description Task 1: Fabrication of the Experimental Setup. Under this task, the test section and liquid holding and viewing chambers will be fabricated. Artificial cylinder cavities will be formed on silicon wafers. A single cavity and two or four cavities, with a prescribed spacing and size, will be formed. The desired nucleation wall superheat will be used to determine the size of the mouth of the cavities. Task 2: Experiments. The basic experiments for flow and temperature field around single and multiple (2 or 4 separated or merged bubbles growing on downward facing or inclined surfaces) will be carried out under normal gravity conditions. The experiments will be conducted at one atmosphere pressure, but liquid subcooling will be varied from 0 to 30C. Water and PF-5050 will be used as test liquids. Task 3: Analytical/Numerical Models. In this task, transient temperature and flow field in vapor and liquid will be determined during growth of a single bubble. Analysis will include the evolution of the vapor-liquid interface and development of microlayer underneath the bubbles. For merged bubbles, detailed calculations of flow and temperature field will be carried out for transient shapes of vapor stems supporting a large bubble and the corresponding evaporation rate. Flow and temperature field for a bubble sliding along a heated wall will also be determined. Microgravity conditions will be simulated and a framework of a numerical tool for prediction of nucleate boiling heat fluxes under microgravity conditions will be developed. Task 4: Experiments in a KC-135. To understand bubble growth and detachment behavior of single or large merged bubbles, boiling experiments will be conducted under low gravity (10-2 g) conditions of the aircraft. In these experiments, 'designed' surfaces will be used. Visual observations and heat transfer data will be taken, but holography will not be used. The apparatus used for laboratory experiments will also be employed for experiments in the aircraft. Task 5: Experiments in the Space Shuttle. Effort will be devoted for defining a boiling experiment to be conducted on a 'designed' surface. The experiment will provide microgravity data on bubble growth and departure. These data are needed for development of a credible model for nucleate boiling heat fluxes under microgravity conditions. The heat transfer data will also be obtained and will be used to validate the models.

Dhir, Vijay K.↗

Interface shape in Czochralski grown crystals - Effect of conduction and radiation

A sequential modular computational scheme has been proposed for the calculation of temperature profiles in the melt and crystal including the melt-crystal interface shape and the crystal pulling rate for a fixed crystal radius. The heat transfer in the melt is assumed to occur by conduction only. The shape of melt-gas meniscus, described by Laplace-Young equation, is incorporated in the model and its effects on interface shape and pulling rate are examined. Further, the model incorporates the detailed radiation interaction among the various surfaces in the puller using the Gebhart enclosure theory which accounts for both direct and reflected radiation. The effects of various process parameters on the interface shape and the pulling rate have been investigated. The results of this study indicate that the shape of the melt-gas meniscus has a significant effect on the pulling rate and the interface shape in conformity with earlier studies. Further, it is shown that a simple radiation model (Stefan's model) is inadequate to model the radiation heat exchange in the Czochralski puller apparatus and the detailed Gebhart analysis is necessary for accurate calculation of both the temperature profile in the crystal as well as the interface shape.

Srivastava, R. K.↗

Coronal loops - Current-based heating processes

Based on new observations, a theoretical model of magnetic-field related heating processes in the solar corona is given. In this model, field-aligned currents are induced along coronal loops in thin current sheaths. Excitation of instabilities involving magnetic reconnection converts the energy associated with the current-related magnetic field directly into particle energy, where the heating process proceeds via short bursts corresponding to an intermittent disruption of the current sheath configuration. Because of the relatively low transverse thermal conduction, only a small fraction of the loop volume is heated to a much higher temperature than the average value. This is consistent with experimental observations of low filling factors of hot plasmas in coronal loops. Thus the model involves a repeated sequence of dynamic events taking into account the observed loop topology, the differential emission measure distribution in the 10 exp 6 - 10 exp 7 K range, the energy balance requirements in the loop, and the probable duty cycles involved in the heating processes.

Beaufume, P.↗

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗

Operating Hardware Impact on the Heat Transfer Properties of Windows

Despite significant advancements in fenestration technology in the last two decades, the thermal transmittance of fenestration products is still significantly higher than that of walls. This corresponds to 60% of the total energy loss of a modern building envelope through the windows. Hence, further development and improvements of fenestration products are necessary. Increasingly stringent codes and standards for fenestration stimulate industry to work on improved solutions. Thus, it is crucial that assessment techniques are able to account for innovations accurately. The thermal effects of non-continuous hardware in window frames are currently ignored by international rating procedures. A preliminary investigation conducted by our team showed significant performance degradation in two of the three out-opening casement profiles caused by the presence of operating hardware. Frames with the structure made of vinyl and fiberglass consist of many air cavities that are penetrated by operating hardware made of highly conductive materials. In these frames, in order to have an accurate assessment, it may be required to employ three-dimensional modeling due to the convective nature of heat transfer within the cavities. However, in this study, we demonstrate that the three-dimensional (3D) effects of non-continuous hardware can be approximated accurately with simpler two-dimensional (2D) simulations. We then develop a simplified model based on weighted average capable of replacing the time- and computation-intensive 3D simulations with 2D simulations and validate it against market available frames and their corresponding hardware. Validation results show that our approximation technique results in discrepancies lower than 0.05 W/(m2K), or 3% of the total thermal transmittance. Thus, we conclude that simplified 2D simulation models may be used for predicting hardware impact in window frames with reasonable accuracy. As windows and glazing structures are becoming ever better thermally insulated, it is becoming even more important to be able to model the impact of the operating hardware on the total thermal performance in order to design the best windows possible and not let the operating hardware ruin an otherwise well-proven design, which is hence addressed in this study.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Gravothermal evolution of dark matter halos with differential elastic scattering

Here, we study gravothermal evolution of dark matter halos in the presence of differential self-scattering that has strong velocity and angular dependencies. We design controlled N-body simulations to model Rutherford and Møller scatterings in the halo, and follow its evolution in both core-expansion and -collapse phases. The simulations show the commonly-used transfer cross section underestimates the effects of dark matter self-interactions, but the viscosity cross section provides an accurate approximation for modeling angular-dependent dark matter scattering. We investigate thermodynamic properties of the halo, and find that the three moments of the Boltzmann equation under the fluid approximation are satisfied. We further propose a constant effective cross section, which integrates over the halo's characteristic velocity dispersion with weighting kernels motivated by kinetic theory of heat conduction. The effective cross section provides a good approximation to differential self-scattering for most of the halo evolution. It indicates that we can map astrophysical constraints on a constant self-interacting cross section to an SIDM model with velocity- and angular-dependent scatterings.

79 ASTRONOMY AND ASTROPHYSICS↗

Effects of the cooling system parameters on heat transfer and performance of the PAFC stack during transient operation

An experimental investigation for the effects of transient operation of a phosphoric acid fuel-cell stack on heat transfer and temperature distribution in the electrodes has been conducted. The proposed work utilized the experimental setup with modifications, which was designed and constructed under NASA Contract No. NCC-3-17(5). The experimental results obtained from this investigation and the mathematical model obtained under NASA Contract No. NCC3-17(4) after modifications, were utilized to develop mathematical models for transient heat transfer coefficient and temperature distribution in the electrode and to evaluate the performance of the cooling - system under unsteady state conditions. The empirical formulas developed were then implemented to modifying the developed computer code. Two incompressible coolants were used to study experimentally the effect of the thermophysical properties of the cool-ants on the transient heat transfer coefficient and the thermal contact resistance during start-up and shut-down processes. Coolant mass flow rates were verified from 16 to 88.2 Kg/hr during the transient process when the electrical power supply was gradually increased or decreased in the range (O to 3000 W/sq m). The effect of the thermal contact resistance with a range of stack pressure from O to 3500 KPa was studied.

Ridha, Rabi M. J.↗

Nonlinear dynamic analysis of solution multiplicity of buoyancy ventilation in a typical underground structure

Buoyancy ventilation is widely used in underground buildings, such as underground hydropower stations. Multiple solutions of buoyancy ventilation may exist in those underground structures. In this study, we developed a transient model comprising an ordinary differential equation system to describe buoyancy ventilation patterns in typical two-zone underground structures. Additionally, the accuracy of the model was validated. Nonlinear dynamical analysis was conducted to study multiple steady-state airflow. According to mathematical derivation, the configuration of one local heat source at the bottom corner introduces two stable solutions. The criterion to determine the stability and existence of solutions for more general scenarios was developed. Using this criterion, we obtained the multiple steady states of any two-zone underground buildings for different stack height ratios and the strength ratios of the heat sources. Furthermore, this criterion can be adopted for the design of buoyancy ventilation or natural smoke ventilation systems. Designers can change the height ratio of the stack or the heat ratio of two zones to induce the desired ventilation patterns. Finally, a case study was conducted with field measurements to demonstrate the use of the nonlinear dynamical analysis method to investigate the multiple steady states of buoyancy ventilation. Through the case study, we validated that the proposed criterion could produce the same result as the nonlinear dynamical analysis.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Hydrogen Risk Assessment Models (HyRAM) (Version 3.0 Technical Reference Manual)

The HyRAM software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen infrastructure and transportation systems. HyRAM is designed to facilitate the use of state-of-the-art science and engineering models to conduct robust, repeatable assessments of hydrogen safety, hazards, and risk. HyRAM includes generic probabilities for hydrogen equipment failures, probabilistic models for the impact of heat flux on humans and structures, and computationally and experimentally validated first-order models of hydrogen release and flame physics. HyRAM integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hydrogen hazards (thermal effects from jet fires, overpressure effects from deflagrations), and assessing impact on people and structures. HyRAM is developed at Sandia National Laboratories for the U.S. Department of Energy to increase access to technical data about hydrogen safety and to enable the use of that data to support development and revision of national and international codes and standards. HyRAM is a research software in active development and thus the models and data may change. This report will be updated at appropriate developmental intervals. This document provides a description of the methodology and models contained in the HyRAM version 3.0. HyRAM 3.0 includes the new ability to model cryogenic hydrogen releases from liquid hydrogen systems, using a different property calculation method and different equations of state. Other changes include modifications to the ignition probability calculations, component leak frequency calculations, and addition of default impulse data.

08 HYDROGEN↗

Hydrogen Risk Assessment Models (HyRAM) (V.3.1) (Technical Reference Manual)

The HyRAM software toolkit provides a basis for conducting quantitative risk assessment and consequence modeling for hydrogen infrastructure and transportation systems. HyRAM is designed to facilitate the use of state-of-the-art science and engineering models to conduct robust, repeatable assessments of hydrogen safety, hazards, and risk. HyRAM includes generic probabilities for hydrogen equipment failures, probabilistic models for the impact of heat flux on humans and structures, and experimentally validated first-order models of hydrogen release and flame physics. HyRAM integrates deterministic and probabilistic models for quantifying accident scenarios, predicting physical effects, and characterizing hydrogen hazards (thermal effects from jet res, overpressure effects from deflagrations), and assessing impact on people and structures. HyRAM is developed at Sandia National Laboratories for the U.S. Department of Energy to increase access to technical data about hydrogen safety and to enable the use of that data to support development and revision of national and international codes and standards. HyRAM is a research software in active development and thus the models and data may change. This report will be updated at appropriate developmental intervals. This document provides a description of the methodology and models contained in HyRAM version 3.1. There have been several impactful updates since version 3.0. HyRAM 3.1 contains a correction to use the volume fraction for two-phase speed of sound calculations; this only affects cryogenic releases in which two-phase ow (vapor and liquid) is predicted in the orifice. Other changes include clarifications that inputs for tank pressure should be given in absolute pressure, not gauge pressure. Additionally, the interface now rejects invalid inputs to probability distributions, and the less accurate single-point radiative source model selection was removed from the interface.

08 HYDROGEN↗

Mathematical modeling of moving boundary problems in thermal energy storage

The capability for predicting the performance of thermal energy storage (RES) subsystems and components using PCM's based on mathematical and physical models is developed. Mathematical models of the dynamic thermal behavior of (TES) subsystems using PCM's based on solutions of the moving boundary thermal conduction problem and on heat and mass transfer engineering correlations are also discussed.

Solomon, A. D.↗

Heat generation in aircraft tires under braked rolling conditions

An analytical model was developed to approximate the internal temperature distribution in an aircraft tire operating under conditions of unyawed braked rolling. The model employs an array of elements to represent the tire cross section and considers the heat generated within the tire to be caused by the change in strain energy associated with cyclic tire deflection. The additional heating due to tire slip and stresses induced by braking are superimposed on the previously developed free rolling model. An extensive experimental program was conducted to verify temperatures predicted from the analytical model. Data from these tests were compared with calculations over a range of operating conditions. The model results were in reasonably good agreement with measured values.

Clark, S. K.↗

A variational formulation for time-dependent climate models

A variational principle for time dependent diffusion problems is presented and is demonstrated by applying it to simple seasonal climate models. Two cases are treated. The first, a North-Coakley-type model with constant properties, is used as a tutorial example for the application of the technique. For the second case, heat capacity and thermal conductivity are considered to be latitude dependent in order to treat the effects of land/ocean distribution on the seasonal temperature distribution over the earth. The variational equations are derived and approximate analytical solutions are developed which delineate the influences of the physical asymmetries of the hemispheres in producing an asymmetric annual cycle.

Smith, G. L.↗