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,351 records · Page 75

"CHON" particles: The interstellar component of cometary dust

Interstellar dust is characterized by strong absorption in the ultraviolet and the mid-IR. Current models of interstellar dust are based on three chemically distinct components: a form of carbon (usually graphite), a silicate, and a blend of polycyclic aromatic hydrocarbons or other carbonaceous material. Previous work using effective medium theories to understand the optical properties of cometary dust suggested that an amalgam of materials could reproduce the observed interstellar and cometary dust features. Recently, Lawler and Brownlee (1992) re-analyzed the PIA and PUMA-1 data sets from the Giotto flyby of P/Halley and discovered that the so-called "CHON" particles were actually composed of a blend of carbon-bearing and silicon-bearing materials. Based on effective medium theories, the absorption spectrum of such a material would display the spectral features of each of the components - strong UV absorption from the carbonaceous component and strong absorption in the IR from the silicate component. To test this idea, vapor-deposited samples were created using two different deposition techniques: sputtering with an argon RF magnetron and deposition from an argon plasma torch. Two different compositions were tested: a blend of graphite and silica in a 7:1 ratio and an amalgam of materials whose approximate composition matches the "CHON"-silicate abundances for the uncompressed PIA data set of Lawler and Brownlee: graphite, iron oxide, magnesium oxide, ammonium sulfate, calcium carbonate, and silica in mass ratios of 6:4.3:4:2.2:1:9. The samples were finely ground and pressed into 2" diameter disks using a 40 ton press. In all, four different experiments were performed: one with each of the compositions (C:SiO and "CHON") in both the RF magnetron and the plasma torch chambers. The RF magnetron created a uniform dark thin film on the substrate surface, and the plasma torch created a coating of small (<100 micron) diameter grey particles. The spectra of all four samples show a strong, broad absorption feature at around 220 nm as well as a strong but narrower absorption peak near 10 microns. The RF magnetron sputtered samples showed some sub-structure in the UV, and the peak of the absorption was shifted toward longer wavelengths. The UV absorption in the plasma torch deposited samples have no sub-structure, and the peak absorption is very near 220 nm. Strong absorption near 9 microns is seen in the spectra from both sample preparation techniques, and is consistent with the IR spectra of some terrestrial silicates. Other features, particularly at 6.2 and 8.6 microns, are seen in the interstellar medium. A strong feature near 2 microns is due to absorbed water in the sample. Based on the results of these experiments, there is evidence that a material with a composition similar to that detected in "CHON" particles in the coma of P/Halley have a spectral signature which reproduces the main absorption features of interstellar dust. This suggests that the "CHON" particles could be the interstellar component of cometary dust.

Lien, David J.↗

Deterministic Stress Modeling of Hot Gas Segregation in a Turbine

Simulation of unsteady viscous turbomachinery flowfields is presently impractical as a design tool due to the long run times required. Designers rely predominantly on steady-state simulations, but these simulations do not account for some of the important unsteady flow physics. Unsteady flow effects can be modeled as source terms in the steady flow equations. These source terms, referred to as Lumped Deterministic Stresses (LDS), can be used to drive steady flow solution procedures to reproduce the time-average of an unsteady flow solution. The goal of this work is to investigate the feasibility of using inviscid lumped deterministic stresses to model unsteady combustion hot streak migration effects on the turbine blade tip and outer air seal heat loads using a steady computational approach. The LDS model is obtained from an unsteady inviscid calculation. The LDS model is then used with a steady viscous computation to simulate the time-averaged viscous solution. Both two-dimensional and three-dimensional applications are examined. The inviscid LDS model produces good results for the two-dimensional case and requires less than 10% of the CPU time of the unsteady viscous run. For the three-dimensional case, the LDS model does a good job of reproducing the time-averaged viscous temperature migration and separation as well as heat load on the outer air seal at a CPU cost that is 25% of that of an unsteady viscous computation.

Busby, Judy↗

