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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

One-Dimensional Hybrid Satellite Track Model for the Dynamics Explorer 2 (DE 2) Satellite

A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. satellite track code but also includes elements of Roble et al. global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the MSIS-90 model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule beating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Deng, Wei↗
NASA NTRS2002

Macroscopic Relationships among Latent Heating, Precipitation, Organized Convection and the Environment

Three-dimensional Cloud Resolving Model (CRM) simulations were conducted to examine the squall line observed on 26 January, 1999 from the Tropical Rainfall Measuring Mission Large Scale Biosphere Atmosphere Experiment in Amazonia (TRMM-LBA) field campaign. The computational domain was 600 kilometers x 180 kilometers x 20 kilometers with a horizontal resolution of 1 kilometer and a vertical resolution of 200 meters. The CRM was initialized from the Abracos Hill and Rebio soundings. Convection was initiated by a surface-based and NW-SE oriented cold pool over a region 60 kilometers in the y-direction and 30 kilometers wide in the x-direction. The cold pool temperature perturbation is a maximum of -6K at the surface, decreasing linearly to zero at 3 kilometers. The simulated convection is in the form of a NW-SE band that moves toward the southwest at a speed of 8 meters per second, and is generally comparable to radar observations.

Moncrieff, Mitchell W.↗
NASA NTRS2007

Transient Numerical Modeling of Catalytic Channels

This paper presents a transient model of catalytic combustion suitable for isolated channels and monolith reactors. The model is a lumped two-phase (gas and solid) model where the gas phase is quasi-steady relative to the transient solid. Axial diffusion is neglected in the gas phase; lateral diffusion, however, is accounted for using transfer coefficients. The solid phase includes axial heat conduction and external heat loss due to convection and radiation. The combustion process utilizes detailed gas and surface reaction models. The gas-phase model becomes a system of stiff ordinary differential equations while the solid phase reduces, after discretization, into a system of stiff ordinary differential-algebraic equations. The time evolution of the system came from alternating integrations of the quasi-steady gas and transient solid. This work outlines the numerical model and presents some sensitivity studies on important parameters including internal transfer coefficients, catalytic surface site density, and external heat-loss (if applicable). The model is compared to two experiments using CO fuel: (1) steady-state conversion through an isothermal platinum (Pt) tube and (2) transient propagation of a catalytic reaction inside a small Pt tube. The model requires internal mass-transfer resistance to match the experiments at lower residence times. Under mass-transport limited conditions, the model reasonably predicted exit conversion using global mass-transfer coefficients. Near light-off, the model results did not match the experiment precisely even after adjustment of mass-transfer coefficients. Agreement improved for the first case after adjusting the surface kinetics such that the net rate of CO adsorption increased compared to O2. The CO / O2 surface mechanism came from a sub-set of reactions in a popular CH4 / O2 mechanism. For the second case, predictions improved for lean conditions with increased external heat loss or adjustment of the kinetics as in the first case. Finally, the results show that different initial surface-species distribution leads to different steady-states under certain conditions. These results demonstrate the utility of a lumped two-phase model of a transient catalytic combustor with detailed chemistry.

Struk, Peter M.↗
NASA NTRS2006

Code for Multiblock CFD and Heat-Transfer Computations

The NASA Glenn Research Center General Multi-Block Navier-Stokes Convective Heat Transfer Code, Glenn-HT, has been used extensively to predict heat transfer and fluid flow for a variety of steady gas turbine engine problems. Recently, the Glenn-HT code has been completely rewritten in Fortran 90/95, a more object-oriented language that allows programmers to create code that is more modular and makes more efficient use of data structures. The new implementation takes full advantage of the capabilities of the Fortran 90/95 programming language. As a result, the Glenn-HT code now provides dynamic memory allocation, modular design, and unsteady flow capability. This allows for the heat-transfer analysis of a full turbine stage. The code has been demonstrated for an unsteady inflow condition, and gridding efforts have been initiated for a full turbine stage unsteady calculation. This analysis will be the first to simultaneously include the effects of rotation, blade interaction, film cooling, and tip clearance with recessed tip on turbine heat transfer and cooling performance. Future plans call for the application of the new Glenn-HT code to a range of gas turbine engine problems of current interest to the heat-transfer community. The new unsteady flow capability will allow researchers to predict the effect of unsteady flow phenomena upon the convective heat transfer of turbine blades and vanes. Work will also continue on the development of conjugate heat-transfer capability in the code, where simultaneous solution of convective and conductive heat-transfer domains is accomplished. Finally, advanced turbulence and fluid flow models and automatic gridding techniques are being developed that will be applied to the Glenn-HT code and solution process.

