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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

A model of the tropical Pacific sea surface temperature climatology

A model for the climatological mean sea surface temperature (SST) of the tropical Pacific Ocean is developed. The upper ocean response is computed using a time dependent, linear, reduced gravity model, with the addition of a constant depth frictional surface layer. The full three-dimensional temperature equation and a surface heat flux parameterization that requires specification of only wind speed and total cloud cover are used to evaluate the SST. Specification of atmospheric parameters, such as air temperature and humidity, over which the ocean has direct influence, is avoided. The model simulates the major features of the observed tropical Pacific SST. The seasonal evolution of these features is generally captured by the model. Analysis of the results demonstrates the control the ocean has over the surface heat flux from ocean to atmosphere and the crucial role that dynamics play in determining the mean SST in the equatorial Pacific. The sensitivity of the model to perturbations in the surface heat flux, cloud cover specification, diffusivity, and mixed layer depth is discussed.

Seager, Richard↗
NASA NTRS1991

Numerical simulation of icing, deicing, and shedding

An algorithm has been developed to numerically model the concurrent phenomena of two-dimensional transient heat transfer, ice accretion, ice shedding and ice trajectory which arise from the use of electrothermal pad. The Alternating Direction Implicit method is used to simultaneously solve the heat transfer and accretion equations occurring in the multilayered body covered with ice. In order to model the phase change between ice and water, a technique was used which assumes a phase for each node. This allows the equations to be linearized such that a direct solution is possible. This technique requires an iterative procedure to find the correct phase at each node. The computer program developed to find this solution has been integrated with the NASA-Lewis flow/trajectory code LEWICE.

Wright, W. B.↗
NASA NTRS1991

Two-dimensional hydrodynamic calculations of the nonlinear development of the Goldreich-Schubert-Fricke instability in a rotating annulus

Two-dimensional nonlinear hydrodynamic calculations are presented which may help assess the effectiveness of the instability in transporting angular momentum in the equatorial zones of stars and planets which are stably stratified with respect to convection. The calculations were made by numerically integrating the 2D axisymmetric Navier-Stokes equations, including viscosity and heat conduction. The instability was followed into the nonlinear regime. The maximum rms velocity amplitude was found to correlate well with the product of the linear growth rate and radial length scale of the instability, consistent with the idea that the instability grows to an amplitude such that an eddy turnover time becomes equal to the growth time defined by the inverse of the growth rate. The time scale for angular momentum to be redistributed to a state of marginal stability was consistent with this picture. The results suggest that in physical situations a state of marginal stability will be maintained, since departures from such a state will be rapidly corrected.

Korycansky, D. G.↗
NASA NTRS1990

Navier-Stokes computations for aerodynamic configurations at high angles of attack

Navier-Stokes computations at high angles of attack over aerodynamic configurations are presented using an implicit finite-volume algorithm. The spatial differencing is upwind-biased for the convective and pressure terms and central for the shear stress and heat transfer terms. The equations are relaxed to steady state with a spatially factored implicit algorithm. In order to treat general geometries, a multiblock patched-grid framework is implemented. Applications and detailed comparisons with experimental data are made for two simple but representative geometric shapes: (1) a highly swept delta wing and (2) a prolate spheroid of 6:1 length-to-diameter. Recent extensions of the algorithm to compute the flow over an F-18 forebody-strake configuration are shown, including comparisons with wind tunnel and flight test results. Comparisons across the range of Reynolds number for all cases indicate that either low Reynolds number or high Reynolds number flows are simulated well, but that flows at intermediate Reynolds number cannot be simulated accurately without a detailed knowledge of the transitional zone between laminar and turbulent flow.

Thomas, J. L.↗
NASA NTRS1992

Numerical analysis of a thermal deicer

An algorithm has been developed to numerically model the concurrent phenomena of two-dimensional transient heat transfer, ice accretion and ice shedding which arise from the use of an electrothermal pad. The Alternating Direction Implicit method is used to simultaneously solve the heat transfer and accretion equations occurring in a multilayered body covered with ice. In order to model the phase change between ice and water, a technique was used which assumes a phase for each node. This allows the equations to be linearized such that a direct solution is possible. This technique requires an iterative procedure to find the correct phase at each node. The computer program developed to find this solution has been integrated with the NASA/Lewis flow/trajectory code LEWICE.