A Modified Linear-Mixing Method for Calculating Atmospheric Path Radiances of Aerosol Mixtures

The top-of-atmosphere (TOA) path radiance generated by an aerosol mixture can be synthesized by linearly adding the contributions of the individual aerosol components, weighted by their fractional optical depths. The method, known as linear mixing, is exact in the single-scattering limit. When multiple scattering is significant, the method reproduces the atmospheric path radiance of the mixture with less than 3% errors for weakly absorbing aerosols up to optical thickness of 0.5. However, when strongly absorbing aerosols are included in the mixture, the errors are much larger. This is due to neglecting the effect of multiple interactions between the aerosol components, especially when the values of the single-scattering albedos of these components are so different that the parameter epsilon = (Sigma)f(sub i) absolute value of bar omega(sub i) - bar omega(sub mix)/bar omega(sub i), is larger than approx. 0.1, where bar omega(sub i) and f(sub i) are the single-scattering albedo and the fractional abundance of the i th component, and bar omega(sub i) is the effective single-scattering albedo of the mixture. We describe an empirical, modified linear-mixing method which effectively accounts for the multiple interactions between aerosol components. The modified and standard methods are identical when epsilon = 0.0 and give similar results when epsilon less than or equal to 0.05. For optical depths larger than approx. 0.5, or when epsilon greater than 0.05, only the modified method can reproduce the radiances within 5% error for common aerosol types up to optical thickness of 2.0. Because this method facilitates efficient and accurate atmospheric path radiance calculations for mixtures of a wide variety of aerosol types, it will be used as part of the aerosol retrieval methodology for the Earth Observing System (EOS) multiangle imaging spectroradiometer (MISR), scheduled for launch into polar orbit in 1998.

Abdou, W. A.↗

A Modified Linear-Mixing Method for Calculating Atmospheric Path Radiances of Aerosol Mixtures

The top-of-atmosphere (TOA) path radiance generated by an aerosol mixture can be synthesized by linearly adding the contributions of the individual aerosol components, weighted by their fractional optical depths. The method, known as linear mixing, is exact in the single-scattering limit. When multiple scattering is significant, the method reproduces the atmospheric path radiance of the mixture with less than 3% errors for weakly absorbing aerosols up to optical thickness of 0.5. However, when strongly absorbing aerosols are included in the mixture, the errors are much larger. This is due to neglecting the effect of multiple interactions between the aerosol components, especially when the values of the single-scattering albedos of these components are so different that the parameter e = the sum of f(sub i)[(bar)omega(sub i) - (bar)omega(sub mix)]/(bar)omega(sub i) is larger than approximately 0.1, where (bar)omega(sub i)and f(sub i) are the single-scattering albedo and the fractional abundance of the ith component, and (bar)omega(sub mix) is the effective single-scattering albedo of the Mixture. We describe an empirical, modified linear-mixing method which effectively accounts for the multiple interactions between aerosol components. The modified and standard methods are identical when epsilon = 0.0 and give similar results when epsilon is less than or equal to 0.05. For optical depths larger than approximately 0.5, or when epsilon is greater than 0.05, only the modified method can reproduce the radiances within 5% error for common aerosol types up to optical thickness of 2.0. Because this method facilitates efficient and accurate atmospheric path radiance calculations for mixtures of a wide variety of aerosol types, it will be used as part of the aerosol retrieval methodology for the Earth Observing System (EOS) multiangle imaging spectroradiometer (MISR), scheduled for launch into polar orbit in 1998.

Abdou, W. A.↗

Experiment 13: The Study of Dopant Segregation Behavior During the Growth of GaAs in Microgravity on USML-2

