ALBERT

Description: The test capabilities of the Stability Wind Tunnel of the Virginia Polytechnic Institute and State University are described, and calibrations for curved and rolling flow techniques are given. Oscillatory snaking tests to determine pure yawing derivatives are considered. Representative aerodynamic data obtained for a current fighter configuration using the curved and rolling flow techniques are presented. The application of dynamic derivatives obtained in such tests to the analysis of airplane motions in general, and to high angle of attack flight conditions in particular, is discussed.

Description: The compressible dynamic stall flowfield over a NACA 0012 airfoil transiently pitching from 0 to 60 deg at a constant rate under compressible flow conditions has been studied using real-time interferometry. A quantitative description of the overall flowfield, including the finer details of dynamic stall vortex formation, growth, and the concomitant changes in the airfoil pressure distribution, has been provided by analyzing the interferograms. For Mach numbers above 0.4, small multiple shocks appear near the leading edge and are present through the initial stages of dynamic stall. Dynamic stall was found to occur coincidentally with the bursting of the separation bubble over the airfoil. Compressibility was found to confine the dynamic stall vortical structure closer to the airfoil surface. The measurements show that the peak suction pressure coefficient drops with increasing freestream Mach number, and also it lags the steady flow values at any given angle of attack. As the dynamic stall vortex is shed, an anti-clockwise vortex is induced near the trailing edge, which actively interacts with the post-stall flow.

Description: The effect of the porous leading edge of an airfoil on the blade-vortex interaction noise, which dominates the far-field acoustic spectrum of the helicopter, is investigated. The thin-layer Navier-Stokes equations are solved with a high-order upwind-biased scheme and a multizonal grid system. The Baldwin-Lomax turbulence model is modified for considering transpiration on the surface. The amplitudes of the propagating acoustic wave in the near field are calculated directly from the computation. The porosity effect on the surface is modeled in two ways: (1) imposition of prescribed transpiration velocity distribution and (2) calculation of transpiration velocity distribution by Darcy's law. Results show leading-edge transpiration can suppress pressure fluctuations at the leading edge during blade-vortex interaction and consequently reduce the amplitude of propagating noise by 30% at a maximum in the near field.

Description: A method has been developed for calculating the viscous flow about airfoils with and without deflected flaps at -90 deg incidence. This method provides for the solution of the unsteady incompressible Navier-Stokes equations by means of an implicit technique. The solution is calculated on a body-fitted computational mesh using a staggered-grid method. The vorticity is defined at the node points, and the velocity components are defined at the mesh-cell sides. The staggered-grid orientation provides for accurate representation of vorticity at the node points and the continuity equation at the mesh-cell centers. The method provides for the noniterative solution of the flowfield and satisfies the continuity equation to machine zero at each time step. The method is evaluated in terms of its stability to predict two-dimensional flow about an airfoil at -90-deg incidence for varying Reynolds number and laminar/turbulent models. The variations of the average loading and surface pressure distribution due to flap deflection, Reynolds number, and laminar or turbulent flow are presented and compared with experimental results. The comparisom indicate that the calculated drag and drag reduction caused by flap deflection and the calculated average surface pressure are in excellent agreement with the measured results at a similar Reynolds number.

Description: High-resolution (0.01/cm) absorption spectra of lean mixtures of CH4 in dry air were recorded with the McMath-Pierce Fourier transform spectrometer (FTS) of the National Solar Observatory on Kitt Peak at various temperatures between 24 and -61 C. The spectra have been analyzed to determine the values at room temperature of pressure-broadened widths and pressure-induced shifts of more than 740 transitions. The temperature dependence of air-broadened widths and pressure-induced shifts was deduced for approx. 370 transitions in the nu(sub 1) + nu(sub 4), nu(sub 3) + nu(sub 4), and nu(sub 2) + nu(sub 3) bands of (12)CH4 located between 4118 and 4615/cm. These results were obtained by analyzing a total of 29 spectra simultaneously using a multi-spectral non-linear least-squares fitting technique. This new technique allowed the determination of correlated spectral line parameters (e.g. intensity and broadening coefficient) better than the procedure of averaging values obtained by fitting the spectra individually. This method also provided a direct determination of the uncertainties in the retrieved parameters due to random errors. For each band analysed in this study the dependence of the various spectral line parameters upon the tetrahedral symmetry species and the rotational quantum numbers of the transitions is also presented.

Description: Rotor noise prediction codes predict the thickness and loading noise produced by a helicopter rotor, given the blade motion, rotor operating conditions, and fluctuating force distribution over the blade surface. However, the criticality of these various inputs, and their respective effects on the predicted acoustic field, have never been fully addressed. This paper examines the importance of these inputs, and the sensitivity of the acoustic predicitions to a variation of each parameter. The effects of collective and cyclic pitch, as well as coning and cyclic flapping, are presented. Blade loading inputs are examined to determine the necessary spatial and temporal resolution, as well as the importance of the chordwise distribution. The acoustic predictions show regions in the acoustic field where significant errors occur when simplified blade motions or blade loadings are used. An assessment of the variation in the predicted acoustic field is balanced by a consideration of Central Processing Unit (CPU) time necessary for the various approximations.

Description: The U.S. National Aeronautics and Space Administration (NASA) Balloon Program has been highly successful since recovering from the catastrophic balloon failure problems of the early to mid 1980s. Balloons have continued to perform at unprecedented success rates. The comprehensive research and development (R&D) effort has continued with advances being made across the spectrum of balloon related disciplines. The long duration balloon project will be transitioning from a development effort to an operational capability this year. Recently, emphasis has been placed on the development and implementation of new support systems and facilities. A new permanent launch facility at Fort Sumner, New Mexico has been established. New ground station support equipment is being implemented, and a new heavy load launch vehicle is scheduled to be implemented in 1992. The progress, status and future plans for these and other aspects of the NASA program will be presented.

Description: The catastrophic balloon failure during the first half of the 1980's identified the need for a comprehensive and continuing balloon research and development (R&D) commitment by NASA. Technical understanding was lacking in many of the disciplines and processes associated with scientific ballooning. A comprehensive balloon R&D plan was developed in 1986 and implemented in 1987. The objectives were to develop the understanding of balloon system performance, limitations, and failure mechanisms. The program consisted of five major technical areas: structures, performance and analysis, materials, chemistry and processing, and quality control. Research activitites have been conducted at NASA/Goddard Space Flight Center (GSFC)-Wallops Flight Facility (WFF), other NASA centers and government facilities, universities, and the balloon manufacturers. Several new and increased capabilities and resources have resulted from this activity. The findings, capabilities, and plan of the balloon R&D program are presented.

Description: Caps have been used to structurally reinforce scientific research balloons since the late 1950's. The scientific research balloons used by the National Aeronautics and Space Administration (NASA) use internal caps. A NASA cap placement specification does not exist since no empirical information exisits concerning cap placement. To develop a cap placement specification, NASA has completed two in-hangar inflation tests comparing the structural contributions of internal caps and external caps. The tests used small scale test balloons designed to develop the highest possible stresses within the constraints of the hangar and balloon materials. An externally capped test balloon and an internally capped test balloon were designed, built, inflated and simulated to determine the structural contributions and benefits of each. The results of the tests and simulations are presented.

Description: The purpose of this Note is to present results from an analytic/experimental study that investigated the potential for passively changing blade twist through the use of extension-twist coupling. A set of composite model rotor blades was manufactured from existing blade molds for a low-twist metal helicopter rotor blade, with a view toward establishing a preliminary proof concept for extension-twist-coupled rotor blades. Data were obtained in hover for both a ballasted and unballasted blade configuration in sea-level atmospheric conditions. Test data were compared with results obtained from a geometrically nonlinear analysis of a detailed finite element model of the rotor blade developed in MSC/NASTRAN.

Description: The paper considers the compressible Rayleigh equation as a model for the Mach wave emission mechanism associated with high-temperature supersonic jets. Solutions to the compressible Rayleigh equation reveal the existence of several families of supersonically convecting instability waves. These waves directly radiate noise to the jet far field. The predicted noise characteristics are compared to previously acquired experimental data for an axisymmetric Mach 2 fully pressure balanced jet operating over a range of jet total temperatures from ambient to 1370 K. The results of this comparison show that the first-order supersonic instability wave and the Kelvin-Hemlhlotz first-, second-, and third-order modes have directional radiation characteristics that are in agreement with observed data. The assumption of equal initial amplitudes for all of the waves leads to the conclusion that the flapping mode of instability dominates the noise radiatio process of supersonic jets. At a jet temperature of 1370 K, supersonic instability waves are predicted to dominate the noise radiated at high frequency at narrow angles to the jet axis.

Description: The objective of the present work is to study the mixing characteristics of a linear array of supersonic rectangular jets under conditions of screech synchronization. The screech synchronization at a fully expanded jet Mach number of 1.61 is achieved by a precise adjustment of the internozzle spacing. To our knowledge, such an experiment on the resonant mixing of screech synchronized multiple rectangular jets has not been reported before. The results are compared with the case where the screech was suppressed in the multijet configuration.

Description: The objective of the present investigation is to assess the effect of the spatial order of accuracy used for the evaluation of the inviscid fluxes on the resolution of higher order quantitites, such as velocity gradients. The viscous terms are computed as second-order accurate with central difference formulas, even though for the explicit part of the algorithm higher order approximations may be used. A viscous/inviscid method is used, and the outer part of the flowfield is computed with the inviscid flow equations. The viscous boundary-layer type flow region close to the body surface is computed with an algebraic eddy viscosity model. Results obtained with the conservative and nonconservative formulations and the viscous/inviscid approach are compared with available experimental data. The effect of grid refinement on the accuracy of the solution is also presented.

Description: The benefits of using a hypersonic waverider for spacecraft trajectory modification are presented. A waverider is a hypersonic vehicle specifically designed so that the undersurface bow shock is attached to the leading edge, which provides for the highest known lift-to-drag ratios achievable at high Mach number flight. Several viable space missions are suggested which could use such configurations for low-drag aero-assisted maneuvers in planetary atmospheres. It is shown that large changes in the spacecraft velocity vector can be accomplished with acceptably small losses in energy due to drag using a waverider aeroshell. The primary advantage of an aero-assist maneuver is suggested by comparison to a traditional gravity-assist trajectory. Some scaling laws are presented for comparing waveriders designed for different planetary atmospheres, and it is shown that the compositional differences between the terrestrial planets has a minimal impact on waverider design.

Description: This paper presents the results of a transient computer simulation that was developed to study phase change energy storage techniques for Space Station Freedom (SSF) solar dynamic (SD) power systems. Such SD systems may be used in future growth SSF configurations. Two solar dynamic options are considered in this paper: Brayton and Rankine. Model elements consist of a single node receiver and concentrator, and takes into account overall heat engine efficiency and power distribution characteristics. The simulation not only computes the energy stored in the receiver phase change material (PCM), but also the amount of the PCM required for various combinations of load demands and power system mission constraints. For a solar dynamic power system in low earth orbit, the amount of stored PCM energy is calculated by balancing the solar energy input and the energy consumed by the loads corrected by an overall system efficiency. The model assumes an average 75 kW SD power system load profile which is connected to user loads via dedicated power distribution channels. The model then calculates the stored energy in the receiver and subsequently estimates the quantity of PCM necessary to meet peaking and contingency requirements. The model can also be used to conduct trade studies on the performance of SD power systems using different storage materials.

