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 847 records · Page 47

Three-Dimensional CFD Investigation of Pre-Spark Heat Release in a Boosted SI Engine

Low-temperature heat release (LTHR) in spark-ignited internal combustion engines is a critical step toward the occurrence of auto-ignition, which can lead to an undesirable phenomenon known as engine knock. As such, correct predictions of LTHR are of utmost importance to improve the understanding of knock and enable techniques aimed at controlling it. While LTHR is typically obscured by the deflagration following the spark ignition, extremely late ignition timings can lead to LTHR occurrence prior to the spark, i.e., pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection SI engine was numerically investigated using three-dimensional computational fluid dynamics (CFD). A hybrid approach was used, based on the G-equation model for representing the turbulent flame front and the multi-zone well-stirred reactor model for tracking the chemical reactions within the unburnt region. A recently developed best practice was also employed which keeps the well-stirred reactor model active throughout the entire simulation. This allowed for correct predictions of the previous cycle trapped residuals which have a considerable effect on the onset of LTHR. Multi-cycle simulations were conducted using Co-Optima alkylate and E30 fuels. The predicted in-cylinder pressure and heat release rate agreed well with the experimental data and served as validation for the CFD model. Following the initial validation, the dynamics of PSHR was discussed and a series of parameters of interest were assessed. First, the effect of exhaust valve temperature was investigated, qualitatively highlighting the importance of boundary conditions uncertainty. Further analyses were carried out on the effects of fuel properties, including laminar flame speed (LFS) and heat of vaporization (HOV). The results indicate that PSHR phasing is slightly advanced with lower LFS as more trapped unburnt fuel is made available for the next cycle. A similar trend was observed with lower HOV as less intense spray cooling led to higher mixture temperatures, i.e., higher mixture reactivity. Finally, a comparison of Co-Optima alkylate and E30 fuels was made using the pressure-temperature trajectory framework. It was shown that the differences between Co-Optima alkylate and E30 in PSHR tendency are correlated with both HOV and chemical effects.

33 ADVANCED PROPULSION SYSTEMS↗

Numerical analysis of particle dispersion and deposition in coal combustion using large-eddy simulation

High-fidelity modeling provides a useful approach to investigate the particle dispersion and deposition mechanism in pulverized coal combustion. To be able to analyze these detailed mechanisms, this work couples detailed multiphysics models for poly-dispersed turbulent, particle-laden flow, particle and gas phase reaction chemistry, convective, conductive and radiative heat transfer, ash formation and deposition with a Large-Eddy Simulation (LES) approach to numerically simulate coal combustion in a downfired self-sustained oxy-fuel combustor (OFC). It is necessary to explicitly capture all but the highest frequency dynamics of the turbulence and its coupling with each of the other physical phenomena using LES. Effects of subgrid-scale unresolved turbulence on the particle motions are also analyzed by a novel Stokes number analysis. Due to the inherent relation ship between coal combustion and ash deposition, this study integrates the improved ash deposition model into the numerical simulation. Overall simulation results are compared with experimentally measured data from the OFC. The simulation and measured data for the averaged gas temperature and deposition rates agree with 5% and 28%. This study shows that high-fidelity LES coupled with other detailed multiphysics models running on a exascale computing facility can provide a good representation of complex coal combustion and deposition in a laboratory-scale furnace.

multiphysics model↗

Compatibility of Segments of Thermoelectric Generators

