Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fitting”

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,621 records · Page 90

Laboratory instrumentation and techniques for characterizing multi-junction solar cells for space applications

Multi-junction solar cells are attractive for space applications because they can be designed to convert a larger fraction of AMO into electrical power at a lower cost than single-junction cells. The performance of multi-junction cells is much more sensitive to the spectral irradiance of the illuminating source than single-junction cells. The design of high efficiency multi-junction cells for space applications requires matching the optoelectronic properties of the junctions to AMO spectral irradiance. Unlike single-junction cells, it is not possible to carry out quantum efficiency measurements using only a monochromatic probe beam and determining the cell short-circuit current assuming linearity of the quantum efficiency. Additionally, current-voltage characteristics can not be calculated from measurements under non-AMO light sources using spectral-correction methods. There are reports in the literature on characterizing the performance of multi junction cells by measuring and convoluting the quantum efficiency of each junction with the spectral irradiance; the technique is of limited value for the characterization of cell performance under AMO power-generating conditions. We report the results of research to develop instrumentation and techniques for characterizing multi junction solar cells for space . An integrated system is described which consists of a standard lamp, spectral radiometer, dual-source solar simulator, and personal computer based current-voltage and quantum efficiency equipment. The spectral radiometer is calibrated regularly using the tungsten-halogen standard lamp which has a calibration based on NIST scales. The solar simulator produces the light bias beam for current-voltage and cell quantum efficiency measurements. The calibrated spectral radiometer is used to 'fit' the spectral irradiance of the dual-source solar simulator to WRL AMO data. The quantum efficiency apparatus includes a monochromatic probe beam for measuring the absolute cell quantum efficiency at various voltage biases, including the voltage bias corresponding to the maximum-power point under AMO light bias. The details of the procedures to 'fit' the spectral irradiance to AMO will be discussed. An assessment of the role of the accuracy of the 'fit' of the spectral irradiance and probe beam intensity on measured cell characteristics will be presented. quantum efficiencies were measured with both spectral light bias and AMO light bias; the measurements show striking differences. Spectral irradiances were convoluted with cell quantum efficiencies to calculate cell currents as function of voltage. The calculated currents compare with measured currents at the 1% level. Measurements on a variety of multi-junction cells will be presented. The dependence of defects in junctions on cell quantum efficiencies measured under light and voltage bias conditions will be presented. Comments will be made on issues related to standards for calibration, and limitations of the instrumentation and techniques. Expeditious development of multi-junction solar cell technology for space presents challenges for cell characterization in the laboratory.

Woodyard, James R.↗

Solar cell angle of incidence corrections

Literature on solar array angle of incidence corrections was found to be sparse and contained no tabular data for support. This lack along with recent data on 27 GaAs/Ge 4 cm by 4 cm cells initiated the analysis presented in this paper. The literature cites seven possible contributors to angle of incidence effects: cosine, optical front surface, edge, shadowing, UV degradation, particulate soiling, and background color. Only the first three are covered in this paper due to lack of sufficient data. The cosine correction is commonly used but is not sufficient when the incident angle is large. Fresnel reflection calculations require knowledge of the index of refraction of the coverglass front surface. The absolute index of refraction for the coverglass front surface was not known nor was it measured due to lack of funds. However, a value for the index of refraction was obtained by examining how the prediction errors varied with different assumed indices and selecting the best fit to the set of measured values. Corrections using front surface Fresnel reflection along with the cosine correction give very good predictive results when compared to measured data, except there is a definite trend away from predicted values at the larger incident angles. This trend could be related to edge effects and is illustrated by a use of a box plot of the errors and by plotting the deviation of the mean against incidence angle. The trend is for larger deviations at larger incidence angles and there may be a fourth order effect involved in the trend. A chi-squared test was used to determine if the measurement errors were normally distributed. At 10 degrees the chi-squared test failed, probably due to the very small numbers involved or a bias from the measurement procedure. All other angles showed a good fit to the normal distribution with increasing goodness-of-fit as the angles increased which reinforces the very small numbers hypothesis. The contributed data only went to 65 degrees from normal which prevented any firm conclusions about extreme angle effects although a trend in the right direction was seen. Measurement errors were estimated and found to be consistent with the conclusions that were drawn. A controlled experiment using coverglasses and cells from the same lots and extending to larger incidence angles would probably lead to further insight into the subject area.