Description: ARC2D is a computational fluid dynamics program developed at the NASA Ames Research Center specifically for airfoil computations. The program uses implicit finite-difference techniques to solve two-dimensional Euler equations and thin layer Navier-Stokes equations. It is based on the Beam and Warming implicit approximate factorization algorithm in generalized coordinates. The methods are either time accurate or accelerated non-time accurate steady state schemes. The evolution of the solution through time is physically realistic; good solution accuracy is dependent on mesh spacing and boundary conditions. The mathematical development of ARC2D begins with the strong conservation law form of the two-dimensional Navier-Stokes equations in Cartesian coordinates, which admits shock capturing. The Navier-Stokes equations can be transformed from Cartesian coordinates to generalized curvilinear coordinates in a manner that permits one computational code to serve a wide variety of physical geometries and grid systems. ARC2D includes an algebraic mixing length model to approximate the effect of turbulence. In cases of high Reynolds number viscous flows, thin layer approximation can be applied. ARC2D allows for a variety of solutions to stability boundaries, such as those encountered in flows with shocks. The user has considerable flexibility in assigning geometry and developing grid patterns, as well as in assigning boundary conditions. However, the ARC2D model is most appropriate for attached and mildly separated boundary layers; no attempt is made to model wake regions and widely separated flows. The techniques have been successfully used for a variety of inviscid and viscous flowfield calculations. The Cray version of ARC2D is written in FORTRAN 77 for use on Cray series computers and requires approximately 5Mb memory. The program is fully vectorized. The tape includes variations for the COS and UNICOS operating systems. Also included is a sample routine for CONVEX computers to emulate Cray system time calls, which should be easy to modify for other machines as well. The standard distribution media for this version is a 9-track 1600 BPI ASCII Card Image format magnetic tape. The Cray version was developed in 1987. The IBM ES/3090 version is an IBM port of the Cray version. It is written in IBM VS FORTRAN and has the capability of executing in both vector and parallel modes on the MVS/XA operating system and in vector mode on the VM/XA operating system. Various options of the IBM VS FORTRAN compiler provide new features for the ES/3090 version, including 64-bit arithmetic and up to 2 GB of virtual addressability. The IBM ES/3090 version is available only as a 9-track, 1600 BPI IBM IEBCOPY format magnetic tape. The IBM ES/3090 version was developed in 1989. The DEC RISC ULTRIX version is a DEC port of the Cray version. It is written in FORTRAN 77 for RISC-based Digital Equipment platforms. The memory requirement is approximately 7Mb of main memory. It is available in UNIX tar format on TK50 tape cartridge. The port to DEC RISC ULTRIX was done in 1990. COS and UNICOS are trademarks and Cray is a registered trademark of Cray Research, Inc. IBM, ES/3090, VS FORTRAN, MVS/XA, and VM/XA are registered trademarks of International Business Machines. DEC and ULTRIX are registered trademarks of Digital Equipment Corporation.

Description: Panel method computer programs are software tools of moderate cost used for solving a wide range of engineering problems. The panel code PMARC_12 (Panel Method Ames Research Center, version 12) can compute the potential flow field around complex three-dimensional bodies such as complete aircraft models. PMARC_12 is a well-documented, highly structured code with an open architecture that facilitates modifications and the addition of new features. Adjustable arrays are used throughout the code, with dimensioning controlled by a set of parameter statements contained in an include file; thus, the size of the code (i.e. the number of panels that it can handle) can be changed very quickly. This allows the user to tailor PMARC_12 to specific problems and computer hardware constraints. In addition, PMARC_12 can be configured (through one of the parameter statements in the include file) so that the code's iterative matrix solver is run entirely in RAM, rather than reading a large matrix from disk at each iteration. This significantly increases the execution speed of the code, but it requires a large amount of RAM memory. PMARC_12 contains several advanced features, including internal flow modeling, a time-stepping wake model for simulating either steady or unsteady (including oscillatory) motions, a Trefftz plane induced drag computation, off-body and on-body streamline computations, and computation of boundary layer parameters using a two-dimensional integral boundary layer method along surface streamlines. In a panel method, the surface of the body over which the flow field is to be computed is represented by a set of panels. Singularities are distributed on the panels to perturb the flow field around the body surfaces. PMARC_12 uses constant strength source and doublet distributions over each panel, thus making it a low order panel method. Higher order panel methods allow the singularity strength to vary linearly or quadratically across each panel. Experience has shown that low order panel methods can provide nearly the same accuracy as higher order methods over a wide range of cases with significantly reduced computation times; hence, the low order formulation was adopted for PMARC_12. The flow problem is solved by modeling the body as a closed surface dividing space into two regions: the region external to the surface in which an unknown velocity potential exists representing the flow field of interest, and the region internal to the surface in which a known velocity potential (representing a fictitious flow) is prescribed as a boundary condition. Both velocity potentials are required to satisfy Laplace's equation. A surface integral equation for the unknown potential external to the surface can be written by applying Green's Theorem to the external region. Using the internal potential and zero flow through the surface as boundary conditions, the unknown potential external to the surface can be solved for. When the internal flow option, which allows the analysis of closed ducts, wind tunnels, and similar internal flow problems, is selected, the geometry is modeled such that the flow field of interest is inside the geometry and the fictitious flow is outside the geometry. Items such as wings, struts, or aircraft models can be included in the internal flow problem. The time-stepping wake model gives PMARC_12 the ability to model both steady and unsteady flow problems. The wake is convected downstream from the wake-separation line by the local velocity field. With each time step, a new row of wake panels is added to the wake at the wake-separation line. Time stepping can start from time t=0 (no initial wake) or from time t=t0 (an initial wake is specified). A wide range of motions can be prescribed, including constant rates of translation, constant rate of rotation about an arbitrary axis, oscillatory translation, and oscillatory rotation about any of the three coordinate axes. Investigators interested in a visual representation of the phenomenon they are studying with PMARC_12 may want to consider obtaining the program GVS (ARC-13361), the General Visualization System. GVS is a Silicon Graphics IRIS program which was created for the purpose of supporting the scientific visualization needs of PMARC_12. GVS is available separately from COSMIC. PMARC_12 is written in standard FORTRAN 77, with the exception of the NAMELIST extension used for input. This makes the code fairly machine independent. A compiler which supports the NAMELIST extension is required. The amount of free disk space and RAM memory required for PMARC_12 will vary depending on how the code is dimensioned using the parameter statements in the include file. The recommended minimum requirements are 20Mb of free disk space and 4Mb of RAM. PMARC_12 has been successfully implemented on a Macintosh II running System 6.0.7 or 7.0 (using MPW/Language Systems Fortran 3.0), a Sun SLC running SunOS 4.1.1, an HP 720 running HP-UX 8.07, an SGI IRIS running IRIX 4.0 (it will not run under IRIX 3.x.x without modifications), an IBM RS/6000 running AIX, a DECstation 3100 running ULTRIX, and a CRAY-YMP running UNICOS 6.0 or later. Due to its memory requirements, this program does not readily lend itself to implementation on MS-DOS based machines. The standard distribution medium for PMARC_12 is a set of three 3.5 inch 800K Macintosh format diskettes and one 3.5 inch 1.44Mb Macintosh format diskette which contains an electronic copy of the documentation in MS Word 5.0 format for the Macintosh. Alternate distribution media and formats are available upon request, but these will not include the electronic version of the document. No executables are included on the distribution media. This program is an update to PMARC version 11, which was released in 1989. PMARC_12 was released in 1993. It is available only for use by United States citizens.

Description: There is no simple and perfect way to measure residual stresses in metal parts that have been welded or deformed to make complex structures such as pressure vessels and aircraft, yet these locked-in stresses can contribute to structural failure by fatigue and fracture. However, one proven and tested technique for determining the internal stress of a metal part is to drill a test hole while measuring the relieved strains around the hole, such as the hole-drilling strain gage method described in ASTM E 837. The program HOLEGAGE processes strain gage data and provides additional calculations of internal stress variations that are not obtained with standard E 837 analysis methods. The typical application of the technique uses a three gage rosette with a special hole-drilling fixture for drilling a hole through the center of the rosette to produce a hole with very small gage pattern eccentricity error. Another device is used to control the drilling and halt the drill at controlled depth steps. At each step, strains from all three strain gages are recorded. The influence coefficients used by HOLEGAGE to compute stresses from relieved hole strains were developed by published finite element method studies of thick plates for specific hole sizes and depths. The program uses a parabolic fit and an interpolating scheme to project the coefficients to other hole sizes and depths. Additionally, published experimental data are used to extend the coefficients to relatively thin plates. These influence coefficients are used to compute the stresses in the original part from the strain data. HOLEGAGE will compute interior planar stresses using strain data from each drilled hole depth layer. Planar stresses may be computed in three ways including: a least squares fit for a linear variation with depth, an integral method to give incremental stress data for each layer, or by a linear fit to the integral data (with some surface data points omitted) to predict surface stresses before strain gage sanding preparations introduced additional residual stresses. Options are included for estimating the effect of hole eccentricity on calculations, smoothing noise from the strain data, and inputting the program data either interactively or from a data file. HOLEGAGE was written in FORTRAN 77 for DEC VAX computers under VMS, and is transportable except for system-unique TIME and DATE system calls. The program requires 54K of main memory and was developed in 1990. The program is available on a 9-track 1600 BPI VAX BACKUP format magnetic tape (standard media) or a TK50 tape cartridge. The documentation is included on the tape. DEC VAX and VMS are trademarks of Digital Equipment Corporation.

Description: This program determines the supersonic flowfield surrounding three-dimensional wing-body configurations of a delta wing. It was designed to provide the numerical computation of three dimensional inviscid, flowfields of either perfect or real gases about supersonic or hypersonic airplanes. The governing equations in conservation law form are solved by a finite difference method using a second order noncentered algorithm between the body and the outermost shock wave, which is treated as a sharp discontinuity. Secondary shocks which form between these boundaries are captured automatically. The flowfield between the body and outermost shock is treated in a shock capturing fashion and therefore allows for the correct formation of secondary internal shocks . The program operates in batch mode, is in CDC update format, has been implemented on the CDC 7600, and requires more than 140K (octal) word locations.

Description: The Solar Space Power Analysis Code, SOSPAC, was developed to examine the solar thermal and photovoltaic power generation options available for a satellite or spacecraft in low earth orbit. SOSPAC is a preliminary systems analysis tool and enables the engineer to compare the areas, weights, and costs of several candidate electric and thermal power systems. The configurations studied include photovoltaic arrays and parabolic dish systems to produce electricity only, and in various combinations to provide both thermal and electric power. SOSPAC has been used for comparison and parametric studies of proposed power systems for the NASA Space Station. The initial requirements are projected to be about 40 kW of electrical power, and a similar amount of thermal power with temperatures above 1000 degrees Centigrade. For objects in low earth orbit, the aerodynamic drag caused by suitably large photovoltaic arrays is very substantial. Smaller parabolic dishes can provide thermal energy at a collection efficiency of about 80%, but at increased cost. SOSPAC allows an analysis of cost and performance factors of five hybrid power generating systems. Input includes electrical and thermal power requirements, sun and shade durations for the satellite, and unit weight and cost for subsystems and components. Performance equations of the five configurations are derived, and the output tabulates total weights of the power plant assemblies, area of the arrays, efficiencies, and costs. SOSPAC is written in FORTRAN IV for batch execution and has been implemented on an IBM PC computer operating under DOS with a central memory requirement of approximately 60K of 8 bit bytes. This program was developed in 1985.

Description: The PYROLASER package is an operating system for the Pyrometer Instrument Company's Pyrolaser. There are 6 individual programs in the PYROLASER package: two main programs, two lower level subprograms, and two programs which, although independent, function predominantly as macros. The package provides a quick and easy way to setup, control, and program a standard Pyrolaser. Temperature and emissivity measurements may be either collected as if the Pyrolaser were in the manual operations mode, or displayed on real time strip charts and stored in standard spreadsheet format for post-test analysis. A shell is supplied to allow macros, which are test-specific, to be easily added to the system. The Pyrolaser Simple Operation program provides full on-screen remote operation capabilities, thus allowing the user to operate the Pyrolaser from the computer just as it would be operated manually. The Pyrolaser Simple Operation program also allows the use of "quick starts". Quick starts provide an easy way to permit routines to be used as setup macros for specific applications or tests. The specific procedures required for a test may be ordered in a sequence structure and then the sequence structure can be started with a simple button in the cluster structure provided. One quick start macro is provided for continuous Pyrolaser operation. A subprogram, Display Continuous Pyr Data, is used to display and store the resulting data output. Using this macro, the system is set up for continuous operation and the subprogram is called to display the data in real time on strip charts. The data is simultaneously stored in a spreadsheet format. The resulting spreadsheet file can be opened in any one of a number of commercially available spreadsheet programs. The Read Continuous Pyrometer program is provided as a continuously run subprogram for incorporation of the Pyrolaser software into a process control or feedback control scheme in a multi-component system. The program requires the Pyrolaser to be set up using the Pyrometer String Transfer macro. It requires no inputs and provides temperature and emissivity as outputs. The Read Continuous Pyrometer program can be run continuously and the data can be sampled as often or as seldom as updates of temperature and emissivity are required. PYROLASER is written using the Labview software for use on Macintosh series computers running System 6.0.3 or later, Sun Sparc series computers running OpenWindows 3.0 or MIT's X Window System (X11R4 or X11R5), and IBM PC or compatibles running Microsoft Windows 3.1 or later. Labview requires a minimum of 5Mb of RAM on a Macintosh, 24Mb of RAM on a Sun, and 8Mb of RAM on an IBM PC or compatible. The Labview software is a product of National Instruments (Austin,TX; 800-433-3488), and is not included with this program. The standard distribution medium for PYROLASER is a 3.5 inch 800K Macintosh format diskette. It is also available on a 3.5 inch 720K MS-DOS format diskette, a 3.5 inch diskette in UNIX tar format, and a .25 inch streaming magnetic tape cartridge in UNIX tar format. An electronic copy of the documentation in Macintosh WordPerfect version 2.0.4 format is included on the distribution medium. Printed documentation is included in the price of the program. PYROLASER was developed in 1992.