An investigation into the segregation behavior of selenium doped gallium arsenide (Se/GaAs) during directional solidification in the microgravity environment was conducted using the Crystal Growth Furnace (CGF) aboard the second United States Microgravity Laboratory (USML-2). Two crystals were successfully processed on USML-2, which lasted from October 20 to November 7, 1995. The first sample was processed for 67 hours, 45 minutes (MET 5/04:53:45-8/00:23:50) and included 19 hours of growth at 0.5 microns/sec which yielded 3.42 cm of sample length, and 5 hours of growth at 1.5 microns/sec which yielded 2.7 cm of sample. During the second experiment, the furnace temperature was adjusted to move the melt-solid interface position towards the hot end of the furnace. The second sample was processed for 50 hours, 10 minutes (MET 8/18:48:49-10/21:58:54) and included 11 hours of growth at 0.5 microns/sec which yielded 1.98 cm of sample, and 1 hour, 25 minutes of growth at 5.0 microns/sec which yielded 2.6 cm of sample. This sample provides an order of magnitude change in growth rate and reproduces one of the growth rates used during USML-1. In contrast to the results from USML-1, no voids were present in either crystal grown on USML-2. The absence of voids in either sample indicates that growth rate changes alone were not responsible for the formation of voids found in the crystals grown on USML-1. Sections of the ground-based and flight crystals grown on USML-2 were cut and polished. All of the interface demarcation lines expected from the current pulse interface demarcation (CPID) system have been identified. These measurements have been analyzed for interface positions, interface shapes, and growth rates. Using a newly developed technique, based on experimental and numerical results, the seeding interface reproducibility from run to run was <= 2.5 mm. The seeding interface position could be controllably moved, with respect to the furnace zones, by adjusting the control set points of the heating zones. The interface shapes flattened slightly as the interface position moved closer to the hot zone but was always an unfavorable concave into the solid shape. The growth rate was found to equal the furnace translation rate, after a 2 -hour transient, for growth rates <= 1.0 microns/sec. Segregation measurements for the ground-based crystals are indicative of complete mixing behavior, as expected. Segregation measurements of the flight crystals are still in progress.

Matthiesen, David H.↗

Boundary Layer Formation in the Magnetotail: Geotail Observations and Comparisons with a Global MHD Simulation

We present Geotail plasma and field observations from the middle magnetotail near X(sub GSE) = -46 R(sub E) for the time period 1400 to 1800 UT on December 14, 1994. During that period, the Wind satellite monitored the solar wind plasma and interplanetary magnetic field (IMF) upstream of the bow shock. The IMF was northward and the plasma parameters near average. Geotail observed slow tailward flows and a northward field. The plasma and field parameters indicate that Geotail is either in the plasma sheet or in a boundary layer. We used the Wind solar wind plasma and IMF data as input for a global simulation of that time interval. Comparison of the simulation results with the observational data show very good overall agreement of the magnitudes of the plasma and field parameters. In particular, the simulation reproduces the slow tailward flows and northward field found at Geotail. Small scale temporal, variations are less well reproduced. The simulation shows the formation of a broad boundary layer (which we call tail flank boundary layer, TFBL) that consists of closed flux which is formed by magnetic magnetic reconnection of IMF and lobe field lines. The simulation results indicate that Geotail is located very close to the TFBL and may have entered the TFBL proper. We show that the TFBL plays an important role in energy transport from the solar wind into the magnetosphere during northward IMF conditions.

Raeder, J.↗

Modeling boot-shaped coronal holes using SOHO-MDI magnetic measurements

The boot-shaped coronal hole observed between the 22 and 30 August 1996, reproduced by rigidly rotating the data computed on 27 August with the rotation rate near the solar equator, is discussed. The data were acquired by the Michelson Doppler imager (MDI) onboard the Solar and Heliospheric Observatory (SOHO). The results suggest that change in size and shape of the holes observed during the period is caused by the projection effect. It was found that the whole solar observational synoptic chart may be used to approximately reproduce boot-shaped holes in solar minimum.

Zhao, X. P.↗

Fitting a 3-D analytic model of the coronal mass ejection to observations

