ALBERT

Description: Helmet streamers on the sun were observed to be the site of coronal mass ejections, dynamic events that eject coronal plasma and magnetic fields into the solar wind. A two dimensional (azimuthally symmetric) helmet streamer configuration was developed by computing solutions of the time dependent magnetohydrodynamic (MHD) equations, for a specified magnetic flux distribution on the sun. The helmet streamer is not symmetric about the equator. The evolution of the configuration, when differential rotation is applied, was investigated. It was found that after many rotations the configuration does not reach a steady state, but disrupts recurrently with the ejection of a plasmoid. These results suggest that differential rotation may be one of the mechanisms by which mass ejections are initiated.

Description: The ground-based observing facilities of the National Solar Observatory (NSO) are reviewed from the perspective of joint observations with SOHO. A specific proposal is presented for observations of the HE-I 1083.0 nm line with the NASA/NSO spectromagnetograph and He 10830 video filtergraph/magnetograph in coordination with ultraviolet sensitive instruments on SOHO. The first task will be to look for associations of low-temperature transition-region lines with He 1083 nm absorption to investigate Andretta's conjecture, i.e. that the He 1083 nm line is formed in two layers where extreme ultraviolet radiation produced both in the low-temperature transition region (the upper layer) and in the surrounding corona products - a lower layer of absorption in the upper chromosphere.

Description: The white light coronograph (WLC) on Skylab provided an opportunity to study the corona at high spatial and temporal resolution. The spatial resolution of the instrument was approximately 25 cm with images taken approximately one per min. One set of images taken over a 10 min period was digitized, providing ten high spatial resolution images for analysis. The progress in data processing techniques available at the time was not sufficient to permit a reliable study of the fine structure in these images. Using current techniques an investigation of the sizes and lifetimes of the smallest scale features in the data was carried out. A preliminary analysis of an area between 2 and 3 Ro was completed. The results show that very narrow rays extend from at least 2 to 3 Ro. The narrowest of these rays has a thickness of approximately 75 cm. The contrast is so low that they are very close to the noise limit of the data. Most of the rays observed become unrecognizable after 10 min, although some remain visible over the entire time. Some notion seems to be detectable in the fine structure rays, but analysis of more frames will be needed to quantify these results.

Description: Filaments, flare sprays, prominences and 'post-flare' loops are familiar to H alpha observers in their frequent appearances 'in absorption', dark against the chromospheric background or plages. Observations of the X-ray corona are generally interpreted as due to emission via optically thin thermal bremsstrahlung. Several cases of X-ray coronal structures in Yohkoh images, due to high opacity, absorbing matter in coronograph loops, are presented. The presence of the absorbing matter, mixed with emitting matter, complicates inference of physical parameters such as emission measures in X-ray sources. In the case of well defined features, absorption provides an opportunity to infer density. Quantitative estimates of the attenuation due to the absorption in example features are presented.

