Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “heat equation”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Three-dimensional MHD modeling of magnetotail dynamics for different polytropic indices

The paper analyzes the 3D resistive tearing instability in a magnetotail configuration by means of a resistive MHD code for various forms of the energy equation, including ohmic and compressional heating. A complete derivation of the set of MHD equations from a kinetic formalism is presented, with all simplifying assumptions noted. Emphasis is placed on the changes in magnetotail evolution caused by a variation of the assumed value of the polytropic exponent gamma, representing the ratio of specific heat in the energy equation which closes the set of MHD equations. Within the gamma range of 2/3-4, a tearing mode is found to develop in all cases. Certain properties, such as the formation and ejection of a plasmoid, associated with the formation of the reversal region of the north-south component of the magnetic field, and strong earthward and tailward flow fields prevail as qualitative features independent of the value of gamma chosen. Ohmic dissipation is found to be the major source of the heating, while adiabatic or nonadiabatic compression may contribute primarily in the closed field region earthward of the reconnection site.

Hesse, Michael↗
NASA NTRS1994

Approximate Solution Methods for Spectral Radiative Transfer in High Refractive Index Layers

Some ceramic materials for high temperature applications are partially transparent for radiative transfer. The refractive indices of these materials can be substantially greater than one which influences internal radiative emission and reflections. Heat transfer behavior of single and laminated layers has been obtained in the literature by numerical solutions of the radiative transfer equations coupled with heat conduction and heating at the boundaries by convection and radiation. Two-flux and diffusion methods are investigated here to obtain approximate solutions using a simpler formulation than required for exact numerical solutions. Isotropic scattering is included. The two-flux method for a single layer yields excellent results for gray and two band spectral calculations. The diffusion method yields a good approximation for spectral behavior in laminated multiple layers if the overall optical thickness is larger than about ten. A hybrid spectral model is developed using the two-flux method in the optically thin bands, and radiative diffusion in bands that are optically thick.

Siegel, R.↗
NASA NTRS1955

Calculations of Laminar Heat Transfer Around Cylinders of Arbitrary Cross Section and Transpiration-Cooled Walls with Application to Turbine Blade Cooling

An approximate method for the development of flow and thermal boundary layers in the laminar region on cylinders with arbitrary cross section and transpiration-cooled walls is obtained by the use of Karman's integrated momentum equation and an analogous heat-flow equation. Incompressible flow with constant property values throughout the boundary layer is assumed. The velocity and temperature profiles within the boundary layer are approximated by expressions composed of trigonometric functions. Shape parameters for these profiles and functions necessary for the solution of the boundary-layer equations are presented as graphs so that the calculation for any specific case is reduced to the solution of two first-order differential equations. This method is applied to determine local heat-transfer coefficients and surface temperatures in the laminar flow region of the transpiration-cooled turbine blades for a given coolant flow rate, or to calculate the coolant flow distribution which is necessary in order to keep the blade temperature uniform along the surface.

Eckert, E R G↗
DOE OSTI2022

Unstable turbulent channel flow response to spanwise-heterogeneous heat fluxes: Prandtl's secondary flow of the third kind

Turbulent secondary flows are defined as Prandtl's secondary flow of the first or second kind, the former produced by stretching and/or tilting of vorticity, the latter produced via spatial heterogeneity of Reynolds stresses. Both mechanisms are instantaneously active within inertia-dominated wall turbulence; Reynolds stress spatial heterogeneity is required for Reynolds-averaged secondary flows. Spanwise-variable surface roughness can induce turbulent stress spatial heterogeneity in the spanwise–wall-normal plane and provide sustenance for streamwise-aligned mean secondary flows. Herein, we demonstrate that turbulent secondary flows can also be sustained by spanwise variability in the surface heat flux in unstably stratified turbulent channels, defined hereafter as Prandtl's secondary flow of the third kind. Support for this mechanism is established with scaling arguments, while large-eddy simulation is used to model inertia-dominated channel turbulence responding to a lower boundary with uniform aerodynamic/hydrodynamic roughness but spanwise-variable surface heat flux. Transport equations for streamwise vorticity and turbulent kinetic energy, k, outline the conditions needed for third-kind production: shear and buoyancy production over the elevated heat flux regions necessitates lateral entrainment of low- k fluid, inducing mean counter-rotating secondary cells aligned such that upwelling and downwelling occur over the high and low heat flux regions, respectively. Buoyancy-driven production of k alters aggregate flow response and thus is a distinctly different mechanism responsible for sustenance of secondary flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1991

