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 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Data correlation and analysis of arc tunnel and wind tunnel tests of RSI joints and gaps, phase 2. Volume 1: Technical report

Heat transfer data measured in gaps representative of those being employed for joints in the space shuttle reusable surface insulation (RSI) thermal protection systems (TPS) were assimilated, analyzed, and correlated. Several types of gap were investigated with emphasis on simple butt joints. Gap widths ranged from 0.0 to 0.76 cm and depths ranged from 1 to 6 cm. Laminar, transitional, and turbulent boundary layer flows over the gap opening were investigated. The angle between gap axis and external flow was varied between 0 and pi/2 radians. The contoured cross section gap performed significantly better than all other wide gaps and slightly better than all other narrow gap geometries. Three dimensional heating variations were observed within gaps in the absence of external flow pressure gradients. Interactions between heating within gaps and heating of adjacent top tile surfaces were observed. Gaps aligned with the flow were observed to promote boundary layer transition. Heat transfer correlation equations were obtained for many of the tests. The TPS thickness requirements with and without gaps were computed for a current shuttle entry trajectory. Experimental data employed in the study are summarized. A description of each test facility, run schedule and test conditions, model descriptive information, and heat flux data are included.

Cristensen, H. E.↗
NASA NTRS2008

A Step Towards CO2-Neutral Aviation

An approximation method for evaluation of the caloric equations used in combustion chemistry simulations is described. The method is applied to generate the equations of specific heat, static enthalpy, and Gibb's free energy for fuel mixtures of interest to gas turbine engine manufacturers. Liquid-phase fuel properties are also derived. The fuels investigated include JP-8, synthetic fuel, and two blends of JP-8 and synthetic fuel. The complete set of fuel property equations for both phases are implemented into a computational fluid dynamics (CFD) flow solver database, and multiphase, reacting flow simulations of a well-tested liquid-fueled combustor are performed. The simulations are a first step in understanding combustion system performance and operational issues when using alternate fuels, at practical engine operating conditions.

Brankovic, Andreja↗
NASA NTRS1986

Theoretical and experimental investigation of the destruction of graphites in a flow of dissociated air

A method is presented for calculating the heating and erosion of blunt bodies made of graphite in a high-enthalpy flow of dissociated air, assuming chemical equilibrium on the surface and taking account of the thermal effects of combustion and sublimation of graphite. The analysis involves the use of a finite difference scheme to solve an equation of unsteady heat conduction. Attention is given to the equilibrium vaporization of C, C2 and C3 molecules. The calculations agree well with experimental data for a wide range of temperatures and stagnation pressures.

Bovina, T. A.↗
NASA NTRS1986

Global diabatic heating during FGGE SOP-1 and SOP-2

With the increase in the observational data provided by FGGE and the use of global circulation models with full physics for the data assimilation, it is now becoming feasible to attempt to estimate globally the atmospheric diabatic heating. The thermodynamic equation in isobaric coordinates and the data generated by the FGGE III-b analysis of the Goddard Laboratory for Atmospheres (GLA) are employed to serve this purpose. The results of the present study generally agree with other previous investigations. However, some important differences are also revealed. (1) The diabatic heating obtained in the tropics in the present study is larger than that obtained elsewhere; (2) the relatively large heating over the mountainous areas shown in other studies does not appear; (3) no significant negative values of diabatic heating are found in the polar regions; and (4) unlike other studies, cooling is noted over parts of Eurasia in the summer.

Chen, Tsing-Chang↗
NASA NTRS2007

A Step Towards CO2-Neutral Aviation

An approximation method for evaluation of the caloric equations used in combustion chemistry simulations is described. The method is applied to generate the equations of specific heat, static enthalpy, and Gibb's free energy for fuel mixtures of interest to gas turbine engine manufacturers. Liquid-phase fuel properties are also derived. The fuels include JP-8, synthetic fuel, and two fuel blends consisting of a mixture of JP-8 and synthetic fuel. The complete set of fuel property equations for both phases are implemented into a computational fluid dynamics (CFD) flow solver database, and multi-phase, reacting flow simulations of a well-tested liquid-fueled combustor are performed. The simulations are a first step in understanding combustion system performance and operational issues when using alternate fuels, at practical engine operating conditions.

Brankovic, Andreja↗
NASA NTRS1954

Heat transfer, diffusion, and evaporation

Although it has long been known that the differential equations of the heat-transfer and diffusion processes are identical, application to technical problems has only recently been made. In 1916 it was shown that the speed of oxidation of the carbon in iron ore depends upon the speed with which the oxygen of the combustion air diffuses through the core of gas surrounding the carbon surface. The identity previously referred to was then used to calculate the amount of oxygen diffusing to the carbon surface on the basis of the heat transfer between the gas stream and the carbon surface. Then in 1921, H. Thoma reversed that procedure; he used diffusion experiments to determine heat-transfer coefficients. Recently Lohrisch has extended this work by experiment. A technically very important application of the identity of heat transfer and diffusion is that of the cooling tower, since in this case both processes occur simultaneously.

