ALBERT

Description: As turbine-engine core operating conditions become more severe, designers must develop more effective means of cooling blades and vanes. In order to design reliable, cooled turbine blades, advanced transient thermal calculation techniques are required. The TACT1 computer program was developed to perform transient and steady-state heat-transfer and coolant-flow analyses for cooled blades, given the outside hot-gas boundary condition, the coolant inlet conditions, the geometry of the blade shell, and the cooling configuration. TACT1 can analyze turbine blades, or vanes, equipped with a central coolant-plenum insert from which coolant-air impinges on the inner surface of the blade shell. Coolant-side heat-transfer coefficients are calculated with the heat transfer mode at each station being user specified as either impingement with crossflow, forced convection channel flow, or forced convection over pin fins. A limited capability to handle film cooling is also available in the program. The TACT1 program solves for the blade temperature distribution using a transient energy equation for each node. The nodal energy balances are linearized, one-dimensional, heat-conduction equations which are applied at the wall-outer-surface node, at the junction of the cladding and the metal node, and at the wall-inner-surface node. At the mid-metal node a linear, three-dimensional, heat-conduction equation is used. Similarly, the coolant pressure distribution is determined by solving the set of transfer momentum equations for the one-dimensional flow between adjacent fluid nodes. In the coolant channel, energy and momentum equations for one-dimensional compressible flow, including friction and heat transfer, are used for the elemental channel length between two coolant nodes. The TACT1 program first obtains a steady-state solution using iterative calculations to obtain convergence of stable temperatures, pressures, coolant-flow split, and overall coolant mass balance. Transient calculations are based on the steady-state solutions obtained. Input to the TACT1 program includes a geometrical description of the blade and insert, the nodal spacing to be used, and the boundary conditions describing the outside hot-gas and the coolant-inlet conditions. The program output includes the value of nodal temperatures and pressures at each iteration. The final solution output includes the temperature at each coolant node, and the coolant flow rates and Reynolds numbers. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 480K of 8 bit bytes. The TACT1 program was developed in 1978.

Description: This is a finite-difference program for calculating the viscous compressible boundary layer flow over either planar or axisymmetric surfaces. The flow may be initially laminar and progress through a transitional zone to a fully turbulent flow, or it may remain laminar, depending on the imposed boundary conditions, laws of viscosity, and numerical solution of the momentum and energy equations. The flow may also be forced into a turbulent flow at a chosen spot by the data input. The input may contain factors of arbitrary Reynolds number, free-stream Mach number, free stream turbulence, wall heating or cooling, longitudinal wall curvature, wall suction or blowing, and wall roughness. The solution may start from an initial Falkner-Skan similarity profile, an approximate equilibrium turbulent profile, or an initial arbitrary input profile. This program has been implemented on the IBM 7094/7044 Direct Couple System. This program is written in FORTRAN IV and was developed in 1974.

Description: A computer program has been developed for the design of sharp-edged throat supersonic nozzles where losses are accounted for by correcting the ideal nozzle geometry for boundary layer displacement thickness. The ideal nozzle is designed by the method of characteristics to produce uniform parallel flow at the nozzle exit in the smallest possible distance. Boundary-layer parameters (displacement and momentum thicknesses) are calculated for the ideal nozzle, and the final nozzle geometry is obtained by adding the displacement thickness to the ideal nozzle coordinates. The boundary layer parameters are also used to calculate the aftermixing conditions downstream of the nozzle assuming the flow mixes to a uniform state. The computer program input consists essentially of the nozzle-exit Mach number, specific-heat ratio, nozzle angle, throat half-height, nozzle subsonic section coordinates and corresponding pressure ratios, total temperature and pressure, gas constant, and initial momentum or displacement thickness. The program gas properties are set up for air; for other gases, changes are required to the program. The computer program output consists of the corrected nozzle coordinates, the principal boundary-layer parameters, and the aftermixing conditions. This program has been implemented on the IBM 7094/7044 Direct Couple System.

Description: A computer program has been developed which analyzes by means of mathematical models the temperature profiles in the contents of a filled propellant tank. In designing space vehicles using cryogenic liquid propellants, it is necessary to know how heat transferred from the tank walls and heat absorbed internally affect the temperature distribution with the tank contents. The mathematical flow model is based on results from small-scale experiments. The results showed that when a subcooled fluid is subject to both nonuniform internal heating and wall heating, two distinct temperature regions are developed. In the lower region, the fluid is thoroughly mixed and maintains a uniform temperature profile. In the upper region, a stratified layer develops, and a temperature gradient is formed from the accumulation of warm fluid from the boundary layer along the tank walls; it also indicated that the temperature profiles in the stratified layer exhibited similarity. This concept was developed primarily for internal heating caused by nuclear radiation. However, the theory and computer program are applicable for any form of internal or bulk heating. This program is written in FORTRAN IV for batch execution and has been implemented on the IBM 7094. This program was developed in 1970.

Description: This program solves the two-dimensional, compressible laminar or turbulent boundary-layer equations in an arbitrary pressure gradient. Cohen and Reshotko's method is used for the laminar boundary layer, Sasman and Cresci's method for the turbulent boundary layer, and the Schlichting-Ulrich-Granville method to predict transition. Transition may also be forced at any point by the user. Separation, if it occurs, is predicted for both laminar and turbulent flow. The user may begin values for displacement thickness and momentum thickness in either laminar or turbulent flow. This program was implemented on the IBM 7094.

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: Newell and Meng (1992) present maps of the occurrence probability of various classifications of particle precipitation as seen in the dayside topside ionosphere. It is argued that these are maps of the magnetospheric regions, a contention with which their critics disagree. The latter conclude that, because of convection, any one population of particles seen at low altitudes will have originated from a wide variety of locations, and particle characteristics cannot be mapped back to those in the magnetosphere without detailed knowledge of both the convection and magnetic field. Steplike boundaries between the regions will arise from nonsteady-state conditions and cannot be envisaged as steady-state magnetospheric boundaries between two plasma populations. In their reply Newell and Meng contend that convection does not move plasma from the LLBL into the cusp. Most of the LLBL plasma comes from the magnetosheath, so the direction of plasma transfer is in the other direction.

Description: Energy of a high harmonic standing wave tends to be distributed equally over the whole wave even in a stratified medium where the wave's peak amplitude can be much larger near the upper boundary than the lower one. This fact is generalized to the many diverse physical problems which solve second-order differential equations of Sturm-Liouville type. For any such solution y(z) whose sign fluctuates along the z-axis, quantities are found which have the same value between any two neighboring zeros of y. One of the equidistributed quantities for an oscillating fluid sphere is similar to kinetic energy but is identical only in limiting cases. The acoustic midpoint of a cavity can be a unique place where some nonlinear perturbations have extra strength. This may apply to the puzzling solar phenomenon called supergranulation.

