ALBERT

Description: Transient solutions were obtained for a square region of heat conducting semitransparent material cooling by thermal radiation. The region is in a vacuum environment, so energy is dissipated only by radiation from within the medium leaving through its boundaries. The effect of heat conduction during the transient is to partially equalize the internal temperature distribution. As the optical thickness of the region is increased, the temperature gradients increase near the boundaries and corners, unless heat conduction is large. The solution procedure must provide accurate temperature distributions in these regions to prevent error in the calculated radiation losses. Two-dimensional numerical Gaussian integration is used to obtain the local radiative source term. A finite difference procedure with variable space and time increments is used to solve the transient energy equation. Variable spacing was used to concentrate grid points in regions with large temperature gradients.

Description: Small hydrocarbon grains in the vicinity of a supernova could be annealed by the absorption of several far-ultraviolet photons to produce the tiny diamonds found in meteorites. These freshly-synthesized diamond grains would be bombarded by the heavy ions and neutrals in the supernovae outflow and would thereby acquire the distinctive noble-gas isotopic signature by which they were first isolated. Only diamonds formed relatively close to supernovae would acquire such a signature, since grains formed farther out would be subjected to a much diluted and less energetic plasma environment.

Description: IUE observations taken during the 1988 January normal outburst of Z Cha are presented and a detailed comparison with the 1987 April superoutburst is made. The most important difference from the superoutburst is that the normal outburst continuum flux shows less than 10 percent orbital variation away from the eclipse, implying that there is no 'cool' bulge on the disk to occult the brighter inner disk periodically. The implications for the outburst mechanism in the types of outburst are discussed. The evolution of the continuum flux distribution and emission-line fluxes, the modulation of the continuum and line fluxes with orbital phase, and the behavior of the mideclipse spectral during normal outburst are investigated.

Description: The Energetic Gamma-Ray Experiment Telescope (EGRET) on the Compton Gamma-Ray Observatory observed the Cygnus region of the Galaxy during the periods 1991 May 30-June 8 and 1991 August 8-15. We report on a periodicity analysis of the gamma rays originating from the direction of Cygnus X-3. This analysis showed no evidence of periodic modulation during these observations of the gamma-ray emission at the approximately 4.8 hr period observed at X-ray wavelengths.

Description: In the present paper, direct numerical methods by which to simulate the spatially developing free shear flows in the transitional region are described and the numerical results of a spatially developing plane wake are presented. The incompressible time-dependent Navier-Stokes equations were solved using Pade finite difference approximations in the streamwise direction, a mapped pseudospectral Fourier method in the cross-stream direction, and a third-order compact Runge-Kutta scheme for time advancement. The unstable modes of the Orr-Sommerfeld equations were used to perturb the inlet of the wake. Statistical analyses were performed and some numerical results were compared with experimental measurements. When only the fundamental mode is forced, the energy spectra show amplification of the fundamental and its higher harmonics. In this case, unperturbed alternate vortices develop in the saturation region of the wake. The phase jitter around the fundamental frequency plays a critical role in generating vortices of random shape and spacing. Large- and small-scale distortions of the fundamental structure are observed. Pairing of vortices of the same sign is observed, as well as vortex coupling of vortices of the opposite sign.

Description: Inflation provides very strong motivation for a flat Universe, Harrison-Zel'dovich (constant-curvature) perturbations, and cold dark matter. However, there are a number of cosmological observations that conflict with the predictions of the simplest such model: one with zero cosmological constant. They include the age of the Universe, dynamical determinations of Omega, galaxy-number counts, and the apparent abundance of large-scale structure in the Universe. While the discrepancies are not yet serious enough to rule out the simplest and most well motivated model, the current data point to a best-fit model with the following parameters: Omega(sub B) approximately equal to 0.03, Omega(sub CDM) approximately equal to 0.17, Omega(sub Lambda) approximately equal to 0.8, and H(sub 0) approximately equal to 70 km/(sec x Mpc) which improves significantly the concordance with observations. While there is no good reason to expect such a value for the cosmological constant, there is no physical principle that would rule out such.

Description: This book first reviews the overall aspects and background information related to thermal radiation heat transfer and incorporates new general information, advances in analytical and computational techniques, and new reference material. Coverage focuses on radiation from opaque surfaces, radiation interchange between various types of surfaces enclosing a vacuum or transparent medium, and radiation including the effects of partially transmitting media, such as combustion gases, soot, or windows. Boundary conditions and multiple layers are discussed with information on radiation in materials with nonunity refractive indices.

Description: The present approach to the prediction of instability generation that is due to the interaction of freestream disturbances with regions of subscale variations in surface boundary conditions can account for the finite Reynolds number effects, while furnishing a framework for the study of receptivity in compressible flow and in 3D boundary layers. The approach is illustrated for the case of Tollmien-Schlichting wave generation in a Blasius boundary layer, due to the interaction of a freestream acoustic wave with a localized wall inhomogeneity. Results are presented for the generation of viscous and inviscid instabilities in adverse pressure-gradient boundary layers, supersonic boundary layer instabilities, and cross-flow vortex instabilities.

Description: The hard X-ray spectrum of AGNs is nonthermal, probably arising from an electron-positron pair cascade, with some emission reflected off relatively cold matter. There has been interest in models on which protons are accelerated and create relativistic electrons on interaction with a local radiation field. It is shown here that a sufficient column density of protons can lead to runaway pair production: photons generated by the relativistic pairs are the targets for the protons to produce more pairs. This can produce X-ray fluxes with the characteristics observed in AGN. The model predicts the maximum ratio of luminosity to source size as well as their spectrum in the early phases. The same mechanism may also be able to create the knots of synchrotron-radiating pair plasma seen in sources such as 3C273.

Description: The contribution of scattered light to the total nebular emission is determined on the basis of linear polarization measurements at 1.25, 1.65, and 2.2 microns of the visual reflection nebulae NGC 7023 and NGC 2023. The percentage polarization of NGC 7023 slowly increases from 0.3 to 1 micron, with peak polarizations of up to 26 percent at 1.25 micron, then rapidly decreases, with values of 4-7 percent at 2.2 microns. This is interpreted as implying that scattered starlight contributes most to the SW emission, while unpolarized emission from small grains or large molecules dominates at longer wavelengths. IR polarization and surface brightness measurements are combined to derive the intensity of scattered light, which is then compared with scattering models. While the near-IR emission of both NGC 2023 and NGC 7023 is dominated by small-grain or large-molecule emission, IR scattered light plays a larger role in NGC 2023 than in NGC 7023.

Description: The existence of the periodic component seen in NGC 6814 with Exosat at 12,000 +/- 100 s is confirmed by a power spectrum and folded light curve analysis of unevenly sampled Ginga data. A comparison of the power spectra produced from simulated light curves with that observed enables the intrinsic shape of the power spectrum of the source to be determined despite the distortions introduced by the window function. The best estimate for the period is 12,132 +/- 3 s, where the error is that derived from simulations. An upper limit to the rate of change of period of about 10 exp -9 is inferred if the light curves are truly phase-coherent, but as this is not required by the data, the conservative upper limit is not greater than 5 x 10 exp -7. The large amount of power in the periodic component and its stability both suggest occultation of the source as its origin.

Description: The collapse of cold, initially spherical systems with varying degrees of central condensation is investigated. The way in which the final shape of a collapsing system depends on the initial density law is examined. For an initial stellar number density rho varies as r exp -n, where n is in the range 0-2.5, the final, nearly prolate shape is given by a/c is approximately equal to 1.28(1 + 0.16 n), where a/c is the ratio of long to short axes of the inertia ellipsoid computed from the moment of inertia tensor of the most tightly bound 80 percent of the mass. The properties associated with the final states in the present computations are also studied. The collapsing systems develop an anisotropic halo dominated by radial orbits surrounding an isotropic core as predicted by Burkert (1990).

