Search NASA⌕ Search

SEARCH · Search NASA

Results for “OIS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

Excitation of the arched filaments near the Galactic Center

We discuss measurements of the far-infrared (FIR) fine structure lines from (S III) (33 microns), (Si II) (35 microns), (O III) (51, 88 microns), (OI) (63 microns), (C II) (158 microns), and the adjacent continua in a strip crossing two of the thermal radio filaments in the Galactic Center 'Arch'. The near spatial coincidence of the line and continuum emission maxima with the radio filaments demonstrates that any excitation mechanism must account for both the line and continuum emission. The peak FIR luminosity and (O III) emission pose difficulties for collisional excitation models; photoionization of molecular cloud edges by a random distribution of stars is the most plausible mechanism proposed.

Colgan, Sean W. J.↗

Turbulence Scales, Rise Times, Caustics, and the Simulation of Sonic Boom Propagation

The general topic of atmospheric turbulence effects on sonic boom propagation is addressed with especial emphasis on taking proper and efficient account of the contributions of the portion oi the turbulence that is associated with extremely high wavenumber components. The recent work reported by Bart Lipkens in his doctoral thesis is reexamined to determine whether the good agreement between his measured rise times with the 1971 theory of the author is fortuitous. It is argued that Lipken's estimate of the distance to the first caustic was a gross overestimate because of the use of a sound speed correlation function shaped like a gaussian curve. In particular, it is argued that the expected distance to the first caustic varies with the kinematic viscosity nu and the energy epsilon dissipated per unit mass per unit time, and the sound speed c as : d(sub first caustic) = nu(exp 7/12) c(exp 2/3)/ epsilon(exp 5/12)(nu x epsilon/c(exp 4))(exp a), where the exponent a is greater than -7/12 and can be argued to be either O or 1/24. In any event, the surprising aspect of the relationship is that it actually goes to zero as the viscosity goes to zero with s held constant. It is argued that the apparent overabundance of caustics can be grossly reduced by a general computational and analytical perspective that partitions the turbulence into two parts, divided by a wavenumber k(sub c). Wavenumbers higher than kc correspond to small-scale turbulence, and the associated turbulence can be taken into account by a renormalization of the ambient sound speed so that the result has a small frequency dependence that results from a spatial averaging over of the smaller-scale turbulent fluctuations. Selection of k(sub c). can be made so large that only a very small number of caustics are encountered if one adopts the premise that the frequency dispersion of pulses is caused by that part of the turbulence spectrum which lies in the inertial range originally predicted by Kolmogoroff. The acoustic propagating wave's dispersion relation has the acoustic wavenumber being of the form k = (omega/c) + F(omega), where c is a spatially averaged sound speed and where, for mechanical turbulence, the extra term F(omega) must depend on only the angular frequency omega, the sound speed c, and the turbulent energy dissipation epsilon per unit fluid mass and per unit time. If the turbulence is weak, then the quantity F(omega) has to be of second order in the portions of the turbulent fluid velocity in the inertial range, so, following Kolmogoroff's reasoning, it must vary with epsilon as epsilon(exp 2/3). Simple dimensional analysis then reveals that F(omega) is K epsilon(exp 2/3) c(exp -7/3) omega(exp l/3), K being a universal dimensionless complex constant.

Pierce, Allan D.↗

Thermospheric O/N2 Based on DE-1 FUV Dayglow Imaging Data

Work performed during the second half of Year 1 of the contract is summarized. The primary objective of the work is to derive global dayside thermospheric oxygen concentrations from DE-1 far ultraviolet imaging data which we are considering under both magnetically quiet and disturbed times. Work to date has been more qualitative in producing maps showing intensity variations beyond those that can be explained by changes in solar zenith angle (SZA) and look angle across an image (Craven et al., 1995; Meier et al., 1994; Gladstone, 1994). In meeting our primary objective, four tasks have been addressed during the reporting period: (1) Investigating the uniqueness of the relationship between the dayglow emission seen using DE-1's 123 filter (dominated by OI 130.4 nm emission) and the column abundance of O relative to N2 referenced to an N2 depth of 10(exp 17) cm(exp 2); (2) Completion of the algorithm for rapid conversion of DE-1 disk dayglow measurements to O/N2 values; (3) Applying the algorithm to a simulation in which a model DE image was constructed using TIGCM atmospheres. The retrieved image of O/N2 was compared to TIGCM O/N2 obtained by integrations of the TIGCM densities; and (4) Applying the algorithm to selected DE-1 data.