Description: A computational code for the vorticity-potential method is developed for a three-dimensional bounded vorticity field. The evaluation of the boundary data for the vector potential in the code is improved so that the numerical solution simulates that in an unbounded domain to a high order. The time evolution of two vortex rings and that of an elliptic ring are investigated with this code. The cut-and-connect phenomena of vortex rings are successfully captured. The results are compared with those of asymptotic theory and the experiment. They also highlight the need for additional theoretical and numerical investigations.

Description: To monitor future global temperature trends, it would be extremely useful if parameters nonlinearly related to surface temperature could be found, thereby amplifying any warming signal that may exist. Evidence that global thunderstorm activity is nonlinearly related to diurnal, seasonal and interannual temperature variations is presented. Since global thunderstorm activity is also well correlated with the earth's ionospheric potential, it appears that variations of ionospheric potential, that can be measured at a single location, may be able to supply valuable information regarding global surface temperature fluctuations. The observations presented enable a prediction that a 1 percent increase in global surface temperatures may result in a 20 percent increase in ionospheric potential.

Description: The largest known effusive eruptions during the Cenozoic and Mesozoic Eras, the voluminous flood basalts, have long been suspected as being associated with major extinctions of biotic species. Despite the possible errors attached to the dates in both time series of events, the significance level of the suspected correlation is found here to be 1 percent to 4 percent. Statistically, extinctions lag eruptions by a mean time interval that is indistinguishable from zero, being much less than the average residual derived from the correlation analysis. Oceanic flood basalts, however, must have had a different biological impact, which is still uncertain owing to the small number of known examples and differing physical factors. Although not all continental flood basalts can have produced major extinction events, the noncorrelating eruptions may have led to smaller marine extinction events that terminated at least some of the less catastrophically ending geologic stages. Consequently, the 26 Myr quasi-periodicity seen in major marine extinctions may be only a sampling effect, rather than a manifestation of underlying periodicity.

Description: The monthly and daily samples of the Ap index for the interval from 1932 through 1982 were studied using the power spectrum technique. Results obtained for Bartel's period (about 27 days), the semiannual period, the dual-peak solar cycle distribution of geomagnetic storms, and certain other medium-scale periodicities are examined in detail. In addition, results on the cumulative occurrence number of storms per decade as a function of the Ap and Dst indices for the storm are presented.

Description: The circulation of polar wind ions from the high-latitude ionosphere to the plasma sheet is investigated. Single-particle trajectory simulations for the geomagnetic tail show, in addition to the breaking of adiabaticity, the possible creation of new high-altitude mirror points. This trajectory feature results from an intense parallel deceleration imparted by the magnetic field rotation during fast ExB transport. This centrifugal deceleration yields a critical parallel velocity which depends on the magnitude of the convection electric field and below which ions remain trapped inside the neutral sheet.

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: High resolution in situ measurements of reactive nitrogen (NO(y)) and O3 were made in the upper troposphere and lower stratosphere at a variety of latitudes and seasons. In the lower stratosphere, NO(y) and O3 are very highly positively correlated at all times and spatial scales sampled. The ratio NO(y)/NO3 is much less variable than either species measured separately. The ratio has a much weaker gradient with altitude than the mixing ratios of O3 or NO(y). The ratio is smaller and decreases more rapidly with altitude in the tropics than at midlatitudes. In the upper troposphere NO(y) and O3 are only weakly correlated. Their ratio in the tropical upper troposphere is about 0.005-0.025 and the ratio in the midlatitude upper troposphere is about 0.004-0.010. The NO(y) in the upper troposhere is probably partly due to lightning.

Description: The chemistry and dynamics of a convective system observed on February 2, 1987 in the EMEX and STEP campaigns are analyzed. Chemical and thermodynamic profiles in undisturbed air near the EMEX 9 system indicate that the troposphere was well mixed by previous convection. As a consequence there was little direct transport from the boundary layer to the upper troposphere and air transported upward was detrained throughout the middle and upper troposphere. There was minimal effect on O3 production. Other convective complexes located 800-900 km upstream produced greater perturbations on trace gas profile immediately below the tropopause. Ozone production was reduced by about 0.25 ppbv/d at these altitudes, representing a reduction in P(O3) of 15-20 percent over the column from 14.5 to 17 km. P(O3) from 12 to 17 km in a region distant from active convection and subject to lightning was 2-3 times higher than it would have been without lightning.

Description: Lateral variations of the temperature field in the lower mantle have been reconstructed using new results in mineral physics and seismic tomographic data. We show that, with the application of high-pressure experimental values of thermal expansivity and of sound velocities, the slow seismic anomalies in the lower mantle under the Pacific and Africa can be converted into realistic-looking plume structures with large dimensions of 0(1000 km). The outer fringes of the plumes have an excess temperature of around 400 K. In the core of the plumes are found tonguelike structures with extremely high thermal anomalies. These values can exceed 1200 K and are too high to be explained on the basis of thermal anomalies alone. We suggest that these major plumes in the deep mantle may be driven by both thermal and chemical buoyancies or that enhanced conductive heat-transfer may be important there.

Description: The largest known Cenozoic impact craters with the most accurately measured ages are found to correlate very closely with geologic stage boundaries. The level of confidence in this result is 98-99 percent even under the most pessimistic assumptions concerning dating errors. One or more large impacts may have led, in at least some cases, to the extinctions and first appearances of biotic species that mark many of the geologic stage boundaries.

Description: Ozone profiles obtained with the Differential Absorption Lidar (DIAL) system at the National Institute for Environmental Studies (NIES) (Tsukuba, Japan) were compared with data provided by the satellite sensor SAGE II. The SAGE II data were selected based on criteria of spatial and temporal differences between the DIAL and the SAGE II measurements: five degrees in latitude and 15 degrees in longitude, within a latitudinal band from 31 deg to 41 deg N, and within one, three and five days after or before the DIAL measurements. Results show very good agreement for the individual and the zonal-mean profiles. The average mean difference between the DIAL and the SAGE II measurements over the altitudes 15-50 km was about 10 percent.