Description: Gamma radiation above 100 MeV in energy has been detected from the radio pulsar PSR1706-44. The gamma emission forms a single broad peak within the pulsar period of 102 ms, in contrast to the two narrow peaks seen in the other three known high-energy gamma-ray pulsars. The emission mechanism in all cases is probably the same, the differences arising from the geometry of the magnetic and rotation axes and the line of sight. Gamma-ray emission accounts for as much as 1 percent of the total neutron star spindown energy in these pulsars, much more than emerges at optical or radio frequencies. Thus, study of this emission is important in understanding pulsar emission and evolution.

Description: The properties and evolution of the major outburst experienced by Comet Halley at 14 AU during February-March 1991 are studied. It is found that the observed halo is a segment of a conical surface populated by solid particles ejected from a localized, temporarily activated region on the sunlit hemisphere of the spinning nucleus. The total mass of the ejecta was at least 10 exp 12 g and CO was probably the prime driver, accelerating the smallest grains to a terminal velocity of about 45 m/s. The mass loading of the gas flow by particulate matter is enormous, with the mass production rate exceeding the expected production rate of CO by a factor of several tens.

Description: A detailed abundance analysis of Cr, Fe, Mn, Ca, and Ti in the Ap star 53 Cam from IUE high-resolution spectra has been performed. It is found that UV lines give significantly lower abundance than visible lines for Cr and Fe. Theoretical spectra of Cr II line calculated with the stratified chromium abundance distributions predicted a priori by the diffusion-mass loss model match well both visible and UV lines. The spectral shape of the Ca II K line and its variation are well accounted for by the diffusion model with mass loss. A spectrum synthesis with a stratified abundance well reproduces the observations. These findings provide convincing evidence for stratification in the photosphere of 53 Cam.

Description: A model for the IR emission and extinction properties of small dust particles is used here to compute the IR emissions in the IRAS photometric bands for a set of isolated, spherical, and nonhomogeneous clouds heated by the local interstellar radiation field. It is predicted that the IR limb brightening (LB) in the IRAS photometric bands caused by selective absorption of UV photons of different energies in the cloud generally happens for central extinctions over 4 mag at 12 microns and greater than 10 mag at 100 microns and intermediate extinctions for the other IRAS bands. The surface brightness in the four IRAS bands is limited to about 0.4, 0.6, 3, and 10 MJy/sr at 12, 25, 60, and 100 microns respectively when no strong density discontinuity is present at the cloud edge.

Description: The growth and development of a horseshoe vortex system in an incompressible, three-dimensional turbulent junction flow were investigated experimentally. A streamlined cylinder mounted with its axis normal to a flat surface was used to generate the junction vortex flow. The flow environment was characterized by a body Reynolds number of 183,000, based on the leading edge diameter of the streamlined cylinder. The study included surface flow visualizations, surface pressure measurements, and mean flow measurements of total pressure, static pressure, and velocity distributions in three planes around the base of the streamlined cylinder, and in two planes in the wake flow. Some characterizations of vortex properties based on the measured mean cross-flow velocity components are presented. The results show the presence of a single large, dominant vortex, with strong evidence of a very small corner vortex in the junction between the cylinder and the flat surface. The center of the dominant vortex drifts away from both the body and the flat surface as the flow develops along and downstream of the body. The growth and development of the core of the large, dominant vortex are documented.

Description: The nonlinear resonant-triad interaction, proposed by Raetz (1959), Craik (1971), and others for a Blasius boundary layer, is analyzed here for an adverse-pressure-gradient boundary layer. We assume that the adverse pressure gradient is in some sense weak and, therefore, that the instability growth rate is small. This ensures that there is a well-defined critical layer located somewhere within the flow and that the nonlinear interaction is effectively confined to that layer. The initial interaction is of the parametric resonance type, even when the modal amplitudes are all of the same order. This means that the oblique instability waves exhibit faster than exponential growth and that the growth rate of the two-dimensional mode remains linear. However, the interaction and the resulting growth rates become fully coupled, once oblique-mode amplitudes become sufficiently large, but the coupling terms are now quartic, rather than quadratic as in the Craik (1971) analysis. More importantly, however, new nonlinear interactions, which were not present in the Craik-type analyses, now come into play. These interactions eventually have a dominant effect on the instability wave development.

Description: Reacting free shear layers are of fundamental importance in many industrial systems including gas turbine combustors and rockets. Efficient propulsion systems are essential for air breathing supersonic ramjets in the high Mach number range. A limiting factor in these engines is the time for fuel and oxidizer to mix in the combustion chamber; for fast mixing, the flow must be vigorously turbulent which requires the laminar flow to be unstable. Understanding the stability characteristics of compressible reacting free shear layers is, therefore, very important and may allow one to control the flow. Low speed shear layers are highly unstable but, as chemical reaction and compressibility effects tend to stabilize them, it is important to investigate the stability of high speed reacting mixing layers. The latter consists of two fluid streams containing fuel and oxidizer respectively, and the conclusions are expected to apply, with quantitative modifications, to other shear flows, e.g., jets. Since low speed reacting cases have been studied earlier, we concentrate on the effects of Mach number and heat release. We are primarily interested in solving the stability problem over a large range of Mach number and heat release. In order to understand the effect of the heat release on the stability of this flow, one must first study the characteristics of the non-reacting flow. Inviscid theory is a reliable guide for understanding stability of compressible shear flows at moderate and large Reynolds numbers and is the basis for this work.

Description: Several direct numerical simulations of high-speed turbulent Couette flow were performed with a new spectral code. Mach numbers up to three and a Reynolds number of 3000 were used. A new time-integration scheme was developed to handle Mach numbers above 1.5, which require greater accuracy and stability than lower Mach numbers. At low Mach number, the large streamwise eddies found by M. J. Lee in high incompressible Couette flow simulations were reproduced. At higher Mach numbers these structures still exist, but they become considerably less organized (although the disorganization may be a function of the spanwise box size). While the same types of vortical structures seen in the incompressible flow are observed at higher Mach numbers, a new structure involving the divergence of the velocity is also observed. This structure is generally associated with low shear areas next to the walls, but it has not been determined whether it is a cause or an effect of the low shear. A 'nonphysical' simulation was performed to determine by what mechanism the Mach number affects the flow. It appears that pressure gradient (acoustic) effects are more important than variable viscosity effects in determining the wall shear, but the size of vortical structures is determined more by the local kinematic viscosity. Low-order mean statistics are provided to help quantify these effects.

Description: Many of the turbulent layers encountered in practical flows develop in adverse pressure gradients; hence, the dynamics of the thickening and possible separation of the boundary layer has important implications for design practices. What are the key physical processes that govern how a turbulent boundary layer responds to an adverse pressure gradient, and how should these processes be modeled? Despite the ubiquity of such flows in engineering and nature, these equations remain largely unanswered. The turbulence closure models presently used to describe these flows commonly use 'wall functions' that have ad hoc corrections for the effects of pressure gradients. There is, therefore, a practical and theoretical need to examine the effects of adverse pressure gradients on wall bounded turbulent flows in order to develop models based on sound physical principle. The evolution of a turbulent boundary layer on a flat wall with an externally imposed pressure gradient is studied.

Description: Advancing the knowledge and understanding of turbulence theory is addressed. Specific problems to be addressed will include studies of subgrid models to understand the effects of unresolved small scale dynamics on the large scale motion which, if successful, might substantially reduce the number of degrees of freedom that need to be computed in turbulence simulation.