Burger, Dale R.↗

Thermospheric nitric oxide from the ATLAS 1 and Spacelab 1 missions

Spectral and spatial images obtained with the Imaging Spectrometric Observatory on the ATLAS 1 and Spacelab 1 missions are used to study the ultraviolet emissions of nitric oxide in the thermosphere. By synthetically fitting the measured NO gamma bands, intensities are derived as a function of altitude and latitude. We find that the NO concentrations inferred from the ATLAS 1 measurements are higher than predicted by our thermospheric airglow model and tend to lie to the high side of a number of earlier measurements. By comparison with synthetic spectral fits, the shape of the NO gamma bands is used to derive temperature as a function of altitude. Using the simultaneous spectral and spatial imaging capability of the instrument, we present the first simultaneously acquired altitude images of NO gamma band temperature and intensity in the thermosphere. The lower thermospheric temperature images show structure as a function of altitude. The spatial imaging technique appears to be a viable means of obtaining temperatures in the middle and lower thermosphere, provided that good information is also obtained at the higher altitudes, as the contribution of the overlying, hotter NO is nonnegligible. By fitting both self-absorbed and nonabsorbed bands of the NO gamma system, we show that the self absorption effects are observable up to 200 km, although small above 150 km. The spectral resolution of the instrument (1.6 A) allows separation of the N(+)(S-5) doublet, and we show the contribution of this feature to the combination of the NO gamma (1, 0) band and the N(+)(S-5) doublet as a function of altitude (less than 10% below 200 km). Spectral images including the NO delta bands support previous findings that the fluorescence efficiency is much higher than that determined from laboratory measurements. The Spacelab 1 data indicate the presence of a significant population of hot NO in the vehicle environment of that early shuttle mission.

Torr, Marsha R.↗

Improved Mechanical Seals For Tubes

Improved tube-to-fitting seals made by application of soft metal pieces to either ends of tubes or to interior of fittings. Metal silver, gold, platinum, tin, or other easily malleable metal selected to be inert with fluid of use. Metal plated, evaporated, or sputter-deposited. Coat of soft metal thin, so no changes required in design or tolerances of standard fittings. Technique suitable for hydraulic or other fluid systems in which leakage from mechanically joined sections of tubes occur.

Babel, Henry W.↗

A Snapshot of the Continuous Emission of the Active Galactic Nucleus in NGC 3783 from Gamma-Ray to Radio Wavelengths

To better understand the physical processes that produce the continuous emission in active galactic nuclei (AGN), a snapshot of the overall continuous energy distribution of NGC 3783, from gamma ray to radio wavelengths, has been obtained within the framework of the World Astronomy Days. The data collected in this campaign are from GRO, ROSAT, Voyager 2, IUE, HST, CTIO, SAAO, and the VLA. Great care has been taken in disentangling the genuine AGN continusous emission from other contributions; depending on the waveband, the latter might be (1) unrelated contaminating sources in cases where the instrument field of view is large (2) components within which the AGN is embedded, such as the stellar bulge population which accounts for a significant fraction of the optical continuum, and free-bound and FE2 blends wich contribute to the ultraviolet flux. After correction for these other contributins, the continuous emission of the isolated AGN appears to be rather flat (i.e., approximately equal energy per unit logarithmic frequency) from soft gamma ray to infrared wavelengths. At high energies (0.1 MeV to 0.1 keV), the AGN continuum can be fitted by a power law F nu approaches Nu(exp -a) with a spectral index of alpha approximately 1. At longer wavelengths, two excesses above this power law ('bumps') appear: in the ultraviolet, the classical big blue bump, which can be interpreted as thermal emission from the accretion disc surrounding a massive black hole, and in the infrared, a second bump which can be ascribed to thermal emission from dust in the vicinity of the AGN, heated by ultraviolet radiation from the central source. By fitting accretion-disk models to the observed AGN spectral energy distribution, we find values for the accretion disk innermost temperature, accretion rate, and black hole mass, with some differences that depend on whether or not we extrapolate the high energy power law up to infrared wavelengths. A fit to the IR bump above the extended alpha equals 1 power law suggests the presence of a dust component covering the region from a distance rho approximately equals 80 light days (hot grains at a temperature of approximately equals 1500 K) to rho approximately equals 60 light years (cool grains at T approximately equals 200 K). The total mass of dust is around 60 solar masses.