Description: The BATSE on-board burst system has been triggered by over 400 terrestrial electron precipitation events. These are the single largest cause of false triggers. Bremsstrahlung is generated as the precipitating electrons interact in the Earth's atmosphere, or in the spacecraft, and this radiation is detected by the Large Area Detectors, often triggering the instrument into burst mode. Several examples of such events are presented here, and the different classes of events are described. A possible correlation of events to strong magnetospheric activity is presented, and the association of a sub-set of events to a powerful VLF transmitter on the western coast of Australia is described.

Description: The Tsyganenko models contain 'modules' for two external current systems, the tail and ring currents; the former is the more detailed. No specific modules, however, cover the magnetopause and the Birkeland current system. The effects of these currents are represented by an all-purpose 'polynomial', and the resulting formulas involve about 30 parameters that specify the model. The dependence on magnetospheric indices is also simplified.

Description: A quantitative model is developed of the magnetic field produced by the electric current system of region 2 Birkeland currents, closed via the partial ring current. The distribution of j-perpendicular is computed from a given axially asymmetric spatial distribution of hot isotropic magnetospheric plasma over an infinitely thin L shell in an axisymmetric purely dipolar geomagnetic field, while the field-aligned current density is found from the continuity of the net electric current. The magnetic field distribution is derived by a Biot-Savart integral over the electric current system. An assumed cosine dependence of the plasma pressure on local time makes it possible to reduce the problem of analytical representation of the B field to two dimensions. The obtained numerical fits for the partial ring current/region 2 Birkeland current magnetic field are relatively simple, continuous, and valid throughout the whole extraterrestrial space from ionospheric heights up to tens of earth radii.

Description: In the present grid-independent approximate Riemann solver for 2D and 3D flows that are governed by the Euler or Navier-Stokes equations, fluxes on grid faces are obtained by wave decomposition; the assumption of information-propagation in the velocity-difference directions leads to a more accurate resolution of shear and shock waves, when these are are oblique to the grid. The model, which yields significantly greater accuracy in both supersonic and subsonic first-order spatially accurate computations, describes the difference in states at each grid interface by the action of five waves.

Description: The action of weak, streamwise vortices on a plane, incompressible, steady mixing layer is examined in the large Reynolds number limit. The outer, inviscid region is bounded by a vortex sheet to which the viscous region is confined. It is shown that the local linear analysis becomes invalid at streamwise distances O(epsilon sup -1), where (epsilon much less than 1) is the crossflow amplitude, and a new nonlinear analysis is constructed for this region. Numerical solutions of the nonlinear problem show that the vortex sheet undergoes an O(1) change in position and that the solution is ultimately terminated by a breakdown in the numerical procedure. The corresponding viscous layer shows downstream thickening, but appears to remain well behaved up to the terminal location.

Description: The effects of multiple phase transitions on mantle convection are investigated by numerical simulations that are based on three-dimensional models. These simulations show that cold sheets of mantle material collide at junctions, merge, and form a strong downflow that is stopped temporarily by the transition zone. The accumulated cold material gives rise to a strong gravitational instability that causes the cold mass to sink rapidly into the lower mantle. This process promotes a massive exchange between the lower and upper mantles and triggers a global instability in the adjacent plume system. This mechanism may be cyclic in nature and may be linked to the generation of superplumes.

Description: This paper investigates the linear stability of confined mixing layers with special emphasis on the effects of heat release and compressibility. The results show that reflection of supersonic disturbances by the walls makes the confined supersonic mixing layer more unstable than the unconfined free shear layer. Decreasing the distance between the walls makes the flow more unstable. However, subsonic disturbances are relatively unaffected by the walls. Heat release and Mach number hardly change the growth rates of supersonic disturbances. The most unstable supersonic disturbances are two-dimensional in rectangular channel flows, but three-dimensional in partially confined flows. Finally, the reactants are not strongly mixed by supersonic instabilities, which mainly disturb one side of the layer.

