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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS1995

Analysis of composite ablators using massively parallel computation

In this work, the feasibility of using massively parallel computation to study the response of ablative materials is investigated. Explicit and implicit finite difference methods are used on a massively parallel computer, the Thinking Machines CM-5. The governing equations are a set of nonlinear partial differential equations. The governing equations are developed for three sample problems: (1) transpiration cooling, (2) ablative composite plate, and (3) restrained thermal growth testing. The transpiration cooling problem is solved using a solution scheme based solely on the explicit finite difference method. The results are compared with available analytical steady-state through-thickness temperature and pressure distributions and good agreement between the numerical and analytical solutions is found. It is also found that a solution scheme based on the explicit finite difference method has the following advantages: incorporates complex physics easily, results in a simple algorithm, and is easily parallelizable. However, a solution scheme of this kind needs very small time steps to maintain stability. A solution scheme based on the implicit finite difference method has the advantage that it does not require very small times steps to maintain stability. However, this kind of solution scheme has the disadvantages that complex physics cannot be easily incorporated into the algorithm and that the solution scheme is difficult to parallelize. A hybrid solution scheme is then developed to combine the strengths of the explicit and implicit finite difference methods and minimize their weaknesses. This is achieved by identifying the critical time scale associated with the governing equations and applying the appropriate finite difference method according to this critical time scale. The hybrid solution scheme is then applied to the ablative composite plate and restrained thermal growth problems. The gas storage term is included in the explicit pressure calculation of both problems. Results from ablative composite plate problems are compared with previous numerical results which did not include the gas storage term. It is found that the through-thickness temperature distribution is not affected much by the gas storage term. However, the through-thickness pressure and stress distributions, and the extent of chemical reactions are different from the previous numerical results. Two types of chemical reaction models are used in the restrained thermal growth testing problem: (1) pressure-independent Arrhenius type rate equations and (2) pressure-dependent Arrhenius type rate equations. The numerical results are compared to experimental results and the pressure-dependent model is able to capture the trend better than the pressure-independent one. Finally, a performance study is done on the hybrid algorithm using the ablative composite plate problem. It is found that there is a good speedup of performance on the CM-5. For 32 CPU's, the speedup of performance is 20. The efficiency of the algorithm is found to be a function of the size and execution time of a given problem and the effective parallelization of the algorithm. It also seems that there is an optimum number of CPU's to use for a given problem.

Shia, David↗

NASA NTRS1987

A numerical simulation to verify the stress-free growth of silicon crystal ribbon

Thermal stresses developed during the growth of silicon crystal ribbon have been shown to be negligible, thus eliminating residual stresses and dislocations, if the temperature profile satisfies a second-order partial differential equation inside the ribbon. This has been numerically verified through a finite element model, an outline of which is presented here. This model shows that, for homogeneous isotropic material with temperature independent thermal expansion coefficients, thermal stresses will vanish if the temperature profile satisfies the Laplacian. A comparison of stresses due to uniform and nonuniform temperature gradients in the plane of the ribbon is also presented. The strategies employed to control the round-off error and to validate the computer model are discussed.

Ray, Sujit K.↗

NASA NTRS1983

Stress studies in EFG

Experimental work in support of stress studies in high speed silicon sheet growth has been emphasized in this quarter. Creep experiments utilizing four-point bending have been made in the temperature range from 1000 C to 1360 C in CZ silicon as well as on EFG ribbon. A method to measure residual stress over large areas using laser interferometry to map strain distributions under load is under development. A fiber optics sensor to measure ribbon temperature profiles has been constructed and is being tested in a ribbon growth furnace environment. Stress and temperature field modeling work has been directed toward improving various aspects of the finite element computing schemes. Difficulties in computing stress distributions with a very high creep intensity and with non-zero interface stress have been encountered and additional development of the numerical schemes to cope with these problems is required. Temperature field modeling has been extended to include the study of heat transfer effects in the die and meniscus regions.

Source record↗

NASA NTRS1994

Superconducting axisymmetric finite elements based on a gauged potential variational principle. Part 1: Formulation

The present work is part of a research program for the numerical simulation of electromagnetic (EM) fields within conventional Ginzburg-Landau (GL) superconductors. The final goal of this research is to formulate, develop and validate finite element (FE) models that can accurately capture electromagnetic thermal and material phase changes in a superconductor. The formulations presented here are for a time-independent Ginzburg-Landau superconductor and are derived from a potential-based variational principle. We develop an appropriate variational formulation of time-independent supercontivity for the general three-dimensional case and specialize it to the one-dimensional case. Also developed are expressions for the material-dependent parameters alpha and beta of GL theory and their dependence upon the temperature T. The one-dimensional formulation is then discretized for finite element purposes and the first variation of these equations is obtained. The resultant Euler equations contain nonlinear terms in the primary variables. To solve these equations, an incremental-iterative solution method is used. Expressions for the internal force vector, external force vector, loading vector and tangent stiffness matrix are therefore developed for use with the solution procedure.