Description: This theoretical aerodynamics program, TAD, was developed to predict the aerodynamic characteristics of vehicles with sounding rocket configurations. These slender, axisymmetric finned vehicle configurations have a wide range of aeronautical applications from rockets to high speed armament. Over a given range of Mach numbers, TAD will compute the normal force coefficient derivative, the center-of-pressure, the roll forcing moment coefficient derivative, the roll damping moment coefficient derivative, and the pitch damping moment coefficient derivative of a sounding rocket configured vehicle. The vehicle may consist of a sharp pointed nose of cone or tangent ogive shape, up to nine other body divisions of conical shoulder, conical boattail, or circular cylinder shape, and fins of trapezoid planform shape with constant cross section and either three or four fins per fin set. The characteristics computed by TAD have been shown to be accurate to within ten percent of experimental data in the supersonic region. The TAD program calculates the characteristics of separate portions of the vehicle, calculates the interference between separate portions of the vehicle, and then combines the results to form a total vehicle solution. Also, TAD can be used to calculate the characteristics of the body or fins separately as an aid in the design process. Input to the TAD program consists of simple descriptions of the body and fin geometries and the Mach range of interest. Output includes the aerodynamic characteristics of the total vehicle, or user-selected portions, at specified points over the mach range. The TAD program is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 123K of 8 bit bytes. The TAD program was originally developed in 1967 and last updated in 1972.

Description: Slush hydrogen, a mixture of the solid and liquid phases of hydrogen, is a possible source of fuel for the National Aerospace Plane (NASP) Project. Advantages of slush hydrogen over liquid hydrogen include greater heat capacity and greater density. However, practical use of slush hydrogen as a fuel requires systems of lines, valves, etc. which are designed to deliver the fuel in slush form with minimal solid loss as a result of pipe heating or flow friction. Engineers involved with the NASP Project developed FLUSH to calculate the pressure drop and slush hydrogen solid fraction loss for steady-state, one-dimensional flow. FLUSH solves the steady-state, one-dimensional energy equation and the Bernoulli equation for pipe flow. The program performs these calculations for each two-node element--straight pipe length, elbow, valve, fitting, or other part of the piping system--specified by the user. The user provides flow rate, upstream pressure, initial solid hydrogen fraction, element heat leak, and element parameters such as length and diameter. For each element, FLUSH first calculates the pressure drop, then figures the slush solid fraction exiting the element. The code employs GASPLUS routines to calculate thermodynamic properties for the slush hydrogen. FLUSH is written in FORTRAN IV for DEC VAX series computers running VMS. An executable is provided on the tape. The GASPLUS physical properties routines which are required for building the executable are included as one object library on the program media (full source code for GASPLUS is available separately as COSMIC Program Number LEW-15091). FLUSH is available in DEC VAX BACKUP format on a 9-track 1600 BPI magnetic tape (standard media) or on a TK50 tape cartridge. FLUSH was developed in 1989.

Description: Scale-space filtering is used to screen information obtained from signals that produce a complex curve (such as that in geographic and thermal analysis) to gain a truer representation of the area under analysis. PSF extends this technique to extract non-periodic hills and valleys from a signal. Because the signal's information is sometimes too complex to determine with certainty if some features are real or artificial, PSF calculates probabilities, with the extracted features corresponding to real events, in order to aid in determining the signal's accuracy. Since the probabilities associated with the features are derived from domain-specific statistics, it is (most likely) necessary to modify the program code to correspond to the user's particular domain. PSF also provides a standard scale-space filtering algorithm for use when the desired features can be identified with certainty or when it is not practical to get the domain-specific statistics. The PSF algorithm is based on Witkin's scale-space filtering theory. The program detects signal variations by finding the points of inflection in the input signal. The number and position of these points are dependent upon the scale of the derivative operators used to detect them. Therefore, instead of assuming any single scale to be correct, PSF identifies points of inflection in a large number of different scales. It then describes the curve according to the groups of points of inflection, across all scales, caused by the same physical process. PSF provides an output table giving the following information: the abscissa of the first inflection of the peak, the type of peak, the distance between the first and second inflection points, the abscissa of the peak, and the probability of the feature corresponding to a real event in the curve. The program will also list points representing a graphical image of the signal and detected peaks. This data can be used with a standard plotting program (not included) to display the signal and its features graphically. PSF is written in C language (49%) and Common LISP (51%) for use on a Sun SPARC workstation running the UNIX operating system. PSF requires 4Mb of RAM. The standard distribution medium for this program is a .25 streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format. PSF was developed in 1991. Sun and SPARC are trademarks of Sun Microsystems, Inc. UNIX is a registered trademark of AT&T Bell Laboratories.

Description: This program, which is called 'AOFA', determines the complete viscous and inviscid flow around a body of revolution at a given angle of attack and traveling at supersonic speeds. The viscous calculations from this program agree with experimental values for surface and pitot pressures and with surface heating rates. At high speeds, lee-side flows are important because the local heating is difficult to correlate and because the shed vortices can interact with vehicle components such as a canopy or a vertical tail. This program should find application in the design analysis of any high speed vehicle. Lee-side flows are difficult to calculate because thin-boundary-layer theory is not applicable and the concept of matching inviscid and viscous flow is questionable. This program uses the parabolic approximation to the compressible Navier-Stokes equations and solves for the complete inviscid and viscous regions of flow, including the pressure. The parabolic approximation results from the assumption that the stress derivatives in the streamwise direction are small in comparison with derivatives in the normal and circumferential directions. This assumption permits the equation to be solved by an implicit finite difference marching technique which proceeds downstream from the initial data point, provided the inviscid portion of flow is supersonic. The viscous cross-flow separation is also determined as part of the solution. To use this method it is necessary to first determine an initial data point in a region where the inviscid portion of the flow is supersonic. Input to this program consists of two parts. Problem description is conveyed to the program by namelist input. Initial data is acquired by the program as formatted data. Because of the large amount of run time this program can consume the program includes a restart capability. Output is in printed format and magnetic tape for further processing. This program is written in FORTRAN IV and has been implemented on a CDC 7600 with a central memory requirement of approximately 35K (octal) of 60 bit words.

Description: The complex environment of the typical research laboratory requires flexible process control. This program provides natural language process control from an IBM PC or compatible machine. Sometimes process control schedules require changes frequently, even several times per day. These changes may include adding, deleting, and rearranging steps in a process. This program sets up a process control system that can either run without an operator, or be run by workers with limited programming skills. The software system includes three programs. Two of the programs, written in FORTRAN77, record data and control research processes. The third program, written in Pascal, generates the FORTRAN subroutines used by the other two programs to identify the user commands with the user-written device drivers. The software system also includes an input data set which allows the user to define the user commands which are to be executed by the computer. To set the system up the operator writes device driver routines for all of the controlled devices. Once set up, this system requires only an input file containing natural language command lines which tell the system what to do and when to do it. The operator can make up custom commands for operating and taking data from external research equipment at any time of the day or night without the operator in attendance. This process control system requires a personal computer operating under MS-DOS with suitable hardware interfaces to all controlled devices. The program requires a FORTRAN77 compiler and user-written device drivers. This program was developed in 1989 and has a memory requirement of about 62 Kbytes.

Description: Accurate simulation of nuclear thermal propulsion systems using computational methods will permit reductions in testing and, thus, the time and cost of achieving a flight ready status for systems utilizing this advanced technology. An accurate simulation must maintain a "balance-of-plant" where the required pump work equals the supplied turbine work. This turbopump assembly balancing must be integrated into the overall system analysis models. TPA was developed to balance turbine and pump work using performance maps. It requires the inlet properties, performance maps, and shaft speed. TPA then computes the exit conditions and work terms. The work terms can then be balanced by varying the input shaft speed. The objective of the pump analysis is to determine the propellant state properties at the pump exit and the pump work. The pump analysis algorithm for liquid flow assumes that the shaft speed, the propellant state properties at the pump entrance, the propellant flow rate, the pump entrance and exit areas, as well as performance curves, are all known. The analysis of both the pump pressure rise and pump efficiency curves is required. The objective of the turbine analysis is to determine the propellant state properties at the turbine exit and the turbine work. The turbine analysis algorithm assumes that the shaft speed, the propellant state properties at the turbine entrance, the propellant flow rate, the turbine root mean square blade diameter, the turbine entrance and exit areas, as well as performance curves, are all known. The analysis also requires the turbine flow parameter curve and the turbine total efficiency curve. TPA is written in FORTRAN 77 to be machine independent. The TPA package includes the NBS+_PH2 code, which is also available separately (LEW-15505). TPA has been successfully implemented on a DEC VAX series computer running VMS, a Sun4 series computer running SunOS, and an IBM PC compatible computer running MS-DOS. Lahey F77L3 EM/32 v. 5.01 or higher is required for compilation on an IBM PC compatible computer; however, a PC executable is included on the distribution diskette. The standard distribution medium for this program is one 5.25 inch 360K MS-DOS format diskette. The diskette's contents have been compressed using PKWARE's archiving tools. The utility to unarchive the file, PKUNZIP.EXE, is included. Alternate distribution media and formats are available upon request. TPA was developed in 1993.

Description: The Comprehensive Analytical Model of Rotorcraft Aerodynamics, CAMRAD, program is designed to calculate rotor performance, loads, and noise; helicopter vibration and gust response; flight dynamics and handling qualities; and system aeroelastic stability. The analysis is a consistent combination of structural, inertial, and aerodynamic models applicable to a wide range of problems and a wide class of vehicles. The CAMRAD analysis can be applied to articulated, hingeless, gimballed, and teetering rotors with an arbitrary number of blades. The rotor degrees of freedom included are blade/flap bending, rigid pitch and elastic torsion, and optionally gimbal or teeter motion. General two-rotor aircrafts can be modeled. Single main-rotor and tandem helicopter and sideby-side tilting proprotor aircraft configurations can be considered. The case of a rotor or helicopter in a wind tunnel can also be modeled. The aircraft degrees of freedom included are the six rigid body motion, elastic airframe motions, and the rotor/engine speed perturbations. CAMRAD calculates the load and motion of helicopters and airframes in two stages. First the trim solution is obtained; then the flutter, flight dynamics, and/or transient behavior can be calculated. The trim operating conditions considered include level flight, steady climb or descent, and steady turns. The analysis of the rotor includes nonlinear inertial and aerodynamic models, applicable to large blade angles and a high inflow ratio, The rotor aerodynamic model is based on two-dimensional steady airfoil characteristics with corrections for three-dimensional and unsteady flow effects, including a dynamic stall model. In the flutter analysis, the matrices are constructed that describe the linear differential equations of motion, and the equations are analyzed. In the flight dynamics analysis, the stability derivatives are calculated and the matrices are constructed that describe the linear differential equations of motion. These equations are analyzed. In the transient analysis, the rigid body equations of motion are numerically integrated, for a prescribed transient gust or control input. The CAMRAD program product is available by license for a period of ten years to domestic U.S. licensees. The licensed program product includes the CAMRAD source code, command procedures, sample applications, and one set of supporting documentation. Copies of the documentation may be purchased separately at the price indicated below. CAMRAD is written in FORTRAN 77 for the DEC VAX under VMS 4.6 with a recommended core memory of 4.04 megabytes. The DISSPLA package is necessary for graphical output. CAMRAD was developed in 1980.