Fabian, John C.↗
NASA NTRS2018

Interactions Between Asian Air Pollution and Monsoon System: South Asia (ROSES-2014 ACMAP)

Asia's rapid economic growth over the past several decades has brought a remarkable increase in air pollution levels in that region. High concentrations of aerosols (also known as particulate matter or PM) from pollution sources pose major health hazards to half of the world population in Asia including South Asia. How do pollution and dust aerosols regulate the monsoon circulation and rainfall via scattering and absorbing solar radiation, changing the atmospheric heating rates, and modifying the cloud properties? We conducted a series of regional model experiments with NASA-Unified Weather Research and Forecast (NUWRF) regional model with coupled aerosol-chemistry-radiation-microphysics processes over South Asia for winter, pre-monsoon, and monsoon seasons to address this question. This study investigates the worsening air quality problem in South Asia by focusing on the interactions between pollution and South Asian monsoon, not merely focusing on the increase of pollutant emissions.

regional model experiments with NASA-Unified Weath↗
DOE OSTI2025

RELAP5-3D Capabilities for High-Temperature Gas-Cooled Reactor Analysis

RELAP5-3D is the premier systems thermal-hydraulic code in the world. The code has been used significantly for high-temperature gas-cooled reactor (HTGR) analyses over the last 10 years, but there has not been a comprehensive summary of RELAP5-3D capabilities and methods for HTGR modeling. In this paper, we provide a guide for how novice users can build models for prismatic HTGRs. This guide includes best practices for developing heat structures and setting up radiation and conduction heat transfer. We also provide a summary of relevant validation studies. These studies are divided into three categories: separate effects tests, integral effects tests, and nuclear reactor studies. We present a comprehensive review of RELAP5-3D HTGR applications for reactor modeling over the last 10 years. These applications include work for the Modular High-Temperature Gas-Cooled Reactor benchmark, design activities for a prismatic HTGR as part of the U.S. Department of Energy’s Advanced Demonstration and Test Reactor Options Study, and studies performed to design gas-cooled microreactors. Following the comprehensive literature review, we present a gap analysis highlighting capabilities of RELAP5-3D that are available and missing for a complete analysis of HTGR thermal hydraulics. We also discuss the limitations of existing code capabilities and potential improvements that could improve analyses. Current research efforts have yet to fully validate RELAP5-3D for a prismatic HTGR analysis; however, some of this is attributable to confounding factors in the available experimental data as opposed to code capabilities. In conclusion, this paper presents the most comprehensive review to date on RELAP5-3D capabilities, limitations, and applications for HTGR analysis.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Development of High-Temperature Emissivity Measurement Capabilities for Molten Salts at PNNL

In this study, an emissivity measurement system for a liquid molten salt was developed on the bases on the Stefan Boltzmann Law. A prototype device operating in ambient atmosphere was constructed an initial result compared well with available literature. The diffuse emissivity of NaCl-KCl system (NaCl, 51mol% NaCl – 49mol% KCl and KCl) was measured from 550 to 850 °C at a wavelength of 1.3 µm. The system appears to be relatively invariant with composition but a clear trend of decreasing emissivity values from ~0.9 at 550 °C to 0.75 at 850 °C was observed. A schematic for an enclosed system for inert atmosphere and actinide salts was drafted. Several of the components for fabrication of this device, improved laser system, furnace, nickel for capsule construction, have been purchased and are ready for assembly. The new laser system will have an improved measurement range of 275-1000 °C. The clam shell furnace has a maximum temperature of 1200 °C and will support the hermetically sealed inert encapsulation vessels. Nickel plate and tubing will create a relatively corrosion resistance capsule for testing that can withstand temperature up to 1000 °C. Emissivity measurements are important for a larger effort to understand radiative heat transfer within molten salts at elevated temperatures (600-1000 °C). Literature review of optical salt properties and radiative heat transfer models indicated that systems within elevated temperature regimes could see upwards of 20% or more of their total heat transferred through photonic mechanisms. Ultimately, future work will measure the total thermal conductivity and break down the molecular and radiative portions. Thermal conductivity, along with other thermal properties of molten salts, both coolant and fuel, are fundamental to modeling and simulations that support MSR design. The results presented here and ongoing work at PNNL will be useful in selecting optimal salt compositions, designing reactor systems, and optimizing heat transfer within a MSR power plant.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

The right conditions for high-precision dynamic temperature and heat capacity measurement via pyrometry and conductivity