Strickland, Douglas J.↗

Upper Thermosphere Winds and Temperatures in the Geomagnetic Polar Cap: Solar Cycle, Geomagnetic Activity, and Interplanetary Magnetic Field Dependencies

Ground-based Fabry-Perot interferometers located at Thule, Greenland (76.5 deg. N, 69.0 deg. W, lambda = 86 deg.) and at Sondre Stromfjord, Greenland (67.0 deg. N, 50.9 deg. W, lambda = 74 deg.) have monitored the upper thermospheric (approx. 240-km altitude) neutral wind and temperature over the northern hemisphere geomagnetic polar cap since 1983 and 1985, respectively. The thermospheric observations are obtained by determining the Doppler characteristics of the (OI) 15,867-K (630.0-nm) emission of atomic oxygen. The instruments operate on a routine, automatic, (mostly) untended basis during the winter observing seasons, with data coverage limited only by cloud cover and (occasional) instrument failures. This unique database of geomagnetic polar cap measurements now extends over the complete range of solar activity. We present an analysis of the measurements made between 1985 (near solar minimum) and 1991 (near solar maximum), as part of a long-term study of geomagnetic polar cap thermospheric climatology. The measurements from a total of 902 nights of observations are compared with the predictions of two semiempirical models: the Vector Spherical Harmonic (VSH) model of Killeen et al. (1987) and the Horizontal Wind Model (HWM) of Hedin et al. (1991). The results are also analyzed using calculations of thermospheric momentum forcing terms from the Thermosphere-ionosphere General Circulation Model TGCM) of the National Center for Atmospheric Research (NCAR). The experimental results show that upper thermospheric winds in the geomagnetic polar cap have a fundamental diurnal character, with typical wind speeds of about 200 m/s at solar minimum, rising to up to about 800 m/s at solar maximum, depending on geomagnetic activity level. These winds generally blow in the antisunward direction, but are interrupted by episodes of modified wind velocity and altered direction often associated with changes in the orientation of the Interplanetary Magnetic Field (IMF). The central polar cap (greater than approx. 80 magnetic latitude) antisunward wind speed is found to be a strong function of both solar and geomagnetic activity. The polar cap temperatures show variations in both solar and geomagnetic activity, with temperatures near 800 K for low K(sub p) and F(sub 10.7) and greater than about 2000 K for high K(sub p) and F(sub 10.7). The observed temperatures are significantly greater than those predicted by the mass spectrometer/incoherent scatter model for high activity conditions. Theoretical analysis based on the NCAR TIGCM indicates that the antisunward upper thermospheric winds, driven by upstream ion drag, basically 'coast' across the polar cap. The relatively small changes in wind velocity and direction within the polar cap are induced by a combination of forcing terms of commensurate magnitude, including the nonlinear advection term, the Coriolis term, and the pressure gradient force term. The polar cap thennospheric thermal balance is dominated by horizontal advection, and adiabatic and thermal conduction terms.

Killeen, T. L.↗

Nocturnal Observations of the Semidiurnal Tide at a Midlatitude Site

