Search NASA⌕ Search

SEARCH · Search NASA

Results for “lumped thermal model”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Numerical comparison of convective heat transfer augmentation devices used in cooling channels of hypersonic vehicles

Hypersonic vehicles are exposed to extreme thermal conditions compared to subsonic aircraft; therefore, some level of thermal management is required to protect the materials used. Normally, hypersonic vehicles experience the highest temperatures in the nozzle throat, and aircraft and propulsion system leading edges. Convective heat transfer augmentation techniques can be used in the thermal management system to increase heat transfer of the cooling channels in those areas. The techniques studied in this report are pin-fin, offset-fin, ribbed and straight roughened channel. A smooth straight channel is used as the baseline for comparing the techniques. SINDA '85, a lumped parameter finite difference thermal analyzer, is used to model the channels. Subroutines are added to model the fluid flow assuming steady one dimensional compressible flow with heat addition and friction. Correlations for convective heat transfer and friction are used in conjunction with the fluid flow analysis mentioned. As expected, the pin-fin arrangement has the highest heat transfer coefficient and the largest pressure drop. All the other devices fall in between the pin-fin and smooth straight channel. The selection of the best heat augmentation method depends on the design requirements. A good approach may be a channel using a combination of the techniques. For instance, several rows of pin-fins may be located at the region of highest heat flux, surrounded by some of the other techniques. Thus, the heat transfer coefficient is maximized at the region of highest heat flux while the pressure drop is not excessive.

Maldonado, Jaime J.↗
NASA NTRS2011

Dynamic Simulation of Human Thermoregulation and Heat Transfer for Spaceflight Applications

Models of human thermoregulation and heat transfer date from the early 1970s and have been developed for applications ranging from evaluating thermal comfort in spacecraft and aircraft cabin environments to predicting heat stress during EVAs. Most lumped or compartment models represent the body as an assemblage cylindrical and spherical elements which may be subdivided into layers to describe tissue heterogeneity. Many existing models are of limited usefulness in asymmetric thermal environments, such as may be encountered during an EVA. Conventional whole-body clothing models also limit the ability to describe local surface thermal and evaporation effects in sufficient detail. A further limitation is that models based on a standard man model are not readily scalable to represent large or small subjects. This work describes development of a new human thermal model derived from the 41-node man model. Each segment is divided into four concentric, constant thickness cylinders made up of a central core surrounded by muscle, fat, and skin, respectively. These cylinders are connected by the flow of blood from a central blood pool to each part. The central blood pool is updated at each time step, based on a whole-body energy balance. Results show the model simulates core and surface temperature histories, sweat evaporation and metabolic rates which generally are consistent with controlled exposures of human subjects. Scaling rules are developed to enable simulation of small and large subjects (5th percentile and 95th percentile). Future refinements will include a clothing model that addresses local surface insulation and permeation effects and developing control equations to describe thermoregulatory effects such as may occur with prolonged weightlessness or with aging.

Miller, Thomas R.↗
NASA NTRS1980

Systems improved numerical differencing analyzer

Program solves physical problems governed by diffusion-type equations, provided that equations can be modeled by lumped-parameter representation. Program is used for thermal analysis, and could be adapted to solve Fourier, Poisson, and Laplace differential equations. Program is in FORTRAN IV and Assembler for execution on UNIVAC 1100-series or CYBER 175.

Source record↗
NASA NTRS2003

Optimal Feedback Control of Thermal Networks

An improved approach to the mathematical modeling of feedback control of thermal networks has been devised. Heretofore software for feedback control of thermal networks has been developed by time-consuming trial-and-error methods that depend on engineers expertise. In contrast, the present approach is a systematic means of developing algorithms for feedback control that is optimal in the sense that it combines performance with low cost of implementation. An additional advantage of the present approach is that a thermal engineer need not be expert in control theory. Thermal networks are lumped-parameter approximations used to represent complex thermal systems. Thermal networks are closely related to electrical networks commonly represented by lumped-parameter circuit diagrams. Like such electrical circuits, thermal networks are mathematically modeled by systems of differential-algebraic equations (DAEs) that is, ordinary differential equations subject to a set of algebraic constraints. In the present approach, emphasis is placed on applications in which thermal networks are subject to constant disturbances and, therefore, integral control action is necessary to obtain steady-state responses. The mathematical development of the present approach begins with the derivation of optimal integral-control laws via minimization of an appropriate cost functional that involves augmented state vectors. Subsequently, classical variational arguments provide optimality conditions in the form of the Hamiltonian equations for the standard linear-quadratic-regulator (LQR) problem. These equations are reduced to an algebraic Riccati equation (ARE) with respect to the augmented state vector. The solution of the ARE leads to the direct computation of the optimal proportional- and integral-feedback control gains. In cases of very complex networks, large numbers of state variables make it difficult to implement optimal controllers in the manner described in the preceding paragraph.