Description: The causal association of major solar particle events seen at earth with coronal mass ejections (CME's), and not with solar flares, is discussed. Evidence that led to the demise of the flare dominated paradigm for major solar energetic particle events are described. The possibility of distinguishing particles from impulsive and gradual events using only observations is described. Particle acceleration at the CME level is discussed. Multi-spacecraft observations of CME events are described. Concerning the interplanetary CME, bidirectional proton events are discussed. Conclusions from progress in understanding the characteristics of solar energetic particles and their relation to the physical mechanisms of acceleration are given.

Description: Two dimensional magnetohydrodynamic simulations of the distortion of a magnetic flux tube, accelerated through ambient solar wind plasma, are presented. Vortices form on the trailing edge of the flux tube, and couple strongly to its interior. If the flux tube azimuthal field is weak, it deforms into an elongated banana-like shape after a few Alfven transit times. A significant azimuthal field component inhibits this distortion. In the case of magnetic clouds in the solar wind, it is suggested that the shape observed at 1 AU was determined by distortion of the cloud in the inner heliosphere. Distortion of the cloud beyond 1 AU takes many days. It is estimated that effective drag coefficients slightly greater than unity are appropriate for modeling flux tube propagation. Synthetic magnetic field profiles as would be seen by a spacecraft traversing the cloud are presented.

Description: Since radio propagation measurements using either natural or spacecraft radio signals are used for probing the solar wind in the vicinity of the sun, they represent a key tool for studying the interplanetary consequences of solar structure and dynamic phenomena. New information on the near sun consequences was obtained from radio scintillation observations of coherent spacecraft signals. The results covering density fluctuations, fractional density fluctuations, coronal streamers, heliospheric current sheets, coronal mass ejections and interplanetary shocks are reviewed. A joint ICE S-band (13 cm wavelength) Doppler scintillation measurement with the SOHO white-light coronograph (LASCO) is described.

Description: An analytical solution to the perturbative multiple collision series of a fragmenting HZE ion beam has limited usefulness since the first collision term has several hundred contributions, the second collision term has tens of thousands of contributions, and each successive collision term progresses to unwieldy computational proportions. Our previous work has revealed the multiple collision terms in the straight-ahead approximation to be simple products of a spatially dependent factor times a linear energy-dependent factor of limited domain and unit normalization. The properties of these forms allow the development of the nonperturbative summation of the series to all orders assuming energy-independent nuclear cross sections as matrix products of a scaled Green's function described herein. This nonperturbative Green's function with multiple scattering correction factors compares well with experiments using 670 MeV/u neon-20 ion beams in thick water targets.

Description: Evidence now suggests that satellite cells constitute a class of myogenic cells that differ distinctly from other embryonic myoblasts. Satellite cells arise from somites and first appear as a distinct myoblast type well before birth. Satellite cells from different muscles cannot be functionally distinguished from one another and are able to provide nuclei to all fibers without regard to phenotype. Thus, it is difficult to ascribe any significant function to establishing or stabilizing fiber type, even during regeneration. Within a muscle, satellite cells exhibit marked heterogeneity with respect to their proliferative behavior. The satellite cell population on a fiber can be partitioned into those that function as stem cells and those which are readily available for fusion. Recent studies have shown that the cells are not simply spindle shaped, but are very diverse in their morphology and have multiple branches emanating from the poles of the cells. This finding is consistent with other studies indicating that the cells have the capacity for extensive migration within, and perhaps between, muscles. Complexity of cell shape usually reflects increased cytoplasmic volume and organelles including a well developed Golgi, and is usually associated with growing postnatal muscle or muscles undergoing some form of induced adaptive change or repair. The appearance of activated satellite cells suggests some function of the cells in the adaptive process through elaboration and secretion of a product. Significant advances have been made in determining the potential secretion products that satellite cells make. The manner in which satellite cell proliferative and fusion behavior is controlled has also been studied. There seems to be little doubt that cellcell coupling is not how satellite cells and myofibers communicate. Rather satellite cell regulation is through a number of potential growth factors that arise from a number of sources. Critical to the understanding of this form of control is to determine which of the many growth factors that can alter satellite cell behavior in vitro are at work in vivo. Little work has been done to determine what controls are at work after a regeneration response has been initiated. It seems likely that, after injury, growth factors are liberated through proteolytic activity and initiate an activation process whereby cells enter into a proliferative phase. After myofibers are formed, it also seems likely that satellite cell behavior is regulated through diffusible factors arising from the fibers rather than continuous control by circulating factors.(ABSTRACT TRUNCATED AT 400 WORDS).

Description: When PCR was used to recover small-subunit (SSU) rRNA genes from a hot spring cyanobacterial mat community, chimeric SSU rRNA sequences which exhibited little or no secondary structural abnormality were recovered. They were revealed as chimeras of SSU rRNA genes of uncultivated species through separate phylogenetic analysis of short sequence domains.

Description: We performed morphologic studies on skin from seven patients with ALS and seven control subjects. By light microscopy, the wide spaces that separated collagen bundles reacted strongly with colloidal iron and alcian blue in ALS patients. Electron microscopy revealed markedly increased amorphous material that was positive for ruthenium red in the ground substance. These findings were not present in controls. Quantitative amino acid analysis showed that the amount of total amino acids (nmoles per mg dry weight) was significantly decreased (p 〈 0.01) in ALS patients compared with that of controls, and there was a significant negative correlation between skin amino acid content and duration of illness in ALS patients (r = -0.83, p 〈 0.001). These morphologic findings and biochemical data indicate that the amorphous material, which is markedly increased in ALS skin, includes glycosaminoglycans.

Description: The present study assessed the contributions of feeding changes and unloading to the overall measured effects of 2-wk hindlimb (Tail) suspension on the mouse femora. Feeding changes were addressed by considering the effects of matched feeding among suspended and control mice. The effects of hind limb unloading were considered by comparing suspended mice to mice equipped identically (though not suspended) and matched-fed. The feeding and unloading aspects of suspension appear to cause distinctly differing effects on the stereotypic modeling of the femora. Matched-feeding was accompanied by increased resorption surface in comparison to suspended mice, while unloading led to reduced bone formation at the mid-diaphysis of the femora. Reduced mineral content was observed in the bones of suspended mice when compared to the other mice groups, but without increased resorption surface. Thus, the unloading aspects of the antiorthostatic suspension protocol apparently causes reduced formation and mineralization in the femur.

Description: Dorsal-ventral patterning in the Xenopus egg becomes established midway through the first cell cycle during a 30 degree rotation of the subcortical yolk mass relative to the egg cortex. This 'rotation of symmetrisation' is microtubule dependent, and its direction is thought to be cued by the usually eccentric sperm centrosome. The fact that parthenogenetically activated eggs also undergo a directed rotation, despite the absence of a sperm centrosome, suggests that an endogenous asymmetry in the unfertilised egg supports the directed polymerisation of microtubules in the vegetal cortex, in the way that an eccentric sperm centrosome would in fertilised eggs. Consistent with this idea, we noticed that the maturation spot is usually located an average of more than 15 degrees from the geometric centre of the pigmented animal hemisphere. In parthenogenetically activated eggs, this eccentric maturation spot can be used to predict the direction of rotation. Although in most fertilised eggs the yolk mass rotates toward the sperm entry point (SEP) meridian, occasionally this relationship is perturbed significantly; in such eggs, the maturation spot is never on the same side of the egg as the SEP. In oocytes tilted 90 degrees from upright during maturation in vitro, the maturation spot developed 15 degrees or more from the centre of the pigmented hemisphere, always displaced towards the point on the equator that was up during maturation. This experimentally demonstrated lability is consistent with an off-axis oocyte orientation during oogenesis determining its eccentric maturation spot position, and, in turn, its endogenous rotational bias.

Description: The mitogen phytohemagglutinin (PHA) works well in both human and cynomolgus monkey (Macaca fascicularis) lymphocyte cultures to stimulate T cell proliferation. T cells from rhesus monkeys (Macaca mulatta) are less responsive than human cells, producing few metaphases when thousands are required, e.g. in biological dosimetry studies. We show that staphylococcal enterotoxin A (SEA), one of the most potent mitogens known, at a concentration of 0.5 microgram/ml stimulated peripheral lymphocytes to grow with a mitotic index (MI) averaging 0.13 metaphases/cell in old, irradiated rhesus macaques. This was significantly greater (p 〈 0.001) than that produced by PHA (MI 〈 0.01) in lymphocytes from the same animals. Whole blood was cultured for 96, 120 and 144 h for five irradiated individuals and for two controls. All cells cultured with SEA produced a high MI with a peak response at 120 h whereas the same cultures showed low MI for each PHA stimulated culture.

Description: No abstract available

Description: The quality of volume-localized magnetic resonance spectroscopy is affected by eddy currents caused by gradient switching. Eddy currents can be reduced with improved gradient systems; however, it has been suggested that the distortion due to eddy currents can be compensated for during postprocessing with a single-frequency reference signal. The authors propose modifying current techniques for acquiring the single-frequency reference signal by using relaxation weighting to reduce interference from components that cannot be eliminated by digital filtering alone. Additional sequences with T1 or T2 weighting for reference signal acquisition are shown to have the same eddy current characteristics as the original signal without relaxation weighting. The authors also studied a new eddy current correction method that does not require a single-frequency reference signal. This method uses two free induction decays (FIDs) collected from the same volume with two sequences with opposite gradients. Phase errors caused by eddy currents are opposite in these two FIDs and can be canceled completely by combining the FIDs. These methods were tested in a phantom. Eddy current distortions were corrected, allowing quantitative measurement of structures such as the -CH = CH- component, which is otherwise undetectable.

Description: Quantitative assessment of the heritable somatic effects of ionizing radiation exposures has relied upon the assumption that radiation-induced lesions were 'fixed' in the DNA prior to the first postirradiation mitosis. Lesion conversion was thought to occur during the initial round of DNA replication or as a consequence of error-prone enzymatic processing of lesions. The standard experimental protocols for the assessment of a variety of radiation-induced endpoints (cell death, specific locus mutations, neoplastic transformation and chromosome aberrations) evaluate these various endpoints at a single snapshot in time. In contrast with the aforementioned approaches, some studies have specifically assessed radiation effects as a function of time following exposure. Evidence has accumulated in support of the hypothesis that radiation exposure induces a persistent destabilization of the genome. This instability has been observed as a delayed expression of lethal mutations, as an enhanced rate of accumulation of non-lethal heritable alterations, and as a progressive intraclonal chromosomal heterogeneity. The genetic controls and biochemical mechanisms underlying radiation-induced genomic instability have not yet been delineated. The aim is to integrate the accumulated evidence that suggests that radiation exposure has a persistent effect on the stability of the mammalian genome.

Description: While studies concerning mitogenic factors have been an important area of research for many years, much less is understood about the mechanisms of action of cell surface growth inhibitors. We have purified an 18 kDa cell surface sialoglycopeptide growth inhibitor (CeReS-18) which can reversibly inhibit the proliferation of diverse cell types. The studies discussed in this article show that three mouse keratinocyte cell lines exhibit sixty-fold greater sensitivity than other fibroblasts and epithelial-like cells to CeReS-18-induced growth inhibition. Growth inhibition induced by CeReS-18 treatment is a reversible process, and the three mouse keratinocyte cell lines exhibited either single or multiple cell cycle arrest points, although a predominantly G0/G1 cell cycle arrest point was exhibited in Swiss 3T3 fibroblasts. The sensitivity of the mouse keratinocyte cell lines to CeReS-18-induced growth inhibition was not affected by the degree of tumorigenic progression in the cell lines and was not due to differences in CeReS-18 binding affinity or number of cell surface receptors per cell. However, the sensitivity of both murine fibroblasts and keratinocytes could be altered by changing the extracellular calcium concentration, such that increased extracellular calcium concentrations resulted in decreased sensitivity to CeReS-18-induced proliferation inhibition. Thus the increased sensitivity of the murine keratinocyte cell lines to CeReS-18 could be ascribed to the low calcium concentration used in their propagation. Studies are currently under way investigating the role of calcium in CeReS-18-induced growth arrest. The CeReS-18 may serve as a very useful tool to study negative growth control and the signal transduction events associated with cell cycling.

Description: Space flight represents a form of sensory stimulus rearrangement requiring modification of established terrestrial response patterns through central reinterpretation. Evidence of sensory reinterpretation is manifested as postflight modifications of eye/head coordination, locomotor patterns, postural control strategies, and illusory perceptions of self or surround motion in conjunction with head movements. Under normal preflight conditions, the head is stabilized during locomotion, but immediately postflight reduced head stability, coupled with inappropriate eye/head coordination, results in modifications of gait. Postflight postural control exhibits increased dependence on vision which compensates for inappropriate interpretation of otolith and proprioceptive inputs. Eye movements compensatory for perceived self motion, rather than actual head movements have been observed postflight. Overall, the in-flight adaptive modification of head stabilization strategies, changes in head/eye coordination, illusionary motion, and postural control are maladaptive for a return to the terrestrial environment.

Description: No abstract available

Description: Bedrest studies of normal subjects provide opportunities to understand physiologic responses to supine posture and inactivity. Furthermore, head-down tilt has been a valuable procedure to investigate adaptation to microgravity and development of countermeasures to maintain the health and well-being of humans during space-flight. Recent bedrest experiments at NASA have ranged in duration from a few hours to 17 weeks. Acute studies of 6 degrees head-down tilt indicate that elevation of capillary blood pressure from 28 to 34 mm Hg and increased capillary perfusion in tissues of the head cause facial and intracranial edema. Intracranial pressure increases from 2 to 17 mm Hg going from upright posture to 6 degrees head-down tilt. Microvessels of the head have a low capacity to constrict and diminish local perfusion. Elevation of blood and tissue fluid pressures/flow in the head may also explain the higher headward bone density associated with long-term head-down tilt. These mechanistic studies of head-down tilt, along with a better understanding of the relative stresses involved with upright posture and lower body negative pressure, have facilitated development of suitable physiologic countermeasures to maintain astronaut health during microgravity. Presently no exercise hardware is available to provide a blood pressure gradient from head to feet in space. However, recent studies in our laboratory suggest that treadmill exercise using a graded lower-body compression suit and 100 mmHg lower body negative pressure provides equivalent or greater physiologic stress than similar upright exercise on Earth. Therefore, exercise within a lower body negative pressure chamber may provide a cost-effective and simple countermeasure to maintain the cardiovascular and neuro-musculoskeletal systems of astronauts during long-duration flight.

Description: There are types of deafness and tinnitus in which ruptures or massive changes in the ionic permeability of the membranes lining the endolymphatic space [e.g., of the reticular lamina (RL)] are believed to allow potassium-rich endolymph to deluge the low [K+] perilymphatic fluid (e.g., in the small spaces of Nuel). This would result in a K+ intoxication of sensory and neural structures. Acute attacks of Meniere's disease have been suggested to be an important example for this event. The present study investigated the effects of transiently elevated [K+] due to the addition of artificial endolymph to the basolateral cell surface of outer hair cells (OHC) in replicating endolymph-induced K+ intoxication of the perilymph in the small spaces of Nuel. The influence of K+ intoxication of the basolateral OHC cell surface on the transduction was then examined. Intoxication resulted in an inhibition of the physiological repolarizing K+ efflux from hair cells. This induced unwanted depolarizations of the hair cells, interfering with mechanoelectrical transduction. A pathological longitudinal OHC shortening was also found, with subsequent compression of the organ of Corti possibly influencing the micromechanics of the mechanically active OHC. Both micromechanical and electrophysiological alterations are proposed to contribute to endolymph leakage induced attacks of deafness and possibly also to tinnitus. Moreover, repeated or long-lasting K+ intoxications of OHC resulted in a chronic and complete loss of OHC motility. This is suggested to be a pathophysiological basis in some patients with chronic hearing loss resulting from Meniere's syndrome.

Description: No abstract available

Description: We have developed a computer monitor of nuclear DNA in the form of chromatin fibre. The fibres are modeled as a ideal solenoid consisting of twenty helical turns with six nucleosomes per turn. The chromatin model, in combination with are Monte Carlo theory of radiation damage induces by charged particles, based on general features of tack structure and stopping power theory, has been used to evaluate the influence of DNA structure on initial damage. An interesting has emerged from our calculations. Our calculated results predict the existence of strong spatial correlations in damage sites associated with the symmetries in the solenoidal model. We have calculated spectra of short fragments of double stranded DNA produced by multiple double strand breaks induced by both high and low LET radiation. The spectra exhibit peaks at multiples of approximately 85 base pairs (the nucleosome periodicity), and approximately 1000 base pairs (solenoid periodicity). Preliminary experiments to investigate the fragment distributions from irradiated DNA, made by B. Rydberg at Lawrence Berkeley Laboratory, confirm the existence of short DNA fragments and are in substantial agreement with the predictions of our theory.

Description: Citrate is an inhibitor of the crystallization of stone-forming calcium salts. Hypocitraturia, frequently encountered in patients with nephrolithiasis, is therefore an important risk factor for stone formation. Potassium citrate provides physiological and physicochemical correction and inhibits new stone formation, not only in hypocitraturic calcium nephrolithiasis but also in uric acid nephrolithiasis. Inhibition of stone recurrence has now been validated by a randomized trial. Ongoing research has disclosed additional causes of hypocitraturia (sodium excess, low intestinal alkali absorption, but not primary citrate malabsorption). Moreover, new insights on potassium citrate action have been shown, notably that some of absorbed citrate escapes oxidation and contributes to the citraturic response, that ingestion with a meal does not sacrifice physiological or physicochemical action, that orange juice mimics but does not completely duplicate its actions, that potassium citrate may have a beneficial bone-sparing effect, that it may reduce stone fragments following ESWL, and that danger of aluminum toxicity is not great in subjects with functioning kidneys. Finally, the research on potassium citrate has led to two promising products, calcium citrate as an optimum calcium supplement and potassium-magnesium citrate which may be superior to potassium citrate in the management of stone disease.

Description: Major strides in the molecular biology of essential hypertension are currently underway. This has tended to obscure the fact that a number of inherited disorders associated with low blood pressure exist and that these diseases may have milder and underrecognized phenotypes that contribute importantly to blood pressure variation in the general population. This review highlights some of the gene products that, if abnormal, could cause hypotension in some individuals. Diseases due to abnormalities in the catecholamine enzymes are discussed in detail. It is likely that genetic abnormalities with hypotensive phenotypes will be as interesting and diverse as those that give rise to hypertensive disorders.

Description: We recorded the horizontal (yaw), vertical (pitch), and torsional (roll) eye movements of two rhesus monkeys with scleral search coils before and after the COSMOS Biosatellite 2229 Flight. The aim was to determine effects of adaptation to microgravity on the vestibulo-ocular reflex (VOR). The animals flew for 11 days. The first postflight tests were 22 h and 55 h after landing, and testing extended for 11 days after reentry. There were four significant effects of spaceflight on functions related to spatial orientation: (1) Compensatory ocular counterrolling (OCR) was reduced by about 70% for static and dynamic head tilts with regard to gravity. The reduction in OCR persisted in the two animals throughout postflight testing. (2) The gain of the torsional component of the angular VOR (roll VOR) was decreased by 15% and 50% in the two animals over the same period. (3) An up-down asymmetry of nystagmus, present in the two monkeys before flight was reduced after exposure to microgravity. (4) The spatial orientation of velocity storage was shifted in the one monkey that could be tested soon after flight. Before flight, the yaw axis eigenvector of optokinetic afternystagmus was close to gravity when the animal was upright or tilted. After flight, the yaw orientation vector was shifted toward the body yaw axis. By 7 days after recovery, it had reverted to a gravitational orientation. We postulate that spaceflight causes changes in the vestibular system which reflect adaptation of spatial orientation from a gravitational to a body frame of reference. These changes are likely to play a role in the postural, locomotor, and gaze instability demonstrated on reentry after spaceflight.

Description: The high molecular weight aggregate (HMWA) fraction was isolated from the water soluble proteins of aged bovine lenses. Its composition and ability to inhibit heat-induced denaturation and aggregation were compared with the lower molecular weight, oligomeric fraction of alpha isolated from the same lens. Although the major components of both fractions were the alpha-A and alpha-B chains, the HMWA fraction possessed a decreased ability to protect other proteins against heat-induced denaturation and aggregation. Immunoelectron microscopy of both fractions demonstrated that alpha particles from the HMWA fraction contained increased amounts of beta and gamma crystallins, bound to a central region of the supramolecular complex. Together, these results demonstrate that alpha crystallins found in the HMWA fraction possess a decreased ability to protect against heat-induced denaturation and aggregation, and suggest that at least part of this decrease could be due to the increased presence of beta and gamma crystallins complexed to the putative chaperone receptor site of the alpha particles.

Description: No abstract available

Description: Experiments conducted on STS-50 indicated that space flight significantly inhibited tumor necrosis factor (TNF)-mediated killing of LM929 cells compared to ground controls. In ground-based studies, activation of protein kinase C (PKC) with phorbol 12-myristate 13-acetate (PMA) also inhibited TNF-mediated killing of LM929 cells. Therefore, we used PKC inhibitors to determine if the inhibitory effects of spaceflight on TNF-mediated cytotoxicity involved the activation of PKC. In experiments conducted onboard space shuttle mission STS-54, we saw that in the presence of the protein kinase C inhibitors H7 and H8, TNF-mediated cytotoxicity was restored to levels of those observed in the ground controls. Subsequent experiments done during the STS-57 mission tested the dose response of two protein kinase inhibitors, H7 and HA1004. We again saw that killing was restored in a dose-dependent manner, with inhibitor concentrations known to inhibit PKC being most effective. These data suggest that space flight ameliorates the action of TNF by affecting PKC in target cells.

Description: Focal sarcomere disruptions were previously observed in adductor longus muscles of rats flown approximately two weeks aboard the Cosmos 1887 and 2044 biosatellite flights. These lesions, characterized by breakage and loss of myofilaments and Z-line streaming, resembled damage induced by unaccustomed exercise that includes eccentric contractions in which muscles lengthen as they develop tension. We hypothesized that sarcomere lesions in atrophied muscles of space flow rats were not produced in microgravity by muscle unloading but resulted from muscle reloading upon re-exposure to terrestrial gravity. To test this hypothesis, we examined temporal changes in sarcomere integrity of adductor longus muscles from rats subjected to 12.5 days of hindlimb suspension unloading and subsequent reloading by return to vivarium cages for 0, 6, 12, or 48 hours of normal weightbearing. Our ultrastructural observations suggested that muscle unloading (0 h reloading) induced myofibril misalignment associated with myofiber atrophy. Muscle reloading for 6 hours induced focal sarcomere lesions in which cross striations were abnormally widened. Such lesions were electron lucent due to extensive myofilament loss. Lesions in reloaded muscles showed rapid restructuring. By 12 hours of reloading, lesions were moderately stained foci and by 48 hours darkly stained foci in which the pattern of cross striations was indistinct at the light and electron microscopic levels. These lesions were spanned by Z-line-like electron dense filamentous material. Our findings suggest a new role for Z-line streaming in lesion restructuring: rather than an antecedent to damage, this type of Z-line streaming may be indicative of rapid, early sarcomere repair.

Description: We assessed the reactivity of purified S-100 antiserum in immuno-electron microscopy by counting the number of gold particles per microns 2 over inner ear tissues embedded in different media. Sections containing predominantly Schwann's cell cytoplasm and nucleus, afferent fiber axoplasm and myelin sheath of chick cochleae were reacted with anti-S-100 IgG, an antibody to a calcium binding protein of neuronal tissues, then labeled with anti-IgG-gold conjugate. This investigation was conducted because previously published procedures, unmodified, did not yield acceptable results. Preparation of all specimens was identical. Only the medium (PolyBed 812, Araldite or Spurr epoxies; and LR White, LR Gold or Lowicryl plastics) was changed. The medium was made the changing variable because antigens available in post-embedding immuno-electron microscopy are decreased by heat, either used and/or released during polymerization of the embedding medium. The results indicate that: (a) none of the embedding media above provided optimal signal-to-noise ratio for all parts of the nerve stained in the same section; (b) aggregation of gold particles over cells was highest in embedding media with high background labeling over areas devoid of tissue (noise); (c) aggregation occurred randomly throughout both cellular and acellular regions; and (d) particles aggregated less and were distributed more evenly in tissues from media yielding good ultrastructural integrity.

Description: A theoretical model is developed to predict the fluid shear stress and streaming potential at the surface of osteocytic processes in the lacunar-canalicular porosity of an osteon when the osteon is subject to mechanical loads that are parallel or perpendicular to its axis. The theory developed in Weinbaum et al. (31) for the flow through a proteoglycan matrix in a canaliculus is employed in a poroelastic model for the osteon. Our formulation is a generalization of that of Petrov et al. (17). Our model predicts that, in order to satisfy the measured frequency dependence of the phase and magnitude of the SGP in macroscopic bone samples, the fiber spacing in the fluid annulus must lie in the narrow range 6-7 nm typical of the spacing of GAG sidechains along a protein monomer. The model predictions for the local SGP profiles in the osteon agree with the experimental observations of Starkebaum et al. (24). The theory predicts that the pore pressure relaxation time, tau d, for a 150-300 microns diameter osteon with the foregoing matrix structure is approximately 0.03-0.13 sec, and that the amplitude of the mean fluid shear stress on the membrane of the osteocytic process at the mean areal radius of the osteon has a maximum at 28 Hz if tau d = 0.06 sec. This maximum, which is independent of the magnitude of the loading, could be important in vivo since the recent experiments of Turner et al. (28) and McLeod et al. (15) have a peak in the strain frequency spectrum between 20 and 30 Hz that also appears to be independent of the type (magnitude) of loading. Numerical predictions for the amplitude of the average fluid shear stress on the osteocytic membrane at the mean areal radius of the osteon show that the fluid shear stress associated with the low amplitude 20-30 Hz spectral strain component is at least as large as the average fluid shear stress associated with the high amplitude 1 Hz stride component, although the latter loading is an order of magnitude larger, and has a magnitude that lies within the middle of the range, 6-30 dynes/cm2, where fluid shear stresses in tissue culture studies with osteoblast monolayers have elicited an intracellular Ca++ response (31). The implications of these results for intracellular electrical communication are discussed.

Description: Class I and class II major histocompatibility complex (MHC) molecules play significant roles in T cell development and immune function. We show that MHCI- and MHCII-deficient mice have low numbers of macrophage precursors and circulating monocytes, as well as abnormal bone marrow cell colony-stimulating factor type 1 secretion and bone composition. We suggest that MHCI and MHCII molecules play a significant role in macrophage development.

Description: Recently, it was observed that noncoding regions of DNA sequences possess long-range power-law correlations, whereas coding regions typically display only short-range correlations. We develop an algorithm based on this finding that enables investigators to perform a statistical analysis on long DNA sequences to locate possible coding regions. The algorithm is particularly successful in predicting the location of lengthy coding regions. For example, for the complete genome of yeast chromosome III (315,344 nucleotides), at least 82% of the predictions correspond to putative coding regions; the algorithm correctly identified all coding regions larger than 3000 nucleotides, 92% of coding regions between 2000 and 3000 nucleotides long, and 79% of coding regions between 1000 and 2000 nucleotides. The predictive ability of this new algorithm supports the claim that there is a fundamental difference in the correlation property between coding and noncoding sequences. This algorithm, which is not species-dependent, can be implemented with other techniques for rapidly and accurately locating relatively long coding regions in genomic sequences.

Description: Organ culture of embryonic mouse lung and pancreas rudiments has been used to investigate development and differentiation, and to assess the effects of microgravity on culture differentiation, during orbital spaceflight of the shuttle Endeavour (mission STS-54). Lung rudiments continue to grow and branch during spaceflight, an initial result that should allow future detailed study of lung morphogenesis in microgravity. Cultured embryonic pancreas undergoes characteristic exocrine acinar tissue and endocrine islet tissue differentiation during spaceflight, and in ground controls. The rudiments developing in the microgravity environment of spaceflight appear to grow larger than their ground counterparts, and they may have differentiated more rapidly than controls, as judged by exocrine zymogen granule presence.

Description: No abstract available

Description: Myoblast cell cultures have been widely employed in conventional (1g) studies of biological processes because characteristics of intact muscle can be readily observed in these cultured cells. We decided to investigate the effects of spaceflight on muscle by utilizing a well characterized myoblast cell line (L8 rat myoblasts) as cultured in the recently designed Space Tissue Loss Flight Module "A" (STL-A). The STL-A is a "state of the art," compact, fully contained, automated cell culture apparatus which replaces a single mid-deck locker on the Space Shuttle. The L8 cells were successfully flown in the STL-A on the Space Shuttle STS-45 mission. Upon return to earth, reculturing of these spaceflown L8 cells (L8SF) resulted in their unexpected failure to fuse and differentiate into myotubes. This inability of the L8SF cells to fuse was found to be a permanent phenotypic alteration. Scanning electron microscopic examination of L8SF cells growing at 1g on fibronectin-coated polypropylene fibers exhibited a strikingly different morphology as compared to control cells. In addition to their failure to fuse into myotubes, L8SF cells also piled up on top of each other. When assayed in fusion-promoting soft agar, L8SF cells gave rise to substantially more and larger colonies than did either preflight (L8AT) or ground control (L8GC) cells. All data to this point indicate that flying L8 rat myoblasts on the Space Shuttle for a duration of 7-10 d at subconfluent densities results in several permanent phenotypic alterations in these cells.

Description: Insulin-like growth factors (IGF) are important regulators of skeletal growth. To determine whether the capacity to produce and respond to these growth factors changes during skeletal development, we measured the protein and mRNA levels for IGF-I, IGF-II, and their receptors (IGF-IR and IGF-IIR, respectively) in the tibia and femur of rats before and up to 28 mo after birth. The mRNA levels remained high during fetal development but fell after birth, reaching a nadir by 3-6 wk. This fall was most pronounced for IGF-II and IGF-IIR mRNA and least pronounced for IGF-I mRNA. However, after 6 wk, both IGF-I and IGF-IR mRNA levels recovered toward the levels observed at birth. In the prenatal bones, the signals for the mRNAs of IGF-II and IGF-IIR were stronger than the signals for the mRNAs of IGF-I and IGF-IR, although the content of IGF-I was three- to fivefold greater than that of IGF-II. IGF-II levels fell postnatally, whereas the IGF-I content rose after birth such that the ratio IGF-I/IGF-II continued to increase with age. We conclude that, during development, rat bone changes its capacity to produce and respond to IGFs with a progressive trend toward the dominance of IGF-I.

Description: In eukaryotic cells, microtubules are 24-nm-diameter tubular structures composed of a class of conserved proteins called tubulin. They are involved in numerous cell functions including ciliary motility, nerve cell elongation, pigment migration, centrosome formation, and chromosome movement. Although cytoplasmic tubules and fibers have been observed in bacteria, some with diameters similar to those of eukaryotes, no homologies to eukaryotic microtubules have been established. Certain groups of bacteria including azotobacters, cyanobacteria, enteric bacteria, and spirochetes have been frequently observed to possess microtubule-like structures, and others, including archaebacteria, have been shown to be sensitive to drugs that inhibit the polymerization of microtubules. Although little biochemical or molecular biological information is available, the differences observed among these prokaryotic structures suggest that their composition generally differs among themselves as well as from that of eukaryotes. We review the distribution of cytoplasmic tubules in prokaryotes, even though, in all cases, their functions remain unknown. At least some tend to occur in cells that are large, elongate, and motile, suggesting that they may be involved in cytoskeletal functions, intracellular motility, or transport activities comparable to those performed by eukaryotic microtubules. In Escherichia coli, the FtsZ protein is associated with the formation of a ring in the division zone between the newly forming offspring cells. Like tubulin, FtsZ is a GTPase and shares with tubulin a 7-amino-acid motif, making it a promising candidate in which to seek the origin of tubulins.

Description: Current therapy for massive venous air embolism (VAE) includes the use of the left lateral recumbent (LLR) position. This recommendation is based on animal studies, conducted 50 yr ago, which looked primarily at survival. Little is known, however, about the concomitant hemodynamic response after VAE in various body positions. The purpose of this study was to investigate the hemodynamic and cardiovascular changes in various body positions after VAE. Twenty-two mechanically ventilated supine mongrel dogs received a venous air infusion of 2.5 mL/kg at a rate of 5 mL/s. One minute after the infusion, 100% oxygen ventilation was commenced and the body position of the dogs was changed to either the LLR (n = 6), the LLR with the head 10 degrees down (LLR-10 degrees; n = 6) or the right lateral recumbent (RLR; n = 5) position. Five dogs were maintained in the supine position (SUP; n = 5). One dog died in every group except in the SUP group, where all the dogs recovered. There were no significant differences among the various body positions in terms of heart rate, mean arterial pressure, pulmonary artery pressure, central venous pressure, left ventricular end-diastolic pressure, or cardiac output. The acute hemodynamic changes occurring during the first 5-15 min after VAE recovered to 80% of control within 60 min. Our data suggest that body repositioning does not influence the cardiovascular response to VAE. Specifically, our data do not support the recommendation of repositioning into the LLR position for the treatment of VAE.

Description: No abstract available

Description: This paper reviews a series of studies that indicate that estrogens play an important role in blood volume regulation. The first study illustrates that the plasma volume (PV) of ambulatory women fluctuates during the menstrual cycle, increasing during periods of elevated estrogens. In the second study, it was shown that exogenous and endogenous elevations in blood estrogens attenuate the decrease in PV during bed rest. In the third study, the hypothesis was tested that women, who naturally have a higher blood estrogen content compared with men, will have a smaller loss of PV during bed rest. Ten men and ten women underwent a 13-day, 6 degrees head-down bed rest. Plasma volume and red cell mass (RCM) were measured before and after bed rest using 125I and 51Cr labeling, respectively. Before bed rest, the men and women had similar blood volume (BV) and PV (mL/kg body weight), but the women had a smaller (P 〈 .01) RCM (22.2 +/- 0.9 versus 26.2 +/- 0.8 mL/kg, mean +/- SE). During bed rest, the decrease in RCM (mL/kg) was similar in men and women. However, the decrease in BV was greater in men (8.0 +/- 0.8 mL/kg versus 5.8 +/- 0.8 mL/kg), because of a greater reduction in PV (6.3 +/- 0.6 mL/kg versus 4.1 +/- 0.6 mL/kg). Because the decline in BV has been proposed to contribute to the cardiovascular deconditioning after bed rest, it is possible that women may experience less cardiac and circulatory strain on reambulation.

Description: Life sciences activities being planned for Space Station Freedom (SSF) as of Fall 1992 are discussed. Planning for these activities is ongoing. Therefore, this description should be viewed as indicative of the prevailing ideas at one particular time in the SSF development cycle. The proposed contributions of the Canadian Space Agency (CSN) the European Space Agency (ESA), Japan, and the United States are all discussed in detail. In each case, the life sciences goals, and the way in which each partner proposes to achieve their goals, are reviewed.

Description: This paper reviews the cardiovascular responses of six healthy male subjects to 6 hours in a 5 degrees head-down bed rest model of weightlessness, and compares these responses to those obtained when subjects were positioned in head-up tilts of 10 degrees, 20 degrees, and 42 degrees, simulating 1/6, 1/3, and 2/3 G, respectively. Thoracic fluid index, cardiac output, stroke volume, and peak flow were measured using impedance cardiography. Cardiac dimensions and volumes were determined from two-dimensional guided M-mode echocardiograms in the left lateral decubitus position at 0, 2, 4, and 6 hours. Cardiovascular response to a stand test were compared before and after bed rest. The impedance values were related to tilt angle for the first 2 hours of tilt; however, after 3 hours, at all four angles, values began to converge, indicating that cardiovascular homeostatic mechanisms seek a common adapted state, regardless of effective gravity level (tilt angle) up to 2/3 G. Echocardiography revealed that left ventricular end-diastolic and end-systolic volume, stroke volume, ejection fraction, heart rate, and cardiac output had returned to control values by hour 6 for all tilt angles. The lack of a significant immediate change in left ventricular end-diastolic volume, despite decrements in stroke volume (P 〈 .05) and heart rate (not significant), indicates that multiple factors may play a role in the adaptation to simulated hypogravity. The echocardiography data indicated that no angle of tilt, whether head-down or head-up for 4 to 6 hours, mimicked exactly the changes in cardiovascular function recorded after 4 to 6 hours of space flight. Changes in left ventricular end-diastolic volume during space flight and tilt may be similar, but follow a different time course. Nevertheless, head-down tilt at 5 degrees for 6 hours mimics some (stroke volume, systolic and diastolic blood pressure, mean arterial blood pressure, and total resistance), but not all, of the changes occurring in an equivalent time of space flight. The magnitude of the change in the mean heart rate response to standing was greater after six hours of tilt at -5 degrees or 10 degrees. Thus, results from the stand test after 6 hours of bed rest at -5 degrees and 10 degrees, but not at 20 degrees or 42 degrees, are similar to those obtained after space flight.

Description: No abstract available

Description: Microbalanced stimuli are dynamic displays which do not stimulate motion mechanisms that apply standard (Fourier-energy or autocorrelational) motion analysis directly to the visual signal. In order to extract motion information from microbalanced stimuli, Chubb and Sperling [(1988) Journal of the Optical Society of America, 5, 1986-2006] proposed that the human visual system performs a rectifying transformation on the visual signal prior to standard motion analysis. The current research employs two novel types of microbalanced stimuli: half-wave stimuli preserve motion information following half-wave rectification (with a threshold) but lose motion information following full-wave rectification; full-wave stimuli preserve motion information following full-wave rectification but lose motion information following half-wave rectification. Additionally, Fourier stimuli, ordinary square-wave gratings, were used to stimulate standard motion mechanisms. Psychometric functions (direction discrimination vs stimulus contrast) were obtained for each type of stimulus when presented alone, and when masked by each of the other stimuli (presented as moving masks and also as nonmoving, counterphase-flickering masks). RESULTS: given sufficient contrast, all three types of stimulus convey motion. However, only one-third of the population can perceive the motion of the half-wave stimulus. Observers are able to process the motion information contained in the Fourier stimulus slightly more efficiently than the information in the full-wave stimulus but are much less efficient in processing half-wave motion information. Moving masks are more effective than counterphase masks at hampering direction discrimination, indicating that some of the masking effect is interference between motion mechanisms, and some occurs at earlier stages. When either full-wave and Fourier or half-wave and Fourier gratings are presented simultaneously, there is a wide range of relative contrasts within which the motion directions of both gratings are easily determinable. Conversely, when half-wave and full-wave gratings are combined, the direction of only one of these gratings can be determined with high accuracy. CONCLUSIONS: the results indicate that three motion computations are carried out, any two in parallel: one standard ("first order") and two non-Fourier ("second-order") computations that employ full-wave and half-wave rectification.

Description: The 3-dimensional crystal structure of glutathione S-transferase (GST) of Schistosoma japonicum (Sj) fused with a conserved neutralizing epitope on gp41 (glycoprotein, 41 kDa) of human immunodeficiency virus type 1 (HIV-1) (Muster T et al., 1993, J Virol 67:6642-6647) was determined at 2.5 A resolution. The structure of the 3-3 isozyme rat GST of the mu gene class (Ji X, Zhang P, Armstrong RN, Gilliland GL, 1992, Biochemistry 31:10169-10184) was used as a molecular replacement model. The structure consists of a 4-stranded beta-sheet and 3 alpha-helices in domain 1 and 5 alpha-helices in domain 2. The space group of the Sj GST crystal is P4(3)2(1)2, with unit cell dimensions of a = b = 94.7 A, and c = 58.1 A. The crystal has 1 GST monomer per asymmetric unit, and 2 monomers that form an active dimer are related by crystallographic 2-fold symmetry. In the binding site, the ordered structure of reduced glutathione is observed. The gp41 peptide (Glu-Leu-Asp-Lys-Trp-Ala) fused to the C-terminus of Sj GST forms a loop stabilized by symmetry-related GSTs. The Sj GST structure is compared with previously determined GST structures of mammalian gene classes mu, alpha, and pi. Conserved amino acid residues among the 4 GSTs that are important for hydrophobic and hydrophilic interactions for dimer association and glutathione binding are discussed.

Description: Allan Hills (ALH) 84001 is the most recently recognized member of a suite of meteorites--the SNCs--that almost certainly originated on Mars. Several factors distinguish ALH84001 from the other SNC meteorites. Preliminary studies suggest that it may be older than other martian meteorites. Moreover, it contains abundant, zoned domains of calcium-iron-magnesium carbonate that are indigenous to the sample and thus may hold important clues regarding near-surface processes on Mars and the evolution of the martian atmosphere. We report here analyses of the carbon and oxygen stable-isotope compositions of the carbonates that place constraints on their formation conditions. Our results imply the presence of at least two chemically distinct carbonates--one Ca,Fe-rich, the other Mg-rich--that are enriched in 13C relative to terrestrial carbonates (delta 13C approximately +41/1000), consistent with martian atmospheric CO2 as the carbon source. The oxygen isotope compositions of the carbonates indicate that they precipitated from a low-temperature fluid in the martian crust. Combined with textural and bulk geochemical considerations, the isotope data suggest that carbonate deposition took place in an open-system environment in which the ambient temperature fluctuated.

Description: A repair/misrepair kinetic model for multiple radiation-induced lesions (mutation inactivation) is coupled to a two-mutation model of initiation-promotion in tissue to provide a parametric description of tumour prevalence in the mouse Harderian gland from high-energy and charge radiations. Track-structure effects are considered using an action-cross section model. Dose-response curves are described for gamma rays and relativistic ions, and good agreement with experiment is found. The effects of nuclear fragmentation are also considered for high-energy proton and alpha-particle exposures. The model described provides a parametric description of age-dependent cancer induction for a wide range of radiation fields. Radiosensitivity parameters found in the model for an initiation mutation (sigma 0 = 7.6 x 10(-10) cm2 and D0 = 148.0 Gy) are somewhat different than previously observed for neoplastic transformation of C3H10T1/2 cell cultures (sigma 0 = 0.7 x 10(-10) cm2 and D0 = 117.0 Gy). We consider the two hypotheses that radiation acts solely as an initiator or as both initiator and promoter and make model calculations for fractionation exposures from gamma rays and relativistic Fe ions. For fractionated Fe exposures, an inverse-dose-rate effect is provided by a promotion hypothesis with an increase of 30% or more, dependent on the dose level and fractionation schedule, using a mutation rate for promotion similar to that of single-gene mutations.

Description: The production of protons in heavy ion collisions through the knockout mechanism (abrasion) is described using the Glauber model. The multiple knockouts from the projectile, including the inelastic collision series with the target, are considered using a closure approximation in treating energy conservation. Calculations for reactions of 12C projectiles with several targets at energies of 1 and 2A GeV arc compared to experiments. For large secondary proton momentum a strong dependence on the target mass is found and attributed to multiple scattering of the projectile knockouts.

Description: The Galileo Ultravilet Spectrometer Experiment (UVS) obtained a partial celestial sphere map of interplanetary Lyman-alpha (IP L alpha) on 13-14 December 1990 during the first Earth encounter. The Galileo spacecraft was near the downwind axis of the local interstellar medium flow. These UVS measurements sampled the downwind, anti-sunward hemisphere. The data were modeled using a hot model of the interplanetary hydrogen density distribution with the goal of studying multiple scattering effects in the inner solar system. The derived ratio in the downwind direction of the observed brightness and a single scattering model brightness, both normalized to unity in the upwind direction, is 1.82 +/- 0.2. This brightness ratio requires a multiple scattering correction which is 36% larger than can be accounted for by theoretical calculations. The hot model may require: (1) a temperature perturbation of the interstellar wind velocity distribution or (2) an additional downstream source of interplanetary hydrogen. However, a more likely exlanation which affects the hot model is the latitude dependence of the radiation pressure. This dependence, based on the known solar L alpha flux latitude variation at solar maximum, causes a downwind brightness enhancement by preferential focusing of H-atoms with trajectory planes containing the solar poles. This result implies that radiation pressure near the solar poles is nearly independent of solar cycle and is insufficient to lead to a net repulsion of hydrogen atoms by the sun, as can occur near the ecliptic plane during the solar maximum. In addition, the UVS performed 13 observations of IP L alpha while in cruise between Venus and the Earth in 3 directions fixed in ecliptic coordinates.

Description: In this paper we generalize earlier gasdynamic analyses of the motion of the heliospheric termination shock in response to upstream disturbances (Barnes, 1993, 1994; Naidu and Barnes, 1994), to include magnetohydrodynamic (MHD) phenomena. We assume that the termination shock is a strong, perpendicular shock and that the initial upstream disturbance is a tangential discontinuity. The resulting configuration after the interaction is very similar to that in the gasdynamic models after an interaction with a contact discontinuity or interplanetary shock, and for an increase (decrease) in dynamic pressure consists of an outward (inward) propagating termination shock and an outward propagating shock (MHD rarefraction wave) that carries the signal of the disturbance into the far downstream plasma. The plasma immediately behind the new termination shock is separated from the downstream signal by a tangential discontinuity. The results of the model show that the speed of the new termination shock depends mainly on the magnitude of the change in dynamic pressure and are typically of order approximately 100 km/s, comparable to the results of the gasdynamic models.

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: The term 'magnetic hole' has been used to denote isolated intervals when the magnitude of the interplanetary magnetic field drops to a few tenths, or less, of its ambient value for a time that corresponds to a linear dimension of tens to a few hundreds of proton gyro-radii. Data obtained by the Ulysses magnetometer and solar wind anlayzer have been combined to study the properties of such magnetic holes in the solar wind between 1 AU and 5.4 AU and to 23 deg south latitude. In order to avoid confusion with decreases in field strength at interplanetary discontinuities, the study has focused on linear holes across which the field direction changed by less than 5 deg. The holes occurred preferentially, but not without exception, in the interaction regions on the leading edges of high-speed solar wind streams. Although the plasma surrounding the holes was generally stable against the mirror instability, there are indications that the holes may have been remnants of mirror-mode structures created upstream of the points of observation. Those indications include the following: (1) For the few holes for which proton of alpha-particle pressure could be measured inside the hole, the ion thermal pressure was always greater than in the plasma adjacent to the holes. (2) The plasma surrounding many of the holes was marginally stable for the mirror mode, while the plasma environment of all holes was significantly closer to mirror instability than was the average solar wind. (3) The plasma containing trains of closely spaced holes was closer to mirror instability than was the plasma containing isolated holes. (4) The near-hole plasma had much higher ion beta (ratio of thermal to magnetic pressure) than did the average solar wind. (5) Near the holes, T(sub perp)/T(sub parallel) tended to be either greater than 1 or larger than in the average wind. (6) The proton and alpha-particle distribution functions measured inside the holes occasionally exhibited the flattened phase-space-density contoures in nu(sub perp)/nu(sub parallel) space found in some numerical simulations of the mirror instability.

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: During the international campaign of June 1991, the active region AR 6659 produced six very large, long-duration flares (X10/12) during its passage across the solar disk. We present the characteristics of four of them (June 4, 6, 9, 15). Precise measurements of the spot motions from Debrecen and Tokyo white-light pictures are used to understand the fragmentation of the main sunspot group with time. This fragmentation leads to a continuous restructuring of the magnetic field pattern while rapid changes are evidenced due to fast new flux emergence (magnetograms of Marshall Space Flight Center (MSFC), Huairou). The first process leads to a shearing of the field lines along which there is energy storage; the second one is the trigger which causes the release of energy by creating a complex topology. We conjecture that these two processes with different time scales are relevant to the production of flares.

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: During the late declining phase of the solar cycle, the tilt of the solar magnetic dipole with respect to the Sun's rotation axis leads to large-scale organization of the solar wind, such that alternating regions of high- and low-speed solar wind are observed in the ecliptic plane. In this paper, we use Doppler scintillation measurements to investigate mass flux of these two types of solar wind in the ecliptic plane and inside 0.3 AU, where in situ measurements have not been possible. To the extent that Doppler scintillation reflects mass flux, we find that mass flux in high-speed streams: (1) is lower (by a factor of approximately 2.2) than the mass flux of the average solar wind in the heliocentric distance range of 0.3-0.5 AU; (2) is lower still (by as much as a factor of about 4) than the mass flux of the slow solar wind associated with the streamer belt; and (3) appears to grow with heliocentric distance. These Doppler scintillation results are consistent with the equator to pole decrease in mass flux observed in earlier spectral broadening measurements, and with trends and differences between high- and low-speed solar wind observed by in situ measurements in the range of 0.3-0.1 AU. The mass flux results suggest that the solar wind flow in high-speed streams is convergent towards the ecliptic near the Sun, becoming less convergent and approaching radial with increasing heliocentric distance beyond 0.3 AU. The variability of mass flux observed within equatorial and polar high-speed streams close to the Sun is strikingly low. This low variability implies that, as Ulysses currently ascends to higher latitudes and spends more time in the south polar high-speed stream after crossing the heliocentric current sheet, it can expect to observe a marked decrease in variations of both mass flux and solar wind speed, a trend that appears to have started already.

Description: Highly structured Langmuir waves, also known as electron plasma oscillations, have been observed in the foreshock of Venus using the plasma wave experiment on the Galileo spacecraft during the gravity assist flyby on February 10, 1990. The Galileo wideband sampling system provides digital electric field waveform measurements at sampling rates up to 201,600 samples per second, much higher than any previous instrument of this type. The main Langmuir wave emission band occurs near the local electron plasma frequency, which was approximately 43 kHz. The Langmuir waves are observed to shift above and below the plasma frequency, sometimes by as much as 20 kHz. The shifts in frequency are closely correlated with the downstream distance from the tangent field line, implying that the shifts are controlled by the electron beam velocity. Considerable fine structure is also evident, with timescales as short as 0.15 ms, corresponding to spatial scales of a few tens of Debye lengths. The frequency spectrum often consists of beat-type waveforms, with beat frequencies ranging from 0.2 to 7 kHz, and in a few cases, isolated wave packets. The peak electric field strengths are approximately 1 mV/m. These field strengths are too small for strongly nonlinear processes to be important. The beat-type waveforms are suggestive of a parametric decay process.

Description: The irregular polygonal pattern of solar granulation is analyzed for size-shape relations using statistical crystallography. In contrast to previous work which has assumed perfectly hexagonal patterns for granulation, more realistic accounting of cell (granule) shapes reveals a broader basis for quantitative analysis. Several features emerge as noteworthy: (1) a linear correlation between number of cell-sides and neighboring shapes (called Aboav-Weaire's law); (2) a linear correlation between both average cell area and perimeter and the number of cell-sides (called Lewis's law and a perimeter law, respectively) and (3) a linear correlation between cell area and squared perimeter (called convolution index). This statistical picture of granulation is consistent with a finding of no correlation in cell shapes beyond nearest neighbors. A comparative calculation between existing model predictions taken from luminosity data and the present analysis shows substantial agreements for cell-size distributions. A model for understanding grain lifetimes is proposed which links convective times to cell shape using crystallographic results.

Description: Dual-frequency ranging and Doppler measurements were conducted in support of the Ulysses Solar Corona Experiment (SCE) at and around the spacecraft's first solar conjunction in 1991 August. The differential group delay time between range codes on the two downlink carrier signals at the wavelengths 13.1 and 3.6 cm, a direct measure of the total electron content between spacecraft and ground station, was used to derive the electron density distribution in the solar corona. Linear power-law representations of the coronal electron density were derived for the range of solar distances from 4 solar radii to 40 solar radii on both sides of the Sun. The corona was found to be very nearly symmetric; the radial falloff exponent being 2.54 +/- 0.05 for occultation ingress (east solar limb) and 2.42 +/- 0.05 for egress (west limb), respectively. The departure of these exponents from the inverse equare relation implies that significant solar wind acceleration is occurring within the radial range of the observations. The electron density level was found to be considerably lower than that observed during the 1988 December solar occultation of Voyager 2. Although the smoothed sunspot number R(sub z) (a standard indicator of solar activity) was almost the same in 1988 December and 1991 August, the mean electron density at 20 solar radii was found to be 1.7 +/- 0.1 x 10(exp 3)/cu cm during the Ulysses conjunction, a decline by almost a factor of 4 from the value obtained during the Voyager conjunction.

Description: One of the long-standing uncertainties in the wave-resonance theory of coronal heating is the stability of the resonance layer. The wave motions in the resonance layer produce highly localized shear flows which vary sinusoidally in time with the resonance period. This configuration is potentially susceptible to the Kelvin-Helmholtz instability (KHI), which can enhance small-scale structure and turbulent broadening of shear layers on relatively rapid ideal timescales. We have investigated numerically the response of a characteristic velocity profile, derived from resonance absorption models, to finite fluid perturbations comparable to photospheric fluctuations. We find that the KHI primarily should affect long (approximately greater than 6 x 10(exp 4) km) loops where higher velocity flows (M approximately greater than 0.2) exist in resonance layers of order 100 km wide. There, the Kelvin-Helmholtz growth time is comparable to or less than the resonance quarter-period, and the potentially stabilizing magnetic effects are not felt until the instability is well past the linear growth stage. Not only is the resonance layer broadened by the KHI, but also the convective energy transport out of the resonance layer is increased, thus adding to the efficiency of the wave-resonance heating process. In shorter loops, e.g., those in bright points and compact flares, the stabilization due to the magnetic field and the high resonance frequency inhibit the growth of the Kelvin-Helmholtz instability beyond a minimal level.

Description: The solar flare and coronal mass ejection (CME) events associated with the large and complex March 1989 active region are discussed. This active region gave us a chance to study the relation of CME with truly major solar flares. The work concentrates on questions of the relation of CMEs and flares to one another and to other types of activity on the Sun. As expected, some major (X-3B class) flares had associated CMEs. However, an unexpected finding is that others did not. In fact, there is strong evidence that the X4-4B flare of March 9th had no CME. This lack of a CME for such an outstanding flare event has important implications to theories of CME causation.Apparently, not all major flares cause CMEs or are caused by CMEs. The relations between CMEs and other types of solar activity are also discussed. No filament disappearances are reported for major CMEs studied here. Comparing these results with other studies, CMEs occur in association with flares and with erupting prominences, but neither are required for a CME. The relation between solar structures showing flaring without filament eruptions and structures showing filament eruptions without flares becomes important. The evolutionary relation between an active flaring sunspot region and extensive filaments without sunspots is reviewed, and the concept of an 'evolving magnetic structure' (EMS) is introduced. It is suggested that all CMEs arise in EMSs and that CMEs provide a major path through which azimuthal magnetic fields escape form the Sun during the solar cycle.

Description: The large gradual solar-energetic-particle (SEP) events, where abundances are commonly measured, are produced when coronal mass ejections (CMEs) drive shock waves through the corona and the interplanetary medium. The shock accelerates particles from the highly-ionized, approximately 1.5 MK, plasma in a manner that depends only weakly upon the Q/A of the ion, except at very high energies. Averaging the approximately 1 MeV/amu abundances over many events compensates for the acceleration effects to produce abundances that appear to correspond directly to those in the coronal source for all observed elements, including H. The resulting abundances reflect the 4 x enhancement of ions with low values of first ionization potential (FIP) arising from ion-neutral fractionation that occurs as the atoms are transported up from the photosphere. A different pattern of fractionation is found for ions that are shock-accelerated from the high speed solar wind emerging from coronal holes.

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: 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: 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: 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 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 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: 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.