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,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Large Eddy Simulation of Low-Reynolds-Number Turbulent Flow of Low-Prandtl-Number Fluid in a Tight Lattice Bundle for Assessment of Reynolds-Averaged Navier-Stokes Turbulence Model

The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2021

Lattice Boltzmann model for conjugate heat transfer across thin walls

A lattice Boltzmann (LB) model with an efficient and accurate interface treatment for conjugate heat transfer across a thin wall between two different media is developed. The proposed interface treatment avoids fine meshing and computation within the thin layer; instead, the energy balance within the thin layer and the conjugate conditions on each interface are utilized to construct explicit updating schemes for the microscopic distribution functions of the LB model at the interior lattice nodes of the two media next to the thin layer. The proposed interface schemes reduce to the standard interface scheme for conjugate conditions in the literature in the limit of zero thickness of the thin layer, and thus it can be considered a more general interface treatment. A simplified version of the interface treatment is also proposed when the heat flux variation along the tangential direction of the thin layer is negligible. Three representative numerical tests are conducted to verify the applicability and accuracy of the proposed interface schemes. The results demonstrate that the intrinsic second-order accuracy of the LB model is preserved with the proposed interface schemes for thin layers with constant tangential fluxes, while for general situations with varying tangential fluxes, first-order accuracy is obtained. Finally, this interface treatment within the LB framework is attractive in conjugate heat transfer modeling involving thin layers for its simplicity, accuracy, and significant reduction in computational resources.

42 ENGINEERING↗
DOE OSTI2021

Direct Observation of the Depth of Active Groundwater Circulation in an Alpine Watershed

Abstract The depth of active groundwater circulation is a fundamental control on stream flows and chemistry in mountain watersheds, yet it remains challenging to characterize and is rarely well constrained. We collected hydraulic conductivity, hydraulic head, temperature, chemical, noble gas, and 3 H/ 3 He groundwater age data from discrete levels in two boreholes 46 and 81 m deep in an alpine watershed, in combination with chemical and age data from shallow groundwater discharge, to discern groundwater flow rates at different depths and directly observe active and inactive groundwater. Vertical head gradients are steep (average of 0.4) and thermal profiles are consistent with typical linear conductive continental geotherms. Groundwater deeper than ∼20 m is distinct from shallow groundwater and creek water in its chemistry, noble gas signature, and age (dominantly >65 years compared to <9 years). Together these results suggest low vertical groundwater flow velocities and a relatively shallow active circulation depth of ∼20 m. This hypothesis is tested with a simple 2‐D numerical fluid flow and heat transport model representing a hillslope transect through the two boreholes. The modeling indicates that the subhorizontally bedded sedimentary rocks underlying the basin are highly anisotropic with low vertical hydraulic conductivity, and at most ∼10% of bedrock recharge (equivalent to <2% of stream baseflow) flows below a depth of 20 m. The study demonstrates the considerable value of discrete‐depth hydrogeologic, chemical, and age data for determining active circulation depth, and illustrates an approach for maximizing the utility of individual boreholes drilled for mountain bedrock aquifer characterization.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2000

Downward Link of Solar Activity Variations Through Wave Driven Equatorial Oscillations (QBO and SAO)

Signatures of the 11-year solar activity/irradiance cycle are observed in the Quasi Biennial Oscillation (QBO) of the lower stratosphere. At these altitudes, the QBO is understood to be the result of "downward control" exerted by the wave mean flow interactions that drive the phenomenon. It is reasonable then to speculate that the QBO is a natural conduit to lower altitudes of solar activity variations in radiance (SAV). To test this hypothesis, we conducted experiments with a 2D version of our Numerical Spectral Model that incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). To account for the SAV, we change the solar heating rate on a logarithmic scale from 0.1% at the surface to 1% at 50 kin to 10% at 100 km. With the same GW flux, we then conduct numerical experiments to evaluate the magnitude of the solar activity irradiance effect (SAE) on the zonal circulation at low latitudes. The numerical results obtained show that, under certain conditions, the SAE is significant in the zonal circulation and does extend to lower altitudes where the SAV is small. The differences in the wind velocities can be as large as 5 m/s at 20 kin. We carried out two numerical experiments with integrations over more than 20 years: 1) With the QBO period "tuned" to be 30 months, of academic interest but instructive, the seasonal cycle in the solar forcing [through the Semi-annual Oscillation (SAO)] acts as a strong pacemaker to produce a firm lock on the period and phase of the QBO. The SAE then shows up primarily as a distinct but relatively weak amplitude modulation. 2) With the QBO period between 30 and 34 (or less than 30, presumably) months, the seasonal phase lock is weak compared with (1). The SAV in the seasonal cycle then causes variations in the QBO period and phase, and this amplifies the SAE to produce relatively large variations in the wind field. We conclude that, under realistic conditions as in (2), the solar seasonal forcing, with activity variations in radiance, causes the phase and period of the QBO to change to produce a relatively large solar activity effect in the zonal circulation at low latitudes and low altitudes.

