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.
Tables and graphs were compiled for solving nonlinear heat transfer constant area fin problems. The differential equation describing one-dimensional steady-state temperature distribution and heat flow under three modes of heat transfer with heat generation was investigated.
The use of classical optimal control methods, in particular variational methods, to solve the airfoil optimization problem, by deriving a set of adjoint (costate) equations and boundary conditions has already been done for inviscid (potential and Euler flows) and two dimensional, steady state, incompressible flow governed by the Navier-Stokes equations. The interior and boundary terms of the volume integral have been derived (in this work) for the steady Navier-Stokes equations in three dimensions for a viscous, compressible heat conducting fluid. This can be used to derive the adjoint equations and numerical boundary conditions for general classes of problems and hence paves the way for a solution to the aerodynamic optimization problem for compressible viscous flows. The next steps to the realization of that goal are projected as below. The usual square integral pressure functional as an objective function is being replaced by a more realistic drag functional subject to a lift constraint. The feasibility of attempting the more difficult time dependent problem is being investigated. It remains to get the full system of adjoint equations and boundary conditions with the new functional. The state and adjoint equations must be discretized and coded. An appropriate optimization program must be used (steepest descents seems inadequate) and various known airfoil shapes should be recovered in test cases of the computer program.
Cost is equally important to power density or efficiency for the adoption of waste heat recovery thermoelectric generators (TEG) in many transportation and industrial energy recovery applications. In many cases the system design that minimizes cost (e.g., the $/W value) can be very different than the design that maximizes the system's efficiency or power density, and it is important to understand the relationship between those designs to optimize TEG performance-cost compromises. Expanding on recent cost analysis work and using more detailed system modeling, an enhanced cost scaling analysis of a waste heat recovery thermoelectric generator with more detailed, coupled treatment of the heat exchangers has been performed. In this analysis, the effect of the heat lost to the environment and updated relationships between the hot-side and cold-side conductances that maximize power output are considered. This coupled thermal and thermoelectric treatment of the exhaust waste heat recovery thermoelectric generator yields modified cost scaling and design optimization equations, which are now strongly dependent on the heat leakage fraction, exhaust mass flow rate, and heat exchanger effectiveness. This work shows that heat exchanger costs most often dominate the overall TE system costs, that it is extremely difficult to escape this regime, and in order to achieve TE system costs of $1/W it is necessary to achieve heat exchanger costs of $1/(W/K). Minimum TE system costs per watt generally coincide with maximum power points, but Preferred TE Design Regimes are identified where there is little cost penalty for moving into regions of higher efficiency and slightly lower power outputs. These regimes are closely tied to previously-identified low cost design regimes. This work shows that the optimum fill factor Fopt minimizing system costs decreases as heat losses increase, and increases as exhaust mass flow rate and heat exchanger effectiveness increase. These findings have profound implications on the design and operation of various thermoelectric (TE) waste heat 3 recovery systems. This work highlights the importance of heat exchanger costs on the overall TEG system costs, quantifies the possible TEG performance-cost domain space based on heat exchanger effects, and provides a focus for future system research and development efforts.
Cost is just as important as power density or efficiency for the adoption of waste heat recovery thermoelectric generators (TEG). Prior work [1] has shown that the system design that minimizes cost (e.g., the $/W value) can be close to the designs that maximize the system’s efficiency or power density, however, it is important to understand the relationship between those designs to optimize TEG performance-cost compromises. Expanding on recent work [1, 2, 3] the impact of heat exchanger conditions on the optimum TEG fill factors and cost scaling of a waste heat recovery thermoelectric generator with a detailed treatment of the hot side exhaust heat exchanger has been investigated further. The effect of the heat lost to the environment and updated relationships between the hot-side and cold-side conductances [4] that maximize power output are considered. The optimum fill factor to minimize TEG energy recovery system costs is strongly dependent on the heat leakage fraction, σ, the mass flow rate of the exhaust, the hot-side heat exchanger effectiveness, heat exchanger UAh, and heat flux. These relationships are explored and characterized for typical exhaust gas-flow conditions to show the inherent design complexities. The heat exchanger costs often dominate the TEG cost equation and it is critical to fully understand the tradeoff between heat exchanger performance, optimum TEG fill factors, and cost to establish potentially optimum design points within the cost-performance design space. This work will explore the design tradeoffs and relationships within the cost-efficiency-power density design space for a typical thermoelectric energy recovery system application. The interplay between optimum TEG fill factors and heat exchanger design can impact system footprint, volume, and mass in weight-sensitive applications. Less-effective, low-cost heat exchangers may outperform higher cost alternatives from a market adoption perspective. This shift of emphasis acknowledging the interdependence of optimum TEG fill factors and heat exchanger performance has significant implications on thermoelectric waste heat recovery systems designs and their operation. In addition, preferred TEG design regimes exist that accommodate reasonable compromises in TE performance and cost. This effort highlights how the optimum fill factor–heat exchanger performance relations couple to these optimum TEG performance-cost domains based on TEG-system-level analyses and provides a focus for future system research and development efforts.
The phonon Boltzmann transport equation (BTE) with the relaxation time approximation (RTA) scattering model is used to calculate transient temperature profiles in graphene, and the results are compared to analogous molecular dynamics (MD) simulations. For the BTE calculations, the phonon dispersion relation and frequency-dependent scattering rates are obtained from a combination of MD data and semi-empirical power-law expressions for the normal and Umklapp phonon lifetimes. The dimensions and initial temperature conditions of graphene are varied to study the size and temperature dependence of thermal transport physics at the mesoscopic scale. Good quantitative agreement to within 5% is found between the BTE and MD results, over a wide range of temperatures and lengthscales of the temperature variation in the graphene sheet. Small differences are attributed to the inaccuracy of the RTA as applied to graphene, and to neglecting four-phonon scattering in the BTE simulations. The present results may further understanding in applications such as the transient heating of nanoelectronics.
The vapor flow in a heat pipe was mathematically modeled and the equations governing the transient behavior of the core were solved numerically. The modeled vapor flow is transient, axisymmetric (or two-dimensional) compressible viscous flow in a closed chamber. The two methods of solution are described. The more promising method failed (a mixed Galerkin finite difference method) whereas a more common finite difference method was successful. Preliminary results are presented showing that multi-dimensional flows need to be treated. A model of the liquid phase of a high temperature heat pipe was developed. The model is intended to be coupled to a vapor phase model for the complete solution of the heat pipe problem. The mathematical equations are formulated consistent with physical processes while allowing a computationally efficient solution. The model simulates time dependent characteristics of concern to the liquid phase including input phase change, output heat fluxes, liquid temperatures, container temperatures, liquid velocities, and liquid pressure. Preliminary results were obtained for two heat pipe startup cases. The heat pipe studied used lithium as the working fluid and an annular wick configuration. Recommendations for implementation based on the results obtained are presented. Experimental studies were initiated using a rectangular heat pipe. Both twin beam laser holography and laser Doppler anemometry were investigated. Preliminary experiments were completed and results are reported.
Methodologies have been developed for modeling both gas dynamics and heat transfer inside the carbon fiber rope (CFR) for applications in the space shuttle reusable solid rocket motor joints. Specifically, the CFR is modeled using an equivalent rectangular duct with a cross-section area, friction factor and heat transfer coefficient such that this duct has the same amount of mass flow rate, pressure drop, and heat transfer rate as the CFR. An equation for the friction factor is derived based on the Darcy-Forschheimer law and the heat transfer coefficient is obtained from pipe flow correlations. The pressure, temperature and velocity of the gas inside the CFR are calculated using the one-dimensional Navier-Stokes equations. Various subscale tests, both cold flow and hot flow, have been carried out to validate and refine this CFR model. In particular, the following three types of testing were used: (1) cold flow in a RSRM nozzle-to-case joint geometry, (2) cold flow in a RSRM nozzle joint No. 2 geometry, and (3) hot flow in a RSRM nozzle joint environment simulator. The predicted pressure and temperature history are compared with experimental measurements. The effects of various input parameters for the model are discussed in detail.
Implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solid objects subjected to nonlinear boundary conditions is shown to generate an accurate and economical approach for the computation of shape sensitivities for this class of problems. This approach involves the employment of analytical derivatives of boundary-element kernel functions with respect to shape design variables. A formulation is presented that can consistently account for both temperature-dependent convection and radiation boundary conditions. Several iterative strategies are presented for the solution of the resulting sets of nonlinear equations and the computational performances examined in detail. Multizone analysis and zone condensation strategies are demonstrated to provide substantive computational economies in this process for models with either localized nonlinear boundary conditions or regions of geometric insensitivity to design variables. A series of nonlinear example problems are presented that have closed-form solutions.
The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.
The flowfield and resulting heat transfer rate over a series of ice accretion shapes is obtained by solving the Navier-Stokes equations. The influence of surface roughness on surface heat transfer is examined by including blockage, form drag, and stagnation heating effects as source terms in the governing equations. The results indicate increases of a factor of three in cooling rates due to distributed roughness compared to smotth surfaces. In addition, droplet impingement efficiencies are studied for the same series of ice accretion shapes using a time-dependent solution procedure.
I will present models of the solar atmospheric response to flare heating. The models solve the equations of non-LTE radiation hydrodynamics with an electron beam added as a flare energy source term. Radiative transfer is solved in detail for many important optically thick hydrogen and helium transitions and numerous optically thin EUV lines making the models ideally suited to study the emission that is produced during flares. I will pay special attention to understanding key EUV lines as well the mechanism for white light production. I will also present preliminary results of how the model solar atmosphere responds to Fletcher & Hudson type flare heating. I will compare this with the results from flare simulations using the standard thick target model.
The adjustment of coronal active-region loops to changes in their heating rate is investigated numerically. The one-dimensional hydrodynamic equations are solved subject to boundary conditions in which heat flux-induced mass exchange between coronal and chromospheric components is allowed. The calculated evolution of physical parameters suggests that (1) mass supplied during chromospheric evaporation is much more effective in moderating coronal temperature excursions than when downward heat flux is dissipated by a static chromosphere, and (2) the method by which the chromosphere responds to changing coronal conditions can significantly influence coronal readjustment time scales. Observations are cited which illustrate the range of possible fluctuations in the heating rates.
On the moon, the available heat sink temperature for a vertical unshaded radiator at the equator is 322 K. A method of reducing this heat sink temperature using a parabolic trough shading device was investigated. A steady state heat balance was performed to predict the available heat sink temperature. The effect of optical surface properties on system performance was investigated. Various geometric configurations were also evaluated. A flexible shade conceptual design is presented which greatly reduces the weight and stowed volume of the system. The concept makes use of the natural catenary shape assumed by a flexible material when supported at two points. The catenary shape is very near parabolic. The lunar radiator parabolic shade design presented integrates the energy collection and rejection of a solar dynamic power cycle with the moderate temperature waste heat rejection of a lunar habitat.
Ordinary chondrites are considered to have experienced thermal metamorphism in small bodies. We are interested in behaviors of volatile elements in such a kind of thermally evolving planetesimals. Volatile elements generally have high vapor pressures at high temperature. In porous bodies, with a high gas permeability, volatile elements are transported efficiently over a long range. Behavior of volatile elements transported by permeable gas flow can be handled by an equation whose form is similar to that of the equation of thermal diffusion. We can follow transportation of heats and volatile elements in planetesimals, when parameters in these equations, initial conditions and chemical behavior of volatile elements are given. Recently, we discovered that nitrogen in equilibrated H-chondrites is mainly trapped in taenite (f.c.c. Fe-Ni), probably dissolved in interstitial sites. Fegley suggests that metallic iron cannot trap nitrogen in the solar nebula gas due to its very low nitrogen partial pressure. Approximately 1 bar of nitrogen pressure is required to explain the nitrogen content in taenite. We may expect high nitrogen gas partial pressure (possibly produced by vaporization of nitrogen-bearing solids such as organic materials) at the interior of thermally evolving planetesimals. Kinetic behavior of nitrogen in taenite suggests that it can easily be equilibrated with the ambient nitrogen gas at temperatures of approximately 500 C or higher. We consider that nitrogen is trapped in taenite through a nitrogen redistribution process occurred during the thermal metamorphic event.
The two-dimensional spatially elliptic Navier-Stokes equations have been used to investigate the radiative interactions in chemically reacting compressible flows of premixed hydrogen and air in an expanding nozzle. The radiative heat transfer term in the energy equation is simulated using the Monte Carlo method (MCM). The nongray model employed is based on the statistical narrow band model with an exponential-tailed inverse intensity distribution. The spectral correlation has been considered in the Monte Carlo formulations. Results obtained demonstrate that the effect of radiation on the flow field is minimal but its effect on the wall heat transfer is significant. Extensive parametric studies are conducted to investigate the effects of equivalence ratio, wall temperature, inlet flow temperature, and the nozzle size on the radiative and conductive wall fluxes.
Experiments were conducted to determine the effects of buoyancy and Coriolis forces on heat transfer in turbine blade internal coolant passages. The experiments were conducted with a large scale, multi-pass, smooth-wall heat transfer model with both radially inward and outward flow. An analysis of the governing flow equations showed that four parameters influence the heat transfer in rotating passages (coolant-to-wall temperature ratio, Rossby number, Reynolds number and radius-to-passage hydraulic diameter ratio). These four parameters were varied over ranges which are typical of advanced gas turbine engine operating conditions. It was found that both Coriolis and buoyancy effects must be considered in turbine blade cooling designs and that the effect of rotation on the heat transfer coefficients was markedly different depending on the flow direction. Local heat transfer coefficients were found to decrease by as much as 60 percent and increase by 250 percent from no rotation levels. Comparisons with a pioneering stationary vertical tube buoyancy experiment showed reasonably good agreement. Correlation of the data is achieved employing dimensionless parameters derived from the governing flow equations.
An experimental and simulation study of warm dense carbon foams at ambient density (n e ~10 21 cm –3 ) is presented. This study of isochorically heated foams is motivated by their potential application in carbon-atmosphere white-dwarf envelopes, where there are modeling uncertainties due to the equation of state. The foams are heated on an approximately picosecond time scale with a laser-accelerated proton beam. The cooling and expansion of the heated foams can be modeled with appropriately initialized radiation-hydrodynamics codes; xRAGE code is used in this work. The primary experimental diagnostic is the streaked optical pyrometer, which images a narrow band of radiation from the rear surface of the heated material. Presented are xRAGE modeling results for both solid aluminum targets and carbonized resorcinol-formaldehyde foam targets, showing that the foam appears to cool slowly on the pyrometer because of partial transparency. So that simulations of cooling foam are processed properly, it is necessary to account for finite optical depth in the photosphere calculation, and the methods for performing that calculation are presented in depth.
Predicting temperature rise accurately during fast charge/discharge of a Li-ion cell is essential to avoid thermal runaway and extend battery life. While modeling the temperature rise, it is necessary to model the heat generation correctly. Heat of mixing, one of the four main sources of heat generation in a Li-ion battery, has often been considered insignificant and therefore excluded from modeling. When included, it is modeled using the expression from a Taylor expansion approximation. In this work, we have shown the conditions when including the heat of mixing becomes important and quantified the error associated with using the Taylor expansion, especially under high charge/discharge conditions. Consequently, we carry out the calculation of the rate of enthalpy change and the heat generation rate from the most fundamental equation for the rate of the total enthalpy change without simplifying assumptions or approximations. The heat generation rate calculated doing so naturally includes irreversible, reversible and mixing heat. We then exclusively separate out the heat of mixing by subtracting the rate of enthalpy change by reaction from the rate of the total enthalpy change. Results show that the contribution of heat of mixing in the total heat generated increases with the charge/discharge rate and is as large as 23% for a 6 C discharge. This result suggests that while modeling heat generation for fast charge/discharge, it is necessary to include the heat of mixing and avoid calculating it using the Taylor expansion approximation.