The pursuit of accurate bulk temperature T under extreme conditions has been a long-standing goal of the high pressure science community, complicated by a lack of data to inform models. To reach these extremely high-pressure, high-temperature (high P − T) conditions, a combination of dynamic and heated static experiments (e.g., diamond or gem anvil cel experiments) are used. For example, in a diamond anvil cell (DAC) experiment, a sample placed in the DAC is first pressurized. Following pressurization, the sample T is increased either by heating the entire DAC (usually using resistive heating, and limited to ∼1000K) or by applying intense laser power to the sample surfaces. In a dynamic experiment, the process of pressurizing the sample also heats it. In the case of shock physics experiments, such heating is substantial, easily reaching thousands of Kelvin; in our work we have seen T ∼17000K. Most methods of measuring temperature at ambient are not compatible with experiments under these high-pressure, high-temperature conditions: thermocouples break, melt, or have conductivity properties that differ from ambient where they are calibrated; thermometers would melt; both are too slow. As a result most methods are based on non-contact techniques such as x-ray diffraction broadening, neutron scattering, or optical methods. Of these, optical methods using the visible and near-infrared region of the spectrum are the most commonly used as the sources and detectors are readily available. In the case of optical methods the optical depth, and therefore the measurement location, is limited to the surface. When a window or anvil material is used, heat flows from the sample into the window/anvil. Likewise, if the sample undergoes a change in thermodynamic state, such as expansion upon release, different T may be expected. As a result, the surface or apparent temperature T app measurement will differ from the bulk or interior temperature that is desired. This surface measurement must be related to the bulk measurement using thermal transport models and material models. While it is tempting to conclude that one should just use x-ray methods that directly probe the interior, even these methods have been shown to depend on thermal transport and material models. Regardless of the method used to create the high P − T condition, therefore, we must understand the role of thermal transport and material models upon our interpretation of the T measurement, as well as the errors and uncertainties associated with the choice of models used in the analysis. This is a substantial area of research and this paper is by no means a complete survey of the relevant sources of uncertainty. For example, we have yet to begin to address alternate transport models in a detailed manner (e.g., Tan-Ahrens), or the many models that use additional layers to approximate melting, turbulence, or epitaxial phenomena). Likewise, we have not explored the impact upon uncertainty of thermal models that use temperature-dependent thermal transport coefficients, or the wide range of material models that can be applied. Instead, this paper focuses on using one simple model, the Urtiew-Grover model, to understand the sources of error in T measurement so that we may identify how best to focus future research efforts to return the best improvements and avoid working on over-optimizing a single type of measurement. To this end, we work through some of the best and worst case scenarios for T measurement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Investigation of CAD-based Geometry Workflows for Multiphysics Fusion Problems Using OpenMC and MOOSE

Fusion system designs are complex and require intricate and accurate meshes to be properly modeled. In this study, we investigate the use of CAD-based geometry workflows in fusion systems multiphysics problems. A simplified tokamak was introduced and modeled in CAD using a multiphysics coupling of OpenMC Monte Carlo transport and MOOSE heat conduction. The meshed geometry was prepared using direct accelerated geometry Monte Carlo (DAGMC) for particle transport, and a volumetric mesh was also prepared to be used in MOOSE and to tally OpenMC results. Cardinal was used to run OpenMC Monte Carlo particle transport within MOOSE framework. The heat source distribution and tritium production were calculated in OpenMC. The data transfer system was used to transfer heat source and temperature distribution between OpenMC and MOOSE. Two computational studies related to mesh refinement were performed: (1) refining the DAGMC and volumetric meshes used for tallying results and solving heat conduction and (2) only refining the DAGMC particle transport mesh. The refinement of the tally mesh has a much larger effect on the runtime compared to the refinement of the DAGMC particle transport surface mesh.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

Evaluation of scramjet nozzle configurations and film cooling for reduction of wall heating

Experiments on relaminarization, or reverse transition of turbulent flow, have been conducted with a scramjet combustor and nozzle model. The Mach-2.7 direct-connect scramjet combustor employed swept-ramp fuel injectors, and the entering flow simulated Mach-6 to Mach-7 flight conditions. Tests were also conducted without combustor fuel to produce a higher nozzle entrance Mach number and lower pressures. Additional tests were conducted to evaluate the effectiveness of film cooling on the cowl region. At each test condition studied, the entrance radius had no significant effect on the measured downstream wall heating rates. Analysis of the data when fuel was injected into the combustor shows good agreement with turbulent theory; for the tests without fuel injection and combustion, the data downstream of the nozzle entrance radius are intermediate to the calculated turbulent and laminar solutions. Results of film-cooling tests, conducted with both hydrogen and nitrogen injectants, showed that, per unit mass, the hydrogen was approximately three to four times more effective than nitrogen in reducing the downstream wall heating rates.

