SEARCH · Search NASA
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.
Convective and radiative heat transfer simulation, describing use of shock tubes and ballistic ranges and combinations thereof
A computational fluid dynamics (CFD) model was built to simulate planned testing of heater assemblies within a canister and overpack for the Hanford Lead Canister (HLC) project. The HLC is a canister storage system that will contain heaters to simulate the decay heat of nuclear material and provide the canister storage system with environmental conditions equivalent to the operating conditions on a dry storage pad. The HLC will be equipped with long-term data collection and monitoring systems to provide an early warning of corrosion, pitting, cracking, or other signs of canister degradation that might threaten the integrity of the containment boundary over the potentially long term of dry storage. An important part of the HLC development is to make pretest numerical predictions for the behavior of the heated canister during the simulated radiolytic decay heat testing, which simulates the dry storage system during loading operations. The simulated radiolytic decay heat test is planned for mid-2024 in a configuration that includes the heater assembly, overpack, and canister, but with the lids removed to allow loading cesium and strontium capsules into the canister. One of the goals of the test is to evaluate the thermal behavior of the canister and overpack assembly in the ambient air of the test facility, which will provide data critical to validating the thermal models and understanding how the HLC will perform as a system once deployed. To best approximate real-world conditions, the CFD model includes the full air volume of the mock-up truck bay the heated canister test will be performed in, enabling detailed investigation of how the heated canister affects airflow around it. Rigorous pre-deployment testing of the complete HLC cask and canister system is intended to be completed before the HLC is deployed in the 2028 timeframe. This study presents the pre-test temperature predictions of the simulated radiolytic decay heat test. A description of the heater assembly, canister, and overpack system is presented. The model was developed with the commercial CFD software STAR-CCM+. An uncertainty analysis was run with the CFD model to determine the uncertainty in the temperature predictions and provide a range over which the predicted temperatures are expected to vary. The uncertainty analysis was preformed by coupling STAR-CCM+ with the software Dakota, which provides advanced parametric analyses, including quantification of margins and uncertainty with computational models. This work is expected to provide insight into SNF canister behavior.
Microwave heating has many potential lunar applications including sintering regolith for lunar surface stabilization and heating regolith for various oxygen production reactors. The microwave properties of lunar simulants must be understood so this technology can be applied to lunar operations. Dielectric properties at microwave frequencies for a common lunar simulant, JSC-1AC, were measured up to 1100 C, which is approximately the melting point. The experimentally determined dielectric properties included real and imaginary permittivity (epsilon', epsilon"), loss tangent (tan delta), and half-power depth, the di stance at which a material absorbs 50% of incident microwave energy. Measurements at 2.45 GHz revealed tan delta of JSC-1A increases from 0.02 at 25 C to 0.31 at 110 C. The corresponding half-power depth decreases from a peak of 286 mm at 110 C, to 13 mm at 1100 C. These data indicate that JSC-1AC becomes more absorbing, and thus a better microwave heater as temperature increases. A half-power depth maximum at 100-200 C presents a barrier to direct microwave heating at low temperatures. Microwave heating experiments confirm the sluggish heating effect of weak absorption below 200 C, and increasingly strong absorption above 200 C, leading to rapid heating and melting of JSC-1AC.
Given the lack of fusion-relevant component test facilities, current estimates of the thermo-fluid performance of plasma-facing components are based for the most part on numerical simulations. A major source of uncertainty in these simulations is the semiempirical turbulence (closure) models for the Reynolds stresses appearing in the governing Reynolds-averaged Navier-Stokes equations, which involve a set of constants that depend upon the flow. The objective of this study is to evaluate Bayesian parameter estimation of turbulence closure constants in ANSYS Fluent to model heat transfer in impinging jets. The Bayesian statistical calibration produces a probability distribution for these constants from experimental data; the maximum a posteriori estimates are then taken to be the calibrated constants, or parameters. The turbulence model constants are calibrated using an experimental study of a submerged jet of air impinging on a flat heated surface at Reynolds numbers Re = O(10 4 ) and impingement distance in jet diameters H/d = 2. Numerical predictions using the calibrated model parameters are then compared with those generated using the default constants. Predictions obtained with model parameters calibrated on datasets of two different sizes are compared to evaluate the effect of the number of calibration samples. Lastly, the extrapolative ability of the calibrated model is examined by predictions at a Re beyond the calibration values.
To mitigate the variation in demand on the electric grid, thermal energy storage (TES) is an alternative to electric batteries and constructing new peaking power plants. The aggressive carbon reduction goals of the United States are accelerating the electrification of building equipment. If space heating and cooling in the US were electrified using conventional heat pump technologies, the winter electric grid peak would approximately double. Integrating TES with heat pump (HP) enables electrification of space cooling and heating devices without overtaxing the grid.This paper evaluates the benefits of a novel phase change material (PCM) integrated heat pump configuration via simulation. A one-dimensional PCM heat exchanger model which discretizes the PCM tank and refrigerant tubes into small control volumes is developed. Each control volume can have different PCM temperatures, PCM properties, and heat transfer coefficients. The PCM tank is charged by a wrapped tank condenser and discharged by an internal refrigerant coil. The PCM heat exchanger model is integrated into DOE/ORNL Heat Pump Design Model for heat pump system simulation.To demonstrate the performance of the PCM integrated heat pump, a case study in Chicago was performed. A Time-of-Use utility structure-based control strategy is implemented to schedule the PCM tank charging and discharging mode switching. Compared with a conventional electric heat pump, the PCM integrated heat pump shows superior performance on load shifting and utility cost reduction. As a result, the proposed system demonstrates 24.6% utility saving for cooling application and 25.8% utility saving for heating application.
Testing Dacron parachute components and assemblies by exposure to ethylene oxide sterilization, simulated package loading, and heat cycle
Heat pumps have a significant impact on the climate, and numerous studies exist discussing the environmental impact of heat pumps quantitatively. Life Cycle Climate Performance (LCCP) is a widely accepted metric and can evaluate the carbon footprint of heat pumps from cradle to grave. LCCP calculation typically requires inputs like annual energy consumption, the material used, refrigerant type, and charge level, which are usually the intermediate or final results from heat pump simulation tools. Thus, a platform linking the heat pump simulation software, local weather conditions, and LCCP evaluation tool can simplify the heat pump design process, especially for complicated systems. This paper presents an LCCP evaluation tool accessing the outputs from the DOE Heat Pump Design Model and the building loads and temperature bin data from AHRI standard 210/240. As for case studies, this paper investigates the LCCP values and life cycle electricity costs (LCEC) of a given cooling-based heat pump with R410A in five typical US cities. In conclusion, this cooling-based heat pump has the highest LCCP value at Chicago and lowest LCCP value at Seattle. In contrast, the heat pump has the highest LCEC at Seattle and the lowest LCEC at Miami. As a future work, a more comprehensive carbon emission and cost analysis of heat pump systems will be conducted through all major cities around the world.
Heat-transfer measurements were made on a simulated glide-rocket shape in free flight at Mach numbers up to 10 and free-stream Reynolds numbers of 2 x 10 based on distance along surface from apex and 3 x 10 based on nominal leading-edge diameter. The model simulated the bottom of a 75 deg delta wing at 8O deg angle of attack. The data indicated that for the test conditions a modified three-dimensional stagnation-point theory will predict to reasonable engineering accuracy the heating on a highly swept wing leading edge, the heating being reduced by sweep by the 3/2 power of the cosine of the sweep angle. The data also indicate that laminar heating rates over the windward surface of a highly swept flat glider wing at moderate angles of attack can be predicted with reasonable engineering accuracy by flat-plate theory using wedge local flow conditions and basing Reynolds numbers on lengths from the wing leading edge parallel to the surface center line.
Two-dimensional three-velocity electrostatic particle simulations were used to determine the nonlinear evolution of the distributions of auroral ion beams and thereby to determine quantitatively signatures in the ion distributions produced by the ion-ion instability for a variety of plasma conditions in the auroral zone. The signatures determined from these simulations were compared with observations from DE 1, making it possible to characterize semiquantitatively the heating of the ionospheric ions, and to investigate the causes of variability seen in the observations of Reiff et al. (1988). A comparison of the simulation with observations showed features consistent with heating via the ion-ion instability including perpendicular heating in the supersonic regime and parallel heating in the subsonic regime, and a change in the heating between these regimes as the ratio of the H(+) beam speed to the local sound speed decreases.
Numerical simulations are performed to explore the heat transfer characteristics of falling films over horizontal round tubes with uniform heat flux imposed for a range of Reynolds numbers spanning the droplet, jet, and sheet regimes. Simulations results agree well with available experimental measurements. The study analyzes the local as well as average heat transfer behavior under the different flow modes. The numerical results show that the local Nusselt number (Nu) distribution depends on the flow features in each mode and varies substantially in all directions for the respective mode. In the droplet mode, the Nu value varies significantly as the droplet impinges and the remnant liquid-bridge retracts (peak instantaneous Nu near 6), followed by wave propagation over the tube surface with peak Nu around 0.25. For the jet modes, the local maximum in the heat transfer occurs off-center to the impingement location with magnitudes of peak Nu = 3.1 for the inline jet mode and Nu = 2.7 for the staggered jet mode, while on the rest of the tube surface, it has an inverse relation with the liquid film thickness. Substantial variations in the heat transfer value are also recorded in the middle of the two impinging jets with Nu = 0.95 in the inline jet mode where the neighboring jets do not interact, and Nu = 0.60 in the crest region of the staggered jet mode where the neighboring jets interact with each other. In the sheet mode, the Nu was seen to depend on the thickness of the liquid waves traversing over the tube surface. Lower Nu values were recorded beneath the crest location of the liquid waves, which increases (1.4–11.6% depending on circumferential location) abruptly in magnitude at the advancing fronts of the waves. The temperature distribution in the liquid film in each of the modes was examined to evaluate the mechanism of heat transfer process. This study also compares the local heat transfer coefficient distribution with the analytical heat transfer models derived to predict heat transfer performance over horizontal tube surfaces.
LEAF-Lite (Laser Enhanced Arc-Jet Facility) is a radiative laser heating facility that has been added to the 60 MW Interaction Heating Facility (IHF) convective plasma arc-jet located at NASA Ames Research Center. Together, these two systems can simulate both convective and radiative heating at heat fluxes reaching 551 W/cm2 by simultaneously combining a highest measured heat flux of 160 W/cm2 convective and 391 W/cm2 radiative heating on a 152-mm x 152-mm wedge model configuration. Adding radiant heating to an existing convective facility better simulates Earth atmospheric entry from hyperbolic lunarreturn speeds. The radiative heat is provided by multiple 50-kW CW IR lasers, which is nearly uniform across the illuminated surface with a total variation less than 6%, while the convective heat is provided by a high enthalpy plasma arc-jet. In a later phase, the facility will expand to test panel test articles of 432-mm x 432-mm and provide 100 W/cm2 of radiative heating in a plasma convective flow environment. The paper describes this new combined heating capability, its current testing conditions, and the unique application of the laser system with respect to the Orion test flight lunar orbits.
This dataset contains data associated with the preprint “A laboratory method for simulating the effects of subsequent fires on pyrogenic organic matter at different exposure depths in a sand matrix” (https://doi.org/10.1101/2024.07.30.605814), which has been accepted for publication in the International Journal of Wildland Fire. The dataset was generated to investigate the effects of subsequent fires on preexisting pyrogenic organic matter (PyOM) at different soil depths. Specifically, we designed a laboratory method to simulate soil heating from above, applying realistic heat flux profiles (High, Low, Control) derived from jack pine (Pinus banksiana Lamb.) log burns. Using a cone calorimeter, these profiles were applied to buried jack pine PyOM (produced at 350°C) in a sand matrix pyrocosm to examine PyOM mass loss at different depths (surface, 1 cm, 5 cm). Each treatment (e.g., High heat flux at 1 cm depth) was applied once to five replicates in a pyrocosm.The dataset includes mass loss fraction of PyOM and time-series temperature profiles (measured with type K thermocouples). Additionally, calibration datasets for high and low heat flux profiles are provided, including programmed heat flux profile data and three sets of calibration results. All datasets are in .csv format and were analyzed in R using code found at https://github.com/MengmengLuo/Reburning-pyrogenic-organic-matter-A-laboratory-method-for-dosing-dynamic-heat-fluxes-from-above. This dataset supports research on PyOM transformations during subsequent fires, and methodological advancements in simulating soil heating. It can be used to explore fire-induced PyOM alterations across different fire intensities and soil depths.
The heating of He(+) ions observed by the GEOS satellite at L of about seven starting from an initial anisotropic proton distribution is studied. A simulation of electromagnetic ion cyclotron waves in a plasma consisting of hot anisotropic protons, a dominant thermal hydrogen plasma, and a minority species of helium ions is utilized. It is shown that the thermal ions are heated to 100 times their initial temperature due to the growth of large-amplitude ion cyclotron waves. He(+) ions are heated mainly in the perpendicular direction, and they are heated more than cold H(+) ions. The He(+) ions are heated in a two-step process. First, they are set into oscillations by the growing wave until some of the ions reach a parallel velocity on the order of the resonant velocity; then, strong heating occurs.
The present work provides an improved fluid flow and heat-transfer modeling of float-zone growth by introducing a RF heating model so that an ad hoc heating temperature profile is not necessary. Numerical simulations were carried out to study the high-temperature float-zone growth of titanium carbide single crystal. The numerical results showed that the thermocapillary convection occurring inside the molten zone tends to increase the convexity of the melt-crystal interface and decrease the maximum temperature of the molten zone, while the natural convection tends to reduce the stability of the molten zone by increasing its height. It was found that the increase of induced heating due to the increase of applied RF voltage is reduced by the decrease of zone diameter. Surface tension plays an important role in controlling the amount of induced heating. Finally, a comparison of the computed shape of the free surface with a digital image obtained during a growth run showed adequate agreement.
The feasibility of simulating in a laboratory the heating environment of the base region of the proposed reusable single-stage-to-orbit vehicle during its ascent is examined. The propellant is assumed to consist of hydrocarbon (RP1), liquid hydrogen (LH2), and liquid oxygen (LO2), which produces CO and H2 as the main combustible components of the exhaust effluent. Since afterburning in the recirculating region can dictate the temperature of the base flowfield and ensuing heating phenomena, laboratory simulation focuses on the thermochemistry of the afterburning. By extrapolating the Saturn V flight data, the Damkohler number, in the base region with afterburning for SSTO vehicle, is estimated to be between 30 and 140. It is shown that a flow with a Damkohler number of 1.8 to 25 can be produced in an impulse ground test facility. Even with such a reduced Damkohler number, the experiment can adequately reproduce the main features of the flight environment.
The Nuclear Thermal Rocket Element Environmental Simulator (NTREES) facility at NASA’s Marshall Space Flight Center (MSFC) subjects materials and components typically found on nuclear thermal propulsion (NTP) systems to non-nuclear simulations of prototypic NTP reactor thermohydraulic environments. With 1.2 MW of inductive power available, the currently employed methods to simulate nuclear heating in NTREES consist of 1) inductively heating a test specimen, which typically cannot be performed on ceramic materials, or 2) heating gas that envelops the test specimen. The latter of these has provided a basis for developing so-called “immersion rigs” to be used in NTREES. Furthermore, immersion rigs are especially useful in the case of smaller, technically challenging test specimens. One such device, the “NTREES-5 Hydrogen Heat Exchanger” (N5HHX-01), designed and deployed to test small material coupons to ~0.5 g/s of 2800 K, 500 psi Hydrogen, is discussed herein. To achieve this desired delivery temperature level, the most suitable material for the N5HHX’s construction was pure tungsten, which imposed budgetary and manufacturing challenges. Also discussed in this paper is the design of the N5HHX rig’s sample container, which became the first iteration of improving NTREES small material coupon test operations by permitting rapid changeout of different material samples.
Artificial heating in plasma simulations is a well-known phenomenon which occurs when, among other things, the Debye length is poorly resolved by the simulation mesh. Here, in this work, the degree to which numerical-heating occurs during a simulation of a nanosecond atmospheric pressure streamer discharge is examined. The streamer is simulated using a two-dimensional finite-element, particle-in-cell code Empire, which uses direct simulation Monte Carlo for binary particle interactions. Initially, an estimate of the numerical-heating rate applied to Empire is performed using a simple plasma model. Second, a positive atmospheric pressure streamer discharge simulation is performed to study the effects of numerical heating on plasma density, electron temperature, and streamer velocity. The nominal Debye length is approximately 1 μm and the amount of numerical heating introduced in the simulation is varied by using mesh sizes ranging from 2 μm to 20 μm. A measurable numerical heating quantity is proposed that can be used to estimate the appropriate element size and quantify the numerical-heating that can be expected over the simulation time for an atmospheric pressure streamer. In conclusion while Δx/λ D violations can be an issue it is not likely to be an issue with streamer discharges that are temporally short and occur in environments where collision frequencies are high. This result validates the rationale of grid size choices for a large amount of previously published works where Δx/λ D violation was not clearly addressed. Primary finding of this work is that numerical heating is of minor concern for plasma simulations where electron–neutral collisions are numerous such that multiple collisions can occur within a single plasma period.
Magnetohydrodynamic (MHD) simulations of electrically exploded aluminum and copper rods demonstrate a technique to validate equations of state (EOS) for rapidly Joule-heated conductors. The balance of internal and magnetic forces at the conductor-insulator interface drives the metal there along the vaporization phase boundary. Variations between critical points and vaporization curves in existing models predict differing densities and temperatures in MHD simulations for these models. Here, the inclusion of Maxwell constructs in the liquid-vapor biphase region of the EOS caused the rod surface to vaporize earlier in time than unmodified tables with van der Waals loops. Velocimetry of recent experiments is used to validate the location of the vaporization curve in existing EOS models and differentiate between the vapor dome treatments. Dielectric coatings applied to the metal surface restricted the conductor’s expansion and diverted the metal into the warm dense matter regime.