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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Implementation of ADI: Schemes on MIMD parallel computers

In order to simulate the effects of the impingement of hot exhaust jets of High Performance Aircraft on landing surfaces a multi-disciplinary computation coupling flow dynamics to heat conduction in the runway needs to be carried out. Such simulations, which are essentially unsteady, require very large computational power in order to be completed within a reasonable time frame of the order of an hour. Such power can be furnished by the latest generation of massively parallel computers. These remove the bottleneck of ever more congested data paths to one or a few highly specialized central processing units (CPU's) by having many off-the-shelf CPU's work independently on their own data, and exchange information only when needed. During the past year the first phase of this project was completed, in which the optimal strategy for mapping an ADI-algorithm for the three dimensional unsteady heat equation to a MIMD parallel computer was identified. This was done by implementing and comparing three different domain decomposition techniques that define the tasks for the CPU's in the parallel machine. These implementations were done for a Cartesian grid and Dirichlet boundary conditions. The most promising technique was then used to implement the heat equation solver on a general curvilinear grid with a suite of nontrivial boundary conditions. Finally, this technique was also used to implement the Scalar Penta-diagonal (SP) benchmark, which was taken from the NAS Parallel Benchmarks report. All implementations were done in the programming language C on the Intel iPSC/860 computer.

Vanderwijngaart, Rob F.↗
DOE OSTI2024

Increased accuracy and signal-to-noise ratio through recent improvements in infra-red video bolometer fabrication and calibration

The infra-red video bolometer (IRVB) is a diagnostic equipped with an infra-red camera that measures the total radiated power in thousands of lines of sight within a large field of view. Recently validated in MAST-U [Fderici et al., Rev. Sci. Instrum. 94, 033502 (2023)], it offers a high spatial resolution map of the radiated power in the divertor region, where large gradients are expected. The IRVB’s sensing element comprises a thin layer of high Z absorbing material, typically platinum, usually coated with carbon to reduce reflections [Peterson et al., Rev. Sci. Instrum. 79, 10E301 (2008)].Here, the possibility of using a relatively inert material such as titanium, is explored that can be produced in layers up to 1 μm compared to 2.5 μm for Pt and then coat it with Pt of the desired thickness (0.3 μm per side here) and carbon. This leads to a higher temperature signal (about 3 times) and better spatial resolution (about 4 times), resulting in higher accuracy in the measured power [Peterson et al., Rev. Sci. Instrum. 79, 10E301 (2008)]. This assembly is also expected to improve foil uniformity, as the Pt layer is obtained via deposition rather than mechanical processes [Mukai et al., Rev. Sci. Instrum. 87, 2014 (2016)].Given its multi-material composition, measuring the thermal properties of the foil assembly is vital. Various methods using a calibrated laser as a heat source have been developed, analyzing the temperature profile shape [Sano et al., Plasma and Fusion Res. 7, 2405039 (2012)] and [Mukai et al., Rev. Sci. Instrum. 89, 10E114 (2018)] or fitting the calculated laser power for different intensities and frequencies [Fderici et al., Rev. Sci. Instrum. 94, 033502 (2023)]. Here, a simpler approach is presented, which relies on analyzing the separate components of the foil heat equation for a single laser exposure in a given area. This can then be iterated over the entire foil to capture local deviations.

Magnetic Fusion Energy↗
NASA NTRS1976

Oscillations and stability of numerical solutions of the heat conduction equation

The mathematical model and results of numerical solutions are given for the one dimensional problem when the linear equations are written in a rectangular coordinate system. All the computations are easily realizable for two and three dimensional problems when the equations are written in any coordinate system. Explicit and implicit schemes are shown in tabular form for stability and oscillations criteria; the initial temperature distribution is considered uniform.

Kozdoba, L. A.↗
NASA NTRS1993

Analysis, approximation, and computation of a coupled solid/fluid temperature control problem

An optimization problem is formulated motivated by the desire to remove temperature peaks, i.e., 'hot spots', along the bounding surfaces of containers of fluid flows. The heat equation of the solid container is coupled to the energy equations for the fluid. Heat sources can be located in the solid body, the fluid, or both. Control is effected by adjustments to the temperature of the fluid at the inflow boundary. Both mathematical analyses and computational experiments are given.

Gunzburger, Max D.↗
NASA NTRS1973

Temperature

The problems in human comfort in heat stress are emphasized, with less emphasis placed upon cold exposure problems. Physiological parameters related to human thermal interactions are discussed, as well as data concerning thermal protective clothing. The energy balance equation, heat transfer equation, thermal comfort, heat stress, and cold stress are also considered. A two node model of human temperature regulation in FORTRAN is appended.

Berenson, P. J.↗
NASA NTRS1992

An engineering aerodynamic heating method for hypersonic flow

A capability to calculate surface heating rates has been incorporated in an approximate three-dimensional inviscid technique. Surface streamlines are calculated from the inviscid solution, and the axisymmetric analog is then used along with a set of approximate convective-heating equations to compute the surface heat transfer. The method is applied to blunted axisymmetric and three-dimensional ellipsoidal cones at angle of attack for the laminar flow of a perfect gas. The method is also applicable to turbulent and equilibrium-air conditions. The present technique predicts surface heating rates that compare favorably with experimental (ground-test and flight) data and numerical solutions of the Navier-Stokes (NS) and viscous shock-layer (VSL) equations. The new technique represents a significant improvement over current engineering aerothermal methods with only a modest increase in computational effort.

Riley, Christopher J.↗
NASA NTRS1992

An Engineering Aerodynamic Heating Method for Hypersonic Flow

A capability to calculate surface heating rates has been incorporated in an approximate three-dimensional inviscid technique. Surface streamlines are calculated from the inviscid solution, and the axisymmetric analog is then used along with a set of approximate convective-heating equations to compute the surface heat transfer. The method is applied to blunted axisymmetric and three-dimensional ellipsoidal cones at angle of attack for the laminar flow of a perfect gas. The method is also applicable to turbulent and equilibrium-air conditions. The present technique predicts surface heating rates that compare favorably with experimental (ground-test and flight) data and numerical solutions of the Navier-Stokes (NS) and viscous shock-layer (VSL) equations. The new technique represents a significant improvement over current engineering aerothermal methods with only a modest increase in computational effort.

Riley, Christopher J.↗
NASA NTRS1994

An axisymmetric analog two-layer convective heating procedure with application to the evaluation of Space Shuttle Orbiter wing leading edge and windward surface heating

A numerical procedure for predicting the convective heating rate of hypersonic reentry vehicles is described. The procedure, which is based on the axisymmetric analog, consists of obtaining the three-dimensional inviscid flowfield solution; then the surface streamlines and metrics are calculated using the inviscid velocity components on the surface; finally, an axisymmetric boundary layer code or approximate convective heating equations are used to evaluate heating rates. This approach yields heating predictions to general three-dimensional body shapes. The procedure has been applied to the prediction of the wing leading edge heating to the Space Shuttle Orbiter. The numerical results are compared with the results of heat transfer testing (OH66) of an 0.025 scale model of the Space Shuttle Orbiter configuration in the Calspan Hypersonic Shock Tunnel (HST) at Mach 10 and angles of attack of 30 and 40 degrees. Comparisons with STS-5 flight data at Mach 9.15 and angle of attack of 37.4 degrees and STS-2 flight data at Mach 12.86 and angle of attack of 39.7 degrees are also given.

Wang, K. C.↗
DOE OSTI2026

Explicit Monotone Stable Super-Time-stepping Methods for Finite Time Singularities

We explore a novel way to numerically resolve the scaling behavior of finite-time singularities in solutions of nonlinear parabolic PDEs. The Runge–Kutta–Legendre (RKL) and Runge–Kutta–Gegenbauer (RKG) super-time-stepping methods were originally developed for nonlinear complex physics problems with diffusion. These are multistage single step second-order, forward-in-time methods with no implicit solves. The advantage is that the time-step size for stability scales with stage number 𝑠 as $\mathcal{O}$⁡(𝑠 2 ). Many interesting nonlinear PDEs have finite-time singularities, and the presence of diffusion often limits one to using implicit or semi-implicit time-step methods for stability constraints. Finite-time singularities are particularly challenging due to the large range of scales that one desires to resolve, often with adaptive spatial grids and adaptive time steps. Here, in this study, we show two examples of nonlinear PDEs for which the self-similar singularity structure has time and space scales that are resolvable using the RKL and RKG methods, without forcing even smaller time steps. Compared to commonly used implicit numerical methods, we achieve a significantly smaller run time while maintaining comparable accuracy. We also prove numerical monotonicity for both the RKL and RKG methods under their linear stability conditions for the constant coefficient heat equation, in the case of infinite domain and periodic boundary condition, leading to a theoretical guarantee of the superiority of the RKL and RKG methods over traditional super-time-stepping methods, such as the Runge-Kutta-Chebyshev and the orthogonal Runge-Kutta-Chebyshev methods. Code can be found at https://github.com/ZT220501/SRK-Singularity.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Migration of the separation point on a deforming cylinder

An iterative scheme of solving the Navier-Stokes equations for unsteady two-dimensional flow about a deforming and translating cylinder is given. In the kth iteration, the convected vorticity appears as the source term in an equation of transient diffusion of vorticity. A novel integral transform is used to reduce this transient vorticity equation to a k-dimensional heat equation. The bounded solution of this equation is obtained with a general method of superposition for problems involving a moving boundary. An equation describing the migration of the separation point on a deforming cylinder in unsteady cross flows is derived from the analytically obtained velocity field. The radial expansion of the cylinder surface is shown to hasten the separation time and to increase the separation angle. The results imply that in the steady flow past a body that is moving forward and sinking at constant rates, the locus of points at which the azimuthal component of skin friction changes sign originates on the leeward ray at a point downstream of the front tip.

Lin, S. P.↗
NASA NTRS1977

Combustion response modeling for composite solid propellants

A computerized mathematical model of the combustion response function of composite solid propellants was developed with particular attention to the contributions of the solid phase heterogeneity. The one-dimensional model treats the solid phase as alternating layers of ammonium perchlorate and binder, with an exothermic melt layer at the surface. Solution of the Fourier heat equation in the solid provides temperature and heat flux distributions with space and time. The problem is solved by conserving the heat flux at the surface from that produced by a suitable model of the gas phase. An approximation of the BDP flame model is utilized to represent the gas phase. By the use of several reasonable assumptions, it is found that a significant portion of the problem can be solved in closed form. A method is presented by which the model can be applied to tetramodal particle size distributions. A computerized steady-state version of the model was completed, which served to validate the various approximations and lay a foundation for the combustion response modeling. The combustion response modeling was completed in a form which does not require an iterative solution, and some preliminary results were acquired.

Source record↗
NASA NTRS1952

Development of a laminar boundary layer behind a suction point

A theoretical investigation is made of the development of a laminar boundary layer behind a suction slot that is assumed to cut off part of the boundary layer without exerting any sink effect. The development, which is approximate, is based on the heat conduction equation. The heat conduction equation enters the analysis through a linearization of the Prandtl-Mises form of the boundary-layer equation.

EQUATIONS - BOUNDARY LAYER↗
NASA NTRS1988

Bistable direction switching in an off-axis pumped continuous wave ruby laser

A report is presented of the observation of hysteretic bistable direction switching in a single-mode CW ruby laser system. This effect is only observed when the pump beam which is focused into the ruby rod is misaligned with respect to the rod end faces. At low pump powers, the ruby lases in a mode nearly collinear with the pump axis. At a higher pump power the ruby switches to a mode that is collinear with the rod end faces and preserves the original polarization. The effect is large enough to switch the beam by an angle equal to twice the diffraction angle. The observations show that under steady-state pumping, a CW ruby laser can exhibit bistable operation in its output direction and power. A calculation using the heat equation with two concentric cylinders with one as a heat source (pump laser) and the outer wall of the other held at 77 K, gives an increase in core temperature of about 0.01 K. Therefore, the increase in temperature is not large enough to change the index of refraction to account for such large macroscopic effects.

Afzal, R. Sohrab↗