Heat flow vs. atmospheric greenhouse on early Mars

Researchers derived a quantitative relationship between the effectiveness of an atmospheric greenhouse and internal heat flow in producing the morphological differences between earlier and later Martian terrains. The derivation is based on relationships previously derived by other researchers. The reasoning may be stated as follows: the CO2 mean residence time in the Martian atmosphere is almost certainly much shorter than the total time span over which early climate differences are thought to have been sustained. Therefore, recycling of previously degassed CO2 quickly becomes more important than the ongoing supply of juvenile CO2. If so, then the atmospheric CO2 pressure, and thereby the surface temperature, may be approximated mathematically as a function of the total degassed CO2 in the atmosphere plus buried material and the ratio of the atmospheric and regolith mean residence times. The latter ratio can also be expressed as a function of heat flow. Hence, it follows that the surface temperature may be expressed as a function of heat flow and the total amount of available CO2. However, the depth to the water table can simultaneously be expressed as a function of heat flow and the surface temperature (the boundary condition). Therefore, for any given values of total available CO2 and regolith conductivity, there exist coupled independent equations which relate heat flow, surface temperature, and the depth to the water table. This means we can now derive simultaneous values of surface temperature and the depth of the water table for any value of the heat flow. The derived relationship is used to evaluate the relative importance of the atmospheric greenhouse effect and the internal regolith thermal gradient in producing morphological changes for any value of the heat flow, and to assess the absolute importance of each of the values of the heat flow which are thought to be reasonable on independent geophysical grounds.

Fanale, F. P.↗
NASA NTRS1989

A one- and two-layer model for estimating evapotranspiration with remotely sensed surface temperature and ground-based meteorological data over partial canopy cover

Surface-air temperature differences are commonly used in a bulk resistance equation for estimating sensible heat flux (H), which is inserted in the one-dimensional energy balance equation to solve for the latent heat flux (LE) as a residual. Serious discrepancies between estimated and measured LE have been observed for partial-canopy-cover conditions, which are mainly attributed to inappropriate estimates of H. To improve the estimates of H over sparse canopies, one- and two-layer resistance models that account for some of the factors causing poor agreement are developed. The utility of the two models is tested with remotely sensed and micrometeorological data for a furrowed cotton field with 20 percent cover and a dry soil surface. It is found that the one-layer model performs better than the two-layer model when a theoretical bluff-body correction for heat transfer is used instead of an empirical adjustment; otherwise, the two-layer model is better.

Kustas, William P.↗
DOE OSTI2021

Data-Driven Reduced-Order Modeling of Convective Heat Transfer in Porous Media

This work presents a data-driven Reduced-Order Model (ROM) for parametric convective heat transfer problems in porous media. The intrusive Proper Orthogonal Decomposition aided Reduced-Basis (POD-RB) technique is employed to reduce the porous medium formulation of the incompressible Reynolds-Averaged Navier–Stokes (RANS) equations coupled with heat transfer. Instead of resolving the exact flow configuration with high fidelity, the porous medium formulation solves a homogenized flow in which the fluid-structure interactions are captured via volumetric flow resistances with nonlinear, semi-empirical friction correlations. A supremizer approach is implemented for the stabilization of the reduced fluid dynamics equations. The reduced nonlinear flow resistances are treated using the Discrete Empirical Interpolation Method (DEIM), while the turbulent eddy viscosity and diffusivity are approximated by adopting a Radial Basis Function (RBF) interpolation-based approach. The proposed method is tested using a 2D numerical model of the Molten Salt Fast Reactor (MSFR), which involves the simulation of both clean and porous medium regions in the same domain. For the steady-state example, five model parameters are considered to be uncertain: the magnitude of the pumping force, the external coolant temperature, the heat transfer coefficient, the thermal expansion coefficient, and the Prandtl number. For transient scenarios, on the other hand, the coastdown-time of the pump is the only uncertain parameter. The results indicate that the POD-RB-ROMs are suitable for the reduction of similar problems. The relative L2 errors are below 3.34% for every field of interest for all cases analyzed, while the speedup factors vary between 54 (transient) and 40,000 (steady-state).

German, Péter (ORCID:0000000307285283)↗
NASA NTRS1951