Description: A new method has been developed to provide a direct test of the error calibrations of gravity models based on actual satellite observations. The basic approach projects the error estimates of the gravity model parameters onto satellite observations, and the results of these projections are then compared with data residual computed from the orbital fits. To allow specific testing of the gravity error calibrations, subset solutions are computed based on the data set and data weighting of the gravity model. The approach is demonstrated using GEM-T3 to show that the gravity error estimates are well calibrated and that reliable predictions of orbit accuracies can be achieved for independent orbits.

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: New research on hypersonic vehicles, such as the National Aero-Space Plane (NASP), has raised concerns about the effects of shock-wave interference on various structural components of the craft. State-of-the-art aerothermal analysis software is inadequate to predict local flow and heat flux in areas of extremely high heat transfer, such as the surface impingement of an Edney-type supersonic jet. EASI revives and updates older computational methods for calculating inviscid flow field and maximum heating from shock wave interference. The program expands these methods to solve problems involving the six shock-wave interference patterns on a two-dimensional cylindrical leading edge with an equilibrium chemically reacting gas mixture (representing, for example, the scramjet cowl of the NASP). The inclusion of gas chemistry allows for a more accurate prediction of the maximum pressure and heating loads by accounting for the effects of high temperature on the air mixture. Caloric imperfections and specie dissociation of high-temperature air cause shock-wave angles, flow deflection angles, and thermodynamic properties to differ from those calculated by a calorically perfect gas model. EASI contains pressure- and temperature-dependent thermodynamic and transport properties to determine heating rates, and uses either a calorically perfect air model or an 11-specie, 7-reaction reacting air model at equilibrium with temperatures up to 15,000 K for the inviscid flowfield calculations. EASI solves the flow field and the associated maximum surface pressure and heat flux for the six common types of shock wave interference. Depending on the type of interference, the program solves for shock-wave/boundary-layer interaction, expansion-fan/boundary-layer interaction, attaching shear layer or supersonic jet impingement. Heat flux predictions require a knowledge (from experimental data or relevant calculations) of a pertinent length scale of the interaction. Output files contain flow-field information for the various shock-wave interference patterns and their associated maximum surface pressure and heat flux predictions. EASI is written in FORTRAN 77 for a DEC VAX 8500 series computer using the VAX/VMS operating system, and requires 75K of memory. The program is available on a 9-track 1600 BPI magnetic tape in DEC VAX BACKUP format. EASI was developed in 1989. DEC, VAX, and VMS are registered trademarks of the Digital Equipment Corporation.

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: Current investigations of the hydrogen-fueled supersonic combustion ramjet engine have delineated several technological problem areas. One area, the analysis of the injection, turbulent mixing, and combusiton of hydrogen, requires the accurate calculation of the supersonic combustion flow fields. This calculation has proven difficult because of an interesting phenomena which makes possible the transition from supersonic to subsonic flow in the combustion field, due to the temperature transitions which occur in the flow field. This computer program was developed to use viscous characteristics theory to analyze supersonic combustion flow fields with imbedded subsonic regions. Intended to be used as a practical design tool for two-dimensional and axisymmetric supersonic combustor development, this program has proven useful in the analysis of such problems as determining the flow field of a single underexpanded hydrogen jet, the internal flow of a gas sampling probe, the effects of fuel-injector strut shape, and the effects of changes in combustor configuration. Both combustion and diffusive effects can significantly alter the wave pattern in a supersonic field and generate significant pressure gradients in both the axial and radial directions. The induced pressure, in turn, substantially influences the ignition delay and reaction times as well as the velocity distribution. To accurately analyze the flow fields, the effects of finite rate chemistry, mixing, and wave propagation must be properly linked to one another. The viscous characteristics theory has been used in the past to describe flows that are purely supersonic; however, the interacting pressure effects in the combustor often allow for the development of shock waves and imbedded subsonic regions. Numerical investigation of these transonic situations has required the development of a new viscous characteristics procedure which is valid within the subsonic region and can be coupled with the standard viscous characteristics procedure in the supersonic region. The basic governing equations used are the 'viscous-inviscid' equations, similar to those employed in higher-order boundary layer analyses, with finite rate chemistry terms included. In addition, the Rankine-Hugoniot and Prandtl-Meyer relations are used to compute shock and expansion conditions. The program can handle up to 20 simultaneous shock waves. Chemistry terms are computed for a 7-species 8-mechanism hydrogen-air reaction scheme. The user input consists of a physical description of the combustor and flow determination parameters. Output includes detail flow parameter values at selected points within the flow field. This computer program is written in FORTRAN IV for batch execution and has been implemented on a CDC CYBER 175 with a central memory requirement of approximately 114K (octal) of 60 bit words. The program was developed in 1978.

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: SINDA, the Systems Improved Numerical Differencing Analyzer, is a software system for solving lumped parameter representations of physical problems governed by diffusion-type equations. SINDA was originally designed for analyzing thermal systems represented in electrical analog, lumped parameter form, although its use may be extended to include other classes of physical systems which can be modeled in this form. As a thermal analyzer, SINDA can handle such interrelated phenomena as sublimation, diffuse radiation within enclosures, transport delay effects, and sensitivity analysis. FLUINT, the FLUid INTegrator, is an advanced one-dimensional fluid analysis program that solves arbitrary fluid flow networks. The working fluids can be single phase vapor, single phase liquid, or two phase. The SINDA'85/FLUINT system permits the mutual influences of thermal and fluid problems to be analyzed. The SINDA system consists of a programming language, a preprocessor, and a subroutine library. The SINDA language is designed for working with lumped parameter representations and finite difference solution techniques. The preprocessor accepts programs written in the SINDA language and converts them into standard FORTRAN. The SINDA library consists of a large number of FORTRAN subroutines that perform a variety of commonly needed actions. The use of these subroutines can greatly reduce the programming effort required to solve many problems. A complete run of a SINDA'85/FLUINT model is a four step process. First, the user's desired model is run through the preprocessor which writes out data files for the processor to read and translates the user's program code. Second, the translated code is compiled. The third step requires linking the user's code with the processor library. Finally, the processor is executed. SINDA'85/FLUINT program features include 20,000 nodes, 100,000 conductors, 100 thermal submodels, and 10 fluid submodels. SINDA'85/FLUINT can also model two phase flow, capillary devices, user defined fluids, gravity and acceleration body forces on a fluid, and variable volumes. SINDA'85/FLUINT offers the following numerical solution techniques. The Finite difference formulation of the explicit method is the Forward-difference explicit approximation. The formulation of the implicit method is the Crank-Nicolson approximation. The program allows simulation of non-uniform heating and facilitates modeling thin-walled heat exchangers. The ability to model non-equilibrium behavior within two-phase volumes is included. Recent improvements to the program were made in modeling real evaporator-pumps and other capillary-assist evaporators. SINDA'85/FLUINT is available by license for a period of ten (10) years to approved licensees. The licensed program product includes the source code and one copy of the supporting documentation. Additional copies of the documentation may be purchased separately at any time. SINDA'85/FLUINT is written in FORTRAN 77. Version 2.3 has been implemented on Cray series computers running UNICOS, CONVEX computers running CONVEX OS, and DEC RISC computers running ULTRIX. Binaries are included with the Cray version only. The Cray version of SINDA'85/FLUINT also contains SINGE, an additional graphics program developed at Johnson Space Flight Center. Both source and executable code are provided for SINGE. Users wishing to create their own SINGE executable will also need the NASA Device Independent Graphics Library (NASADIG, previously known as SMDDIG; UNIX version, MSC-22001). The Cray and CONVEX versions of SINDA'85/FLUINT are available on 9-track 1600 BPI UNIX tar format magnetic tapes. The CONVEX version is also available on a .25 inch streaming magnetic tape cartridge in UNIX tar format. The DEC RISC ULTRIX version is available on a TK50 magnetic tape cartridge in UNIX tar format. SINDA was developed in 1971, and first had fluid capability added in 1975. SINDA'85/FLUINT version 2.3 was released in 1990.

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: Spectral local isotropy tests are applied to direct numerical simulation data, mainly at the centerline of a fully developed turbulent channel flow. Despite the small Reynolds number of the simulation, the high-wavenumber behavior of velocity and vorticity spectra is consistent with local isotropy. This consistency is verified by the relationship between streamwise wavenumber spectra and spanwise wavenumber spectra. The high-wavenumber behavior of the pressure spectrum is also consistent with local isotropy and compares favorably with the calculation of Batchelor (1951), which assumes isotropy and joint normality of the velocity field at two points in space. The latter assumption is validated by the shape but not the magnitude of the quadruple correlation of the streamwise velocity fluctuation at small separations. There is only partial support for local spectral isotropy away from the centerline as the magnitude of the mean strain rate increases.

Description: The stability of TIROS-N Microwave Sounding Unit (MSU) channel 4, which monitors the deep-layer averaged temperature of the lower stratosphere, is tested by intercalibrating MSU channel 4 data from the NIROS-N series of NOAA satellites during 1979-1991. The monthly gridpoint anomalies are validated with 10 years of radiosonde data during 1979-1988. The results demonstrated that the satellite sensors are very stable in their calibration and that the previously reported uncertainties in the stratospheric temperature information produced by NOAA are not the result of calibration changes in the MSUs. It was found that the largest globally averaged temperature variations during 1979-1991 occur after the El Chichon (1982) and Pinatubo (1991) volcanic eruptions. These warm events are superimposed upon a net downward trend in temperatures during the period. The cooling trend is strongest in polar regions and the Northern Hemisphere middle latitudes.