Fabry-Perot interferometer observations of the mesospheric hydroxyl emission and the lower thermospheric OI (5577A) emission have been conducted from an airglow observatory at a dark field site in southeastern Michigan for the past several years. The primary functions of the observatory are to provide a database for correlative observations with the UARS satellite and to provide a synoptic measurement program for the coupling energetics and dynamics of atmospheric regions effort, An intensive operational effort between May 1993 and July 1994 has resulted in a substantial data set from which neutral winds have been determined from the bifilter acquisition sequence. A 'best fit' analysis in the least squares sense of the simultaneous measurements of the neutral winds to a 12-hour periodicity has provided amplitude and phase parameters for the semidiurnal tide as well as a measure of the mean wind. The measured tidal amplitude is greater at the higher altitude, though the seasonal behavior at both altitudes is similar with greater amplitudes during August/September and April/May. Both meridional and zonal wind components are consistent with a semidiurnal tidal description during the entire observational sequence except for the May to July 1993 period. The mean winds show annual variation in the meridional flow, being equatorward from May to October and poleward during the winter. The zonal flow is primarily eastward during the entire observational window with higher speed flows during May/June at the higher attitude and June/July at the lower altitude. A comparison with a semidiurnal tidal model indicates that the measured tidal amplitudes are a factor of 2 times greater, while the phases show similar equinoctial transitions.

Niciejewski, R. J.↗

Studies for the Loss of Atomic and Molecular Species from Io

The general objective of this project has been to advance our theoretical understanding of Io's atmosphere and how various atomic and molecular species are lost from this atmosphere and are distributed in the circumplanetary environment of Jupiter. The scientific objectives of the larger collaborative program between AER, Inc., and the University of Michigan have been to undertake theoretical modeling studies to simulate the distributions of the exospheric gases in Io's corona and extended clouds, to investigate the importance of the various physical processes that shape their relative abundances, and with these tools to analyze observations of O, S and Na obtained by four observers: M.A. McGrath of the Space Telescope Science Institute and G.E. Ballester of the University of Michigan who each have obtained Hubble Space Telescope observations of O and S near Io, F. Scherb who continues an effort to obtain 6300 A OI observations as part of the University of Wisconsin Fabry-Perot program, and N.M. Schneider of the University of Colorado who obtained an extensive set of spectral and spatial observations of the Na emission near Io in the D-lines.

Combi, Michael R.↗

A Survey of Large-Scale Variations in Thermospheric Oxygen Column Density with Magnetic Activity as Inferred from Observations of the FUV Dayglow

Brightness of the terrestrial far-ultraviolet (FUV) dayglow is dominated by the 130.4-nm emission of neutral atomic oxygen, OI, and variations in the brightness observed from altitudes high above the emitting region reflect variations in thermospheric oxygen density. This paper summarizes the results of an initial survey of the Dynamics Explorer 1 observations of the FUV dayglow through a presentation of 13 representative events selected to demonstrate the spatial extent and short-term temporal stability of the brightness perturbations. The emphasis here is on the morning sector of local time and the polar cap for observations obtained in the time interval from September 23, 1981, through January 19, 1982. An analytic expression is derived for the average response of the FUV photometer to the dayglow during periods of high-latitude magnetic quiescence. The remaining observations in this time interval are then analyzed for their deviations from the established quiet time values. Deviations of -40% to +30% are found following intervals of increased magnetic activity. The most significant decreases (-30% to -40%) are observed equatorward of the instantaneous aurora] oval only after sustained periods (approx. 6 hours) of intense magnetic activity (average AE greater than approx. 700 nT). Decreases extend equatorward from the aurora to geographic latitudes as low as approx. 30 deg N. Decreases of lesser magnitude that do not extend as far equatorward are associated with sustained periods of more moderate activity in which the average value of AE is smaller (approx. 300-400 nT). Also, the spatial extent and magnitude of the decreases in the morning sector appear greater when the IMF B(sub y) component is positive. In both cases, decreases are readily observed within the polar cap. Localized enhancements of +20% to +30% occur much less frequently and are detected at the middle latitudes, well equator-ward of the auroral oval.

Nicholas, A. C.↗

Using AppletMagic(tm) to Implement an Orbit Propagator: New Life for Ada Objects

