SEARCH · Engineering Papers
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.
Computer program evaluates heat transfer modes and calculates either the transient or steady-state temperature distributions throughout an object of complex shape when heat sources are applied to specified points on the object. It uses an electrothermal model to simulate the conductance, heat capacity, and temperature potential of the object.
A thermal computer model of the NSTAR (Nasa Solar Electric Propulsion Technology Applications Readiness) xenon ion thruster has been produced using a lumped parameter thermal nodal network scheme. This model contains 104 nodes on the thruster and was implemented using SINDA (Systems Improved Numerical Differencing Analyzer) and TRASYS (Thermal Radiation Analyzer System) on various UNIX workstations. The model includes radiation and conduction heat transfer, the effect of plasma interaction on the thruster, and an account for finely perforated surfaces. The model was developed in conjunction with an NSTAR thruster outfitted with approximately 20 thermocouples for thermal testing at the NASA Lewis Research Center. The results of these experiments were used to calibrate and confirm the computer model first without and then with the plasma interaction. The calibrated model was able to predict discharge chamber temperatures to within 10 C of measured temperatures. To demonstrate the ability of the model under various circumstances the heat flux was examined for a thruster operating in the environment of space.
The storage of solar energy in a solid form, referred to as a “solar fuel”, can be achieved through a process known as endothermic solar thermochemistry. This process transforms the absorbed solar energy into a stable and retrievable form that can be stored for extended periods of time. This paper presents a low–order heat transfer model of a counter–current tubular falling bed reactor designed to produce thermally reduced magnesium manganese oxide pellets for long duration thermochemical energy storage. The energy required for the endothermic reduction was supplied by concentrated solar energy or renewable electricity via indirect heating of the gas and solid reactants flowing in a ceramic tube. The counter-current gas flow enhances the mixing of the solid particles with the heat recuperation zone, allowing the gas and particles to enter and exit the tubular reactor close to room temperature. Further, the reactor was vertically oriented and was heated circumferentially by an adjustable level heat flux along a finite segment of its length. The temperature distribution of the reactor in response to transient changes along the tube was modeled by considering conduction, convection, and radiation heat transfer. Governing equations for the heat transfer model were solved by discretizing the reactor tube into a finite number of control volumes and using an energy balance for the heat exchange between the reactor wall, gas, and particles within the control volume. The energy absorbed during this endothermic reaction was modeled numerically by fitting the data of the chemical conversion rate with the corresponding temperature of particles in the heating zone. The numerical model has been experimentally validated using a reactor prototype made of a 121.92 cm alumina tube heated by a 7kW electric tube–furnace. The alumina tube receives magnesium manganese oxide pellets of 3.66±0.516 mm in diameter from the top, and a counter–current gas flow from the bottom. The reactor wall temperature was monitored by six thermocouples installed along the reactor tube length. The experimental procedure was numerically simulated, and the temperature variation along the reactor tube was compared with a matrix of experimental runs for a range of particles mass flowrates (0.75–1.25g/s) and corresponding gas flowrates (36–65 SLPM). The reactor system was heated gradually from room temperature to a steady state temperature of 1673K, and then cooled down to room temperature. The heating and cooling processes were simulated, and the numerical and experimental results were compared throughout processes. The numerical model showed similar trends to the experimental results, with an error of 0.69 to 7.9% for the particle inlet and 0.7 to 7.9% for the gas inlet during steady-state operation. The proposed numerical model can be implemented as a simplified physical model to design a feedback control system to regulate reactor temperature.
Heat exchangers are critical components of aerospace systems and improve efficiency of operation by providing waste heat recovery. Supercritical CO has emerged as a promising candidate working fluid due to its potential in increasing the power density of heat exchangers. In this paper, a generalized costing model is developed to estimate the capital costs incurred to manufacture microtube shell-and-tube heat exchangers. This model is utilized in conjunction with an accurate and efficient 2D numerical shell-and-tube heat exchanger performance prediction model to conduct optimization studies with two key objectives – minimization of cost and maximization of heat exchanger power density – on supercritical CO microtube heat exchangers utilizing superalloy Haynes 282 as the solid material. A methodology is then demonstrated to optimize these heat exchangers using Particle Swarm Optimization for aerospace applications and highly compact and cost-effective optimal designs with power density around 20 and cost per conductance less than 5 $ ∙ K/W are obtained.
The solution of the inverse problem of nonsteady heat conduction is discussed, based on finding the coefficient of the heat conduction and the coefficient of specific volumetric heat capacity. These findings are included in the equation used for the electrical model of this phenomenon.
This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.
A primary concern in modern fluid dynamics research is the experimental verification of computational aerothermodynamic codes. This research requires high precision and detail in the test model employed. Ceramic materials are used for these models because of their low heat conductivity and their survivability at high temperatures. To fabricate such models, slip casting techniques were developed to provide net-form, precision casting capability for high-purity ceramic materials in aqueous solutions. In previous slip casting techniques, block, or flask molds made of plaster-of-paris were used to draw liquid from the slip material. Upon setting, parts were removed from the flask mold and cured in a kiln at high temperatures. Casting detail was usually limited with this technique -- detailed parts were frequently damaged upon separation from the flask mold, as the molded parts are extremely delicate in the uncured state, and the flask mold is inflexible. Ceramic surfaces were also marred by 'parting lines' caused by mold separation. This adversely affected the aerodynamic surface quality of the model as well. (Parting lines are invariably necessary on or near the leading edges of wings, nosetips, and fins for mold separation. These areas are also critical for flow boundary layer control.) Parting agents used in the casting process also affected surface quality. These agents eventually soaked into the mold, the model, or flaked off when releasing the case model. Different materials were tried, such as oils, paraffin, and even an algae. The algae released best, but some of it remained on the model and imparted an uneven texture and discoloration on the model surface when cured. According to the present invention, a wax pattern for a shell mold is provided, and an aqueous mixture of a calcium sulfate-bonded investment material is applied as a coating to the wax pattern. The coated wax pattern is then dried, followed by curing to vaporize the wax pattern and leave a shell mold of the calcium sulfate-bonded investment material. The shell mold is cooled to room temperature, and a ceramic slip is poured therein. After a ceramic shell of desired thickness has set up in the shell mold, excess ceramic slip is poured out. While still wet, the shell mold is peeled from the ceramic shell to expose any delicate or detailed parts, after which the ceramic shell is cured to provide a complete, detailed, precision ceramic article without parting lines.
Temperatures recorded during two KOSI comet nucleus simulation experiments strongly suggest that heat transport by vapor flow into the interior of the sample is very important. Two comet nucleus simulation experiments have been done by the KOSI team in a big space simulator. The thermal evolution of the sample during insolation and the results of simplified thermal evolution calculations are discussed. The observed thermal histories cannot be explained by a simple model with heat transferred by heat conduction at a constant conductivity, so a coupled heat and mass transfer problem was considered. The porous ice matrix was assumed to have a constant thermal conductivity and to be in thermal equilibrium with vapor in the pores, the internal pressure being the vapor pressure. The vapor was modelled as an ideal gas because, at the temperatures relevant to the problem, the mean free path length of the vapor molecules is large in comparison with the pore dimensions. The heat capacity at constant volume per unit mass of the two phase mixture was also assumed constant. The vapor was allowed to flow and transfer heat in response to an internal pressure gradient.
We present simulations of heat flow relevant to gaspipe experiments on the National Ignition Facility to investigate kinetic effects on transport phenomena. D 2 and neopentane (C 5 H 12 ) filled targets are used to study the laser preheat stage of a MagLIF scheme where an axial magnetic field is sometimes applied to the target. Simulations were done with the radiation-MHD code HYDRA with a collision-dominated fluid model and the SNB nonlocal electron thermal conduction model. Using the SNB model to evolve the electron temperature increased the heat front propagation of neopentane gas targets compared to a local model by limiting radial heat flow. This increases electron temperature near the axis, which decreases laser absorption. We find that the effect of heat flow models on temperature profiles and laser propagation is modest. Beyond the SNB model, we utilize HYDRA to initialize plasma conditions for the Vlasov–Fokker–Planck K2 code. We run K2 until a quasi-steady state is reached and examine the impact of kinetic effects on heat transport. Although axial heat flow is well predicted by fluid models, the fluid model consistently overpredicts radial heat flow up to 150% in regions with the largest temperature gradient of D 2 filled gaspipes. On the other hand, the SNB nonlocal electron conduction model is found to be adequate for capturing kinetic heat flow in gaspipes.
The paper reports on results of heat-transfer tests conducted on a 1/29-scale model of the X-24C-12I hypersonic research aircraft configuration in a Mach 6 tunnel at a Reynolds number of thirteen million using the phase-change heat transfer technique. Sequences of phase-change heat transfer pattern photographs are presented showing windward side and leeward side heating processes. Theoretical predictions of dimensionless heat transfer coefficients along a data line on lower fuselage and on fuselage side bracket the experimental values. A turbulent heating theory gives good agreement with data when shifted to a new virtual origin.
Using the De Vries models for thermal conductivity and heat capacity, thermal inertia was determined as a function of soil moisture for 12 classes of soil types ranging from sand to clay. A coupled heat and moisture balance model was used to describe the thermal behavior of the top soil, while microwave remote sensing was used to estimate the soil moisture content of the same top soil. Soil hydraulic parameters are found to be very highly correlated with the combination of soil moisture content and thermal inertia at the same moisture content. Therefore, a remotely sensed estimate of the thermal behavior of the soil from diurnal soil temperature observations and an independent remotely sensed estimate of soil moisture content gives the possibility of estimating soil hydraulic properties by remote sensing.
Creep strength-enhanced ferritic (CSEF) steels such as Grade 91 (9Cr-1Mo-V) and Grade 92 (Fe-9Cr-2W-0.5Mo) steels are widely used in the fossil-fuel-fired and nuclear power plants. The weld integrity of these steels is crucial for power plants' safe and reliable operations. Due to harsh service conditions, the steel weld can become susceptible to environmental degradation. Field welding repair is used to restore the degraded weld’s performance where a controlled temper-bead welding technique is commonly used to temper the freshly formed martensite during welding. However, knowledge of weld repairability is limited and experimental trial and error optimization to achieve desired microstructure and joint properties is expensive and time-consuming. Many existing computational models, e.g., finite element models, are limited to solving heat conduction equation and ignoring convective heat transfer due to molten metal flow. These models can result in over-prediction of peak temperatures of weld pool and heat-affected zone (HAZ), which in turn can affect the accuracy of tempering prediction. Moreover, these finite element models require an input of the deposit profiles in advance and thus limits the usability of these models. Here, a molten pool-based, multi-pass multi-layer model has been developed based on computational fluid dynamics (CFD) approach with the Volume of Fluid (VOF) method. The model calculates the bead formation, thereby eliminating the need for pre-determined bead profiles required by finite element models. For computational efficiency, a coordinate system attached to the moving heat source is utilized. A subroutine is developed to convert the temperature profiles in the reference frame stationary to the heat source to that stationary to the workpiece. The converted thermal cycles are then imported into a microstructure model to compute the tempering kinetics and resultant hardness using a Johnson-Mehl-Avrami-Kolmogorov (JMAK), and modified Grange-Baughman parameter. The modeling approach is first developed and validated on single- and multi-pass deposition of stainless steel filler metal onto a SA-533 high strength steel substrate. The models are then applied to a multi-pass V-groove repair weld of Grade 91 steel plate as well as directed energy deposition of Grade 92 steel. Non-destructive characterization of microstructures was performed on Grade 91 and 92 steel welds. Two welding processes, cold metal transfer (CMT) and flux-cored arc welding (FCAW), were investigated for the Grade 91 steel weld samples. For the Grade 92 weld samples, three different heat inputs (low, medium, and high) of gas tungsten arc welding (GTAW) were utilized to replicate traditional field welding processes. The non-destructive evaluation (NDE) method used for this research was immersion ultrasonic testing (UT) using a micro-resolution ultrasonic imaging methodology specifically designed to operate in the through-transmission configuration operating at 20 MHz of frequency. The system used a focused ultrasonic beam spot size diameter between 250-300 μm, and a 6 μm laser vibrometer spot size for detection, to produce highly defined images with longitudinal and mode-converted shear waves. From the micro-resolution ultrasonic C-scan images, three microstructural regions, i.e., weld metal (WM), HAZ, and base metal (BM), were clearly identifiable. Various levels of ultrasonic amplitudes distributed over the three regions were correlated with electron beam backscattered diffraction (EBSD) images using grain size, grain boundaries, and dislocation densities. The results showed that areas with relatively higher ultrasonic amplitude levels were associated with smaller grains and higher dislocation densities, while areas with lower amplitude levels were associated with larger grains and lower dislocation densities. In addition, ultrasonic velocity data obtained across the three different weld microstructural regions of Grade 91 test samples were correlated with optical metallographic images and hardness measurements. The results showed distinctive decreases in ultrasonic velocity and hardness over the HAZ region, where weld failures often occur during service.
A new, thermal-flow simulation code, called SFLOW. has been developed to model the gas dynamics, heat transfer, as well as O-ring and flow path erosion inside the space shuttle solid rocket motor joints by combining SINDA/Glo, a commercial thermal analyzer. and SHARPO, a general-purpose CFD code developed at Thiokol Propulsion. SHARP was modified so that friction, heat transfer, mass addition, as well as minor losses in one-dimensional flow can be taken into account. The pressure, temperature and velocity of the combustion gas in the leak paths are calculated in SHARP by solving the time-dependent Navier-Stokes equations while the heat conduction in the solid is modeled by SINDA/G. The two codes are coupled by the heat flux at the solid-gas interface. A few test cases are presented and the results from SFLOW agree very well with the exact solutions or experimental data. These cases include Fanno flow where friction is important, Rayleigh flow where heat transfer between gas and solid is important, flow with mass addition due to the erosion of the solid wall, a transient volume venting process, as well as some transient one-dimensional flows with analytical solutions. In addition, SFLOW is applied to model the RSRM nozzle joint 4 subscale hot-flow tests and the predicted pressures, temperatures (both gas and solid), and O-ring erosions agree well with the experimental data. It was also found that the heat transfer between gas and solid has a major effect on the pressures and temperatures of the fill bottles in the RSRM nozzle joint 4 configuration No. 8 test.
Hot spots are local regions of high temperature that are widely considered to govern explosive initiation. Hot spot dynamics rests on a delicate balance between heat generation due to chemical reactions and heat loss through thermal conduction, making accurate determinations of the conductivity under extreme conditions a key component of predictive explosive models. We develop here an approach to directly determine the thermal transport properties of explosive hot spots with realistic initial structures through a combination of molecular dynamics (MD) and diffusive heat equation (HEq) modeling. Effective thermal conductivity values are determined by fitting HEq models to MD predictions of long timescale hot spot relaxation. The approach is applied to model hot spots in the molecular crystalline explosive 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) for a range of shock strengths and two limiting cases for impact orientation. Isotropic and anisotropic HEq models yield similar results, despite TATB exhibiting some of the largest and most anisotropic thermal conductivity values for explosive near normal conditions. The conductivity is found to be a strong function of density, which parametrically captures dependence on temperature, pressure, and material state. The associated root-mean-square errors of the fitted HEq models are approximately 5% of MD predicted final equilibrium temperatures. Here, the conductivity values determined here for TATB hot spots are considerably larger than those used in a prior hot spot criticality study, which may significantly impact predictions for critical hot spot sizes. The approach provides a convenient foundation for determining the effective thermal conductivity for hot spot problems in other explosives and directly yields information on reasonable approximations that might be taken in higher-level models for those materials.
This report provides an initial study for producing reduced-order models (ROMs) of pebble-bed high temperature gas reactor (HTGR) models for the purposes of design optimization. As an initial study, this work is meant to be exploratory---identifying useful workflows and methods for ROM generation---and not meant to be a catch-all analysis of HTGR ROM generation and usage for optimization. This report summarizes three tasks performed in Fiscal Year 2022: 1) the creation of HTGR model, 2) the sensitivity analysis of model design parameters, and 3) an introduction to ROM generation techniques. The representative HTGR model created in this work is a multiphysics equilibrium-core using the BlueCRAB (comprehensive reactor analysis bundle) reactor analysis application, coupling four physical phenomena: neutronics, streamline depletion, porous flow thermal hydraulics, and pebble heat conduction. Part of the model creation was identifying some design parameters and quantities of interest that are relevant in an optimization analysis and adjustable in the model. The sensitivity analysis utilized a polynomial chaos expansion methodology to compute global sensitivity metrics. This analysis showed that thermal hydraulics parameters and quantities of interest had a relatively small impact on simulation results. Finally, the ROM generation work involved exploring three different ROM methodologies: polynomial regression, a Gaussian process, and artificial neural networks. Using a cross-validation technique to characterize ROM performance, the Gaussian process and single-layer artificial neural networks showed the most promising results. Overall, this study was insightful and the lessons learned will be invaluable for the eventual development of an HTGR design optimization workflow.
A coupled heat and moisture balance model was used to determine the thermal inertia of a grass covered top soil under different meteorological conditions. Relations between thermal inertia and soil moisture were established using the De Vries models for thermal conductivity and heat capacity to relate soil moisture and thermal inertia as a function of soil type. A sensitivity study of the surface roughness length and thermal inertia on diurnal surface temperature shows the necessity of focusing on the night time surface temperature rather than on the day time surface temperature, in order to estimate the soil moisture content of the top soil.
Explore the source record for details and available documents.
During the Arctic winter, the conductive heat flux through the sea ice and snow balances the radiative and turbulent heat fluxes at the surface. Snow on sea ice is a thermal insulator that reduces the magnitude of the conductive flux. The thermal conductivity of snow, that is, how readily energy is conducted, is known to vary significantly in time and space from observations, but most forecast and climate models use a constant value. This work begins with a demonstration of the importance of snow thermal conductivity in a regional coupled forecast model. Varying snow thermal conductivity impacts the magnitudes of all surface fluxes, not just conduction, and their responses to atmospheric forcing. Given the importance of snow thermal conductivity in models, we use observations from sea ice mass balance buoys installed during the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition to derive the profiles of thermal conductivity, density, and conductive flux. From 13 sites, median snow thermal conductivity ranges from 0.33 W m −1 K −1 to 0.47 W m −1 K −1 with a median from all data of 0.39 W m −1 K −1 from October to February. In terms of surface energy budget closure, estimated conductive fluxes are generally smaller than the net atmospheric flux by as much as 20 W m −2 , but the average residual during winter is −6 W m −2 , which is within the uncertainties. The spatial variability of conductive heat flux is highest during clear and cold time periods. Higher surface temperature, which often occurs during cloudy conditions, and thicker snowpacks reduce temporal and spatial variability. These relationships are compared between observations and the coupled forecast model, emphasizing both the importance and challenge of describing thermodynamic parameters of snow cover for modeling the Arctic as a coupled system.