ALBERT

Description: Wind tunnel tests have been conducted on an NACA 2412 airfoil section at Reynolds number of 2.2 x 10(exp 6) and Mach number of 0.13. Detailed measurements of flow fields associated with turbulent boundary layers have been obtained at angles of attack of 12.4 degrees, 14.4 degrees, and 16.4 degrees. Pre- and post-separated velocity and pressure survey results over the airfoil and in the associated wake are presented. Extensive force, pressure, tuft survey, hot-film survey, local skin friction, and boundary layer data are also included. Pressure distributions and separation point locations show good agreement with theory for the two layer angles of attack. Boundary layer displacement thickness, momentum thickness, and shape factor agree well with theory up to the point of separation. There is considerable disparity between extent of flow reversal in the wake as measured by pressure and hot-film probes. The difference is attributed to the intermittent nature of the flow reversal.

Description: We have studied the magnetic structure in AR 7150 (S09E06) observed on 29 April 1992 by the Soft X-Ray Telescope (SXT) on Yohkoh. The observed X-ray images are compared with force-free magnetic fields with different values of alpha, extrapolated from the Marshall Space Flight Center (MSFC) photospheric magnetogram observed at the same time. The results show that the magnetic field of the active region is not potential. Different groups of loops are characterized by different values of alpha. Bright loops correlation between the brightness of individual loops with the amount of twist. Further investigation of the magnetic state of the loop structure requires accurate nonlinear force-free calculations.

Description: After its fly-by of the planet Jupiter in February 1992, the Ulysses spacecraft is now in a highly inclined heliocentric orbit that will bring it above the south polar regions of the Sun in September 1994. The high-latitude phenomena observed to date have been strongly influenced by the near-minimum solar activity conditions encountered during this phase of the mission. In late April 1993, when Ulysses was at approximately 29 deg S heliographic latitude, the recurrent high speed solar wind stream that had been observed at the location of the spacecraft for 11 consecutive solar rotation underwent a dramatic change. The wind speed in the valleys between successive peaks increased in a single step from approximately 420 km/s to aopproximately 560 km/s. This change in solar wind flow was accompanied by the disappearance at the spacecraft of the magnetic sector structure that had been observed until then. Both these finding are consistent with Ulysses having climbed beyond the latitude of the coronal streamer belt in which is embedded the heliospheric current sheet (HCS). In its subsequent poleward journey, no further evidence for an encounter with the HCS has been seen at Ulysses. Other phenomena observed include the evolution with latitude of corotating interaction region (CIRs) and their influence on the acceleration of energetic particles, and the characteristics of the solar wind flows emanating from the south polar coronal hole. In this paper, we present details of the above observations. Finally, while the polar passes of the prime mission will take place near solar minimum, an extended mission will bring Ulysses back over the poles near the maximum of the next cycle. A summary of scientific goals for Ulysses at solar maximum is given.

Description: Extensive hard X-ray (HXR)/gamma-ray (GR) observations of solar flares, performed during solar cycles 21 and 22 have led to important new discoveries. These data, combined with observations obtained in other parts of the electromagnetic spectrum (soft X-ray, Hard X-ray, optical, and radio) largley contributed to get a better understanding and to develop new ideas on particle acceleration and transport during solar flares. This review presents new observational facts relevant to hard X-ray/gamma-ray producing flares. Among these are the frequent presence of sub-second time structure in the hard X-ray emission, the variability in hard X-ray and radio spatial distributions during a flare and from flare to flare, the evidence for strong gamma-ray line emission from the Corona and the existence of extended phases of the gamma-ray emission lasting for several hours after the flare onset. This ensemble of observations indicates that particle acceleration takes place at different sites in a complex and dynamic magnetic field environment.

Description: A flood of new observations of the solar corona have been made with high spatial resolution, good temporal coverage and resolution, and large linear dynamic range by the Soft X-ray Telescope (SXT) on Yohkoh. These data are changing our fundamental understanding of how solar magnetic fields emerge, interact, and dissipate. This paper reviews some of the results from Yohkoh in the context of earlier results from the Solar Maximum Mission (SMM) and in comjunction with ground-based optical and radio observations.

Description: We have studied the relation between flux emergence and flare activity in the active region NOAA 7260, using images from the Soft X-ray Telescope (SXT) aboard the Yohkoh spacecraft and other supporting ground-based data. It is found that microflares start around the time of flux emergence as recorded in white-light data, which generally precedes a major flare by several hours. We interpret the microflares as due to fast reconnection that takes place intermittently in the slow reconnection stage while more energy is accumulated in preparation for a larger flare.

Description: The implementation of a two-equation k-omega turbulence model into the NPARC flow solver is described. Motivation for the selection of this model is given, major code modifications are outlined, new imputs to the code are described, and results are presented for several validation cases: an incompressible flow over a smooth flat plate, a subsonic diffuser flow, and a shock-induced separated flow. Comparison of results with the k-epsilon model indicate that the k-omega model predicts simple flows equally well whereas, for adverse pressure gradient flows, the k-omega model outperforms the other turbulence models in NPARC.

Description: The NPARC Alliance is a partnership between the NASA Lewis Research Center (LeRC) and the USAF Arnold Engineering Development Center (AEDC) dedicated to the establishment of a national CFD capability, centered on the NPARC Navier-Stokes computer program. The three main tasks of the Alliance are user support, code development, and validation. The present paper is a status report on the validation effort. It describes the validation approach being taken by the Alliance. Representative results are presented for laminar and turbulent flat plate boundary layers, a supersonic axisymmetric jet, and a glancing shock/turbulent boundary layer interaction. Cases scheduled to be run in the future are also listed. The archive of validation cases is described, including information on how to access it via the Internet.

Description: An approach for solving the compressible Euler and Navier-Stokes equations upon meshes composed of nearly arbitrary polyhedra is described. Each polyhedron is constructed from an arbitrary number of triangular and quadrilateral face elements, allowing the unified treatment of tetrahedral, prismatic, pyramidal, and hexahedral cells, as well the general cut cells produced by Cartesian mesh approaches. The basics behind the numerical approach and the resulting data structures are described. The accuracy of the mixed volume grid approach is assessed by performing a grid refinement study upon a series of hexahedral, tetrahedral, prismatic, and Cartesian meshes for an analytic inviscid problem. A series of laminar validation cases are made, comparing the results upon differing grid topologies to each other, to theory, and experimental data. A computation upon a prismatic/tetrahedral mesh is made simulating the laminar flow over a wall/cylinder combination.

Description: The paper lists US solar missions, both those planned and funded by NASA alone as well as those carried out in collaboration with other space agencies. Soe of the missions are now in operation, the others are either planned and approved or under active discussion. The paper also describes the principal scientific objects of the missions and gives some orbital characteristics.

Description: Interplanetary scintillation (IPS) measurements of the 'disturbance factor' g, obtained with the Cambridge (UK) array can be used to explore the heliospheric density structure. We have used these data to construct synoptic (Carrington) maps, representing the large-scale enhancements of the g-factor in the inner heliosphere. These maps emphasize the stable corotating, rather than the transient heliospheric density enhancements. We have compared these maps with Carrington maps of Fe XIV observations National Solar Observatory ((NSO), Sacramento Peak) and maps based on Yohkoh Soft X-Ray Telescope (SXT) X-ray observations. Our results indicate that the regions of enhanced g tend to map to active regions rather than the current sheet. The implication is that act ve regions are the dominant source of the small-scale (approximately equal 200 km) density variations present in the quiet solar wind.