GASES - DIFFUSION↗
DOE OSTI2025

Advancing Understanding of Deep Convective Anvil Clouds

In this project we use ground-based DOE ARM observations, satellite observations, and Weather Research and Forecasting (WRF) model simulations to arrive at a more complete understanding (or, new theory) of how atmospheric heating profiles in thunderstorm cores drives how quickly air rises in the storms. The rate at which air motion increases or decreases with altitude directly relates to how much cloud mass is ejected laterally into the adjacent atmosphere, which then influences the development and maintenance of “stratiform raining cloud anvils”. Our work has led to new equations connecting the heating in these thunderstorms to these stratiform clouds, with the connection modulated by the stability structure of the atmosphere (i.e., if more stable, then a more rapid change in thunderstorm heating with altitude results in more stratiform cloud, independent of whether over land or ocean, or how organized the thunderstorms are). Improving knowledge of stratiform clouds is important, since they contribute to substantial accumulated rainfall (important to agriculture and local hydrology) as well as radiation, given their large areal extents and lifetime. The radiation impacts likely impact simulated cloud feedback in all Earth System models (ESMs).

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1951

Blade-to-coolant heat-transfer results and operating data from a natural-convection water-cooled single-stage turbine

Blade-to-coolant heat-transfer data and operating data were obtained with a natural-convection water-cooled turbine over range of turbine speeds and inlet-gas temperatures. The convective coefficients were correlated by the general relation for natural-convection heat transfer. The turbine data were displaced from a theoretical equation for natural convection heat transfer in the turbulent region and from natural-convection data obtained with vertical cylinders and plates; possible disruption of natural convection circulation within the blade coolant passages was thus indicated. Comparison of non dimensional temperature-ratio parameters for the blade leading edge, midchord, and trailing edge indicated that the blade cooling effectiveness is greatest at the midchord and least at the trailing edge.

COOLING - GAS TURBINE SYSTEMS↗
NASA NTRS1974

A model of heat transfer in immersed man

An equation representing man's thermal balance under water is considered. The equation states that the body thermal loading from metabolic heat production and artificial heat input must be offset by respiratory and environmental heat exchange to maintain a constant body temperature. Critical body regions are affected by cold-water thermal stress. A model of the thermoregulatory system may be divided into the physical-controlled system and the dynamic controlling system. The thermal model is simulated by computer programs.

Montgomery, L. D.↗
NASA NTRS1974

Variation of the ratio of specific heats across a detached bow shock

Equations are derived which allow the ratio of specific heats behind the earth's bow shock to be evaluated if several pre-shock parameters (the specific-heat ratio, the Alfvenic Mach number, the sonic Mach number, and the angle between the shock normal at the stagnation point and the magnetic field) and the density jump across the shock are known. Numerical examples show that the dependence of the post-shock ratio on the pre-shock ratio is weak.

Chao, J. K.↗
DOE OSTI2025

Development for Integrated System-Level Analysis Capabilities in SAM for Molten Salt Reactors

In recent years, there has been renewed interest in Molten Salt Reactors (MSRs) for their potential advantages compared to reactors that rely on solid fuel. In response to such interest, many methods and codes have been developed to capture the unique features of MSRs. Among them, the System Analysis Module (SAM) is a modern system analysis tool that provides fast-running, modest-fidelity, whole-plant transient analysis capabilities, essential for fast-turnaround design scoping and engineering analyses of advanced reactor concepts. For liquid-fuel MSRs, the complex physics and chemistry involved in MSR operation—such as reactor kinetics, fluid flow, heat transfer, and salt composition dynamics—pose significant challenges for system-level modeling. Specific modeling capabilities are needed for system-level transient simulation. This paper presents recent advancements in SAM capability enhancements for system-level modeling of MSRs, focusing on improved simulation fidelity, computational efficiency, and multi-physics integration. Key enhancements include the development of species transport, Delayed Neutron Precursor (DNP) drift, modified Point Kinetics Equations (PKE), decay heat modeling, key fission product behavior, salt corrosion, and thermal-hydraulic coupling, as well as code robustness and performance enhancements for MSR applications. The code enhancement allows for better predictive accuracy in safety analysis, transient behavior, and operational optimization, thus supporting the design and licensing of next-generation MSRs. Results from case studies are presented to demonstrate the benefits of these enhancements in accurately capturing key reactor transient behaviors.

Hu, Rui (ORCID:0000000237712920)↗
NASA NTRS2004