Alloin, D.↗

The 11 micron Silicon Carbide Feature in Carbon Star Shells

Silicon carbide (SiC) is known to form in circumstellar shells around carbon stars. SiC can come in two basic types - hexagonal alpha-SiC or cubic beta-SiC. Laboratory studies have shown that both types of SiC exhibit an emission feature in the 11-11.5 micron region, the size and shape of the feature varying with type, size and shape of the SiC grains. Such a feature can be seen in the spectra of carbon stars. Silicon carbide grains have also been found in meteorites. The aim of the current work is to identity the type(s) of SiC found in circumstellar shells and how they might relate to meteoritic SiC samples. We have used the CGS3 spectrometer at the 3.8 m UKIRT to obtain 7.5-13.5 micron spectra of 31 definite or proposed carbon stars. After flux-calibration, each spectrum was fitted using a chi(exp 2)-minimisation routine equipped with the published laboratory optical constants of six different samples of small SiC particles, together with the ability to fit the underlying continuum using a range of grain emissivity laws. It was found that the majority of observed SiC emission features could only be fitted by alpha-SiC grains. The lack of beta-SiC is surprising, as this is the form most commonly found in meteorites. Included in the sample were four sources, all of which have been proposed to be carbon stars, that appear to show the SiC feature in absorption.

Speck, A. K.↗

Hubble Space Telescope Images of Nearby Luminous Quasars: Results for Eight Quasars and Tests of the Detection Sensitivity - 2

Observations with the Wide-Field Camera of the Hubble Space Telescope (HST) are presented for eight intrinsically luminous quasars with redshifts between 0.16 and 0.29. These observations, when combined with a similar HST study of the quasar PKS 2349-014, show that luminous nearby quasars exist in a variety of environments. Seven companion galaxies brighter than M(V) = 16.5 (H(sub 0) = 100 km s(sup -1) Mpc(sup -1), Omega(sub 0) = 1.0) lie within a projected distance of 25 kpc of the quasars; three of the companions are located closer than 3'' (6 kpc projected distance) from the quasars, well within the volume that would be enclosed by a typical L* host galaxy. The observed association of quasars and companion galaxies is statistically significant and may he an important element in the luminous-quasar phenomenon. Apparent host galaxies are detected for three of the quasars: PG 1116+215, 3C 273, and PG 1444+407; the hosts have an average absolute magnitude of about 0.6 mag brighter than L*. The agreement between the previously published major-axis directions in ground-based images and in the present HST images of 3C 273 and PG 1444+407 constitutes important evidence supporting the reality of these candidate host galaxies. Upper limits are placed on the visual-band brightnesses of representative galactic hosts for all the quasars. These limits are established by placing galaxy images obtained with HST underneath the quasars and measuring at what faintness level the known galaxies are detected. On average, the HST spirals would have been detected if they were as faint as 1 mag below L*, and the early-type galaxies could have been detected down to a brightness level of about L*, where L* is the Schechter characteristic luminosity of field galaxies. Smooth, featureless galaxy models (exponential disks or de Vaucouleurs profiles) are fitted to the residual light after a best-fitting point source is subtracted from the quasar images. The results show that smooth spiral galaxies brighter than, on average, about L*, would have been detected. These upper limits, or possible detections, are consistent with, for example, the eight luminous quasars studied in this paper, occurring in host galaxies that have a Shechter luminosity function with a lower cutoff in the range 0.01-0.1 L*. Tests are performed to determine if our failure to detect, in some cases, luminous host galaxies could be an artifact caused by our analysis procedures. These tests include comparing the measured point-spread function (PSF) for our HST observations with the PSFs used in previous ground-based studies of host galaxies, measuring the fluctuations in the sky signals that were subtracted from the quasar images, evaluating empirically the effects of using different stellar PSFs in the analysis, carrying out the subtraction of the stellar (nuclear) source in different ways, creating and analyzing artificial active galactic nuclei (AGNs) with known surface brightnesses, and fitting the observed quasar light to an analytic model that includes a host galaxy.

Bahcall, John N.↗

The alpha Centauri Line of Sight: D/H Ratio, Physical Properties of Local Interstellar Gas, and Measurement of Heated Hydrogen (The 'Hydrogen Wall') Near the Heliopause