Wright, W. B.↗
NASA NTRS1993

Nongray gas analyses for reflecting walls utilizing a flux technique

A flux formulation for a planar slab of molecular gas radiation bounded by diffuse reflecting walls is developed. While this formulation is limited to the planar geometry, it is useful for studying approximations necessary in modeling nongray radiative heat transfer. The governing equations are derived by considering the history of multiple reflections between the walls. Accurate solutions are obtained by explicitly accounting for a finite number of reflections and approximating the spectral effects of the remaining reflections. Four approximate methods are presented and compared using a single absorption band of H2O. All four methods reduce to an identical zeroth-order formulation, which accounts for all reflections approximately but does handle nonreflected radiation correctly. A single absorption band of CO2 is also considered using the best-behaved approximation for higher orders. A zeroth-order formulation is sufficient to predict the radiative transfer accurately for many cases considered. For highly reflecting walls, higher order solutions are necessary for better accuracy. Including all the important bands of H2O, the radiative source distributions are also obtained for two different temperature and concentration profiles.

Menart, J. A.↗
NASA NTRS1972

An analysis of the radiation field beneath a bank of tubular quartz lamps

Equations governing the incident heat flux distribution beneath a lamp-reflector system were developed. Analysis of a particular radiant heating facility showed good agreement between theory and experiment when a lamp power loss correction was used. In addition, the theory was employed to estimate thermal disruption in the radiation field caused by a protruding probe.

Ash, Robert L.↗
NASA NTRS1994

An analysis code for the Rapid Engineering Estimation of Momentum and Energy Losses (REMEL)

Nonideal behavior has traditionally been modeled by defining efficiency (a comparison between actual and isentropic processes), and subsequent specification by empirical or heuristic methods. With the increasing complexity of aeropropulsion system designs, the reliability of these more traditional methods is uncertain. Computational fluid dynamics (CFD) and experimental methods can provide this information but are expensive in terms of human resources, cost, and time. This report discusses an alternative to empirical and CFD methods by applying classical analytical techniques and a simplified flow model to provide rapid engineering estimates of these losses based on steady, quasi-one-dimensional governing equations including viscous and heat transfer terms (estimated by Reynold's analogy). A preliminary verification of REMEL has been compared with full Navier-Stokes (FNS) and CFD boundary layer computations for several high-speed inlet and forebody designs. Current methods compare quite well with more complex method results and solutions compare very well with simple degenerate and asymptotic results such as Fanno flow, isentropic variable area flow, and a newly developed, combined variable area duct with friction flow solution. These solution comparisons may offer an alternative to transitional and CFD-intense methods for the rapid estimation of viscous and heat transfer losses in aeropropulsion systems.

Dechant, Lawrence J.↗
NASA NTRS1995

Simulated dynamic response of a multi-stage compressor with variable molecular weight flow medium

A mathematical model of a multi-stage compressor with variable molecular weight flow medium is derived. The modeled system consists of a five stage, six cylinder, double acting, piston type compressor. Each stage is followed by a water cooled heat exchanger which serves to transfer the heat of compression from the gas. A high molecular weight gas (CFC-12) mixed with air in varying proportions is introduced to the suction of the compressor. Condensation of the heavy gas may occur in the upper stage heat exchangers. The state equations for the system are integrated using the Advanced Continuous Simulation Language (ACSL) for determining the system's dynamic and steady state characteristics under varying operating conditions.

Babcock, Dale A.↗
NASA NTRS1999

Dynamics and Structure of Dusty Reacting Flows: Inert Particles in Strained, Laminar, Premixed Flames

A detailed numerical study was conducted on the dynamics and thermal response of inert, spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport, For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.

Egolfopoulos, Fokion N.↗
NASA NTRS1999

Dynamics and Structure of Dusty Reacting Flows: Inert Particles in Strained, Laminar, Premixed Flames

A detailed numerical study was conducted on the dynamics and thermal response of inert spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport. For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.

Egolfopoulos, Fokion N.↗
NASA NTRS1998

Multi-Component Diffusion with Application To Computational Aerothermodynamics

The accuracy and complexity of solving multicomponent gaseous diffusion using the detailed multicomponent equations, the Stefan-Maxwell equations, and two commonly used approximate equations have been examined in a two part study. Part I examined the equations in a basic study with specified inputs in which the results are applicable for many applications. Part II addressed the application of the equations in the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) computational code for high-speed entries in Earth's atmosphere. The results showed that the presented iterative scheme for solving the Stefan-Maxwell equations is an accurate and effective method as compared with solutions of the detailed equations. In general, good accuracy with the approximate equations cannot be guaranteed for a species or all species in a multi-component mixture. 'Corrected' forms of the approximate equations that ensured the diffusion mass fluxes sum to zero, as required, were more accurate than the uncorrected forms. Good accuracy, as compared with the Stefan- Maxwell results, were obtained with the 'corrected' approximate equations in defining the heating rates for the three Earth entries considered in Part II.