Description: ARC2D is a computational fluid dynamics program developed at the NASA Ames Research Center specifically for airfoil computations. The program uses implicit finite-difference techniques to solve two-dimensional Euler equations and thin layer Navier-Stokes equations. It is based on the Beam and Warming implicit approximate factorization algorithm in generalized coordinates. The methods are either time accurate or accelerated non-time accurate steady state schemes. The evolution of the solution through time is physically realistic; good solution accuracy is dependent on mesh spacing and boundary conditions. The mathematical development of ARC2D begins with the strong conservation law form of the two-dimensional Navier-Stokes equations in Cartesian coordinates, which admits shock capturing. The Navier-Stokes equations can be transformed from Cartesian coordinates to generalized curvilinear coordinates in a manner that permits one computational code to serve a wide variety of physical geometries and grid systems. ARC2D includes an algebraic mixing length model to approximate the effect of turbulence. In cases of high Reynolds number viscous flows, thin layer approximation can be applied. ARC2D allows for a variety of solutions to stability boundaries, such as those encountered in flows with shocks. The user has considerable flexibility in assigning geometry and developing grid patterns, as well as in assigning boundary conditions. However, the ARC2D model is most appropriate for attached and mildly separated boundary layers; no attempt is made to model wake regions and widely separated flows. The techniques have been successfully used for a variety of inviscid and viscous flowfield calculations. The Cray version of ARC2D is written in FORTRAN 77 for use on Cray series computers and requires approximately 5Mb memory. The program is fully vectorized. The tape includes variations for the COS and UNICOS operating systems. Also included is a sample routine for CONVEX computers to emulate Cray system time calls, which should be easy to modify for other machines as well. The standard distribution media for this version is a 9-track 1600 BPI ASCII Card Image format magnetic tape. The Cray version was developed in 1987. The IBM ES/3090 version is an IBM port of the Cray version. It is written in IBM VS FORTRAN and has the capability of executing in both vector and parallel modes on the MVS/XA operating system and in vector mode on the VM/XA operating system. Various options of the IBM VS FORTRAN compiler provide new features for the ES/3090 version, including 64-bit arithmetic and up to 2 GB of virtual addressability. The IBM ES/3090 version is available only as a 9-track, 1600 BPI IBM IEBCOPY format magnetic tape. The IBM ES/3090 version was developed in 1989. The DEC RISC ULTRIX version is a DEC port of the Cray version. It is written in FORTRAN 77 for RISC-based Digital Equipment platforms. The memory requirement is approximately 7Mb of main memory. It is available in UNIX tar format on TK50 tape cartridge. The port to DEC RISC ULTRIX was done in 1990. COS and UNICOS are trademarks and Cray is a registered trademark of Cray Research, Inc. IBM, ES/3090, VS FORTRAN, MVS/XA, and VM/XA are registered trademarks of International Business Machines. DEC and ULTRIX are registered trademarks of Digital Equipment Corporation.

Description: Panel methods are moderate cost tools for solving a wide range of engineering problems. PMARC (Panel Method Ames Research Center) is a potential flow panel code that numerically predicts flow fields around complex three-dimensional geometries. PMARC's predecessor was a panel code named VSAERO which was developed for NASA by Analytical Methods, Inc. PMARC is a new program with many additional subroutines and a well-documented code suitable for powered-lift aerodynamic predictions. The program's open architecture facilitates modifications or additions of new features. Another improvement is the adjustable size code which allows for an optimum match between the computer hardware available to the user and the size of the problem being solved. PMARC can be resized (the maximum number of panels can be changed) in a matter of minutes. Several other state-of-the-art PMARC features include internal flow modeling for ducts and wind tunnel test sections, simple jet plume modeling essential for the analysis and design of powered-lift aircraft, and a time-stepping wake model which allows the study of both steady and unsteady motions. PMARC is a low-order panel method, which means the singularities are distributed with constant strength over each panel. In many cases low-order methods can provide nearly the same accuracy as higher order methods (where the singularities are allowed to vary linearly or quadratically over each panel). Low-order methods have the advantage of a shorter computation time and do not require exact matching between panels. The flow problem is solved by assuming that the body is at rest in a moving flow field. The body is modeled as a closed surface which divides space into two regions -- one region contains the flow field of interest and the other contains a fictitious flow. External flow problems, such as a wing in a uniform stream, have the external region as the flow field of interest and the internal flow as the fictitious flow. This arrangement is reversed for internal flow problems where the internal region contains the flow field of interest and the external flow field is fictitious. In either case it is assumed that the velocity potentials in both regions satisfy Laplace's equation. PMARC has extensive geometry modeling capabilities for handling complex, three-dimensional surfaces. As with all panel methods, the geometry must be modeled by a set of panels. For convenience, the geometry is usually subdivided into several pieces and modeled with sets of panels called patches. A patch may be folded over on itself so that opposing sides of the patch form a common line. For example, wings are normally modeled with a folded patch to form the trailing edge of the wing. PMARC also has the capability to automatically generate a closing tip patch. In the case of a wing, a tip patch could be generated to close off the wing's third side. PMARC has a simple jet model for simulating a jet plume in a crossflow. The jet plume shape, trajectory, and entrainment velocities are computed using the Adler/Baron jet in crossflow code. This information is then passed back to PMARC. The wake model in PMARC is a time-stepping wake model. The wake is convected downstream from the wake separation line by the local velocity flowfield. With each time step, a new row of wake panels is added to the wake at the wake separation line. PMARC also allows an initial wake to be specified if desired, or, as a third option, no wakes need be modeled. The effective presentation of results for aerodynamics problems requires the generation of report-quality graphics. PMAPP (ARC-12751), the Panel Method Aerodynamic Plotting Program, (Sterling Software), was written for scientists at NASA's Ames Research Center to plot the aerodynamic analysis results (flow data) from PMARC. PMAPP is an interactive, color-capable graphics program for the DEC VAX or MicroVAX running VMS. It was designed to work with a variety of terminal types and hardcopy devices. PMAPP is available separately from COSMIC. PMARC was written in standard FORTRAN77 using adjustable size arrays throughout the code. Redimensioning PMARC will change the amount of disk space and memory the code requires to be able to run; however, due to its memory requirements, this program does not readily lend itself to implementation on MS-DOS based machines. The program was implemented on an Apple Macintosh (using 2.5 MB of memory) and tested on a VAX/VMS computer. The program is available on a 3.5 inch Macintosh format diskette (standard media) or in VAX BACKUP format on TK50 tape cartridge or 9-track magnetic tape. PMARC was developed in 1989.

Description: The Transonic Airfoil analysis computer code, TAIR, was developed to employ a fast, fully implicit algorithm to solve the conservative full-potential equation for the steady transonic flow field about an arbitrary airfoil immersed in a subsonic free stream. The full-potential formulation is considered exact under the assumptions of irrotational, isentropic, and inviscid flow. These assumptions are valid for a wide range of practical transonic flows typical of modern aircraft cruise conditions. The primary features of TAIR include: a new fully implicit iteration scheme which is typically many times faster than classical successive line overrelaxation algorithms; a new, reliable artifical density spatial differencing scheme treating the conservative form of the full-potential equation; and a numerical mapping procedure capable of generating curvilinear, body-fitted finite-difference grids about arbitrary airfoil geometries. Three aspects emphasized during the development of the TAIR code were reliability, simplicity, and speed. The reliability of TAIR comes from two sources: the new algorithm employed and the implementation of effective convergence monitoring logic. TAIR achieves ease of use by employing a "default mode" that greatly simplifies code operation, especially by inexperienced users, and many useful options including: several airfoil-geometry input options, flexible user controls over program output, and a multiple solution capability. The speed of the TAIR code is attributed to the new algorithm and the manner in which it has been implemented. Input to the TAIR program consists of airfoil coordinates, aerodynamic and flow-field convergence parameters, and geometric and grid convergence parameters. The airfoil coordinates for many airfoil shapes can be generated in TAIR from just a few input parameters. Most of the other input parameters have default values which allow the user to run an analysis in the default mode by specifing only a few input parameters. Output from TAIR may include aerodynamic coefficients, the airfoil surface solution, convergence histories, and printer plots of Mach number and density contour maps. The TAIR program is written in FORTRAN IV for batch execution and has been implemented on a CDC 7600 computer with a central memory requirement of approximately 155K (octal) of 60 bit words. The TAIR program was developed in 1981.

Description: The Vibration Pattern Imager (VPI) system was designed to control and acquire data from laser vibrometer sensors. The PC computer based system uses a digital signal processing (DSP) board and an analog I/O board to control the sensor and to process the data. The VPI system was originally developed for use with the Ometron VPI Sensor (Ometron Limited, Kelvin House, Worsley Bridge Road, London, SE26 5BX, England), but can be readily adapted to any commercially available sensor which provides an analog output signal and requires analog inputs for control of mirror positioning. VPI's graphical user interface allows the operation of the program to be controlled interactively through keyboard and mouse-selected menu options. The main menu controls all functions for setup, data acquisition, display, file operations, and exiting the program. Two types of data may be acquired with the VPI system: single point or "full field". In the single point mode, time series data is sampled by the A/D converter on the I/O board at a user-defined rate for the selected number of samples. The position of the measuring point, adjusted by mirrors in the sensor, is controlled via a mouse input. In the "full field" mode, the measurement point is moved over a user-selected rectangular area with up to 256 positions in both x and y directions. The time series data is sampled by the A/D converter on the I/O board and converted to a root-mean-square (rms) value by the DSP board. The rms "full field" velocity distribution is then uploaded for display and storage. VPI is written in C language and Texas Instruments' TMS320C30 assembly language for IBM PC series and compatible computers running MS-DOS. The program requires 640K of RAM for execution, and a hard disk with 10Mb or more of disk space is recommended. The program also requires a mouse, a VGA graphics display, a Four Channel analog I/O board (Spectrum Signal Processing, Inc.; Westborough, MA), a break-out box and a Spirit-30 board (Sonitech International, Inc.; Wellesley, MA) which includes a TMS320C30 DSP processor, 256Kb zero wait state SRAM, and a daughter board with 8Mb one wait state DRAM. Please contact COSMIC for additional information on required hardware and software. In order to compile the provided VPI source code, a Microsoft C version 6.0 compiler, a Texas Instruments' TMS320C30 assembly language compiler, and the Spirit 30 run time libraries are required. A math co-processor is highly recommended. A sample MS-DOS executable is provided on the distribution medium. The standard distribution medium for this program is one 5.25 inch 360K MS-DOS format diskette. The contents of the diskettes are compressed using the PKWARE archiving tools. The utility to unarchive the files, PKUNZIP.EXE, is included. VPI was developed in 1991-1992.

Description: The TAWFIVE program calculates transonic flow over a transport-type wing and fuselage. Although more complex Euler and Navier-Stokes methods are available, TAWFIVE combines a multi-grid acceleration technique in the iterative solution of the potential equation with the use of integral-form boundary-layer equations to provide a computationally efficient and sufficiently accurate design tool. TAWFIVE simplifies the solution process by breaking the problem into a loosely coupled set of modified equations. The inviscid method, using standard inviscid equations (nonlinear full potential), is valid in the "outer" region away from the wing, whereas the boundary-layer equations are valid in the thin region near the solid surface of the wing. The two types of equations are coupled by a technique of modifying surface boundary conditions for the inviscid equations. This interaction process starts with a solution of the outer flow field. Pressures are computed at the wing surface and are used to calculate the boundary layer. The boundary-layer and wake properties are then computed using a three-dimensional integral method, and the computed displacement thickness is added to the surface of the "hard" geometry. This new displaced wing surface is then regridded and the inviscid flowfield is recomputed. New values of the inviscid pressures are then used by the boundary-layer method to predict a new displacement thickness distribution. An under-relaxed update of the previously predicted displacement thickness is then made to obtain a new displacement thickness correction that is added to the "hard" geometry. These global iterations are continued until suitable convergence is obtained. Input to TAWFIVE is limited to geometric definition of the configuration, free-stream flow quantities, and iteration control parameters. The geometric input consists of the definition of a series of airfoil sections to define the wing and a series of fuselage cross sections to model the fuselage. High-aspect-ratio wings are modeled more accurately than low-aspect-ratio wings since no special provisions are made to accurately model the wing-fuselage juncture or the wingtip region. The user can specify the solution either in terms of lift or in terms of angle of attack. TAWFIVE can produce tabular output and input files for PLOT3D (COSMIC program number ARC-12779). TAWFIVE is written in FORTRAN 77 for CRAY series computers running UNICOS. The main memory requirement is 2.7Mb for execution. This program is available on a 9-track 1600 BPI UNIX tar format magnetic tape. TAWFIVE was under development from 1979 to 1989 and first released by COSMIC in 1991. CRAY and UNICOS are registered trademarks of Cray Research, Inc.