This paper will discuss the use of the Intermetrics AppletMagic tool to build an applet to display a satellite ground track on a world map. This applet is the result of a prototype project that was developed by the Goddard Space Flight Center's Flight Dynamics Division (FDD), starting in June of 1996. Both Version 1 and Version 2 of this applet can be accessed via the URL http://fdd.gsfc.nasa.gov/Java.html. This paper covers Version 1, as Version 2 did not make radical changes to the Ada part of the applet. This paper will briefly describe the design of the applet, discuss the issues that arose during development, and will conclude with lessons learned and future plans for the FDD's use of Ada and Java. The purpose of this paper is to show examples of a successful project using Oi AppletMagic, and to highlight some of the pitfalls that occurred along the way. It is hoped that this discussion will be useful both to users of AppletMagic and to organizations such as Intermetrics that develop new technology.

Stark, Michael E.↗

A 2 THz Heterodyne Array Receiver for SOFIA

We proposed to perform a comprehensive design study of a 16-element heterodyne array receiver for SOFIA. The array was designed to utilize hot-electron bolometers in an efficient, low-cost waveguide mount to achieve low noise performance between approx. 1500 and 2400 GHz. Due to the prevailing physical conditions in the interstellar medium, this frequency range is one of the richest in the FIR portion of the spectrum. An array designed for this wavelength range will make excellent use of the telescope and the available atmospheric transmission, and will provide a new perspective on stellar, chemical, and galaxy evolution in the present as well as past epochs. A few of the most important molecular and atomic species which the instrument will sample are CII, OI, CO, OH, NII, and CH. The system used the most sensitive detectors available in an efficient optical system. The local oscillator was a compact CO2 pumped far-infrared laser currently under development for SOFIA. The backend spectrometer was an array acousto-optic spectrometer (aAOS). The spectrometer utilizes proven hardware and technologies to provide broadband performance (greater than or equal to 1 GHz per AOS channel) and high spectral resolution (1 MHz) with the maximum sensitivity and minimum complexity and cost. The proposed instrument would be the fastest and most sensitive heterodyne receiver ever to operate in the 1.5 - 2.4 THz band. One of the key technologies developed for the proposed instrument is the laser micromachining of waveguide structures. These structures provide both the optical link between the instrument and the telescope (via an array of efficient feedhorns) and the impedance transformation between the detectors and free space. With the assistance of funds provided from this grant, we were able to fabricate and test the world's first laser micromachined feedhorns. Figure 1 is a photograph of the 2 THz double feedhorn structure designed and constructed under the auspices of this grant. The quality of the waveguide structure is far better than that obtainable using any other fabrication technique. Figure 2 is the beam pattern obtained from it. The beam parameters are an excellent match to what is expected from theory. The success of this experiment demonstrates the viability of using laser micromachined components in the development of high performance, large format array receivers. Figure 3 is an illustration of the focal plane array concept we developed lot- SOFIA base upon this work. Un- fortunately, our instrument proposal was not selected as a first generation SOFIA instrument. However, we have continued our development efforts and will propose to build a 2nd generation instrument based on the same design concepts. Our work under NAG 2-1057 was very rewarding and we appreciate the opportunity provided by NASA to pursue it. The research led directly to publishing 2 papers (listed below) and the award of an NSF grant to our group to construct a laser micromachining system on the University of Arizona campus.

Walker, Christopher K.↗

Stress Measurements on Blair High School Gymnasium: A Demonstration of Space Technology Transfer

This Report describes an actual demonstration of transfer to non-space use of technologies developed for space programs applications. Techniques used in assessing static and dynamic characteristics of the Blair High School gymnasium involved data acquisition by continuous scanning of strain gauge data acquired over a time of wide-temperature range, and analysis by a computer routine developed by Jet Propulsion Laboratory five years ago. The advantage of this method over conventional structural testing of uniquely designed structures was proved. More importantly, the process of demonstration was shown to be of great assistance to, and extension of, normal methods of disseminating information of new technologies. It is felt that significant benefit will derive from this improved mode oi concept transfer.

Kastel, Dean↗

Scattering Polarization in the Chromosphere