Mengel, J. G.↗
DOE OSTI2022

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2016

Weather and Climate Change Impacts on Human Mortality in Bangladesh

Weather and climate profoundly affect human health. Several studies have demonstrated a U-, V-, or J-shaped temperature-mortality relationship with increasing death rates at the lower and particularly upper end of the temperature distribution. The objectives of this study were (1) to analyze the relationship between temperature and mortality in Bangladesh for different subpopulations and (2) to project future heat-related mortality under climate change scenarios. We used (non-)parametric Generalized Additive Models adjusted for trend, season and day of the month to analyze the effect of temperature on daily mortality. We found a decrease in mortality with increasing temperature over a wide range of values; between the 90th and 95th percentile an abrupt increase in mortality was observed which was particularly pronounced for the elderly above the age of 65 years, for males, as well as in urban areas and in areas with a high socio-economic status. Daily historical and future temperature values were obtained from the NASA Earth Exchange Global Daily Downscaled Projections (NEX-GDDP) dataset. This dataset is comprised of downscaled climate scenarios for the globe that are derived from the General Circulation Model (GCM) runs conducted under the Coupled Model Intercomparison Project Phase 5 (CMIP5). The derived dose-response functions were used to estimate the number of heat-related deaths occurring during the 1990s (1980-2005), the 2020s (2011-2040) and the 2050s (2041-2070). We estimated that excess deaths due to heat will triple from the 1990s to the 2050s, with an annual number of 0.5 million excess deaths in 1990 to and expected number of 1.5 millions in 2050.

mortality↗
NASA NTRS2011

Design for On-Sun Evaluation of Evaporator Receivers

A heat pipe designed for operation as a solar power receiver should be optimized to accept the solar energy flux and transfer this heat into a reactor. Optical properties of the surface, thermal conductance of the receiver wall, contact resistance of the heat pipe wick, and other heat pipe wick properties ultimately define the maximum amount of power that can be extracted from the concentrated sunlight impinging on the evaporator surface. Modeling of solar power receivers utilizing optical and physical properties provides guidance to their design. On-sun testing is another important means of gathering information on performance. A test rig is being designed and built to conduct on-sun testing. The test rig is incorporating a composite strip mirror concentrator developed as part of a Small Business Innovative Research effort and delivered to NASA Glenn Research Center. In the strip concentrator numerous, lightweight composite parabolic strips of simple curvature were combined to form an array 1.5 m x 1.5 m in size. The line focus of each strip is superimposed in a central area simulating a point of focus. A test stand is currently being developed to hold the parabolic strip concentrator, track the sun, and turn the beam downward towards the ground. The hardware is intended to be sufficiently versatile to accommodate on-sun testing of several receiver concepts, including those incorporating heat pipe evaporators. Characterization devices are also being developed to evaluate the effectiveness of the solar concentrator, including a receiver designed to conduct calorimetry. This paper describes the design and the characterization devices of the on-sun test rig, and the prospect of coupling the concentrated sunlight to a heat pipe solar power receiver developed as part of another Small Business Innovative Research effort.

Jaworske, Donald A.↗
NASA NTRS1998

Effects of Ice Accretion on Aircraft Aerodynamics

