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 811 records · Page 45

Material parameter measurements at high temperatures

Alternate fixtures of techniques for the measurement of the constitutive material parameters at elevated temperatures are presented. The technique utilizes scattered field data from material coated cylinders between parallel plates or material coated hemispheres over a finite size groundplane. The data acquisition is centered around the HP 8510B Network Analyzer. The parameters are then found from a numerical search algorithm using the Newton-Ralphson technique with the measured and calculated fields from these canonical scatters. Numerical and experimental results are shown.

Dominek, A.↗

Three dimensional thermal stresses in angle-ply composite laminates

The room temperature stress distributions and shapes of a family of angle ply graphite/epoxy laminates have been obtained using a three-dimensional linear finite element analysis. The sensitivity of the corners to fiber angle variations is examined, in addition to the errors introduced by assuming planes of symmetry which do not exist in angle-ply laminates. The results show that angle ply laminates with 'clustered' plies will tend to delaminate at diagonally opposite corners, and that matrix cracks in this family of laminates will be initiated in the laminate interior.

Griffin, O. Hayden, Jr.↗

COI Structural Analysis Presentation

This report discusses the structural analysis of the Next Generation Space Telescope Mirror System Demonstrator (NMSD) developed by Composite Optics Incorporated (COI) in support of the Next Generation Space Telescope (NGST) project. The mirror was submitted to Marshall Space Flight Center (MSFC) for cryogenic testing and evaluation. Once at MSFC, the mirror was lowered to approximately 40 K and the optical surface distortions were measured. Alongside this experiment, an analytical model was developed and used to compare to the test results. A NASTRAN finite element model was provided by COI and a thermal model was developed from it. Using the thermal model, steady state nodal temperatures were calculated based on the predicted environment of the large cryogenic test chamber at MSFC. This temperature distribution was applied in the structural analysis to solve for the deflections of the optical surface. Finally, these deflections were submitted for optical analysis and comparison to the interferometer test data.

Cline, Todd↗

Reduced Equations for Calculating the Combustion Rates of Jet-A and Methane Fuel

Simplified kinetic schemes for Jet-A and methane fuels were developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) that is being developed at Glenn. These kinetic schemes presented here result in a correlation that gives the chemical kinetic time as a function of initial overall cell fuel/air ratio, pressure, and temperature. The correlations would then be used with the turbulent mixing times to determine the limiting properties and progress of the reaction. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentration of carbon monoxide as a function of fuel air ratio, pressure, and temperature. The NASA Glenn GLSENS kinetics code calculates the reaction rates and rate constants for each species in a kinetic scheme for finite kinetic rates. These reaction rates and the values obtained from the equilibrium correlations were then used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide, and NOx were obtained for both Jet-A fuel and methane.

Molnar, Melissa↗

Design/Analysis of Metal/Composite Bonded Joints for Survivability at Cryogenic Temperatures

A major design and analysis challenge for the JWST ISM structure is the metal/composite bonded joints that will be required to survive down to an operational ultra-low temperature of 30K (-405 F). The initial and current baseline design for the plug-type joint consists of a titanium thin walled fitting (1-3mm thick) bonded to the interior surface of an M555/954-6 composite truss square tube with an axially stiff biased lay-up. Metallic fittings are required at various nodes of the truss structure to accommodate instrument and lift-point bolted interfaces. Analytical experience and design work done on metal/composite bonded joints at temperatures below liquid nitrogen are limited and important analysis tools, material properties, and failure criteria for composites at cryogenic temperatures are virtually nonexistent. Increasing the challenge is the difficulty in testing for these required tools and parameters at 30K. A preliminary finite element analysis shows that failure due to CTE mismatch between the biased composite and titanium or aluminum is likely. Failure is less likely with Invar, however an initial mass estimate of Invar fittings demonstrates that Invar is not an automatic alternative. In order to gain confidence in analyzing and designing the ISM joints, a comprehensive joint development testing program has been planned and is currently running. The test program is designed for the correlation of the analysis methodology, including tuning finite element model parameters, and developing a composite failure criterion for the effect of multi-axial composite stresses on the strength of a bonded joint at 30K. The testing program will also consider stress mitigation using compliant composite layers and potential strength degradation due to multiple thermal cycles. Not only will the finite element analysis be correlated to the test data, but the FEA will be used to guide the design of the test. The first phase of the test program has been completed and the preliminary analysis has been revisited based on the test data In this work, we present an overview of the test plan, results today, and resulting design improvements.