Scattering polarization from the photosphere observed close to the solar limb has recently become of interest to study turbulent magnetic fields, abundances, and radiative transfer effects. We extend these studies by measuring the scattering polarization off the limb, i.e. in the chromosphere. However, instrumental effects are much more pronounced and more complicated than those affecting on-disk measurements. In particular, scattered light from the telescope mirrors leads to a new type of instrumental polarization that we describe in detail. The differences between the linearly polarized spectra on the disk and off the limb are often very substantial. Here we show the profiles of HeI D(sub 3), the OI triplet at 777 nm, and the Nal D lines. The change in the latter is in reasonable agreement with the recent modeling efforts of atomic polarization in the lower level by Landi Degl'Innocenti (1998).

Keller, C. U.↗

A 2 THz Heterodyne Array Receiver for SOFIA: Summary of Research

We proposed to perform a comprehensive design study of a 16-element heterodyne array receiver for SOFIA. The array was designed to utilize hot-electron bolometers in an efficient, low-cost waveguide mount to achieve low noise performance between approximately 1500 and 2400 GHz. Due to the prevailing physical conditions in the interstellar medium, this frequency range is one of the richest in the Far-Infra Red (FIR) portion of the spectrum. An array designed for this wavelength range will make excellent use of the telescope and the available atmospheric transmission, and will provide a new perspective on stellar, chemical, and galactic evolution in the present as well as past epochs. A few of the most important molecular and atomic species which the instrument will sample are CII, OI, CO, OH, NII, and CH. The system used the most sensitive detectors available in an efficient optical system. The local oscillator was a compact CO2 pumped far-infrared laser currently under development for SOFIA. The backend spectrometer was an array acousto-optic spectrometer (aAOS). The spectrometer utilizes proven hardware and technologies to provide broadband performance (> 1 GHz per AOS channel) and high spectral resolution (1 MHz) with the maximum sensitivity and minimum complexity and cost. The proposed instrument would be the fastest and most sensitive heterodyne receiver ever to operate in the 1.5 - 2.4 THz band. One of the key technologies developed for the proposed instrument is the laser micromachining of waveguide structures. These structures provide both the optical link between the instrument and the telescope (via an array of efficient feedhorns) and the impedance transformation between the detectors and free space. With the assistance of funds provided from this grant, we were able to fabricate and test the world's first laser micromachined feedhorns. The quality of the waveguide structure is far better than that obtainable using any other fabrication technique. The beam parameters are an excellent match to what is expected from theory. The success of this experiment demonstrates the viability of using laser micromachined components in the development of high performance, large format array receivers. We have continued our development efforts and will propose to build a 2nd generation instrument based on the same design concepts.

Walker, Christopher K.↗

Mid-Infrared Observations of the Galactic Center

Under this grant, Susan Stolovy completed her thesis work and performed an analysis of the galactic center. For her thesis Susan analyzed observations of the galactic center obtained with the KAO using the KEGS spectrograph, built at Cornell. These observations present a study of the distribution and kinematics of the atomic gas in the inner few parsecs of the Galaxy as traced by the forbidden [SiII] line at 34.814 microns. The integrated [SiII] emission peaks near Sgr A* and extends past the inner edge of the Circumnuclear Disk (CND), passing through a gap in the dense molecular material to the northwest. The [SiII] maps have a spatial resolution of 15" and a spectral resolution of 50 km/s. The spectra, which are characterized by broad linewidths of order 100 km/s, are kinematically consistent with the CND rotation to the southwest but not to the north. The northern extension may be experiencing shocks and is likely to be infalling along the Northern Arm. Observations of high [Sill]/ [OI] and [SiII]/dust continuum ratios support the conjecture that turbulent motions and shocks in the inner few parsecs of the Galaxy are destroying dust grains, thus elevating the abundance of atomic silicon.

Stolovy, Susan↗

Studies for the Loss of Atomic and Molecular Species from Io

