Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite temperature”

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 919 records · Page 51

Marangoni instability in a liquid layer with two free surfaces

The onset of the Marangoni instability in a liquid layer with two free nearly insulating surfaces heated from below is studied. Linear stability analysis yields a condition for the emergence of a longwave or a finite wavelength instability from the quiescent equilibrium state. Using the method of asymptotic expansions, a weakly nonlinear evolution equation describing the spatiotemporal behavior of the velocity and temperature fields at the onset of the longwave instability is derived. The latter is given by delta(M) = 24, delta(M) being the difference between the upper and the lower Marangoni numbers. It is shown that in some parametric range one convective cell forms across the layer, while in other parametric domains two convective cells emerge between the two free surfaces.

Deissler, Robert J.↗

Models of thermal/chemical boundary layer convection: Potential application to Venus

The upper boundary layer of Venus is comprised of at least two distinct chemical components, mantle and crust. Fluid dynamical models of convection within Venus' mantle were primarily of the thermal boundary layer type. Models assessing the ability of convective mantle flows to deform the crust were undertaken, but models exploring the effects of a variable thickness crust on mantle convection were largely lacking. A Venusian crust of variable thickness could couple back into, and alter, the mantle flow patterns that helped create it, leading to deformation mechanisms not predicted by purely thermal boundary layer convection models. This possibility is explored through a finite element model of thermal/chemical boundary layer convection. Model results suggest that a crust of variable thickness can serve as a mantle flow driver by perturbing lateral temperature gradients in the upper mantle. Resulting mantle flow is driven by the combination of free convective and nonuniform crustal distribution. This combination can lead to a flow instability manifest in the occurrence of episodic mantle lithosphere subduction initiated at the periphery of a crustal plateau. The ability of a light, near surface, chemical layer to potentially alter mantle flow patterns suggest that mantle convection and the creation and/or deformation of such a chemical layer may be highly nonseparable problems on time scales of 10(exp 8) years.

Lenardic, A.↗

Comparison of numerical simulation and experimental data for steam-in-place sterilization

A complex problem involving convective flow of a binary mixture containing a condensable vapor and noncondensable gas in a partially enclosed chamber was modelled and results compared to transient experimental values. The finite element model successfully predicted transport processes in dead-ended tubes with inside diameters of 0.4 to 1.0 cm. When buoyancy driven convective flow was dominant, temperature and mixture compositions agreed with experimental data. Data from 0.4 cm tubes indicate diffusion to be the primary air removal method in small diameter tubes and the diffusivity value in the model to be too large.

Young, Jack H.↗

Computational Fluid Dynamic (CFD) analysis of axisymmetric plume and base flow of film/dump cooled rocket nozzle

Film/dump cooling a rocket nozzle with fuel rich gas, as in the National Launch System (NLS) Space Transportation Main Engine (STME), adds potential complexities for integrating the engine with the vehicle. The chief concern is that once the film coolant is exhausted from the nozzle, conditions may exist during flight for the fuel-rich film gases to be recirculated to the vehicle base region. The result could be significantly higher base temperatures than would be expected from a regeneratively cooled nozzle. CFD analyses were conduced to augment classical scaling techniques for vehicle base environments. The FDNS code with finite rate chemistry was used to simulate a single, axisymmetric STME plume and the NLS base area. Parallel calculations were made of the Saturn V S-1 C/F1 plume base area flows. The objective was to characterize the plume/freestream shear layer for both vehicles as inputs for scaling the S-C/F1 flight data to NLS/STME conditions. The code was validated on high speed flows with relevant physics. This paper contains the calculations for the NLS/STME plume for the baseline nozzle and a modified nozzle. The modified nozzle was intended to reduce the fuel available for recirculation to the vehicle base region. Plumes for both nozzles were calculated at 10kFT and 50kFT.

Tucker, P. K.↗

Efficient Use of Distributed Systems for Scientific Applications

Distributed computing has been regarded as the future of high performance computing. Nationwide high speed networks such as vBNS are becoming widely available to interconnect high-speed computers, virtual environments, scientific instruments and large data sets. One of the major issues to be addressed with distributed systems is the development of computational tools that facilitate the efficient execution of parallel applications on such systems. These tools must exploit the heterogeneous resources (networks and compute nodes) in distributed systems. This paper presents a tool, called PART, which addresses this issue for mesh partitioning. PART takes advantage of the following heterogeneous system features: (1) processor speed; (2) number of processors; (3) local network performance; and (4) wide area network performance. Further, different finite element applications under consideration may have different computational complexities, different communication patterns, and different element types, which also must be taken into consideration when partitioning. PART uses parallel simulated annealing to partition the domain, taking into consideration network and processor heterogeneity. The results of using PART for an explicit finite element application executing on two IBM SPs (located at Argonne National Laboratory and the San Diego Supercomputer Center) indicate an increase in efficiency by up to 36% as compared to METIS, a widely used mesh partitioning tool. The input to METIS was modified to take into consideration heterogeneous processor performance; METIS does not take into consideration heterogeneous networks. The execution times for these applications were reduced by up to 30% as compared to METIS. These results are given in Figure 1 for four irregular meshes with number of elements ranging from 30,269 elements for the Barth5 mesh to 11,451 elements for the Barth4 mesh. Future work with PART entails using the tool with an integrated application requiring distributed systems. In particular this application, illustrated in the document entails an integration of finite element and fluid dynamic simulations to address the cooling of turbine blades of a gas turbine engine design. It is not uncommon to encounter high-temperature, film-cooled turbine airfoils with 1,000,000s of degrees of freedom. This results because of the complexity of the various components of the airfoils, requiring fine-grain meshing for accuracy. Additional information is contained in the original.