Schuler, James J.↗

NASA NTRS1977

Finite-element solutions for geothermal systems

Vector potential and scalar potential are used to formulate the governing equations for a single-component and single-phase geothermal system. By assuming an initial temperature field, the fluid velocity can be determined which, in turn, is used to calculate the convective heat transfer. The energy equation is then solved by considering convected heat as a distributed source. Using the resulting temperature to compute new source terms, the final results are obtained by iterations of the procedure. Finite-element methods are proposed for modeling of realistic geothermal systems; the advantages of such methods are discussed. The developed methodology is then applied to a sample problem. Favorable agreement is obtained by comparisons with a previous study.

Chen, J. C.↗

NASA NTRS1999

Oxygen Diffusion and Reaction Kinetics in Continuous Fiber Ceramic Matrix Composites

Previous stressed oxidation tests of C/SiC composites at elevated temperatures (350 C to 1500 C) and sustained stresses (69 MPa and 172 MPa) have led to the development of a finite difference cracked matrix model. The times to failure in the samples suggest oxidation occurred in two kinetic regimes defined by the rate controlling mechanisms (i.e. diffusion controlled and reaction controlled kinetics). Microstructural analysis revealed preferential oxidation along as-fabricated, matrix microcracks and also suggested two regimes of oxidation kinetics dependent on the oxidation temperature. Based on experimental results, observation, and theory, a finite difference model was developed. The model simulates the diffusion of oxygen into a matrix crack bridged by carbon fibers. The model facilitates the study of the relative importance of temperature, the reaction rate constant, and the diffusion coefficient on the overall oxidation kinetics.

Halbig, Michael C.↗

NASA NTRS2001

Finite Element Analysis of Active and Sensory Thermopiezoelectric Composite Materials

Analytical formulations are developed to account for the coupled mechanical, electrical, and thermal response of piezoelectric composite materials. The coupled response is captured at the material level through the thermopiezoelectric constitutive equations and leads to the inherent capability to model both the sensory and active responses of piezoelectric materials. A layerwise laminate theory is incorporated to provide more accurate analysis of the displacements, strains, stresses, electric fields, and thermal fields through-the-thickness. Thermal effects which arise from coefficient of thermal expansion mismatch, pyroelectric effects, and temperature dependent material properties are explicitly accounted for in the formulation. Corresponding finite element formulations are developed for piezoelectric beam, plate, and shell elements to provide a more generalized capability for the analysis of arbitrary piezoelectric composite structures. The accuracy of the current formulation is verified with comparisons from published experimental data and other analytical models. Additional numerical studies are also conducted to demonstrate additional capabilities of the formulation to represent the sensory and active behaviors. A future plan of experimental studies is provided to characterize the high temperature dynamic response of piezoelectric composite materials.

Lee, Ho-Jun↗

NASA NTRS1990

Use of satellite data and modeling to assess the influence of stratospheric processes on the troposphere

Current research focused on the examination of the effects of ozone heating and Newtonian cooling on the linear stability of transient planetary waves. A finite difference numerical model was developed in which the basic state wind, temperature, and ozone distributions were constructed using ground-based and satellite data. Preliminary results indicate that ozone dynamics interaction may play a more important role in the stability and maintenance of planetary waves than previously thought, particularly in summer when the mean solar zenith angle is smallest and thus ozone heating largest. Examination of the finite amplitude interactions among radiation, ozone, and dynamics was carried out. Self-consistent, coupled equations governing the weakly nonlinear interactions are currently being analyzed to provide a better understanding of wave-mean flow interactions and ozone transport in a continuously stratified model of the troposphere-stratosphere coupled system.

Nathan, Terrence↗

NASA NTRS2013

MHD Simulations of the Plasma Flow in the Magnetic Nozzle

