Finite transform solution of the temperature of a plate heated by a moving discrete source.
Thin plate temperature distribution from solving conduction equation, considering boundary convection losses and moving discrete heat source
THIN PLATE↗
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.
Particle-in-cell Monte Carlo-collision modeling of non-ideal effects in wave-heated dense microplasmas
A computational model for non-ideal plasma effects during the time evolution of a second-stage laser-heated discharge at high pressures is presented. The model extends a classical one-dimensional particle-in-cell Monte Carlo-collision (PIC-MCC) approach coupled with Maxwell's equations for the laser-heating process of a xenon plasma at 300 K temperature and 10–100 bar pressure. Plasma non-ideality resulting from Coulomb coupling at high plasma densities is manifested as a depression in the effective ionization potential of atoms and enhanced collision cross sections. These non-ideal effects are represented using the Ecker–Kröll model in the context of the PIC-MCC approach. We find that full ionization of the plasma is obtained on the picosecond timescale, starting from the skin layer and quickly expanding throughout the domain through an anomalous extension of the skin depth. More critically, we show that the inclusion of the non-ideal plasma effects results in more rapid ionization compared to an ideal plasma, especially at higher pressures. The ionization delay reduction is of the order of a fraction of a picosecond, corresponding to a 16% decrease at 100 bar. As the amplitude of the wave field is lowered, the ionization rate is lowered, making the plasma non-ideality effects more prominent.
Solmaz, Evrim (ORCID:0000000211830881)↗
An approximate analysis of the diffusing flow in a self-controlled heat pipe.
Constant-density two-dimensional axisymmetric equations are presented for the diffusing flow of a class of self-controlled heat pipes. The analysis is restricted to the vapor space. Condensation of the vapor is related to its mass fraction at the wall by the gas kinetic formula. The Karman-Pohlhausen integral method is applied to obtain approximate solutions. Solutions are presented for a water heat pipe with neon control gas.
Somogyi, D.↗
Planetary-scale circulations in the presence of climatological and wave-induced heating
Interaction between the large-scale circulation and the convective pattern is investigated in a coupled system governed by the linearized primitive equations. Convection is represented in terms of two components of heating: A 'climatological component' is prescribed stochastically to represent convection that is maintained by fixed distributions of land and sea and sea surface temperature (SST). An 'induced component' is defined in terms of the column-integrated moisture flux convergence to represent convection that is produced through feedback with the circulation. Each component describes the envelope organizing mesoscale convective activity. As SST on the equator is increased, induced heating amplifies in the gravest zonal wavenumbers at eastward frequencies, where positive feedback offsets dissipation. Under barotropic stratification, a critical SST of 29.5 C results in positive feedback exactly cancelling dissipation in wavenumber 1 for an eastward phase speed of 6 m/s. Sympathetic interaction between the circulation and the induced heating is the basis for 'frictional wave-Conditional Instability of the Second Kind (CISK)', which is distinguished from classical wave-CISK by rendering the gravest zonal dimensions most unstable. Under baroclinic stratification, the coupled system exhibits similar behavior. The critical SST is only 26.5 C for conditions representative of equinox, but in excess of 30 C for conditions representative of solstice. Having the form of an unsteady Walker circulation, the disturbance produced by frictional wave-CISK compares favorably with the observed life cycle of the Madden-Julian oscillation (MJO). SST above the critical value produces an amplifying disturbance in which enhanced convection coincides with upper-tropospheric westerlies and is positively correlated with temperature and surface convergence. Conversely, SST below the critical value produces a decaying disturbance in which enhanced convection coincides with upper-tropospheric easterlies and is nearly in quadrature with temperature and surface convergence. While sharing essential features with the MJO in the Eastern Hemisphere, frictional wave-CISK does not explain observed behavior in the Western Hemisphere, where the convective signal is largely absent. Comprised of Kelvin structure with the same frequency, observed behavior in the Western Hemisphere can be understood as a propagating response that is excited in and radiates away from the fluctuation of convection in the Eastern Hemisphere.
Salby, Murry L↗
A boundary integral method for an inverse problem in thermal imaging
An inverse problem in thermal imaging involving the recovery of a void in a material from its surface temperature response to external heating is examined. Uniqueness and continuous dependence results for the inverse problem are demonstrated, and a numerical method for its solution is developed. This method is based on an optimization approach, coupled with a boundary integral equation formulation of the forward heat conduction problem. Some convergence results for the method are proved, and several examples are presented using computationally generated data.
Bryan, Kurt↗
An Analysis of Exact and Approximate Equations for the Temperature Distribution in an Insulated Thick Skin Subjected to Aerodynamic Heating
The problem of calculating the temperature distribution in an insulated slab is investigated. Exact and approximate solutions are obtained, and the results are compared to determine the ranges of applicability of the approximations. The approximations are found to be within 5 percent of the exact solution when the ratio of the thermal capacitance of the metal to that of the insulation and the ratio of the conductance of the metal to that of the insulation are sufficiently large. The roots of the characteristic equation of the exact solution are generally applicable to the two-slab heat-transfer problem and are tabulated up to the first nine roots.
Harris, Robert S., Jr.↗
Direct numerical simulation of supercritical oxy-methane mixing layers with CO 2 substituted counterparts
Direct numerical simulations (DNS) of temporally developing, three dimensional, CH 4 /CO 2 , CH 4 /O 2 and CO 2 /O 2 mixing layers, are conducted at a supercritical pressure of 300 atm. To effectively model the supercritical regime, the employed formulation includes the compressible form of the governing equations, the cubic Peng-Robinson equation of state and a generalized formulation for heat and mass flux vectors derived from non-equilibrium thermodynamics and fluctuation theory. A linear inviscid stability analysis is also performed for each case, to determine its most unstable wavelength. Flow visualizations reveal the presence of high density gradient magnitude regions for all three mixing layers, with conditional averages indicating increased presence of heavier fluid species within these regions. No significant departures are observed from perfect gas behavior, with compressibility factors very close to unity for all three mixing cases. Applicability of presumed probability density function methods (PDF) is examined for the three supercritical mixing layers. An a priori analysis is also conducted to investigate various simplifying assumptions employed in modeling various subgrid scale (SGS) flux models. Two additional terms are identified in the large eddy simulations (LES) equations, the gradient of SGS contribution of pressure in the momentum equation and the gradient of SGS contribution of heat flux in energy equation, whose magnitudes are similar and comparable with their respective resolved terms. The performance of the scale similarity model to represent these additional terms is investigated. Lastly, the performance of Smagorinsky, gradient and scale similarity models are also investigated.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Moist wind relationships
Equations describing the temporal and spatial behavior of the kinematic moisture and heat flux are introduced. In these nonlinear equations, the contribution by diabatic processes to the large-scale flux is composed of two parts. One part is associated with a Rayleigh damping term, while the other arises from temporal and spatial changes in the pressure gradient term. The influence of diabatic processes on the large-scale moisture fluxes depends greatly on the degree of balance between forcing and damping terms in the governing equations. The existence of a near balance requires a reduction in the large-scale horizontal geostrophic wind speed. From a scale analysis of the moisture flux equations it is argued that reductions in the large-scale horizontal wind speed, observed within major cumulus cloud systems, help conserve large-scale moisture fluxes. The deviation of the wind from geostrophic conditions is easily estimated. This wind modification induces secondary vertical circulations that contribute to the convergence, creating or supporting long-lived mesoscale flows. In the tropics the wind modification has an antitriptic relationship. These diagnostic findings suggest possible modifications to the wind field in the application of cumulus parameterization, and may be important in diabatic initialization of numerical weather prediction models.
Raymond, William H.↗
Moist wind relationships
The equations describing the temporal and spatial behavior of the kinematic moisture and heat flux are described. In these nonlinear equations, the contribution by diabatic processes to the large-scale flux is composed of two parts. One part is associated with a Rayleigh damping term while the other arises from temporal and spatial changes in the pressure gradient term. It was found that the influence of diabatic processes on large-scale moisture fluxes depends greatly on the degree of balance between forcing and damping terms in the governing kinematic flux equations. The existence of a near balance requires a reduction in the large-scale horizontal geostrophic wind speed. Based on an examination of the moisture flux equations, it is argued that reductions in the large-scale horizontal wind speed observed within major cumulus cloud systems help conserve large-scale moisture fluxes. The deviation of the wind from geostrophic conditions is easily estimated for a near balanced state. This wind modification induces secondary vertical circulations that contribute to convergence, creating or supporting long-lived mesoscale flows. We believe this process to be a major supporter of the mesoscale circulations observed in severe storms and squall lines. In the tropics the wind modification has an antitriptic relationship. These diagnostic findings suggest possible modifications to the wind field in the application of a cumulus parameterization, and may be important in diabatic initialization of numerical weather prediction models.
Raymond, William H.↗
On the Divergence of the Velocity Vector in Real-Gas Flow
A theoretical study was performed addressing the degree of applicability or inapplicability, to a real gas, of the occasionally stated belief that for an ideal gas, incompressibility is synonymous with a zero or very low Mach number. The measure of compressibility used in this study is the magnitude of the divergence of the flow velocity vector [V(bar) (raised dot) u (where u is the flow velocity)]. The study involves a mathematical derivation that begins with the governing equations of flow and involves consideration of equations of state, thermodynamics, and fluxes of heat, mass, and the affected molecular species. The derivation leads to an equation for the volume integral of (V(bar) (raised dot) u)(sup 2) that indicates contributions of several thermodynamic, hydrodynamic, and species-flux effects to compressibility and reveals differences between real and ideal gases. An analysis of the equation leads to the conclusion that for a real gas, incompressibility is not synonymous with zero or very small Mach number. Therefore, it is further concluded, the contributions to compressibility revealed by the derived equation should be taken into account in simulations of real-gas flows.
Bellan, Josette↗
Compressible Laminar Boundary Layer over a Yawed Infinite Cylinder with Heat Transfer and Arbitrary Prandtl Number
The equations are presented for the development of the compressible laminar boundary layer over a yawed infinite cylinder. For compressible flow with a pressure gradient the chordwise and spanwise flows are not independent. Using the Stewartson transformation and a linear viscosity-temperature relation yields a set of three simultaneous ordinary differential equations in a form yielding similar solutions. These equations are solved for stagnation-line flow for surface temperatures from zero to twice the free-stream stagnation temperature and for a wide range of yaw angle and free-stream Mach number. The results indicate that the effect of yaw on the heat-transfer coefficient at the stagnation line depends markedly on the free-stream Mach number. An unusual result of the solutions is that for large yaw angles and stream Mach numbers the chordwise velocity within the boundary layer exceeds the local external chordwise velocity, even for a highly cooled wall.
Reshotko, Eli↗
Development of OPTIMOD - The new Space Shuttle Orbiter/Payload Thermal Integration Model
The development of the Orbiter/Payload Thermal Integration Model (OPTIMOD), a simplified Space Shuttle Orbiter thermal math model to be used for payload thermal integration analysis, is described. The Orbiter lower fuselage math model has been reduced to the minimum level of detail consistent with accuracy requirements. Instead, an 'environment radiation sensor', a simple closed-end cylinder oriented in a fixed position with respect to the Orbiter, is introduced. Equations are presented which relate the cylinder surface heating rates to external fuselage steady-state temperatures. These equations were derived from a regression correlation analysis of steady-state Orbiter thermal solutions for 2372 points on orbit, representing a complete range of vehicle-earth-sun orientations. An exponential equation is developed to relate the steady-state data to time-variant temperatures for the particular transient attitude timeline being simulated. The accuracy of OPTIMOD is shown to be excellent in the payload bay and to compare favorably with previous models in Orbiter elements outside the bay.
Russell, David J.↗
Anisotropic ion heating and parallel O(+) acceleration in regions of rapid E x B convection
A numerical solution to the 20-moment set of transport equations has been found in order to study subauroral ionospheric outflows during periods of enhanced perpendicular ion drifts. The numerical model solves the time-dependent O(+) density, momentum, and both the parallel and perpendicular energy and heat flow equations in the 200-6000 km altitude range. Assuming perpendicular drifts of 3 km/s relative to the neutral atmosphere, we have found that anisotropic heating of O(+) (a result of ion-neutral collisions) leads to a temperature anisotropy, with perpendicular temperatures exceeding 8000 K and parallel temperatures greater than 5000 K (near 200 km altitude). Above approximately 2000 km, transport processes dominate the effects of collisions and wavelike oscillations in O(+) velocity, temperature and heat flux were noted.
Korosmezey, A.↗
NASTRAN solutions of problems described by simultaneous parabolic differential equations
NASTRAN solution techniques are shown for a numerical analysis of a class of coupled vector flow processes described by simultaneous parabolic differential equations. To define one physical problem type where equations of this form arise, the differential equations describing the coupled transfers of heat and mass in mechanical equilibrium with negligible mass average velocity are presented and discussed. Also shown are the equations describing seepage when both electrokinetic and hydrodynamic forces occur. Based on a variational statement of the general problem type, the concepts of scalar transfer elements and parallel element systems are introduced. It is shown that adoptation of these concepts allows the direct use of NASTRAN's existing Laplace type elements for uncoupled flow (the heat transfer elements) for treating multicomponent coupled transfer. Sample problems are included which demonstrate the application of these techniques for both steady-state and transient problems.
Mason, J. B.↗
Boundary formulations for shape sensitivity of temperature dependent conductivity problems
Used in concert with the Kirchhoff transformation, implicit differentiation of the discretized boundary integral equations governing the conduction of heat in solids with temperature dependent thermal conductivity is shown to generate an accurate and economical approach for computation of shape sensitivities. For problems with specified temperature and heat flux boundary conditions, a linear problem results for both the analysis and sensitivity analysis. In problems with either convection or radiation boundary conditions, a nonlinear problem is generated. 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 is presented that have closed form solutions. Exact analytical expressions for the shape sensitivities associated with these problems are developed and these are compared with the sensitivities computed using the boundary element formulation.
Kane, James H.↗
Thermodynamic Properties of Nitrogen Including Liquid and Vapor Phases from 63K to 2000K with Pressures to 10,000 Bar
Tables of thermodynamic properties of nitrogen are presented for the liquid and vapor phases for temperatures from the freezing line to 2000K and pressures to 10,000 bar. The tables include values of density, internal energy, enthalpy, entropy, isochoric heat capacity, isobaric heat capacity velocity of sound, the isotherm derivative, and the isochor derivative. The thermodynamic property tables are based on an equation of state, P=P (p,T), which accurately represents liquid and gaseous nitrogen for the range of pressures and temperatures covered by the tables. Comparisons of property values calculated from the equation of state with measured values for P-p-T, heat capacity, enthalpy, latent heat, and velocity of sound are included to illustrate the agreement between the experimental data and the tables of properties presented here. The coefficients of the equation of state were determined by a weighted least squares fit to selected P-p-T data and, simultaneously, to isochoric heat capacity data determined by corresponding states analysis from oxygen data, and to data which define the phase equilibrium criteria for the saturated liquid and the saturated vapor. The vapor pressure equation, melting curve equation, and an equation to represent the ideal gas heat capacity are also presented. Estimates of the accuracy of the equation of state, the vapor pressure equation, and the ideal gas heat capacity equation are given. The equation of state, derivatives of the equation, and the integral functions for calculating derived thermodynamic properties are included.
Thermodynamic properties↗
Temperatures of individual ion species and heating due to charge exchange in the ionosphere of Venus
The coupled electron and multispecies ion energy equations were solved for daytime conditions in the Venus ionosphere. The heating rates due to charge exchange between hot oxygen atoms and thermal oxygen ions were calculated and incorporated into the energy equations. The combination of the traditional EUV heating and this hot oxygen energy source leads to calculated electron and individual ion temperatures significantly lower than the measured values during solar cycle maximum conditions. Calculations were also carried out for solar cycle minimum conditions, which led to considerably lower temperatures; no data are available which would allow direct comparisons of these results with measurements. In order to obtain calculated temperature values consistent with the observed ones, for solar cycle maximum conditions, topside heat inflows into the ion and electron gases have to be introduced or the thermal conductivity must be reduced by considering the effect of steady and fluctuating magnetic fields, as was done in previous studies. The addition of hot oxygen heating leads to minor increases in the calculated ion temperatures except for the case of reduced thermal conductivities. Separate temperatures were calculated for each ion species for a number of different conditions and in general the differences were found to be relatively small.
Kim, Jhoon↗
Pronghorn Theory Manual
This manual is intended to provide a thorough description of the governing equations; the relations used for material, heat transfer, and fluid flow parameters; and the numerical method utilized in Pronghorn. Detailed derivations of the governing equations are given, beginning from first principles in order to be as explicit as possible about assumptions made along the way and guide any future methods extensions. Additional explanatory or tangential information that is not required for detailed understanding of the governing equations used in Pronghorn is shown in gray boxes in the text. The convention used in this document is to represent all class names (or file names) in typewriter font.
42 ENGINEERING↗