Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reproducibility”

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

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

At least 1,387 records · Page 77
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

A Composite View of Lower Stratospheric Ozone Developed Using a Chemistry Transport Model and Observations from Airborne Lidar and Satellite

An ozone simulation from the Goddard three-dimensional chemistry and transport model for the 1995-96 northern hemisphere winter is compared with ozone observations from airborne Differential Absorption Lidar (DIAL), from the Polar Ozone and Aerosol Measurement (POAM), from the Microwave Limb Sounder (MLS), and from the Halogen Occultation Experiment (HALOE). The 3D model uses winds from the Goddard Data Assimilation System. The 3D model reproduces the latitude dependence of the horizontal and vertical ozone gradients of the subtropical DIAL observations. Comparisons with subtropical satellite observations, which lack the spatial resolution of DIAL but provide near continuous coverage throughout the subtropics, show that the model also reproduces longitude and temporal dependence in the tropical-midlatitude boundary. At polar latitudes, observations from DIAL flights on December 9 and January 30, and POAM and MLS between late December and late January are compared with the 3D model. Data from the three platforms consistently show that the observed ozone has a negative trend relative to the modeled ozone, and that the trend is uniform in time between early and mid winter, with no obvious dependence on proximity to the vortex edge.

Douglass, A. R.↗
NASA NTRS1999

Ocean Turbulence I: One-Point Closure Model Momentum and Heat Vertical Diffusivities

Since the early forties, one-point turbulence closure models have been the canonical tools used to describe turbulent flows in many fields. In geophysics, Mellor and Yamada applied such models using the 1980 state-of-the art. Since then, no improvements were introduced to alleviate two major difficulties: 1) closure of the pressure correlations, which affects the correct determination of the critical Richardson number Ri(sub cr) above which turbulent mixing is no longer possible and 2) the need to express the non-local third-order moments (TOM) in terms of lower order moments rather than via the down-gradient approximation as done thus far, since the latter seriously underestimates the TOMs. Since 1) and 2) are still being dealt with adjustable parameters which weaken the credibility of the models, alternative models, not based on turbulence modeling, have been suggested. The aim of this paper is to show that new information, partly derived from the newest 2-point closure model discussed, can be used to solve these shortcomings. The new one-point closure model, which in its simplest form is algebraic and thus simple to implement, is first shown to reproduce a variety of data. Then, it is used in a Ocean-General Circulation Model (O-GCM) where it reproduces well a large variety of ocean data. While phenomenological models are specifically tuned to ocean turbulence, the present model is not. It is first tested against laboratory data on stably stratified flows and then used in an O-GCM. It is more general, more predictive and more resilient, e.g., it can incorporate phenomena like wave-breaking at the surface, salinity diffusivity, non-locality, etc. One important feature that naturally comes out of the new model is that the predicted Richardson critical value Ri(sub cr) is Ri (sub cr approx. = 1) in agreement with both Large Eddy Simulations (LES) and empirical evidence while all previous models predicted Ri (sub cr approx. = 0.2) which led to a considerable underestimate of the extent of turbulent mixing and thus to an incorrect mixed layer depth. The predicted temperature and salinity profiles (vs. depth) are presented and compared with those of the Kolmogorov-Petruvsky-Piskunuv (KPP) model and Levitus data.

Canuto, V. M.↗
NASA NTRS1998

EOS Interdisciplinary Investigation: Observational and Modeling Studies of Radiative, Chemical, and Dynamical Interactions in the Earth's Atmosphere