Description: Remote sensing of aerosol optical thickness from space is difficult over continental surfaces. There are two retrieval algorithms, one based on the use of dark targets and a second based on contrast reduction between selected pixels. Improvements in the contrast reduction method are reported. A procedure is developed for using the satellite image to evaluate whether conditions for applying the structure method are met. The theoretical background is discussed, and the usefulness of the structure functions is demonstrated. The method is applied to NOAA Advanced Very High Resolution Radiometer (AVHRR) imagery where simultaneous ground measurements are used for validation.

Description: An MHD simulation was performed to obtain a self-consistent model of magnetic field and plasma density near the X point reconnection region. The MHD model was used to perform extensive ray tracing calculations in order to clarify the propagation characteristics of the plasma waves near the X point reconnection region. The dynamic wave spectra possibly observed by the Geotail spacecraft during a typical cross-tail trajectory are reconstructed. By comparing the extensive ray tracing calculations with the plasma wave data from Geotail, it is possible to perform a kind of 'remote sensing' to identify the location and structure of potential X point reconnection regions.

Description: A previously derived analytical model for the small-scale structure of turbulence is reformulated in such a way that the energy spectrum may be computed. The model is an ensemble of two-dimensional (2D) vortices with internal spiral structure, each stretched by an axially symmetric strain flow. Stretching and differential rotation produce an energy cascade to smaller scales in which the stretching represents the effect of instabilities and the spiral structure is the source of dissipation at the end of the cascade. The energy spectrum of the resulting flow may be expressed as a time integration involving only the enstrophy spectrum of the time evolving 2D cross section flow, which may be obtained numerically. Examples are given in which a k exp -5/3 spectrum is obtained by this method. The k exp -5/3 inertial range spectrum is shown to be related to the existence of a self-similar enstrophy preserving range in the 2D enstrophy spectrum. The results are found to be insensitive to time dependence of the strain rate, including even intermittent on-or-off strains.

Description: A k-epsilon model is proposed for wall bonded turbulent flows. In this model, the eddy viscosity is characterized by a turbulent velocity scale and a turbulent time scale. The time scale is bounded from below by the Kolmogorov time scale. The dissipation equation is reformulated using this time scale and no singularity exists at the wall. The damping function used in the eddy viscosity is chosen to be a function of R(sub y) = (k(sup 1/2)y)/v instead of y(+). Hence, the model could be used for flows with separation. The model constants used are the same as in the high Reynolds number standard k-epsilon model. Thus, the proposed model will be also suitable for flows far from the wall. Turbulent channel flows at different Reynolds numbers and turbulent boundary layer flows with and without pressure gradient are calculated. Results show that the model predictions are in good agreement with direct numerical simulation and experimental data.

Description: The 1992 global average total ozone, measured by the Total Ozone Mapping Spectrometer (TOMS) on the Nimbus-7 satellite, was 2 to 3 percent lower than any earlier year observed by TOMS (1979 to 1991). Ozone amounts were low in a wide range of latitudes in both the Northern and Southern Hemispheres, and the largest decreases were in the regions from 10 to 20 deg S and 10 to 60 deg N. Global ozone in 1992 is at least 1.5 percent lower than would be predicted by a statistical model that includes a linear trend and accounts for solar cycle variation and the quasi-biennial oscillation. These results are confirmed by comparisons with data from other ozone monitoring instruments: the SBUV/2 instrument on the NOAA-11 satellite, the TOMS instrument on the Russian Meteor-3 satellite, the World Standard Dobson Instrument 83, and a collection of 22 ground-based Dobson instruments.

Description: Results are presented of an experimental and numerical study of a turbulent boundary layer with pressure gradients conducted using the recent 'fringe method' with its numerical advantages and good inflow quality. After an inflow transient good agreement is observed; the differences, of up to 13 percent, are discussed. Moderate deviations from the law of the wall are found in the velocity profiles of the simulation. They are fully correlated with the pressure gradient, are in fair quantitative agreement with the experimental results of Nagano et al. (1992), and are roughly the opposite of uncorrected mixing-length-model predictions. Large deviations from the wall scaling are observed for other quantities, notably for the turbulence dissipation rate. The a(1) structure parameter drops mildly in the upper layer with adverse pressure gradient.

Description: The paper formulates a tensorially consistent near-wall second-order closure model. Redistributive terms in the Reynolds stress equations are modeled by an elliptic relaxation equation in order to represent strongly nonhomogeneous effects produced by the presence of walls; this replaces the quasi-homogeneous algebraic models that are usually employed, and avoids the need for ad hoc damping functions. The model is solved for channel flow and boundary layers with zero and adverse pressure gradients. Good predictions of Reynolds stress components, mean flow, skin friction, and displacement thickness are obtained in various comparisons to experimental and direct numerical simulation data. The model is also applied to a boundary layer flowing along a wall with a 90-deg, constant-radius, convex bend.

Description: A computer code for solving the Reynolds-averaged full Navier-Stokes equations has been developed and applied using H- and C-type grids. The Baldwin-Lomax eddy-viscosity model is used for turbulence closure. The integration in time is based on an explicit four-stage Runge-Kutta scheme. Local time stepping, variable coefficient implicit residual smoothing, and a full multigrid method have been implemented to accelerate steady-state calculations. A grid independence analysis is presented for a transonic rotor blade. Comparisons with experimental data show that the code is an accurate viscous solver and can give very good blade-to-blade predictions for engineering applications.

Description: Accounts are given of the development status and prospective efficacy of the Geoscience Environmental Data Display, the Space Environment Laboratory Data Acquisition and Display System, and the Geospace Environment Modeling program, which are all concerned with the 3D definition of the earth's magnetosphere. Attention is given to current and prospective improvements in the integration of all these data-gathering systems.