Description: Commercially available hot wires/films were used to measure the velocities of evaporated hydrogen or helium gas during cryogenic mixing experiments. Hot wires were found to be too delicate to use in this harsh environment. Hot films were rugged enough to use at cryogenic temperatures even though they failed after a number of thermal cycles. Since the hot films have small aspect ratios, 13.4 and 20, they are quite sensitive to the thermal loading, Tw/Tg, even with a correction for the conduction end loss. In general, although the increase of the Nusselt number with Reynolds number at low temperatures was similar to that at room temperature, there was also a pronounced variation with Tw/Tg over the large range of 1.2 to 12 investigated.

Description: A summary is presented of vortex control applications and current techniques for the control of longitudinal vortices produced by bodies, leading edges, tips and intersections. Vortex control has up till now been performed by many approaches in an empirical fashion, assisted by the essentially inviscid nature of much of longitudinal vortex behavior. Attention is given to Reynolds number sensitivities, vortex breakdown and interactions, vortex control on highly swept wings, and vortex control in juncture flows.

Description: The flow-field within an axial flow turbomachine, such as a turbine or compressor, is extremely complex because of three-dimensional features such as hub-corner stall, tip-leakage flows, and airfoil wakes. These flow features interact with each other and with rotor and stator airfoils inducing time-varying forces on the airfoils. These complicated rotor-stator interactions must be understood in order to design turbomachines that are light and compact as well as reliable and efficient. Two codes, STAGE-2 and STAGE-3, have been developed to compute these unsteady rotor-stator interaction flows in multistage turbomachines. An implicit, thin-layer Euler/Navier-Stokes zonal algorithm is used to compute the unsteady flow-field within both turbine and compressor configurations. Results include surface pressures and wake profiles for two-dimensional turbine and compressor configurations and surface pressures for a three-dimensional single-stage turbine configuration. The results compare well with experimental data and other unsteady computations.

Description: This paper demonstrates that optical microlithography can be used to produce a crossed grating which diffracts light into multiple orders sufficient to record moire interferograms with sensitivities ranging from 2.0 to 0.285 micron/fringe. The grating profile produced by the method is analyzed to establish the diffraction efficiency in each diffraction order, and generalized expressions are given for variable sensitivity moire interferometry. Experimental tests are conducted to verify analytical arguments. In one of these tests, two different diffraction order pairs are used simultaneously to verify that surface displacement can be measured at different sensitivities.

Description: Jet noise and jet-induced structural loads have become key issues in the design of commercial and military aircraft. Computational Fluid Dynamics (CFD) can be of use in predicting the underlying jet shear-layer instabilities and, in conjunction with classical acoustic theory, jet noise. The computational issues involved in the resolution of high Reynolds number unsteady jet flows are addressed in this paper. Once these jet flows can be accurately resolved, it should be possible to use acoustic theory to extract, for example, the far-field jet noise. An assessment of future work and computational resources required for directly computing far-field jet noise is also presented.

Description: We report the use of a short-length, multimode sapphire rod as an extension to a Michelson configuration, but operated as a low-finesse Fabry-Perot cavity. We demonstrate the performance of such a device as an interferometric sensor, where the interference between the reflections from the sapphire-air interface and an air-metallic surface is observed for microdisplacement of the metallic surface which is placed close to the sapphire endface. We describe in detail the fabrication procedure and present results obtained from the detection of temperature changes, applied strain, and surface acoustic waves.

Description: This computer program is designed to calculate the flow fields in two-dimensional and three-dimensional axisymmetric supersonic inlets. The method of characteristics is used to compute arrays of points in the flow field. At each point the total pressure, local Mach number, local flow angle, and static pressure are calculated. This program can be used to design and analyze supersonic inlets by determining the surface compression rates and throat flow properties. The program employs the method of characteristics for a perfect gas. The basic equation used in the program is the compatibility equation which relates the change in stream angle to the change in entropy and the change in velocity. In order to facilitate the computation, the flow field behind the bow shock wave is broken into regions bounded by shock waves. In each region successive rays are computed from a surface to a shock wave until the shock wave intersects a surface or falls outside the cowl lip. As soon as the intersection occurs a new region is started and the previous region continued only in the area in which it is needed, thus eliminating unnecessary calculations. The maximum number of regions possible in the program is ten, which allows for the simultaneous calculations of up to nine shock waves. Input to this program consists of surface contours, free-stream Mach number, and various calculation control parameters. Output consists of printed and/or plotted results. For plotted results an SC-4020 or similar plotting device is required. This program is written in FORTRAN IV to be executed in the batch mode and has been implemented on a CDC 7600 with a central memory requirement of approximately 27k (octal) of 60 bit words.

Description: This program was developed to predict turbine stage performance taking into account the effects of complex passage geometries. The method uses a quasi-3D inviscid-flow analysis iteratively coupled to calculated losses so that changes in losses result in changes in the flow distribution. In this manner the effects of both the geometry on the flow distribution and the flow distribution on losses are accounted for. The flow may be subsonic or shock-free transonic. The blade row may be fixed or rotating, and the blades may be twisted and leaned. This program has been applied to axial and radial turbines, and is helpful in the analysis of mixed flow machines. This program is a combination of the flow analysis programs MERIDL and TSONIC coupled to the boundary layer program BLAYER. The subsonic flow solution is obtained by a finite difference, stream function analysis. Transonic blade-to-blade solutions are obtained using information from the finite difference, stream function solution with a reduced flow factor. Upstream and downstream flow variables may vary from hub to shroud and provision is made to correct for loss of stagnation pressure. Boundary layer analyses are made to determine profile and end-wall friction losses. Empirical loss models are used to account for incidence, secondary flow, disc windage, and clearance losses. The total losses are then used to calculate stator, rotor, and stage efficiency. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370/3033 under TSS with a central memory requirement of approximately 4.5 Megs of 8 bit bytes. This program was developed in 1985.

Description: Turbomachinery components are often connected by ducts, which are usually annular. The configurations and aerodynamic characteristics of these ducts are crucial to the optimum performance of the turbomachinery blade rows. The ANDUCT computer program was developed to calculate the velocity distribution along an arbitrary line between the inner and outer walls of an annular duct with axisymmetric swirling flow. Although other programs are available for duct analysis, the use of the velocity gradient method makes the ANDUCT program fast and convenient while requiring only modest computer resources. A fast and easy method of analyzing the flow through a duct with axisymmetric flow is the velocity gradient method, also known as the stream filament or streamline curvature method. This method has been used extensively for blade passages but has not been widely used for ducts, except for the radial equilibrium equation. In ANDUCT, a velocity gradient equation derived from the momentum equation is used to determine the velocity variation along an arbitrary straight line between the inner and outer wall of an annular duct. The velocity gradient equation is used with an assumed variation of meridional streamline curvature. Upstream flow conditions may vary between the inner and outer walls, and an assumed total pressure distribution may be specified. ANDUCT works best for well-guided passages and where the curvature of the walls is small as compared to the width of the passage. The ANDUCT program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 60K of 8 bit bytes. The ANDUCT program was developed in 1982.

Description: The Panel Code for Planar Cascades was developed as an aid for the designer of turbomachinery blade rows. The effective design of turbomachinery blade rows relies on the use of computer codes to model the flow on blade-to-blade surfaces. Most of the currently used codes model the flow as inviscid, irrotational, and compressible with solutions being obtained by finite difference or finite element numerical techniques. While these codes can yield very accurate solutions, they usually require an experienced user to manipulate input data and control parameters. Also, they often limit a designer in the types of blade geometries, cascade configurations, and flow conditions that can be considered. The Panel Code for Planar Cascades accelerates the design process and gives the designer more freedom in developing blade shapes by offering a simple blade-to-blade flow code. Panel, or integral equation, solution techniques have been used for several years by external aerodynamicists who have developed and refined them into a primary design tool of the aircraft industry. The Panel Code for Planar Cascades adapts these same techniques to provide a versatile, stable, and efficient calculation scheme for internal flow. The code calculates the compressible, inviscid, irrotational flow through a planar cascade of arbitrary blade shapes. Since the panel solution technique is for incompressible flow, a compressibility correction is introduced to account for compressible flow effects. The analysis is limited to flow conditions in the subsonic and shock-free transonic range. Input to the code consists of inlet flow conditions, blade geometry data, and simple control parameters. Output includes flow parameters at selected control points. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 590K of 8 bit bytes. This program was developed in 1982.

Description: An exact, full-potential-equation model for the steady, irrotational, homoentropic, and homoenergetic flow of a compressible, inviscid fluid through a two-dimensional planar cascade together with its appropriate boundary conditions has been derived. The CAS2D computer program numerically solves an artificially time-dependent form of the actual full-potential-equation, providing a nonrotating blade-to-blade, steady, potential transonic cascade flow analysis code. Comparisons of results with test data and theoretical solutions indicate very good agreement. In CAS2D, the governing equation is discretized by using type-dependent, rotated finite differencing and the finite area technique. The flow field is discretized by providing a boundary-fitted, nonuniform computational mesh. This mesh is generated by using a sequence of conformal mapping, nonorthogonal coordinate stretching, and local, isoparametric, bilinear mapping functions. The discretized form of the full-potential equation is solved iteratively by using successive line over relaxation. Possible isentropic shocks are captured by the explicit addition of an artificial viscosity in a conservative form. In addition, a four-level, consecutive, mesh refinement feature makes CAS2D a reliable and fast algorithm for the analysis of transonic, two-dimensional cascade flows. The results from CAS2D are not directly applicable to three-dimensional, potential, rotating flows through a cascade of blades because CAS2D does not consider the effects of the Coriolis force that would be present in the three-dimensional case. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 200K of 8 bit bytes. The CAS2D program was developed in 1980.

Description: A computer program, QSONIC, has been developed for calculating the full potential, transonic quasi-three-dimensional flow through a rotating turbomachinery blade row. The need for lighter, more efficient turbomachinery components has led to the consideration of machines with fewer stages, each with blades capable of higher speeds and higher loading. As speeds increase, the numerical problems inherent in the transonic regime have to be resolved. These problems include the calculation of imbedded shock discontinuities and the dual nature of the governing equations, which are elliptic in the subcritical flow regions but become hyperbolic for supersonic zones. QSONIC provides the flow analyst with a fast and reliable means of obtaining the transonic potential flow distribution on a blade-to-blade stream surface of a stationary or rotating turbomachine blade row. QSONIC combines several promising transonic analysis techniques. The full potential equation in conservative form is discretized at each point on a body-fitted period mesh. A mass balance is calculated through the finite volume surrounding each point. Each local volume is corrected in the third dimension for any change in stream-tube thickness along the stream tube. The nonlinear equations for all volumes are of mixed type (elliptic or hyperbolic) depending on the local Mach number. The final result is a block-tridiagonal matrix formulation involving potential corrections at each grid point as the unknowns. The residual of each system of equations is solved along each grid line. At points where the Mach number exceeds unity, the density at the forward (sweeping) edge of the volume is replaced by an artificial density. This method calculates the flow field about a cascade of arbitrary two-dimensional airfoils. Three-dimensional flow is approximated in a turbomachinery blade row by correcting for stream-tube convergence and radius change in the through flow direction. Several significant assumptions were made in developing the QSONIC program, including: (1) the flow is inviscid and adiabatic, (2) the flow relative to the blade is steady, (3) the fluid is a perfect gas with constant specific heat, (4) the flow is isentropic and any discontinuities (shocks) are weak enough to be approximated as isentropic jumps, (5) there is no velocity component normal to the stream surface, and (6) the flow relative to a fixed frame in space (absolute velocity) is completely irrotational. These assumptions place some limitations on the application of QSONIC. Sharp leading edges at high incidence and high-Mach-number turbine blade trailing edges with substantial deviation will both cause large velocity peaks on the blade. In addition, the program may have difficulty converging if the passage is nearly choked. Input to QSONIC consists of case control parameters, a geometry description, upstream boundary conditions, and a rotor description. Output includes solution scheme parameters and flow field parameters. A data file is also output which contains data on the solution mesh, surface Mach numbers, surface static pressures, isomachs, and the velocity vector field. This data may be used for further processing or for plotting. The QSONIC is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 500K of 8 bit bytes. QSONIC was developed in 1982.