Calculations of Laminar Heat Transfer Around Cylinders of Arbitrary Cross Section and Transpiration-Cooled Walls with Application to Turbine Blade Cooling

An approximate method for development of flow and thermal boundary layers in laminar regime on cylinders with arbitrary cross section and transpiration-cooled walls is obtained by use of Karman's integrated momentum equation and an analogous heat-flow equation. Incompressible flow with constant property values throughout boundary layer is assumed. Shape parameters for approximated velocity and temperature profiles and functions necessary for solution of boundary-layer equations are presented as charts, reducing calculations to a minimum. The method is applied to determine local heat-transfer coefficients and surface temperature-cooled turbine blades for a given flow rate. Coolant flow distributions necessary for maintaining uniform blade temperatures are also determined.

TURBINE COOLING↗
NASA NTRS1989

Convergence of infinite dimensional sampled LQR problems - Theory and numerical results

A theory is developed for the convergence of the closed-loop solution to infinite-dimensional discrete-time linear-quadratic regulator (LQR) problems on the infinite time interval to the solution of a corresponding continuous-time LQR problem as the length of the sampling interval tends toward zero. Convergence of solutions to the operator algebraic Riccati equation and corresponding optimal feedback control gains is guaranteed under appropriate uniform stabilizability and detectability conditions and consistent sampling. Also presented are numerical results involving the optimal LQ control of a heat or diffusion equation, a hereditary or delay differential equation, and a hybrid system of ordinary and partial differential equations describing the transverse vibration of a cantilevered Voigt-Kelvin viscoelastic beam with tip mass.

Rosen, I. G.↗
DOE OSTI2022

Computational Fluid Dynamics and Heat Transfer Modeling of a Dimpled Heat Exchanger

Multiphysics and analytical calculations were conducted for a heat exchanger with passive, natural circulation flow. A glycol/water working fluid convects the heat to a dimpled heat exchanger shell, which subsequently transfers the heat to the soil, which acts as the ultimate heat sink. Because the system is fully-passive, it is not subject to the expenses, maintenance, and mechanical breakdowns associated with moving parts. Density, heat capacity, and thermal conductivity material properties were measured for various soil samples, and subsequently included as input for the soil heat conduction model. The soil model was coupled to a computational fluid dynamics (CFD) heat exchanger model that included the dynamic Smagorinsky large eddy simulation and k- omega turbulence models. The analysis showed that the fluid dynamics and heat transfer models worked properly, albeit at a slow pace. Nevertheless, the coupled CFD/heat conduction simulation ran long enough to determine a key parameter—the amount of heat conducted from the heat exchanger to the ground. This unique performance value, along with experimental data, was used as input for stand-alone, fast-running CFD models, as well as boundaries to obtain solutions to partial differential equations for soil heat conduction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
NASA NTRS1978

Interpolation and Extrapolation of Creep Rupture Data by the Minimum Commitment Method. Part 3: Analysis of Multiheats

The Minimum Commitment Method was applied to two sets of data for which multiple heat information was available. For one alloy, a 304 stainless steel studied in Japan, data on nine well characterized heats were used, while for a proprietary low alloy carbon steel studied in the United Kingdom data were available on seven heats - in many cases to very long rupture times. For this preliminary study no instability factors were used. It was discovered that heat-to-heat variations would be accounted for by introducing heat identifiers in the form A + B log sigma where sigma is the stress and the constants A and B depend only on the heat. With these identifiers all the data could be collapsed onto a single master curve, even though there was considerable scatter among heats. Using these identifiers together with the average behavior of all heats made possible the determination of an accurate constitutive equation for each individual heat. Two basic approaches are discussed for applying the results of the analysis.

Manson, S. S.↗
NASA NTRS1993

Analyzing Martian winds and tracer concentrations using Mars Observer data