Bartoszyk, Andrew E.↗

Multiscale Modeling of Grain Boundaries in ZrB2: Structure, Energetics, and Thermal Resistance

A combination of ab initio, atomistic and finite element methods (FEM) were used to investigate the structures, energetics and lattice thermal conductance of grain boundaries for the ultra high temperature ceramic ZrB2. Atomic models of idealized boundaries were relaxed using density functional theory. Information about bonding across the interfaces was determined from the electron localization function. The Kapitza conductance of larger scale versions of the boundary models were computed using non-equilibrium molecular dynamics. The interfacial thermal parameters together with single crystal thermal conductivities were used as parameters in microstructural computations. FEM meshes were constructed on top of microstructural images. From these computations, the effective thermal conductivity of the polycrystalline structure was determined.

Lawson, John W.↗

Finite-difference fluid dynamics computer mathematical models for the design and interpretation of experiments for space flight

Numerical methods are used to design a spherical baroclinic flow model experiment of the large scale atmosphere flow for Spacelab. The dielectric simulation of radial gravity is only dominant in a low gravity environment. Computer codes are developed to study the processes at work in crystal growing systems which are also candidates for space flight. Crystalline materials rarely achieve their potential properties because of imperfections and component concentration variations. Thermosolutal convection in the liquid melt can be the cause of these imperfections. Such convection is suppressed in a low gravity environment. Two and three dimensional finite difference codes are being used for this work. Nonuniform meshes and implicit iterative methods are used. The iterative method for steady solutions is based on time stepping but has the options of different time steps for velocity and temperature and of a time step varying smoothly with position according to specified powers of the mesh spacings. This allows for more rapid convergence. The code being developed for the crystal growth studies allows for growth of the crystal as the solid-liquid interface. The moving interface is followed using finite differences; shape variations are permitted. For convenience in applying finite differences in the solid and liquid, a time dependent coordinate transformation is used to make this interface a coordinate surface.

Roberts, G. O.↗

Development of Individually Addressable Micro-Mirror-Arrays for Space Applications

We have been developing a 32 x 32 prototype array of individually addressable Micro-Mirrors capable of operating at cryogenic temperature for Earth and Space Science applications. Micro-Mirror-Array technology has the potential to revolutionize imaging and spectroscopy systems for NASA's missions of the 21st century. They can be used as programmable slits for the Next Generation Space Telescope, as smart sensors for a steerable spectrometer, as neutral density filters for bright scene attenuation etc. The, entire fabrication process is carried out in the Detector Development Laboratory at NASA, GSFC. The fabrication process is low temperature compatible and involves integration of conventional CMOS technology and surface micro-machining used in MEMS. Aluminum is used as the mirror material and is built on a silicon substrate containing the CMOS address circuit. The mirrors are 100 microns x l00 microns in area and deflect by +/- 10 deg induced by electrostatic actuation between two parallel plate capacitors. A pair of thin aluminum torsion straps allow the mirrors to tilt. Finite-element-analysis and closed form solutions using electrostatic and mechanical torque for mirror operation were developed and the results were compared with laboratory performance. The results agree well both at room temperature and at cryogenic temperature. The development demonstrates the first cryogenic operation of two-dimensional Micro-Mirrors with bi-state operation. Larger arrays will be developed meeting requirements for different science applications. Theoretical analysis, fabrication process, laboratory test results and different science applications will be described in detail.

Dutta, Sanghamitra B.↗

A microstructural lattice model for strain oriented problems: A combined Monte Carlo finite element technique

A specialized, microstructural lattice model, termed MCFET for combined Monte Carlo Finite Element Technique, was developed which simulates microstructural evolution in material systems where modulated phases occur and the directionality of the modulation is influenced by internal and external stresses. In this approach, the microstructure is discretized onto a fine lattice. Each element in the lattice is labelled in accordance with its microstructural identity. Diffusion of material at elevated temperatures is simulated by allowing exchanges of neighboring elements if the exchange lowers the total energy of the system. A Monte Carlo approach is used to select the exchange site while the change in energy associated with stress fields is computed using a finite element technique. The MCFET analysis was validated by comparing this approach with a closed form, analytical method for stress assisted, shape changes of a single particle in an infinite matrix. Sample MCFET analytical for multiparticle problems were also run and in general the resulting microstructural changes associated with the application of an external stress are similar to that observed in Ni-Al-Cr alloys at elevated temperature.