Description: This computer program, WIND, was developed to numerically solve the exact, full-potential equation for three-dimensional, steady, inviscid flow through an isolated wind turbine rotor. The program automatically generates a three-dimensional, boundary-conforming grid and iteratively solves the full-potential equation while fully accounting for both the rotating and Coriolis effects. WIND is capable of numerically analyzing the flow field about a given blade shape of the horizontal-axis type wind turbine. The rotor hub is assumed representable by a doubly infinite circular cylinder. An arbitrary number of blades may be attached to the hub and these blades may have arbitrary spanwise distributions of taper and of the twist, sweep, and dihedral angles. An arbitrary number of different airfoil section shapes may be used along the span as long as the spanwise variation of all the geometeric parameters is reasonably smooth. The numerical techniques employed in WIND involve rotated, type-dependent finite differencing, a finite volume method, artificial viscosity in conservative form, and a successive overrelaxation combined with the sequential grid refinement procedure to accelerate the iterative convergence rate. Consequently, WIND is cabable of accurately analyzing incompressible and compressible flows, including those that are locally transonic and terminated by weak shocks. Along with the three-dimensional results, WIND provides the results of the two-dimensional calculations to aid the user in locating areas of possible improvement in the aerodynamic design of the blade. Output from WIND includes the chordwise distribution of the coefficient of pressure, the Mach number, the density, and the relative velocity components at spanwise stations along the blade. In addition, the results specify local values of the lift coefficient and the tangent and axial aerodynamic force components. These are also given in integrated form expressing the total torque and the total axial force acting on the shaft. WIND can also be used to analyze the flow around isolated aircraft propellers and helicopter rotors in hover as long as the relative oncoming flow is subsonic. The WIND program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 253K of 8 bit bytes. WIND was developed in 1980.

Description: This computer program calculates the flow field in the supersonic portion of a mixed-compression aircraft inlet at non-zero angle of attack. This approach is based on the method of characteristics for steady three-dimensional flow. The results of this program agree with those produced by the two-dimensional method of characteristics when axisymmetric flow fields are calculated. Except in regions of high viscous interaction and boundary layer removal, the results agree well with experimental data obtained for threedimensional flow fields. The flow field in a variety of axisymmetric mixed compression inlets can be calculated using this program. The bow shock wave and the internal shock wave system are calculated using a discrete shock wave fitting procedure. The internal flow field can be calculated either with or without the discrete fitting of the internal shock wave system. The influence of molecular transport can be included in the calculation of the external flow about the forebody and in the calculation of the internal flow when internal shock waves are not discretely fitted. The viscous and thermal diffussion effects are included by treating them as correction terms in the method of characteristics procedure. Dynamic viscosity is represented by Sutherland's law and thermal conductivity is represented as a quadratic function of temperature. The thermodynamic model used is that of a thermally and calorically perfect gas. The program assumes that the cowl lip is contained in a constant plane and that the centerbody contour and cowl contour are smooth and have continuous first partial derivatives. This program cannot calculate subsonic flow, the external flow field if the bow shock wave does not exist entirely around the forebody, or the internal flow field if the bow flow field is injected into the annulus. Input to the program consists of parameters to control execution, to define the geometry, and the vehicle orientation. Output consists of a list of parameters used, solution planes, and a description of the shock waves. This program is written in FORTRAN IV for batch execution and has been implemented on a CDC 6000 series machine with a central memory requirement of 110K (octal) of 60 bit words when it is overlayed. This flow analysis program was developed in 1978.

Description: This computer program was developed for calculating the subsonic or transonic flow on the hub-shroud mid-channel stream surface of a single blade row of a turbomachine. The design and analysis of blades for compressors and turbines ideally requires methods for analyzing unsteady, three-dimensional, turbulent viscous flow through a turbomachine. Since an exact solution is impossible at present, solutions on two-dimensional surfaces are calculated to obtain a quasi-three dimensional solution. When three-dimensional effects are important, significant information can be obtained from a solution on a cross-sectional surface of the passage normal to the flow. With this program, a solution to the equations of flow on the meridional surface can be carried out. This solution is chosen when the turbomachine under consideration has significant variation in flow properties in the hubshroud direction, especially when input is needed for use in blade-to-blade calculations. The program can also perform flow calculations for annular ducts without blades. This program should prove very useful in the design and analysis of any turbomachine. This program calculates a solution for two-dimensional, adiabatic shockfree flow. The flow must be essentially subsonic, but there may be local areas of supersonic flow. To obtain the solution, this program uses both the finite difference and the quasi-orthogonal (velocity gradient) methods combined in a way that takes maximum advantage of both. The finite-difference method solves a finite-difference equation along the meridional stream surface in a very efficient manner but is limited to subsonic velocities. This approach must be used in cases where the blade aspect ratios are above one, cases where the passage is curved, and cases with low hub-tip-ratio blades. The quasi-orthogonal method solves the velocity gradient equation on the meridional surface and is used if it is necessary to extend the range of solutions into the transonic regime. In general the blade row may be fixed or rotating and the blades may be twisted and leaned. The flow may be axial, radial, or mixed. The upstream and downstream flow conditions can vary from hub to shroud with provisions made for an approximate correction for loss of stagnation pressure. Also, viscous forces are neglected along solution mesh lines running from hub to tip. The capabilities of this program include handling of nonaxial flows without restriction, annular ducts without blades, and specified streamwise loss distributions. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 700K of 8 bit bytes. This core requirement can be reduced depending on the size of the problem and the desired solution accuracy. This program was developed in 1977.

Description: A computer program has been developed for the design of supersonic rotor blades where losses are accounted for by correcting the ideal blade geometry for boundary layer displacement thickness. The ideal blade passage is designed by the method of characteristics and is based on establishing vortex flow within the passage. Boundary-layer parameters (displacement and momentum thicknesses) are calculated for the ideal passage, and the final blade geometry is obtained by adding the displacement thicknesses to the ideal nozzle coordinates. The boundary-layer parameters are also used to calculate the aftermixing conditions downstream of the rotor blades assuming the flow mixes to a uniform state. The computer program input consists essentially of the rotor inlet and outlet Mach numbers, upper- and lower-surface Mach numbers, inlet flow angle, specific heat ratio, and total flow conditions. The program gas properties are set up for air. Additional gases require changes to be made to the program. The computer output consists of the corrected rotor blade coordinates, the principal boundary-layer parameters, and the aftermixing conditions. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 7094. This program was developed in 1971.

Description: This program obtains a transonic flow solution on a blade-to-blade surface between blades of a turbomachine. The flow must be essentially subsonic, but there may be locally supersonic flow. The solution is two-dimensional, isentropic, and shock free. The blades may be fixed or rotating. The flow may be axial, radial, or mixed, and there may be a change in stream-channel thickness in the through-flow direction. A loss in relative stagnation pressure may be accounted for. The program input consists of blade and stream-channel geometry, stagnation flow conditions, inlet and outlet flow angles, and blade-to-blade stream-channel weight flow. The output includes blade surface velocities, velocity magnitude and direction at all interior mesh points in the blade-to-blade passage, and streamline coordinates throughout the passage. The transonic solution is obtained by a combination of a finite-difference, stream-function solution and a velocity-gradient solution. The finite-difference solution at a reduced weight flow provides information needed to obtain a velocity-gradient solution. This program is written in FORTRAN IV for batch execution and has been implemented on the IBM 360 computer with a central memory requirement of approximately 36K of 8 bit bytes. This program was developed in 1969 and last updated in 1979.

Description: This program is a revision of an existing program for blade-to-blade aerodynamic analysis of turbomachine blades and it is a simpler program while consistent with related programs. The analysis is for two-dimensional, subsonic, compressible (or incompressible), nonviscous flow in a circular or straight infinite cascade of blades, which may be fixed or rotating. The flow may be axial, radial, or mixed, and the stream channel thickness may change in the through-flow direction. The program input consists of blade and stream channel geometry, total flow conditions, inlet and outlet flow angles, and blade-to-blade stream channel weight flow. The output includes blade surface velocities, velocity magnitude and direction at all interior mesh points in the blade-to-blade passage, and streamline coordinates throughout the passage. This program was developed on an IBM 7094/7044 DCS.

Description: This computer program gives the blade-to-blade solution of the two-dimensional, subsonic, compressible (or incompressible), nonviscous flow problem for a circular or straight infinite cascade of tandem or slotted turbomachine blades. The blades may be fixed or rotating. The flow may be axial, radial , or mixed. The method of solution is based on the stream function using an iterative solution of nonlinear finite-difference equations. These equations are solved using two major levels of iteration. The inner iteration consists of the solution of simultaneous linear equations by successive over-relaxation, using an estimated optimum over-relaxation factor. The outer iteration then changes the coefficients of the simultaneous equations to correct for compressibility. The program input consists of the basic blade geometry, the meridional stream channel coordinates, fluid stagnation conditions, weight flow and flow split through the slot, and inlet and outlet flow angles. The output includes blade surface velocities, velocity magnitude and direction throughout the passage, and the streamline coordinates.

Description: This FORTRAN IV computer program which incorporates the method of characteristics was written to assist in the design of supersonic inlets. There were two objectives: (1) to study a greater variety of supersonic inlet configurations and (2) to reduce the time required for trial-and-error procedures to arrive at optimum inlet design. The computer program was written with the intention of being able to construct a variety of inlet configurations by interchanging specific subroutines. In this manner, greater flexibility of choice was attained, and the time required to program a specific inlet configuration was greatly reduced. The second objective was accomplished by a reformulation of the boundary value problem for hyperbolic equations. By this reformulation of the boundary data, the engineering design quantities, throat Mach number and flow angle, were introduced as direct input quantities to the computer program. As a consequence of introducing the engineering parameters as input, the computer program will calculate the surface contours required to satisfy the specific throat conditions. Inviscid flow is assumed and the method used to calculate the inlet contour results in minimum distortion to the flow in the throat. This program was developed on an IBM 7094.

Description: This program represents a subsonic aerodynamic method for determining the mean camber surface of trimmed noncoplaner planforms with minimum vortex drag. With this program, multiple surfaces can be designed together to yield a trimmed configuration with minimum induced drag at some specified lift coefficient. The method uses a vortex-lattice and overcomes previous difficulties with chord loading specification. A Trefftz plane analysis is used to determine the optimum span loading for minimum drag. The program then solves for the mean camber surface of the wing associated with this loading. Pitching-moment or root-bending-moment constraints can be employed at the design lift coefficient. Sensitivity studies of vortex-lattice arrangements have been made with this program and comparisons with other theories show generally good agreement. The program is very versatile and has been applied to isolated wings, wing-canard configurations, a tandem wing, and a wing-winglet configuration. The design problem solved with this code is essentially an optimization one. A subsonic vortex-lattice is used to determine the span load distribution(s) on bent lifting line(s) in the Trefftz plane. A Lagrange multiplier technique determines the required loading which is used to calculate the mean camber slopes, which are then integrated to yield the local elevation surface. The problem of determining the necessary circulation matrix is simplified by having the chordwise shape of the bound circulation remain unchanged across each span, though the chordwise shape may vary from one planform to another. The circulation matrix is obtained by calculating the spanwise scaling of the chordwise shapes. A chordwise summation of the lift and pitching-moment is utilized in the Trefftz plane solution on the assumption that the trailing wake does not roll up and that the general configuration has specifiable chord loading shapes. VLMD is written in FORTRAN for IBM PC series and compatible computers running MS-DOS. This program requires 360K of RAM for execution. The Ryan McFarland FORTRAN compiler and PLINK86 are required to recompile the source code; however, a sample executable is provided on the diskette. The standard distribution medium for VLMD is a 5.25 inch 360K MS-DOS format diskette. VLMD was originally developed for use on CDC 6000 series computers in 1976. It was originally ported to the IBM PC in 1986, and, after minor modifications, the IBM PC port was released in 1993.