A method of calculating (usually for the purpose of maximizing) the power-conversion efficiency of a segmented thermoelectric generator is based on equations derived from the fundamental equations of thermoelectricity. Because it is directly traceable to first principles, the method provides physical explanations in addition to predictions of phenomena involved in segmentation. In comparison with the finite-element method used heretofore to predict (without being able to explain) the behavior of a segmented thermoelectric generator, this method is much simpler to implement in practice: in particular, the efficiency of a segmented thermoelectric generator can be estimated by evaluating equations using only hand-held calculator with this method. In addition, the method provides for determination of cascading ratios. The concept of cascading is illustrated in the figure and the definition of the cascading ratio is defined in the figure caption. An important aspect of the method is its approach to the issue of compatibility among segments, in combination with introduction of the concept of compatibility within a segment. Prior approaches involved the use of only averaged material properties. Two materials in direct contact could be examined for compatibility with each other, but there was no general framework for analysis of compatibility. The present method establishes such a framework. The mathematical derivation of the method begins with the definition of reduced efficiency of a thermoelectric generator as the ratio between (1) its thermal-to-electric power-conversion efficiency and (2) its Carnot efficiency (the maximum efficiency theoretically attainable, given its hot- and cold-side temperatures). The derivation involves calculation of the reduced efficiency of a model thermoelectric generator for which the hot-side temperature is only infinitesimally greater than the cold-side temperature. The derivation includes consideration of the ratio (u) between the electric current and heat-conduction power and leads to the concept of compatibility factor (s) for a given thermoelectric material, defined as the value of u that maximizes the reduced efficiency of the aforementioned model thermoelectric generator.

Snyder, G. Jeffrey↗

Skylab extravehicular mobility unit thermal simulator

The analytical methods, thermal model, and user's instructions for the Skylab Extravehicular Mobility Unit (SEMU) routine are presented. This digital computer program was developed for detailed thermal performance predictions of the SEMU on the NASA-JSC Univac 1108 computer system. It accounts for conductive, convective, and radiant heat transfer as well as fluid flow and special component characterization. The program provides thermal performance predictions for a 967 node thermal model in one thirty-sixth (1/36) of mission time when operated at a calculating interval of three minutes (mission time). The program has the operational flexibility to: (1) accept card or magnetic tape data input for the thermal model describing the SEMU structure, fluid systems, crewman and component performance, (2) accept card and/or magnetic tape input of internally generated heat and heat influx from the space environment, and (3) output tabular or plotted histories of temperature, flow rates, and other parameters describing system operating modes.

Hixon, C. W.↗

Candle Flames in Non-Buoyant Atmospheres

This paper addresses the behavior of a candle flame in a long-duration, quiescent microgravity environment both on the space Shuttle and the Mir Orbiting Station. On the Shuttle, the flames became dim blue after an initial transient where there was significant yellow (presumably soot) in the flame. The flame lifetimes were typically less than 60 seconds. The safety-mandated candlebox that contained the candle flame inhibited oxygen transport to the flame and thus limited the flame lifetime. The flames on the Mir were similar, except that the yellow luminosity persisted longer into the flame lifetime because of a higher initial oxygen concentration, The Mir flames burned for as long as 45 minutes. The difference in the flame lifetime between the Shuttle and Mir flames was primarily the redesigned candlebox that did not inhibit oxygen transport to the flame. In both environments, the flame intensity and the height-to-width ratio gradually decreased as the ambient oxygen content in the sealed chamber slowly decreased. Both sets of experiments showed spontaneous, axisymmetric flame oscillations just prior to extinction. The paper also presents a numerical model of a candle flame. The formulation is two-dimensional and time-dependent in the gas phase with constant specific heats, thermal conductivity and Lewis number (although different species can have different Lewis numbers), one-step finite-rate kinetics, and gas-phase radiative losses from CO2 and H2O. The treatment of the liquid/wick phase assumes that the, fuel evaporates from a constant diameter sphere connected to an inert cone. The model predicts a steady flame with a shape and size quantitatively similar to the Shuttle and Mir flames. The computation predicts that the flame size will increase slightly with increasing ambient oxygen mole fraction. The model also predicts pre-extinction flame oscillations if the rate of decrease in ambient oxygen is small enough, such as that which would occur for a flame burning in a sealed ambient.

Dietrich, D. L.↗

Review of thermophysical property methods applied to fueled and un-fueled molten salts