Baker, N. R.↗
NASA NTRS1976

Results from a convective heat transfer rate distribution test on a 0.0175 scale model (22-0) of the Rockwell International vehicle 4 space shuttle configuration in the AEDC-VKF tunnel B (OH49B), volume 1

The tests were conducted in a hypersonic wind tunnel at Mach number 8 to investigate reentry mode convective heat--transfer rates to the vehicle 4 shuttle orbiter. The thin skin thermocouple technique was used to obtain the heat transfer rate measurements. A complete set of tabulated data is presented.

Herrera, B. J.↗
NASA NTRS1976

Aeroheating (pressure) characteristics on a 0.10-scale version of the vehicle 3 space shuttle configuration (26-OTS) in the Langley Research Center 4-foot wind tunnel (IH4)

Results of wind tunnel tests, conducted at the Langley Research Center Unitary Plan Wind Tunnel, are presented. The model tested was an 0.010-scale version of the Vehicle 3 Space Shuttle Configuration. Pressure measurements were made on the launch configuration, Orbiter alone, external tank alone, and solid rocket booster alone, to provide heat transfer pressure data. The tests were conducted for a Mach number range from 2.36 to 4.6 and Reynolds number range from 1.2 to 5 million per foot. The model was tested at angles of attack from -10 to 20 deg for a sideslip angle range from -5 to +5 deg, and at sideslip angles from -5 to 48 deg for 0 deg angle of attack. Tabulated data are given and photographs of the test configuration are shown.

Kingsland, R. B.↗
NASA NTRS1974

Pressure and heat flux results from the space shuttle/external fuel tank interaction test at Mach numbers 16 and 19

Heat transfer rates and pressures were measured on a 0.0175-scale model of the space shuttle external tank (ET), model MCR0200. Tests were conducted with the ET model separately and while mated with a 0.0175-scale model of the orbiter, model 21-OT (Grumman). The tests were conducted in the AEDC-VKF Hypervelocity Wind Tunnel (F) at Mach numbers 16 and 19. The primary data consisted of the interaction heating rates experienced by the ET while mated with the orbiter in the flight configuration. Data were taken for a range of Reynolds numbers from 50,000 to 65,000 under laminar flow conditions.

Brewer, E. B.↗
NASA NTRS2000

Numerical Modeling and Test Data Comparison of Propulsion Test Article Helium Pressurization System

A transient model of the propulsion test article (PTA) helium pressurization system was developed using the generalized fluid system simulation program (GFSSP). The model included pressurization lines from the facility interface to the engine purge interface and liquid oxygen (lox) and rocket propellant-1 (RP-1) tanks, the propellant tanks themselves including ullage space, and propellant feed lines to their respective pump interfaces. GFSSP's capability was extended to model a control valve to maintain ullage pressure within a specified limit and pressurization processes such as heat transfer between ullage gas, propellant, and the tank wall as well as conduction in the tank wall. The purpose of the model is to predict the flow system characteristics in the entire pressurization system during 80 sec of lower feed system priming, 420 sec of fuel and lox pump priming, and 150 sec of engine firing.

Holt, Kimberly↗
NASA NTRS1984

The friable sponge model of a cometary nucleus

The mantle/core model of cometary nuclei, first suggested by Whipple and subsequently developed by Mendis and Brin, is modified and extended. New terms are added to the heat conduction equation for the mantle, which is solved in order to obtain the temperature distribution in the mantle and the gas production rate as a function of mantle thickness and heliocentric distance. These results are then combined with some specific assumptions about the mantle structure (the friable sponge model) in order to make predictions for the variation of gas production rate and mantle thickness as functions of heliocentric distance for different comets. A solution of the time-dependent heat conduction equation is presented in order to check some of the assumptions.

Horanyi, M.↗
DOE OSTI2022

Progress on Demonstration of a MOOSE-Based Coupled Capability for Hot Channel Factors in Fast Reactors