A 3-dimensional model of atmospheric dynamics and photochemistry has been developed from the primitive equations in isentropic coordinates. The model extends from the upper troposphere through the middle atmosphere and is driven by observed dynamical and chemical structure in the troposphere. The model's formulation is entirely spectral: Horizontal structure is represented in terms of vector Hough functions, which explicitly partition the motion into its rotational and divergent components. This formulation enables the model's computational performance to be increased dramatically by filtering high-frequency gravity waves, which do not affect PV conservation. Vertical structure is represented in terms of eigenfunctions that follow from the primitive equations in isentropic coordinates. The model's fully spectral formulation enables scale-selective dissipation, necessary for numerical stability, to be applied at 6th order-in all 3 coordinate directions. This feature leaves all but the shortest vertical scales undamped and, consequently, allows potential vorticity and chemical tracers to be conserved quite accurately, distinctly better than in the model's layered counterpart. These dynamical features are complemented by a basic but fairly complete treatment of gas phase photochemistry, which accounts for some 48 chemical species, diurnally-varying SW absorption by O2 and O3, and LW cooling calculated via a detailed band calculation of radiative transfer. Driven by observed tropospheric structure, the 3D model reproduces observed structure in the middle atmosphere, inclusive of transport by the Brewer-Dobson circulation and accompanying chemical variations. Calculated residual motion is consistent with diabatic cooling rates, poleward transport, and descent in the polar night derived from satellite measurements of chemical tracers made from LIMS, SAGE, and UARS. The model also reproduces the observed structure of chemical species like O3, HNO3, and ClO. At the same time, it provides a respectable 3D simulation of the wintertime increase of midlatitude total ozone, as well as its geographical structure observed by Nimbus-7 TOMS. The 3D model has been used to explore the origin of diabatic upwelling in the tropical lower stratosphere, which has been related to planetary wave absorption in the winter- time stratosphere, as well as to influences from tropospheric convection at the tropical tropopause.

Salby, Murry L.↗
NASA NTRS1968

Use of a Magnetic Azimuth-Indicator System Developed for Balloon Payloads of the Planetary Entry Parachute Program

An azimuth-sensing system was used for continuous ground monitoring of the azimuth orientation of the balloon-borne Planetary Entry Parachute Program (PEPP) spacecraft. The system utilized two magnetic field sensors located in the spacecraft so as to produce a unique set of voltage outputs for any azimuth. Electronics onboard the spacecraft encoded the magnetometer outputs and fed the encoded signal to a C-band modulator which imposed the intelligence onto the radar tracking beacon pulses. A conveniently located ground radar received the modulated pulses and the magnetometer outputs were reproduced after demodulation. For fast, direct readout, an X-Y plotter was calibrated to cross plot the reproduced signal of the two magnetometers on a combination rectangular-polar graph indicating the correct true azimuth in real time. The method was used successfully on three of the PEPP balloon-borne spacecraft to determine when the spacecraft was pointed in the most advantageous direction for release from the balloon, at 130,000 ft altitude. The use of such a system of azimuth monitoring is suitable to balloon payloads which are virtually stable with respect to the horizontal but may be rotating about the vertical axis. Angular accuracies within plus or minus 10 deg. could generally be expected.

Darnell, Wayne L.↗
NASA NTRS1994

Theoretical Investigation of Wave-Vector-Dependent Analytical and Numerical Formulations of the Interband Impact-Ionization Transition Rate for Electron in Bulk Silicon and GaAs

