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 955 records · Page 53

Electron Diffusion Regions in Magnetotail Reconnection Under Varying Guide Fields

Kinetic structures of electron diffusion regions (EDRs) under finite guide fields in magnetotail reconnection are reported. The EDRs with guide fields 0.14–0.5 (in unit of the reconnecting component) are detected by the Magnetospheric Multiscale spacecraft. The key new features include the following: (1) cold inflowing electrons accelerated along the guide field and demagnetized at the magnetic field minimum while remaining a coherent population with a low perpendicular temperature, (2) wave fluctuations generating strong perpendicular electron flows followed by alternating parallel flows inside the reconnecting current sheet under an intermediate guide field, and (3) gyrophase bunched electrons with high parallel speeds leaving the X‐line region. The normalized reconnection rates for the three EDRs range from 0.05 to 0.3. The measurements reveal that finite guide fields introduce new mechanisms to break the electron frozen‐in condition.

Chen, L.‐J.↗

A finite element method for the thermochemical decomposition of polymeric materials. II - Carbon phenolic composites

A previously developed formulation for modeling the thermomechanical behavior of chemically decomposing, polymeric materials is verified by simulating the response of carbon phenolic specimens during two high temperature tests: restrained thermal growth and free thermal expansion. Plane strain and plane stress models are used to simulate the specimen response, respectively. In addition, the influence of the poroelasticity constants upon the specimen response is examined through a series of parametric studies.

Sullivan, R. M.↗

Finite element large deflection random response of thermally buckled plates

The governing nonlinear equations of motion are derived for a rectangular isotropic plate subjected to a combined thermal-acoustic loading. The thermal loading is taken to be a steady state temperature distribution, and the acoustic loading is considered to be a stationary Gaussian pressure with zero mean and uniform magnitude and phase over the plate surface. With these types of loads, the thermal-acoustic analysis can be separated into two problems: (1) thermal post-buckling; and (2) random vibration of the thermally buckled plate. Since both analyses are nonlinear, iterative techniques are used to solve each. For the random vibration analysis, the method of equivalent linearization is used; both auto- and cross-correlation terms are included. Results obtained using the present formulation are compared with previous classical equivalent linearization results.

Locke, J. E.↗

Thermal Stability of a 4 Meter Primary Reflector for the Scanning Microwave Limb Sounder

The Scanning Microwave Limb Sounder (SMLS) is a space-borne heterodyne radiometer which will measure pressure, temperature and atmospheric constituents from thermal emission in [180,680] GHz. SMLS, planned for the NRC Decadal Survey's Global Atmospheric Composition Mission, uses a novel toric Cassegrain antenna to perform both elevation and azimuth scanning. This provides better horizontal and temporal resolution and coverage than were possible with elevation-only scanning in the two previous MLS satellite instruments. SMLS is diffraction-limited in the vertical plane but highly astigmatic in the horizontal (beam aspect ratio approx. 1:20). Nadir symmetry ensures that beam shape is nearly invariant over plus or minus 65 deg azimuth. A low-noise receiver FOV is swept over the reflector system by a small azimuth-scanning mirror. We describe the fabrication and thermal-stability test of a composite demonstration primary reflector, having full 4m height and 1/3 the width planned for flight. Using finite-element models of reflectors and structure, we evaluate thermal deformations and optical performance for 4 orbital environments and isothermal soak. We compare deformations with photogrammetric measurements made during soak tests in a chamber. The test temperature range exceeds predicted orbital ranges by large factors, implying in-orbit thermal stability of 0.21 micron rms (root mean square)/C, which meets SMLS requirements.

thermal distortion↗

Simplified Two-Time Step Method for Calculating Combustion and Emission Rates of Jet-A and Methane Fuel With and Without Water Injection

A simplified kinetic scheme for Jet-A, and methane fuels with water injection was developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) or even simple FORTRAN codes. The two time step method is either an initial time averaged value (step one) or an instantaneous value (step two). The switch is based on the water concentration in moles/cc of 1x10(exp -20). The results presented here results in a correlation that gives the chemical kinetic time as two separate functions. This two time step method is used as opposed to a one step time averaged method previously developed to determine the chemical kinetic time with increased accuracy. The first time averaged step is used at the initial times for smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used with higher water concentrations, gives the chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure and temperature (T4). The simple correlations would then be compared to the turbulent mixing times to determine the limiting rates of the reaction. 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 are used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide and NOx are obtained for Jet-A fuel and methane with and without water injection to water mass loadings of 2/1 water to fuel. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, and pressure.