Gayda, J.↗

Micromechanics, Fracture Mechanics and Gas Permeability of Composite Laminates for Cryogenic Storage Systems

A micromechanics method is developed to investigate microcrack propagation in a liquid hydrogen composite tank at cryogenic temperature. The unit cell is modeled using square and hexagonal shapes depends on fiber and matrix layout from microscopic images of composite laminates. Periodic boundary conditions are applied to the unit cell. The temperature dependent properties are taken into account in the analysis. The laminate properties estimated by the micromechanics method are compared with empirical solutions using constituent properties. The micro stresses in the fiber and matrix phases based on boundary conditions in laminate level are calculated to predict the formation of microcracks in the matrix. The method is applied to an actual liquid hydrogen storage system. The analysis predicts micro stresses in the matrix phase are large enough to cause microcracks in the composite. Stress singularity of a transverse crack normal to a ply-interface is investigated to predict the fracture behavior at cryogenic conditions using analytical and finite element analysis. When a transverse crack touches a ply-interface of a composite layer with same fiber orientation, the stress singularity is equal to 1/2. When the transverse crack propagates to a stiffer layer normal to the ply-direction, the singularity becomes less than 1/2 and vice versa. Finite element analysis is performed to predict the fracture toughness of a laminated beam subjected to fracture loads measured by four-point bending tests at room and cryogenic temperatures. As results, the fracture load at cryogenic temperature is significantly lower than that at room temperature. However, when thermal stresses are taken into consideration, for both cases of room and cryogenic temperatures, the difference of the fracture toughness becomes insignificant. The result indicates fracture toughness is a characteristic property, which is independent to temperature changes. The experimental analysis is performed to investigate the effect of cryogenic cycling on permeability for various composite material systems. Textile composites have lower permeability than laminated composites even with increasing number of cryogenic cycle. Nano-particles dispersed in laminated composites do not show improvement on permeability. The optical inspection is performed to investigate the microcrack propagation and void content in laminated composites and compared the microscopic results before and after cryogenic cycling.

Choi, Sukjoo↗

Progress Towards Modeling a Rapid Cycle Adsorption Pump for CO2 Compression

A Rapid Cycle Adsorption Pump (RCAP) is a competitive technology for capturing and pressuring CO2 within a Martian In-Situ Resource Utilization (ISRU) system. In an ISRU plant, CO2 from the Martian atmosphere at ~0.69-0.925kPa must first be pressured to ~101-500kPa to produce O2 and/or CH4. A RCAP pressurizes CO2 by imposing fast temperature swings on an adsorbent bed low pressure CO2 is adsorbed onto the cooled bed, and higher pressure CO2 is desorbed from the heated bed. To aid the design of a RCAP for NASA's Advanced Exploration Systems (AES) ISRU project, a finite difference thermal model of a single stack RCAP was developed in Thermal Desktop. The stack consists of one gas passage sandwiched between two sorbent beds and two cold plates (for heating/cooling each bed). The model implements adsorption/desorption physics via a linear driving force approximation in order to predict both temperature and pressure swings in the pump. The modeling approach is presented along with a discussion of its results and the current design. The model was also used to trade cooling speed when constructing the RCAP with 3D printed high thermal conductivity copper (GRCop-84) verses 3D printed aluminum (AlSi10mg). A wide assembly was modeled to predict the performance of multiple stacks in parallel. Major performance drivers were identified to be 1) the contact heat transfer to the sorbent bed, and 2) the pump's thermal mass.

Modeling↗

An Approach to Solving Enclosure Radiation Problems in A Multi-Physics Context