Recent investigation into designs of molten salt reactors (MSR) requires accurate data of candidate salts. Review of the literature shows lacking thermophysical (TP) property data for coolants and fuel-bearing salts of Advanced High Temperature Reactors (AHTR). Measurement uncertainties of current data are even more deficient, and this affects ongoing research where MSR modeling and testing are being performed. Melting temperature, viscosity, density, specific heat capacity, and thermal and electrical conductivity measurements are often not arranged together or discussed in great detail, especially for fuel-bearing salts. Raw data from key findings spanning over 70 years for LiF, LiF-BeF 2 , LiF-NaF-KF, KCl-MgCl 2 , KF-ZrF 4 , and NaF-KF-MgF 2 have been organized and presented here. Brief description of methods used to obtain these results, in addition to the TP data of some of these salts as fuel carriers, has also been discussed. Thermal conductivity data are particularly weak due to older, high-error measurement techniques. More recent measurement systems to evaluate salt thermal conductivity show promise for obtaining consistent absolute values. This work concludes that measurement systems such as transient hot-wire method, forced Rayleigh scattering, and laser flash analysis provide the most promise to fill these data gaps and reduce current uncertainty in thermal conductivity measurements of both fueled and un-fueled molten salts.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evolution of magnetic field induced ordering in the layered quantum Heisenberg triangular-lattice antiferromagnet Ba 3 CoSb 2 O 9

Quantum fluctuations in the effective spin-12 layered triangular-lattice quantum Heisenberg antiferromagnet Ba 3 CoSb 2 O 9 lift the classical degeneracy of the antiferromagnetic ground state in magnetic field, producing a series of novel spin structures for magnetic fields applied within the crystallographic ab plane, including a celebrated collinear “up-up-down” spin ordering with magnetization equal to 1/3 of the saturation magnetization over an extended field range. Theoretically unresolved, however, are the effects of interlayer antiferromagnetic coupling and transverse magnetic fields on the ground states of this system. Additional magnetic field induced phase transitions are theoretically expected and in some cases have been experimentally observed, but details regarding their number, location, and physical character appear inconsistent with the predictions of existing models. Conversely, an absence of experimental measurements as a function of magnetic-field orientation has left other key predictions of these models untested. To address these issues, in this study we have used specific heat, neutron diffraction, thermal conductivity, and magnetic torque measurements to map out the phase diagram as a function of magnetic field intensity and orientation relative to the crystallographic ab plane. For H||ab, we have discovered an additional magnetic field induced phase transition at low temperature and an unexpected tetracritical point in the high-field phase diagram, which coupled with the apparent second-order nature of the phase transitions eliminates several theoretically proposed spin structures for the high-field phases. Our calorimetric measurements as a function of magnetic field orientation are in general agreement with theory for field-orientation angles close to plane parallel (H||a) but diverge at angles near plane perpendicular; a predicted convergence of two phase boundaries at finite angle and a corresponding change in the order of the field induced phase transition are not observed experimentally. Our results emphasize the role of interlayer coupling in selecting and stabilizing field induced phases, provide guidance on the nature of the magnetic order in each phase, and reveal the need for new physics to account for the nature of magnetic ordering in this archetypal two-dimensional spin-12 triangular-lattice quantum Heisenberg antiferromagnet.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electrification Options for Multi-Family Water Heating in Cold Climates - Final Report

In multi-family buildings, large water storage tanks in centralized domestic hot water (DHW) systems can serve as thermal energy storage (TES) batteries to mitigate grid impact. These systems offer demand shift and efficiency benefits, significantly reducing peak power consumption, particularly in cold climates. This study evaluates the load-shifting benefits of a centralized heat pump water heater (HPWH) system equipped with a CO2 heat pump in multi-family buildings through simulation. The heat pump system and water storage tank are sized using design-day sizing. A finite-element-based stratified tank model and CO2 heat pump performance map from a commercial DHW product are used. Annual simulations are conducted to assess the benefits of the centralized DHW system for energy efficiency improvements, load shifting, and emission reductions. These simulations incorporate utility tariffs and marginal grid emission data from Los Angeles and Chicago. In Los Angeles, using a water tank as a thermal battery achieves 7.4% utility cost savings and 10.2% emission reduction. In Chicago, compared to HPWH conventional operation without preheating, TES-enabled central HPWH provides 15% utility cost savings and 13% emission reduction. The case study demonstrates that the demand reduction potential of central CO2 HPWHs is significant in cold climate regions.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Calculations of the moon's heat history at different concentrations of radioactive elements taking account of the material differentiation with melting