Taylor, Valerie↗

Residual Stresses in Thermal Barrier Coatings for a Cu-8Cr-4Nb Substrate System

Analytical calculations were conducted to determine the thermal stresses developed in a coated copper-based alloy, Cu-8%(at.%)Cr-4%Nb (designated as GRCop-84), after plasma spraying and during heat-up in a simulated rocket engine environment. Finite element analyses were conducted for two coating systems consisting of a metallic top coat, a pure copper bond coat and the GRCop-84. The through thickness temperature variations were determined as a function of coating thickness for two metallic coatings, a Ni-17%(wt%)Cr-6%Al-0.5%Y alloy and a Ni-50%(at.%)Al alloy. The residual stresses after low-pressure plasma spraying of the NiCrAlY and NiAl coatings on GRCop-84 substrate were also evaluated. These analyses took into consideration a 50.8 mm copper bond coat and the effects of an interface coating roughness. The through the thickness thermal stresses developed in coated liners were also calculated after 15 minutes of exposure in a rocket environment with and without an interfacial roughness.

Ghosn, Louis J.↗

Analysis and Sizing for Transient Thermal Heating of Insulated Aerospace Vehicle Structures

An analytical solution was derived for the transient response of an insulated structure subjected to a simplified heat pulse. The solution is solely a function of two nondimensional parameters. Simpler functions of these two parameters were developed to approximate the maximum structural temperature over a wide range of parameter values. Techniques were developed to choose constant, effective thermal properties to represent the relevant temperature and pressure-dependent properties for the insulator and structure. A technique was also developed to map a time-varying surface temperature history to an equivalent square heat pulse. Equations were also developed for the minimum mass required to maintain the inner, unheated surface below a specified temperature. In the course of the derivation, two figures of merit were identified. Required insulation masses calculated using the approximate equation were shown to typically agree with finite element results within 10%-20% over the relevant range of parameters studied.

Blosser, Max L.↗

Mixed Mode Fracture of Plasma Sprayed Thermal Barrier Coatings: Effects of Anisotropy and Heterogeneity

The combined mode I-mode II fracture behavior of anisotropic ZrO2-8wt%Y2O3 thermal barrier coatings was determined in asymmetric flexure loading at both ambient and elevated temperatures. A fracture envelope of KI versus KII was determined for the coating material at ambient and elevated temperatures. Propagation angles of fracture as a function of KI/KII were also determined. The mixed-mode fracture behavior of the microsplat coating material was modeled using Finite Element approach to account for anisotropy and micro cracked structures, and predicted in terms of fracture envelope and propagation angle using mixed-mode fracture theories.

Zhu, Dongming↗

Finite Element Simulation of Solid Rocket Booster Separation Motors During Motor Firing

One of the toughest challenges facing Solid Rocket Booster (SRB) engineers is to ensure that any design changes made to the Shuttle-Derived Booster Separation Motors (BSM) for future space exploration vehicles is able to withstand the increasingly hostile motor firing environment without cracking its critical component - the graphite throat. This paper presents a critical analysis methodology and techniques for assessing effects of BSM design changes with great accuracy and precision. For current Space Shuttle operation, the motor firing occurs at SRB separation - approximately 125 seconds after Shuttle launch at an altitude of about 28 miles. The motor operation event lasts about two seconds, however, the surface temperature of the graphite throat increases approximately 3400 F in less than one second with a corresponding increase in surface pressure of approximately 2200 pounds per square inch (psi) in less than one-tenth of a second. To capture this process fully and accurately, a two-phase sequentially coupled thermal-mechanical finite element approach was developed. This method allows the time- and location-dependent pressure fields to interact with the spatial-temporal thermal fields throughout the operation. The material properties of graphite throat are orthotropic and temperature-dependent. The analysis involves preload and multiple body contacts.

Yu. Weiping↗

Influence of Cooling Channel Geometry on the Thermal Response in Silicon Nitride Plates Studied