Sutton, Kenneth↗
NASA NTRS2007

Gravity Scaling of a Power Reactor Water Shield

A similarity analysis on a water-based reactor shield examined the effect of gravity on free convection between a reactor shield inner and outer vessel boundaries. Two approaches established similarity between operation on the Earth and the Moon: 1) direct scaling of Rayleigh number equating gravity-surface heat flux products, 2) temperature difference between the wall and thermal boundary layer held constant. Nusselt number for natural convection (laminar and turbulent) is assumed of form Nu = CRa(sup n).

Reid, Robert S.↗
NASA NTRS2012

Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure

The simulation of high-pressure turbulent flows, where the pressure, p, is larger than the critical value, p(sub c), for the species under consideration, is relevant to a wide array of propulsion systems, e.g. gas turbine, diesel, and liquid rocket engines. Most turbulence models, however, have been developed for atmospheric-p turbulent flows. The difference between atmospheric-p and supercritical-p turbulence is that, in the former situation, the coupling between dynamics and thermodynamics is moderate to negligible, but for the latter it is very significant, and can dominate the flow characteristics. The reason for this stems from the mathematical form of the equation of state (EOS), which is the perfect-gas EOS in the former case, and the real-gas EOS in the latter case. For flows at supercritical pressure, p, the large eddy simulation (LES) equations consist of the differential conservation equations coupled with a real-gas EOS. The equations use transport properties that depend on the thermodynamic variables. Compared to previous LES models, the differential equations contain not only the subgrid scale (SGS) fluxes, but also new SGS terms, each denoted as a correction. These additional terms, typically assumed null for atmospheric pressure flows, stem from filtering the differential governing equations, and represent differences between a filtered term and the same term computed as a function of the filtered flow field. In particular, the energy equation contains a heat-flux correction (q-correction) that is the difference between the filtered divergence of the heat flux and the divergence of the heat flux computed as a function of the filtered flow field. In a previous study, there was only partial success in modeling the q-correction term, but in this innovation, success has been achieved by using a different modeling approach. This analysis, based on a temporal mixing layer Direct Numerical Simulation database, shows that the focus in modeling the q-correction should be on reconstructing the primitive variable gradients rather than their coefficients, and proposes the approximate deconvolution model (ADM) as an effective means of flow field reconstruction for LES heat flux calculation. Further, results for a study conducted for temporal mixing layers initially containing oxygen in the lower stream, and hydrogen or helium in the upper stream, show that, for any LES, including SGS-flux models (constant-coefficient Gradient or Scale-Similarity models, dynamic-coefficient Smagorinsky/Yoshizawa or mixed Smagorinsky/Yoshizawa/Gradient models), the inclusion of the q-correction in the LES leads to the theoretical maximum reduction of the SGS heat-flux difference. The remaining error in modeling this new subgrid term is thus irreducible.

Bellan, Josette↗
NASA NTRS2017

Stabilization of Hypersonic Boundary Layers by Linear and Nonlinear Optimal Perturbations