We analyze high-resolution spectra of the nearby (1.34 pc) stars alpha Cen A (G2 V) and alpha Cen B (K1 V), which were obtained with the Goddard High Resolution Spectrograph on the Hubble Space Telescope. The observations consist of echelle spectra of the Mg II 2800 A and Fe II 2599 A resonance lines and the Lyman-alpha lines of hydrogen and deuterium. The interstellar gas has a velocity (v = - 18.0 +/- 0.2 km/s) consistent with the local flow vector proposed for this line of sight by Lailement & Berlin (1992). The temperature and nonthermal velocity inferred from the Fe II, Mg II, and D I line profiles are T = 5400 +/- 500 K and xi = 1.20 +/- 0.25 km/s, respectively. However, single-component fits to the H I Lyman-alpha lines yield a Doppler parameter (b(sub HI) = 11.80 km/s) that implies a significantly warmer temperature of 8350 K, and the velocity of the H I absorption (v = - 15.8 +/- 0.2 km/s) is redshifted by about 2.2 km/s with respect to the Fe II, Mg II, and D I lines. The one-component model of the interstellar gas suggests natural logarithm N base HI = 18.03 +/- 0.01 and D/H = (5.7 +/- 0.2) x 10(exp -6) . These parameters lead to a good fit to the observed spectra, but this model does not explain the higher temperature and redshift of H I relative to the other interstellar lines. The most sensible way to resolve the discrepancy between H(I) and the other lines is to add a second absorption component to the H(I) lines. This second component is hotter (T approx. equals 30,000 K), is redshifted relative to the primary component by 2-4 km/s, and has a column density too low to be detected in the Fe(II), Mg(II), and D(I) lines. We propose that the gas responsible for this component is located near the heliopause, consisting of the heated H I gas from the interstellar medium that is compressed by the solar wind. This so-called 'hydrogen wall' is predicted by recent multifluid gasdynamical models of the interstellar gas and solar wind interaction. Our data provide the first measurements of the temperature and column density of H(I) in the hydrogen wall. After considering the effects that a corresponding hydrogen wall around alpha Cen would have on our analysis, our best estimates for the parameters of the solar hydrogen wall are natural log N(sup (2))(H(I)) = 14.74 +/- 0.24, b(sup (2))(H(I)) = 21.9 +/- 1.7 km/s (corresponding to T = 29,000 +/- 5000 K), and v(sup (2))(H(I)) greater than -16km/s. Unfortunately, the existence of this heated H(I) reduces our ability to compute the H(I) column density of the interstellar medium accurately because, with slight alterations to our assumed stellar Lyman-alpha profiles, we discovered that acceptable two-component fits also exist with natural log N(H(I))approx. 17.6. We, therefore, quote large error bars for the H I column density along the alpha Cen line of sight, natural log N(H(I)) = 17.80 +/- 0.30. For this range in N(H(I)), n(H(I)) = 0.15 /cu.cm (+/- a factor of 2) and D/H = (0.5-1.9) x 10(exp -5). This is the first direct measurement of the H(I) density in a local cloud and allows us to predict the distance from the Sun to the edge of the local cloud along various lines of sight. This range in D/H is consistent with the value D/H = 1.6 x 10(exp -5) previously derived for the Capella and Procyon lines of sight. We cannot tell whether D/H ratio varies or is constant in the local interstellar medium, but we do find that the D(I)/Mg(II) ratio for the alpha Cen line of sight is about 4 times smaller than for the Capella and Procyon lines of sight. Therefore, either D/H or the Mg depletion varies significantly over distance scales of only a few parsecs.

Linsky, Jeffrey L.↗

Simultaneous X-ray and Far-Ultraviolet Spectra of AGN with ASCA and HUT