Higher-Order Finite Elements for Computing Thermal Radiation

Two variants of the finite-element method have been developed for use in computational simulations of radiative transfers of heat among diffuse gray surfaces. Both variants involve the use of higher-order finite elements, across which temperatures and radiative quantities are assumed to vary according to certain approximations. In this and other applications, higher-order finite elements are used to increase (relative to classical finite elements, which are assumed to be isothermal) the accuracies of final numerical results without having to refine computational meshes excessively and thereby incur excessive computation times. One of the variants is termed the radiation sub-element (RSE) method, which, itself, is subject to a number of variations. This is the simplest and most straightforward approach to representation of spatially variable surface radiation. Any computer code that, heretofore, could model surface-to-surface radiation can incorporate the RSE method without major modifications. In the basic form of the RSE method, each finite element selected for use in computing radiative heat transfer is considered to be a parent element and is divided into sub-elements for the purpose of solving the surface-to-surface radiation-exchange problem. The sub-elements are then treated as classical finite elements; that is, they are assumed to be isothermal, and their view factors and absorbed heat fluxes are calculated accordingly. The heat fluxes absorbed by the sub-elements are then transferred back to the parent element to obtain a radiative heat flux that varies spatially across the parent element. Variants of the RSE method involve the use of polynomials to interpolate and/or extrapolate to approximate spatial variations of physical quantities. The other variant of the finite-element method is termed the integration method (IM). Unlike in the RSE methods, the parent finite elements are not subdivided into smaller elements, and neither isothermality nor other unrealistic physical conditions are assumed. Instead, the equations of radiative heat transfer are integrated numerically over the parent finite elements by use of a computationally efficient Gaussian integration scheme.

Gould, Dana C.↗
NASA NTRS2004

Experiences on p-Version Time-Discontinuous Galerkin's Method for Nonlinear Heat Transfer Analysis and Sensitivity Analysis

The focus of this research is on the development of analysis and sensitivity analysis equations for nonlinear, transient heat transfer problems modeled by p-version, time discontinuous finite element approximation. The resulting matrix equation of the state equation is simply in the form ofA(x)x = c, representing a single step, time marching scheme. The Newton-Raphson's method is used to solve the nonlinear equation. Examples are first provided to demonstrate the accuracy characteristics of the resultant finite element approximation. A direct differentiation approach is then used to compute the thermal sensitivities of a nonlinear heat transfer problem. The report shows that only minimal coding effort is required to enhance the analysis code with the sensitivity analysis capability.

Hou, Gene↗
NASA NTRS1987

Application of turbulence modeling to predict surface heat transfer in stagnation flow region of circular cylinder

A theoretical analysis and numerical calculations for the turbulent flow field and for the effect of free-stream turbulence on the surface heat transfer rate of a stagnation flow are presented. The emphasis is on the modeling of turbulence and its augmentation of surface heat transfer rate. The flow field considered is the region near the forward stagnation point of a circular cylinder in a uniform turbulent mean flow. The free stream is steady and incompressible with a Reynolds number of the order of 10 to the 5th power and turbulence intensity of less than 5 percent. For this analysis, the flow field is divided into three regions: (1) a uniform free-stream region where the turbulence is homogeneous and isotropic; (2) an external viscid flow region where the turbulence is distorted by the variation of the mean flow velocity; and, (3) an anisotropic turbulent boundary layer region over the cylinder surface. The turbulence modeling techniques used are the kappa-epsilon two-equation model in the external flow region and the time-averaged turbulence transport equation in the boundary layer region. The turbulence double correlations, the mean velocity, and the mean temperature within the boundary layer are solved numerically from the transport equations. The surface heat transfer rate is calculated as functions of the free-stream turbulence longitudinal microlength scale, the turbulence intensity, and the Reynolds number.

Wang, Chi R.↗
DOE OSTI2024

Model-based optimization strategies for direct hydrogenation of carbon dioxide to dimethyl ether

Here, this study discusses model-based optimization strategies for CO 2 hydrogenation to dimethyl ether (DME) over aCuZnZr (CZZ) and ferrite (FER) mixed catalyst system in a packed-bed reactor configuration. A two-dimensional axisymmetric, nonisothermal packed-bed reactor model was developed using COMSOL Multiphysics 6.2 software. The model solves two-dimensional (radial and axial) heat and mass transport equations in the packed-bed and integrates intraparticle diffusion and heat transfer in a 1D approach. This powerful feature differs from a traditional porous media approach and takes into account any heat and mass transfer limitations that may exist. Analysis shows that the heat transfer limitations are negligible, but strong internal mass transfer limitations were observed at 10 ≤ WHSV ≤ 90 h -1 on the FER catalyst and at 240 °C. The optimum catalyst composition (i.e., mixing ratio) varies depending on the operating regime. The FER catalyst weight in the mixture can be as low as 5 wt.%, but the ideal composition de-pends on the internal mass transfer limitation and its relationship with the operating regime (i.e., weight hourly space velocity, temperature). A catalyst composition of 80 wt.% CZZ and 20 wt.% FER was suggested; this composition can provide high CO 2 conversion and DME production rates at a wide range of temperatures and flow rates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Using Ensemble Data Assimilation to Estimate Transient Hydrologic Exchange Flow Under Highly Dynamic Flow Conditions