Engine manufacturers are continually attempting to improve the performance and efficiency of internal combustion engines. Usually they raise the operating temperature or reduce the cooling air requirement for the hot section turbine components. However, the success of these attempts depends on finding materials that are lightweight, are strong, and can withstand high temperatures. Ceramics are among the top candidate materials considered for such harsh applications. They hold low-density, high-temperature strength, and thermal conductivity, and they are undergoing investigation as potential materials for replacing nickel-base alloys and superalloys that are currently used for engine hot-section components. Ceramic structures can withstand higher operating temperatures and a harsh combustion environment. In addition, their low densities relative to metals help reduce component mass. The long-term objectives of the High Temperature Propulsion Components (HOTPC) Project are to develop manufacturing technology, thermal and environmental barrier coatings (TBC/EBC), and the analytical modeling capability to predict thermomechanical stresses in minimally cooled silicon nitride turbine nozzle vanes under simulated engine conditions. Two- and three-dimensional finite element analyses with TBC were conducted at the NASA Glenn Research Center. Nondestructive evaluation was used to determine processing defects. The study included conducting preliminary parametric analytical runs of heat transfer and stress analyses under steady-state conditions to demonstrate the feasibility of using cooled Si3N4 parts for turbine applications. The influence of cooling-channel shapes (such as circular, square, and ascending-order cooling channels) on cooling efficiency and thermal stresses was investigated. Temperature distributions were generated for all cases considered under both cooling and no-cooling conditions, with air being the cooling medium. The table shows the magnitude of the maximum and minimum temperature obtained for the plates under air cooling. Each channel's cross-sectional shape delivered a different temperature; however, the two-dimensional analyses for circular and square or equal-side rectangular holes produced close results. Moreover, the model of the panel with ascending order cooling channels experienced the lowest temperature. A difference of near 260 C was found among the three cooling-hole configurations investigated. The ascending-order cooling channels arrangement showed superior performance by attaining the lowest temperature (1077 C) in comparison to the circular (1379 C) and square (1343 C) channels for the same cooling-hole size. This indicates that the panel with ascending-order cooling channels is the most suitable configuration regardless of the complexity involved in its manufacture. More details pertaining to this study are reported.

Abdul-Aziz, Ali↗

Inelastic Strain Analysis of Solder Joint in NASA Fatigue Specimen

The solder fatigue specimen designed by NASA-GSFC/UNISYS is analyzed in order to obtain the inelastic strain history during two different representative temperature cycles specified by UNISYS. In previous reports (dated July 25, 1990, and November 15, 1990), results were presented of the elastic-plastic and creep analysis for delta T = 31 C cycle, respectively. Subsequent results obtained during the current phase, from viscoplastic finite element analysis of the solder fatigue specimen for delta T = 113 C cycle are summarized. Some common information is repeated for self-completeness. Large-deformation continuum formulations in conjunction with a standard linear solid model is utilized for modeling the solder constitutive creep-plasticity behavior. Relevant material properties are obtained from the literature. Strain amplitudes, mean strains, and residual strains (as well as stresses) accumulated due to a representative complete temperature cycle are obtained as a result of this analysis. The partitioning between elastic strains, time-independent inelastic (plastic) strains, and time-dependent inelastic (creep) strains is also explicitly obtained for two representative cycles. Detailed plots are presented for two representative temperature cycles. This information forms an important input for fatigue damage models, when predicting the fatigue life of solder joints under thermal cycling

Dasgupta, Abhijit↗

The effect of energetically produced O2/+/ on the ion temperatures of the Martian thermosphere

The heating rates and ion temperatures of the Martian thermosphere resulting from the thermalization of the energetic O2(+) ions produced by the reactions of the solar ionization products CO2(+) and O(+) with neutral particles are calculated and the effects of small magnetic fields on the ion thermal balance are investigated. The energy transfer and transport of energetic ions is modelled by solving the continuity equation for the ions over a series of finite energy cells yielding the equilibrium densities as a function of energy and altitude. Particles which are able to proceed upwards without collisions are dealt with separately from the continuity equation. It is shown that the thermalization of the energetic O2(+) ions can greatly increase ion temperatures above 200 km compared to those calculated for only ambient electron heating. Current solar wind interaction models predict that small horizontal magnetic fields act to restrict the ion thermal conductivity and to increase upper altitude ion temperatures. The combined effects of these processes provide a partial agreement with measurements made by Viking 1.

Rohrbaugh, R. P.↗

Emissions Prediction and Measurement for Liquid-Fueled TVC Combustor with and without Water Injection