A summary discussion of research undertaken in this project is presented and is related to six published papers attached in the appendix. The discussion is divided into six sections describing a variety of studies for the loss of atomic and molecular species from Io. They include studies for: (1) atomic sodium, (2) SO2, (3) O and S, (4) spectacular high-spatial resolution ultraviolet image observations of O, S, and possibly H in Io's atmosphere and/or corona acquired by the Space Telescope Imaging Spectrometer (STIS) of the Hubble Space Telescope (HST), (5) spectacular high-spatial resolution visible Io eclipse image observations acquired by the Solid State Imager (SSI) of Galileo spacecraft, (6) ground-based observations acquired by the Solid State Imager (SSI) of Galileo spacecraft, and (7) ground-based observations of Io's neutral cloud in [OI] 6300 angstrom emission. New source rates at Io's exobase for SO2, O, and H are given and a variety of interesting implications for Io's atmosphere and for the Io plasma torus are discussed. Appendices that are comprised of articles published during the contract are also presented.

Smyth, William H.↗

A 2 THz Heterodyne Array Receiver for SOFIA

We proposed to perform a comprehensive design study of a 16-element heterodyne array receiver for SOFIA. The array was designed to utilize hot-electron bolometers in an efficient, low-cost waveguide mount to achieve low noise performance between approx. 1500 and 2400 GHz. Due to the prevailing physical conditions in the interstellar medium, this frequency range is one of the richest in the FIR portion of the spectrum. An array designed for this wavelength range will make excellent use of the telescope and the available atmospheric transmission, and will provide a new perspective on stellar, chemical, and galaxy evolution in the present as well as past epochs. A few of the most important molecular and atomic species which the instrument will sample are CII, OI, CO, OH, NII, and CH. The system used the most sensitive detectors available in an efficient optical system. The local oscillator was a compact CO2 pumped far-infrared laser currently under development for SOFIA. The backend spectrometer was an array acousto-optic spectrometer (aAOS). The spectrometer utilizes proven hardware and technologies to provide broadband performance (greater than or equal to 1 GHz per AOS channel) and high spectral resolution (1 MHz) with the maximum sensitivity and minimum complexity and cost. The proposed instrument would be the fastest and most sensitive heterodyne receiver ever to operate in the 1.5 - 2.4 THz band. One of the key technologies developed for the proposed instrument is the laser micromachining of waveguide structures. These structures provide both the optical link between the instrument and the telescope (via an array of efficient feedhorns) and the impedance transformation between the detectors and free space. With the assistance of funds provided from this grant, we were able to fabricate and test the world's first laser micromachined feedhorns. The quality of the waveguide structure is far better than that obtainable using any other fabrication technique. The beam parameters are an excellent match to what is expected from theory. The success of this experiment demonstrates the viability of using laser micromachined components in the development of high performance, large format array receivers. Unfortunately, our instrument proposal was not selected as a first generation SOFIA instrument. However, we have continued our development efforts and will propose to build a 2nd generation instrument based on the same design concepts.

Walker, Christopher K.↗

Processing and Properties of SiC/MoSi2-SiC Composites Fabricated by Melt Infiltration

Hi-Nicalon SiC fiber reinforced MoSi2-SiC matrix composites (SiC/MoSi2-SiC) have been fabricated by the melt infiltration approach. The composite consists of approximately 60 vol%, 2-D woven BN/SiC coated Hi-Nicalon SiC fibers and approximately 40 vol% MoSi2-SiC matrix. The room temperature tensile properties and thermal conductivity of the SiC/MoSi2-SiC composites were measured and compared with those of the melt infiltrated SiC/SiC composites. The influence oi fiber architecture on tensile properties was also evaluated. Results indicate that the primary modulus, stress corresponding to deviation from linearity, and transverse thermal conductivity values for the SiC/MoSi2-SiC composites are significantly lower than those for the SiC/SiC composites. Microcracking of the matrix due to the large difference in thermal expansion between MoSi2 and SiC appears to be the reason for the lower matrix dominated properties of SiC/MoSi2-SiC composites.

Bhatt, Ramakrishna T.↗

Near Earth Asteroid Rendezvous (NEAR) Revised Eros Orbit Phase Trajectory Design