Abstract Quantifying dynamic hydrologic exchange flows (HEFs) within river corridors that experience high‐frequency flow variations caused by dam regulations is important for understanding the biogeochemical processes at the river water and groundwater interfaces. Heat has been widely used as a tracer to infer steady‐state flow velocities through analytical solutions of heat transport defined by the diurnal temperature signals. Under sub‐daily dynamic flow conditions, however, such analytical solutions are not applicable due to the violation of their fundamental assumptions. In this study, we developed a data assimilation‐based approach to estimate the sub‐daily flux under highly dynamic flow conditions using multi‐depth temperature observations at a 5‐min resolution. If the hydraulic gradient is measured, Darcy's law was used to calculate the flux with permeability estimated from temperature responses below the riverbed. Otherwise, flux was estimated directly by assimilating multi‐depth temperature data at 1‐ or 2‐hr time intervals assuming one‐dimensional flow and heat transport governing equation. By comparing estimated fluxes with model‐generated synthetic truth, we demonstrated that both schemes have robust performance in estimating fluxes under highly dynamic flow conditions. This data assimilation‐based flux estimation method was able to capture the vertical sub‐daily fluxes using multi‐depth high‐resolution temperature data alone, even in the presence of multi‐dimensional flow. This approach has been successfully applied to real field temperature data collected at the Hanford site, which experiences highly dynamic HEFs. Our study shows the promise of adopting distributed 1‐D temperature monitoring to capture spatial and temporal exchange dynamics in river corridors at a watershed scale or beyond.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Comparison of hydrodynamic and semi-kinetic treatments for plasma flow along closed field lines

Hydrodynamic and semi-kinetic treatments of plasma flow along closed geomagnetic field lines are compared. The hydrodynamic treatment is based on a simplified 16-moment set of transport equations as the equations for the heat flows are not solved; the heat flows are treated heuristically. The semi-kinetic treatment is based on a particle code. The comparison deals with the distributions of the plasma density, flow velocity, and parallel and perpendicular temperatures as obtained from the two treatments during the various stages of the flow. In the kinetic treatment, the appropriate boundary condition is the prescription of the velocity distribution functions for the particles entering the flux tubes at the ionospheric boundaries; those particles leaving the system are determined by the processes occurring in the flux tube. The prescribed distributions are half-Maxwellian with temperature T(sub 0) and density n(sub 0). In the hydrodynamic model, the prescribed boundary conditions are on density (n(sub 0)), flow velocity (V(sub 0)) and temperature (T(sub 0). It was found that results from the hydrodynamic treatment critically depend on V(sub 0); for early stages of the flow this treatment yields results in good agreement with those from the kinetic treatment, when V(sub 0) = square root of (kT(sub 0)/2 (pi)m), which is the average velocity of particles moving in a given direction for a Maxwellian distribution. During this early stage, the flows developing form the conjugate ionospheres show some distinct transitions. For the first hour or so, the flows are highly supersonic and penetrate deep into the opposite hemispheres, and both hydrodynamics and kinetic treatments yield almost similar features. It is found that during this period heatflow effects are negligibly small. When a flow penetrates deep into the opposite hemisphere, the kinetic treatment predicts reflection and setting up of counterstreaming. In contrast, the hydrodynamic treatment yields a shock in the flow. The reasons for this difference in the two treatments is discussed, showing that in view of the relatively warm ions, the coupling of ion beams and the consequent shock formation in the offequatorial region are not likely due to the enhancements in the beam temperatures. The counterstreaming in the kinetic treatment and the shock in the hydrodynamic treatment first advance upward to the equator and then downward to the ionospheric boundary from where the flow originated. The transit time for this advancement is found to be about 1 hour for the respective models. After 2 hours or so, both models predict that the flows from the ionospheric boundaries are generally subsonic with respect to the local ion-sound speed. At late stages of the flow, when a substantial fraction of ions entering the flux tube begin to return back in the kinetic treatment, the hydrodynamic treatment with the boundary condition V(sub 0) = square root of (kT(sub 0)/2(pi)m) yields an over-refilling, and the choice of V(sub 0) becomes uncertain.

Singh, Nagendra↗