The electron interband impact-ionization rate for both silicon and gallium arsenide is calculated using an ensemble Monte Carlo simulation with the expressed purpose of comparing different formulations of the interband ionization transition rate. Specifically, three different treatments of the transition rate are examined: the traditional Keldysh formula, a new k-dependent analytical formulation first derived by W. Quade, E Scholl, and M. Rudan, and a more exact, numerical method of Y. Wang and K. F. Brennan. Although the completely numerical formulation contains no adjustable parameters and as such provides a very reliable result, it is highly computationally intensive. Alternatively, the Keldysh formula, although inherently simple and computationally efficient, fails to include the k dependence as well as the details of the energy band structure. The k-dependent analytical formulation of Quade and co-workers overcomes the limitations of both of these models but at the expense of some new parameterization. It is found that the k-dependent analytical method of Quade and co-workers produces very similar results to those obtained with (he completely numerical model for some quantities. Specifically, both models predict that the effective threshold for impact ionization in GaAs and silicon is quite soft, that the majority of ionization events originate from the second conduction band in both materials, and that the transition rate is k dependent. Therefore, it is concluded that the k-dependent analytical model can qualitatively reproduce results similar to those obtained with the numerical model yet with far greater computational efficiency. Nevertheless, there exist some important drawbacks to the k-dependent analytical model of Quade and co-workers: These are that it does not accurately reproduce the quantum yield data for bulk silicon, it requires determination of a new parameter, related physically to (he overlap integrals of the Bloch state which can only be adjusted by comparison to experiment, and fails to account for any wave-vector dependence of the overlap integrals. As such [he transition rate may be overestimated at those points for which "near vertical," small change in k, transitions occur.

Kolnik, Jan↗
NASA NTRS2000

On the Spatial Distribution of High Velocity Al-26 Near the Galactic Center

We present results of simulations of the distribution of 1809 keV radiation from the decay of Al-26 in the Galaxy. Recent observations of this emission line using the Gamma Ray Imaging Spectrometer (GRIS) have indicated that the bulk of the AL-26 must have a velocity of approx. 500 km/ s. We have previously shown that a velocity this large could be maintained over the 10(exp 6) year lifetime of the Al-26 if it is trapped in dust grains that are reaccelerated periodically in the ISM. Here we investigate whether a dust grain velocity of approx. 500 km/ s will produce a distribution of 1809 keV emission in latitude that is consistent with the narrow distribution seen by COMPTEL. We find that dust grain velocities in the range 275 - 1000 km/ s are able to reproduce the COMPTEL 1809 keV emission maps reconstructed using the Richardson-Lucy and Maximum Entropy image reconstruction methods while the emission map reconstructed using the Multiresolution Regularized Expectation Maximization algorithm is not well fit by any of our models. The Al-26 production rate that is needed to reproduce the observed 1809 keV intensity yields in a Galactic mass of Al-26 of approx. 1.5 - 2 solar mass which is in good agreement with both other observations and theoretical production rates.

Sturner, Steven J.↗
NASA NTRS2000

Photochemistry Saturn's Atmosphere: Effects of an Influx of External Oxygen - 2

We use a one-dimensional diurnally averaged model of photochemistry and diffusion in Saturn's stratosphere to investigate the influence of extraplanetary debris on atmospheric chemistry. In particular, we consider the effects of an influx of oxygen from micrometeoroid ablation or from ring-particle diffusion; the contribution from cometary impacts, satellite debris, or ring vapor is deemed to be less important. The photochemical model results are compared directly with Infrared Space Observatory (ISO) observations to constrain the influx of extraplanetary oxygen to Saturn. From the ISO observations, we determine that the column densities of CO2 and H2O above 10 mbar in Saturn's atmosphere are (6.3 +/- 1) x 10(exp 14) and (1.4 +/- 0.4) x 10(exp 15)/ square cm, respectively; our models indicate that a globally averaged oxygen influx of (4+/-2) x 10(exp 6) O atoms /sq cm/s is required to explain these observations. Models with a locally enhanced influx of H20 operating over a small fraction of the projected area do not provide as good a fit to the ISO H2O observations. If volatile oxygen compounds comprise one-third to one-half of the exogenic source by mass, then Saturn is currently being bombarded with (3 +/- 2) x 10(exp -16) g/square cm/s of extraplanetary material. To reproduce the observed CO2/H2O ratio in Saturn's stratosphere, some of the exogenic oxygen must arrive in the form of a carbon-oxygen bonded species such as CO or CO2. An influx consistent with the composition of cometary ices fails to reproduce the high observed CO2/H2O ratio, suggesting that (i) the material has ices that are slightly more carbon-rich than is typical for comets, (ii) a contribution from an organic-rich component is required, or (iii) some of the hydrogen-oxygen bonded material is converted to carbon-oxygen bonded material without photochemistry (e.g., during the ablation process). We have also reanalyzed the 5-micron CO observations of Noll and Larson and determine that the CO lines are most sensitive to the 100- to 400-mbar column density for which we derive a range of (0.7-1.5) x 10(exp 17)/square cm; the CO observations do not allow us to distinguish between an external or internal source of CO on Saturn. If we assume that all the extraplanetary oxygen derives from a micrometeoroid source, then the unfocused dust flux at 9.5 AU must be (i) (1 +/- 0.7) x 10(exp -16) g/square cm/s if interstellar grains are the source of the external oxygen on Saturn, (ii) (4 +/- 3) x 10 (exp -17) g/ sq. cm./s if IDPs on randomly inclined, highly eccentric orbits (e.g., Halley-type comet grains) are the source of the external oxygen, or (iii) (2+/- 1.4) x 10(exp -18)/g/sq cm/s if IDPs on low inclination, low eccentricity orbits (e.g., Kuiper-belt grains) are the source of the external oxygen. These estimates can be used in combination with future Cassini dust detection data to determine the ultimate source of the dust at Saturn's distance from the Sun.

Moses, Julianne I.↗
NASA NTRS2000

Analysis of the NASA White Sands Test Facility (WSTF) Test System for Friction-Ignition of Metallic Materials

Friction is a known ignition source for metals in oxygen-enriched atmospheres. The test system developed by the NASA White Sands Test Facility in response to ASTM G-94 has been used successfully to determine the relative ignition from friction of numerous metallic materials and metallic materials pairs. These results have been ranked in terms of a pressure-velocity product (PV) as measured under the prescribed test conditions. A high value of 4.1(exp 8) watts per square meter for Inconel MA 754 is used to imply resistance to friction ignition, whereas a low value of 1.04(exp 8) watts per square meter for stainless steel 304 is taken as indicating material susceptible to friction ignition. No attempt has been made to relate PV values to other material properties. This work reports the analysis of the WSTF friction-ignition test system for producing fundamental properties of metallic materials relating to ignition through friction. Three materials, aluminum, titanium, and nickel were tested in the WSTF frictional ignition instrument system under atmospheres of oxygen or nitrogen. Test conditions were modified to reach a steady state of operation, that is applied, the force was reduced and the rotational speed was reduced. Additional temperature measurements were made on the stator sample. The aluminum immediately galled on contact (reproducible) and the test was stopped. Titanium immediately ignited as a result of non-uniform contact of the stator and rotor. This was reproducible. A portion of the stator sampled burned, but the test continued. Temperature measurements on the stator were used to validate the mathematical model used for estimating the interface (stator/rotor) temperature. These interface temperature measurements and the associate thermal flux into the stator were used to distinguish material-phase transitions, chemical reaction, and mechanical work. The mechanical work was used to analyze surface asperities in the materials and to estimate a coefficient of fiction. The coefficient of fiction was analyzed in terms of material properties that is, hardness, Young's modulus and elasticity/plasticity of the material.

Shoffstall, Michael S.↗
NASA NTRS2000

The Non-Stellar Infrared Continuum of Seyfert Galaxies

JHKL'M (1 - 5 micrometers) imaging of a sample of Seyfert 2 galaxies is presented. We have performed an accurate estimate of the near-infrared non-stellar nuclear fluxes. We confirm that the near-infrared nuclear continuum between 1 and 2.2microns of some Seyfert 2s is dominated by stellar emission, whereas the continuum emission at longer wavelengths (lambda = 3 - 5 micrometers) is almost entirely non-stellar in origin. The non-stellar spectral energy distributions (SED) in the infrared (up to 15 micrometers) of Seyfert galaxies show a variety of shapes, and they are well reproduced with the tapered disk models of Efstathiou & Rowan-Robinson (1995). We have used two models, one including an optically thin cone component found to fit the SED of NGC 1068, and a coneless model. Although our modelling of the SEDs does not allow us to favor either model to account for all the observed SEDs, we find that the viewing angle towards the central source is well constrained by both models. The galaxies in our sample have fitted values of the viewing angle in the range Theta(sub V) = 0 deg - 64 deg, for the assumed model parameters. We have also investigated non-stellar color-color diagrams (L' - M vs. H - M and L' - M vs. H - L'). The colors of the Seyfert galaxies with viewing angles Theta(sub v) less than 30 deg are better reproduced with the cone model. These diagrams provide a good means to separate Seyfert 2s with moderate obscuration (A(sub V) approx. less than 20 mag from hard X-ray observations) from those with high obscuration. The ground-based 4.8 microns and ISO 9.6 microns luminosities are well correlated with the hard X-ray luminosities of Seyfert ls and 2s. These continuum emissions appear as a good indicator of the AGN luminosity, at least in the cases of hard X-ray Compton-thin Seyfert galaxies (N(sub H) less than or = 10(exp 24)/sq cm). We finally stress the finding that some Compton thick galaxies show bright non-stellar emission at 5 microns This suggests that the near-infrared emission in Seyfert galaxies is produced in an extended component illuminated by the central source, that is more visible from all viewing angles, providing a good explanation for the differing N(sub H)/A(sub V) ratios found in some Seyfert 2s. We discuss possible implications of mid-infrared surveys for the search of counterparts of highly obscured hard X-ray sources.

Alonso-Herrero, Almudena↗
NASA NTRS2001

Bridgman Growth of Detached GeSi Crystals

The growth of detached crystals by the Bridgman technique, in which the growing crystal is not in contact with the crucible wall, has been observed both on earth and in microgravity conditions. Such detachment has particularly been in evidence in microgravity experiments, where the pressure head of the molten sample is absent. At present, the mechanisms contributing to the detachment are not completely understood and until recently detachment has not been reproducibly obtained. Key parameters which must be considered are the contact angle between the melt and the crucible and the growth angle. Another essential parameter is the pressure difference between the annular gap around the solid below the melt and the volume above the melt. Here we present results of the growth of Ge(0.98)Si(0.02) using adjustments in the applied temperature profile to control the pressure difference between the bottom and top of the melt. The technique is less susceptible to sample contamination than controlling the pressure by connecting the crucible to external gas sources. Using this technique, a pressure difference is created by decreasing the temperature in the volume above the melt while the sample is molten but prior to growth. A maximum pressure difference approximately equal to the pressure head of the molten sample can thus be obtained. Several GeSi crystals were grown in pyrolitic boron nitride crucibles. When a pressure difference was applied, samples were reproducibly grown mostly detached. For comparison, samples were also grown in a configuration in which gas could pass freely between the gap below the melt and the volume above the melt and no pressure difference could be established. These samples were initially attached. Existence of detachment was determined both by measuring the radius of the samples with a profilometer and by observations of the sample surfaces with optical and electron microscopy. The gap thickness between the crucible and detached crystal was on the order of 10 micron. The surfaces of the attached areas of the crystals had the same shape as the interior crucible wall whereas in the detached areas the crystal facets could usually be observed.

Volz, M. P.↗
NASA NTRS2001

Climate and Ozone Response to Increased Stratospheric Water Vapor

Stratospheric water vapor abundance affects ozone, surface climate, and stratospheric temperatures. From 30-50 km altitude, temperatures show global decreases of 3-6 K over recent decades. These may be a proxy for water vapor increases, as the Goddard Institute for Space Studies (GISS) climate model reproduces these trends only when stratospheric water vapor is allowed to increase. Observations suggest that stratospheric water vapor is indeed increasing, however, measurements are extremely limited in either spatial coverage or duration. The model results suggest that the observed changes may be part of a global, long-term trend. Furthermore, the required water vapor change is too large to be accounted for by increased production within the stratosphere, suggesting that ongoing climate change may be altering tropospheric input. The calculated stratospheric water vapor increase contributes an additional approximately equals 24% (approximately equals 0.2 W/m(exp 2)) to the global warming from well-mixed greenhouse gases over the past two decades. Observed ozone depletion is also better reproduced when destruction due to increased water vapor is included. If the trend continues, it could increase future global warming and impede stratospheric ozone recovery.

Shindell, Drew T.↗
NASA NTRS2001

Ground Based Experiments in Support of Microgravity Research Results-Vapor Growth of Organic Nonlinear Optical Thin Film

This work is specifically focused on explaining previous results obtained for the crystal growth of an organic material in a reduced gravity environment. On STS-59, in April 1994, two experiments were conducted with N,N-dimethyl-p-(2,2-dicyanovinyl) aniline (DCVA), a promising nonlinear optical (NLO) material. The space experiments were set to reproduce laboratory experiments that yielded small, bulk crystals of DCVA. The results of the flight experiment, however, were surprising. Rather than producing a bulk single crystal, the result was the production of two high quality, single crystalline thin films. This result was even more intriguing when it is considered that thin films are more desirable for NLO applications than are bulk single crystals. Repeated attempts on the ground to reproduce these results were fruitless. A second set of flight experiments was conducted on STS-69 in September 1995. This time eight DCVA experiments were flown, with each of seven experiments containing a slight change from the first reference experiment. The reference experiment was programmed with growth conditions identical to those of the STS-59 mission. The slight variations in each of the other seven were an attempt to understand what particular parameter was responsible for the preference of thin film growth over bulk crystal growth in microgravity. Once again the results were surprising. In all eight cases thin films were grown again, albeit with varying quality. So now we were faced with a phenomenon that not only takes place in microgravity, but also is very robust, resisting all attempts to force the growth of bulk single crystals.

Zugrav, M. Ittu↗
NASA NTRS2000

Ammonia Formation by the Reduction of Nitrite/Nitrate by FeS: Ammonia Formation Under Acidic Conditions

FeS reduces nitrite to, ammonia at pHs lower than the corresponding reduction by aqueous Fe+2. The reduction follows a reasonable first order decay, in nitrite concentration, with a half life of about 150 min (room temperature, CO2, pH 6.25). The highest ammonia product yield measured was 53%. Under CO2, the product yield decreases from pH 5.0 to pH 6.9. The increasing concentration of bicarbonate at higher pH interferes with the reaction. Bicarbonate interference is shown by comparing runs under N2 and CO2. The reaction proceeds well in the presence of such species as chloride, sulfate, and phosphate though the yield drops significantly with phosphate. FeS also reduces nitrate and, unlike with Fe+2, the reduction shows more reproducibility. Again, the product yield decreases with increasing pH, from 7% at pH 4.7 to 0% at pH 6.9. It appears as if nitrate is much more sensitive to the presence of added species, perhaps not competing as well for binding sites on the FeS surface. This may be the cause of the lack of reproducibility of nitrate reduction by Fe+2 (which also can be sensitive to binding by certain species).

Summers, David P.↗
NASA NTRS2001

An Improved Model for the Turbulent PBL

Second order turbulence models of the Mellor and Yamada type have been widely used to simulate the PBL. It is however known that these models have several deficiencies. For example, they all predict a critical Richardson number which is about four times smaller than the Large Eddy Simulation (LES) data, they are unable to match the surface data, and they predict a boundary layer height lower than expected. In the present model, we show that these difficulties are all overcome by a single new physical input: the use of the most complete expression for both the pressure-velocity and the pressure-temperature correlations presently available. Each of the new terms represents a physical process that, was not accounted for by previous models. The new model is presented in three different levels according to Mellor and Yamada's terminology, with new, ready-to-use expressions for the turbulent, moments. We show that the new model reproduces several experimental and LES data better than previous models. As far as the PBL is concerned, we show that the model reproduces both the Kansas data as analyzed by Businger et al. in the context of Monin-Obukhov similarity theory for smaller Richardson numbers, as well as the LES and laboratory data up to Richardson numbers of order unity. We also show that the model yields a higher PBL height than the previous models.

Cheng, Y.↗
NASA NTRS1994

Towards Understanding the Fluid Dynamic Phenomenon of Interest to Rocket Base Heating: A Review

The significance of the base heating problem for rockets during ascent is due to the complex interaction between the rocket nozzle plumes and the external-flow which can change the flow characteristics in the base region dramatically. At lower altitudes the external-flow merges with the plume-flow, without the formation of a large separated flow region, and the cooler external-flow promotes convective cooling of the base wall. Under these conditions the majority of the base heating is due to radiative heating from the shock heated plume gases. At higher altitudes, however, the process of base heating is not so straightforward. The plume and the base flow expands dramatically and separated flow regions occur in the base area. Hot exhaust gases from the rocket nozzle will be entrained into the separated flow regions and produce a convective component to the base wall heating. Further, if the rocket exhaust-gas contains soot, the soot can increase the emission from the gas and dramatically increase the wall absorption coefficient for radiative heating if it is deposited on the walls . In addition, if the rocket exhaust gas is fuel rich, the fuel can bum in the separated flow regions and further increase the base heating. The base burning phenomenon, and the increased base heating caused by it at higher altitudes, have been observed for the Space Shuttle and Saturn Rocket. Under these conditions, the total heating is significantly higher than the heating without separated flow in the base region, and the increase in heating is directly attributable to the fluid dynamic complexity of the base region. Realistic simulation of the base heating requires that the calculated flow environment reproduce the fluid dynamic flow features accurately. Thus, it will be necessary to introduce into the CFD codes the capability for the flow to respond to the complex vehicle geometry, the effect of turbulence, the ability to accurately reproduce the plume shock/shear layer structures and capture such phenomenon as shock induced base separation and base-burning phenomenon. A survey of experimental, theoretical and computational work that details the fluid dynamics of the base flow environment will be presented in the proposed paper. CFD simulations of rocket base flows using standard CFD codes such as OVERFLOW or GASP will be explored to capture these phenomenon accurately. Merits and limitations of these codes for base flow environment predictions will be explored.

Venkatapathy, E.↗
NASA NTRS2002

Numerical Simulation of Screech Tones from Supersonic Jets: Physics and Prediction

The objectives of this project are to: (1) perform a numerical simulation of the jet screech phenomenon; and (2) use the data of the simulations to obtain a better understanding of the physics of jet screech. The original grant period was for three years. This was extended at no cost for an extra year to allow the principal investigator time to publish the results. We would like to report that our research work and results (supported by this grant) have fulfilled both objectives of the grant. The following is a summary of the important accomplishments: (1) We have now demonstrated that it is possible to perform accurate numerical simulations of the jet screech phenomenon. Both the axisymmetric case and the fully three-dimensional case were carried out successfully. It is worthwhile to note that this is the first time the screech tone phenomenon has been successfully simulated numerically; (2) All four screech modes were reproduced in the simulation. The computed screech frequencies and intensities were in good agreement with the NASA Langley Research Center data; (3) The staging phenomenon was reproduced in the simulation; (4) The effects of nozzle lip thickness and jet temperature were studied. Simulated tone frequencies at various nozzle lip thickness and jet temperature were found to agree well with experiments; (5) The simulated data were used to explain, for the first time, why there are two axisymmetric screech modes and two helical/flapping screech modes; (6) The simulated data were used to show that when two tones are observed, they co-exist rather than switching from one mode to the other, back and forth, as some previous investigators have suggested; and (7) Some resources of the grant were used to support the development of new computational aeroacoustics (CAA) methodology. (Our screech tone simulations have benefited because of the availability of these improved methods.)

Tam, Christopher K. W.↗
NASA NTRS2001

Modeling Tides, Planetary Waves, and Equatorial Oscillations in the MLT

Applying Hines Doppler Spread Parameterization for gravity waves (GW), our 3D model reproduces some essential features that characterize the observed seasonal variations of tides and planetary waves in the upper mesosphere. In 2D, our model also reproduces the large Semi-Annual Oscillation (SAO) and Quasi Biennial Oscillation (QBO) observed in this region at low latitudes. It is more challenging to describe these features combined in a more comprehensive self consistent model, and we give a progress report that outlines the difficulties and reports some success. In 3D, the GW's are partially absorbed by tides and planetary waves to amplify them. Thus the waves are less efficient in generating the QBO and SAO at equatorial latitudes. Some of this deficiency is compensated by the fact that the GW activity is observed to be enhanced at low latitudes. Increasing the GW source has the desired effect to boost the QBO, but the effect is confined primarily to the stratosphere. With increasing altitude, the meridional circulation becomes more important in redistributing the momentum deposited in the background flow by the GW's. Another factor involved is the altitude at which the GW's originate, which we had originally chosen to be the surface. Numerical experiments show that moving this source altitude to the top of the troposphere significantly increases the efficiency for generating the QBO without affecting much the tides and planetary waves in the model. Attention to the details in which the GW source comes into play thus appears to be of critical importance in modeling the phenomenology of the MLT. Among the suite of numerical experiments reported, we present a simulation that produced significant variations of tides and planetary waves in the upper mesosphere. The effect is related to the QBO generated in the model, and GW filtering is the likely cause.

Mengel, J. G.↗
NASA NTRS2001

Evidence for the Widespread Presence of Liquid-Phase Particles During the 1999-2000 Arctic Winter

In-situ Multi-Angle Spectrometer Probe (MASP) particle measurements have been analyzed to determine the typical behaviour of sulphate particles during the SOLVE campaign. The analysis has focussed on variations in the total particle concentration measured by MASP, which is a diagnostic that is primarily sensitive to the growth of small particles (those which are initially smaller than 0.2 micrometers in radius). At all levels of the stratosphere, the total MASP concentration varies continuously with temperature. This behaviour is well-reproduced by assuming that the sulphate aerosols are liquid ternary solutions, but cannot be reproduced if the aerosol is assumed to be frozen. This liquid-like behaviour is apparent for all measurements made during SOLVE, both inside and outside the vortex, and even at the coldest temperatures sampled during the campaign. The only anomalous measurements were made during the flight of January 14th, 2001; however, this midlatitude flight is very unlikely to contain SAT particles based on the recent warm temperatures experienced by the air. At the levels with the coldest measured temperatures, at least 90% of the particles grow as liquids. Therefore, the freezing that occurred during the 1999-2000 Arctic winter was selective, causing most of the particles to remain liquid even in the presence of a small number of frozen particles.

Drdla, K↗