Trajectory design of the orbit phase of the NEAR mission involves a new process that departs significantly from those procedures used in previous missions. In most cases, a precise spacecraft ephemeris is designed well in advance of arrival at the target body. For NEAR, the uncertainty in the dynamic environment around Eros does not allow the luxury of a precise spacecraft trajectory to be defined in advance. The principal cause of this uncertainty is the limited knowledge oi' the gravity field a,-id rotational state of Eros. As a result, the concept for the NEAR trajectory design is to define a number of rules for satisfying spacecraft, mission, and science constraints, and then apply these rules to various assumptions for the model of Eros. Nominal, high, and low Eros mass models are used for testing the trajectory design strategy and to bracket the ranges of parameter variations that are expected upon arrival at the asteroid. The final design is completed after arrival at Eros and determination of the actual gravity field and rotational state. As a result of the unplanned termination of the deep space rendezvous maneuver on December 20, 1998, the NEAR spacecraft passed within 3830 km of Eros on December 23, 1998. This flyby provided a brief glimpse of Eros, and allowed for a more accurate model of the rotational parameters and gravity field uncertainty. Furthermore, after the termination of the deep space rendezvous burn, contact with the spacecraft was lost and the NEAR spacecraft lost attitude control. During the subsequent gyrations of the spacecraft, hydrazine thruster firings were used to regain attitude control. This unplanned thruster activity used Much of the fuel margin allocated for the orbit phase. Consequently, minimizing fuel consumption is now even more important.

Helfrich, J↗

A Tool for Automatic Data Distribution for CFD Applications on Structured Grids

Development of HPF versions of NPB and ARC3D has shown that HPF provides an efficient, concise way to express parallelism and to organize data traffic. The use of HPF, as noted in the papers, requires an intimate knowledge of the applications and a detailed analysis of data affinity, data movement, and data granularity. To simplify and accelerate the task of developing HPF versions of existing CFD applications we have designed and implemented ADAPT (Automatic Data Alignment and Placement Tool). ADAPT analyzes a CFD application working on a single structured grid and generates HPF TEMPLATE, (RE)DISTRIBUTION, ALIGNMENT, and INDEPENDENT directives. The directives can be generated on the nest level, subroutine level, application level, or on the application interface level. ADAPT annotates an existing CFD FORTRAN application, performing computations on single or multiple grids. On each grid the application is considered as a sequence of operators, each applied to a set of variables defined in a particular grid domain. ADAPT automatically detects implicit operators (i.e., having data dependences) and explicit operators (without data dependences). For parallelization of an explicit operator ADAPT creates a template for the operator domain, aligns arrays used in the operator with the template, distributes the template, and declares the loops over the distributed dimensions as INDEPENDENT. For parallelization of an implicit operator, the distribution of the operator's domain should be consistent with the operator's dependences. Any dependence between sections distributed on different processors would preclude parallelization if the compiler does not have an ability to pipeline computations. If a data distribution is "orthogonal" to the dependences of an implicit operator, then the loop which implements the operator can be declared as INDEPENDENT. ADAPT starts with an analysis of array index expressions of the loop nests. For each pair of arrays referenced in an assignment statement, it generates an arc in the alignment graph and annotates it with an affinity relation. The template, alignment, and distribution directives for a particular loop nest are then derived from a transitive closure of the affinity relation. A compromise of data distributions in different nests and subroutines is achieved by merging annotated alignment graphs for adjacent nests/stibroutine calls in the nest/call graph of the application in the process called distribution lifting. ADAPT has been implemented as a C++ program running in conjunction with a parallelization tool called CAPTools. ADAPT uses the parse tree, interprocedural analysis and application database generated by CAPTools. It also uses the Directed Graph class, initially implemented in p2d2 (parallel debugger oi distributed programs), and some other classes supporting symbolic computations. ADAPT uses data distribution techniques described. ADAPT was tested with ARC3D and the FT benchmark and has demonstrated a code performance within a factor of 1.5 of handwritten versions.

Frumkin, Michael↗