The primary objective of this research was to support the development of a new ice accretion model by improving our physical understanding of the ice accretion process through experimental measurements. The focus was on the effect of the initial ice roughness (smooth/rough boundary) on the accretion process. This includes understanding the boundary-layer development over the roughness and especially its effect on the heat transfer which is fundamental to the ice accretion process. The research focused on acquiring the experimental data needed to formulate a new ice accretion physical model. Research was conducted to analyze boundary-layer data taken on a NACA 0012 airfoil with roughness to simulate the smooth/rough boundary. The effect of isolated roughness on boundary-layer transition was studied experimentally to determine if the classical critical roughness Reynolds number criteria could be applied to transition in the airfoil leading-edge area. The effect of simulated smooth/rough boundary roughness on convective heat transfer was studied to complete the study. During the course of this research the effect of free-stream wind tunnel turbulence on the boundary layer was measured. Since this quantity was not well known, research to accurately measure the wind tunnel turbulence in an icing cloud was undertaken. Preliminary results were attained and the final data were acquired, reduced and presented under a subsequent grant.

Bragg, Michael B.↗
NASA NTRS1994

The effect of gravity waves on the global mean temperature and composition structure of the upper atmosphere

Formulas are presented that parameterize the heating rate and coefficient of turbulent heat conduction produced by saturated internal gravity waves (IGW) in the upper atmosphere. Estimates of these values are made using observational data. The parameterization of IGW influences are introduced into a one-dimensional model of global mean thermal and composition balances of the upper atmosphere. Computations are performed for different values of IGW energy fluxes entering into the upper atmosphere from below. It is shown that realistic vertical profiles of the global mean temperature can be obtained using different values of IGW energy flux into the upper atmosphere. Increasing the IGW intensity leads not only to an increase of the heating rate due to wave enery dissipation, but also to an increase of the heating rate due to wave energy dissipation, but also to an increase in the coefficient of turbulent heat conduction and cooling rate produced by turbulence generated by the wave. So, near an altitude of 100 km the main part of solar heating is compensated by infrared cooling on one hand, and the main part of wave dissipation heating is compensated by turbulent cooling on the other hand. These quasi-balances generally hold for different values of IGW intensity.

Gavrilov, Nikolai M.↗
NASA NTRS1981

Satellite estimation of the surface energy balance, moisture availability and thermal inertia

A method for inferring the distribution of surface heat and evaporative fluxes and the ground moisture availability and thermal inertia (ground conductive capacity) is used to analyze two urbanized areas, Los Angeles and St. Louis. The technique employs infrared satellite temperature measurements in conjunction with a one-dimensional boundary-layer model. Results show that there is a marked reduction of evaporation and moisture availability and a corresponding elevation of sensible heat flux over urbanized areas and over cropped areas with low vegetative cover. Conversely, low heat flux and high evaporation characterize vegetated and, especially, forested areas. Warm urban centers appear directly related to a reduction in vegetation, which normally allows for a greater fraction of available radiant energy to be converted into latent heat flux. The distribution of thermal inertia was surprisingly ill-defined and its variation between urban and rural areas was quite small. Thus, the increased heat storage within the urban fabric, which has been proposed as the underlying cause of the nocturnal heat island, may be caused mainly by enhanced daytime surface heating which occurs because of surface dryness, rather than by large spatial variations in the conductivity of the surface.

Carlson, T. N.↗
DOE OSTI2026

Working Fluid Characterization and Performance Assessment of Subcritical Organic Rankine Cycles Based on the Lee–Kesler Approach for Energy Recovery

Here, a generalized model using the Lee–Kesler approach based on the corresponding states principle is developed to assess the performance of subcritical Organic Rankine Cycles operating with different working fluids. Each fluid is characterized by five parameters: the acentric factor, critical temperature, critical pressure, molar mass, and the ideal-gas ratio of specific heats at the critical temperature. The model was developed using the compressibility factor modified version of the Benedict–Webb–Rubin equation proposed by Lee and Kesler and the enthalpy and entropy functions to calculate thermodynamic state properties. The model was validated by comparing the results calculated with the model and working fluid thermodynamic properties obtained with the CoolProp database. This comparison was conducted for 91 working fluids, obtaining a relative error below 5% for 88 out of the 91 fluids (∼97%). A generalized parametric study was conducted to determine the influence of the pinch point and each fluid parameter on the performance of Organic Rankine Cycle (ORC) systems. It was found that efficiency increases with critical temperature, ideal-gas ratio of specific heats at the critical temperature, and acentric factor, reaching up to 13%. The developed model enables the evaluation of ORC system performance for existing working fluids. It also allows the formulation and evaluation of new fluids to enhance the performance of the ORC while retrieving energy from any kind of source; and likewise, the methodology can be applied to other power generation cycles.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS1996

