Engineering topics
Publications and source records attributed to Howser, L. M..
Transient response of thermal-protection material to heat applied to surface calculated using CHAP III computer program. Used to analyze chemical kinetics of pyrolysis gas in detail and examine pyrolysis reactions-in-depth. Analysis includes deposition of solid products produced by chemical reactions in gas phase. Uses modeling technique that approximates wide range of ablation problems. CHAP III written in FORTRAN IV.
A representative hydrodynamics problem, the shock initiated flow over a flat plate, was used for exploring data organizations and program structures needed to exploit the STAR-100 vector processing computer. A brief description of the problem is followed by a discussion of how each portion of the computational process was vectorized. Finally, timings of different portions of the program are compared with equivalent operations on serial machines. The speed up of the STAR-100 over the CDC 6600 program is shown to increase as the problem size increases. All computations were carried out on a CDC 6600 and a CDC STAR 100, with code written in FORTRAN for the 6600 and in STAR FORTRAN for the STAR 100.
A reexamination of some numerical methods is considered in light of the new class of computers which use vector streaming to achieve high computation rates. A study has been made of the effect on the relative efficiency of several numerical methods applied to a particular fluid flow problem when they are implemented on a vector computer. The method of Brailovskaya, the alternating direction implicit method, a fully implicit method, and a new method called partial implicitization have been applied to the problem of determining the steady state solution of the two-dimensional flow of a viscous imcompressible fluid in a square cavity driven by a sliding wall. Results are obtained for three mesh sizes and a comparison is made of the methods for serial computation.
Two representative algorithms used on a serial computer and presently executed on the Control Data Corporation 6000 computer were adapted to execute efficiently on the Control Data STAR-100 computer. Gaussian elimination for the solution of simultaneous linear equations and the Gauss-Legendre quadrature formula for the approximation of an integral are the two algorithms discussed. A description is given of how the programs were adapted for STAR and why these adaptations were necessary to obtain an efficient STAR program. Some points to consider when adapting an algorithm for STAR are discussed. Program listings of the 6000 version coded in 6000 FORTRAN, the adapted STAR version coded in 6000 FORTRAN, and the STAR version coded in STAR FORTRAN are presented in the appendices.
Two algorithms are presented for smoothing arbitrary sets of data. They are the explicit variable algorithm and the parametric variable algorithm. The former would be used where large gradients are not encountered because of the smaller amount of calculation required. The latter would be used if the data being smoothed were double valued or experienced large gradients. Both algorithms use a least-squares technique to obtain a cubic spline fit to the data. The advantage of the spline fit is that the first and second derivatives are continuous. This method is best used in an interactive graphics environment so that the junction values for the spline curve can be manipulated to improve the fit.
The differential equations governing the transient response of the char layer of an ablating axisymmetric body, internal pyrolysis gas flow effects being considered, have been derived. These equations have been expanded into finite difference form and programed for numerical solution on a digital computer. Numerical results compare favorably with simplified exact solutions. The complete numerical analysis was used to obtain solutions for two representative body shapes subjected to a typical entry heating environment. Pronounced effects of the lateral flow of pyrolysis gases on the mass flow field within the char layer and the associated surface and pyrolysis interface recession rates are shown.
A Fe-57 Mossbauer spectrometer, equally efficient in transmission and reflection geometries, is described. The radiation detector consists of a 1.524- by 5.08 by 5.08-cm rectangular NaI(Tl) crystal with a hole 1.524 cm in diameter. The front and back faces of the crystal are covered with beryllium windows 0.0127 cm thick and 3.81 cm in diameter. The energy of the radiation accepted for counting ranges from 6.3 keV conversion X-rays to the 14.4 keV reemitted gamma rays. The spectrometer was used to measure various types of low frequency (10 Hz) and low amplitude (0.254 mm) periodic motion of steel specimens.
A computer program to analyze Mossbauer data is presented in detail. The least-squares curve fitting techniques described apply to single line spectra, single hyperfine spectra, or when the constituent spectra are separated well enough to let the individual absorption peaks stand alone. The present program is not adapted for complex spectra resulting from the existence of several local environments in the absorber iron alloy. Sample problems are presented to aid the user in setting up and running the program. The program is written in FORTRAN 4 language for the Control Data 6000 series digital computer with the SCOPE 3.0 operating system and requires approximately 115,000 octal locations of core storage. A typical case with one absorption peak runs in 20 seconds, and a typical problem with six absorption peaks requires 50 seconds.
Governing differential equation, boundary conditions for the analysis on which computer program is based, and method of solution of resulting finite difference equations are discussed in the documentation.
Computer program for analyzing transient response of ablating axisymmetric bodies including effects of shape change
Differential equations governing transient response of ablating axisymmetric orthotropic bodies and effect of shape change