The effect of stationary, finite-amplitude, linear and nonlinear optimal perturbations on the modal disturbance growth in a Mach 6 axisymmetric flow over a 7 deg. half-angle cone with 0:126 mm nose radius and 0:305 m length is investigated. The freestream parameters (M = 6, Re(exp 1) = 18 x 10(exp. 6) /m) are selected to match the flow conditions of a previous experiment in the VKI H3 hypersonic tunnel. Plane-marching parabolized stability equations are used in conjunction with a partial-differential equation based planar eigenvalue analysis to characterize the boundary layer instability in the presence of azimuthally periodic streaks. The streaks are observed to stabilize nominally planar Mack mode instabilities, although oblique Mack mode and first-mode disturbances are destabilized. Experimentally measured transition onset in the absence of any streaks correlates with an amplification factor of N = 6 for the planar Mack modes. For high enough streak amplitudes, the transition threshold of N = 6 is not reached by the Mack mode instabilities within the length of the cone; however, subharmonic first-mode instabilities, which are destabilized by the presence of the streaks, do reach N = 6 near the end of the cone. The highest stabilization is observed at streak amplitudes of approximately 20 percent of the freestream velocity. Because the use of initial disturbance profiles based on linear optimal growth theory may yield suboptimal control in the context of nonlinear streaks, the computational predictions are extended to nonlinear optimal growth theory. Results show that by using nonlinearly optimal perturbation leads to slightly enhanced stabilization of plane Mack mode disturbances as well as reduced destabilization of subharmonic first-mode disturbances.

Incompressible boundary layer↗
NASA NTRS2019

Enhanced Habitability on High Obliquity Bodies near the Outer Edge of the Habitable Zone of Sun-Like Stars

High obliquity planets represent potentially extreme limits of terrestrial climate, as they exhibit large seasonality, a reversed annual-mean pole-to-equator gradient of stellar heating, and novel cryospheres. A suite of 3D global climate model simulations is performed for low and high obliquity planets with various stellar fluxes, CO2 concentrations, and initial conditions to explore the propensity for high obliquity climates to undergo global glaciation. We also simulate planets with thick CO2 or H2 atmospheres, such as those expected to develop near or beyond the outer edge of the habitable zone. We show that high obliquity planets are hotter than their low obliquity counterparts due to ice-albedo feedbacks for cold climates, and water vapor in warm climates. We suggest that the water vapor greenhouse trapping is greater on high obliquity bodies for a given global-mean temperature due to the different dynamical regimes that occur between the two states. While equatorial ice belts are stable at high obliquity in some climate regimes, it is substantially harder to achieve global glaciation than for a low obliquity planet. Temperate polar conditions can be present at high obliquity at forcings for which low obliquity planets would be in a hard snowball state. Furthermore, open ocean can persist even in the winter hemisphere and when global-mean temperatures are well below freezing. However, the influence of obliquity diminishes for dense atmospheres, in agreement with calculations from 1D energy balance models.

Planet and satellite atmospheres↗
NASA NTRS2019

Enhanced Habitability on High Obliquity Bodies Near the Outer Edge of the Habitable Zone of Sun-Like Stars

High obliquity planets represent potentially extreme limits of terrestrial climate, as they exhibit large seasonality, a reversed annual-mean pole-to-equator gradient of stellar heating, and novel cryospheres. A suite of 3-D global climate model simulations with a dynamic ocean is performed with Earthlike atmospheres for low and high obliquity planets with various stellar fluxes, CO2 concentrations, and initial conditions to explore the propensity for high obliquity climates approaching the outer edge of the Habitable Zone to undergo global glaciation. We also simulate planets with thick CO2 or H2 atmospheres, such as those expected to develop near or beyond the outer edge of the Habitable Zone. We show that high obliquity planets are hotter than their low obliquity counterparts due to ice-albedo feedbacks for cold climates, and water vapor in warm climates. We suggest that the water vapor greenhouse trapping is greater on high obliquity bodies due to the different dynamical regimes that occur between the two states. While equatorial ice-belts are stable at high obliquity in some climate regimes, it is harder to achieve global glaciation than for a low obliquity planet. Temperate polar conditions can be present at high obliquity at forcings for which low obliquity planets would be in a hard snowball state. We suggest the conditions on high obliquity planets are likely to be more favorable for a robust biosphere to develop approaching the outer edge of the HZ. However, the influence of obliquity diminishes for dense atmospheres, in agreement with calculations from 1-D Energy Balance Models.

planets and satellites: atmospheres↗