An investigation is performed to evaluate the performance of a computational fluid dynamics (CFD) tool for the prediction of the reacting flow in a liquid-fueled combustor that uses water injection for control of pollutant emissions. The experiment consists of a multisector, liquid-fueled combustor rig operated at different inlet pressures and temperatures, and over a range of fuel/air and water/fuel ratios. Fuel can be injected directly into the main combustion airstream and into the cavities. Test rig performance is characterized by combustor exit quantities such as temperature and emissions measurements using rakes and overall pressure drop from upstream plenum to combustor exit. Visualization of the flame is performed using gray scale and color still photographs and high-frame-rate videos. CFD simulations are performed utilizing a methodology that includes computer-aided design (CAD) solid modeling of the geometry, parallel processing over networked computers, and graphical and quantitative post-processing. Physical models include liquid fuel droplet dynamics and evaporation, with combustion modeled using a hybrid finite-rate chemistry model developed for Jet-A fuel. CFD and experimental results are compared for cases with cavity-only fueling, while numerical studies of cavity and main fueling was also performed. Predicted and measured trends in combustor exit temperature, CO and NOx are in general agreement at the different water/fuel loading rates, although quantitative differences exist between the predictions and measurements.

Brankovic, A.↗

Loading tests of a wing structure for a hypersonic aircraft

Room-temperature loading tests were conducted on a wing structure designed with a beaded panel concept for a Mach 8 hypersonic research airplane. Strain, stress, and deflection data were compared with the results of three finite-element structural analysis computer programs and with design data. The test program data were used to evaluate the structural concept and the methods of analysis used in the design. A force stiffness technique was utilized in conjunction with load conditions which produced various combinations of panel shear and compression loading to determine the failure envelope of the buckling critical beaded panels The force-stiffness data did not result in any predictions of buckling failure. It was, therefore, concluded that the panels were conservatively designed as a result of design constraints and assumptions of panel eccentricities. The analysis programs calculated strains and stresses competently. Comparisons between calculated and measured structural deflections showed good agreement. The test program offered a positive demonstration of the beaded panel concept subjected to room-temperature load conditions.

Fields, R. A.↗

Thermal elastohydrodynamic lubrication of spur gears

An analysis and computer program called TELSGE were developed to predict the variations of dynamic load, surface temperature, and lubricant film thickness along the contacting path during the engagement of a pair of involute spur gears. The analysis of dynamic load includes the effect of gear inertia, the effect of load sharing of adjacent teeth, and the effect of variable tooth stiffness which are obtained by a finite-element method. Results obtained from TELSGE for the dynamic load distributions along the contacting path for various speeds of a pair of test gears show patterns similar to that observed experimentally. Effects of damping ratio, contact ratio, tip relief, and tooth error on the dynamic load were examined. In addition, two dimensionless charts are included for predicting the maximum equilibrium surface temperature, which can be used to estimate directly the lubricant film thickness based on well established EHD analysis.

Wang, K. L.↗

Thermally driven flow in a rotating spherical shell Axisymmetic states

A spherical analogue of the rotating annulus experiments modeling atmospheric motion, in which a liquid is contained between two rigid, corotating and concentric hemispheres upon both of which thermal gradients are imposed, is presently studied by means of numerical models. Temperatures are lower on the inner than on the outer sphere, and decrease towards the pole. Using Navier-Stokes equations which assume symmetry about the polar axis, finite difference numerical models yield steady-state solutions to the equations. Hydrostatic and nonhydrostatic solutions are compared for cylindrical and spherical cases, and it is found in the case of the spherical shell that the differences between hydrostatic and nonhydrostatic solutions are small and largely confined to the regions near the pole and equator. It is suggested that nonhydrostatic effects on the axisymmetric state will not affect the flow's baroclinic stability.

Miller, T. L.↗

Structural sensitivity analysis: Methods, applications, and needs

Some innovative techniques applicable to sensitivity analysis of discretized structural systems are reviewed. These techniques include a finite-difference step-size selection algorithm, a method for derivatives of iterative solutions, a Green's function technique for derivatives of transient response, a simultaneous calculation of temperatures and their derivatives, derivatives with respect to shape, and derivatives of optimum designs with respect to problem parameters. Computerized implementations of sensitivity analysis and applications of sensitivity derivatives are also discussed. Finally, some of the critical needs in the structural sensitivity area are indicated along with Langley plans for dealing with some of these needs.

Adelman, H. M.↗

Free-edge delamination: Laminate width and loading conditions effects

The width and loading conditions effects on free-edge stress fields in composite laminates are investigated using a three-dimensional finite element analysis. This analysis includes a special free-edge region refinement or superelement with progrssive substructuring (mesh refinement) and finite thickness interply layers. The different loading conditions include in-plane and out-of-plane bending, combined axial tension and in-plane shear, twisting, uniform temperature and uniform moisture. Results obtained indicate that: axial tension causes the smallest magnitude of interlaminar free edge stress compared to other loading conditions; free-edge delamination data obtained from laboratory specimens cannot be scaled to structural components; and composite structural components are not likely to delaminate.

Murthy, P. L. N.↗