Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “correspondence”

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 397 records · Page 22

Measurement of diurnal and semidiurnal rotational variations and tidal parameters of Earth

We discuss the determination of diurnal and semidiurnal variations in the rotation rate and the direction of rotation axis of Earth from the analysis of 8 years of very long baseline interferometry (VLBI) data. This analysis clearly show that these variations are largely periodic and tidally driven; that is, the periods of the variations correspond to the periods of the largest lunar and solar tides. For rotation rate variations, expressed in terms of changes in universal time (UT), the tidal lines with the largest observed signals are O1 (amplitude 23.5 microseconds in time (microseconds), period 25.82 solar hours); KL (18.9 microseconds, 23.93 hours); M2 (17.9 microseconds, 12.54 hours); and S2 (8.6 microseconds, 12.00 hours). For variations in the direction of the rotation axis (polar motion), significant signals exist in the retrograde semidiurnal band at the M2 and S2 tides (amplitudes 265 and 119 microarc seconds (microarc seconds, respectively); the prograde diurnal band at the O1, K1, and P1 tides (amplitudes 199, 152, and 60 microarc seconds, respectively); and the prograde semidiurnal band at the M2 and K2 tides (amplitudes 58 and 39 microarc seconds, respectively). Variations in the retrograde diurnal band are represented by corrections with previous estimates except that a previously noted discrepancy in the 13.66-day nutation (corresponding to the O1 tide) is largely removed in this new analysis. We estimate that the standard deviations of these estimates are 1.0 microseconds for the UT1 variations and 14-16 microarc seconds for the polar motion terms. These uncertainties correspond to surface displacements of approximately 0.5 mm. From the analysis of atmospheric angular momentum data we conclude that variations in UT1 excited by the atmosphere with subdaily periods are small (approximately 1 microsecond). We find that the average radial tidal displacements of the VLBI sites in the diurnal band are largely consistent with known deficiencies in current tidal models, i.e., deficiencies of up to 0.9 mm in the treatment of the free core nutation resonance. In the semidiurnal band, our analysis yields estimates of the second-degree harmonic radial Love number h(sub 2) at the M2 tide of 0.604 + i0.005 +/- 0.002. The most likely explanation for the rotational variations are the effects of ocean tides, but there may also be some contributions from atmospheric tides, the effects of triaxiality of Earth, and the equatorial second-degree-harmonic components of the core- mantle boundary.

Herring, Thomas A.↗

Neutron-capture Cl-36, Ca-41, Ar-36, and Sm-150 in large chondrites: Evidence for high fluences of thermalized neutrons

We have measured significant concentrations of Cl-36, Ca-41, Ar-36 from decay of Cl-36, and Sm-150 produced from the capture of thermalized neutrons in the large Chico L6 chondrite. Activities of Cl-36 and Ca-41, corrected for a high-energy spallogenic component and a terrestrial age of approximately 50 ka, give average neutron-capture production rates of 208 atoms/min/g-Cl and 1525 atoms/min/kg-Ca, which correspond to thermal neutron (n) fluxes of 6.2 n/sq cm/s and 4.3 n/sq cm/s, respectively. If sustained for the approximately 65 Ma single-stage, cosmic ray exposure age of Chico, these values correspond to thermal neutron fluences of approximately 1.3 x 10(exp 16) and 0.8 x 10(exp 16) n/sq cm for Cl-36 and Ca-41, respectively. Stepwise temperature extraction of Ar in Chico impact melt shows Ar-36/Ar-38 ratios as large as approximately 9. The correlation of high Ar-36/Ar-38 with high Cl/Ca phases in neutron-irradiated Chico indicates that the excess Ar-36 above that expected from spallation is due to decay of neutron-produced Cl-36. Excess Ar-36 in Chico requires a thermal neutron fluence of 0.9-1.7 x 10(exp 16) n/sq cm. Decreases in Sm-149/Sm-152 due to neutron-capture by Sm-149 correlate with increases in Sm-150/Sm-152 for three samples of Chico, and one of the Torino H-chondrite. The 0.08% decrease in Sm-149 shown by Chico corresponds to a neutron fluence of 1.23 x 10(exp 16) n/sq cm. This fluence derived from Sm considers capture of epithermal neutrons and effects of chemical composition on the neutron energy distribution. Excess Ar-36 identified in the Arapahoe, Bruderheim, and Torino chondrites and the Shallowater aubrite suggest exposure to neutron fluences of approximately 0.2-0.2 x 10(exp 16) n/sq cm. Depletion of Sm-149 in Torino and the LEW86010 angrite suggest neutron fluences of 0.8 x 10(exp 16) n/sq cm and 0.25 x 10(exp 16) n/sq cm, respectively. Neutron fluences of approximately 10(exp 16) n/sq cm in Chico are almost as large as those previously observed for some lunar soils. Consideration of exposure ages suggests that the neutron flux in Chico may have been greater than that in many lunar soils.

Bogard, D. D.↗

Direct calculations of waves in fluid flows using high-order compact difference scheme

The solution of the unsteady Euler equations by a sixth-order compact difference scheme combined with a fourth-order Runge-Kutta method is investigated. Closed-form expression for the amplification factors and their corresponding dispersion correlations are obtained by Fourier analysis of the fully discretized, two-dimensional Euler equations. The numerical dissipation, dispersion, and anisotropic effects are assessed. It is found that the Courant-Friedrichs-Lewy (CFL) limit for stable calculations is about 0.8. For a CFL number equal to 0.6, the smallest wavelength which is resolved without numerical damping is about six - eight grid nodes. For phase speeds corresponding to acoustic waves, the corresponding time period is resolved by about 200 - 300 time steps. Three numerical examples of waves in compressible flow are included: (1) sound propagation in a duct with linear shear, (2) linear wave growth in a compressible free shear layer, and (3) vortex pairing in a compressible free shear layer perturbed at two frequencies.

Yu, Sheng-Tao↗

Fast temporal neural learning using teacher forcing

A neural network is trained to output a time dependent target vector defined over a predetermined time interval in response to a time dependent input vector defined over the same time interval by applying corresponding elements of the error vector, or difference between the target vector and the actual neuron output vector, to the inputs of corresponding output neurons of the network as corrective feedback. This feedback decreases the error and quickens the learning process, so that a much smaller number of training cycles are required to complete the learning process. A conventional gradient descent algorithm is employed to update the neural network parameters at the end of the predetermined time interval. The foregoing process is repeated in repetitive cycles until the actual output vector corresponds to the target vector. In the preferred embodiment, as the overall error of the neural network output decreasing during successive training cycles, the portion of the error fed back to the output neurons is decreased accordingly, allowing the network to learn with greater freedom from teacher forcing as the network parameters converge to their optimum values. The invention may also be used to train a neural network with stationary training and target vectors.

Toomarian, Nikzad↗

Multiple Coarse Grid Multigrid Methods for Solving Elliptic Problems

In this paper we describe some classes of multigrid methods for solving large linear systems arising in the solution by finite difference methods of certain boundary value problems involving Poisson's equation on rectangular regions. If parallel computing systems are used, then with standard multigrid methods many of the processors will be idle when one is working at the coarsest grid levels. We describe the use of Multiple Coarse Grid MultiGrid (MCGMG) methods. Here one first constructs a periodic set of equations corresponding to the given system. One then constructs a set of coarse grids such that for each grid corresponding to the grid size h there are four grids corresponding to the grid size 2*h. Multigrid operations such as restriction of residuals and interpolation of corrections are done in parallel at each grid level. For suitable choices of the multigrid operators the MCGMG method is equivalent to the Parallel Superconvergent MultiGrid (PSMG) method of Frederickson and McBryan. The convergence properties of MCGMG methods can be accurately analyzed using spectral methods.

Xiao, Shengyou↗

Experimental Investigation of Unsteady Shock Wave Turbulent Boundary Layer Interactions About a Blunt Fin

A series of experiments were performed to investigate the effects of Mach number variation on the characteristics of the unsteady shock wave/turbulent boundary layer interaction generated by a blunt fin. A single blunt fin hemicylindrical leading edge diameter size was used in all of the experiments which covered the Mach number range from 2.0 to 5.0. The measurements in this investigation included surface flow visualization, static and dynamic pressure measurements, both on centerline and off-centerline of the blunt fin axis. Surface flow visualization and static pressure measurements showed that the spatial extent of the shock wave/turbulent boundary layer interaction increased with increasing Mach number. The maximum static pressure, normalized by the incoming static pressure, measured at the peak location in the separated flow region ahead of the blunt fin was found to increase with increasing Mach number. The mean and standard deviations of the fluctuating pressure signals from the dynamic pressure transducers were found to collapse to self-similar distributions as a function of the distance perpendicular to the separation line. The standard deviation of the pressure signals showed initial peaked distribution, with the maximum standard deviation point corresponding to the location of the separation line at Mach number 3.0 to 5.0. At Mach 2.0 the maximum standard deviation point was found to occur significantly upstream of the separation line. The intermittency distributions of the separation shock wave motion were found to be self-similar profiles for all Mach numbers. The intermittent region length was found to increase with Mach number and decrease with interaction sweepback angle. For Mach numbers 3.0 to 5.0 the separation line was found to correspond to high intermittencies or equivalently to the downstream locus of the separation shock wave motion. The Mach 2.0 tests, however, showed that the intermittent region occurs significantly upstream of the separation line. Power spectral densities measured in the intermittent regions were found to have self-similar frequency distributions when compared as functions of a Strouhal number for all Mach numbers and interaction sweepback angles. The maximum zero-crossing frequencies were found to correspond with the peak frequencies in the power spectra measured in the intermittent region.

Barnhart, Paul J.↗

An Assessment of the Axial and Radial Dilation of a DPIMS Tantalum Cartridge for Space Shuttle Flight Experiments

Ground-based heat treatment tests are planned on an argon gas-filled tantalum cartridge developed as pan of a Diffusion Processes in Molten Semiconductors (DPIMS) experiment conducted on NASA's Space Shuttle. The possibility that the cartridge may creep during testing and touch the furnace walls is of real concern in this program. The present paper discusses the results of calculations performed to evaluate this possibility. Deformation mechanism maps were constructed using literature data in order to identify the creep mechanism dominant under the appropriate stresses and temperatures corresponding to the test conditions. These results showed that power-law creep was dominant when the grain size of the material exceeded 55 gm but Coble creep was the important mechanism below this value of grain size. Finite element analysis was used to analyze the heat treatment tWs assuming a furnace run away condition (which is a worst case scenario) using the appropriate creep parameters corresponding to grain sizes of 1 and 100 gm. Calculations were also conducted to simulate the effect of an initial 3 tilt of the cartridge assembly, the maximum possible tilt angle. The von Mises stress and su-ain distributions were calculated assuming that the cartridge was fixed at one end as it was heated from ambient temperature to 1823 K in 1.42 h, maintained at 1823 K for 9.5 h and then further heated to an over temperature condition of 2028 K in 0.3 h. The inelastic axial and radial displacements of the cartridge walls were evaluated by resolving the von Mises strain along the corresponding directions. These calculations reveal that the maximum axial and radial displacements are expected to be about 2.9 and 0.25 mm, respectively, for both fine and coarse-grained materials at 2028 K. It was determined that these displacements occur during heat-up to temperature and creep of the cartridge is likely to be relatively insignificant irrespective of grain size. Furthermore, with a 3' tilt of the cartridge, the deflection is increased by only 0.39 gm which is negligible. Since the gap between the furnace heating elements and the cartridge is about 7.5 mm and less than the maximum radial dilation of 0.25 mm at 2028 K, it is concluded that the cartridge is unlikely to touch the furnace walls during the experiments.

Raj, S.V.↗

Assessment Study of the State of the Art in Adaptive Control and its Applications to Aircraft Control

Many papers relevant to reconfigurable flight control have appeared over the past fifteen years. In general these have consisted of theoretical issues, simulation experiments, and in some cases, actual flight tests. Results indicate that reconfiguration of flight controls is certainly feasible for a wide class of failures. However many of the proposed procedures although quite attractive, need further analytical and experimental studies for meaningful validation. Many procedures assume the availability of failure detection and identification logic that will supply adequately fast, the dynamics corresponding to the failed aircraft. This in general implies that the failure detection and fault identification logic must have access to all possible anticipated faults and the corresponding dynamical equations of motion. Unless some sort of explicit on line parameter identification is included, the computational demands could possibly be too excessive. This suggests the need for some form of adaptive control, either by itself as the prime procedure for control reconfiguration or in conjunction with the failure detection logic. If explicit or indirect adaptive control is used, then it is important that the identified models be such that the corresponding computed controls deliver adequate performance to the actual aircraft. Unknown changes in trim should be modelled, and parameter identification needs to be adequately insensitive to noise and at the same time capable of tracking abrupt changes. If however, both failure detection and system parameter identification turn out to be too time consuming in an emergency situation, then the concepts of direct adaptive control should be considered. If direct model reference adaptive control is to be used (on a linear model) with stability assurances, then a positive real or passivity condition needs to be satisfied for all possible configurations. This condition is often satisfied with a feedforward compensator around the plant. This compensator must be robustly designed such that the compensated plant satisfies the required positive real conditions over all expected parameter values. Furthermore, with the feedforward only around the plant, a nonzero (but bounded error) will exist in steady state between the plant and model outputs. This error can be removed by placing the compensator also in the reference model. Design of such a compensator should not be too difficult a problem since for flight control it is generally possible to feedback all the system states.

Kaufman, Howard↗

Synchronous Parallel Emulation and Discrete Event Simulation System with Self-Contained Simulation Objects and Active Event Objects

The present invention is embodied in a method of performing object-oriented simulation and a system having inter-connected processor nodes operating in parallel to simulate mutual interactions of a set of discrete simulation objects distributed among the nodes as a sequence of discrete events changing state variables of respective simulation objects so as to generate new event-defining messages addressed to respective ones of the nodes. The object-oriented simulation is performed at each one of the nodes by assigning passive self-contained simulation objects to each one of the nodes, responding to messages received at one node by generating corresponding active event objects having user-defined inherent capabilities and individual time stamps and corresponding to respective events affecting one of the passive self-contained simulation objects of the one node, restricting the respective passive self-contained simulation objects to only providing and receiving information from die respective active event objects, requesting information and changing variables within a passive self-contained simulation object by the active event object, and producing corresponding messages specifying events resulting therefrom by the active event objects.

Steinman, Jeffrey S.↗

Hydrodynamic Instability and Thermal Coupling in a Dynamic Model of Liquid-Propellant Combustion

For liquid-propellant combustion, the Landau/Levich hydrodynamic models have been combined and extended to account for a dynamic dependence of the burning rate on the local pressure and temperature fields. Analysis of these extended models is greatly facilitated by exploiting the realistic smallness of the gas-to-liquid density ratio rho. Neglecting thermal coupling effects, an asymptotic expression was then derived for the cellular stability boundary A(sub p)(k) where A(sub p) is the pressure sensitivity of the burning rate and k is the disturbance wavenumber. The results explicitly indicate the stabilizing effects of gravity on long-wave disturbances, and those of viscosity and surface tension on short-wave perturbations, and the instability associated with intermediate wavenumbers for critical negative values of A(sub p). In the limit of weak gravity, hydrodynamic instability in liquid-propellant combustion becomes a long-wave, instability phenomenon, whereas at normal gravity, this instability is first manifested through O(1) wavenumbers. In addition, surface tension and viscosity (both liquid and gas) each produce comparable effects in the large-wavenumber regime, thereby providing important modifications to the previous analyses in which one or more of these effects was neglected. For A(sub p)= O, the Landau/Levich results are recovered in appropriate limiting cases, although this typically corresponds to a hydrodynamically unstable parameter regime for p << 1. In addition to the classical cellular form of hydrodynamic stability, there exists a pulsating form corresponding to the loss of stability of steady, planar burning to time-dependent perturbations. This occurs for negative values of the parameter A(sub p), and is thus absent from the original Landau/Levich models. In the extended model, however, there exists a stable band of negative pressure sensitivities bounded above by the Landau type of instability, and below by this pulsating form of hydrodynamic instability. Indeed, nonsteady modes of combustion have been observed at low pressures in hydroxylammonium nitrate (HAN)-based liquid propellants, which often exhibit negative pressure sensitivities. While nonsteady combustion may correspond to secondary and higher-order bifurcations above the cellular boundary, it may also be a manifestation of this pulsating type of hydrodynamic instability. In the present work, a nonzero temperature sensitivity is incorporated into our previous asymptotic analyses. This entails a coupling of the energy equation to the previous purely hydrodynamic problem, and leads to a significant modification of the pulsating boundary such that, for sufficiently large values of the temperature-sensitivity parameter, liquid-propellant combustion can become intrinsically unstable to this alternative form of hydrodynamic instability. For simplicity, further attention is confined here to the inviscid version of the problem since, despite the fact that viscous and surface-tension effects are comparable, the qualitative nature of the cellular boundary remains preserved in the zero-viscosity limit, as does the existence of the pulsating boundary. The mathematical model adopts the classical assumption that there is no distributed reaction in either the liquid or gas phases, but now the reaction sheet, representing either a pyrolysis reaction or an exothermic decomposition at the liquid/gas interface, is assumed to depend on local conditions there.

Margolis, S. B.↗

C-Band Backscatter Measurements of Winter Sea-Ice in the Weddell Sea, Antarctica

During the 1992 Winter Weddell Gyre Study, a C-band scatterometer was used from the German ice-breaker R/V Polarstern to obtain detailed shipborne measurement scans of Antarctic sea-ice. The frequency-modulated continuous-wave (FM-CW) radar operated at 4-3 GHz and acquired like- (VV) and cross polarization (HV) data at a variety of incidence angles (10-75 deg). Calibrated backscatter data were recorded for several ice types as the icebreaker crossed the Weddell Sea and detailed measurements were made of corresponding snow and sea-ice characteristics at each measurement site, together with meteorological information, radiation budget and oceanographic data. The primary scattering contributions under cold winter conditions arise from the air/snow and snow/ice interfaces. Observations indicate so e similarities with Arctic sea-ice scattering signatures, although the main difference is generally lower mean backscattering coefficients in the Weddell Sea. This is due to the younger mean ice age and thickness, and correspondingly higher mean salinities. In particular, smooth white ice found in 1992 in divergent areas within the Weddell Gyre ice pack was generally extremely smooth and undeformed. Comparisons of field scatterometer data with calibrated 20-26 deg incidence ERS-1 radar image data show close correspondence, and indicate that rough Antarctic first-year and older second-year ice forms do not produce as distinctively different scattering signatures as observed in the Arctic. Thick deformed first-year and second-year ice on the other hand are clearly discriminated from younger undeformed ice. thereby allowing successful separation of thick and thin ice. Time-series data also indicate that C-band is sensitive to changes in snow and ice conditions resulting from atmospheric and oceanographic forcing and the local heat flux environment. Variations of several dB in 45 deg incidence backscatter occur in response to a combination of thermally-regulated parameters including sea-ice brine volume, snow and ice complex dielectric properties, and snow physical properties.

Drinkwater, M. R.↗

Determining the Molecular Growth Mechanisms of Protein Crystal Faces by Atomic Force Microscopy

A high resolution atomic force microscopy (AFM) study had shown that the molecular packing on the tetragonal lysozyme (110) face corresponded to only one of two possible packing arrangements, suggesting that growth layers on this face were of bimolecular height. Theoretical analyses of the packing also indicated that growth of this face should proceed by the addition of growth units of at least tetramer size corresponding to the 43 helices in the crystal. In this study an AFM linescan technique was devised to measure the dimensions of individual growth units on protein crystal faces as they were being incorporated into the lattice. Images of individual growth events on the (110) face of tetragonal lysozyme crystals were observed, shown by jump discontinuities in the growth step in the linescan images as shown in the figure. The growth unit dimension in the scanned direction was obtained from these images. A large number of scans in two directions on the (110) face were performed and the distribution of lysozyme growth unit sizes were obtained. A variety of unit sizes corresponding to 43 helices, were shown to participate in the growth process, with the 43 tetramer being the minimum observed size. This technique represents a new application for AFM allowing time resolved studies of molecular process to be carried out.

Nadarajah, Arunan↗

Microphysical Retrievals Over Stratiform Rain Using Measurements from an Airborne Dual-Wavelength Radar-Radiometer

The need to understand the complementarity of the radar and radiometer is important not only to the Tropical Rain Measuring Mission (TRMM) program but to a growing number of multi-instrumented airborne experiment that combine single or dual-frequency radars with multichannel radiometers. The method of analysis used in this study begins with the derivation of dual-wavelength radar equations for the estimation of a two-parameter drop size distribution (DSD). Defining a "storm model" as the set of parameters that characterize snow density, cloud water, water vapor, and features of the melting layer, then to each storm model there will usually correspond a set of range-profiled drop size distributions that are approximate solutions of the radar equations. To test these solutions, a radiative transfer model is used to compute the brightness temperatures for the radiometric frequencies of interest. A storm model or class of storm models is considered optimum if it provides the best reproduction of the radar and radiometer measurements. Tests of the method are made for stratiform rain using simulated storm models as well as measured airborne data. Preliminary results show that the best correspondence between the measured and estimated radar profiles usually can be obtained by using a moderate snow density (0.1-0.2 g/cu cm), the Maxwell-Garnett mixing formula for partially melted hydrometeors (water matrix with snow inclusions), and low to moderate values of the integrated cloud liquid water (less than 1 kg/sq m). The storm-model parameters that yield the best reproductions of the measured radar reflectivity factors also provide brightness temperatures at 10 GHz that agree well with the measurements. On the other hand, the correspondence between the measured and modeled values usually worsens in going to the higher frequency channels at 19 and 34 GHz. In searching for possible reasons for the discrepancies, It is found that changes in the DSD parameter Mu, the radar constants, or the path-integrated attenuation can affect the high frequency channels significantly. In particular, parameters that cause only modest increases in the median mass diameter of the snow, and which have a minor effect on the radar returns or the low frequency brightness temperature, can produce a strong cooling of the 34 GHz brightness temperature.

Meneghini, Robert↗

Lunar Meteorites and Implications for Compositional Remote Sensing of the Lunar Surface

Lunar meteorites (LMs) are rocks found on Earth that were ejected from the Moon by impact of an asteroidal meteoroid. Three factors make the LMs important to remote-sensing studies: (1) Most are breccias composed of regolith or fragmental material; (2) all are rocks that resided (or breccias composed of material that resided) in the upper few meters of the Moon prior to launch and (3) most apparently come from areas distant from the Apollo sites. How Many Lunar Locations? At this writing (June 1999), there are 18 known lunar meteorite specimens. When unambiguous cases of terrestrial pairing are considered, the number of actual LMs reduces to 13. (Terrestrial pairing is when a single piece of lunar rock entered Earth's atmosphere, but multiple fragments were produced because the meteoroid broke apart on entry, upon hitting the ground or ice, or while being transported through the ice.) We have no reason to believe that LMs preferentially derive from any specific region(s) of the Moon; i.e., we believe that they are samples from random locations. However, we do not know how many different locations are represented by the LMs; mathematically, it could be as few as 1 or as many as 13. The actual maximum is < 13 because in some cases a single impact appears to have yielded more than one LM. Yamato 793169 and Asuka 881757 are considered "source-crater paired" or "launch paired" because they are compositionally and petrographically similar to each other and distinct from the others, and both have similar cosmic-ray exposure (CRE) histories. The same can be said of QUE 94281 and Y 793274. Thus the 13 meteorites probably represent a maximum of 11 locations on the Moon. The minimum number of likely source craters is debated and in flux as new data for different isotopic systems are obtained. Conservatively, considering CRE data only, a minimum of about 5 impacts is required. Compositional and petrographic data offer only probabilistic constraints. An extreme, but not unreasonable viewpoint, is that such data offer no constraint. For example, if one were to cut up the Apollo 17 landing site (which was selected for its diversity) into softball-sized pieces, some of those pieces (e.g., sample 70135) would be crystalline mare basalts like Y 793169 whereas others (e.g., sample 73131 would be feldspathic regolith breccias like MAC 88104/ 88105. However, nature is not so devious. Warren argues that LMs come from craters of only a few kilometers in diameter. If so, even though CRE data allow, for example, that ALHA 81005 and Y 791197) were launched simultaneously from the same crater, the probability is nevertheless low because the two meteorites are compositionally and mineralogically distinct. Thus, within the allowed range (5-11) for the number of locations represented by the LMs, values at the high end of the range are probably more likely. Mare Meteorites: Three LMs consist almost entirely of mare basalt. Two, Y 793169 and Asuka 881757, are unbrecciated, low-Ti, crystalline rocks that are compositionally and mineralogically similar (but not identical) to each other; they probably derive from a single lunar-mare location. The third, EET 87521/96008, is a fragmental breccia consisting predominantly of VLT mare basalt. Thus, these LMs probably represent only two lunar mare locations. The basaltic LMs have mineral and bulk compositions distinct from Apollo mare basalts. The petrography of Calcalong Creek has not been described in detail, but compositionally it is unique in that it corresponds to a mixture (breccia) of about one-half feldspathic material (i.e., the mean composition of the feldspathic lunar meteorites, below), one-fourth KREEP norite, one-fourth VLT mare basalt (like EET 87521), and 1% CI chondrite. With 4 micro g/g Th and correspondingly high concentrations of other incompatible elements, it is the only lunar meteorite that is likely to have come from within the Procellarum KREEP Terrane (PKT). Yamato 793274 and QUE 94281 are together distinct in being fragmental breccias containing subequal parts of feldspathic highland material and VLT mare basalt. Jolliff et al. estimate a mare to highland ratio of 54:46 for QUE 94281 and 62:38 for Y 793274; this difference is well within the range observed for soils collected only centimeters apart (in cores) at interface site like Apollo 15 and 17 [11]. Although the two meteorites were found on opposite sides of Antarctica, they are probably launch-paired. The strongest evidence is that the pyroclastic glass spherules that occur in both are of two compositional groups and the two groups are essentially the same in both meteorites. Yamato 791197 is nominally a feldspathic lunar meteorite (below), but among FLMs, it probably contains the highest abundance of clasts and glasses of mare derivation. As a consequence, its composition is at the high-Fe, low-Mg end of the range for FLMs and is not included in the FLM average of Table 1. Its composition is consistent with about 10% mare-derived material. Similarly, the two small (Y 82) pieces of Y 82192/82193186032 are more mafic than the large (Y 86) piece, probably as a result of about 7% mare-derived material. All Apollo missions went to areas in or near the PKT, and, consequently, all Apollo regolith samples are contaminated with Th-rich material from the PKT. At the nominally "typical" highland site, Apollo 16, about 30% of the regolith (<1-mm fines) is Th-rich ejecta from the Imbrium impact and about 6% is mare material probably derived from mare basins. Thus Apollo 16 regolith is not typical of the highlands. Among Apollo rocks, the compositions of the FLMs correspond most closely to the feldspathic granulitic breccias of Apollo 16 and 17. (Additional information is contained in original)

Korotev, R. L.↗

Shipboard Sunphotometer Measurements of Aerosol Optical Depth Spectra and Columnar Water Vapor During ACE-2 and Comparison with Selected Land, Ship, Aircraft, and Satellite Measurements

Analyses of aerosol optical depth (AOD) and colurnmn water vapor (CWV) measurements acquired with NASA Ames Research Center's 6-channel Airborne Tracking Sunphotometer (AATS-6) operated aboard the R/V Professor Vodyanitskiy during the 2nd Aerosol Characterization Experiment (ACE-2) are discussed. Data are compared with various in situ and remote measurements for selected cases. The focus is on 10 July, when the Pelican airplane flew within 70 km of the ship near the time of a NOAA-14/AVHRR satellite overpass and AOD measurements with the 14-channel Ames Airborne Tracking Sunphotometer (AATS-14) above the marine boundary layer (MBL) permitted calculation of AOD within the MBL from the AATS-6 measurements. A detailed column closure test is performed for MBL AOD on 10 July by comparing the AATS-6 MBL AODs with corresponding values calculated by combining shipboard particle size distribution measurements with models of hygroscopic growth and radiosonde humidity profiles (plus assumptions on the vertical profile of the dry particle size distribution and composition). Large differences (30-80% in the mid-visible) between measured and reconstructed AODs are obtained, in large part because of the high sensitivity of the closure methodology to hygroscopic growth models, which vary considerably and have not been validated over the necessary range of particle size/composition distributions. The wavelength dependence of AATS-6 AODs is compared with the corresponding dependence of aerosol extinction calculated from shipboard measurements of aerosol size distribution and of total scattering mearured by a shipboard integrating nephelometer for several days. Results are highly variable, illustrating further the great difficulty of deriving column values from point measurements. AATS-6 CWV values are shown to agree well with corresponding values derived from radiosonde measurements during 8 soundings on 7 days and also with values calculated from measurements taken on 10 July with the AATS-14 and the University of Washington Passive Humidigraph aboard the Pelican.

Livingston, John M.↗

The Relationship Between Magnetotail Dynamics and Substorm Onset Longitudes Determined from Spacecraft Images

Geotail plasma and magnetic field observations of plasmoids between 25 and 30 Re have been compared to Polar UVI observations of auroral brightenings. Plasmoids almost always corresponded to brightenings but the brightenings were sometimes weak and spatially limited and did not always grow to a global substorm. Even a case where a plasmoid event occurred with fast post-plasmoid flow corresponded to a weak brightening but no substorm. Some brightenings did not correspond to plasmoids, but these events were observed away from the longitude of Geotail. The plasmoids are observed 0-2 min after the brightenings in most cases. It seems likely that formation of a near-Earth neutral line causes each brightening in the polar ionosphere, but these formations do not have a one-to-one relationship with a substorm onset. What causes development of the full, large-scale substorm remains an open question.

Ieda, A.↗

Pulsating Hydrodynamic Instability and Thermal Coupling in an Extended Landau/Levich Model of Liquid-Propellant Combustion

Hydrodynamic (Landau) instability in combustion is typically associated with the onset of wrinkling of a flame surface, corresponding to the formation of steady cellular structures as the stability threshold is crossed. In the context of liquid-propellant combustion, such instability has recently been shown to occur for critical values of the pressure sensitivity of the burning rate and the disturbance wavenumber, significantly generalizing previous classical results for this problem that assumed a constant normal burning rate. Additionally, however, a pulsating form of hydrodynamic instability has been shown to occur as well, corresponding to the onset of temporal oscillations in the location of the liquid/gas interface. In the present work, we consider the realistic influence of a nonzero temperature sensitivity in the local burning rate on both types of stability thresholds. It is found that for sufficiently small values of this parameter, there exists a stable range of pressure sensitivities for steady, planar burning such that the classical cellular form of hydrodynamic instability and the more recent pulsating form of hydrodynamic instability can each occur as the corresponding stability threshold is crossed. For larger thermal sensitivities, however, the pulsating stability boundary evolves into a C-shaped curve in the disturbance-wavenumber/ pressure-sensitivity plane, indicating loss of stability to pulsating perturbations for all sufficiently large disturbance wavelengths. It is thus concluded, based on characteristic parameter values, that an equally likely form of hydrodynamic instability in liquid-propellant combustion is of a nonsteady, long-wave nature, distinct from the steady, cellular form originally predicted by Landau.

Margolis, S. B.↗

The Role of Clear Sky Identification in the Study of Cloud Radiative Effects: Combine Analysis from ISCCP and the Scanner of Radiation Budget (ScaRaB)

Since the effect of clouds on the earth's radiation balance is often estimated as the difference of net radiative fluxes at the top of the atmosphere between all situations and monthly averaged clear sky situations of the same regions, a reliable identification of clear sky is important for the study of cloud radiative effects. The Scanner for Radiation Balance (ScaRaB) radiometer on board the Russian Meteor-3/7 satellite provided earth radiation budget observations from March 1994 to February 1995 with two ERBE-Re broad-band longwave and shortwave channels. Two narrow-band channels, in the infrared atmospheric window and in the visible band, have been added to the ScaRaB instrument to improve the cloud scene identification. The International Satellite Cloud Climatology Project (ISCCP) method for cloud detection and determination of cloud and surface properties uses the same narrow-band channels as ScaRaB, but is employed to a collection of measurements at a better spatial resolution of about 5 km. By applying the original ISCCP algorithms to the ScaRaB data, the clear sky frequency is about 5% lower than the one over quasi-simultaneous original ISCCP data, an indication that the ISCCP cloud detection is quite stable. However, one would expect an about 10 to 20% smaller clear sky occurrence over the larger ScaRaB pixels. Adapting the ISCCP algorithms to the reduced spatial resolution of 60 km and to the different time sampling of the ScaRaB data leads therefore to a reduction of a residual cloud contamination. A sensitivity study with time-space collocated ScaRaB and original ISCCP data at a spatial resolution of 1deg longitude x 1deg latitude shows that the effect of clear sky identification method plays a higher role on the clear sky frequency and therefore on the statistics than on the zonal mean values of the clear sky fluxes. Nevertheless, the zonal outgoing longwave fluxes corresponding to ERBE clear sky are in general about 2 to 10 W/sq m higher than those obtained from the ScaRaB adapted ISCCP clear sky identifications. The latter are close to (about 1 W/sq m higher) fluxes corresponding to clear sky regions from original ISCCP data, whereas ScaRaB clear sky LW fluxes obtained with the original ISCCP identification lie about 1 to 2 W/sq m below. Especially in the tropics where water vapor abundance is high, the ERBE clear sky LW fluxes seem to be systematically overestimated by about 4 W/sq m, and SW fluxes are lower by about 5 to 10 W/sq m. However, the uncertainty in the analysis of monthly mean zonal cloud radiative effects is also produced by the low frequency of clear sky occurrence, illustrated when averaging over pixels or even over regions of 4deg longitude x 5deg latitude, corresponding to the spatial resolution of General Circulation Models. The systematic bias in the clear sky fluxes is not reflected in the zonal cloud radiative effects, because the clear sky regions selected by the different algorithms can occur in different geographic regions with different cloud properties.

Rossow, W. B.↗