Description: The increase in the range of length scales with increasing Reynolds number limits the direct simulation of turbulent flows to relatively simple geometries and low Reynolds numbers. However, since most flows of engineering interest occur at much higher Reynolds number than is currently within the capabilities of full simulation, prediction of these flow fields can only be obtained by solving some suitably-averaged set of governing equations. In the traditional Reynolds-averaged approach, the Navier-Stokes equations are averaged over time. This in turn yields correlations between various turbulence fluctuations. It is these terms, e.g. the Reynolds stresses, for which a turbulence model must be derived. Turbulence modeling of incompressible flows has received a great amount of attention in the literature. An area of research that has received comparatively less attention is the modeling of compressible turbulent flows. An approach to simulating compressible turbulence at high Reynolds numbers is through the use of Large-Eddy Simulation (LES). In LES the dependent variables are decomposed into a large-scale (resolved) component and a sub-grid scale component. It is the small-scale components of the velocity field which are presumably more homogeneous than the large scales and, therefore, more easily modeled. Thus, it seems plausible that simpler models, which should be more universal in character than those employed in second-order closure schemes, may be developed for LES of compressible turbulence. The objective of the present research, therefore, is to explore models for the Large-Eddy Simulation of compressible turbulent flows. Given the recent successes of Zeman in second order closure modeling of compressible turbulence, model development was guided by principals employed in second-order closures.

Description: With the recent revitalization of high speed flow research, compressibility presents a new set of challenging problems to turbulence researchers. Questions arise as to what extent compressibility affects turbulence dynamics, structures, the Reynolds stress-mean velocity (constitutive) relation, and the accompanying processes of heat transfer and mixing. In astrophysical applications, compressible turbulence is believed to play an important role in intergalactic gas cloud dynamics and in accretion disk convection. Understanding and modeling of the compressibility effects in free shear flows, boundary layers, and boundary layer/shock interactions is discussed.

Description: I discuss the most recent model of inflation. In first-order inflation the inflationary epoch is associated with a first-order phase transition, with the most likely candidate being GUT symmetry breaking. The transition from the false-vacuum inflationary phase to the true-vacuum radiation-dominated phase proceeds through the nucleation and percolation of true-vacuum bubbles. The first successful and simplest model of first-order inflation, extended inflation, is discussed in some detail: evolution of the cosmic-scale factor, reheating, density perturbations, and the production of gravitational waves both from quantum fluctuations and bubble collisions. Particular attention is paid to the most critical issue in any model of first-order inflation: the requirements on the nucleation rate to ensure a graceful transition from the inflationary phase to the radiation-dominated phase.

Description: Using the metric in the local inertial frame of the Earth, we calculate relativistic effects on the lunar orbit with the synodic month period. It is shown that such perturbations arise entirely from the gravitomagnetic components of the local metric which exist because of the relative motion of the sun with respect to the Earth. In the case of general relativity, the net perturbation has an amplitude of 3 cm for the lunar range.

Description: Because of the breakdown of the Einstein equivalence principle in the nonsymmetric gravitational theory (NGT) of Moffat, orthogonally polarized electromagnetic waves can propagate at different velocities in a gravitational field. Moffat has proposed that galactic dark matter, in the form of cosmions, may act as a significant source of gravity in the NGT. We discuss how observations of the highly polarized radiation from distant pulsars could provide significant limits on the strength of the coupling of cosmions in the NGT.

Description: Two approaches to accelerating the method of complete linearization for calculating NLTE model stellar atmospheres are suggested. The first one, the so-called Kantorovich variant of the Newton-Raphson method, consists of keeping the Jacobi matrix of the system fixed, which allows us to calculate the costly matrix inversions only a few times and then keep them fixed during the subsequent computations. The second method is an application of the Ng acceleration. Both methods are extremely easy to implement with any model atmosphere code based on complete linearization. It is demonstrated that both methods, and especially their combination, yield a rapidly and globally convergent algorithm, which takes 2 to 5 times less computer time, depending on the model at hand and the required accuracy, than the ordinary complete linearization. Generally, the time gain is more significant for more complicated models. The methods were tested for a broad range of atmospheric parameters, and in all cases they exhibited similar behavior. Ng acceleration applied on the Kantorovich variant thus offers a significant improvement of the standard complete-linearization method, and may now be used for calculating relatively involved NLTE model stellar atmospheres.

Description: We have constructed high resolution IRAS maps of the B209 star forming region in the Taurus complex. By combining survey scans with pointed observations using advanced techniques specially developed to handle IRAS data in crowded fields, we have constructed images that reveal two new premain-sequence (PMS) objects. One of these, IRAS 04114+2757G, has FIR colors close to those of typical T Tauri stars and is probably associated with a faint stellar object visible on the POSS prints. The other new PMS source, IRAS 04111+2800G, is a deeply embedded object not detected at 12 microns and not optically visible. The object is projected close to the center of an ammonia core. Our detection of compact blue and red shifted (C-12)O emission peaks around the position of IRAS 04111+2800G confirms the PMS nature and youth of this object. We have examined the reasons for the exclusion of IRAS 04111+2800G and IRAS 04114+2757G from the IRAS Point Source Catalog and estimated the incompleteness of the currently known sample of embedded PMS stars in Taurus. The implications for the determination of evolutionary time scales in low mass star formation are discussed.

Description: The Giotto spacecraft carried two different instruments - the JPA and the IMS - for the observation of hot ions in the coma of P/Halley. Although there are many similarities in the time and distance profiles of the plasma flow parameters (bulk velocity, number density, and temperature) computed from the two data sets, there are also some significant differences, especially at cometocentric distances less than 500,000 km. The principal discrepancies between the JPA results presented by Formisano et al. (1990) and the IMS observations are: (1) the IMS did not detect the levelling off of the speed and temperature profiles that Formisano et al. interpreted as flow stabilization; (2) the IMS detected differential north-south flow between the solar wind and cometary ions for only a brief interval when the magnetic field was oriented nearly southward, whereas Formisano et al. reported more extensive differential north-south flow that was independent of the direction of the field; (3) the JPA ion densities were factors of 2 to 4 higher than the IMS ion densities which, in turn, were an order of magnitude greater than theoretical values.

Description: We present new theoretical calculations for the red wing of the Lyman-alpha profile. Close collisions with neutral and ionized hydrogen lead to the formation of the pseudomolecules H-H and H-H(+) with the appearance of satellite features near 1600 and 1400 A. The calculations include multiperturber effects, which are responsible for the formation of H3 and H3 with features near 1950 and 2600 A. The theoretical absorption profiles are included in stellar atmosphere codes and used to predict synthetic spectra for DA white dwarfs of intermediate temperatures (20,000 to 8000 K). These new calculations offer a unique opportunity to determine accurate effective temperatures and surface gravities for the variable ZZ Ceti stars.

Description: We present new observations of the H-alpha line profile for five main-sequence A-type stars, where very high SNRs were achieved. A search for weak asymmetries was carried out to detect a stellar wind, but the many telluric absorption lines in this range prevent us from taking the full benefit of the high SNRs. This situation was improved by modeling the telluric absorptions to remove them from the observed spectra, but the H-alpha profiles were nevertheless found to be quite symmetric. We calculated the profile of the H-alpha line for a grid of model atmospheres of an A-type star including a weak wind, in order to assess an upper limit on the mass loss rate. The asymmetry is found to be sensitive to the velocity law and to the turbulent velocity of the wind, but its first moment depends much less on them. The upper limits deduced on the mass loss rate are between 1 and 2 x 10 exp -10 solar mass/yr. An improvement up to a factor 10 could be anticipated if new observations could be secured from a much drier location, achieving similar SNRs.

Description: The LBV characteristics of the new LBV candidate WRA 751 suggested by Hu et al. (1990) are studied in greater detail. A comparison is made with observational characteristics of the well-known LBVs Ag Car and HR Car using photometric and spectroscopic data of Hu et al. supplemented with an IUE spectrum and direct CCD images in H-alpha, blue, and red. It is found that the spectra of WRA 751 and HR Car discussed in this paper are similar to the minimum-phase spectrum of AG Car and therefore resemble the Of/WN 9 spectral type. The characteristics of the optical, near- and far-UV (IUE) spectra of the Fe II spectrum in these wavelength regions is due to continuum fluorescence. The H, He I, Fe II and forbidden Fe II emission lines are formed in different parts of the stars' envelopes, corresponding to their expanding velocities. It is concluded that WRA 751 has properties which fit well the LBV properties at minimum-brightness phase, and possesses circumstellar material with observational properties similar to those of the known LBVs AG Car and HR Car.