Papalexandris, Miltiadis↗
NASA NTRS1983

Large-area sheet task: Advanced dendritic-web-growth development

Thermally generated stresses in the growing web crystal were reduced. These stresses, which if too high cause the ribbon to degenerate, were reduced by a factor of three, resulting in the demonstrated growth of high-quality web crystals to widths of 5.4 cm. This progress was brought about chiefly by the application of thermal models to the development of low-stress growth configurations. A new temperature model was developed which can analyze the thermal effects of much more complex lid and top shield configurations than was possible with the old lumped shield model. Growth experiments which supplied input data such as actual shield temperature and melt levels were used to verify the modeling results. Desirable modifications in the melt level-sensing circuitry were made in the new experimental web growth furnace, and this furnace has been used to carry out growth experiments under steady-state conditions. New growth configurations were tested in long growth runs at Westinghouse AESD which produced wider, lower stress and higher quality web crystals than designs previously used.

Duncan, C. S.↗
NASA NTRS1992

Preliminary investigation of a low power pulsed arcjet thruster

A type of pulsed arcjet is examined which operates with gaseous helium propellant at powers from 100 to 1500 W and pulse rats from 360 to 6000 pulses per second. During the pulse, peak power is 50 to 250 kW generating chamber pressures of 10 to 50 atmospheres. The high operating pressure substantially reduces ionization in the 2.5 mm diameter x 12.5 mm long capillary and lowers frozen flow losses in the nozzle. The thruster thermal efficiency, measured calorimetrically, is 43 percent at 440 W and 8 mg/sec helium mass flow rate. Thruster performance trends are predicted by a time-dependent lumped-parameter model which includes heat exchange between the propellant and the wall. The model substantially underpredicts the experimental thermal efficiency values. Alternative thruster configurations yielding significant improvements in thrust efficiency and specific impulse are discussed.

Taylor, Russell D.↗
NASA NTRS1971

Thermal modeling of a cryogenic turbopump for space shuttle applications.

Thermal modeling of a cryogenic pump and a hot-gas turbine in a turbopump assembly proposed for the Space Shuttle is described in this paper. A model, developed by identifying the heat-transfer regimes and incorporating their dependencies into a turbopump system model, included heat transfer for two-phase cryogen, hot-gas (200 R) impingement on turbine blades, gas impingement on rotating disks and parallel plate fluid flow. The ?thermal analyzer' program employed to develop this model was the TRW Systems Improved Numerical Differencing Analyzer (SINDA). This program uses finite differencing with lumped parameter representation for each node. Also discussed are model development, simulations of turbopump startup/shutdown operations, and the effects of varying turbopump parameters on the thermal performance.

Knowles, P. J.↗
NASA NTRS1990

Analysis of the thermal performance of heat pipe radiators

A comprehensive mathematical model and computational methodology are presented to obtain numerical solutions for the transient behavior of a heat pipe radiator in a space environment. The modeling is focused on a typical radiator panel having a long heat pipe at the center and two extended surfaces attached to opposing sides of the heat pipe shell in the condenser section. In the set of governing equations developed for the model, each region of the heat pipe - shell, liquid, and vapor - is thermally lumped to the extent possible, while the fin is lumped only in the direction normal to its surface. Convection is considered to be the only significant heat transfer mode in the vapor, and the evaporation and condensation velocity at the liquid-vapor interface is calculated from kinetic theory. A finite-difference numerical technique is used to predict the transient behavior of the entire radiator in response to changing loads.

Boo, J. H.↗
NASA NTRSDate not supplied

Single and Multi-Node Modeling of Direct, Submerged, and Self-Pressurization of A Cryogenic Propellant Tank Using Nodal Tools