We obtained ASCA spectra of the Seyfert 1 galaxy NGC 3516 in March 1995. Simultaneous far-UV observations were obtained with the Hopkins Ultraviolet Telescope on the Astro-2 shuttle mission. The ASCA spectrum shows a lightly absorbed power law of energy index 0.78. The low energy absorbing column is significantly less than previously seen. Prominent 0 VII and 0 VIII absorption edges are visible, but, consistent with the much lower total absorbing column, no Fe K absorption edge is detectable. A weak, narrow Fe K(alpha) emission line from cold material is present as well as a broad Fe K(alpha) line. These features are similar to those reported in other Seyfert 1 galaxies. A single warm absorber model provides only an imperfect description of the low energy absorption. In addition to a highly ionized absorber with ionization parameter U = 1.66 and a total column density of 1.4 x 10(exp 22)/sq cm, adding a lower ionization absorber with U = 0.32 and a total column of 6.9 x 10(exp 21)/sq cm significantly improves the fit. The contribution of resonant line scattering to our warm absorber models limits the Doppler parameter to less than 160 km/s at 90% confidence. Turbulence at the sound speed of the photoionized gas provides the best fit. None of the warm absorber models fit to the X-ray spectrum can match the observed equivalent widths of all the UV absorption lines. Accounting for the X-ray and UV absorption simultaneously requires an absorbing region with a broad range of ionization parameters and column densities.

Kriss, Gerard A.↗

Advanced Models of Accretion Disk Atmospheres and Spectra for Close Binary Stars

This work led to the development of code for fitting models to data, and to an understanding of the nature of the models which enabled a more rapid search of 'parameter space' for optimal fits to spectral data sets. The code was used to find optimal fits to IUE spectra of quiescent dwarf novae that have been reported to show evidence for the white dwarf. The models consisted of a white dwarf component and an accretion disk with boundary conditions appropriate for the choice of the white dwarf. The preliminary work has strengthened the initial impression that accretion disk spectra can mimic the appearance of white dwarf spectra in the short-wavelength ultraviolet, so that additional constraints (such as distance) are needed to distinguish to two cases.

Wade, Richard A.↗

Incorporating the APS Catalog of the POSS I and Image Archive in ADS

The primary purpose of this contract was to develop the software to both create and access an on-line database of images from digital scans of the Palomar Sky Survey. This required modifying our DBMS (called Star Base) to create an image database from the actual raw pixel data from the scans. The digitized images are processed into a set of coordinate-reference index and pixel files that are stored in run-length files, thus achieving an efficient lossless compression. For efficiency and ease of referencing, each digitized POSS I plate is then divided into 900 subplates. Our custom DBMS maps each query into the corresponding POSS plate(s) and subplate(s). All images from the appropriate subplates are retrieved from disk with byte-offsets taken from the index files. These are assembled on-the-fly into a GIF image file for browser display, and a FITS format image file for retrieval. The FITS images have a pixel size of 0.33 arcseconds. The FITS header contains astrometric and photometric information. This method keeps the disk requirements manageable while allowing for future improvements. When complete, the APS Image Database will contain over 130 Gb of data. A set of web pages query forms are available on-line, as well as an on-line tutorial and documentation. The database is distributed to the Internet by a high-speed SGI server and a high-bandwidth disk system. URL is http://aps.umn.edu/IDB/. The image database software is written in perl and C and has been compiled on SGI computers with MIX5.3. A copy of the written documentation is included and the software is on the accompanying exabyte tape.

Humphreys, Roberta M.↗

X-Ray Spectral Variability in NGC 7469

We present analyses of two Ginga observations and two observations from the ROSAT database of NGC 7469, focusing on the spectral variability observed on timescales of days and longer. During the 1988 Ginga observation, the hardness ratio (8-21 keV/3.4-5.7 keV) increased significantly as the total flux decreased by 30%. As the spectrum is well fit by the reflection model and since the spectra variability dominates the higher energy band, this could be explained by either a variation in the power law index or in the effective covering fraction of the reflecting material. This ambiguity is inherent in reflection modeling of Ginga spectra from moderate flux Seyfert 1 galaxies. Assuming that the power law index did not change, we find that the reflected flux is consistent with being constant, suggesting that much of the reflecting material may be located more than 3 light-days from the continuum source with the molecular torus being a plausible site. This scenario is also supported by the report of a narrow rather than broad iron K-alpha line in the ASCA data by Guainazzi et al. NGC 7469 was faint during the 1989 Ginga observation, but variability was observed with doubling timescale of 5 hr, and the spectrum was harder. A reflection component could not be constrained, and the change in the spectrum could be explained by an increase in neutral absorption. The brighter of two ROSAT spectra was significantly softer, and in both spectra there was evidence of spectral complexity, as has been previously reported by Turner, George, & Mushotzky and Brandt et al. The spectrum could be fit by a variety of two-component models, including a warm absorber model, an ionized disk model, and a thermal model with single-component blackbody spectrum, but joint fitting of the 1988 average Ginga spectrum and the nonsimultaneous ROSAT spectra favored thermal models, and other models required an anomalously high reflection ratio. This model is supported by the observation of a soft excess component and the lack of ionized absorption edges in the ASCA spectrum by Guainazzi et al. The long-term spectral variability could be explained by relative variability between the power-law and soft excess component normalizations, perhaps implying that hard X-ray reprocessing in thermal material does not dominate on long timescales.