Description: An account is given of interface-driven motions of drops and bubbles. It is shown that even in the simplest cases, theory predicts exotic flow topologies. Attention is given to several unsolved problems that must be addressed both theoretically and experimentally.

Description: The first results of the STARDUST project, aimed at producing and analyzing cosmic-dust analog materials in microgravity conditions, are summarized. The discussion covers the purpose of the investigation, cosmic-dust formation and properties, previous simulations of cosmic-dust formation, the current approach, the microgravity experimental apparatus, and potential advantages of studying dust formation under microgravity conditions.

Description: Three molecular clouds of Sgr B2 have been mapped in the 1(10) - 1(01) orthotransition of C3H2 at 18.3 GHz. The distribution of C3H2 and the C3H2 optical depth of the three clouds, centered near 50, 65, and 80 km/s, are reported.

Description: The synchrotron self-Compton (SSC) emission model for gamma-ray bursts is capable of producing narrow line features for a variety of modelled field strengths, primary electron injection distributions and burst luminosities. Multiple resonant scattering with synchrotron continuum photons efficiently traps and cools pairs in the ground state to an average energy where the Compton energy loss rate is zero. Annihilation between pairs in these cooled distributions can be very efficient. For isotropic injection of primary electrons, there is an anti-correlation of the intensity of the angular-binned emission between the portions of the continuum spectrum below the cyclotron peak and above the annihilation line feature. Small-angle emission dominates the continuum above the annihilation line feature and is smooth through the line up to the pair production cut-off, which can be above several MeV for small enough angles. The angle-averaged annihilation line is made of components which can have peak emission centered at energies away from 511 keV, due to Doppler shifting. For beamed injection,the annihilation line breaks up into relatively narrow components in the angular emission, which sum into a broad feature in the angle-averaged emission.

Description: With improved data from BATSE and other instruments, it is important to develop a range of diagnostic tools to link gamma-ray burst observations with theory. I will review some of the physical processes which may take place to form the spectrum of gamma-ray burst sources, assuming that the bursts originate on strongly magnetized neutron stars. The important diagnostics that these processes provide to probe the emission region and how they might be used to interpret observed spectra will also be discussed.

Description: Data on the velocity distributions of pickup ions in the vicinity of Comet Halley that were obtained by the Giotto spacecraft are presented. Combination of data from two different instruments allows a comparative study of the pitch-angle scattering rates and velocity diffusion rates of protons and water-group ions.

Description: It has been suggested that the presence of red giants in the vicinity of an active nucleus can contribute significantly to the line emission by the reprocessing of the continuum radiation onto the red giant's surface or in the stellar wind that usually accompanies these stars. We present herein the line profiles and the corresponding transfer functions for this model for the line emission, using a realistic phase space distribution function for the stars in the vicinity of the active nucleus.

Description: We review recent hard X-ray and gamma-ray observations of Cen A, NGC 4151, 3C273, and 3C279. The new data showed remarkable variability of these objects in the high-energy gamma-ray regime. Recent observations by GRANAT show an unexpected sharp break at about 50 keV in the 1990 spectrum of NGC 4151. Recent observations by the Compton Observatory show intense GeV emission from 3C279 in June 1991. We present composite energy spectra of these objects from the radio to gamma-ray energies.

Description: The general properties of the high energy AGN continuum are reviewed with particular emphasis on the connection between radio loud and radio quiet AGN in the broad context of the dynamics of accretion onto a black hole. Arguments are provided indicating the possibility that the emission from radio loud AGN originates at distances much larger than a few Schwarzschild radii from the black hole and hence the role of pairs in defining the observed spectrum is limited. A connection is made between the radio and high energy continuum of AGN.

Description: Contributions of amateur astronomers to research on dwarf novae, which are based on carefully monitoring the outburst behavior of these objects, are reviewed. These contributions range from scheduling of observations to the observational basis for research on the dwarf nova outburst mechanism. It is suggested, that, with better equipment, observations of orbital light variations in dwarf novae might be performed by amateur astronomers.

Description: We discuss the results of a survey of Quasar radio structures over redshifts from 0.6 to 3.7. There are clear evolutionary trends in size and luminosity, which suggest that the duty cycle of individual Quasars has increased over cosmic time. This affects source count statistics and gives clues on the evolution of Quasar environments.

Description: Direct numerical simulations were made to examine the local structure of the reaction zone for a moderately fast reaction between unmixed species in decaying, homogeneous turbulence and in a homogeneous turbulent shear flow. Pseudospectral techniques were used in domains of 64 exp 3 and higher wavenumbers. A finite-rate, single step reaction between non-premixed reactants was considered, and in one case temperature-dependent Arrhenius kinetics was assumed. Locally intense reaction rates that tend to persist throughout the simulations occur in locations where the reactant concentration gradients are large and are amplified by the local rate of strain. The reaction zones are more organized in the case of a uniform mean shear than in isotropic turbulence, and regions of intense reaction rate appear to be associated with vortex structures such as horseshoe vortices and fingers seen in mixing layers. Concentration gradients tend to align with the direction of the most compressive principal strain rate, more so in the isotropic case.

Description: The objective of this work is to develop a space-time accurate numerical method for the solution of incompressible Navier-Stokes equations in generalized coordinates. The resulting code is to be used for direct and large-eddy simulation of turbulence in complex geometries. In a previous paper, the system of Navier-Stokes equations in general curvilinear coordinates was solved by a second-order accurate finite-difference scheme. Satisfactory results were obtained for several flows in two and three dimensions. The system of Navier-Stokes for the fluxes are given in Orlandi (1989). The main deficiency of the numerical scheme was the large CPU time required for the solution of the Poisson equation for the 'pressure' field. The point SOR relaxation, in conjunction with a multigrid scheme, was used for the Poisson equation. In some cases, particularly with very fine grids, it was impossible to obtain a divergent-free flow. A preliminary attempt is made to compute the spatially evolving flow of Swearingen & Blackwelder. To reduce the streamwise distance, the inflow was at a distance x = 60 cm from the leading edge.

Description: Vortex filaments in superfluids such as helium 2 may provide new insights into very high Reynolds number flows. The behavior of a superfluid vortex ring interacting with a normal fluid shear flow, specifically channel flow, is simulated. The vortex ring evolves into a stable horseshoe configuration which propagates without further change of form. In this simulation, a boundary layer behavior in a superfluid through the coupling of the superfluid and the normal fluid is demonstrated.

Description: Explicit solutions of the stationary Hopf equation are discussed and their computational possibilities are explored. The motivation is to circumvent the infinite hierarchy of coupled equations for the velocity moments and obtain an exact closure of the steady-state 3D Navier-Stokes equations, without modeling assumptions or truncation. The Hopf formulation of the Navier-Stokes equation is reviewed. A stationary homogeneous solution for 2D flow is displayed and discussed. It is shown how depletion of nonlinearity may arise for 3D forced homogeneous flow. The general 3D forced case is considered and a method for closing the 3D unforced equations with arbitrary boundary conditions is derived.