Description: This code was developed to aid design engineers in the selection and evaluation of aerodynamically efficient wing-canard and wing-horizontal-tail configurations that may employ simple hinged-flap systems. Rapid estimates of the longitudinal aerodynamic characteristics of conceptual airplane lifting surface arrangements are provided. The method is particularly well suited to configurations which, because of high speed flight requirements, must employ thin wings with highly swept leading edges. The code is applicable to wings with either sharp or rounded leading edges. The code provides theoretical pressure distributions over the wing, the canard or horizontal tail, and the deflected flap surfaces as well as estimates of the wing lift, drag, and pitching moments which account for attainable leading edge thrust and leading edge separation vortex forces. The wing planform information is specified by a series of leading edge and trailing edge breakpoints for a right hand wing panel. Up to 21 pairs of coordinates may be used to describe both the leading edge and the trailing edge. The code has been written to accommodate 2000 right hand panel elements, but can easily be modified to accommodate a larger or smaller number of elements depending on the capacity of the target computer platform. The code provides solutions for wing surfaces composed of all possible combinations of leading edge and trailing edge flap settings provided by the original deflection multipliers and by the flap deflection multipliers. Up to 25 pairs of leading edge and trailing edge flap deflection schedules may thus be treated simultaneously. The code also provides for an improved accounting of hinge-line singularities in determination of wing forces and moments. To determine lifting surface perturbation velocity distributions, the code provides for a maximum of 70 iterations. The program is constructed so that successive runs may be made with a given code entry. To make additional runs, it is necessary only to add an identification record and the namelist data that are to be changed from the previous run. This code was originally developed in 1989 in FORTRAN V on a CDC 6000 computer system, and was later ported to an MS-DOS environment. Both versions are available from COSMIC. There are only a few differences between the PC version (LAR-14458) and CDC version (LAR-14178) of AERO2S distributed by COSMIC. The CDC version has one main source code file while the PC version has two files which are easier to edit and compile on a PC. The PC version does not require a FORTRAN compiler which supports NAMELIST because a special INPUT subroutine has been added. The CDC version includes two MODIFY decks which can be used to improve the code and prevent the possibility of some infrequently occurring errors while PC-version users will have to make these code changes manually. The PC version includes an executable which was generated with the Ryan McFarland/FORTRAN compiler and requires 253K RAM and an 80x87 math co-processor. Using this executable, the sample case requires about four hours to execute on an 8MHz AT-class microcomputer with a co-processor. The source code conforms to the FORTRAN 77 standard except that it uses variables longer than six characters. With two minor modifications, the PC version should be portable to any computer with a FORTRAN compiler and sufficient memory. The CDC version of AERO2S is available in CDC NOS Internal format on a 9-track 1600 BPI magnetic tape. The PC version is available on a set of two 5.25 inch 360K MS-DOS format diskettes. IBM AT is a registered trademark of International Business Machines. MS-DOS is a registered trademark of Microsoft Corporation. CDC is a registered trademark of Control Data Corporation. NOS is a trademark of Control Data Corporation.

Description: This program provides a wing design algorithm based on modified linear theory which takes into account the effects of attainable leading-edge thrust. A primary objective of the WINGDES2 approach is the generation of a camber surface as mild as possible to produce drag levels comparable to those attainable with full theoretical leading-edge thrust. WINGDES2 provides both an analysis and a design capability and is applicable to both subsonic and supersonic flow. The optimization can be carried out for designated wing portions such as leading and trailing edge areas for the design of mission-adaptive surfaces, or for an entire planform such as a supersonic transport wing. This program replaces an earlier wing design code, LAR-13315, designated WINGDES. WINGDES2 incorporates modifications to improve numerical accuracy and provides additional capabilities. A means of accounting for the presence of interference pressure fields from airplane components other than the wing and a direct process for selection of flap surfaces to approach the performance levels of the optimized wing surfaces are included. An increased storage capacity allows better numerical representation of those configurations that have small chord leading-edge or trailing-edge design areas. WINGDES2 determines an optimum combination of a series of candidate surfaces rather than the more commonly used candidate loadings. The objective of the design is the recovery of unrealized theoretical leading-edge thrust of the input flat surface by shaping of the design surface to create a distributed thrust and thus minimize drag. The input consists of airfoil section thickness data, leading and trailing edge planform geometry, and operational parameters such as Mach number, Reynolds number, and design lift coefficient. Output includes optimized camber surface ordinates, pressure coefficient distributions, and theoretical aerodynamic characteristics. WINGDES2 is written in FORTRAN V for batch execution and has been implemented on a CDC CYBER computer operating under NOS 2.7.1 with a central memory requirement of approximately 344K (octal) of 60 bit words. This program was developed in 1984, and last updated in 1990. CDC and CYBER are trademarks of Control Data Corporation.

Description: Progress in developing fiber-optic interferometric sensors for aeroacoustic measurements in wind tunnels, performed under the NASA program, is reported. Preliminary results show that the fiber-optic interferometer sensor array is a powerful instrument for solving complex acoustic measurement problems in wind tunnels, which cannot be resolved with the conventional transducer technique.

Description: Multiresponse imaging is a process that acquires A images, each with a different optical response, and reassembles them into a single image with an improved resolution that can approach 1/sq rt A times the photodetector-array sampling lattice. Our goals are to optimize the performance of this process in terms of the resolution and fidelity of the restored image and to assess the amount of information required to do so. The theoretical approach is based on the extension of both image restoration and rate-distortion theories from their traditional realm of signal processing to image processing which includes image gathering and display.

Description: Unsteady flow behavior and load characteristics of a 2D VR-7 airfoil with and without a leading-edge slat were studied in the water tunnel of the Aeroflightdynamics Directorate, NASA Ames Research Center. Both airfoils were oscillated sinusoidally between 5 and 25 deg at Re = 200,000 to obtain the unsteady lift, drag, and pitching moment data. A fluorescent dye was released from an orifice located at the leading edge of the airfoil for the purpose of visualizing the boundary layer and wake flow. The flowfield and load predictions of an incompressible Navier-Stokes code based on a velocity-vorticity formulation were compared with the test data. The test and predictions both confirm that the slatted VR-7 airfoil delays both static and dynamic stall as compared to the VR-7 airfoil alone.

Description: The Far Infrared Limb Observing Spectrometer (FILOS) is an instrument designed to measure chemical species in the upper atmosphere using limb emission in the FIR region of the spectrum. FILOS uses three Fabry-Perot etalons in series to obtain a resolution of 0.0017/cm near 100/cm (100 microns). It is compact and has low power and low data rate requirements so that it may be flown as an auxiliary balloon payload with larger instruments. FILOS has two 0.05/cm bandwidth channels which are currently tuned to a HCl line at 104.2/cm and a pair of OH lines at 101.3/cm. The instrument is described in further detail and results are presented from two recent balloon flights in which OH was measured as a function of time on one hour centers from sunrise to sunset.

Description: Experimental results are reported for submerged injection pressurization and expulsion tests of a 4.89 cu m liquid hydrogen tank. The pressurant injector was positioned near the bottom of the test vessel to simulate liquid engulfment of the pressurant gas inlet, a condition that may occur in low-gravity conditions. Results indicate a substantial reduction in pressurization efficiency with pressurant gas requirements approximately five times greater than ideal amounts. Consequently, submerged vapor injection should be avoided as a low-gravity autogenous pressurization method whenever possible. The work presented herein validates that pressurant requirements are accurately predicted by a homogeneous thermodynamic model when the submerged injection technique is employed.

Description: In an effort to improve the signal to noise in an interference experiment, we have developed a method to remove systematic phase drift between data sets acquired over long time intervals. Using this technique, it is possible to average repeatedly acquired phase measurements and improve the phase estimate without sacrificing spatial resolution. Results from tests using real-time phase stepping holographic interferometry applied to cantilever bending of a piezoelectric bimorph indicate that white noise has been reduced from 3 to less than 1 deg (lambda/360) by averaging 36 phase compensated data sets before object bending and 36 data sets after bending.

Description: The potential for using high-temperature superconductive elements, screen-printed onto ceramic substrates, as thermal bridges to replace the currently employed manganin wires is studied at NASA-LaRC. Substrate selection is considered to be the most critical parameter in device production. Due to the glass-like thermal behavior of yttria-stabilized-zirconia (YSZ) and fused silica substrates, these materials are found to reduce the heat load significantly. The estimated thermal savings for superconductive leads printed onto YSZ or fused silica substrates range from 6 to 14 percent.

Description: The problem of the hypersonic double ellipse in rarefied flow is treated by a particle method using the collision model first described by McDonald (1988). In the approach used here, the computational overhead is reduced by using simple cubic cells. The problem of the definition of complex geometries is addressed by developing an algorithm to define the relation of a body surface to the network of cells.

Description: Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D viscous hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The code has been applied to Problem 8.2 of this workshop for the case of thermochemical nonequilibrium flow through a nozzle. Chemical reaction rates are defined with the model of Park (1987). Thermal nonequilibrium is modeled using a two-temperature approximation in which the vibrational energies of all molecules are assumed to be in equilibrium at a single temperature which is generally different from the translational-rotational temperature. Two grids were used to define the flow for the original problem, with a stagnation temperature of 6500 K. A third case with a stagnation temperature of 10,000 K is also presented. The solution domain includes the converging nozzle, subsonic flow domain in which the gas is substantially in thermochemical equilibrium and the diverging nozzle, hypersonic flow domain in which the gas is substantially in thermochemical nonequilibrium.

Description: Solutions have been computed and results are presented for Problem 1, the case of Mach 9 transitional flow past a 7 deg half-angle cone at zero incidence. The solutions were computed using a code developed for the integration of the parabolized Navier-Stokes equations. The algorithm employed in the code is based on a Roe-type flux-difference-splitting scheme applied following a finite-volume approach. The basic algorithm has been modified to make it implicit and second-order accurate in the crossflow directions. Results are presented in terms of surface pressure and heat transfer as well as boundary layer profiles of pitot pressure, Mach number, and tangential velocity. The case was recalculated several times in an effort to determine sensitivities to such parameters as grid density, wall temperature, turbulence model parameters, as well as freestream expansion. Comparisons with the experimental data are presented and discussed.

Description: A development status evaluation is presented for the theoretical understanding and design conceptualization of boundary layer control (BLC) systems applicable to supersonic transports, such as the currently envisioned NASA High Speed Civil Transport. By reducing fuel burned, supersonic BLC techniques could expand ranges to Pacific-crossing scales, while lowering sonic boom effects and upper-atmosphere pollution and even reducing skin friction temperature. The critical consideration for supersonic BLC is the presence of wave effects.

Description: Attention is drawn to the influence of preexisting finite-amplitude instabilities on the growth of other disturbances; current design tools for LFC take no notice of this kind of interaction. When a rational accounting is accomplished for the evolution of incoming disturbances in finite-amplitude solutions of the equations of motion, future transition-prediction methods will need to take these wave interactions into account. Attention is given here to interactions in the presence of crossflow vortices and interactions involving Goertler vortices.

Description: A development history and a development-trends evaluation are presented for laminar flow controlled airfoil technologies and design concepts, including the search for 'natural' laminar flow and actively controlled flow via suction through small pores on the airfoil surface. While most NASA activities in this field have been concerned with subsonic aircraft, it has been projected that the control of boundary layer turbulence may be even more critical to the aerodynamic efficiency of supersonic aircraft. Developmental programs for these techniques have been conducted with several modified conventional aircraft.

Description: An account is given of the development history of natural laminar-flow (NLF) airfoil profiles under guidance of an experimentally well-verified theoretical method for the design of airfoils suited to virtually all subcritical applications. This method, the Eppler Airfoil Design and Analysis Program, contains a conformal-mapping method for airfoils having prescribed velocity-distribution characteristics, as well as a panel method for the analysis of potential flow about given airfoils and a boundary-layer method. Several of the NLF airfoils thus obtained are discussed.

Description: The instantaneous velocity fields of time-dependent flows, or of a collection of objects moving with spatially varying velocities, can be measured by means of digital image velocimetry (DIV). DIV overcomes several shortcomings of such existing techniques as laser-speckle or particle-image velocimetry. Attention is presently given to numerically generated images representing objects in uniform motion which are then used for the experimental validation of DIV.