Leighly, Karen↗

The Local ISM and its Interaction with the Winds of Nearby Late-type Stars

We present new Goddard High-Resolution Spectrograph (GHRS) observations of the Ly-alpha and Mg II absorption lines seen toward the nearby stars 61 Cyg A and 40 Eri A. We use these data to measure interstellar properties along these lines of sight and to search for evidence of circumstellar hydrogen walls, which are produced by collisions between the stellar winds and the Local InterStellar Medium (LISM). We were able to model the Ly-alpha lines of both stars without hydrogen-wall absorption components, but for 61 Cyg A the fit required a stellar Ly-alpha, line profile with an improbably deep self-reversal, and for 40 Eri A the fit required a very low deuterium-to-hydrogen ratio that is inconsistent with previous GHRS measurements. Since these problems could be rectified simply by including stellar hydrogen-wall components with reasonable attributes, our preferred fits to the data include these components. We have explored several ways in which the hydrogen-wall properties measured here and in previous work can be used to study stellar winds and the LISM. We argue that the existence of a hydrogen wall around 40 Eri A and a low H I column density along that line of sight imply that either the interstellar density must decrease toward 40 Eri A or the hydrogen ionization fraction (chi) must increase. We find that hydrogen-wall temperatures are larger for stars with faster velocities through the LISM. The observed temperature-velocity relation is consistent with the predictions of hydromagnetic shock jump conditions. More precise comparison of the data and the jump conditions suggests crude upper limits for both chi and the ratio of magnetic to thermal pressure in the LISM (alpha): chi less than 0.6 and alpha less than 2. The latter upper limit corresponds to a limit on the LISM magnetic field of B less than 5 micro G. These results imply that the plasma Mach number of the interstellar wind flowing into the heliosphere is M(sub A) greater than 1.3, which indicates that the collision is supersonic and that there should therefore be a bow shock outside the heliopause in the upwind direction. Finally, we estimate stellar wind pressures (P sub wind) from the measured hydrogen-wall column densities. These estimates represent the first empirical measurements of wind properties for late-type main-sequence stars. The wind pressures appear to be correlated with stellar X-ray surface fluxes, F(x), in a manner consistent with the relation P(wind) varies as F(x)(exp -1/2), a relation that is also consistent with the variations of P(sub wind) and F(sub x) observed during the solar activity cycle. If this relation can in fact be generalized to solar-like stars, as is suggested by our data, then it is possible to estimate stellar wind properties simply by measuring stellar X-rays. One implication of this is that stellar wind pressures and mass-loss rates are then predicted to increase with time, since F(sub x) is known to decrease with stellar age.

Wood, Brian E.↗

Development of a Computer Architecture to Support the Optical Plume Anomaly Detection (OPAD) System

