SEARCH · Engineering Papers
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.
The symmetrical problem of two transverse cracks in an elastic strip with reinforced surfaces is formulated in terms of a singular integral equation. The special cases of one central crack or two edge cracks are discussed. Numerical methods for solving the problems with internal cracks are outlined and stress intensity factors are presented for various geometrics and degrees of surface reinforcement.
A detailed formulation of radiative heat transfer problems involving the use of bidirectional reflectance is discussed. Except for very simple systems this formulation is very intricate. For this reason and the fact that bidirectional data are scarce, such a formulation is not in common use. However, with the development of the digital computer there have been numerical methods developed for detailed radiative investigations using the bidirectional reflectance. Since computations using bidirectional reflectance are coming into use for spacecraft radiative studies, it is necessary to have reflectance data on the materials involved. Materials that are frequently studied are coatings which are used in controlling the thermal environment of spacecrafts. This study is primarily concerned with the bidirectional reflectance of zinc oxide which is one constituent commonly used for coatings.
A technique consisting of freezing photo-elasticity and a numerical method was used to obtain stress intensity factors for natural cracks emanating from the corner at which a hole intersects a plate surface. Geometries studied were: (1) crack depth to thickness ratios of approximately 0.2, (2) 0.5 and 0.75; (3) crack depth to crack length ratios of approximately 1.0 to 2.0; and (4) crack length to hole radius ratios of about 0.5 to 2.0. All final crack geometries were grown under monotonic loading and growth was not self similar, with most of the growth occuring through the thickness under remote extension. Stress intensity factors were determined at the intersection of the flaw border.
The intensity of turbulence and the Lagrangian correlation coefficient for a gaseous rocket combustion chamber have been determined from the experimental measurements of the tracer gas diffusion. A combination of Taylor's turbulent diffusion theory and Spalding's numerical method for solving the conservation equations of fluid mechanics was used to calculate these quantities. Taylor's theory was extended to consider the inhomogeneity of the turbulence field in the axial direction of the combustion chamber. An exponential function was used to represent the Lagrangian correlation coefficient. The results indicate that the maximum value of the intensity of turbulence is about 15% and the Lagrangian correlation coefficient drops to about 0.12 in one inch of the chamber length.
Surface pressure and heat transfer, flow separation, flow field, and oil flow patterns on the leeward side of a space shuttle orbiter model are investigated at a free stream Mach number of 6. The free stream Reynolds numbers are between 1.64 times 10 to the 7th power and 1.31 times 10 to the 8th power per meter, and the angle of attack is varied between 0 deg and 40 deg for the present experiments. The stagnation temperatures for the tests are approximately 500 K and the wall temperature is maintained at 290 K. Existing numerical methods of three-dimensional inviscid supersonic flow theory and compressible boundary layer theory are used to predict the present experimental measurements. Results of the present tests indicate two distinct types of flow separation and surface peak heating depending on the angle of attack.
The unsteady turbulent boundary layer and potential flow about a pitching airfoil are analyzed using numerical methods to determine the effect of pitch rate on the delay in forward movement of the rear flow reversal point. An explicit finite difference scheme is used to integrate the unsteady boundary layer equations, which are coupled at each instant of time to a fully unsteady and nonlinear potential flow analysis. A substantial delay in forward movement of the reversal point is demonstrated with increasing pitch rate, and it is shown that the delay results partly from the alleviation of the gradients in the potential flow, and partly from the effects of unsteadiness in the boundary layer itself. The predicted delay in flow-reversal onset, and its variation with pitch rate, are shown to be in reasonable agreement with experimental data relating to the delay in dynamic stall. From the comparisons it can be concluded (a) that the effects of time-dependence are sufficient to explain the failure of the boundary layer to separate during the dynamic overshoot, and (b) that there may be some link between forward movement of the reversal point and dynamic stall.
The postbuckling behavior of shear webs divided into rectangular panels by stiffeners (uprights) was studied numerically, using the STAGS program in which two-dimensional finite differences are used to solve buckling and nonlinear problems. Universal nondimensional parameters, suggested by linear buckling data, are found to hold for the postbuckling range. Results indicate that the postbuckling stiffness of shear webs with isotropic panels is roughly two thirds of the prebuckling stiffness. The postbuckling behavior of shear webs with isotropic and orthotropic material properties is compared. Some practical aspects of solving nonlinear problems of this type by the numerical method employed are examined.
Numerical procedures for accurate determination of elastic stress intensity factors for the general two-dimensional crack problem are reviewed. The elastic perfectly plastic state of crack tip deformation is studied by a finite element procedure. Elastic-plastic fields in the immediate vicinity of a crack tip are determined numerically by finite element procedures based on asymptotic studies of crack tip singularities in plastic materials. The small-scale yielding problem is modeled, and expressions for crack tip opening displacement, shear singularity amplitude, and plastic zone extent are derived. A finite element solution to the large-scale yielding of a circumferentially cracked round tension bar is obtained. The three-dimensional aspects of flawed structures and numerical methods of treating them are studied. Ductile fracture mechanisms, in particular crack tip fracture on the microscale, are discussed.
The coupled unsteady turbulent boundary layer and potential flow about a pitching airfoil are analyzed using numerical methods to determine the effect of pitch rate on the delay in forward movement of the rear flow reversal point. An explicit finite difference scheme is used to integrate the unsteady boundary layer equations, which are coupled at each instant of time to a fully unsteady and nonlinear potential flow analysis. A substantial delay in forward movement of the reversal point is demonstrated with increasing pitch rate, and it is shown that the delay results partly from the alleviation of the gradients in the potential flow, and partly from the effects of unsteadiness in the boundary layer itself.
Many future NASA programs require very high accurate pointing stability. These pointing requirements are well beyond anything attempted to date. This paper suggests a control system which has the capability of meeting these requirements. An optimal control law for the suggested system is specified. However, since no direct method of solution is known for this complicated system, a computation technique using successive approximations is used to develop the required solution. The method of calculus of variations is applied for estimating the changes of index of performance as well as those constraints of inequality of state variables and terminal conditions. Thus, an algorithm is obtained by the steepest descent method and/or conjugate gradient method. Numerical examples are given to show the optimal controls.
Definitions are provided of the basic concepts occurring in the solution of multiple scattering problems involving planetary atmospheres and attention is given to aspects of problem characterization. Approaches are considered for finding the answer to a particular problem without the performance of detailed calculations. The characteristics of albedos are investigated, taking into account semiinfinite atmospheres and finite atmospheres. Questions of surface illumination are discussed along with aspects related to energy deposition in the atmosphere, intensity, and polarization. Precise numerical methods are examined and analytical solutions are presented.
With applications to geophysical subsurface probings, electromagnetic fields due to a horizontal electric dipole laid on the surface of a two-layer medium are solved by a combination of analytic and numerical methods. Interference patterns are calculated for various layer thickness. The results are interpreted in terms of normal modes, and the accuracies of the methods are discussed.
The theory and computer programs, based on electromagnetic theory, for the analysis and design of echelle gratings were developed. The gratings are designed for instruments that operate in the ultraviolet portion of the spectrum. The theory was developed so that the resulting computer programs will be able to analyze deep (up to 30 wavelengths) gratings by including as many as 100 real or homogeneous diffraction orders. The program calculates the complex amplitude coefficient for each of the diffracted orders. A check on the numerical method used to solve the integral equations is provided by a conservation of energy calculation.
Multiple scattering problems, radiative transfer problems in planetary atmospheres within extended visible portion of the spectrum, are examined. The direct and inverse problems and the extinction coefficient are defined, along with other scattering characteristics. Albedos in semi-infinite and finite atmospheres are considered, as well as surface illumination, energy deposition, and polarization. The Eddington approximation figures prominently in the calculations. Precise numerical methods and analytical solutions are included.
The development of analytical techniques and numerical methods for the prediction of the stability behavior of liquid propellant rocket combustors is discussed. The studies involved nonlinear wave propagation effects and nonlinear combustion zone, liner, and nozzle responses. A second primary extension of the basic theory was directed at the inclusion of the effects of injector face baffles. Attention was restricted to the linear (small amplitude) problem, finite Mach numbers, three dimensional oscillations, and distributed combustion effects.
The numerical method of obtaining three field components from total field measurements, using double Fourier series expansion, is presented. The expressions for moments of the anomalous field components over a finite area are given. The magnitude and direction of the magnetization vector indicate that the vertical component of the magnetic field calculated from total field observations is more accurate at higher geomagnetic latitudes than at lower latitudes. The opposite is true for the horizontal components. The error in determining the magnetization vector directions are significantly large over most of the range of variation of declination and inclination of the vector, demonstrating the practical limitations of computing field components from total field data even under the best of conditions.
Steady-state and dynamic performance of a gas-lubricated, self-acting face seal was determined using numerical methods based on a variable grid, finite-difference, time-transient procedure. Results were obtained for a gas turbine main shaft seal operating at 206.9 newton per square centimeter (300 psi) sealed air pressure and 152.4 meters per second (500 ft/sec) sliding velocity. Analysis of the seal dynamics revealed that the response of the seal nosepiece to runout of the seat face is markedly affected by secondary seal friction and by nosepiece inertia. The nosepiece response was determined for various levels of secondary seal friction and seat face runout magnitudes.
Analytical investigations of fluid dynamics problems of relevance to the gaseous core nuclear reactor program are presented. The vortex type flow which appears in the nuclear light bulb concept is analyzed along with the fluid flow in the fuel inlet region for the coaxial flow gaseous core nuclear reactor concept. The development of numerical methods for the solution of the Navier-Stokes equations for appropriate geometries is extended to the case of rotating flows and almost completes the gas core program requirements in this area. The investigations demonstrate that the conceptual design of the coaxial flow reactor needs further development.