Description: Selective deposition and abrasion, as well as etching in atomic oxygen or reduced-pressure air, have been used to prepare patterned polycrystalline diamond films which, on further processing by anisotropic Si etching, yield the microstructures of such devices as flow sensors and accelerometers. Both types of sensor have been experimentally tested in the respective functions of hot-wire anemometer and both single- and double-hinged accelerometer.

Description: AVIRIS is a facility consisting of a flight system, a ground data system, a calibration facility, and a full-time operations team. The facility was developed by JPL under funding from NASA. NASA also provides funding for operations and maintenance. The flight system is a whisk-broom imager that acquires data in 224 narrow, contiguous spectral bands covering the solar reflected portion of the electromagnetic spectrum. It is flown aboard the NASA high altitude ER-2 research aircraft. The ground data system is a facility dedicated to the processing and distribution of data acquired by AVIRIS. It operates year round at JPL. The calibration facility consists of a calibration laboratory at JPL and a suite of field instruments and procedures for performing inflight calibration of AVIRIS. A small team of engineers, technicians, and scientists supports a yearly operations schedule that includes 6 months of flight operations, 6 months of routine ground maintenance of the flight system, and year-round data processing and distribution. Details of the AVIRIS system, its performance history, and future plans are described.

Description: The solar corona, supernova remnants, the hot diffuse interstellar gas in the Galaxy, galactic halos, and the hot intracluster gas in rich clusters of galaxies, are examples of extended astrophysical plasmas which emit line-rich spectra in the X-ray spectral range from 1.5 to 25 A. These phenomena represent a significant fraction of the baryonic matter in the universe. The study of the composition, structure and dynamics of these astrophysical plasmas requires observations with both high spectral and spatial resolution simultaneously. The Objective Double Crystal Spectrometer, coupled with a grazing incidence X-ray telescope, represents a stigmatic instrument which is highly efficient for the study of such sources. We describe the configuration and performance (spatial resolution, spectral resolution and efficiency) of the Objective Double Crystal spectrometer.

Description: The present treatment of vector magnetic field measurement in coronas by means of the Hanle effect of the Lyman-alpha line uses data from all-reflecting imaging coronagraph/polarimeters. The polarization sensitivity, bandpass, and spatial resolution of these instruments are defined through a modeling of the Hanle-effect signature in Lyman-alpha emission from coronal magnetic loops; the line-of-sight integration through an inhomogeneous coronal medium is taken into account. The use of the Hanle effect to measure solar corona vector magnetic fields is verified.

Description: Multilayer optics operated at normal incidence offer a powerful new technology for the study of the solar spectrum in the XUV. The spectra of most cosmic X-ray sources are strongly extinguished at wavelengths above 40 A due to absorption and scattering by interstellar grains. We describe a number of configurations which allow multilayer optics to be used at nonnormal angles of incidence in conjunction with grazing incidence optics to analyze the spectra of cosmic X-ray sources in the wavelength interval between 1.5 and 40 A. These optical configurations utilize both multilayer mirrors and gratings, and permit the efficient observation of extended sources using stigmatic spectrographs. The response of the instruments described to typical cosmic X-ray sources is also discussed.

Description: The rocketborne Multi-Spectral Solar Telescope Array (MSSTA) uses an array of Ritchey-Chretien, Cassegrain, and Herschelian telescopes to produce ultrahigh-resolution full-disk images of the sun within the soft X-ray, EUV, and FUV ranges. Such imaging of the solar disk and corona out to several solar radii placed great demands on the MSSTA's data storage capabilities; in addition, its photographic films required very low outgassing rates. Results are presented from calibration tests conducted on the MSSTA's emulsions, based on measurements at NIST's synchrotron facility.

Description: The spherical Schwarzschild microscope for soft X-ray applications in microscopy and projection lithography consists of two concentric spherical mirrors configured such that the third-order spherical aberration and coma are zero. Since multilayers are used on the mirror substrates for X-ray applications, it is desirable to have only two reflecting surfaces in a microscope. To reduce microscope aberrations and increase the field of view, generalized mirror surface profiles are here considered. Based on incoherent and sine wave modulation transfer function calculations, the object plane resolution of a microscope has been analyzed as a function of the object height and numerical aperture (NA) of the primary for several spherical Schwarzschild, conic, and aspherical Head reflecting two-mirror microscope configurations. The Head microscope with a NA of 0.4 achieves diffraction limited performance for objects with a diameter of 40 microns. Thus, it seems possible to record images with a feature size less than 100 A with a 40x microscope when using 40 A radiation.

Description: Consideration is given to a cryogenic multichannel electronically scanned pressure (ESP) module developed and tested over an extended temperature span from -184 to +50 C and a pressure range of 0 to 5 psig. The ESP module consists of 32 pressure sensor dice, four analog 8 differential-input multiplexers, and an amplifier circuit, all of which are packaged in a physical volume of 2 x 1 x 5/8 in with 32 pressure and two reference ports. Maximum nonrepeatability is measured at 0.21 percent of full-scale output. The ESP modules have performed consistently well over 15 times over the above temperature range and continue to work without any sign of degradation. These sensors are also immune to repeated thermal shock tests over a temperature change of 220 C/sec.

Description: The paper reports the development and initial testing of a digital resolver to replace existing analog signal processing instrumentation. Radiometers, mounted directly on one of the fully articulated blades, are electrically connected through a slip ring to analog signal processing circuitry. The measured signals are periodic with azimuth angle and are resolved into harmonic components, with 0 deg over the tail. The periodic nature of the helicopter blade motion restricts the frequency content of each flapping and yaw signal to the fundamental and harmonics of the rotor rotational frequency. A minicomputer is employed to collect these data and then plot them graphically in real time. With this and other information generated by the instrumentation, a helicopter test pilot can then adjust the helicopter model's controls to achieve the desired aerodynamic test conditions.

Description: Multi-component strain-gage force transducer design requires the designer to determine the spring constant of the numerous beams or flexures incorporated in the transducer. The classical beam deflection formulae that are used in calculating these spring constants typically assume that the beam has a uniform moment of inertia along the entire beam length. In practice all beams have a radius at the end where the beam interfaces with the shoulder of the transducer, and on short beams in particular this increases the beam spring constant considerably. A Basic computer program utilizing numerical integration is presented to determine this effect.

Description: The objective is to describe continuing efforts to develop methods for measuring surface heat flux, gauge active surface temperature, and heat transfer coefficient quantities. The methodology involves inventing a procedure for fabricating improved plug-type heat flux gauges and also for formulating inverse heat conduction models and calculation procedures. These models and procedures are required for making indirect measurements of these quantities from direct temperature measurements at gauge interior locations. Measurements of these quantities were made in a turbine blade thermal cycling tester (TBT) located at MSFC. The TBT partially simulates the turbopump turbine environment in the Space Shuttle Main Engine. After the TBT test, experiments were performed in an arc lamp to analyze gauge quality.

Description: Presented here are results of a test program undertaken to further define the response of the solar dynamic radiator to hypervelocity impact (HVI). Tests were conducted on representative radiator panels (under ambient, nonoperating conditions) over a range of velocity. Target parameters are also varied. Data indicate that analytical penetration predictions are conservative (i.e., pessimistic) for the specific configuration of the solar dynamic radiator. Test results are used to define the solar dynamic radiator reliability with respect to HVI more rigorously than previous studies. Test data, reliability, and survivability results are presented.

Description: A ground-based optical telescope system has been constructed with the capability to locate fast optical transients that may be associated with gamma-ray bursts (GRBs). The instrument has been integrated and operated during a shakedown period at GSFC, Maryland. Results of 35 hours of 'state mode' data are presented. The telescope has the proven capability to slew to any point on the night sky within 1.0 sec, track that position with better than one arcsecond stability, and image a 9 x 12 arcmin field of view with 1 arcsec angular resolution with 1.5 sec time resolution. The telescope-CCD camera system has a sensitivity of 13th magnitude for transients and 14th mag for field stars. In the 35 hr of operation many single frame transients of instrumental and optical origin have been observed; no two-sequential frame astrophysical transients have been identified. The combined rate of instrumental transients (predominantly sea-level muons) is 7.2/hr and of optical transients (satellite glints, airplane strobe lights, and meteors) is 5.1/hr. The RMT will operate in conjunction with the MIT Explosive Transient Camera survey instrument at Kitt Peak National Observatory, Tucson. The RMT is now being installed at Kitt Peak. Full operation will begin this summer.

Description: Attention is given to a space-borne engine plume experiment study to fly an experiment which will both verify and quantify the reduced contamination from advanced rhenium-iridium earth-storable bipropellant rockets (hot rockets) and provide a correlation between high-fidelity, in-space measurements and theoretical plume and surface contamination models. The experiment conceptual design is based on survey results from plume and contamination technologists throughout the U.S. With respect to shuttle use, cursory investigations validate Hitchhiker availability and adaptability, adequate remote manipulator system (RMS) articulation and dynamic capability, acceptable RMS attachment capability, adequate power and telemetry capability, and adequate flight altitude and attitude/orbital capability.

Description: A new procedure, dubbed the Munich Method, has been proposed recently for the modeling of rocket engine performance. The author of the Munich Method claims it to be an extension and improvement of the thermodynamic procedures used to model rocket engines in the NASA-Lewis chemical equilibrium program. An examination of the Munich Method shows that it contains several flaws. If these defects are corrected then the Munich Method will produce results identical to those generated by the NASA-Lewis Code.

Description: This report describes the calibration of a nonnulling, conical, seven-hole pressure probe over a large range of flow onset angles. The calibration procedure is based on the use of differential pressures to determine the three components of velocity. The method allows determination of the flow angle and velocity magnitude to within an average error of 1.0 deg and 1.0 percent, respectively. Greater accuracy can be achieved by using high-quality pressure transducers. Also included is an examination of the factors which limit the use of the probe, a description of the measurement chain, an error analysis, and a typical experimental result. In addition, a new general analytical model of pressure probe behavior is described, and the validity of the model is demonstrated by comparing it with experimentally measured calibration data for a three-hole yaw meter and a seven-hole probe.

Description: This paper describes a recently completed electrooptical camera flying onboard the NASA ER-2 high altitude aircraft. The device includes a six-position filter wheel which can be fitted with a combination of polarizing and/or spectral filters. An alternate configuration will include a polarizing filter which can be rotated to any angle under computer control. The camera mount in the nose of the ER-2 can tilt forward or aft up to 40 degrees, both for bidirectional reflectance studies and for image motion compensation (the aircraft moves 34 meters between frame acquisitions). The ground resolution is nominally 5 meters from and altitude of 20 km. Spectral responsivity is that of the silicon imaging array (Kodak KAF-1400). Initial data sets were acquired in support of the International Satellite Cloud Climatology Program Regional Experiment of November, 1991, and will be used to study cirrus cloud properties.

Description: Airloads measured on a two-bladed helicopter rotor in flight, from the Tip Aerodynamic and Acoustic Test, are compared with calculations from a comprehensive helicopter analysis (CAMRAD/JA), and the pressures compared with calculations from a full-potential rotor code (FPR). The flight test results cover an advance ratio range from 0.19 to 0.38. The lowest speed case is characterized by the presence of significant blade-vortex interactions. Good correlation of peak-to-peak vortex-induced loads and the corresponding pressures is obtained. The results of the correlation for this two-bladed rotor are substantially similar to the results for three- and four-bladed rotors, concerning the tip vortex core size for best correlation, calculation of the peak-to-peak loads on the retreating side, and calculation of vortex-induced loads on inboard radial stations.

Description: A new CFD potential code, FPX (eXtended Full-Potential), has been developed for application to both helicopters and tilt-rotors. The code solves the unsteady, three-dimensional full potential equation and is an extension of the rotor code, FPR. Both entropy and viscosity corrections are included to enhance the physical modeling capabilities. A number of efficiency related modifications have yielded a factor of two speed-up in the code. An axial flow capability has been added to treat tilt-rotor in forward flight (cruise mode). In order to employ streamwise periodicity and accurately solve for the propagation of acoustic signals in the tip region, an H-H topology has been added to the basic O-H grid system. Computations are performed for the XV-15 Standard and ATB blades at high-speed conditions. Comparisons are made for the blade aerodynamics and the induced fuselage cabin pressure for a range of Mach numbers. Grid generation, wake treatment, and far-field wall treatment are identified as problem areas with recommendations for future research.