Hot channel factors (HCFs) are computed values that account for the impact on predicted peak fuel, cladding, and coolant temperatures due to uncertainties in the as-built reactor’s material properties and geometry as well as uncertainties due to modeling approximations. Reduction in computed HCF values via reduction or elimination of modeling approximations may translate to significant economic savings if the reactor power can be raised due to the extra temperature margin gained. While limited historical datasets exist for sodium-cooled fast reactors (SFRs), there are no available HCF data for lead-cooled fast reactors (LFRs) outside of work generated previously within NEAMS. The computation of HCFs involves insights from reactor physics, thermal fluids and heat conduction calculations to determine how the peak temperatures respond to various uncertainties in the design. Due to the significant advantages for multi-physics coupling offered by the MOOSE framework, Griffin (MOOSE-based reactor physics code), MOOSE Heat Conduction Module, and Cardinal (MOOSE-wrapped multi-physics application which includes the NekRS thermal fluids code) are being coupled together using the MOOSE MultiApp System to develop a highfidelity multi-physics modeling capability for HCF simulations. This high-fidelity coupling workflow may also be beneficial for other fast reactor applications in the future. In previous work, Griffin and NekRS were individually assessed to ensure the necessary capabilities were in place. This work describes initial efforts to couple the codes (including folding in the MOOSE Heat Conduction Module) and determining the workflow for the perturbed calculations which will leverage the Stochastic Tools Module (STM). To our knowledge, this is the first coupling of Griffin and NekRS as well as the first exploratory use of Stochastic Tools Module for Cardinal. In this report, the neutronics code Griffin, the heat conduction solver in MOOSE, and the MOOSE-wrapped application containing NekRS (Cardinal) are linked together to demonstrate the coupled capability. Griffin and Cardinal are linked dynamically by specifying shared libraries. Different coupling hierarchies are tested for selecting the most appropriate coupling strategy. A coupling scheme is selected based on the efficiency of calculation and ease of data communication. Multiple tests are performed to choose suitable mesh structure, model configurations, scheme setup and boundary conditions to avoid loss of energy due to data interpolation between different modules or weak imposition of fluxes in finite element codes. Computational experiments are performed to study the tolerance control of each type of iteration to avoid false convergence. The coupled capability is demonstrated in both single pin and 7-pin models based on LFR materials and geometry. The study finds that the use of too large a time step size in the heat conduction module can lead to temperature oscillation even though the heat conduction equation does not have a time-derivative kernel, but only the time-dependent boundary condition. A 7-pin model without duct region achieved good convergence in the coupled calculation while a 7 pin model with duct region experienced data communication issues which need to be resolved.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1977

The SRB heat shield: Aeroelastic stability during reentry

Wind tunnel tests of a 3% scale model of the aft portion of the SRB equipped with partially scaled heat shields were conducted for the purpose of measuring fluctuating pressure levels in the aft skirt region. During these tests, the heat shields were observed to oscillate violently, the oscillations in some instances causing the heat shields to fail. High speed films taken during the tests reveal a regular pattern of waves in the fabric starting near the flow stagnation point and progressing around both sides of the annulus. The amplitude of the waves was too great, and their pattern too regular, for them to be attributed to the fluctuating pressure levels measured during the tests. The cause of the oscillations observed in the model heat shields, and whether or not similar oscillations will occur in the full scale SRB heat shield during reentry were investigated. Suggestions for modifying the heat shield so as to avoid the oscillations are provided, and recommendations are made for a program of vibration and wind tunnel tests of reduced-scale aeroelastic models of the heat shield.

Ventres, C. S.↗
NASA NTRSDate not supplied

Thermal Modeling of a Compressor for CO2 Removal in the International Space Station

Contaminant removal technology, such as the Carbon Dioxide Removal Assembly (CDRA) on the International Space Station (ISS), is critical to Environmental Control and Life Support Systems (ECLSS), which enable humans to live and work in outer space. The CDRA includes an adsorbent bed, which absorbs CO2 from the cabin air and releases it to a Sabatier reactor for water production when thermally cycled. An effective system would heat and cool the adsorbent quickly and uniformly to maximize the amount of CO2 removed from the cabin and delivered to the Sabatier reactor. Air-Cooled Temperature Swing Adsorption Compressors (AC-TSAC) are promising for use downstream of the CDRA because they are a simple alternative to more mechanically-complex compressors, which are currently in use on the ISS. A thermal model of the AC-TSAC was created in COMSOL and validated by comparing the model outputs to experimental data for the heating phase of the cycle. Subsequently, several design trades and performance sensitivities were conducted, including bed geometry (rectangular vs. cylindrical), bed structural material (i.e., the material that makes up the walls and shelves), thermal conductivity of the adsorbent material and bed structural materials, and the input power for heating. The results from these studies will inform the design of the next generation of the AC-TSAC, and the thermal modeling results will be further validated through testing.

Hannah Alpert↗