Description: The solar nebula, from which the planets in our solar system formed, featured a disk of gas and dust grains in rapid, differential rotation, and at some stage was likely to have been unstable to thermal convection. This situation is suspected by many to lead to significant turbulent Reynolds stress production and angular momentum transport in such systems, and estimates of transport rates have been attempted from unsubstantiated phenomenological models. In order to determine the circumstances and physical conditions under which our own planetary system formed and to explain recent observations of young stellar systems, it is necessary to develop realistic models of heat and angular momentum transport for such flows. Developing an understanding of complicated flows featuring thermal convection, rotation, and shear is also of wide interest in stellar astrophysics and in planetary and terrestrial atmospheric studies. The ultimate objective is to develop workable models based on the numerical simulations for constructing global solar nebula models; viz., relatively simple prescriptions for heat and angular momentum fluxes from given system parameters (e.g., ratios of rotation, shear, and convective lapse rates) are characterized, quantified, and developed. Toward this end, our program has been to attempt to understand the behavior of the direct numerical simulations of Boussinesq convection, which, despite the complexity of the results, is still an overly simplified approximation to the real system and should be more amenable to analysis. These results are also intended to be tested against turbulence models, especially those designed for atmospheric boundary layers, and may provide a basis for subgrid-scale models. In order to make the numerical simulations more realistic with regard to the solar nebula problem, a fully compressible code that will allow incorporation of large density stratifications and realistic thermodynamic and radiative properties is developed. In order to explore the properties of these flows at the very high values of Re found in natural systems and the very low values of Pr found in most astrophysical contexts, we will need to employ large-eddy simulations for which we want to determine the most appropriate subgrid-scale model to incorporate.

Description: Turbulent plane Couette flow was numerically simulated at a Reynolds number (U(sub w)h/nu) of 6000, where U(sub w) is the relative wall speed and h is half the channel-height. Unlike in Poiseuille flow, where the mean shear rate changes its sign at the centerline, the sign of mean shear rate in plane Couette flow remains the same across the whole channel. This difference is expected to yield several differences between the two flows, especially in the core region. The most significant and dramatic difference observed was the existence of large-scale structures in the core region of the plane Couette flow. The large eddies are extremely long in the flow direction and fill the entire channel (i.e., their vertical extent is 2h). The large-scale structures have the largest contribution from the wavenumber (k(sub x)h,k(sub z)h) = (0, plus or minus 1.5), corresponding to a wavelength lambda(sub z)/h is approximately equal to 4. The secondary motion associated with the k(sub x)h = 0 mode consists of the large-scale vortices. The large eddies contribute about 30 percent of turbulent kinetic energy.

Description: The National Aerospace Plane (NASP) will require thermal insulation systems which are consistent with cryogenic fluids, high thermal loads, and design restrictions such as weight and volume. Test sections of the proposed system have been constructed and evaluated. In this paper we discuss the components of the insulation system, the application of the insulation system to the NASP liquid hydrogen fuel tank system, and thermal conductivity measurements performed on test sections of the system. Both steady-state and transient thermal measurements are presented.

Description: The surface distribution of star clusters in the LMC has revealed the existence of two elliptical systems superimposed. A small surface density outer ellipse outlining the large system and a higher density elliptical inner system. The major axis of the two subsystems form an angle of 60 deg. The central subsystem contains all youngest populous globular star clusters from both the stellar and dynamical point of view. The large outer system seems to contain all old globular clusters, showing a dumpy distribution and outlining two arms at the northeast and southwest. The small LMC clusters were found to have masses less than 1000 solar masses and are distributed all over the large elliptical system.

Description: On the basis of optical and IUE observations we have investigated the possibility that some of the very young stellar systems of the LMC, such as SL360, may be regarded as globular clusters at the very early phases of their dynamical evolution after the gas cloud expulsion.

Description: A tentative detection of the J = 1 - 0 emission line of (C-13)O has been obtained with SEST from a 24.4 hour integration. The velocity resolution used was 0.23 km/s and the FWHP beamwidth was 45 arcsec. If the (C-13)O line data are conservatively interpreted as an upper limit, the (C-12)O/(C-13)O ratio is not less than 60. Our result supports the previous determination of a large value of the isotope ratio in this cloud, made using radio emission lines with a 1.6-arcmin beam, and extends the ratio based on emission lines to a smaller region. When interpreted as a lower limit, our data is consistent with the ratio obtained from UV absorption line data for (C-12)O and (C-13)O.

Description: We present measurements of the integrated radio continuum flux density of M33 at frequencies between 22 and 610 MHz and discuss the radio continuum spectrum of M33 between 22 MHz and 10 GHz. This spectrum has a turnover between 500 and 900 MHz, depending on the steepness of the high frequency radio spectrum of M33. Below 500 MHz the spectrum is relatively flat. We discuss possible mechanisms to explain this spectral shape and consider efficient free-free absorption of nonthermal emission by a cool (not greater than 1000 K) ionized gas to be a very likely possibility. The surface filling factor of both the nonthermal and the thermal material appears to be small (of order 0.001), which could be explained by magnetic field/density fluctuations in the M 33 interstellar medium. We briefly speculate on the possible presence of a nuclear radio source with a steep spectrum.

Description: We present aperture synthesis observations of H I (21 cm) line radiation and continuum emission at 408 and 1420 MHz towards a field centered on the molecular cloud B5. The H I emission shows an extended atomic halo around the molecular cloud. The opacity of the halo is derived using H I absorption toward several background sources and a simple source model is presented. The model indicates that the halo is not gravitationally bound to the molecular cloud and that it is in fact expanding away from it. Approximately 350 solar masses are contained in the H I halo. Flux densities and spectral indices for the sources detected in both of the continuum bands are given.

Description: The long-period variable Mira Omicron Ceti has been observed at 800 nm wavelength with the Mk III Optical Interferometer at photometric phases 0 = 0.96, 0.05, and 0.14 in 1990; some additional data were taken in 1989. The star is not spherically symmetric, and temporal variations of the size and the position angle of the asymmetry are detected. The visibility data can be represented by two uniform elliptical disks, which could correspond to different layers in Mira's atmosphere. If the size variations are interpreted in terms of physical motion of the emitting material, an infall velocity of about 15 km/s around phi = 0 is derived.

Description: The 1(10)-1(11) transition of ortho-H2D(+) at 372 GHz has been sought in several dark clouds. The transition was not detected; the best upper limits obtained are about 0.3 K (3 sigma). We derive upper limits for the ortho-H2D(+) column density and briefly discuss their meaning in comparison with a simple chemical model we have developed (Pagani et al., 1992).

Description: We present additional photometric observations (Stromgren y filter) of the Herbig Ae Star HR 5999. This new set of data, composed of 282 data points covering April 1983 to August 1989, was used in conjunction with the 362 data points compiled by Baade and Stahl (1989). Our aim was to detect single or multiple periods in the integrated set of data spanning nearly 20 yr. Through the use of series analysis techniques, and by removing any linear component present in the raw data, we were able to detect weak peaks in the power spectra (in the order of intensity the strongest peaks are at 301 and 113 d) which were not self-evident in the phase diagrams. The overall picture of the data shows a small linear component, which is more noticeable in the new set of data, indicating that the recent maxima appear brighter. By a detailed analysis of the full light curve we were able to detect the presence of well-defined 'pulses' or bursts. Further Gaussian fits of these bursts indicated that the pulses are relatively rapid, of the order of 10 d, with some of them being closely spaced. No periodicity was found for the pulses. A possible explanation of these aperiodic outbursts is that they detect a flow of matter accompanied by magnetic field disturbances originated from the interior of the star, unlike those originating externally such as binary perturbations, which tend to be periodic.

Description: On May 3, 1991, the Energetic Gamma Ray Experiment Telescope on the Compton Gamma Ray Observatory detected a gamma-ray burst both in the energy measurement subsystem and independently in the spark chamber assembly. Six individual photons were detected in the spark chamber, allowing a determination of the burst arrival direction which was l(II) = 171.9 deg +/- 1.3 deg, b(II) = 5.3 deg +/- 1.1 deg. Three energy spectra were measured from 1 to 200 MeV; they were measured during the first second after the Burst and Transient Sources Experiment trigger, the next two seconds, and the subsequent four seconds. The first two spectra exhibit a similar differential spectra index of about -2.2 with no apparent high-energy cut-off. By the time of the third spectrum, an additional soft component is evident.