Description: Activities at the NASA Langley Research Center's distributed active archive centers (DAACs) intended to capitalize on existing centers of scientific expertise and to prevent a single point of failure are described. A Version 0 Langley DAAC, a prototype of an Earth Observing System Data and Information System, started archiving and distributing existing datasets on the earth's radiation budget, clouds, aerosols, and tropospheric chemistry in late 1992. The major goals of the LaRC Version 0 effort include to enhance scientific use of existing data; to develop institutional expertise in maintaining and distributing data; to encourage cooperative interagency and international involvement with datasets and research; and to use institutional capability for processing data from previous missions to prepare for processing the future EOS satellite data.

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: Space geodetic data from very long baseline interferometry (VLBI) was used to estimate velocity relative to the plate interiors of two sites on the deforming leading edge at the Japan trench. Elastic models of interseismic deformation and results obtained were used to put constraints on the slip rate along the main thrust of the Japan subduction zone. Observed velocities reflect the sum of permanent west-northwest shortening in Honshu, elastic deformation due to locking of the main thrust fault at the Japan trench, and deformation associated with the subducting Phillipine plate. These velocities limit the locked segment of the main thrust at the Japan trench to 27 km vertically and 100 km along the dip. This indicates that the main Pacific plate thrust fault is not strongly coupled and probably does not generate strong earthquakes.

Description: Monthly values of the J2 and J3 earth gravitational coefficients were estimated using LAGEOS satellite laser ranging data collected between 1980 and 1989. Monthly variations in gravitational coefficients caused by atmospheric mass redistribution were calculated using measurements of variations in surface atmospheric pressure. Results for correlation studies of the two time series are presented. The LAGEOS and atmospheric J2 time series agree well and it appears that variations in J2 can be attributed to the redistribution of atmospheric mass. Atmospheric and LAGEOS estimates for J3 show poorer agreement, J3 estimates appear to be very sensitive to unmodeled forces acting on the satellite. Results indicate that the LAGEOS data can be used to detect small variations in the gravitational field.

Description: A numerical study was conducted to investigate how bleed through a two-dimensional slot affects shock-wave induced, boundary-layer separation on a flat plate. This study is based on the ensemble-averaged, compressible, Navier-Stokes equations closed by the Baldwin-Lomax, algebraic turbulence model. The algorithm used to obtain solutions was the implicit, partially split, two-factored scheme of Steger. This study examined the effects of the following parameters in controlling shock-wave induced flow separation: location of slot in relation to where the incident shock wave impinged on the boundary layer, size of slot in relation to the boundary-layer thickness, number of slots, spacings between slots, and strength of the incident shock wave. This study also showed the nature of the very complex flowfield about the slot or slots and how the plenum affects the bleed process. The results of this study are relevant to problems where bleed is used to control shock-wave induced, boundary-layer separation (e.g., inside jet engine inlets and wind tunnels).

Description: In his paper, Aikin investigates the hypothesis that high latitude O3 distributions near 1 mbar are modified by the precipitation of relativistic electrons. The authors of this comment believe that Aikin's results are based on an incomplete analysis and contain misperceptions regarding studies done by the authors. Aikin's reply follows this comment.

Description: HALOE observations of O3, CH4, HF, H2O, NO, NO2, and HCl collected during the October 1991 Antarctic spring period are reported. The data show a constant CH4 mixing ratio of about 0.25 ppmv for the altitude range from 65 km down to about 25 km at the position of minimum wind speed in the vortex: i.e., the vortex center, and depressions in pressure versus longitude contours of NO, NO2, HF, and HCl in this same region. Water vapor, HF, and HCl enhancement are also observed in the vortex center region above about 25 km. Between 10 and 20 km, the expected mixing ratio signatures exist within the vortex, i.e., low ozone and dehydration. The water vapor increased by 50 percent, and the ozone level doubled inside the vortex between October 11 and 24 in the 15-20 km layer. These changes imply a time constant for recovery from ozone hole conditions of 19 and 30 days for O3 and H2O, respectively. The data further show the presence of air inside the vortex between 3 and 30 mb which has mixing ratios characteristic of midlatitudes.

Description: Turbulent flow in a 2D channel with repeated rectangular rib roughness was numerically simulated using a low Reynolds number form of the k-epsilon turbulence model. Friction factors and average Stanton numbers were calculated for various pitch to rib height ratios and bulk Reynolds numbers. Comparisons with experiment were generally adequate, with the predictions of friction superior to those for heat transfer. The effect of variable properties for channel flow was investigated, and the results showed a greater effect for friction than for heat transfer. Comparison with experiment yielded no clear conclusions. The turbulence model was also validated for a related problem, that of flow downstream of an abrupt pipe expansion.

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: The purpose of this note is to construct a local solution that eliminates a residual velocity discontinuity in the inviscid portion of a solution obtained in a recent paper by Goldstein, Leib and Cowley (1992). This result is of importance because it shows that the solution obtained in that paper is entirely non-singular outside the viscous wall boundary layer and that any singularity in the problem will have to arise in the usual way through a breakdown in the viscous boundary layer.

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: The paper develops a simple representation of the global circuit of Birkeland currents on the basis of a representation of the current density (j) in terms of Euler potentials (alpha, chi). The underlying magnetic field, which shares the potential alpha with j, is assumed to be dipolar, making the model applicable chiefly to region 2 Birkeland currents. A form of j is chosen that gives a current sheet with peak outflow at dawn and peak inflow at dusk (or vice versa), connected across a flat polar cap sheet. The superposition of harmonics of the same type, centered at different 'foci', are found to provide a flexible and powerful representation of harmonic functions, accurate within less than 1 percent. An interpolation formula by which current sheets of finite width could be consistently represented is developed.

Description: We use a set of equations, sometimes referred to as the 'delta-equations', to approximate the two-dimensional inviscid motion of an initially circular vortex sheet released from rest in a cross-flow. We present numerical solutions of these equations for the case with delta-square = 0 (for which the equations are exact) and for delta-square greater than 0. For small values of the smoothing parameter delta, a spectral filter must be used to eliminate spurious instabilities due to round-off error. Two singularities appear simultaneously in the vortex sheet when delta-square = 0 at a critical time t(c). After t(c), the solutions do not converge as the computational mesh is refined. With delta-square greater than 0, converged solutions were found for all values of delta-square when t is less than t(c), and for all but the two smallest values of delta-square used when t is greater than t(c). Our results show that, when delta-square is greater than 0, the vortex sheet deforms into two doubly branched spirals some time after t(c). The limiting solution as delta approaching 0 clearly exists and equals the delta = 0 solution when t is less than t(c).