Description: The large-scale structure of the solar corona is investigated using synoptic maps produced from Fe XIV (530.3 nm), Fe X (637.4 nm) and Ca XV (569.4 nm) data obtained at National Solar Observatory (NSO/SP), Yohkoh/Soft X-ray Telescope (SXT) X-ray data and Wilcox Solar Observatory (WSO) 'source surface' maps. We find that the Fe XIV data are an excellent proxy for spatially-average Yohkoh/SXT data. Isolated emission features and large-scale structures are nearly identical in SXT and Fe XIV maps. In addition, coronal holes and other low-emission regions are very similar. Synoptic temperature maps, calculated from the Fe X/Fe XIV ratio, show a tendency for the highest temperatures to occur where the large-scale magnetic fields change polarity at high latitudes, while lower-latitude features, including active regions, have lower apparent tempertures. Regions of enhanced temperature generally follow the helisopheric current sheet (HCS) as defined by the WSO maps. Further, emission in Ca XV (formed at T is approximately equal to 3 MK), generally occurs only over low-latitude regions that are bright in both FE X (T approximately equal to 1 MK) and Fe XIV (T approximately equal to 2 MK). Thus, there is evidence for low (approximately equal to 1 MK), moderate (approximately equal to 2 MK) and high (approximately 3 MK) temperatures in close proximity in the low corona.