Description: The effects of convection on diffusive-convective physical vapor transport process are examined computationally. We analyze conditions ranging from typical laboratory conditions to conditions achievable only in a low gravity environment. This corresponds to thermal Rayleigh numbers Ra(T) ranging from 1.80 to 1.92 x 10 exp 6. Our results indicate that the effect of the sublimation and condensation fluxes at the boundaries is 10 increase the threshold of instability. For typical ground based conditions time dependent oscillatory convection can occur. This results in nonuniform temperature and concentration gradients at the crystal interface. Spectral analysis of the flow field shows regions of both periodic and quasi-periodic states. Low gravity conditions can effectively reduce convective effects, thus resulting in uniform temperature and concentration gradients at the interface, a desirable condition for crystal growth.

Description: A model is developed for describing the interaction of vortex-drop clusters in a flowing gaseous jet, convecting downstream from an injection location. Results are presented for a stationary case representing the situation when identical clusters are continuously injected and the injection rate is constant. The results indicate that, in a rich mixture high-drop-number density regime, the mass evaporated from the drops controls the velocity of the cluster-in-vortex as it propagates downstream.

Description: An algorithm is presented for unsteady two-dimensional incompressible Navier-Stokes calculations. This algorithm is based on the fourth order partial differential equation for incompressible fluid flow which uses the streamfunction as the only dependent variable. The algorithm is second order accurate in both time and space. It uses a multigrid solver at each time step. It is extremely efficient with respect to the use of both CPU time and physical memory. It is extremely robust with respect to Reynolds number.

Description: A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.

Description: The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

Description: An implicit numerical algorithm for the time accurate solution of the compressible Navier-Stokes equations is described. Results for steady flow past a finite flat plate are presented, together with preliminary results for the temporal simulation of second mode instability in a flat plate boundary layer at Mach 4.5.

Description: The main objective of this research is to address two important but unresolved problems: (1) the measurement of vertical and transverse length scales via space correlations for all Reynolds stress components and velocity-temperature correlations, both in the free stream and within the boundary layer using the existing triple and quad-wire probes; and (2) to relate the character of the free stream turbulence to the character of the turbulence within the boundary layer in order to determine the effect on surface heat transfer.

Description: An experimental research project aimed at obtaining quantitative data on the behavior of the secondary vortex structure in a turbulent mixing layer at moderate Reynolds numbers (Re(sub delta) = 2.9 x 10(exp 4)) is discussed. This project was terminated before all the contemplated measurements could be made, and data were obtained only on the spatially stationary part of the secondary structure. Nonetheless, these data reveal some interesting facets of mixing layer behavior which are discussed.

Description: The wavelet-transformed Navier-Stokes equations are used to define quantities such as the transfer of kinetic energy and the flux of kinetic energy through scale r at position x. Direct numerical simulations of turbulent shear flow reveal that although their mean spatial values agree with their traditional counterparts in Fourier space, their spatial variability at every scale is very large, exhibiting non-Gaussian statistics. The local flux of energy involving scales smaller than some r also exhibits large spatial intermittency, and it is negative quite often, indicative of local inverse cascades.

Description: Large-scale structures in turbulent and transitional wall-bounded flows make a significant contribution to the Reynolds stress and turbulent energy. The behavior of these structures is examined. Small perturbations are introduced into a laminar and a turbulent boundary layer to trigger the formation of large-scale features. Both flows use the same inlet unit Reynolds number, and they experience the same pressure gradient history, i.e. a favorable pressure gradient (FPG) followed by an adverse pressure gradient (APG). The perturbation consists of a small short duration flow repetitively introduced through a hole in the wall located at the C(sub p) minimum. Hot-wire data are averaged on the basis of the phase of the disturbance, and automation of the experiment was used to obtain measurements on large spatially dense grids. In the turbulent boundary, the perturbation evolves into a vortex loop which retains its identity for a considerable streamwise distance. In the laminar layer, the perturbation decays to a very small magnitude before growing rapidly and triggering the transition process in the APG. The 'time-like' animations of the phase-averaged data are used to gain insight into the naturally occurring physical mechanisms in each flow.

Description: The dependence of product generation on Peclet and Reynolds numbers in a numerically simulated, reacting, two dimensional, temporally growing mixing layer is related theoretically to the fractal dimension of the passive scalar interfaces. This relation is verified using product generation measurements and dimensions derived from a standard box counting technique. A transition from a low initial dimension to a higher one of approximately 5/3 is identified and shown to be associated to the kinematic distortion on the flow field during the first pairing interaction. It is suggested that the structures responsible for this transition are non-deterministic, non-random, inhomogeneous fractals. Only the large scales are involved. No further transitions, either in the spectra of the vorticity field or in the mixing behavior, are found for Reynolds numbers up to 90,000.

Description: Studies of supersonic mixing which were accomplished over the last year with Center for Turbulence Research (CTR) support are described. During this period, a Nd:YAG laser, optical components, and data acquisition computer were obtained. This allowed detailed visualizations of the flow structure to be performed at a rapid rate, representing a significant improvement over our previous attempts. Aspects of the flow structure are described below. In addition, preliminary findings on a possible mixing enhancement strategy are also shown using the flow visualization technique.

Description: A formulation has been developed to describe the evaporation of dense or dilute clusters of binary-fuel drops. The binary fuel is assumed to be made of a solute and a solvent whose volatility is much lower than that of the solute. Convective flow effects, inducing a circulatory motion inside the drops, are taken into account, as well as turbulence external to the cluster volume. Results obtained with this model show that, similar to the conclusions for single isolated drops, the evaporation of the volatile is controlled by liquid mass diffusion when the cluster is dilute. In contrast, when the cluster is dense, the evaporation of the volatile is controlled by surface layer stripping, that is, by the regression rate of the drop, which is in fact controlled by the evaporation rate of the solvent. These conclusions are in agreement with existing experimental observations. Parametric studies show that these conclusions remain valid with changes in ambient temperature, initial slip velocity between drops and gas, initial drop size, initial cluster size, initial liquid mass fraction of the solute, and various combinations of solvent and solute. The implications of these results for computationally intensive combustor calculations are discussed.

Description: It is found through two-dimensional temporal simulations of high-speed free shear layers that mean flow distortion is significantly increased when supersonic disturbances are introduced as initial conditions. The shear layer exhibits no subharmonic growth or roll-up, but rather a spectral broadening as energy is distributed into higher harmonics. Increasing the velocity of one side of the mixing layer (u2) to roughly 1/5 the speed of the high speed side (u1), allows a slight subharmonic growth at a very slow rate for two-dimensional modes. A first look at three-dimensional free shear flows is also presented for M = 2. No effect is seen for incompressible flow; however, stabilization is seen with respect to maximum temporal growth rates as the transverse velocity is increased. A much stronger, but similar effect is seen if u2 is increased. The wave direction of maximum growth for u2 is greater than 0.4 is found to be the direction of the faster stream (u1) over a broad range of transverse velocities.

Description: An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's symmetric total variation diminishing scheme. Newton relaxation of the fully coupled equation set is employed on a cell-to-cell basis. Under-relaxation of the inviscid and over-relaxation of the viscous contributions to the residual are implemented. Computational work is easily partitioned among many processors in an asynchronous, dynamic mode for convergence acceleration. An overview of the physical models employed herein for thermochemical nonequilibrium is included. Several test cases and comparisons with experimental data are presented involving hypersonic flow over blunt bodies which illustrate the qualitative and quantitative capabilities of this approach.

Description: The outburst of Halley comet that occurred in February 1991 when the comet was over 14 AU from the sun, between the orbits of Saturn and Uranus, is described. This region of the solar system was known to be so cold that a comet nucleus would be totally inert there. On February 12th, 1991 Halley brightened by a factor of about 300, and a dust coma was observed around the nucleus. Different explanations of the comet's strange behavior are suggested.