The magnetohydrodynamic (MHD) flow of plasma through a magnetic nozzle is simulated by solving the governing equations for the plasma flow in the presence of an static magnetic field representing the applied nozzle. This work will numerically investigate the flow and behavior of the plasma as the inlet plasma conditions and magnetic nozzle field strength are varied. The MHD simulations are useful for addressing issues such as plasma detachment and to can be used to gain insight into the physical processes present in plasma flows found in thrusters that use magnetic nozzles. In the model, the MHD equations for a plasma, with separate temperatures calculated for the electrons and ions, are integrated over a finite cell volume with flux through each face computed for each of the conserved variables (mass, momentum, magnetic flux, energy) [1]. Stokes theorem is used to convert the area integrals over the faces of each cell into line integrals around the boundaries of each face. The state of the plasma is described using models of the ionization level, ratio of specific heats, thermal conductivity, and plasma resistivity. Anisotropies in current conduction due to Hall effect are included, and the system is closed using a real-gas equation of state to describe the relationship between the plasma density, temperature, and pressure.A separate magnetostatic solver is used to calculate the applied magnetic field, which is assumed constant for these calculations. The total magnetic field is obtained through superposition of the solution for the applied magnetic field and the self-consistently computed induced magnetic fields that arise as the flowing plasma reacts to the presence of the applied field. A solution for the applied magnetic field is represented in Fig. 1 (from Ref. [2]), exhibiting the classic converging-diverging field pattern. Previous research was able to demonstrate effects such as back-emf at a super-Alfvenic flow, which significantly alters the shape of the magnetic field in both the near- and far-field regions. However, in that work the downstream domain was constrained to a channel of constant cross-sectional area. In the present work we seek to address this issue by modeling the downstream region with a domain that permits free expansion of the plasma, permitting a better evaluation of the downstream effects the applied field has on the plasma. The inlet boundary conditions and applied magnetic field values will also be varied to determine the effect the initial plasma energy content and applied magnetic field energy density have on the near- and far-field plasma properties on the MHD code. This will determine the effect of inlet boundary conditions on the results downstream and address issues related to the restrictive numerical domain previously used.

Smith, T. E. R.↗

NASA NTRS1982

A theoretical study of a two-wavelength lidar technique for the measurement of atmospheric temperature profiles

The theory of differential absorption lidar measurements for lines with a Voigt profile is given and applied to a two-wavelength technique for measuring the atmospheric temperature profile using a high J line in the oxygen A band. Explicit expressions for the temperature and pressure dependence of the absorption coefficient are developed for lines with a Voigt profile. An iteration procedure for calculating the temperature for narrow laser bandwidths is described which has an accuracy better than 0.2 K for bandwidths less than 0.01/cm. To reduce the errors in lidar measurements due to uncertainties in pressure, a method for estimating the pressure from the temperature profile is described. A procedure for extending the differential absorption technique to the case of finite laser bandwidth with good accuracy is also described. Simulation results show that a knowledge of the laser frequency is needed to 0.005/cm for accurate temperature measurements. Evaluation of the sensitivity for both ground- and Shuttle-based measurements shows accuracies generally better than 1 K. This technique allows up to an order of magnitude improvement in sensitivity compared to other differential absorption lidar techniques.

Korb, C. L.↗

NASA NTRS1993

Brief summary of the evolution of high-temperature creep-fatigue life prediction models for crack initiation

The evolution of high-temperature, creep-fatigue, life-prediction methods used for cyclic crack initiation is traced from inception in the late 1940's. The methods reviewed are material models as opposed to structural life prediction models. Material life models are used by both structural durability analysts and by material scientists. The latter use micromechanistic models as guidance to improve a material's crack initiation resistance. Nearly one hundred approaches and their variations have been proposed to date. This proliferation poses a problem in deciding which method is most appropriate for a given application. Approaches were identified as being combinations of thirteen different classifications. This review is intended to aid both developers and users of high-temperature fatigue life prediction methods by providing a background from which choices can be made. The need for high-temperature, fatigue-life prediction methods followed immediately on the heels of the development of large, costly, high-technology industrial and aerospace equipment immediately following the second world war. Major advances were made in the design and manufacture of high-temperature, high-pressure boilers and steam turbines, nuclear reactors, high-temperature forming dies, high-performance poppet valves, aeronautical gas turbine engines, reusable rocket engines, etc. These advances could no longer be accomplished simply by trial and error using the 'build-em and bust-em' approach. Development lead times were too great and costs too prohibitive to retain such an approach. Analytic assessments of anticipated performance, cost, and durability were introduced to cut costs and shorten lead times. The analytic tools were quite primitive at first and out of necessity evolved in parallel with hardware development. After forty years more descriptive, more accurate, and more efficient analytic tools are being developed. These include thermal-structural finite element and boundary element analyses, advanced constitutive stress-strain-temperature-time relations, and creep-fatigue-environmental models for crack initiation and propagation. The high-temperature durability methods that have evolved for calculating high-temperature fatigue crack initiation lives of structural engineering materials are addressed. Only a few of the methods were refined to the point of being directly useable in design. Recently, two of the methods were transcribed into computer software for use with personal computers.

Halford, Gary R.↗

NASA NTRS1983

Thermal and elastohydrodynamic analysis of reciprocating rod seals in the Stirling engine