The application of an analytic magnetohydrodynamic model is presented to observations of the time-dependent explusion of 3D coronal mass ejections (CMEs) out of the solar corona. This model relates the white-light appearance of the CME to its internal magnetic field, which takes the form of a closed bubble, filled with a partly anchored, twisted magnetic flux rope and embedded in an otherwise open background field. The density distribution frozen into the expanding CME expanding field is fully 3D, and can be integrated along the line of sight to reproduce observations of scattered white light. The model is able to reproduce the three conspicuous features often associated with CMEs as observed with white-light coronagraphs: a surrounding high-density region, an internal low-density cavity, and a high-density core. The model also describes the self-similar radial expansion of these structures. By varying the model parameters, the model can be fitted directly to observations of CMEs. It is shown how the model can quantitatively match the polarized brightness contrast of a dark cavity emerging through the lower corona as observed by the HAO Mauna Loa K-coronameter to within the noise level of the data.

Gibson, S. E.↗

Variability in shell models of GRBs

Many cosmological models of gamma-ray bursts (GRBs) assume that a single relativistic shell carries kinetic energy away from the source and later converts it into gamma rays, perhaps by interactions with the interstellar medium or by internal shocks within the shell. Although such models are able to reproduce general trends in GRB time histories, it is difficult to reproduce the high degree of variability often seen in GRBs. The authors investigate methods of achieving this variability using a simplified external shock model. Since the model emphasizes geometric and statistical considerations, rather than the detailed physics of the shell, it is applicable to any theory that relies on relativistic shells. They find that the variability in GRBs gives strong clues to the efficiency with which the shell converts its kinetic energy into gamma rays.

Sumner, M. C.↗

On the Relative Importance of Convection and Temperature on the Behavior of the Ionosphere in North American during January 6-12, 1997

Measurements from a network of digisondes and an incoherent scatter radar In Eastern North American For January 6-12, 1997 have been compared with the Field Line Interhemispheric Plasma (FLIP) model which now includes the effects of electric field convective. With the exception of Bermuda, the model reproduces the daytime electron density very well most of the time. As is typical behavior for winter solar minimum on magnetically undisturbed nights, the measurements at Millstone Hill show high electron temperatures before midnight followed by a rapid decay, which is accompanied by a pronounced density enhancement in the early morning hours. The FLIP model reproduces the nighttime density enhancement well, provided the model is constrained to follow the topside electron temperature and the flux tube is full. Similar density enhancements are seen at Goose Bay, Wallops Island and Bermuda. However, the peak height variation and auroral images indicate the density enhancements at Goose Bay are most likely due to particle precipitation. Contrary to previously published work we find that the nighttime density variation at Millstone Hill is driven by the temperature behavior and not the other way around. Thus, in both the data and model, the overall nighttime density is lowered and the enhancement does not occur if the temperature remains high all night. Our calculations show that convections of plasma from higher magnetic latitudes does not cause the observed density maximum but it may enhance the density maximum if over-full flux tubes are convected over the station. On the other had, convection of flux tubes with high temperatures and depleted densities may prevent the density maximum from occurring. Despite the success in modeling the nighttime density enhancements, there remain two unresolved problems. First, the measured density decays much faster than the modeled density near sunset at Millstone Hill and Goose Bay though not at lower latitude stations. Second, we cannot fully explain the large temperatures before midnight nor the sudden decay near midnight.

Richards, P. G.↗

Streamwise Vorticity Generation in Laminar and Turbulent Jets

Complex streamwise vorticity fields are observed in the evolution of non-circular jets. Generation mechanisms are investigated via Reynolds-averaged (RANS), large-eddy (LES) and direct numerical (DNS) simulations of laminar and turbulent rectangular jets. Complex vortex interactions are found in DNS of laminar jets, but axis-switching is observed only when a single instability mode is present in the incoming mixing layer. With several modes present, the structures are not coherent and no axis-switching occurs, RANS computations also produce no axis-switching. On the other hand, LES of high Reynolds number turbulent jets produce axis-switching even for cases with several instability modes in the mixing layer. Analysis of the source terms of the mean streamwise vorticity equation through post-processing of the instantaneous results shows that, complex interactions of gradients of the normal and shear Reynolds stresses are responsible for the generation of streamwise vorticity which leads to axis-switching. RANS computations confirm these results. k - epsilon turbulence model computations fail to reproduce the phenomenon, whereas algebraic Reynolds stress model (ASM) computations, in which the secondary normal and shear stresses are computed explicitly, succeeded in reproducing the phenomenon accurately.