A mathematical procedure for analyzing the heat conductivity of the lunar surface is discussed. The solution is based on homogeneous and laminated moon models and considers the effects of radioactive elements conveyed to the lunar surface by melting. The various parameters which introduce uncertainties into the numerical analysis are identified. The application of data obtained from radio astronomy and from analyses of lunar samples returned by the Apollo flights is explained. Tables of data are included to show the types and amounts of radioactive materials which have been identified.

Arnatskaya, O. I.↗

Coronal plasma parameters in a long-duration X-ray event observed by Skylab

Physical conditions in the solar coronal X-ray-emitting plasma are analyzed for an event observed by the S-056 X-ray telescope on Skylab in August 1973. The event was located at the west limb and consisted of a series of high looplike structures that underwent dramatic changes in structure and intensity for more than 50 hours. A total X-ray energy (8 to 20 A) of about 5 by 10 to the 28th power erg is obtained by integrating over the lifetime and is shown to be equivalent to that reported previously for a 'typical' flare. It is found that a first major peak in total X-ray flux and thermal energy density was due to an increase in temperature, a major secondary peak in these parameters corresponded to the time of maximum in density, and conductivity energy losses near maximum intensity were much greater than radiative losses. Models dealing with sources of heating and energy balance are evaluated.

Vorpahl, J. A.↗

X-ray emission from the winds of hot stars

A phenomenological theory is proposed for the structure of the unstable line-driven winds of early-type stars. These winds are conjectured to break up into a population of blobs that are being radiatively driven through, and confined by ram pressure of an ambient gas that is not itself being radiatively driven. Radiation from the bow shocks preceding the blobs can account for the X-ray luminosity of zeta Puppis. The theory breaks down when used to model the much lower density wind of tau Scorpii, for then the blobs are destroyed by heat conduction from shocked gas. This effect explains why the profiles of this star's UV resonance lines depart from classical P Cygni form.

Lucy, L. B.↗

Preliminary results of a study of the relationship between free stream turbulence and stagnation region heat transfer

The mechanism that causes free stream turbulence to increase heat transfer in the stagnation region of turbine vanes and blades was studied. The work is being conducted in a wind tunnel at atmospheric conditions to facilitate measurements of turbulence and heat transfer. The model size is scaled up to simulate Reynolds numbers (based on leading edge diameter) that are to be expected on a turbine blade leading edge. Reynolds numbers from 13,000 to 177,000 were run in the present tests. Spanwise averaged heat transfer measurement with high and low turbulence were made with rough and smooth surface stagnation regions. Results of these measurements show that the boundary layer remains laminar in character even in the presence of free stream turbulence at the Reynolds numbers tested. If roughness is added the boundary layer becomes transitional as evidenced by the heat transfer increase with increasing distance from the stagnation line. Hot wire measurements near the stagnation region downstream of an array of parallel wires has shown that vorticity in the form of mean velocity gradients is amplified as flow approaches the stagnation region.

Vaanfossen, G. J., Jr.↗

Preliminary results of a study of the relationship between free stream turbulence and stagnation region heat transfer

The mechanism that causes free stream turbulence to increase heat transfer in the stagnation region of turbine vanes and blades was studied. The work is being conducted in a wind tunnel at atmospheric conditions to facilitate measurements of turbulence and heat transfer. The model size is scaled up to simulate Reynolds numbers (based on leading edge diameter) that are to be expected on a turbine blade leading edge. Reynolds numbers from 13,000 to 177,000 were run in the present tests. Spanwise averaged heat transfer measurement with high and low turbulence were made with rough and smooth surface stagnation regions. Results of these measurements show that the boundary layer remains laminar in character even in the presence of free stream turbulence at the Reynolds numbers tested. If roughness is added the boundary layer becomes transitional as evidenced by the heat transfer increase with increasing distance from the stagnation line. Hot wire measurements near the stagnation region downstream of an array of parallel wires have shown that vorticity in the form of mean velocity gradients is amplified as flow approaches the stagnation region.