Sliding seals and pumping rings for use in Stirling engines are analyzed from an elastohydrodynamic point of view. The oil film thickness and pressure distribution are found by a finite element method and then used to determine the operating temperature of sliding seals. Thermal aspects of dry seals (cap seals) are also discussed. A parametric study has been made and the results summarized in a set of curves.

Yang, Y.↗

NASA NTRS1984

Thermocapillary convection in a liquid bridge

Results from experiments on thermocapillary flow in a vertical liquid bridge are reported. Data on surface velocities are presented for four different silicone oils at several values of the temperature difference across the liquid bridge and for three different average temperatures. Also, axial velocity measured as a function of radial position in the bridge is reported. A theoretical model of the velocity and temperature fields in the system is formulated and solved by finite differences. The experimental data on surface velocities are nondimensionalized, and are found to collapse on a single line when plotted against the Marangoni number. Scaled velocities are relatively insensitive to changes in the value of the Marangoni number from 20 to 290, and the Prandtl number from 1460 to 7120. The data are in reasonable accord with predictions from theory.

Mcneil, Thomas J.↗

NASA NTRS1991

Algorithm Computes Hypersonic Flow Of Air

Algorithm simulates three-dimensional hypersonic flow of air or another mixture of gases around blunt body. Accounts for nonequilibrium thermochemical effects via two-temperature mathematical model of molecular excitations and equations for finite-rate chemical reactions. These equations fully coupled to fluid dynamical equations.

Palmer, Grant↗

NASA NTRS1991

Interaction of the sea breeze with a river breeze in an area of complex coastal heating

The interaction of the sea-breeze circulation with a river-breeze circulation in an area of complex coastal heating (east coast of Florida) was studied using a 3D finite-element mesoscale model. The model simulations are compared with temperature and wind fields observed on a typical fall day during the Kennedy Space Center Atmospheric Boundary Layer Experiment. The results from numerical experiments designed to isolate the effect of the river breeze indicate that the convergence in the sea-breeze front is suppressed when it passes over the cooler surface of the rivers.

Zhong, Shiyuan↗

NASA NTRS1992

Stagnation-point heat-transfer rate predictions at aeroassist flight conditions

The results are presented for the stagnation-point heat-transfer rates used in the design process of the Aeroassist Flight Experiment (AFE) vehicle over its entire aeropass trajectory. The prediction methods used in this investigation demonstrate the application of computational fluid dynamics (CFD) techniques to a wide range of flight conditions and their usefulness in a design process. The heating rates were computed by a viscous-shock-layer (VSL) code at the lower altitudes and by a Navier-Stokes (N-S) code for the higher altitude cases. For both methods, finite-rate chemically reacting gas was considered, and a temperature-dependent wall-catalysis model was used. The wall temperature for each case was assumed to be radiative equilibrium temperature, based on total heating. The radiative heating was estimated by using a correlation equation. Wall slip was included in the N-S calculation method, and this method implicitly accounts for shock slip. The N-S/VSL combination of projection methods was established by comparison with the published benchmark flow-field code LAURA results at lower altitudes, and the direct simulation Monte Carlo results at higher altitude cases. To obtain the design heating rate over the entire forward face of the vehicle, a boundary-layer method (BLIMP code) that employs reacting chemistry and surface catalysis was used. The ratio of the VSL or N-S method prediction to that obtained from the boundary-layer method code at the stagnation point is used to define an adjustment factor, which accounts for the errors involved in using the boundary-layer method.

Gupta, Roop N.↗

NASA NTRS1991

Residual thermal strains and stresses in nickel aluminide matrix composites

Thermally induced residual strains and stresses developed during postfabrication cooling in Saphikon/NiAl and tungsten/NiAl high-temperature composites are investigated through three-dimensional elastoplastic finite-element analyses. Average axial and transverse strains in the matrix are found to be tensile and compressive, respectively, and similar for both Saphikon and W-fiber-reinforced NiAl composites. It is suggested that the residual matrix stresses and strains are controlled more by the low-matrix yield stress than by the fiber/matrix expansion mismatch. Residual thermal strains in the matrix of these composites are measured by using a neutron-diffraction technique; the measured axial and transverse strains in the matrix are found to be in agreement with the computed values.

Saigal, A.↗

NASA NTRS1992

One-dimensional transient finite difference model of an operational salinity gradient solar pond

This paper describes the modeling approach used to simulate the transient behavior of a salinity gradient solar pond. A system of finite difference equations are used to generate the time dependent temperature and salinity profiles within the pond. The stability of the pond, as determined by the capacity of the resulting salinity profile to suppress thermal convection within the primary gradient region of the pond, is continually monitored and when necessary adjustments are made to the thickness of the gradient zone. Results of the model are then compared to measurements taken during two representative seasonal periods at the University of Texas at El Paso's (UTEP's) research solar pond.

Hicks, Michael C.↗