Description: AGN spectrum and spectral features, polarization, inclination, and X-ray line and continuum reflection features are discussed in a critical way in order to determine the ones that are the least model-dependent. The sign and strength of absorption and emission edges are found to be model-dependent, and relativistic broadening and shifting makes them hard to detect. The presence or absence of the predicted Lyman edge polarization feature may be used as a decisive test for thin, bare AGN disks. Other good model-independent tests are several inclination-related line and continuum correlations in big AGN samples. It is shown that electron temperature near the surface of the disk can greatly exceed the disk equilibrium temperature, which causes deviations from LTE. This effect must be incorporated into realistic disk models.

Description: Active galaxies are believed to consist of a compact nucleus, the standard model for which is a massive black hole or a cluster of black holes. A different paradigm is considered here, deriving from quark confinement theory in QCD. It is an 'astrophysical bag', modelled after the 'hadron bags' of particle physics which have already been studied in astrophysics as quark stars. Another interpretation of the cosmological constant in general relativity, and possibly a new quasar redshift formula, are introduced. As a highly-energetic object, this model may resolve the baryonic matter problem for fuelling AGN accretion processes which black hole paradigms cannot account for. Here, baryons, cosmic rays, and neutrinos are free.

Description: Although there are some flaws in each of the methods to constrain the black hole mass, which are described in the above paper, all of them unambiguously indicate that the mass of a putative black hole at the Galactic center is much smaller than a widely quoted value, about 10 exp 6 M(solar). Possibly, the mass of the black hole, if any, does not exceed 100 M(solar) or so. Further work is needed in order to distinguish between such a modest-mass black hole and its potential rival such as a very massive star of about the same mass. Anyway, the Galactic nucleus seems to be a tracer of a past star burst rather than a remnant of a Seyfert galaxy activity.

Description: Attention is given to an empirical model for transition to turbulence in oscillatory flows in straight tubes. Designed after a correlation for transition of a boundary layer on a flat plate, the model yields the laminar flow momentum thickness Reynolds number that must be met before transition to turbulence will occur. The transition point is located by comparing this to the actual momentum thickness Reynolds number. A scheme is proposed for estimating the momentum thickness Reynolds number in terms of the position within the cycle, the maximum value of the diameter Reynolds within the cycle, Re(max), and the dimensionless frequency, Valensi number. Results from an experimental study of oscillatory flow in a tube are employed to develop the model. When the flow is determined to be turbulent, it is proposed that a fully-developed, steady flow friction coefficient be applied. When the flow is laminar, the assumption of fully developed flow cannot be made; thus, a method is suggested for estimating the friction factor.

Description: The effects of rotation on grid turbulence were studied in a low speed tunnel. Results are presented for both axisymmetric and non-axisymmetric flows. As observed by previous workers, the cascade process was seen to be effectively blocked by rotation leading to reduced dissipation. In the non-axisymmetric case, rotation caused the flow to tend rapidly towards axisymmetry; however, for the length studied, the flow tended to move away from isotropy. Based on the results, some recommendations are made on the design of a larger facility for studying similar flows.

Description: It is shown that the ideas of selective amplification and direct resonance, based on linear theory, can not provide an explanation for the well-defined streak spacing of about 100 wall units (referred to as 100(+) hereafter) in wall-bounded turbulent shear flows. In addition, for the direct resonance theory, the streaks are created by the non-linear self-interaction of the vertical velocity rather than of the directly forced vertical vorticity. In view of the failure of these approaches, it is then proposed that the selection mechanism must be inherently non-linear and correspond to a self-sustaining mechanism. The 100(+) value should thus be considered as a critical Reynolds number for that mechanism. Indeed, in the case of Poiseuille flow, this 100(+) criterion for transition to turbulence corresponds to the usually quoted value of 1000 based on the half-width and the centerline velocity. In Couette flow, it corresponds to a critical Reynolds number of about 400 based on the half width and half velocity difference.

Description: The development of an automated technique for the eduction of 3-D spatial patterns in vector or scalar diagnostics was completed. The method is based on an iterative convolution between a trial pattern and the data field. It was applied to the analysis of low Reynolds number turbulent channel flow and homogeneous shear flow. The results yielded new information on the dominant flow structures in these flows, particularly with respect to the spatial relationships between various forms of organized motion. A particular application of the pattern eduction method, which is tentatively referred to as an 'adaptive wavelet transformation', is proposed with the objective of investigating the way turbulence structure changes with scale. Preliminary results using data from homogeneous turbulent shear flow simulations are presented. At the low Reynolds numbers of the simulations, there is no evidence of scale similarity. The small scales appear to be associated with the edges of the larger scale vortical structures.

Description: The existence and importance of large-scale spanwise vortical structures for 2-D straight mixing layers has been well documented in the last decade. Computer models and simulations have sought to reproduce these vortical structures associated with the Kelvin-Helmholtz (K-H) instability mode which is due to the shear per se. Secondary streamwise vortical structures for the same flows were also seen experimentally and have recently been given importance in computational efforts. Curved mixing layers can be characterized as stable (the high-speed stream is placed on the outside of the longitudinal bend), leading to a suppression of the Taylor-Gortler (T-G) instability, and unstable (high-speed stream on the inside of the bend), leading to an enhancement of the T-G instability. The T-G instability is associated with the centripetal acceleration that the curvature imparts. Thus, curvature superimposed on 2-D shear layer flows provides a way for studying the importance of streamwise vorticity, its competition with spanwise vorticity, and changes to entrainment and mixing. Furthermore, the outcome of the competition of a relatively enhanced or suppressed T-G instability with the K-H instability offers the possibility of achieving passive mixing enhancement. As a first step in understanding the competition between the K-H and the T-G instabilities and the resulting changes to the structure of the flow, highly resolved visualizations of the flow structure for the stable and the unstable configurations are provided. The straight layer is also visualized for comparison with earlier works.

Description: Computer flow simulation aided by dynamical systems analysis is used to investigate the kinematics of time-periodic vortex shedding past a two-dimensional circular cylinder in the context of the following general questions: (1) Is a dynamical systems viewpoint useful in the understanding of this and similar problems involving time-periodic shedding behind bluff bodies; and (2) Is it indeed possible, by adopting such a point of view, to complement previous analyses or to understand kinematical aspects of the vortex shedding process that somehow remained hidden in previous approaches. We argue that the answers to these questions are positive. Results are described.

Description: An ongoing research effort designed to apply the best possible second-moment-closure model to simulate complex hypersonic flows is presented. The baseline model under consideration is the Launder-Reece-Rodi Reynolds stress transport turbulence model. Two add-ons accounting for wall effects, namely, the Launder-Shima low-Reynolds-number model and the compressible wall-function technique, are tested. Results are reported for flow over a flat plate, both adiabatic-wall and cooled-wall cases. It has been found that further improvements of the existing models are necessary to achieve accurate prediction in high Mach number flow range.

Description: Scalar fields undergoing random advection have attracted much attention from researchers in both the theoretical and practical sectors. Research interest spans from the study of the small scale structures of turbulent scalar fields to the modeling and simulations of turbulent reacting flows. The probability density function (PDF) method is an effective tool in the study of turbulent scalar fields, especially for those which involve chemical reactions. It has been argued that a one-point, joint PDF approach is the one to choose from among many simulation and closure methods for turbulent combustion and chemically reacting flows based on its practical feasibility in the foreseeable future for multiple reactants. Instead of the multi-point PDF, the joint PDF of a scalar and its gradient which represents the roles of both scalar and scalar diffusion is introduced. A proper closure model for the molecular diffusion term in the PDF equation is investigated. Another direction in this research is to study the mapping closure method that has been recently proposed to deal with the PDF's in turbulent fields. This method seems to have captured the physics correctly when applied to diffusion problems. However, if the turbulent stretching is included, the amplitude mapping has to be supplemented by either adjusting the parameters representing turbulent stretching at each time step or by introducing the coordinate mapping. This technique is still under development and seems to be quite promising. The final objective of this project is to understand some fundamental properties of the turbulent scalar fields and to develop practical numerical schemes that are capable of handling turbulent reacting flows.