Demuren, Aodeji O.↗

Search for Unique Organic Biomarkers in ALH84001

Four goals were outlined for this project. These were: [1] to reproduce the measurement of polycyclic aromatic hydrocarbons (PAHS) in ALH84001 with both a higher spatial resolution and sensitivity than has been previously reported; [2] to extend such measurements to include other members of the Martian SNC (Shergotties, Nahklites, and Chassigny) meteorite clan, in particular the Antarctic Martian meteorite EETA79001; [3] to address issues of potential organic contamination, because at present very little is known about the effect of terrestrial weathering in the Antarctic environment as it pertains to perturbing an indigenous organic distribution within a meteoritic matrix; and [4] to diversify the range of organic compounds studied to include species that can serve as unique biological markers - "molecular fossils" - derived from once living organisms. In order to achieve this, three specific goals were outlined for the funding period 06/01/97 to 02/28/98. They were: [1] to investigate the effects of terrestrial weathering and organic contamination of meteoritic samples collected from Antarctica; [2] to reproduce and extend upon the measurements of PAHs in ALH84001 with the aim of establishing or refuting the indigeneity of these species; and [3] to extend the analysis of organic compounds in ALH84001 and EETA79001 to address compounds that are considered to be more biologically relevant than PAHS. All three were successfully accomplished, as detailed in the previous performance report. In brief, however, the results achieved were to establish that the PAHs found in ALH84001 were indigenous and not due to contamination, and to determine that a novel and sensitive technique in meteoritic work, capillary zone electrophoresis (CE), could indeed detect amino acids, a potential class of biomarker.

Zare, Richard N.↗

Properties of QBO and SAO Generated by Gravity Waves

We present an extension for the 2D (zonal mean) version of our Numerical Spectral Mode (NSM) that incorporates Hines' Doppler spread parameterization (DSP) for small scale gravity waves (GW). This model is applied to describe the seasonal variations and the semi-annual and quasi-biennial oscillations (SAO and QBO). Our earlier model reproduced the salient features of the mean zonal circulation in the middle atmosphere, including the QBO extension into the upper mesosphere inferred from UARS measurements. In the present model we incorporate also tropospheric heating to reproduce the upwelling at equatorial latitudes associated with the Brewer-Dobson circulation that affects significantly the dynamics of the stratosphere as Dunkerton had pointed out. Upward vertical winds increase the period of the QBO observed from the ground. To compensate for that, one needs to increase the eddy diffusivity and the GW momentum flux, bringing the latter closer to values recommended in the DSP. The QBO period in the model is 30 months (mo), which is conducive to synchronize this oscillation with the seasonal cycle of solar forcing. Multi-year interannual oscillations are generated through wave filtering by the solar driven annual oscillation in the zonal circulation. Quadratic non-linearities generate interseasonal variations to produce a complicated pattern of variability associated with the QBO. The computed temperature amplitudes for the SAO and QBO are in substantial agreement with observations at equatorial and extratropical latitudes. At high latitudes, however, the observed QBO amplitudes are significantly larger, which may be a signature of propagating planetary waves not included in the present model. The assumption of hydrostatic equilibrium not being imposed, we find that the effects from the vertical Coriolis force associated with the equatorial oscillations are large for the vertical winds and significant for the temperature variations even outside the tropics but are relatively small for the zonal winds.

Mayr, H. G.↗

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.↗

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.↗

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.↗

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.↗

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↗