Molnar, Melissa↗

Thermal cycling of tungsten-fibre-reinforced superalloy composites

The thermal cycling of a tungsten-fiber-reinforced superalloy (TFRS) composite is typical of its application in high-temperature engine environments. The mismatch in thermal expansion coefficients between fiber and matrix causes substantial longitudinal (0 deg) stresses in the composite, which can produce inelastic damage-producing matrix strains. The case of thermal fatigue is explored as a "worst case" of the possible matrix damage, in comparison with specimens which are also mechanically loaded in tension. The thermally generated cyclic stresses and the attendant matrix plasticity may be estimated using a nonlinear finite-element program, by proposing a physical analog to the micromechanics equations. A damage metric for the matrix is proposed using the Coffin-Manson criterion, which metric can facilitate comparisons of damage among different candidate materials, and also comparisons for a given material subjected to different temperature cycles. An experimental program was carried out for thermal cycling of a 37 vol pct TFRS composite to different maximum temperatures. The results confirm the prediction that thermal cycling produces matrix degradation and composite strength reduction, which become more pronounced with increasing maximum cyclic temperature. The strength of the fiber is shown to be identical for the as-fabricated and thermally cycled specimens, suggesting that the reduction in composite strength is due to the loss of matrix contribution and also to notching effects of the matrix voids on the fiber.

Wetherhold, Robert C.↗

Structural Anomalies Detected in Ceramic Matrix Composites Using Combined Nondestructive Evaluation and Finite Element Analysis (NDE and FEA)

Most reverse engineering approaches involve imaging or digitizing an object and then creating a computerized reconstruction that can be integrated, in three dimensions, into a particular design environment. The rapid prototyping technique builds high-quality physical prototypes directly from computer-aided design files. This fundamental technique for interpreting and interacting with large data sets is being used here via Velocity2 (an integrated image-processing software, ref. 1) using computed tomography (CT) data to produce a prototype three-dimensional test specimen model for analyses. A study at the NASA Glenn Research Center proposes to use these capabilities to conduct a combined nondestructive evaluation (NDE) and finite element analysis (FEA) to screen pretest and posttest structural anomalies in structural components. A tensile specimen made of silicon nitrite (Si3N4) ceramic matrix composite was considered to evaluate structural durability and deformity. Ceramic matrix composites are being sought as candidate materials to replace nickel-base superalloys for turbine engine applications. They have the unique characteristics of being able to withstand higher operating temperatures and harsh combustion environments. In addition, their low densities relative to metals help reduce component mass (ref. 2). Detailed three-dimensional volume rendering of the tensile test specimen was successfully carried out with Velocity2 (ref. 1) using two-dimensional images that were generated via computed tomography. Subsequent, three-dimensional finite element analyses were performed, and the results obtained were compared with those predicted by NDE-based calculations and experimental tests. It was shown that Velocity2 software can be used to render a three-dimensional object from a series of CT scan images with a minimum level of complexity. The analytical results (ref. 3) show that the high-stress regions correlated well with the damage sites identified by the CT scans and the experimental data. Furthermore, modeling of the voids collected via NDE offered an analytical advantage that resulted in more accurate assessments of the material s structural strength. The top figure shows a CT scan image of the specimen test section illustrating various hidden structural entities in the material and an optical image of the test specimen considered in this study. The bottom figure represents the stress response predicted from the finite element analyses (ref .3 ) for a selected CT slice where it clearly illustrates the correspondence of the high stress risers due to voids in the material with those predicted by the NDE. This study is continuing, and efforts are concentrated on improving the modeling capabilities to imitate the structural anomalies as detected.

Abdul-Aziz, Ali↗

New Reduced Two-Time Step Method for Calculating Combustion and Emission Rates of Jet-A and Methane Fuel With and Without Water Injection