Description: Coronal mass ejections (CME's) are thought to result from the loss of stability within a magnetically confined coronal structure leading to its radial expansion into interplanetary space. As the CME expands into the corona current sheets will form between the expanding CME and surrounding field lines in the ambient wind. This configuration may lead to reconnection between the CME and adjacemt field lines. Such reconnection may produce double ion beams as has been observed in the terrrestrial magnetosphere. We examine all 24 distinct signatures of CME's observed by Ulysses during the in-ecliptic portion of the mission. In 5 of these 24 cases the ion spectra were not clear and thus the presence of double ion beams could not be determined. In 13 of the remaining 19 CME's double ion beams were found on the leading and/or trailing edge of the CME but not in the interior of the CME. In 3 of the CME's double ion beams were found throughout the CME while in the remaining 3 CME's double ion beams were not present near or just inside of the CME. In contrast in a control sample of 19 randomly chosen intervals, double ion beams were present at the leading and/or trailing edges of the random intervals in only 3 of the 19 cases. There appears to be no correlation between probability of occurrence of double ion beams and a magnetic cloud or non-cloud configuration of the CME and no correlation between the presence of the double ion beams at the edges of the CME and the CME being a fast or slow CME.

Description: Imaging observations of solar flare hard X-ray sources with the Hard X-ray Telescope (HXT) aboard the Yohkoh satellite have revealed that hard X-ray emissions (greater than 30 ke V) originate most frequently from double sources. The double sources are located on both sides of the magnetic neutral line, suggesting that the bulk of hard X-rays is emitted from footpoints of flaring magnetic loops. We also found that hard X-rays from the double sources are emitted simultaneously within a fraction of second and that the weaker source tends to be located in the stronger magnetic field region, showing a softer spectrum. Physcial implications on the observed characteristics of the hard X-ray double sources are discussed.

Description: The scaling properties of a time series of Doppler images obtained in good visibility conditions are studied. A 28 cm vacuum telescope and a vacuum spectroheliograph in video spectra-spectroheliograph mode, are used. Sixty line-of-sight Doppler images of an area of the quiet sun are investigated. They were taken at 60 sec intervals over a one hour span and have a 2 arcsec resolution. After the removal of the five-minute oscillations, the time-spatial spectrum is calculated. To study the turbulence of photospheric flows, two scaling parameters in the spectra, are estimated: the exponent of the spatial part of the power spectrum, and the exponent governing the scaling of time correlations. The implied diffusive behavior is discussed. This includes the estimation of a diffusion coefficient and the type of diffusion involved.

Description: The nonlocal non-diffusive transport of passive scalars in turbulent magnetohydrodynamic (MHD) convection is investigated using transilient matrices. These matrices describe the probability that a tracer particle beginning at one position in a flow will be advected to another position after some time. A method for the calculation of these matrices from simulation data which involves following the trajectories of passive tracer particles and calculating their transport statistics, is presented. The method is applied to study the transport in several simulations of turbulent, rotating, three dimensional compressible, penetrative MDH convection. Transport coefficients and other diagnostics are used to quantify the transport, which is found to resemble advection more closely than diffusion. Some of the results are found to have direct relevance to other physical problems, such as the light element depletion in sun-type stars. The large kurtosis found for downward moving particles at the base of the convection zone implies several extreme events.

Description: Long uninterrupted sequences of solar magnetograms from the global oscillations network group (GONG) network and from the solar and heliospheric observatory (SOHO) satellite will provide the opportunity to study the proper motions of magnetic features. The possible use of multiscale regularization, a scale-recursive estimation technique which begins with a prior model of how state variables and their statistical properties propagate over scale. Short magnetogram sequences are analyzed with the multiscale regularization algorithm as applied to optical flow. This algorithm is found to be efficient, provides results for all the spatial scales spanned by the data and provides error estimates for the solutions. It is found that the algorithm is less sensitive to evolutionary changes than correlation tracking.

Description: Ulysses has collected data between 1 and 5 AU during, and just following solar maximum, when the heliospheric current sheet (HCS) can be thought of as reaching its maximum tilt and being subject to the maximum amount of turbulence in the solar wind. The Ulysses solar wind plasma instrument measures the vector velocity and can be used to estimate the flow speed and direction in turbulent 'eddies' in the solar wind that are a fraction of an astronomical unit in size and last (have either a turnover or dynamical interaction time of) several hours to more than a day. Here, in a simple exercise, these solar wind eddies at the HCS are characterized using Ulysses data. This character is then used to define a model flow field with eddies that is imposed on an ideal HCS to estimate how the HCS will be deformed by the flow. This model inherently results in the complexity of the HCS increasing with heliocentric distance, but the result is a measure of the degree to which the observed change in complexity is a measure of the importance of solar wind flows in deforming the HCS. By comparison with randomly selected intervals not located on the HCS, it appears that eddies on the HCS are similar to those elsewhere at this time during the solar cycle, as is the resultant deformation of the interplanetary magnetic field (IMF). The IMF deformation is analogous to what is often termed the 'random walk' of interplanetary magnetic field lines.

Description: We study the onset of a pure Marangoni convection in a liquid layer with two deformable interfaces in the no-gravity environment. Both oscillatory and stationary instabilities are considered for a wide range of parameters. It is shown that only stationary instability is possible when surface tension at the colder interface is lower than that at the hotter one. Oscillatory instability tends to disappear and to be replaced by the stationary instability with increase of the Prandtl number and decrease of surface tension at the colder interface.

Description: Several physical and observational effects contribute to the significant imbalances of magnetic flux that are often observed in active regions. We consider an effect not previously treated: the influence of electric currents in the photosphere. Electric currents can cause a line-of-sight flux imbalance because of the directionality of the magnetic field they produce. Currents associated with magnetic flux tubes produce larger imbalances than do smoothly-varying distributions of flux and current. We estimate the magnitude of this effect for current densities, total currents, and magnetic geometry consistent with observations. The expected imbalances lie approximately in the range 0-15%, depending on the character of the current-carying fields and the angle from which they are viewed. Observationally, current-induced flux imbalances could be indicated by a statistical dependence of the imbalance on angular distance from disk center. A general study of magnetic flux balance in active regions is needed to determine the relative importance of other- probably larger- effects such as dilute flux (too weak to measure or rendered invisible by radiative transfer effects), merging with weak background fields, and long-range connections between active regions.

Description: Measurements of wing buffeting, using root strain gages, were made in the NASA Langley 0.3 m cryogenic wind tunnel to refine techniques which will be used in larger cryogenic facilities such as the United States National Transonic Facility (NTF) and the European Transonic Wind Tunnel (ETW). The questions addressed included the relative importance variations in frequency parameter and Reynolds number, the choice of model material (considering both stiffness and damping) and the effects of static aeroelastic distortion. The main series of tests was made on three half models of slender 65 deg delta wings with a sharp leading edge. The three delta wings had the same planform but widely differing bending stiffnesses and frequencies (obtained by varying both the material and the thickness of the wings). It was known that the steady flow on this configuration would be insensitive to variations in Reynolds number. On this wing at vortex breakdown the spectrum of the unsteady excitation is unusual, having a sharp peak at particular frequency parameter. Additional tests were made on one unswept half-wing of aspect ratio 1.5 with an NPL 9510 aerofoil section, known to be sensitive to variations in Reynolds number at transonic speeds. The test Mach numbers were M = 0.21 and 0.35 for the delta wings and to M = 0.30 for the unswept wing. On this wing the unsteady excitation spectrum is fairly flat (as on most wings). Hence correct representation of the frequency parameter is not particularly important.

Description: A consistent solution of the radiative transfer equation characterizing photon transport in a semi-infinite medium of refractive index greater than or equal to one is obtained following the method of Sobolev. Fresnel specular reflection, Snell's law and isotropic scattering are assumed. An algorithm is developed and its accuracy is demonstrated. A numerical Laplace transform inversion leads to an efficient evaluation for the interior flux and source function distributions.

Description: A semianalytic method is derived for dealing simultaneously with large numbers of linear stellar oscillation modes trapped in a cavity (a shell) of fluid which is rotating and convecting. A simple generalization of mixing-length theory shows how convection is modulated by weak rotational effects and by the horizontal wind fields of linear r-mode oscillations. The modulated convection is then used to compute the energy lost to turbulent viscosity by a family of nondegenerate oscillations. Viscosity terms of fourth degree in the wind shear can be included if they are a perturbation affecting only a small portion of the r-mode. Viscous energy loss strenghthens convection in a narrow layer near the base of the H and He ionization zone. In the Sun, this layer is about 7 Mm thick and centered at 0.932 of a solar radius where convection cells have a typical size of about 20 Mm and a lifetime of 0.3 Ms, both similar to what is observed in supergranules. If the rms velocity of r-modes at the surface exceeds 5 m/s, then energy is deposited inside the Sun at a sufficient rate to power the supergranulation and impose on it a weak latitude dependence.

Description: We present observational evidence that eruptions of quiescent filaments and associated coronal mass ejections (CMEs) occur as a consequence of the destabilization of large-scale coronal arcades due to interactions between these structures and new and growing active regions. Both statistical and case studies have been carried out. In a case study of a 'bulge' observed by the High-Altitude Observatory Solar Maximum Mission coronagraph, the high-resolution magnetograms from the Big Bear Solar Observatory show newly emerging and rapidly changing flux in the magnetic fields that apparently underlie the bugle. For other case studies and in the statistical work the eruption of major quiescent filaments was taken as a proxy for CME eruption. We have found that two thirds of the quiescent-filament-associated CMEs occurred after substantial amounts of new magnetic flux emerged in the vicinity of the filament. In addition, in a study of all major quiescent filaments and active regions appearing in a 2-month period we found that 17 of the 22 filaments that were associated with new active regions erupted and 26 of the 31 filaments that were not associated with new flux did not erupt. In all cases in which the new flux was oriented favorably for reconnection with the preexisting large-scale coronal arcades; the filament was observed to erupt. The appearance of the new flux in the form of new active regions begins a few days before the eruption and typically is still occurring at the time of the eruption. A CME initiation scenario taking account of these observational results is proposed.

Description: Total solar irradiance measurements from the 1984-1993 Earth Radiation Budget Satellite (ERBS) active cavity radiometer and 1978-1993 Nimbus 7 transfer cavity radiometer spacecraft experiments are analyzed to detect the presence of 11-, 22-, and 80-year irradiance variability components. The analyses confirmed the existence of a significant 11-year irradiance variability component, associated with solar magnetic activity and the sunspot cycle. The analyses also suggest the presence of a 22- or 80-year variability component. The earlier Nimbus 7 and Solar Maximum Mission (SMM) spacecraft irradiance measurements decreased approximately 1.2 and 1.3 W/sq m, respectively, between 1980 and 1986. The Nimbus 7 values increased 1.2 W/sq m between 1986 and 1989. The ERBS irradiance measurements increased 1.3 W/sq m during 1986-1989, and then decreased 0.4 W/sq m (at an annual rate of 0.14 W/sq. m/yr) during 1990-1993. Considering the correlations between ERBS, Nimbus 7, and SMM irradiance trends and solar magnetic activity, the total solar irradiance should decrease to minimum levels by 1997 as solar activity decreases to minimum levels, and then increase to maximum levels by the year 2000 as solar activity rises. The ERBS measurements yielded 165.4 +/- 0.7 W/sq m as the mean irradiance value with measurement accuracies and precisions of 0.2% and 0.02%, respectively. The ERBS mean irradiance value is within 0.2% of the 1367.4, 1365.9, and 1366.9 W/sq m mean values for the SMM, Upper Atmosphere Research Satellite (UARS), and Space Shuttle Atmospheric Laboratory for Applications and Science (ATLAS 1) Solar Constant (SOLCON) active cavity radiometer spacecraft experiments, respectively. The Nimbus 7 measurements yielded 1372.1 W/sq m as the mean value with a measurement accuracy of 0.5%. Empirical irradiance model fits, based upon 10.7 -cm solar radio flux (F10) and photometric sunspot index (PSI), were used to assess the quality of the ERBS, Numbus 7, SMM, and the UARS irradiance data sets and to identify irradiance variability trends which may be caused by drifts or shifts in the spacecraft sensor responses. Comparisons among the fits and measured irradiances indicate that the Nimbus 7 radiometer response shifted by a total of 0.8 W/sq m between September 1989 and April 1990 and that the ERBS and UARS radiometers each drifted approximately 0.5 W/sq m during the first 5 months in orbit.

Description: Results are reported of the Surface Tension Driven Convection Experiment (STDCE) aboard USML-1 Spacelab. Steady and transient thermocapillary flows were investigated in a 10 cm dia. circular container filled with 10 Cs silicone oil. The velocity and temperature fields were studied in detail under various conditions. It is shown in this paper how the Marangoni number affects the velocity field. A numerical analysis of the flows was also conducted and its results were compared to the experimental data. Good agreement is shown.

Description: Compressible MHD simulations in one dimension with three-dimensional vectors are used to investigate a number of processes relevant to problems in interplanetary physics. The simulations indicate that a large-amplitude nonequilibrium (e.g., linearly polarized) Alfvenic wave, which always starts with small relative fluctuations in the magnitude B of the magnetic field, typically evolves to flatten the magnetic profile in most regions. Under a wide variety of conditions B and the density rho become anticorrelated on average. If the mean magnetic field is allowed to decrease in time, the point where the transverse magnetic fluctuation amplitude delta B(sub T) is greater than the mean field B(sub 0) is not special, and large values of delta B(sub T)/B(sub 0) do not cause the compressive thermal energy to increase remarkably or the wave energy to dissipate at an unusually high rate. Nor does the 'backscatter' of the waves that occurs when the sound speed is less than the Alfven speed result, in itself, in substantial energy dissipation, but rather primarily in a phase change between the magnetic and velocity fields. For isolated wave packets the backscatter does not occur for any of the parameters examined; an initial radiation of acoustic waves away from the packet establishes a stable traveling structure. Thus these simulations, although greatly idealized compared to reality, suggest a picture in which the interplanetary fluctuations should have small deltaB and increasingly quasi-pressure balanced compressive fluctuations, as observed, and in which the dissipation and 'saturation' at delta B(sub T)/B(sub 0) approximately = 1 required by some theories of wave acceleration of the solar wind do not occur. The simulations also provide simple ways to understand the processes of nonlinear steepening and backscattering of Alfven waves and demonstrate the existence of previously unreported types of quasi-steady MHD states.

Description: Frequency shifts of high frequency p-modes during the solar cycle are calculated for a non-magnetic polytrope convection zone model. An isothermal chromospheric atmosphere threaded by a uniform horizontal magnetic field is correlated to this model. The relevant observations of such frequency changes are discussed. The calculated simultaneous changes in the field strength and chromospheric temperature result in the frequency shifts that are similar to those of the observations.

Description: The specific attraction and, in large part, the significance of solar magnetograms lie in the fact that they give the most important data on the electric currents and the nonpotentiality of active regions. Using the vector magnetograms from the Marshall Space Flight Center (MSFC), we employ a unique technique in the area of data analysis for resolving the 180 deg ambiguity in order to calculate the spatial structure of the vertical electric current density. The 180 deg ambiguity is resolved by applying concepts from the nonlinear multivariable optimization theory. The technique is shown to be of particular importance in very nonpotential active regions. The characterization of the vertical electric current density for a set of vector magnetograms using this method then gives the spatial scale, locations, and magnitude of these current systems. The method, which employs an intermediate parametric function which covers the magnetogram and which defines the local `preferred' direction, minimizes a specific functional of the observed transverse magnetic field. The specific functional that is successful is the integral of the square of the vertical current density. We find that the vertical electric current densities have common characteristics for the extended bipolar (beta) (gamma) (delta)-regions studied. The largest current systems have j(sub z)'s which maximizes around 30 mA/sq m and have a linear decreasing distribution to a diameter of 30 Mn.

Description: Relative abundances of oxygen, neon, and magnesium have been derived for a sample of nine solar active regions, flares, and an erupting prominance by combining plots of the ion differential emission measures. The observations were photographed in the 300-600 A range by the Naval Research Laboratory (NRL) spectroheliograph on Skylab. Methods for deriving the Mg/Ne abundance ratio-which measures the separation between the low- first ionization potential (FIP) and high-FIP abundnace plateaus-have been described in previous papers. In this paper we describe the spectroscopic methods for deriving the O/Ne abundance ratio, which gives the ratio between two high-FIP elements. The plot of the O/Ne ratio versus the Mg/Ne ratio in the sample of nine Skylab events is shown. The variation in the Mg/Ne ratio by a factor of 6 is associated with a much smaller range in the O/Ne ratio. This is broadly consistent with the presence of the standard FIP pattern of abundances in the outer atmosphere of the Sun. However, a real change in the relative abundances of oxygen and neon by a factor of 1.5 cannot be excluded.

Description: We use the ray description of acoustic-gravity modes to calculate time-distance diagrams for the quiet Sun and for regions in the vicinity of a sunspot with a monolithic flux-tube structure. Time-distance curves for the quiet Sun match the observations of Duvall et al. In the vicinity of a sunspot these quiet Sun curves split into a family of closely spaced curves. The structure of this bandlike feature is found to be sensitive to the sunspot model and can be a diagnostic of the subsurface geometry of the sunspot flux tube.

Description: The signature of the solar cycle appears in helioseismic frequencies and splittings. It is known that the changing outer superadiabatic region of the sun is responsible for this. The deeper solar-cycle mechanism from the surface changes, and, in particular, how magnetic fields perturb the global modes, the solar irradiance and the luminosity, is discussed. The irradiance and helioseismic changes are described. The interpretation of seismic and photometric data is discussed, considering current one-dimensional models and phenomenology. It is discussed how the long term solar-cycle luminosity effect could be caused by changes occurring near the base of the convection zone (CZ). It is shown that a thin toroidal flux sheath at the top of the radiative zone changed the thermal stratification immediately below the CZ over a solar-cycle timescale in two ways: the temperature of the magnetized fluid becomes hotter than the surrounding fluid, and the temperature gradient steepens above the magnetized region. The testing of CZ dynamics and extension of numerical experiments to global scales are considered.

Description: The chromospheric bright points are the sites where intense heating occurs in three minute period waves. The bright points are grouped into three classes depending on the amount of intensity enhancement and the pattern of their dynamical evolution. A 35-minute time series of photographic spectra in the Ca(II) H line on a quiet region ofthe center of the solar disk was used to show that the period of intensity oscillations seen at sites of the bright points is independent of their intensity enhancements. The series was also used to show that the period may not depend on the strength of the magnetic fields with which they are associated. A linear regression equation was fitted to a curve representing the variation of the period of intensity oscillations with the peak value of I(sub H2V). The correlation coefficient was found to be 0.19.

Description: The preliminary results of the photometry of CaII K spectroheliograms are presented. From the spectrograms for 1992, plages, the magnetic network, intranetwork elements and the chromospheric background were separated using the histogram method. The intensity and area of these separated features, as well as the full disk intensity, were derived. The spatial K index was compared to the spectral CaII K index derived from line profiles. It was found that the spatial K index and intensity of plages, the network elements and the intranetwork and background regions were highly correlated with the MgII h and k c/w ratio.

Description: The observations made in July 1994 on the impact of fragment A of the comet P/Shoemaker-Levy 9 with Jupiter are described. The instrumentation used was a magneto-optical filter, acting as a two-channel filter. The data showed a double-peak transient which occurred after the impact, and whose general properties indicated a true jovian origin. The peaks appear in absorption. A numerical simulation can explain the main characteristics of the observed signal where the two peaks have the same polarity and appear only in the channel at shorter wavelengths. The simulation carried out appeared to indicate that the observed signal could be produced by the combination of shock waves and the expanding material with a velocity of 13 +/- 8 km/s. This implies that two separate impacts may have been observed. The developed simulation can be extended to predict long term effects.

Description: Both weak magnetic fields and latitudinally dependent acoustic perturbations remove the degeneracy of the azimuthal quantum number, m, of acoustic modes of otherwise spherically symmetrical solar model. In the case of acoustic perturbations, the degeneracy is removed because the range of latitude in which a mode propagates depends on m, and therefore modes of like principle order n and degree l sample the aspherical scalar sound speed distribution differently. In the magnetic case, the removal of the degeneracy is caused by the same geometrical effect, and is influenced by the anisotropy of the Lorentz forces. Asymptotic analysis is used to show that the frequency splittings cannot be unambiguously attributed to the direct effect of a magnetic field, and that the effect of such a field on frequency splittings can be reproduced by a perturbation to the sound speed.

Description: The central intensities of Na(I) D1 and D2 linear profiles at the sites of the chromospheric bright points in the interior of the supergranulation cells were derived from photographic spectra. The observation scheme sampled spectra simultaneously in seven lines at a repetition rate of 12 sec. It is shown that the Na(I) D1 and D2 lines exhibit a four minute periodicity in their intensity oscillations. It is seen that the period of intensity oscillations decreases outwardly from the photosphere to the corona. It is surmised that the spatial and temporal relationships between intensity and/or velocity in the photosphere and chromosphere may explain the physical mechanisms of the underlying oscillations.

Description: The dependence of the brightness of chromospheric network elements on latitude was investigated for quiet solar regions. Calibrated photographic CaII K-spectroheliograms were used to compare the variation in brightness at the center of the disk with higher latitude of chromospheric network elements in a quiet region as a function of solar activity. It was found that there was no significant difference in brightness between the center of the solar disk and higher latitude. It is concluded that the brightness of the chromospheric network elements in a quiet region does not depend on the latitude, but that the variation in the intensity enhancement is related to the level of solar activity.

Description: The preliminary analysis of a 69 day observation run taken at the JPL using the magneto-optical filter is presented. The aim is to estimate the rotational splitting of l = 1 modes. A value of Delta nu = 0.44 +/- 0.09 micro-Hz is found. In a second, more accurate analysis, it is planned to investigate the low frequency part of the power spectrum. The observational statistics are presented.

Description: The plane-wave decomposition of the acoustic-gravity wave effects observed in the photosphere provides a computationally efficient technique that probes the structure of the upper convective zone and boundary. In this region, the flat sun approximation is considered as being reasonably accurate. A technique to be used for the systematic plane-wave analysis of Michelson Doppler imager data, as part of the solar oscillations investigation, is described. Estimates of sensitivity are presented, and the effects of using different planar mappings are discussed. The technique is compared with previous approaches to the three dimensional plane-wave problem.

Description: It was observed that the p-mode power is substantially suppressed in magnetic regions. One possible explanation is that the upper turning point, the acoustic cut-off point of the solar p-modes is lowered in the presence of a magnetic field. A related possibility is that the attenuation length scale in the evanescent region is reduced in the presence of a magnetic field. It is likely that the observations sample a different position in the evanescent tails of the eigenfunctions in magnetic regions because of different temperature structures in these regions. A model is used to quantify the first of these effects.

Description: The fluctuations in magnetic field and plasma velocity in solar wind, which possess many features of fully developed magnetohydrodynamic (MHD) turbulence, are discussed. Direct spacecraft observations from 0.3 to over 20 AU, remote sensing radio scintillation observations, numerical simulations, and various models provide complementary methods that show that the fluctuations in the wind parameters undergo significant dynamical evolution independent of whatever turbulence might exist in the solar photosphere and corona. The Cluster mission, with high time resolution particle and field measurements and its variable separation strategies, should be able to provide data for answering many questions on MHD turbulence.

Description: Double ion beams are often observed in the solar wind, but little work has been done in relating these beams to structures within the solar wind. Double ion beams are observed as beams of a given ion species and charge state occurring at two different energies. We use the three-dimensional ion plasma instrument on board the Ulysses spacecraft to look for evidence of such beams associated with the heliospheric current sheet. In a subset chosen independently of plasma parameters consisting of 8 of cover 47 crossings of the current sheet made during the inecliptic phase of the Ulysses mission we find that these double ion beams are always present on either side of the current sheet. The double beams are present in both the proton and helium species. The secondary beam typically has a higher helium abundance, which suggests that these beams are formed in the helium-rich corona rather than in interplanetary space. The double beams are not present in the interior of the current sheet. Neither collisions nor effects of plasma beta can account for the disappearance of the double beams inside the current sheet in all eight cases. We postulate that these beams are formed by reconnection occurring near the Sun in the boundary region between the open field lines of the coronal holes and the closed field line region of the heliospheric current sheet. Such a scenario would be consistent with previous X ray measurements which suggect that reconnection is occurring in this region.

Description: A global resistive, two-dimensional, time-dependent magnetohydrodynamic (MHD) model is used to introduce and support the hypothesis that the quiet solar middle chromosphere is heated by resistive dissipation of large-scale electric currents which fill most of its volume. The scale height and maximum magnitude of the current density are 400 km and 31.3 m/sq m, respectively. The associated magnetic field is almost horizontal, has the same scale height as the current density, and has a maximum magnitude of 153 G. The current is carried by electrons flowing across magnetic field lines at 1 m/s. The resistivity is the electron contribution to the Pedersen resitivity for a weakly ionized, strongly magnetized, hydrogen gas. The model does not include a driving mechanism. Most of the physical quantities in the model decrease exponentially with time on a resistive timescale of 41.3 minutes. However, the initial values and spatial; dependence of these quantities are expected to be essentially the same as they would be if the correct driving mechanism were included in a more general model. The heating rate per unit mass is found to be 4.5 x 10(exp 9) ergs/g/s, independent of height and latitude. The electron density scale height is found to be 800 km. The model predicts that 90% of the thermal energy required to heat the middle chromosphere is deposited in the height range 300-760 km above the temperature minimum. It is shown to be consistent to assume that the radiation rate per unit volume is proportional to the magnetic energy density, and then it follows that the heating rate per unit volume is also proportional to the energy from the photosphere into the overlying chromosphere are briefly discussed as possible driving mechanisms for establishing and maintaining the current system. The case in which part of or all of the current is carried by protons and metal ions, and the contribution of electron-proton scattering to the current are also considered, with the conclusion that these effects do not change the qualitative prediction of the model, but probably change the quantitative predictions slightly, mainly by increasing the maximum magntiude of the current density and magnetic field to at most approximately 100 mA/m and approximately 484 G, respectively. The heating rate per unit mass, current density scale height, magnetic field scale height, temperatures, and pressures are unchanged or are only slightly changed by including these additional effects due to protons and ions.

Description: A number of different solar constant observations all made from space during the ATLAS 2 mission have been gathered and compared to each other. The Sun did not have a single sunspot during several days. As eight of the radiometric channels were all within 0.1%, the mean of the observations has been used to determine a set of adjustment factors providing de facto the definition of the Space Absolute Radiometric Reference (SARR). The differential absolute radiometers of Solar Constant (SOLCON) experiment and the Solar Variability-1 (SOVA 1) experiment, as well as the SOVA 2 and Active Cavity Radiometer (ACR) radiometers that have been brought back to the Earth may, if used in the same conditions, reproduce and maintain the SARR for the future.

Description: Eleven tests were carried out in DLR's high enthalpy tunnel in Goettingen (HEG), Germany, for reservoir conditions ranging from 10 to 23 MJ/kg and a free stream Mach number of approximately 10. A blunted cone model with a cylindrical afterbody (sting) was investigated. To obtain information on the influence of defined parameters on the body back flow, especially of the reacting gas, the heat transfer rate along the body contour was measured with fast response surface thermocouples on the forebody and sensitive thin film heat transfer gauges on the base and sting of the model. Flow visualization with a holographic interferometry system was provided. Tests are described, and a preliminary interpretation of the observed flow effects are given.

Description: Get Away Special (GAS) payload G-093, also called VORTEX (VOrtex Ring Transit EXperiment), is an investigation of the propagation of a vortex ring through a liquid-gas interface in microgravity. This process results in the formation of one or more liquid droplets similar to earth based liquid atomization systems. In the absence of gravity, surface tension effects dominate the drop formation process. The Shuttle's microgravity environment allows the study of the same fluid atomization processes as using a larger drop size than is possible on Earth. This enables detailed experimental studies of the complex flow processes encountered in liquid atomization systems. With VORTEX, deformations in both the vortex ring and the fluid surface will be measured closely for the first time in a parameters range that accurately resembles liquid atomization. The experimental apparatus will record images of the interactions for analysis after the payload has been returned to earth. The current design of the VORTEX payload consists of a fluid test cell with a vortex ring generator, digital imaging system, laser illumination system, computer based controller, batteries for payload power, and an array of housekeeping and payload monitoring sensors. It is a self-contained experiment and will be flown on board the Space Shuttle in a 5 cubic feet GAS canister. The VORTEX Project is entirely run by students at the University of Michigan but is overseen by a faculty advisor acting as the payload customer and the contact person with NASA. This paper summarizes both the technical and programmatic aspects of the VORTEX Project.

Description: This paper focuses on the flight results of the Cryogenic Two-Phase Flight Experiment (CRYOTP), which was a Hitchhiker based experiment that flew on the space shuttle Columbia in March of 1994 (STS-62). CRYOTP tested two new technologies for advanced cryogenic thermal control; the Space Heat Pipe (SHP), which was a constant conductance cryogenic heat pipe, and the Brilliant Eyes Thermal Storage Unit (BETSU), which was a cryogenic phase-change thermal storage device. These two devices were tested independently during the mission. Analysis of the flight data indicated that the SHP was unable to start in either of two attempts, for reasons related to the fluid charge, parasitic heat leaks, and cryocooler capacity. The BETSU test article was successfully operated with more than 250 hours of on-orbit testing including several cooldown cycles and 56 freeze/thaw cycles. Some degradation was observed with the five tactical cryocoolers used as thermal sinks, and one of the cryocoolers failed completely after 331 hours of operation. Post-flight analysis indicated that this problem was most likely due to failure of an electrical controller internal to the unit.

Description: The Capillary Pumped Loop Flight Experiment (CAPL) employs a passive two-phase thermal control system that uses the latent heat of vaporization of ammonia to transfer heat over long distances. CAPL was designed as a prototype of the Earth Observing System (EOS) instrument thermal control systems. The purpose of the mission was to provide validation of the system performance in micro-gravity, prior to implementation on EOS. CAPL was flown on STS-60 in February, 1994, with some unexpected results related to gravitational effects on two-phase systems. Flight test results and post flight investigations will be addressed, along with a brief description of the experiment design.

Description: This paper presents an overview of the Cryogenic Test Bed (CTB) experiments including experiment results, integration techniques used, and lessons learned during integration, test and flight phases of the Cryogenic Heat Pipe Flight Experiment (STS-53) and the Cryogenic Two Phase Flight Experiment (OAST-2, STS-62). We will also discuss the Cryogenic Flexible Diode Heat Pipe (CRYOFD) experiment which will fly in the 1996/97 time frame and the fourth flight of the CTB which will fly in the 1997/98 time frame. The two missions tested two oxygen axially grooved heat pipes, a nitrogen fibrous wick heat pipe and a 2-methylpentane phase change material thermal storage unit. Techniques were found for solving problems with vibration from the cryo-collers transmitted through the compressors and the cold heads, and mounting the heat pipe without introducing parasitic heat leaks. A thermally conductive interface material was selected that would meet the requirements and perform over the temperature range of 55 to 300 K. Problems are discussed with the bi-metallic thermostats used for heater circuit protection and the S-Glass suspension straps originally used to secure the BETSU PCM in the CRYOTP mission. Flight results will be compared to 1-g test results and differences will be discussed.

Description: With the recent focus on the needs of design and applications CFD, research groups have begun to address the traditional bottlenecks of grid generation and surface modeling. Now, a host of emerging technologies promise to shortcut or dramatically simplify the simulation process. This paper discusses the current status of these emerging technologies. It will argue that some tools are already available which can have positive impact on portions of the design cycle. However, in most cases, these tools need to be integrated into specific engineering systems and process cycles to be used effectively. The rapidly maturing status of unstructured and Cartesian approaches for inviscid simulations makes suggests the possibility of highly automated Euler-boundary layer simulations with application to loads estimation and even preliminary design. Similarly, technology is available to link block structured mesh generation algorithms with topology libraries to avoid tedious re-meshing of topologically similar configurations. Work in algorithmic based auto-blocking suggests that domain decomposition and point placement operations in multi-block mesh generation may be properly posed as problems in Computational Geometry, and following this approach may lead to robust algorithmic processes for automatic mesh generation.

Description: The last decade has witnessed a vigorous and sustained research effort on unstructured methods for computational fluid dynamics. Unstructured mesh generators and flow solvers have evolved to the point where they are now in use for design purposes throughout the aerospace industry. In this paper we survey the various mesh types, structured as well as unstructured, and examine their relative strengths and weaknesses. We argue that unstructured methodology does offer the best prospect for the next generation of computational fluid dynamics algorithms.

Description: A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.

Description: Grid related issues of the Chimera overset grid method are discussed in the context of a method of solution and analysis of unsteady three-dimensional viscous flows. The state of maturity of the various pieces of support software required to use the approach is considered. Current limitations of the approach are identified.

Description: This paper presents 'A view from the trenches' on CFD grid generation from a Pratt & Whitney perspective. We anticipate that other organizations have similar views. We focus on customer expectations and the consequent requirements. We enunciate a vision for grid generation, discuss issues that developers must recognize.

Description: This paper presents a perspective on the requirements that Computational Fluid Dynamics (CFD) technology must meet for its effective use in aerospace design. General observations are made on current aerospace design practices and deficiencies are noted that must be rectified for the U.S. aerospace industry to maintain its leadership position in the global marketplace. In order to rectify deficiencies, industry is transitioning to an integrated product and process development (IPPD) environment and design processes are undergoing radical changes. The role of CFD in producing data that design teams need to support flight vehicle development is briefly discussed. An overview of the current state of the art in CFD is given to provide an assessment of strengths and weaknesses of the variety of methods currently available, or under development, to produce aerodynamic data. Effectiveness requirements are examined from a customer/supplier view point with design team as customer and CFD practitioner as supplier. Partnership between the design team and CFD team is identified as an essential requirement for effective use of CFD. Rapid turnaround, reliable accuracy, and affordability are offered as three key requirements that CFD community must address if CFD is to play its rightful role in supporting the IPPD design environment needed to produce high quality yet affordable designs.

Description: The efforts in geometry modeling and grid generation at the NASA Lewis Research Center, as applied to the computational fluid dynamic (CFD) analysis of aeropropulsion systems, are presented. The efforts are mainly characterized by a focus on the analysis of components of an aeropropulsion system, which involve turbulent viscous flow with heat transfer and chemistry. Thus, this discussion will follow that characterization and will sequence through the components of typical propulsion systems consisting of inlets, compressors, combustors, turbines, and nozzles. For each component, some applications of CFD analysis will be presented to show how CFD is used to compute the desired performance information, how geometry modeling and grid generation are performed, and what issues have developed related to geometry modeling and grid generation. The discussion will illustrate the following needs related to geometry modeling and grid generation as observed in aeropropulsion analysis: (1) accurate and efficient resolution of turbulent viscous and chemically-reacting flowfields; (2) easy-to-use interfaces with CAD data for automated grid generation about complex geometries; and (3) automated batch grid generation software for use with design and optimization software.

Description: This viewgraph presentation discusses key features of current combustion system CFD modeling capabilities at GE Aircraft Engines provided by the CONCERT code; CONCERT development history; modeling applied for designing engine combustion systems; modeling applied to improve fundamental understanding; CONCERT3D results for current production combustors; CONCERT3D model of NASA/GE E3 combustor; HYBRID CONCERT CFD/Monte-Carlo modeling approach; and future modeling directions.

Description: Analytical investigation of dynamic stall on HAWT (horizontal-axis wind turbines) rotor loads was conducted. Dynamic stall was modeled using the Gormont approach on the MOD-2 rotor, treating the blade as a rigid body teetering about a fixed axis. Blade flapwise bending moments at station 370 were determined with and without dynamic stall for spatial variations in local wind speed due to wind shear and yaw. The predicted mean flapwise bending moments were found to be in good agreement with test results. Results obtained with and without dynamic stall showed no significant difference for the mean flapwise bending moment. The cyclic bending moments calculated with and without dynamic stall effects were substantially the same. None of the calculated cyclic loads reached the level of the cyclic loads measured on the MOD-2 using the Boeing five-minute-average technique.

Description: A procedure was developed to improve the turn-around time for computational fluid dynamics (CFD) simulations of an inlet-bleed problem involving oblique shock-wave/boundary-layer interactions on a flat plate with bleed into a plenum through one or more circular holes. This procedure is embodied in a preprocessor called AUTOMAT. With AUTOMAT, once data for the geometry and flow conditions have been specified (either interactively or via a namelist), it will automatically generate all input files needed to perform a three-dimensional Navier-Stokes simulation of the prescribed inlet-bleed problem by using the PEGASUS and OVERFLOW codes. The input files automatically generated by AUTOMAT include those for the grid system and those for the initial and boundary conditions. The grid systems automatically generated by AUTOMAT are multi-block structured grids of the overlapping type. Results obtained by using AUTOMAT are presented to illustrate its capability.

Description: The role of Computational Fluid Dynamics (CFD) at Ames Research Center has expanded to address a broad range of aeronautical problems, including wind tunnel support, flight test support, design, and analysis. Balancing the requirements of each new problem against the available resources - software, hardware, time, and expertise - is critical to the effective use of CFD. Several case studies of recent applications highlight the depth of CFD capability at Ames, the tradeoffs involved in various approaches, and lessons learned in the use of CFD as an engineering tool.

Description: When computing the flow around complex three dimensional configurations, the generation of the mesh is the most time consuming part of any calculation. With some meshing technologies this can take of the order of a man month or more. The requirement for a number of design iterations coupled with ever decreasing time allocated for design leads to the need for a significant acceleration of this process. Of the two competing approaches, block-structured and unstructured, only the unstructured approach will allow fully automatic mesh generation directly from a CAD model. Using this approach coupled with the techniques described in this paper, it is possible to reduce the mesh generation time from man months to a few hours on a workstation. The desire to closely couple a CFD code with a design or optimization algorithm requires that the changes to the geometry be performed quickly and in a smooth manner. This need for smoothness necessitates the use of Bezier polynomials in place of the more usual NURBS or cubic splines. A two dimensional Bezier polynomial based design system is described.

Description: The results of a first order perfect gas correction for the effects of the boundary layer formation within expansion tubes with nozzles are presented. The analytical model developed to describe the boundary layer formation within the expansion tube and an expansion nozzle located at the end of the acceleration tube is based on the Karman integral equations. The results of this analytical model are compared with experimental data from an expansion diffuser. The model provides a useful tool for the preliminary design of nozzles for such facilities.

Description: In this paper the capabilities of an automated CFD system which is currently available at NLR are demonstrated. Transonic flow around the AS28G wing/body configuration and hypersonic flow through a generic three-dimensional mixed-compression airbreathing inlet are simulated. An assessment of the level of automation of the current CFD-system is made. The problem-turnaround time lies within the order of a week for both applications.

Description: Detailed simulations of viscous flows in complicated geometries pose a significant challenge to current capabilities of Computational Fluid Dynamics (CFD). To enable routine application of CFD to this class of problems, advanced methodologies are required that employ (a) automated grid generation, (b) adaptivity, (c) accurate discretizations and efficient solvers, and (d) advanced software techniques. Each of these ingredients contributes to increased accuracy, efficiency (in terms of human effort and computer time), and/or reliability of CFD software. In the long run, methodologies employing structured grid systems will remain a viable choice for routine simulation of flows in complex geometries only if genuinely automatic grid generation techniques for structured grids can be developed and if adaptivity is employed more routinely. More research in both these areas is urgently needed.

Description: A coordinated effort has been underway over the past four years to elevate unstructured-grid methodology to a mature level. The goal of this endeavor is to provide a validated capability to non-expert users for performing rapid aerodynamic analysis and design of complex configurations. The Euler component of the system is well developed, and is impacting a broad spectrum of engineering needs with capabilities such as rapid grid generation and inviscid flow analysis, inverse design, interactive boundary layers, and propulsion effects. Progress is also being made in the more tenuous Navier-Stokes component of the system. A robust grid generator is under development for constructing quality thin-layer tetrahedral grids, along with a companion Navier-Stokes flow solver. This paper presents an overview of this effort, along with a perspective on the present and future status of the methodology.

Description: A discussion of the strengths and weaknesses of overset composite grid and solution technology is given, along with a sampling of current work in the area. Major trends are identified, and the observation is made that generalized and hybridized overset methods provide a natural framework for combining disparate mesh types and physics models. Because of this, the author concludes that overset methods will be the foundation for the general purpose computational fluid dynamics programs of the future.

Description: The status of CFD methods based on the use of block-structured grids for analyzing viscous flows over complex configurations is examined. The objective of the present study is to make a realistic assessment of the usability of such grids for routine computations typically encountered in the aerospace industry. It is recognized at the very outset that the total turnaround time, from the moment the configuration is identified until the computational results have been obtained and postprocessed, is more important than just the computational time. Pertinent examples will be cited to demonstrate the feasibility of solving flow over practical configurations of current interest on block-structured grids.

Description: This viewgraph presentation discusses turbulence modeling requirements, development philosophy, and approach; two major areas of concentration (high speed and low speed turbulence modeling); high speed turbulence modeling; compressibility effects; turbulence models adapted to USA code; M = 9.2 flat plate flow; Mach 7.05 flow over axisymmetric flare; Mach 8.6 flow over cold wall edge; low speed turbulence modeling; turbulence models being assessed; turbulence model deck structure and integration with Navier-Stokes solver; nonlinear algebraic-stress model; rotation modified k-epsilon model; and Reynolds stress model.

Description: This viewgraph presentation covers gas turbine combustor flow physics, turbulence model investigations, turbulent combustion modeling, and present status and future needs.

Description: This viewgraph presentation discusses geometry and flow configuration, effect of y+ on heat transfer computations, standard and extended k-epsilon turbulence model results with wall function, low-Re model results (the Lam-Bremhorst model without wall function), a criterion for flow reversal in a radially rotating square duct, and a summary.

Description: This viewgraph presentation demonstrates that computationally efficient k-l and k-kl turbulence models have been developed and implemented at Lockheed Fort Worth Company. Many years of experience have been gained applying two equation turbulence models to complex three-dimensional flows for design and analysis.

Description: The objective of this viewgraph presentation is to evaluate turbulence models for integrated aircraft components such as the forebody, wing, inlet, diffuser, nozzle, and afterbody. The one-equation models have replaced the algebraic models as the baseline turbulence models. The Spalart-Allmaras one-equation model consistently performs better than the Baldwin-Barth model, particularly in the log-layer and free shear layers. Also, the Sparlart-Allmaras model is not grid dependent like the Baldwin-Barth model. No general turbulence model exists for all engineering applications. The Spalart-Allmaras one-equation model and the Chien k-epsilon models are the preferred turbulence models. Although the two-equation models often better predict the flow field, they may take from two to five times the CPU time. Future directions are in further benchmarking the Menter blended k-w/k-epsilon and algorithmic improvements to reduce CPU time of the two-equation model.

Description: This viewgraph presentation summarizes some CFD experience at GE Aircraft Engines for flows in the primary gaspath of a gas turbine engine and in turbine blade cooling passages. It is concluded that application of the standard k-epsilon turbulence model with wall functions is not adequate for accurate CFD simulation of aerodynamic performance and heat transfer in the primary gas path of a gas turbine engine. New models are required in the near-wall region which include more physics than wall functions. The two-layer modeling approach appears attractive because of its computational complexity. In addition, improved CFD simulation of film cooling and turbine blade internal cooling passages will require anisotropic turbulence models. New turbulence models must be practical in order to have a significant impact on the engine design process. A coordinated turbulence modeling effort between NASA centers would be beneficial to the gas turbine industry.

Description: Program goals at the Center for Modeling of Turbulence and Transition (CMOTT), NASA Lewis Research Center, are (1) to develop reliable turbulence (including bypass transition) and combustion models for complex flows in propulsion systems and (2) to integrate developed models into deliverable CFD tools for propulsion systems in collaboration with industry. This viewgraph presentation covers the following topics: development of turbulence and combustion models; collaboration with industry and technology transfer; isotropic eddy viscosity models; algebraic Reynolds stress models; scalar turbulence models; second order closure models; multiple scale k-epsilon models; and PDF modeling of turbulent reacting flows.

Description: An assessment of two unstructured methods is presented in this paper. A tetrahedral unstructured method USM3D, developed at NASA Langley Research Center is compared to a Cartesian unstructured method, SPLITFLOW, developed at Lockheed Fort Worth Company. USM3D is an upwind finite volume solver that accepts grids generated primarily from the Vgrid grid generator. SPLITFLOW combines an unstructured grid generator with an implicit flow solver in one package. Both methods are exercised on three test cases, a wing, and a wing body, and a fully expanded nozzle. The results for the first two runs are included here and compared to the structured grid method TEAM and to available test data. On each test case, the set up procedure are described, including any difficulties that were encountered. Detailed descriptions of the solvers are not included in this paper.

Description: The past ten years have seen steady progress in surface modeling procedures, and wholesale changes in grid generation technology. Today, it seems fair to state that a satisfactory grid can be developed to model nearly any configuration of interest. The issues at present focus on operational concerns such as cost and quality. Continuing evolution of the engineering process is placing new demands on the technologies of surface modeling and grid generation. In the evolution toward a multidisciplinary analysis-bascd design environment, methods developed for Computational Fluid Dynamics are finding acceptance in many additional applications. These two trends, the normal evolution of the process and a watershed shift toward concurrent and multidisciplinary analysis, will be considered in assessing current capabilities and needed technological improvements.

Description: This viewgraph presentation discusses what models are used in this package and what their advantages and disadvantages are, how the probability density function (PDF) model is implemented and the features of the program, and what can be expected in the future from the NASA Lewis PDF code.

Description: This viewgraph presentation discusses an extension of the probability density function (PDF) method to the modeling of spray flames to evaluate the limitations and capabilities of this method in the modeling of gas-turbine combustor flows. The comparisons show that the general features of the flowfield are correctly predicted by the present solution procedure. The present solution appears to provide a better representation of the temperature field, particularly, in the reverse-velocity zone. The overpredictions in the centerline velocity could be attributed to the following reasons: (1) the use of k-epsilon turbulence model is known to be less precise in highly swirling flows and (2) the swirl number used here is reported to be estimated rather than measured.

Description: Viewgraphs are presented on computation of turbulent combustion, governing equations, closure problem, PDF modeling of turbulent reactive flows, validation cases, current projects, and collaboration with industry and technology transfer.

Description: Viewgraphs are presented on modeling turbulent reacting flows, regimes of turbulent combustion, regimes of premixed and regimes of non-premixed turbulent combustion, chemical closure models, flamelet model, conditional moment closure (CMC), NO(x) emissions from turbulent H2 jet flames, probability density function (PDF), departures from chemical equilibrium, mixing models for PDF methods, comparison of predicted and measured H2O mass fractions in turbulent nonpremixed jet flames, experimental evidence of preferential diffusion in turbulent jet flames, and computation of turbulent reacting flows.

Description: This viewgraph presentation discusses project description, turbulence models, and computational engine and results for second-order closures for compressible turbulence.

Description: This viewgraph presentation gives a profile of Advanced Scientific Computing (ASC) Ltd., applications, clients and clients' needs, ASC's directions, and how the Center for Modeling of Turbulence and Transition (CMOTT) can help.

Description: This viewgraph presentation discusses the following: STAR-CD computational features; STAR-CD turbulence models; common features of industrial complex flows; industry-specific CFD development requirements; applications and experiences of industrial complex flows, including flow in rotating disc cavities, diffusion hole film cooling, internal blade cooling, and external car aerodynamics; and conclusions on turbulence modeling needs.

Description: This viewgraph presentation provides a brief review of two-equation eddy-viscosity models (TEM's) from the perspective of applied CFD. It provides objective assessment of both well-known and newer models, compares model predictions from various TEM's with experiments, identifies sources of modeling error and gives historical perspective of their effects on model performance and assessment, and recommends directions for future research on TEM's.

Description: This viewgraph presentation discusses the following: introduction to CFD Research Corporation; experiences with two-equation models - models used, numerical difficulties, validation and applications, and strengths and weaknesses; and answers to three questions posed by the workshop organizing committee - what are your customers telling you, what are you doing in-house, and how can NASA-CMOTT (Center for Modeling of Turbulence and Transition) help.

Description: This viewgraph presentation concludes that a Monte Carlo probability density function (PDF) solution successfully couples with an existing finite volume code; PDF solution method applied to turbulent reacting flows shows good agreement with data; and PDF methods must be run on parallel machines for practical use.

Description: This viewgraph presentation discusses (1) turbulence modeling: challenges in turbulence modeling, desirable attributes of turbulence models, turbulence models in FLUENT, and examples using FLUENT; and (2) combustion modeling: turbulence-chemistry interaction and FLUENT equilibrium model. As of now, three turbulence models are provided: the conventional k-epsilon model, the renormalization group model, and the Reynolds-stress model. The renormalization group k-epsilon model has broadened the range of applicability of two-equation turbulence models. The Reynolds-stress model has proved useful for strongly anisotropic flows such as those encountered in cyclones, swirlers, and combustors. Issues remain, such as near-wall closure, with all classes of models.

Description: This viewgraph presentation discusses joint velocity-scalar PDF method; turbulent combustion modeling issues for gas turbine combustors; PDF calculations for a recirculating flow; stochastic dissipation model; joint PDF calculations for swirling flows; spray calculations; reduced kinetics/manifold methods; parallel processing; and joint PDF focus areas.

Description: An algorithm has been developed for time-dependent forced convective diffusion-reaction having convection by a recirculating flow field within the drop that is hydrodynamically coupled at the interface with a convective external flow field that at infinity becomes a uniform free-streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet, or reactions can take place in both phases. The algorithm has been implemented, and for comparison results are shown here for the case of no reaction in either phase and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Description: The application of spectral methods, using a Chebyshev collocation scheme, to solve hydrodynamic stability problems is demonstrated on the Benard problem. Implementation of the Chebyshev collocation formulation is described. The performance of the spectral scheme is compared with that of a 2nd order finite difference scheme. An exact solution to the Marangoni-Benard problem is used to evaluate the performance of both schemes. The error of the spectral scheme is at least seven orders of magnitude smaller than finite difference error for a grid resolution of N = 15 (number of points used). The performance of the spectral formulation far exceeded the performance of the finite difference formulation for this problem. The spectral scheme required only slightly more effort to set up than the 2nd order finite difference scheme. This suggests that the spectral scheme may actually be faster to implement than higher order finite difference schemes.

Description: The effects of different turbulence boundary conditions were examined for two classical flows: a turbulent plane free shear layer and a flat plate turbulent boundary layer with zero pressure gradient. The flow solver used was DTNS, an incompressible Reynolds averaged Navier-Stokes solver with k-epsilon turbulence modeling, developed at the U.S. Navy David Taylor Research Center. Six different combinations of turbulence boundary conditions at the inflow boundary were investigated: In case 1, 'exact' k and epsilon profiles were used; in case 2, the 'exact' k profile was used, and epsilon was extrapolated upstream; in case 3, both k and epsilon were extrapolated; in case 4, the turbulence intensity (I) was 1 percent, and the turbulent viscosity (mu(sub t)) was equal to the laminar viscosity; in case 5, the 'exact' k profile was used and mu(sub t) was equal to the laminar viscosity; in case 6, the I was 1 percent, and epsilon was extrapolated. Comparisons were made with experimental data, direct numerical simulation results, or theoretical predictions as applicable. Results obtained with DTNS showed that turbulence boundary conditions can have significant impacts on the solutions, especially for the free shear layer.

Description: The exchange of models is one of the most serious problems currently encountered in the practice of spacecraft thermal analysis. Essentially, the problem originates in the diversity of computing environments that are used across different sites, and the consequent proliferation of native tool formats. Furthermore, increasing pressure to reduce the development's life cycle time has originated a growing interest in the so-called spacecraft concurrent engineering. In this context, the realization of the interdependencies between different disciplines and the proper communication between them become critical issues. The use of a neutral format represents a step forward in addressing these problems. Such a means of communication is adopted by consensus. A neutral format is not directly tied to any specific tool and it is kept under stringent change control. Currently, most of the groups promoting exchange formats are contributing with their experience to STEP, the Standard for Exchange of Product Model Data, which is being developed under the auspices of the International Standards Organization (ISO 10303). This paper presents the different efforts made in Europe to provide the spacecraft thermal analysis community with a Thermal Neutral Format (TNF) based on STEP. Following an introduction with some background information, the paper presents the characteristics of the STEP standard. Later, the first efforts to produce a STEP Spacecraft Thermal Application Protocol are described. Finally, the paper presents the currently harmonized European activities that follow up and extend earlier work on the area.

Description: Using global interpolation functions (GIF's) boundary element solutions are obtained for two-dimensional laminar flows. Two schemes are proposed for handling the convective terms. The first treats convection as a forcing function, and converts the flow equations to pseudo-Poisson equations. In the second scheme, some convective effect is incorporated into the fundamental solution used in constructing the pertinent integral equations. The lid-driven cavity flow is selected as the benchmark problem.

Description: A two dimensional finite volume method is used to predict the film coefficients in the transitional flow region (laminar or turbulent) for the radiator panel tubes. The code used to perform this analysis is CAST (Computer Aided Simulation of Turbulent Flows). The information gathered from this code is then used to augment a Sinda85 model that predicts overall performance of the radiator. A final comparison is drawn between the results generated with a Sinda85 model using the Sinda85 provided transition region heat transfer correlations and the Sinda85 model using the CAST generated data.

Description: This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

Description: This paper presents a numerical scheme, based on the finite element method, to solve strongly coupled fluid flow and heat transfer problems. The surface radiation effect for gray, diffuse and isothermal surfaces is considered. A procedure for obtaining the view factors between the radiating surfaces is discussed. The overall solution strategy is verified by comparing the available results with those obtained using this approach. An analysis of a thermosyphon is undertaken and the effect of considering the surface radiation is clearly explained.

Description: An algorithm has been developed for the forced convective diffusion-reaction problem for convection inside and outside a droplet by a recirculating flow field hydrodynamically coupled at the droplet interface with an external flow field that at infinity becomes a uniform streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet or reactions can take place in both phases. The algorithm has been implemented and results are shown here for the case of no reaction and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Description: The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon, RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.

Description: The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.