The NASA OPAD spectrometer system relies heavily on extensive software which repetitively extracts spectral information from the engine plume and reports the amounts of metals which are present in the plume. The development of this software is at a sufficiently advanced stage where it can be used in actual engine tests to provide valuable data on engine operation and health. This activity will continue and, in addition, the OPAD system is planned to be used in flight aboard space vehicles. The two implementations, test-stand and in-flight, may have some differing requirements. For example, the data stored during a test-stand experiment are much more extensive than in the in-flight case. In both cases though, the majority of the requirements are similar. New data from the spectrograph is generated at a rate of once every 0.5 sec or faster. All processing must be completed within this period of time to maintain real-time performance. Every 0.5 sec, the OPAD system must report the amounts of specific metals within the engine plume, given the spectral data. At present, the software in the OPAD system performs this function by solving the inverse problem. It uses powerful physics-based computational models (the SPECTRA code), which receive amounts of metals as inputs to produce the spectral data that would have been observed, had the same metal amounts been present in the engine plume. During the experiment, for every spectrum that is observed, an initial approximation is performed using neural networks to establish an initial metal composition which approximates as accurately as possible the real one. Then, using optimization techniques, the SPECTRA code is repetitively used to produce a fit to the data, by adjusting the metal input amounts until the produced spectrum matches the observed one to within a given level of tolerance. This iterative solution to the original problem of determining the metal composition in the plume requires a relatively long period of time to execute the software in a modern single-processor workstation, and therefore real-time operation is currently not possible. A different number of iterations may be required to perform spectral data fitting per spectral sample. Yet, the OPAD system must be designed to maintain real-time performance in all cases. Although faster single-processor workstations are available for execution of the fitting and SPECTRA software, this option is unattractive due to the excessive cost associated with very fast workstations and also due to the fact that such hardware is not easily expandable to accommodate future versions of the software which may require more processing power. Initial research has already demonstrated that the OPAD software can take advantage of a parallel computer architecture to achieve the necessary speedup. Current work has improved the software by converting it into a form which is easily parallelizable. Timing experiments have been performed to establish the computational complexity and execution speed of major components of the software. This work provides the foundation of future work which will create a fully parallel version of the software executing in a shared-memory multiprocessor system.

Katsinis, Constantine↗

A Numerical Method for Calculating Stellar Occultation Light Curves from an Arbitrary Atmospheric Model

We present a method for speeding up numerical calculations of a light curve for a stellar occultation by a planetary atmosphere with an arbitrary atmospheric model that has spherical symmetry. This improved speed makes least-squares fitting for model parameters practical. Our method takes as input several sets of values for the first two radial derivatives of the refractivity at different values of model parameters, and interpolates to obtain the light curve at intermediate values of one or more model parameters. It was developed for small occulting bodies such as Pluto and Triton, but is applicable to planets of all sizes. We also present the results of a series of tests showing that our method calculates light curves that are correct to an accuracy of 10(exp -4) of the unocculted stellar flux. The test benchmarks are (i) an atmosphere with a l/r dependence of temperature, which yields an analytic solution for the light curve, (ii) an atmosphere that produces an exponential refraction angle, and (iii) a small-planet isothermal model. With our method, least-squares fits to noiseless data also converge to values of parameters with fractional errors of no more than 10(exp -4), with the largest errors occurring in small planets. These errors are well below the precision of the best stellar occultation data available. Fits to noisy data had formal errors consistent with the level of synthetic noise added to the light curve. We conclude: (i) one should interpolate refractivity derivatives and then form light curves from the interpolated values, rather than interpolating the light curves themselves; (ii) for the most accuracy, one must specify the atmospheric model for radii many scale heights above half light; and (iii) for atmospheres with smoothly varying refractivity with altitude, light curves can be sampled as coarsely as two points per scale height.

Chamberlain, D. M.↗

Multitaper Spectral Analysis and Wavelet Denoising Applied to Helioseismic Data

Estimates of solar normal mode frequencies from helioseismic observations can be improved by using Multitaper Spectral Analysis (MTSA) to estimate spectra from the time series, then using wavelet denoising of the log spectra. MTSA leads to a power spectrum estimate with reduced variance and better leakage properties than the conventional periodogram. Under the assumption of stationarity and mild regularity conditions, the log multitaper spectrum has a statistical distribution that is approximately Gaussian, so wavelet denoising is asymptotically an optimal method to reduce the noise in the estimated spectra. We find that a single m-upsilon spectrum benefits greatly from MTSA followed by wavelet denoising, and that wavelet denoising by itself can be used to improve m-averaged spectra. We compare estimates using two different 5-taper estimates (Stepian and sine tapers) and the periodogram estimate, for GONG time series at selected angular degrees l. We compare those three spectra with and without wavelet-denoising, both visually, and in terms of the mode parameters estimated from the pre-processed spectra using the GONG peak-fitting algorithm. The two multitaper estimates give equivalent results. The number of modes fitted well by the GONG algorithm is 20% to 60% larger (depending on l and the temporal frequency) when applied to the multitaper estimates than when applied to the periodogram. The estimated mode parameters (frequency, amplitude and width) are comparable for the three power spectrum estimates, except for modes with very small mode widths (a few frequency bins), where the multitaper spectra broadened the modest compared with the periodogram. We tested the influence of the number of tapers used and found that narrow modes at low n values are broadened to the extent that they can no longer be fit if the number of tapers is too large. For helioseismic time series of this length and temporal resolution, the optimal number of tapers is less than 10.