Description: The evolution of three-dimensional temporally evolving plane mixing layers through as many as three pairings has been simulated numerically. All simulations were begun from a few low-wavenumber disturbances, usually derived from linear stability theory, in addition to the mean velocity. Three-dimensional perturbations were used with amplitudes ranging from infinitesimal to large enough to trigger a rapid transition to turbulence. Pairing is found to inhibit the growth of infinitesimal three-dimensional disturbances, and to trigger the transition to turbulence in highly three-dimensional flows. The mechanisms responsible for the growth of three-dimensionality and onset of transition to turbulence are described. The transition to turbulence is accompanied by the formation of thin sheets of spanwise vorticity, which undergo secondary rollups. The post-transitional simulated flow fields exhibit many properties characteristic of turbulent flows.

Description: An analytical solution was obtained by Siegel (1991, 1992) for local boundary heat fluxes by a radiating medium at uniform temperature in a 2D rectangular region. It is shown here that, after local fluxes from the medium to the walls have been evaluated, it is very easy to compute local fluxes arriving from the adjacent and opposite walls. This extends the previous analysis and provides convenient relations to include radiation from a black boundary, each side of the rectangle being at a different uniform temperature. The final expressions are helpful in performing spectral calculations that must be made for many spectral bands.

Description: Effects of helicity in three-dimensional incompressible inhomogeneous turbulence are examined with the aid of a two-scale direct-interaction approximation (DIA). The turbulent helicity gives a measure of the reflectional asymmetry in a turbulent flow and its inhomogeneity contributes to the sustainment of large-scale vorticity field in a three-dimensional mean flow. The importance of helicity effects is discussed in the context of flows in a rotating system and swirling flows in a pipe. A three-equation model with the turbulent helicity incorporated is proposed using the theoretical results. The validity of the model is confirmed quantitatively through the application to a decaying swirling flow in a pipe.

Description: The use of satellite-derived imagery for measuring subresolution horizontal terrain displacments associated with present-day earthquakes is discussed with reference to data from the French SPOT satellite whose sensor array provides 10-m panchromatic imagery. The measured terrain displacements can be up to several meters, but usually no more than 6-8 m even for major earthquakes. The general approach is to spatially match the after image to the before image at each point on a half-kilometer grid by iteratively interpolating one and testing its correlation with the other. The discussion covers the basic algorithm, results of initial tests, error types and limitations, and future work.

Description: A new approach to stratigraphic analysis is described which uses photogeologic and spectral interpretation of multispectral remote sensing data combined with topographic information to determine the attitude, thickness, and lithology of strata exposed at the surface. The new stratigraphic procedure is illustrated by examples in the literature. The published results demonstrate the potential of spectral stratigraphy for mapping strata, determining dip and strike, measuring and correlating stratigraphic sequences, defining lithofacies, mapping biofacies, and interpreting geological structures.

Description: The transient response of a thermal protection material to heat applied to the surface can be calculated using the CHAP III computer program. CHAP III can be used to analyze pyrolysis gas chemical kinetics in detail and examine pyrolysis reactions-indepth. The analysis includes the deposition of solid products produced by chemical reactions in the gas phase. CHAP III uses a modelling technique which can approximate a wide range of ablation problems. The energy equation used in CHAP III incorporates pyrolysis (both solid and gas reactions), convection, conduction, storage, work, kinetic energy, and viscous dissipation. The chemically reacting components of the solid are allowed to vary as a function of position and time. CHAP III employs a finite difference method to approximate the energy equations. Input values include specific heat, thermal conductivity, thermocouple locations, enthalpy, heating rates, and a description of the chemical reactions expected. The output tabulates the temperature at locations throughout the ablator, gas flow within the solid, density of the solid, weight of pyrolysis gases, and rate of carbon deposition. A sample case is included, which analyzes an ablator material containing several pyrolysis reactions subjected to an environment typical of entry at lunar return velocity. CHAP III is written in FORTRAN IV for batch execution and has been implemented on a CDC CYBER 170 series computer operating under NOS with a central memory requirement of approximately 102K (octal) of 60 bit words. This program was developed in 1985.

Description: This program performs a one-dimensional numerical analysis of the transient thermal response of multi-layer insulative systems. The analysis can determine the temperature distribution through a system consisting of from one to four layers, one of which can be an air gap. Concentrated heat sinks at any interface can be included. The computer program based on the analysis will determine the thickness of a specified layer that will satisfy a temperature limit criterion at any point in the insulative system. The program will also automatically calculate the thickness at several points on a system and determine the total system mass. This program was developed as a tool for designing thermal protection systems for high-speed aerospace vehicles but could be adapted to many areas of industry involved in thermal insulation systems. In this package, the equations describing the transient thermal response of a system are developed. The governing differential equation for each layer and boundary condition are put in finite-difference form using a Taylor's series expansion. These equations yield an essentially tridiagonal matrix of unknown temperatures. A procedure based on Gauss' elimination method is used to solve the matrix. This program is written in FORTRAN IV for the CDC RUN compiler and has been implemented on a CDC 6000 series machine operating under SCOPE 3.0. This program requires a minimum of 44K (octal) of 60 bit words of memory.

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 computer program has been developed to analyze the transient response of an ablating axisymmetric body, including the effect of shape change. The governing differential equation, the boundary conditions for the analysis on which the computer program is based, and the method of solution of the resulting finite-difference equations are discussed in the documentation. Some of the features of the analysis and the associated program are (1) the ablation material is considered to be orthotropic with temperature-dependent thermal properties; (2) the thermal response of the entire body is considered simultaneously; (3) the heat transfer and pressure distribution over the body are adjusted to the new geometry as ablation occurs; (4) the governing equations and several boundary-condition options are formulated in terms of generalized orthogonal coordinates for fixed points in a moving coordinate system; (5) the finite-difference equations are solved implicitly; and (6) other instantaneous body shapes can be displayed with a user-supplied plotting routine. The physical problem to be modeled with the analysis is described by FORTRAN input variables. For example, the external body geometry is described in the W, Z coordinates; material density is given; and the stagnation cold-wall heating rate is given in a time-dependent array. Other input variables are required which control the solution, specify boundary conditions, and determine output from the program. The equations have been programmed so that either the International System of Units or the U. S. Customary Units may be used. This program is written in FORTRAN IV for batch execution and has been implemented on a CDC 6000 Series computer. This program was developed in 1972.

