A special investigation to develop a general method for three-dimensional photoelastic stress analysis
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Explore the source record for details and available documents.
The method of strain measurement after annealing is reviewed and found to be satisfactory for the materials available in this country. A new general method is described for the photoelastic determination of the principal stresses at any point of a general body subjected to arbitrary load. The method has been applied to a sphere subjected to diametrical compressive loads. The results show possibilities of high accuracy.
The strategic routing of aircraft through navigable and controlled airspace to minimize adverse noise impact over sensitive areas is critical in the proper management and planning of the U.S. based airport system. A major objective of this phase of research is to identify, inventory, characterize, and analyze the various environmental, land planning, and regulatory data bases, along with potential three dimensional software and hardware systems that can be potentially applied for an impact assessment of any existing or planned air route. There are eight data bases that have to be assembled and developed in order to develop three dimensional aircraft route impact methodology. These data bases which cover geographical information systems, sound metrics, land use, airspace operational control measures, federal regulations and advisories, census data, and environmental attributes have been examined and aggregated. A three dimensional format is necessary for planning, analyzing space and possible noise impact, and formulating potential resolutions. The need to develop this three dimensional approach is essential due to the finite capacity of airspace for managing and planning a route system, including airport facilities. It appears that these data bases can be integrated effectively into a strategic aircraft noise routing system which should be developed as soon as possible, as part of a proactive plan applied to our FAA controlled navigable airspace for the United States.
The purpose of this study is to investigate the performance of the integral equation computations using numerical source field-panel method in a massively parallel processing (MPP) environment. A comparative study of computational performance of the MPP CM-5 computer and conventional Cray-YMP supercomputer for a three-dimensional flow problem is made. A serial FORTRAN code is converted into a parallel CM-FORTRAN code. Some performance results are obtained on CM-5 with 32, 62, 128 nodes along with those on Cray-YMP with a single processor. The comparison of the performance indicates that the parallel CM-FORTRAN code near or out-performs the equivalent serial FORTRAN code for some cases.
The designers of the Jupiter Icy Moons Orbiter (JIMO) are investigating the potential of nuclear powered-electric propulsion technology to provide deep space propulsion. In one design scenario a closed-Brayton-cycle power converter is used to convert thermal energy from a nuclear reactor to electrical power for the spacecraft utilizing an inert gas as the working fluid to run a turboalternator as described in L.S. Mason, "A Power Conversion for the Jupiter Icy Moons Orbiter," Journal of Propulsion and Power, vol. 20, no. 5, pp. 902-910. A key component in the turboalternator is the radial flow turbine wheel which may be fabricated from a cast superalloy. This turbine wheel is envisioned to run continuously over the life of the mission, which is anticipated to be about ten years. This scenario places unusual material requirements on the turbine wheel. Unlike the case of terrestrial turbine engines, fatigue, associated with start-up and shut-down of the engine, foreign-object damage, and corrosion issues are insignificant and thus creep issues become dominate. The purpose of this paper is to present estimates for creep growth of a prototypical JIMO turbine wheel over a ten year life. Since an actual design and bill of materials does not exist, the results presented in this paper are based on preliminary concepts which are likely to evolve over time. For this reason, as well as computational efficiency, a simplified 2-D, in lieu of a 3-D, viscoelastic, finite element model of a prototypical turbine wheel will be utilized employing material properties for the cast superalloy MAR-M247. The creep data employed in this analysis are based on preliminary data being generated at NASA Glenn Research Center.
Computational analyses have been conducted on the Wall-mounted Glauert-Goldschmied type body ("hump" model) with the Full Unstructured Navier-Stokes 2-D (FUN2D) flow solver developed at NASA LaRC. This investigation uses the time-accurate Reynolds-averaged Navier- Stokes (RANS) approach to predict aerodynamic performance of the active flow control experimental database for the hump model. The workshop is designed to assess the current capabilities of different classes of turbulent flow solution methodologies, such as RANS, to predict flow fields induced by synthetic jets and separation control geometries. The hump model being studied is geometrically similar to that previously tested both experimentally and computationally at NASA LaRC [ref. 1 and 2, respectively].
The objective of this study is to investigate the effect of nozzle out-of-roundness on the transient startup side loads. The out-of-roundness could be the result of asymmetric loads induced by hardware attached to the nozzle, asymmetric internal stresses induced by previous tests and/or deformation, such as creep, from previous tests. The rocket engine studied encompasses a regeneratively cooled thrust chamber and a film cooled nozzle extension with film coolant distributed from a turbine exhaust manifold. The computational methodology is based on an unstructured-grid, pressure-based computational fluid dynamics formulation, and a transient inlet history based on an engine system simulation. Transient startup computations were performed with the out-of-roundness achieved by four degrees of ovalization of the nozzle: one perfectly round, one slightly out-of-round, one more out-of-round, and one significantly out-of-round. The computed side load physics caused by the nozzle out-of-roundness and its effect on nozzle side load are reported and discussed.
The objective of this study is to investigate the effect of nozzle out-of-roundness on the transient startup side loads at a high altitude, with an anchored computational methodology. The out-of-roundness could be the result of asymmetric loads induced by hardware attached to the nozzle, asymmetric internal stresses induced by previous tests, and deformation, such as creep, from previous tests. The rocket engine studied encompasses a regeneratively cooled thrust chamber and a film cooled nozzle extension with film coolant distributed from a turbine exhaust manifold. The computational methodology is based on an unstructured-grid, pressure-based computational fluid dynamics formulation, and a transient inlet history based on an engine system simulation. Transient startup computations were performed with the out-of-roundness achieved by four different degrees of ovalization: one perfectly round, one slightly out-of-round, one more out-of-round, and one significantly out-of-round. The results show that the separation-line-jump is the peak side load physics for the round, slightly our-of-round, and more out-of-round cases, and the peak side load increases as the degree of out-of-roundness increases. For the significantly out-of-round nozzle, however, the peak side load reduces to comparable to that of the round nozzle and the separation line jump is not the peak side load physics. The counter-intuitive result of the significantly out-of-round case is found to be related to a side force reduction mechanism that splits the effect of the separation-line-jump into two parts, not only in the circumferential direction and most importantly in time.
Detailed two-dimensional finite element analyses of the cross-sections of a model CVI (chemical vapor infiltrated) SiC/SiC (silicon carbide fiber in a silicon carbide matrix) ceramic matrix composites are performed. High resolution images of the cross-section of this composite material are generated using serial sectioning of the test specimens. These images are then used to develop very detailed finite element models of the cross-sections using the public domain software OOF2 (Object Oriented Analysis of Material Microstructures). Examination of these images shows that these microstructures have significant variability and irregularity. How these variabilities manifest themselves in the variability in effective properties as well as the stress distribution, damage initiation and damage progression is the overall objective of this work. Results indicate that even though the macroscopic stress-strain behavior of various sections analyzed is very similar, each section has a very distinct damage pattern when subjected to in-plane tensile loads and this damage pattern seems to follow the unique architectural and microstructural details of the analyzed sections.
A method for computing three-dimensional turbulent subsonic flow in curved ducts is described. An approximate set of governing equations is given for viscous flows which have a primary flow direction. The derivation is coordinate invariant, and the resulting equations are expressed in terms of tensors. General tube-like coordinates were developed for a general class of geometries applicable to many internal flow problems. The coordinates are then particularized to pipes having superelliptic cross sections whose shape can vary continuously between a circle and a near rectangle. The analysis is applied to a series of relevant aerodynamic problems including transition from nearly square to round pipes and flow through a pipe with an S-shaped bend.
A three-dimensional solution procedure for the approximate form of the Navier-Stokes equation was exercised in the two- and three-dimensional modes to compute the unsteady turbulent boundary layer on a flat plate corresponding to the data of Karlsson. The procedure is based on the use of a consistently split Linearized Block Implicit technique in conjunction with a QR operator scheme. New time-dependent upstream boundary conditions were developed that yielded realistic solutions for the interior in the vicinity of the upstream boundary. Comparisons of the computation employing these boundary conditions with the data indicate that both qualitative and quantitative agreement was obtained for the mean velocity and the in phase and out of phase components of the first harmonic of the velocity. In addition, the calculation gave results for the skin friction phase angle that had expected physical behavior for large distances downstream of the inflow boundary. For the three-dimensional case, the two-dimensional data of Karlsson was considered, but in a coordinate system skewed at 45 deg to the free stream direction. The results of the calculations were in excellent agreement with the data and the two-dimensional computations.
The three-dimensional boundary-layer equations in the limit as the normal coordinate tends to infinity are called the surface Euler equations. The present paper describes an accurate method for generating edge conditions for three-dimensional boundary-layer codes using these equations. The inviscid pressure distribution is first interpolated to the boundary-layer grid. The surface Euler equations are then solved with this pressure field and a prescribed set of initial and boundary conditions to yield the velocities along the two surface coordinate directions. Results for typical wing and fuselage geometries are presented. The smoothness and accuracy of the edge conditions obtained are found to be superior to the conventional interpolation procedures.
A compressible stability analysis computer code is developed. The code uses a matrix finite-difference method for local eigenvale solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used inital-value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed that may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black-box stability code. Sample stability calculations are presented for the boundary layer profiles of an LFC swept wing.
Variational methods (VM) sensitivity analysis employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.
The paper develops an efficient free-vibration analysis procedure of two-dimensional structures. This is achieved by employing a discretization technique based on a recently developed concept of finite dynamic elements, involving higher order dynamic correction terms in the associated stiffness and inertia matrices. A plane rectangular dynamic element is developed in detail. Numerical solution results of free-vibration analysis presented herein clearly indicate that these dynamic elements combined with a suitable quadratic matrix eigenproblem solution technique effect a most economical and efficient solution for such an analysis when compared with the usual finite element method.
The FORTRAN IV Program developed to analyze the flow field associated with scramjet exhaust systems is presented. The instructions for preparing input and interpreting output are described. The program analyzes steady three dimensional supersonic flow by the reference plane characteristic technique. The governing equations and numerical techniques employed are presented in Volume 1 of this report.
An approximate analysis is presented for calculating three-dimensional, low Mach number, laminar viscous flows in curved passages with large secondary flows and corner boundary layers. The analysis is based on the decomposition of the overall velocity field into inviscid and viscous components with the overall velocity being determined from superposition. An incompressible vorticity transport equation is used to estimate inviscid secondary flow velocities to be used as corrections to the potential flow velocity field. A parabolized streamwise momentum equation coupled to an adiabatic energy equation and global continuity equation is used to obtain an approximate viscous correction to the pressure and longitudinal velocity fields. A collateral flow assumption is invoked to estimate the viscous correction to the transverse velocity fields. The approximate analysis is solved numerically using an implicit ADI solution for the viscous pressure and velocity fields. An iterative ADI procedure is used to solve for the inviscid secondary vorticity and velocity fields. This method was applied to computing the flow within a turbine vane passage with inlet flow conditions of M = 0.1 and M = 0.25, Re = 1000 and adiabatic walls, and for a constant radius curved rectangular duct with R/D = 12 and 14 and with inlet flow conditions of M = 0.1, Re = 1000, and adiabatic walls.