Komm, R. W.↗

Interactions in Undersaturated and Supersaturated Lysozyme Solutions: Static and Dynamic Light Scattering Results

We have performed multiangle static and dynamic light scattering studies of lysozyme solutions at pH=4.7. The Rayleigh ratio R(sub g) and the collective diffusion coefficient D(sub c) were determined as function of both protein concentration c(sub p) and salt concentration c(sub s) with two different salts. At low salt concentrations, the scattering ratio K(sub c)(sub p)/R(sub theta) and diffusivity increased with protein concentration above the values for a monomeric, ideal solution. With increasing salt concentration this trend was eventually reversed. The hydrodynamic interactions of lysozyme in solution, extracted from the combination of static and dynamic scattering data, decreased significantly with increasing salt concentration. These observations reflect changes in protein interactions, in response to increased salt screening, from net repulsion to net attraction. Both salts had the same qualitative effect, but the quantitative behavior did not scale with the ionic strength of the solution. This indicates the presence of salt specific effects. At low protein concentrations, the slopes of K(sub c)(sub p)/R(sub theta) and D(sub c) vs c(sub p) were obtained. The dependence of the slopes on ionic strength was modeled using a DLVO potential for colloidal interactions of two spheres, with the net protein charge Z(sub e) and Hamaker constant A(sub H) as fitting parameters. The model reproduces the observed variations with ionic strength quite well. Independent fits to the static and dynamic data, however, led to different values of the fitting parameters. These and other shortcomings suggest that colloidal interaction models alone are insufficient to explain protein interactions in solutions.

Muschol, Martin↗

Modelling of Collision Induced Absorption Spectra Of H2-H2 Pairs for the Planetary Atmospheres Structure: The Second Overtone Band

The main objective of the proposal was to model the collision induced, second overtone band of gaseous hydrogen at low temperatures. The aim of this work is to assist planetary scientists in their investigation of planetary atmospheres, mainly those of Uranus and Neptune. The recently completed extended database of collision induced dipole moments of hydrogen pairs allowed us, for the first time, to obtain dipole moment matrix elements responsible for the roto-vibrational collision induced absorption spectra of H2-H2 in the second overtone band. Despite our numerous attempts to publish those data, the enormous volume of the database did not allow us to do this. Instead, we deposited the data on a www site. The final part of this work has been partially supported by NASA, Division for Planetary Atmospheres. In order to use our new data for modelling purpose, we first needed to test how well we can reproduce the existing experimental data from theory, when using our new input data. Two papers resulted from this work. The obtained agreement between theoretical results and the measurements appeared to be within 10-30%. The obviously poorer agreement than observed for the first H2 overtone, the fundamental, and the rototranslational bands can be attributed to the fact that dipole moments responsible for the second overtone are much weaker, therefore susceptible to larger numerical uncertainties. At the same time, the intensity of the second overtone band is much weaker and therefore it is much harder to be measured accurately in the laboratory. We need to point out that until now, no dependable model of the 2nd overtone band was available for modelling of the planetary atmospheres. The only one, often referred to in previous works on Uranian and Neptune's atmospheres, uses only one lineshape, with one (or two) parameter(s) deduced at the effective temperature of Uranus (by fitting the planetary observation). After that, the parameter(s) was(were) made temperature dependent according to some very simple relation. Summarizing, no reliable temperature-dependent model has been available yet. Our approach was a bit different from similar attempts done earlier, on account of the poorer agreement of theory with experiment. We needed to resort to some semi-empirical procedure. While we were in a favourable position to be able to rely on the physical input data, these, apparently, did not supply the most dependable predictions (simply because the results did not agree well enough with experimental data). On the other hand, the relative deviations between the theory and experiment were comparable at 77 and at 298 K. That fact indicated that theory is capable of predicting the temperature dependence of the absorption spectra well. We have thus chosen the "middle way". We have fitted the existing measurements with many 3- parameter lineshapes, in order to achieve the closest fit.

Borysow, Aleksandra↗