Description: The GTRAN program was developed to solve transient, as well as steady state, problems for gas piping systems. GTRAN capabilities allow for the analysis of a variety of system configurations and components. These include: multiple pipe junctions; valves that change position with time; fixed restrictions (orifices, manual valves, filters, etc.); relief valves; constant pressure sources; and heat transfer for insulated piping and piping subjected to free or forced convection. In addition, boundary conditions can be incorporated to simulate specific components. The governing equations of GTRAN are the one-dimensional transient gas dynamic equations. The three equations for pressure, velocity, and density are reduced to numerical equations using an implicit Crank-Nicholson finite difference technique. Input to GTRAN includes a description of the piping network, the initial conditions, and any events (e.g. valve closings) occuring during the period of analysis. Output includes pressure, velocity, and density versus time. GTRAN is written in FORTRAN 77 for batch execution and has been implemented on a DEC VAX series computer. GTRAN was developed in 1983.

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: The JPL MkIV interferometer, a Fourier Transform Infra-Red (FTIR) spectrometer designed specifically for atmospheric remote sensing, made measurements of the composition of the Arctic stratosphere in January, February and March 1992. These measurements were made from the NASA DC-8 aircraft as part of the AASE2 campaign. The data reveal that despite 5 to 6 km of subsidence inside the vortex, which more than doubled the vertically integrated column amounts (burdens) of HF and HNO3 with respect to outside the vortex, considerable losses of NO2, HCl and ClNO3 were evident by mid-January. Temporary freeze-out of HNO3 was observed only on one occasion, Jan. 19, and was accompanied by substantial reductions in HCl and ClNO3. During February and March, ClNO3 and NO2 amounts increased dramatically. HCl also recovered but at a much slower rate, so that by March ClNO3 was the major reservoir of inorganic chlorine, at times exceeding HCl by a factor 2.

Description: A Runge-Kutta formula in time is presently used to advance schemes in which central differences are used to solve the time-dependent Euler equations; a second difference is added near shocks as an artificial viscosity to reduce the given scheme to a first-order upwind one at shocks. A matrix-valued dissipation is introduced and compared with the scalar viscosity; a connection is shown between this artificial viscosity and flux limiters. The use of various flux limiters for this central difference scheme is compared.

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: The joint hypocentral determination method is used to relocate deep seismicity reported in the International Seismological Center catalog for earthquakes deeper than 400 km in the Honshu, Bonin, Mariannas, Java, Banda, and South America subduction zones. Each deep seismic zone is found to display planar features of seismicity parallel to the Harvard centroid-moment tensor nodal planes, which are identified as planes of shear failure. The sense of displacement on these planes is one of resistance to deeper penetration.

Description: Stratospheric aerosol collections by wire impactors, taken mainly over the western U.S. from early 1984 to late 1986, show the diminishing effects of El Chichon's 1982 eruption, and provide a set of data for judging subsequent volcanic effects. Decrease in sulfate burden during the study time is due to preferential gravitational settling-out of large (above 0.5-micron diameter) sulfate aerosol droplets. As a result of settling from higher levels, lower altitude (12.2-15.2 km) air early in 1984 tends to contain more sulfate than higher level (19.8-21.3 km) air. As of late 1986, however, high- and low-altitude sulfate contents have decreased and are similar, suggesting large-particle settling has been completed. The later sulfate collection size distributions resemble unimodal background spectra, whereas earlier ones are bimodal. Average sulfate load for all altitudes decreases during the period of study from 0.4 to 0.08 microgram/cu m. The latter value is somewhat higher than a volcanically unenriched pre-El Chichon level, suggesting that even as of 1987, stratospheric background had not been obtained.

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: Using a comprehensive ionospheric data set comprised of all available ion composition and plasma temperature measurements from satellites, the vertical distributions of ion composition and plasma temperatures are defined from middle latitudes up into the polar cap for summer conditions for altitudes below about 1200 km. These data are sufficient to allow a numerical estimation of the latitudinal variation of the light ion outflows from within the plasmasphere to the polar wind regions. The altitude at which significant light ion outflow begins is found to be lower during solar minimum conditions than during solar maximum. The H(+) outward speeds are of the order of 1 km/s near 1100 km during solar maximum but attain several km/s speeds for solar minimum. He(+) shows a similar altitude development of flow but attains polar cap speeds much less than 1 km/s at altitudes below 1100 km, particularly under solar maximum conditions. Outward flows are also found in the topside F-region for noontime magnetic flux tubes within the plasmasphere.

Description: The results of limitation studies performed with the UARS MLS are presented. A consistent set of algorithms allows the extraction of the spectral coefficients in time and longitude from asynoptically sampled satellite data and the subsequent reconstruction of synoptic maps from that spectral information. In addition to providing synoptic maps, the asynoptic technique allows the use of standard spectral analysis tools such as autocorrelation and cross correlation.

Description: Airglow and chemical processes in the terrestrial mesosphere and lower thermosphere are reviewed, and initial parameterizations of the processes applicable to multidimensional models are presented. The basic processes by which absorbed solar energy participates in middle atmosphere energetics for absorption events in which photolysis occurs are illustrated. An approach that permits the heating processes to be incorporated in numerical models is presented.

Description: Attention is given to a program to compare stratospheric models and measurements to calculate the atmospheric effects of emissions of a projected fleet of High Speed Civil Transports on the stratosphere is presented. The goals of the effort is to establish a standard set of atmospheric measurements that can be used to test the reliability of atmospheric chemistry models, develop a method for evaluating model/data comparisons, and direct the first major international stratospheric model/data comparison.

Description: Multispectra aerosol extinction data for the fall and spring of 1987 measured by the SAGE II sensor are employed to determine the physical characteristics of aerosols within the springtime Antarctic polar vortex. Attention is given to the physical processes that give rise to the apparent springtime 'cleansing' of the Antarctic stratosphere. The inferred vertical and radial structure compare favorably with in situ measurements but yield a previously unavailable 2D structure to the distribution of aerosols within the polar vortex. The springtime 'cleansing' of the Antarctic stratosphere is found to be a result of both large-scale subsidence and the preferential removal of large particles by the nucleation and subsequent sedimentation of polar stratospheric clouds.

Description: The NASA Langley 3D GCM chemical transport model is used to investigate the distribution of atmospheric N2O up to 60 km altitude. The transport characteristics of the model is evaluated without the complications of a detailed chemical formulation for all of the relevant stratospheric minor constituents. Interpretation of the yearly average zonal mean N2O distribution in terms of transport by the yearly averaged meridional circulation and stratospheric photochemical loss indicates large regions in the Northern Hemisphere stratosphere where dynamical mixing apparently plays a large role in maintaining the N2O distribution. In these regions, slopes of the N2O mixing ratio isopleths are maintained by competition between advection by the meridional circulation acting to steepen and dynamical mixing acting to flatten the slopes.