Thermal protection system analysis of complex features or damage sites can sometimes require modeling of high temperature enclosures. Implementing efficient and accurate view-factor algorithms required to model such problems is complex. The current work leverages the Non-equilibrium Radiation (NERO) software, which solves the radiation transport equation in a finite-volume scheme, to alleviating challenges often faced with view-factor calculations. By assuming heat transfer occurs only between grey bodies and that the medium is non-participating, computational cost of the method is significantly reduced. The enclosure physics are modeled through emitting and reflecting boundary conditions in NERO. The emitted radiative flux is dependent on the wall temperature which is a solution to the material response, obtained from Icarus, in this context. The Ares framework manages the time-advancement and exchange of the necessary data between the solvers. The surface energy balance is modified to account for the enclosure terms within the material response boundary condition. The methodology was verified against analytical solutions including radiating parallel plates, a hollow cylinder (shown in Fig. 1), and a hemisphere. Application of the methodology to inform the design of components of the Dragonfly system will be shown.

Ablation↗

An implicit-iterative solution of the heat conduction equation with a radiation boundary condition

For the problem of predicting one-dimensional heat transfer between conducting and radiating mediums by an implicit finite difference method, four different formulations were used to approximate the surface radiation boundary condition while retaining an implicit formulation for the interior temperature nodes. These formulations are an explicit boundary condition, a linearized boundary condition, an iterative boundary condition, and a semi-iterative boundary method. The results of these methods in predicting surface temperature on the space shuttle orbiter thermal protection system model under a variety of heating rates were compared. The iterative technique caused the surface temperature to be bounded at each step. While the linearized and explicit methods were generally more efficient, the iterative and semi-iterative techniques provided a realistic surface temperature response without requiring step size control techniques.

Williams, S. D.↗

Elevated temperature crack growth

Critical gas turbine engine hot section components such as blades, vanes, and combustor liners tend to develop minute cracks during early stages of operations. The ability of currently available path-independent (P-I) integrals to correlate fatigue crack propagation under conditions that simulate the turbojet engine combustor liner environment was determined. To date, an appropriate specimen design and a crack displacement measurement method were determined. Alloy 718 was selected as the analog material based on its ability to simulate high temperature behavior at lower temperatures in order to facilitate experimental measurements. Available P-I integrals were reviewed and the best approaches are being programmed into a finite element post processor for eventual comparison with experimental data. The experimental data will include cyclic crack growth tests under thermomechanical conditions, and, additionally, thermal gradients.

Kim, K. S.↗

A two-dimensional finite difference program for thermal analysis of rocket thrust chambers

A two-dimensional finite difference computer model for thermal analysis of rocket thrust chambers has been developed. The model uses an iterative scheme for calculating the temperature distribution within the chamber wall and implements a successive overrelaxation formula for a quick convergence. The inputs of the models are the dimensions of the thrust chamber wall, types of materials used, heat transfer coefficients and temperatures of the hot gas and the coolant. The resulting output of the program consists of the nodal temperature distribution, heat transfer to the coolant and heat transfer from the hot gas.

Naraghi, Mohammad H.↗

Transient conduction-radiation analysis of an absolute active cavity radiometer using finite elements

A NASA-developed finite element-based model of a generic active cavity radiometer (ACR) has been developed in order to study the dependence on operating temperature of the closed-loop and open-loop transient response of the instrument. Transient conduction within the sensing element is explored, and the transient temperature distribution resulting from the application of a time-varying radiative boundary condition is calculated. The results verify the prediction that operation of an ACR at cryogenic temperatures results in large gains in frequency response.

Mahan, J. R.↗

Temperature induced distortions in space telescope mirrors

In this paper, it is illustrated how measured instantaneous coefficients of thermal expansion (CTE) can be accurately taken into account when modeling the structural behavior of space based optical systems. In particular, the importance of including CTE spatial variations in the analysis of optical elements is emphasized. A comparison is made between the CTE's of three optical materials commonly used in the construction of space mirrors (ULE, Zerodur, and beryllium). The overall impact that selection of any one of these materials has on thermal distortions is briefly discussed. As an example of how temperature dependent spatial variations in thermal strain can be accurately incorporated in the thermo-structural analysis of a precision optical system, a finite element model is developed, which is used to estimate the thermally induced distortions in the Hubble Space Telescope's (HST) primary mirror. In addition to the structural analysis, the optical aberrations due to thermally induced distortions are also examined. These calculations indicate that thermal distortions in HST's primary mirror contribute mainly to defocus error with a relatively small contribution to spherical aberration.

Nied, H. F.↗

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