Description: The immediate goal is to understand and validate the Yakhot-Orszag model of the velocity-derivative skewness and model equation for the rate of energy dissipation epsilon. A summary of a more detailed manuscript in preparation is presented. The purpose is to clarify some limitations of the theory by careful examination of key assumptions and approximations, and thereby to encourage its improvement.

Description: Investigations carried out for revisiting homogeneous turbulent flows in the presence of mean shear, rotation, or external compression are summarized. The simplest and most concise RDT (Rapid Distortion Theory) formulation, which includes a comprehensive linear stability analysis, is used for this purpose. Such a linear approach could be extended by a generalized EDQNM (Eddy Damped Quasi-Normal Markovian) to two point closure in order to model non-linear interactions, especially when pure Coriolis effects are present. The results are discussed in connection with databases obtained by DNS (Direct Numerical Simulations), including previous Center for Turbulence Research (CTR) results and new calculations in progress. The main goal is to contribute to and significantly improve on the rational one-point closure models in progress at the CTR and at Ecole Centrale de Lyon (ECL).

Description: The near-wall region plays an essential role in turbulent boundary layers: it is a region of high shear; the peak rate of production and peak intensity of turbulence occurs there; and the peak rate of dissipation occurs right at the wall. Nevertheless, this region has received less attention from modelers than have more nearly homogeneous flows. One reason for this is that when the boundary layer is near equilibrium, experimental data can be used to prescribe the flow in the wall layer. Another reason is that most turbulence models are developed under assumptions of near homogeneity. This is a poor approximation in the wall region. A single-point moment closure model for the strongly non-homogeneous A turbulent flow near a rigid boundary is developed.

Description: This paper presents a 2D axisymmetric combined conduction and radiation model of a multizone crystal growth furnace. The model is based on a programmable multizone furnace (PMZF) designed and built at NASA Lewis Research Center for growing high quality semiconductor crystals. A novel feature of this model is a control algorithm which automatically adjusts the power in any number of independently controlled heaters to establish the desired crystal temperatures in the furnace model. The control algorithm eliminates the need for numerous trial and error runs previously required to obtain the same results. The finite element code, FIDAP, used to develop the furnace model, was modified to directly incorporate the control algorithm. This algorithm, which presently uses PID control, and the associated heat transfer model are briefly discussed. Together, they have been used to predict the heater power distributions for a variety of furnace configurations and desired temperature profiles. Examples are included to demonstrate the effectiveness of the PID controlled model in establishing isothermal, Bridgman, and other complicated temperature profies in the sample. Finally, an example is given to show how the algorithm can be used to change the desired profile with time according to a prescribed temperature-time evolution.

Description: The cooling arrangement of the Space Shuttle Main Engine High Pressure Oxidizer Turbopump (HPOTP) incorporates two jet rings, each of which produces 19 high-velocity coolant jets. At some operating conditions, the frequency of excitation associated with the 19 jets coincides with the natural frequency of the turbine blades, contributing to fatigue cracking of blade shanks. In this paper, an alternate turbine disk cooling arrangement, applicable to disk faces of zero hub radius, is evaluated, which consists of a single coolant jet impinging at the center of the turbine disk. Results of the CFD analysis show that replacing the jet ring with a single central coolant jet in the HPOTP leads to an acceptable thermal environment at the disk rim. Based on the predictions of flow and temperature fields for operating conditions, the single central jet cooling system was recommended for implementation into the development program of the Technology Test Bed Engine at NASA Marshall Space Flight Center.

Description: Horizontal-branch (HB) sequences have been computed for different rates of convective overshooting at the edge of the convective core. We find that the convective core undergoes a series of 'breathing pulses' during most of the HB phase. The average behavior of these sequences closely mimics the behavior of sequences computed with canonical semiconvection.

Description: The interplanetary magnetic field (IMF) behaves in a reasonably well-understood manner between the Sun and the heliospheric termination shock. At the shock, the azimuthal field is amplified by a factor of four (for a strong shock) and undergoes secular amplification in the heliosheath until the flow is fully turned into the downstream direction and has reached its asymptotic state in the distant heliotail. This amplification may lead to important MHD effects that can cause the shock to be closer to the Sun than otherwise expected. Here we further examine whether there are important MHD effects in the heliosheath. We do this by calculating the kinematic compression of the magnetic field in the heliosheath using an analytic incompressible flow model of the dynamics downstream of the shock. We conclude that it is likely that MHD effects are important in the heliosheath in a narrow cone about the upstream direction.

Description: The suitability of the acetate replication method for monitoring the growth of small cracks is discussed. Applications of this technique are shown for cracks growing at the notch root in semicircular-edge-notch specimens of a variety of aluminum alloys and one steel. The calculated crack growth rate versus Delta K relationship for small cracks was compared to that for large cracks obtained from middle-crack-tension specimens. The primary advantage of this techinque is that it provides an opportunity, at the completion of the test, to go backward in time towards the crack initiation event and 'zoom in' on areas of interest on the specimen surface with a resolution of about 0.1 micron. The primary disadvantage is the inability to automate the process. Also, for some materials, the replication process may alter the crack-tip chemistry or plastic zone, thereby affecting crack growth rates.

Description: We analyze the dynamics of a rotating cloudy gaseous medium relative to 2D perturbations having wavelengths shorter than the disk thickness. The viscosity due to cloud-cloud collisions is confronted with the derived critical viscosity, below which a collective mode of perturbations in the disk is able to grow until going into the saturation regime where fully developed turbulence is established. It is shown that for the assumed parameters of AGN accretion disks consistent with the observational data, the viscosity due to collective modes usually dominates over that due to cloud-cloud collisions.

Description: Archival EXOSAT and HEAO1-A2 data from Cyg X-1 show the 'high energy excess' above 10 keV seen in X-ray observations of AGN. Using a likelihood ratio test, we are for the first time able to distinguish conclusively in favor of Compton reflection rather than partial covering as the origin of the high energy excess. This supports the idea of an X-ray illuminated accretion disk in Cyg X-1, but the line equivalent width is smaller by a factor of 2-3 than that expected from such a disk. While the larger optical depth required for reflection as opposed to line emission admit the possibility of seeing line without reflection, the converse is not possible. To see a reflection spectrum, including the strong iron absorption edge, implies that strong iron emission must be observed as the line and edge are causally linked.

Description: A review is given of methods to constrain the mass of a putative black hole at the center of the Milky Way galaxy. These approaches are based on: (1) tidal disruption of stars by a massive black hole; (2) displacement between Sgr A* and either of the IRS 16 components; (3) electron-positron pair production by a black hole, and (4) mass outflow centered on IRS 16/Sgr A* complex.

Description: The origin of the cosmic X-ray and gamma-ray backgrounds is explained via the mechanism of AGN spectral-luminosity evolution. The spectral evolution of precursor active galaxies into AGN, and Newton-Raphson input and output parameters are discussed.

Description: We present results of monitoring a sample of five OVV qausars in an effort to determine if these types of AGNs exhibit microvariability similar to that detected for BL Lacertae objects. Three of the five quasars investigated clearly exhibited microvariability. The timescales for the variations for both BL Lacertae objects and OVV quasars are similar which suggests that the radiation is generated in similar spatial volumes. These variations are consistent with the results of models which explain the observed microvariability as the result of excess emission produced by flares or hot spots randomly appearing and disappearing on the accretion disks around supermassive black holes.