The pressurization of cryogenic propellant tanks will always be an important process so long as cryogenic liquids are being considered as fuel sources or used for other in-space applications. Pressure control of the tank ullage is necessary for managing propellant flowrates to an engine or a receiver tank, and modeling of the process is used to predict the pressurant requirements and the amount of propellant boiloff. Direct ullage pressurization is the more traditional approach to tank pressurization, as the physics are straight-forward, and ample test (flight) data have been collected and analyzed over the past several decades. Submerged injection pressurization is an alternate method for tank pressurization and has been shown to reduce pressurant requirements, subcool the propellant, and reduce the risk of ullage collapse. Additionally, the pressurant gas entering the ullage is usually much colder when using the submerged pressurization approach, resulting in reduced propellant boiloff. These benefits are at the expense of vaporizing a small percentage of the propellent. Both tank pressurization methods are viable options for current and future space missions, and it is important to have the capability of analyzing the tank ullage conditions for both approaches. Our previous work has demonstrated the development of a Generalized Fluid System Simulation Program (GFSSP) model, which contains a thermodynamic equilibrium heat and mass transfer subroutine capable of effectively analyzing both direct and submerged pressurization systems [1-2]. This subroutine has most recently been enhanced to include the non-equilibrium effect of pressurant dissolution into the propellant. To date the ullage has always been represented as a single node, and although the simulated single-node temperatures have good comparison with the volume-averaged ullage temperatures computed from test data, the physics of the thermal stratification in the ullage were never captured, and adjustment factors in the model were required. The purpose of this paper is to introduce the development of a multi-node ullage model using GFSSP and to discuss the improvements of the simulated ullage temperature distribution and its resulting effects on ullage heat transfer processes. Test data from the Cryogenic Propellant Storage and Transfer Engineering Developmental Unit (CPST EDU) was used for model validation. For additional comparison, a Thermal Desktop (TD) model was also developed to analyze the CPST EDU direct ullage pressurization tests using both a single node and multi-node approach. The model includes the direct pressurant line, vent line, fill/drain line, and a TD FloCAD Compartment. The TD FloCAD Compartment is employed to represent the liquid and ullage as single volumes inside the tank, to include a liquid/vapor interface, and to generate network level objects such as lumps (analogous to nodes in GFSSP), paths, and ties between the fluid and thermal elements. An established heat load on the model tank was leveraged from a pre-existing higher-fidelity model correlated to CPST EDU test data.

pressurization↗
NASA NTRSDate not supplied

Predicting Maximum Thermal Response in a Li-ion Cell as a Thermal Runaway Predictor for a UAV Fight