Vanfossen, G. J., Jr.↗

Laser anemometry measurements of natural circulation flow in a scale model PWR reactor system

The natural circulation of a single phase fluid in a scale model of a pressurized water reactor system during a postulated grade core accident is analyzed. The fluids utilized were water and SF6. The design of the reactor model and the similitude requirements are described. Four LDA tests were conducted: water with 28 kW of heat in the simulated core, with and without the participation of simulated steam generators; water with 28 kW of heat in the simulated core, with the participation of simulated steam generators and with cold upflow of 12 lbm/min from the lower plenum; and SF6 with 0.9 kW of heat in the simulated core and without the participation of the simulated steam generators. For the water tests, the velocity of the water in the center of the core increases with vertical height and continues to increase in the upper plenum. For SF6, it is observed that the velocities are an order of magnitude higher than those of water; however, the velocity patterns are similar.

Kadambi, J. R.↗

The normal modes of the thermosphere

The linearized momentum, energy, and continuity equations for the thermosphere can be reduced to a form that gives the vertical structure for each horizontal wave mode. The vertical structure equation can be described in terms of the normal modes, or eigenmodes, of the thermosphere. The latter are obtained by using a 27-layer model that includes a realistic temperature profile and the effects of the Lorentz force, viscosity, and heat conduction. The normal modes have one real eigenfrequency for every two complex conjugate eigenfrequency values. The real modes have a dominant rotational wind component and are nonpropagating. The complex modes have comparable divergent and rotational wind components. The complex eigenvalues give vertically propagating modes, primarily associated with the transient response to forcing, and are significantly affected by the dissipation in the upper E region and F region. Results show that the rotational wind component dominates in the steady state when the forcing is due to the two-cell convection pattern at high latitudes and that the normal modes explain the large shears and large winds speeds that are typically observed in the high-latitude E region. The vertical energy flux for the normal modes is also calculated. The results show that the flux is upward above 130 km but downward in the lower E region for the total solution. The downward energy flux is a contribution from the real eigenmode structure.

Larsen, M. F.↗

An analysis of the response of the thermospheric normal modes to temporally varying convection at high latitudes

The response of the high-latitude thermosphere to temporal variations in the forcing due to the plasma convection pattern is investigated using a normal modes model that includes realistic temperature and electron density profiles as well as the effects of viscosity, heat conduction, and Lorentz forcing. The calculations show that the rotational wind component is only slightly modified when the plasma rotation is varied, while the divergent wind component responds more strongly to variations of the rotation rate, with the effect being greatest in the E region. A variation of the rotation rate then results in a greater magnitude and phase variation with the height of the divergent wind component in the E region, indicating that vertically propagating gravity waves are of importance in that height range.

Clark, M. A.↗

Test Data of Flow Field of Shuttle SRM Nozzle Joint with Bond Defects, Using Unheated Air

The nozzle-to-case joint on the Shuttle SRM (as redesigned after the Challenger accident) features an adhesive sealant filling and bonding the joint, with a wiper O-ring to prevent the adhesive from reaching and disabling the closure O-ring. Flawless implementation of that joint design would ensure that hot, corrosive propellant combustion gases never reach the closure O-ring. However, understanding the flow field related to bonding defects is prudent. A comprehensive test program was conducted to quantify such flow fields and associated heating environments. A two-dimensional, full-scale model represented 65 inches of the nozzle joint, using unheated air as the test medium, in a blowdown mode. Geometry variations modeled RSRM assembly tolerances, and two types of bonding defects: pullaways and blowholes. A range of the magnitude of each type defect was tested. Also a range of operational parameters was tested, representative of the RSRM flow environment, including duplication of RSRM Mach and Reynolds numbers. Extensive instrumentation was provided to quantify pressures, heat rates, and velocities. The resulting data established that larger geometric defects cause larger pressure and larger heating, at the closure O-ring region. Velocity trends were not so straight-forward. Variations in assembly tolerances did not generally affect flow fields or heating. Operational parameters affected flow fields and heating as might be expected, increasing density or velocity increased heating. Complete details of this test effort are presented.

Hair, Leroy M.↗