During the courses of a day, the Mars Observer spacecraft will acquire globally distributed profiles of the martian atmosphere. It is highly desirable that this data be assembled into synoptic weather maps (complete specifications of the atmospheric pressure, temperature, and winds at a given time), which can in turn be used as starting points in the study of many meteorological phenomena. Unfortunately, the special nature of the Mars Observer data presents several challenges above and beyond the usual difficult problem of data initialization. Mars Observer atmospheric data will consist almost exclusively of asynoptic vertical profiles of temperatures (or radiances) and pressures, whereas winds are generally in balance with horizontal gradients of these quantities (which will not be observed). It will therefore be necessary to resort to dynamical models to analyze the wind fields. As a rule, data assimilation into atmospheric models can result in the generation of spurious gravity waves, so special steps must be taken to suppress these. In addition, the asynoptic nature of the data will require a four-dimensional (space and time) data assimilation scheme. The problem is to find a full set of meteorological fields (winds and temperatures) such that, when marched forward in time in the model, they achieve a best fit (in the weighted least-squares sense) to the data. The proposed solution is to develop a model especially for the Mars Observer data assimilation problem. Gravity waves are filtered from the model by eliminating all divergence terms from the prognostic divergence equation. This leaves a diagnostic gradient wind relation between the rotational wind and the temperature field. The divergent wind is diagnosed as the wind required to maintain the gradient wind balance in the presence of the diabatic heating. The primitive equations of atmospheric dynamics (with three principal dependent variables) are thus reduced to a simpler system with a single prognostic equation for temperature - the variable that will be best observed. (This balance system was apparently first derived by Charney as a first-order Rossby number expansion of the equations of motion). Experience with a full primitive equation model of the martian atmosphere indicates that a further simplification is possible: at least for short-term integrations, the model can be linearized about the zonally symmetric basic state.

Houben, Howard C.↗
NASA NTRS1961

ANALYSIS OF QUARTZ AND TEFLON SHIELDS FOR A PARTICULAR RE-ENTRY MISSION

The transient performance of ablation type heat protection shields is treated herein for the surface of a vehicle returning from outer space to the earth. The vehicle weighs 8640 kg, has a ballistic factor of 500 lb/ft2, re-enters with a speed of 11 km/sec at 120 km altitude, has a lift over drag ratio of -0.5, and is subjected to a maximum deceleration of 7.7 times the gravity constant. By use of well-known equations for the heat transfer and the mass transfer at a heated surface, a numerical calculation method is derived which, for the investigated ablation processes, yields exact transient solutions of the fundamental system of partial differential equations. The method is applied to various quartz shields and to one Teflon shield, which all evaporate so readily under the conditions of the problem at hand that practically no flow of molten shield material exists. The solutions also show comparatively small temperature changes parallel to the surface. The results show that the nose of the vehicle is cooled predominantly by the evaporation of the quartz or the Teflon; the rest of the vehicle's surface is cooled by radiation of the quartz or evaporation of the Teflon. The large mass transfer effects on the nose of the vehicle are detrimental since the resulting low surface temperatures prevent the radiative heat transfer out of the shield, which does not involve any mass loss, from being the desirable governing cooling factor.

ABLATING MATERIAL↗
NASA NTRS2000

Thermal-Flow Code for Modeling Gas Dynamics and Heat Transfer in Space Shuttle Solid Rocket Motor Joints

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.

Wang, Qunzhen↗
NASA NTRS1961

Analysis of Heat Transfer and Pressure Drop for a Gas Flowing Through a set of Multiple Parallel Flat Plates at High Temperatures

Equations were derived representing heat transfer and pressure drop for a gas flowing in the passages of a heater composed of a series of parallel flat plates. The plates generated heat which was transferred to the flowing gas by convection. The relatively high temperature level of this system necessitated the consideration of heat transfer between the plates by radiation. The equations were solved on an IBM 704 computer, and results were obtained for hydrogen as the working fluid for a series of cases with a gas inlet temperature of 200 R, an exit temperature of 5000 0 R, and exit Mach numbers ranging from 0.2 to O.8. The length of the heater composed of the plates ranged from 2 to 4 feet, and the spacing between the plates was varied from 0.003 to 0.01 foot. Most of the results were for a five- plate heater, but results are also given for nine plates to show the effect of increasing the number of plates. The heat generation was assumed to be identical for each plate but was varied along the length of the plates. The axial variation of power used to obtain the results presented is the so-called "2/3-cosine variation." The boundaries surrounding the set of plates, and parallel to it, were assumed adiabatic, so that all the power generated in the plates went into heating the gas. The results are presented in plots of maximum plate and maximum adiabatic wall temperatures as functions of parameters proportional to f(L/D), for the case of both laminar and turbulent flow. Here f is the Fanning friction factor and (L/D) is the length to equivalent diameter ratio of the passages in the heater. The pressure drop through the heater is presented as a function of these same parameters, the exit Mach number, and the pressure at the exit of the heater.

Einstein, Thomas H.↗