As the energy storage devices continue to "pack" more energy in a small space, any damage, battery component failure, manufacturing defect, or electrically abusing the battery can lead to catastrophic thermal runaway events. A catastrophic thermal event in a cell leads to high temperature, in some instances to spewing of battery materials due to gas development from side reactions initiated due to high internal temperatures. Moreover, a thermal runaway can propagate from a single "failed" cell to a whole battery pack, resulting in a more serious event. The commercial and automotive sectors need to mitigate thermal runaway events. An electric aircraft (or air taxi) has no alternative in the event of a failure, so preventing such events is paramount. Currently, battery prognostics algorithms only predict the state-of-charge (SOC) and end-of-life (EOL) of a Li-ion battery in a UAV (unmanned air vehicle), but do not predict the maximum temperature during a flight or the likelihood of thermal runaway. Our current work focuses on adding the thermal model (based on Birkl's OCV model, in Ref. [1]) to the existing hybrid electrochemical model (single-particle model with lumped parameters as in Ref. [2]) to determine end-of-discharge, end-of-life, and maximum temperature during discharge.

Thermal runaway↗
NASA NTRS1993

Program Computes Flows Of Fluids And Heat

SINDA'85/FLUINT incorporates lumped-parameter-network and one-dimensional-flow mathematical models. System enables analysis of mutual influences of thermal and flow phenomena. Offers two finite-difference numerical solution techniques: forward-difference explicit approximation and Crank-Nicholson approximation. Enables simulation of nonuniform heating and facilitates mathematical modeling of thin-walled heat exchangers. Ability to model nonequilibrium behavior within two-phase volumes included. Recent changes in program improve modeling of real evaporator pumps and other capillary-assist evaporators. Written in FORTRAN 77.

Cullimore, Brent↗
NASA NTRS1992

A simplified generic cryostat thermal model for predicting cryogen mass and lifetime

This study presents a simple generic cryostat thermal model developed for predicting cryogen mass and the lifetime of cryogenic space telescopes. The model is based on a lumped parameter representation of eight nodes and over 35 conductors. The major telescope Dewar components represented as nodes are the main cryogen tank, the three vapor-cooled shields, the outer shell, barrel baffle, and the barrel baffle heat exchanger. The input to the model consists of the vapor-cooled shield, support straps and instrument cable geometry, and the cryostat heat loads from the instrument and through the aperture. The predictions from the model were compared with those from the more detailed models of IRAS, COBE, and SIRTF. Correlation of the helium flow rates of the IRAS and COBE Dewars was 12 percent above the actual flight or ground test data and for SIRTF it was less than 2 percent of the prediction from the detailed model.

Birur, G. C.↗
NASA NTRS1980

On the performance of explicit and implicit algorithms for transient thermal analysis

The status of an effort to increase the efficiency of calculating transient temperature fields in complex aerospace vehicle structures is described. The advantages and disadvantages of explicit and implicit algorithms are discussed. A promising set of implicit algorithms, known as the GEAR package is described. Four test problems, used for evaluating and comparing various algorithms, have been selected and finite element models of the configurations are discribed. These problems include a space shuttle frame component, an insulated cylinder, a metallic panel for a thermal protection system and a model of the space shuttle orbiter wing. Calculations were carried out using the SPAR finite element program, the MITAS lumped parameter program and a special purpose finite element program incorporating the GEAR algorithms. Results generally indicate a preference for implicit over explicit algorithms for solution of transient structural heat transfer problems when the governing equations are stiff. Careful attention to modeling detail such as avoiding thin or short high-conducting elements can sometimes reduce the stiffness to the extent that explicit methods become advantageous.

Adelman, H. M.↗
NASA NTRSDate not supplied

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

Additive Manufacturing↗
NASA NTRSDate not supplied

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

pragmatic method↗
NASA NTRS1992

Semi-coupled flow and thermal analysis of the field joint during rapid pressurization of the redesigned solid rocket motor

A transient, semi-coupled, multi-dimensional thermal and flow analysis methodology was developed to predict the thermal/gas dynamic conditions in the field joint region of the Redesigned Solid Rocket Motor (RSRM). Transient temperature response, pressure history, and combustion gas flow rate (within the field joint region), were of principle interest, in the course of this study. The thermal environment in the field joint was modeled using SINDA, a finite difference based thermal network analyzer. The combustion gas flow boundary condition was generated using the FLAP code; this code performs a transient, lumped-parameter, control volume analysis to solve the mass, momentum, and energy conservation equations. The FLAP computer code was modified to account for erosion of the NBR insulation material, following ignition. An independent grid sensitivity study was conducted to determine an appropriate grid distribution near the wall. The predicted results, obtained using an optimum grid distribution and computer generated flow boundary condition, were compared with subscale test data.

Ghaffarian, Benny↗
NASA NTRS2018

Thermal Modeling of Zero Boil Off Tank Experiment

Understanding fluid behavior in microgravity environments is essential to further development of cryogenic storage in space environments. The Zero Boil Off Tank (ZBOT) experiment was designed to investigate two-phase pressurization and depressurization of a tank in a microgravity environment. The test fluid was the refrigerant Perfluoro-normal-Pentane (PNP). Thermal modeling for the ZBOT model was conducted using Thermal Desktop and SINDA/FLUINT. The temperature distribution within the fluid of the tank is of particular interest. This particular work is centered on ascertaining the thermal behavior of the refrigerant in order to build more complete models of fluid in microgravity. Separate cases were run modeling experiments that were conducted both on the ground and on the International Space Station (ISS) to compare 1g and microgravity environments. The microgravity modeling cases consisted of a fluid lump representing the vapor ullage suspended in a solid to represent the liquid. Mass flow between the liquid and vapor was modeled using the Schrage equation for mass flow. Initial results indicate that the pressure rise and temperature increase within the fluid closely align with the experimental data by matching initial conditions of the experiment. This work is ongoing and will yield further insights into the thermal behavior of fluid mixing in microgravity.

low-gravity fluid modeling↗
NASA NTRS2001

Thermal Analysis of a Carbon Fiber Rope Barrier for Use in the Reusable Solid Rocket Motor Nozzle Joint-2

This study provides for development and verification of analysis methods used to assess performance of a carbon fiber rope (CFR) thermal barrier system that is currently being qualified for use in Reusable Solid Rocket Motor (RSRM) nozzle joint-2. Modeled geometry for flow calculations considers the joint to be vented with the porous CFR barriers placed in the "open' assembly gap. Model development is based on a 1-D volume filling approach where flow resistances (assembly gap and CFRs) are defined by serially connected internal flow and the porous media "Darcy" relationships. Combustion gas flow rates are computed using the volume filling code by assuming a lumped distribution total joint fill volume on a per linear circumferential inch basis. Gas compressibility, friction and heat transfer are included in the modeling. Gas-to-wall heat transfer is simulated by concurrent solution of the compressible flow equations and a large thermal 2-D finite element (FE) conduction grid. The derived numerical technique loosely couples the FE conduction matrix with the compressible gas flow equations, Free constants that appear in the governing equations are calibrated by parametric model comparison to hot fire subscale test results. The calibrated model is then used to make full-scale motor predictions using RSRM aft dome environments. Model results indicate that CFR thermal barrier systems will provide a thermally benign and controlled pressurization environment for the RSRM nozzle joint-2 primary seal activation.

Clayton, J. Louie↗