Heat Transfer Experiments in the Internal Cooling Passages of a Cooled Radial Turbine Rotor

An experimental study was conducted (1) to experimentally measure, assess and analyze the heat transfer within the internal cooling configuration of a radial turbine rotor blade and (2) to obtain heat transfer data to evaluate and improve computational fluid dynamics (CFD) procedures and turbulent transport models of internal coolant flows. A 1.15 times scale model of the coolant passages within the NASA LERC High Temperature Radial Turbine was designed, fabricated of Lucite and instrumented for transient beat transfer tests using thin film surface thermocouples and liquid crystals to indicate temperatures. Transient heat transfer tests were conducted for Reynolds numbers of one-fourth, one-half, and equal to the operating Reynolds number for the NASA Turbine. Tests were conducted for stationary and rotating conditions with rotation numbers in the range occurring in the NASA Turbine. Results from the experiments showed the heat transfer characteristics within the coolant passage were affected by rotation. In general, the heat transfer increased and decreased on the sides of the straight radial passages with rotation as previously reported from NASA-HOST-sponsored experiments. The heat transfer in the tri-passage axial flow region adjacent to the blade exit was relatively unaffected by rotation. However, the heat transfer on one surface, in the transitional region between the radial inflow passage and axial, constant radius passages, decreased to approximately 20 percent of the values without rotation. Comparisons with previous 3-D numerical studies indicated regions where the heat transfer characteristics agreed and disagreed with the present experiment.

Johnson, B. V.↗
NASA NTRS2023

Importance of Atmospheric Feedbacks in Simulating the Seasonal Cycle of the Antarctic Sea Ice and Its Response to Perturbations

The seasonal cycle of the Antarctic sea ice extent is largely controlled by the evolution of the insolation received at the top of the atmosphere. However, sea ice processes and feedbacks with the ocean and the atmosphere can modulate this seasonal cycle. Here, the atmospheric feedbacks are quantified in a series of idealized sensitivity experiments performed with an eddy-permitting (1/4°) NEMO-LIM3 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-CLM 2 regional atmospheric model. In these 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/release upon sea ice growth/melt, ice dynamics and surface albedo. We show that the changes in surface heat fluxes are very different between the configurations driven by the reanalysis and those coupled to the atmosphere. Atmospheric feedbacks enhance the response of the modelled winter ice extent 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. By contrast, the atmospheric feedbacks can damp the impact of the perturbations affecting the heat conductivity through sea ice.

Atmospheric feedbacks↗
NASA NTRS1977

Results of phase change heat transfer test OH51 using 0.006-scale space shuttle orbiter models 46-0 and 90-0 and partial wing 0.0175-scale model 64-0 in the LaRC 31-inch CFHT

Test OH51 was a phase change paint test conducted in the LaRC 31-inch CFHT utilizing models 46-0, 64-0, and 90-0. Model 46-0 represented the space shuttle configuration 139 Orbiter. Model 90-0 represented the configuration 140 Orbiter. Model 64-0 represented the forward 45% portion of the Orbiter wing. The partial wing was tested with a shock generator located at various positions relative to the wing. The test was conducted at Mach 10.0, angles of attack from 27.5 deg through 37.5 deg, and Reynolds numbers of 0.5 and 1.5 million per foot.

Cummings, J. W.↗
DOE OSTI2021

Specifics of Thermal Transport in Graphene Composites: Effect of Lateral Dimensions of Graphene Fillers