A simplified kinetic scheme for Jet-A, and methane fuels with water injection was developed to be used in numerical combustion codes, such as the National Combustor Code (NCC) or even simple FORTRAN codes that are being developed at Glenn. The two time step method is either an initial time averaged value (step one) or an instantaneous value (step two). The switch is based on the water concentration in moles/cc of 1x10(exp -20). The results presented here results in a correlation that gives the chemical kinetic time as two separate functions. This two step method is used as opposed to a one step time averaged method previously developed to determine the chemical kinetic time with increased accuracy. The first time averaged step is used at the initial times for smaller water concentrations. This gives the average chemical kinetic time as a function of initial overall fuel air ratio, initial water to fuel mass ratio, temperature, and pressure. The second instantaneous step, to be used with higher water concentrations, gives the chemical kinetic time as a function of instantaneous fuel and water mole concentration, pressure and temperature (T4). The simple correlations would then be compared to the turbulent mixing times to determine the limiting properties of the reaction. 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 were then used to calculate the necessary chemical kinetic times. Chemical kinetic time equations for fuel, carbon monoxide and NOx were obtained for Jet-A fuel and methane with and without water injection to water mass loadings of 2/1 water to fuel. A similar correlation was also developed using data from NASA's Chemical Equilibrium Applications (CEA) code to determine the equilibrium concentrations of carbon monoxide and nitrogen oxide as functions of overall equivalence ratio, water to fuel mass ratio, pressure and temperature (T3). The temperature of the gas entering the turbine (T4) was also correlated as a function of the initial combustor temperature (T3), equivalence ratio, water to fuel mass ratio, and pressure.

Molnar, Melissa↗

Low-cycle fatigue analysis of a cooled copper combustion chamber

A three-dimensional finite element elastoplastic strain analysis was performed for the throat section of regeneratively cooled rocket engine combustion chamber. The analysis included thermal and pressure loads, and the effects of temperature dependent material properties, to determine the strain range corresponding to the engine operating cycle. The strain range was used in conjunction with OFHC copper isothermal fatigue test data to predict engine low-cycle fatigue life. The analysis was performed for chamber configuration and operating conditions corresponding to a hydrogen-oxygen chamber which was fatigue tested to failure at the NASA Lewis Research Center.

Miller, R. W.↗

An investigation of temperature distribution in cooled guide vanes

A numerical study to determine the temperature distribution in the guide vane blades of a radial inflow turbine is presented. A computer program was developed which permits the temperature distribution to be calculated when the blade is cooled internally using a combination of impingement and film cooling techniques. The study is based on the use of the finite difference method in a two dimensional heat conduction problem. The results are then compared to determine the best cooling configuration for a certain coolant to primary mass flow ratio.

Kotwal, R.↗

Techniques for thermal sensitivity analysis

Techniques for computing the sensitivity of temperatures to changes in design variables needed for designing thermal protection systems are considered. It is shown that the choice of the most efficient technique depends on the ratio of number of temperature constraints to the number of design variables, as well as on the thermal analysis method employed. The analysis is specialized to the case of a structure modeled by finite elements, and an example of an insulation panel is used to demonstrate the techniques.

Haftka, R. T.↗

On making large nonlinear problems small

Reduction methods for solving large-scale nonlinear problems are considered. Attention is given to: (1) the selection of basis vectors for steady-state problems; (2) the identification and determination of bifurcation and limit points, including tracing post-limit-point and post-bifurcation-point paths using reduction methods; (3) the application of reduction methods to nonlinear problems with prescribed nonzero values of the fundamental unknowns; and (4) the use of reduction methods in conjunction with multifield (mixed) finite element models. The effectiveness of using reduction methods is demonstrated on the basis of several numerical examples including two-dimensional steady conduction in a square plate with temperature dependent thermal conductivity and a shallow spherical cap subjected to a central-point load and a ring load.

Noor, A. K.↗

Computation of consistent boundary quantities in finite element thermal-fluid solutions

The consistent boundary quantity method for computing derived quantities from finite element nodal variable solutions is investigated. The method calculates consistent, continuous boundary surface quantities such as heat fluxes, flow velocities, and surface tractions from nodal variables such as temperatures, velocity potentials, and displacements. Consistent and lumped coefficient matrix solutions for such problems are compared. The consistent approach may produce more accurate boundary quantities, but spurious oscillations may be produced in the vicinity of discontinuities. The uncoupled computations of the lumped approach provide greater flexibility in dealing with discontinuities and provide increased computational efficiency. The consistent boundary quantity approach can be applied to solution boundaries other than those with Dirichlet boundary conditions, and provides more accurate results than the customary method of differentiation of interpolation polynomials.