In this work, we report on the investigation of thermal transport in noncured silicone composites with graphene fillers of different lateral dimensions. Graphene fillers are comprised of few-layer graphene flakes with lateral sizes in the range from 400 to 1200 nm and the number of atomic planes from 1 to ~100. The distribution of the lateral dimensions and thicknesses of graphene fillers has been determined via atomic force microscopy statistics. It was found that in the examined range of the lateral dimensions, the thermal conductivity of the composites increases with increasing size of the graphene fillers. The observed difference in thermal properties can be related to the average gray phonon mean free path in graphene, which has been estimated to be around ~800 nm at room temperature. The thermal contact resistance of composites with graphene fillers of 1200 nm lateral dimensions was also smaller than that of composites with graphene fillers of 400 nm lateral dimensions. The effects of the filler loading fraction and the filler size on the thermal conductivity of the composites were rationalized within the Kanari model. The obtained results are important for the optimization of graphene fillers for applications in thermal interface materials for heat removal from high-power-density electronics.

36 MATERIALS SCIENCE↗
NASA NTRS1980

Dynamic spectral characteristics of thermal models for solar hard X-ray bursts

It is proposed that the source for solar hard X-ray bursts may consist of a distribution of many small impulsively-heated kernels, each cooled by anomalous conduction, with lifetimes shorter than current burst data temporal resolution. The dynamic spectra of bursts are governed by the dynamic evolution of the kernel production process, such as magnetic-field dissipation in the tearing mode. An integral equation is formulated, the solution of which yields information on this kernel production process from dynamic burst spectra, for any kernel model.

Brown, J. C.↗
NASA NTRS2016

Hall Thruster Thermal Modeling and Test Data Correlation

HERMeS - Hall Effect Rocket with Magnetic Shielding. Developed through a joint effort by NASA/GRC and the Jet Propulsion Laboratory (JPL). Design goals: High power (12.5 kW) high Isp (3000 sec), high efficiency (> 60%), high throughput (10,000 kg), reduced plasma erosion and increased life (5 yrs) to support Asteroid Redirect Robotic Mission (ARRM). Further details see "Performance, Facility Pressure Effects and Stability Characterization Tests of NASAs HERMeS Thruster" by H. Kamhawi and team. Hall Thrusters (HT) inherently operate at elevated temperatures approx. 600 C (or more). Due to electric magnetic (E x B) fields used to ionize and accelerate propellant gas particles (i.e., plasma). Cooling is largely limited to radiation in vacuum environment.Thus the hardware components must withstand large start-up delta-T's. HT's are constructed of multiple materials; assorted metals, non-metals and ceramics for their required electrical and magnetic properties. To mitigate thermal stresses HT design must accommodate the differential thermal growth from a wide range of material Coef. of Thermal Expansion (CTEs). Prohibiting the use of some bolted/torqued interfaces.Commonly use spring loaded interfaces, particularly at the metal-to-ceramic interfaces to allow for slippage.However most component interfaces must also effectively conduct heat to the external surfaces for dissipation by radiation.Thus contact pressure and area are important.

Model↗
DOE OSTI2022

COMPASS-U Global Heat Balance Calculations

COMPASS-U is a medium size, high magnetic field experimental tokamak (R = 0.9 m, B t = 5 T, and I p = 2 MA), built at the Institute of Plasma Physics, Czech Academy of Sciences (IPP-CAS). This global heat balance calculation for COMPASS-U was done at Princeton Plasma Physics Laboratory (PPPL). Based on our previous experience of building global thermal model for National Spherical Torus Experiment-Upgrade (NSTX-U) at PPPL, this 2-D global thermal model geometry represents a typical cross section of COMPASS-U machine. The model includes thermal radiation, conduction, and convection among components and also between the machine and outer environment. Helium gas heating/cooling was modeled with fluid element and surface element. This model was used to calculate component temperatures and heat distribution during heat up, cryogenic cool down, normal operation, and fault operation scenarios. Definition of 14 main thermal scenarios was provided by IPP-CAS. First nine are normal operation scenarios. Scenarios #10–#14 are fault scenarios. Results of thermal scenario #5 will be given and discussed in this article, including peak temperatures, temperature ratcheting, energy distribution, and required cooling power.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