Thornton, E. A.↗

Elevated temperature crack growth

The purpose is to determine the ability of currently available P-I integrals to correlate fatigue crack propagation under conditions that simulate the turbojet engine combustor liner environment. The utility of advanced fracture mechanics measurements will also be evaluated during the course of the program. To date, an appropriate specimen design, a crack displacement measurement method, and boundary condition simulation in the computational model of the specimen were achieved. Alloy 718 was selected as an analog material based on its ability to simulate high temperature behavior at lower temperatures. Tensile and cyclic tests were run at several strain rates so that an appropriate constitutive model could be developed. Suitable P-I integrals were programmed into a finite element post-processor for eventual comparison with experimental data.

Malik, S. N.↗

Vapor condensation at the free surface of an axisymmetric liquid mixed by a laminar jet

This paper presents numerical solutions of jet-induced mixing in a partially full cryogenic tank. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to slove the nondimensional form of steady-state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of previously stated parameters on the condensation Nusselt and Stanton numbers that characterize the steady-state interface condensation process are investigated. Detailed analysis is performed to gain a better understanding of the fundamentals of fluid mixing and interface condensation.

Lin, Chin-Shun↗

Development of an integrated BEM approach for hot fluid structure interaction: BEST-FSI: Boundary Element Solution Technique for Fluid Structure Interaction

As part of the continuing effort at NASA LeRC to improve both the durability and reliability of hot section Earth-to-orbit engine components, significant enhancements must be made in existing finite element and finite difference methods, and advanced techniques, such as the boundary element method (BEM), must be explored. The BEM was chosen as the basic analysis tool because the critical variables (temperature, flux, displacement, and traction) can be very precisely determined with a boundary-based discretization scheme. Additionally, model preparation is considerably simplified compared to the more familiar domain-based methods. Furthermore, the hyperbolic character of high speed flow is captured through the use of an analytical fundamental solution, eliminating the dependence of the solution on the discretization pattern. The price that must be paid in order to realize these advantages is that any BEM formulation requires a considerable amount of analytical work, which is typically absent in the other numerical methods. All of the research accomplishments of a multi-year program aimed toward the development of a boundary element formulation for the study of hot fluid-structure interaction in Earth-to-orbit engine hot section components are detailed. Most of the effort was directed toward the examination of fluid flow, since BEM's for fluids are at a much less developed state. However, significant strides were made, not only in the analysis of thermoviscous fluids, but also in the solution of the fluid-structure interaction problem.

Dargush, G. F.↗

Sonic boom interaction with turbulence

A recently developed transonic small-disturbance model is used to analyze the interactions of random disturbances with a weak shock. The model equation has an extended form of the classic small-disturbance equation for unsteady transonic aerodynamics. It shows that diffraction effects, nonlinear steepening effects, focusing and caustic effects and random induced vorticity fluctuations interact simultaneously to determine the development of the shock wave in space and time and the pressure field behind it. A finite-difference algorithm to solve the mixed-type elliptic hyperbolic flows around the shock wave is presented. Numerical calculations of shock wave interactions with various deterministic vorticity and temperature disturbances result in complicate shock wave structures and describe peaked as well as rounded pressure signatures behind the shock front, as were recorded in experiments of sonic booms running through atmospheric turbulence.

Zvi Rusak↗

Modeling of 3-D Woven Ceramic Matrix Composites

Three different approaches are being pursued at the NASA Glenn Research Center to predict the nanostructural behavior of three-dimensional woven ceramic matrix composites. These are: a micromechanics-based approach using W-CEMCAN (Woven Ceramic Matrix Composite Analyzer), a laminate analogy method and a structural frame approach (based on the finite element method). All three techniques are applied to predict the thermomechanical properties of a three-dimensional woven angle interlock C/SiC composite. The properties are predicted for room temperature and 1100 C and the predicted properties are compared to measurements. General observations regarding the three approaches for three-dimensional composite modeling are discussed.

Murthy, Pappu L. N.↗