Search NASA⌕ Search

SEARCH · Search NASA

Results for “Sparse approximation”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

TOI-908: a planet at the edge of the Neptune desert transiting a G-type star

ABSTRACT We present the discovery of an exoplanet transiting TOI-908 (TIC-350153977) using data from TESS sectors 1, 12, 13, 27, 28, and 39. TOI-908 is a T = 10.7 mag G-dwarf (Teff = 5626 ± 61 K) solar-like star with a mass of 0.950 ± 0.010 M⊙ and a radius of 1.028 ± 0.030 R⊙. The planet, TOI-908 b, is a 3.18 ± 0.16 R⊕ planet in a 3.18 d orbit. Radial velocity measurements from HARPS reveal TOI-908 b has a mass of approximately 16.1 ± 4.1 M⊕, resulting in a bulk planetary density of $2.7^{+0.2}_{-0.4}$ g cm−3. TOI-908 b lies in a sparsely populated region of parameter space known as the Neptune desert. The planet likely began its life as a sub-Saturn planet before it experienced significant photoevaporation due to X-rays and extreme ultraviolet radiation from its host star, and is likely to continue evaporating, losing a significant fraction of its residual envelope mass.

Hawthorn, Faith (ORCID:000000028675182X)↗
NASA NTRS1991

Solution of three-dimensional afterbody flow using reduced Navier-Stokes equations

The flow over afterbody geometries was investigated using the reduced Navier-Stokes (RNS) approximation. Both pressure velocity flux-split and composites velocity primitive variable formulations were considered. Pressure or pseudopotential relaxation procedures are combined with sparse matrix or coupled strongly implicit algorithms to form a three-dimensional solver for general non-orthogonal coordinates. Three-dimensional subsonic and transonic viscous/inviscid interacting flows were evaluated. Solutions with and without regions of recirculation were obtained.

Almahroos, H. M. H.↗
NASA NTRS2013

Sweeping Jet Actuators - A New Design Tool for High Lift Generation

Active Flow Control (AFC) experiments performed at the Caltech Lucas Wind Tunnel on a generic airplane vertical tail model proved the effectiveness of sweeping jets in improving the control authority of a rudder. The results indicated that a momentum coefficient (C(sub u)) of approximately 2% increased the side force in excess of 50% at the maximum conventional rudder deflection angle in the absence of yaw. However, sparsely distributed actuators providing a collective C(sub u) approx. = 0.1% were able to increase the side force in excess of 20%. This result is achieved by reducing the spanwise flow along the swept back rudder and its success is attributed to the large sweep back angle of the vertical tail. This current effort was sponsored by the NASA Environmentally Responsible Aviation (ERA) project.

Graff, Emilio↗
DOE OSTI2019

Sparse Bayesian mass mapping with uncertainties: local credible intervals

ABSTRACT Until recently, mass-mapping techniques for weak gravitational lensing convergence reconstruction have lacked a principled statistical framework upon which to quantify reconstruction uncertainties, without making strong assumptions of Gaussianity. In previous work, we presented a sparse hierarchical Bayesian formalism for convergence reconstruction that addresses this shortcoming. Here, we draw on the concept of local credible intervals (cf. Bayesian error bars) as an extension of the uncertainty quantification techniques previously detailed. These uncertainty quantification techniques are benchmarked against those recovered via Px-MALA – a state-of-the-art proximal Markov chain Monte Carlo (MCMC) algorithm. We find that, typically, our recovered uncertainties are everywhere conservative (never underestimate the uncertainty, yet the approximation error is bounded above), of similar magnitude and highly correlated with those recovered via Px-MALA. Moreover, we demonstrate an increase in computational efficiency of $\mathcal {O}(10^6)$ when using our sparse Bayesian approach over MCMC techniques. This computational saving is critical for the application of Bayesian uncertainty quantification to large-scale stage IV surveys such as LSST and Euclid.

Price, M. A.↗
NASA NTRS2004

Precision X-ray Timing of RX J0806.3+1527 with CHANDRA: Evidence for Gravitational Radiation from an Ultracompact Binary

RX J0806.3+1527 is a candidate double degenerate binary with possibly the shortest known orbital period. The source shows an approximately equal to 100% X-ray intensity modulation at the putative orbital frequency of 3.11 mHz (321.5 s). If the system is a detached, ultracompact binary gravitational radiation should drive spin-up with a magnitude of nu(sup dot) approximately 10(exp -16) Hz per second. Efforts to constrain the X-ray frequency evolution to date have met with mixed success, principally due to the sparseness of earlier observations. Here we describe the results of the first phase coherent X-ray monitoring campaign on RX J0806.3+1527 with Chandra. We obtained a total of 70 ksec of exposure in 6 epochs logarithmically spaced over 320 days. With these data we conclusively show that the X-ray frequency is increasing at a rate of 3.77 plus or minus 0.8 x 10(exp -16) Hz per second. Using the ephemeris derived from the new data we are able to phase up all the earlier Chandra and ROSAT data and show they are consistent with a constant nu(sup dot) = 3.63 plus or minus 0.06 x 10(exp -16) Hz per second over the past decade. This value appears consistent with that recently derived by Israel et al. largely from monitoring of the optical modulation, and is in rough agreement with the solutions reported initially by Hakala et al., based on ground-based optical observations. The large and stable nu(sup dot) over a decade is consistent with gravitational radiation losses driving the evolution. An intermediate polar (IP) scenario where the observed X-ray period is the spin period of an accreting white dwarf appears less tenable because the observed nu(sup dot) requires an m(sup dot) approximately equal to 4 x 10 (exp -8) solar mass yr(sup -l), that is much larger than that inferred from the observed X-ray luminosity (although this depends on the uncertain distance and bolometric corrections), and it is difficult to drive such a high m(sup dot) in a binary system with parameters consistent with all the multiwavelength data. If the ultracompact scenario is correct, then the X-ray flux cannot be powered by stable accretion which would drive the components apart, suggesting a new type of energy source (perhaps electromagnetic) may power the X-ray flux.

Strohymayer, Tod E.↗
DOE OSTI2025

Continuous surface-to-distributed acoustic sensor snapshots explain reactivation of individual natural fractures during an unconventional reservoir stimulation

ABSTRACT Fiber-optic sensing technologies allow petroleum engineering teams to detect hydraulic fracture interaction with boreholes during unconventional reservoir stimulation. In combination with high-repeatability seismic sources, the same distributed acoustic sensors (DASs) enable vertical seismic profiling (VSP) of the fracture evolution away from the boreholes. We discovered clear signatures of seismic scattering on activated fractures during nine days of continuous seismic monitoring of the fracturing stages at the Austin Chalk/Eagle Ford Field Laboratory. The present study applies a novel approach for quantitative analysis of the scattering events in terms of the evolution of the geometry and elastic stiffness of individual fractures. Our characterization strategy sequentially refines the fracture models: from a stack of 1D soft layers to 3D rectangular inclusions. First, we estimate the number of fracture locations and reflectivity using a modified sparse-spike deconvolution of the stacked VSP traces. The fracture set consists of five fractures spaced by 15–30 m with a reflectivity of approximately 1%. Then, we develop a scattering integral method to refine these estimates along with an inversion of the fracture top and bottom for each monitoring vintage. We find that, initially, some of the fractures are located above the monitoring fiber with the height of approximately 100 m. Then we integrate the seismic interpretation with the low-frequency DAS and pressure and microseismic monitoring to reconstruct the activation process of the fractures. Most likely, some of the natural fractures slowly grew downward to the monitoring fiber as a result of fluid injections in the stimulated well. This led to bright strain anomalies but did not trigger seismicity. The top of the fractures remained almost constant and were limited by a lithologic boundary/stress barrier. To our knowledge, this is the first time VSP data enabled tracking of the fracture evolution with such high spatial and temporal resolution, which was previously only available for crosswell surveys and at a much smaller scale.

Glubokovskikh, Stanislav↗
NASA NTRS2000

Chromosomal instability induced by heavy ion irradiation

PURPOSE: To establish the dose-response relationship for the induction of chromosomal instability in GM10115 cells exposed to high-energy iron ions (1 GeV/nucleon, mean LET 146 keV/microm) and gold ions (11 GeV/nucleon, mean LET 1450 keV/microm). Past work has established that sparsely ionizing X-rays can induce a long-lived destabilization of chromosomes in a dose-dependent manner at an incidence of approximately 3% per gray. The present investigation assesses the capacity of High-Z and High-energy (HZE) particles to elicit this same endpoint. MATERIALS AND METHODS: Clonal populations derived from single progenitor cells surviving heavy-ion irradiation were analyzed cytogenetically to identify those clones showing a persistent destablization of chromosomes. RESULTS: Dose-response data, with a particular emphasis at low dose (< 1.0 Gy), indicate a frequency of approximately 4% per gray for the induction of chromosomal instability in clones derived from single progenitor cells surviving exposure to iron ions. The induction of chromosomal instability by gold ions was, however, less responsive to applied dose, as the observed incidence of this phenotype varied from 0 to 10% over 1-8 Gy. Both iron and gold ions gave dose-dependent increases in the yield of chromosomal aberrations (both chromosome- and chromatid-type) measured at the first mitosis following irradiation, as well as shoulderless survival curves having D0=0.87 and 1.1 Gy respectively. CONCLUSIONS: Based on the present dose-response data, the relative biological effectiveness of iron ions is 1.3 for the induction of chromosomal instability, and this indicates that heavy ions are only slightly more efficient than X-rays at eliciting this delayed phenotype.

Non-NASA Center↗
NASA NTRS2017

Crater Morphometry and Crater Degradation on Mercury: Mercury Laser Altimeter (MLA) Measurements and Comparison to Stereo-DTM Derived Results

Two types of measurements of Mercury's surface topography were obtained by the MESSENGER (MErcury Surface Space ENvironment, GEochemisty and Ranging) spacecraft: laser ranging data from Mercury Laser Altimeter (MLA) [1], and stereo imagery from the Mercury Dual Imaging System (MDIS) camera [e.g., 2, 3]. MLA data provide precise and accurate elevation meaurements, but with sparse spatial sampling except at the highest northern latitudes. Digital terrain models (DTMs) from MDIS have superior resolution but with less vertical accuracy, limited approximately to the pixel resolution of the original images (in the case of [3], 15-75 m). Last year [4], we reported topographic measurements of craters in the D=2.5 to 5 km diameter range from stereo images and suggested that craters on Mercury degrade more quickly than on the Moon (by a factor of up to approximately 10×). However, we listed several alternative explanations for this finding, including the hypothesis that the lower depth/diameter ratios we observe might be a result of the resolution and accuracy of the stereo DTMs. Thus, additional measurements were undertaken using MLA data to examine the morphometry of craters in this diameter range and assess whether the faster crater degradation rates proposed to occur on Mercury is robust.

morphometry↗
DOE OSTI2022

Assessing entropy for catalytic processes at complex reactive interfaces

When chemical reactions are accelerated by a catalyst, entropy differences between reactants and their transient intermediates can be the driving force behind the promotion or inhibition of desired and parasitic chemical pathways. Understanding and controlling catalytic processes therefore requires both a fundamental and practicable understanding of entropy in addition to enthalpy. In unstructured media such as the vapor phase equilibrated with sparsely covered surfaces, entropy can be adequately accounted for by well-established approaches based on translational, rotational, and harmonic vibrational partition functions. However, these approximations become inadequate in more complex condensed phase environments, e.g., solid liquid interfaces of confined reaction spaces. In this chapter, we provide an overview of the state-of-art in the computational quantification of entropy and its known ramifications on catalysis. The fundamental roles of thermodynamics and kinetics in catalysis are covered in enough detail to appreciate and contextualize the computational methods employed to compute chemically accurate estimates of entropy. These methods are discussed in appropriate detail and range from the ubiquitous harmonic oscillator approximation where entropy unrelated to high frequency oscillations is typically underestimated, to enhanced free energy sampling with molecular dynamics where the desired accuracy must be weighed against the associated computational cost of obtaining it. The rising importance of machine learning and artificial intelligence in accelerating methodological progress in this field is touched upon, as well. Finally, applications, successes, and pitfalls of using these methods are provided to showcase past and present accomplishments while clarifying where improvements in both understanding and methodology are still needed.

Kollias, Loukas↗
DOE OSTI2021

Identifying Entangled Physics Relationships Through Sparse Matrix Decomposition to Inform Plasma Fusion Design

We report a sustainable burn platform through inertial confinement fusion (ICF) has been an ongoing challenge for over 50 years. Mitigating engineering limitations and improving the current design involves an understanding of the complex coupling of physical processes. While sophisticated simulation codes are used to model ICF implosions, these tools contain necessary numerical approximation but miss physical processes that limit predictive capability. Identification of relationships between controllable design inputs to ICF experiments and measurable outcomes (e.g., neutron yield, neutron velocity, areal density) from performed experiments can help guide the future design of experiments and development of simulation codes, to potentially improve the accuracy of the computational models used to simulate ICF experiments. We use sparse matrix decomposition methods to identify clusters of a few related design variables. Sparse principal component analysis (SPCA) identifies groupings that are related to the physical origin of the variables (laser, hohlraum, and capsule). A variable importance analysis finds that in addition to variables highly correlated with neutron yield, such as picket power and laser energy, variables that represent a dramatic change of the ICF design, such as number of pulse steps, are also very important. The obtained sparse components are then used to train a random forest (RF) regression surrogate for predicting total yield. The RF performance on the training and testing data compares with the performance of the RF trained using all the design variables considered. This work is intended to inform design changes in future ICF experiments by augmenting the expert intuition and simulation results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2016

TEM Analyses of Itokawa Regolith Grains and Lunar Soil Grains to Directly Determine Space Weathering Rates on Airless Bodies

Samples returned from the moon and Asteroid Itokawa by NASA's Apollo Missions and JAXA's Hayabusa Mission, respectively, provide a unique record of their interaction with the space environment. Space weathering effects result from micrometeorite impact activity and interactions with the solar wind. While the effects of solar wind interactions, ion implantation and solar flare particle track accumulation, have been studied extensively, the rate at which these effects accumulate in samples on airless bodies has not been conclusively determined. Results of numerical modeling and experimental simulations do not converge with observations from natural samples. We measured track densities and rim thicknesses of three olivine grains from Itokawa and multiple olivine and anorthite grains from lunar soils of varying exposure ages. Samples were prepared for analysis using a Leica EM UC6 ultramicrotome and an FEI Quanta 3D dual beam focused ion beam scanning electron microscope (FIB-SEM). Transmission electron microscope (TEM) analyses were performed on the JEOL 2500SE 200kV field emission STEM. The solar wind damaged rims on lunar anorthite grains are amorphous, lack inclusions, and are compositionally similar to the host grain. The rim width increases as a smooth function of exposure age until it levels off at approximately 180 nm after approximately 20 My (Fig. 1). While solar wind ion damage can only accumulate while the grain is in a direct line of sight to the Sun, solar flare particles can penetrate to mm-depths. To assess whether the track density accurately predicts surface exposure, we measured the rim width and track density in olivine and anorthite from the surface of rock 64455, which was never buried and has a surface exposure age of 2 My based on isotopic measurements. The rim width from 64455 (60-70nm) plots within error of the well-defined trend for solar wind amorphized rims in Fig. 1. Measured solar flare track densities are accurately reflecting the surface exposure of the grains. Track densities correlate with the amorphous rim thicknesses. While the space-weathered rims of anorthite grains are amorphous, the space-weathered rims on both Itokawa and lunar olivine grains show solar wind damaged rims that are not amorphous. Instead, the rims are nanocrystalline with high dislocation densities and sparse inclusions of nanophase Fe metal. The rim thicknesses on the olivine grains also correlate with track density. The Itokawa olivine grains have track densities that indicate surface exposures of approximately 10(exp 5) years. Longer exposures (up to approximately 10(exp 7) years) do not amorphize the rims, as evidenced by lunar soil olivines with high track densities (approximately 10(exp 11) cm(exp -2)). From the combined data, shown in Fig. 1, it is clear that olivine is damaged (but not amorphized) more rapidly by the solar wind compared to anorthite. The olivine damaged rim forms quickly (in approximately 10(exp 6) y) and saturates at approximately 120nm with longer exposure time. The anorthite damaged rims form more slowly, amorphize, and grow thicker than the olivine rims. This is in agreement with numerical modeling data which predicts that solar wind damaged rims on anorthite will be thicker than olivine. However, the models predict that both olivine and anorthite rims will amorphize and reach equilibrium widths in less than 10(exp 3) y, in contrast to what is observed for natural samples. Laboratory irradiation experiments, which show rapid formation of fully amorphous and blistered surfaces from simulated solar wind exposures are also in contrast to observations of natural samples. These results suggest that there is a flux dependence on the type and extent of irradiation damage that develops in olivine. This flux dependence suggests that great caution be used in extrapolating between high-flux laboratory experiments and the natural case, as demonstrated by. We constrain the space weathering rate through analysis of returned samples. Provided that the track densities and the solar wind damaged rim widths exhibited by the Itokawa grains are typical of the fine-grained regions of Itokawa, then the space weathering rate is on the order of 10(exp 5) y. Space weathering effects in lunar soils saturate within a few My of exposure while those in Itokawa regolith grains formed in approximately 10(exp 5) y. Olivine and anorthite respond differently to solar wind irradiation. The space weathering effects in olivine are particularly difficult to reconcile with laboratory irradiation studies and numerical models. Additional measurements, experiments, and modeling are required to resolve the discrepancies among the observations and calculations involving solar wind amorphization of different minerals on airless bodies.

Berger, Eve L.↗
DOE OSTI2026

A score-based diffusion model approach for adaptive learning of stochastic partial differential equation solutions

In this paper, we propose a novel framework for adaptively learning the time-evolving solutions of stochastic partial differential equations (SPDEs) using score-based diffusion models within a recursive Bayesian inference setting. SPDEs play a central role in modeling complex physical systems under uncertainty, but their numerical solutions often suffer from model errors and reduced accuracy due to incomplete physical knowledge and environmental variability. To address these challenges, we encode the governing physics into the score function of a diffusion model using simulation data and incorporate observational information via a likelihood-based correction in a reverse-time stochastic differential equation. This enables adaptive learning through iterative refinement of the solution as new data becomes available. To improve computational efficiency in high-dimensional settings, we introduce the ensemble score filter, a training-free approximation of the score function designed for real-time inference. Numerical experiments on benchmark SPDEs demonstrate the accuracy and robustness of the proposed method under sparse and noisy observations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Forward sensitivity approach for estimating eddy viscosity closures in nonlinear model reduction

In this paper, we propose a variational approach to estimate eddy viscosity using forward sensitivity method (FSM) for closure modeling in nonlinear reduced order models. FSM is a data assimilation technique that blends model's predictions with noisy observations to correct initial state and/or model parameters. We apply this approach on a projection based reduced order model (ROM) of the one-dimensional viscous Burgers equation with a square wave defining a moving shock, and the two-dimensional vorticity transport equation formulating a decay of Kraichnan turbulence. We investigate the capability of the approach to approximate an optimal value for eddy viscosity with different measurement configurations. Specifically, we show that our approach can sufficiently assimilate information either through full field or sparse noisy measurements to estimate eddy viscosity closure to cure standard Galerkin reduced order model (GROM) predictions. Therefore, our approach provides a modular framework to correct forecasting error from a sparse observational network on a latent space. We highlight that the proposed GROM-FSM framework is promising for emerging digital twin applications, where real-time sensor measurements can be used to update and optimize surrogate model's parameters.

42 ENGINEERING↗
NASA NTRS1986

Canopy reflectance modeling in a tropical wooded grassland

Geometric/optical canopy reflectance modeling and spatial/spectral pattern recognition are used to study the form and structure of savanna in West Africa. An invertible plant canopy reflectance model is tested for its ability to estimate the amount of woody vegetation cover in areas of sparsely wooded grassland from remotely sensed data. Dry woodlands and wooded grasslands, commonly referred to as savannas, are important ecologically and economically in Africa, and cover approximately forty percent of the continent by some estimates. The Sahelian and Sudanian savanna make up the important and sensitive transition zone between the tropical forests and the arid Saharan region. The depletion of woody cover, used for fodder and fuel in these regions, has become a very severe problem for the people living there. LANDSAT Thematic Mapper (TM) data is used to stratify woodland and wooded grassland into areas of relatively homogeneous canopy cover, and then by applying an invertible forest canopy reflectance model to estimate directly the height and spacing of the trees in the stands. Since height and spacing are proportional to biomass in some cases, a successful application of the segmentation/modeling techniques will allow direct estimation of woody biomass, as well as cover density, over significant areas of these valuable and sensitive ecosystems. Sahelian savanna sites in the Gourma area of Mali being used by the NASA/GIMMS project (Global Inventory Modeling and Monitoring System, at Goddard Space Flight Center), in conjunction with CIPEA/Mali (Centre International pour l'Elevage en Afrique) will be used for testing the canopy model. The model will also be tested in a Sudanian zone crop/woodland area in the Region of Segou, Mali.

Simonett, D.↗
NASA NTRS1994

Superhumps in cataclysmic binaries. 3: V795 Herculis

We report photometry of the cataclysmic variable V795 Herculis during 1990-1994. Comparison of recent with previously published data shows that the 2.78-hr periodicity has essentially disappeared; the amplitude has diminished by a factor of at least 10. In its place, there is a low-amplitude signal near the radial-velocity (orbital?) period of 2.60 hr. During 1983-1989, the star did show an obvious 2.78-hr signal. We study the published data on this signal, together with our own data under the (questionable) assumption that the signal still exists with an amplitude too small to survive the rigors of unbiased period detection. The data are really too sparse and too poorly distributed to decide the issue of phase stability. But since the photometric signal is transient and exceeds the likely orbital period by approximately 8%, the photometric waves are likely to be yet another example of 'superhumps', the phenomenon made famous by dwarf novae in superoutburst. At high frequencies, the star shows a quasi-periodic signal with P approximately 1160 s, and possibly also a stable signal with P = 1310.2 s (or 1330.4 s, the one-day alias). The beat period between these signals could be the orbital period, or that of the vanished superhump.

Patterson, Joseph↗
NASA NTRS2008

First Science with the Keck Interferometer Nuller High Spatial Resolution N-Band Observations of the Recurrent Nova RS Ophiuchi

We report observations of the nova RS 0phiuchi using the Keck Interferometer Nuller (KIN) taken approximately 3.8 days following the most recent outburst that occurred on 2006 February 12. The KIN operates in N-band from approx. 8 to 12.5 microns in a nulling mode - the sparse aperture equivalent of the conventional coronagraphic technique used in filled aperture te1escopes. In this mode the stellar light is suppressed by a destructive fringe, effectively enhancing the contrast of the circumstellar material located near the star. In a second, constructive-fringe mode, the instrument detects primarily the light from the central, bright source. These are the outer and inner spatial regimes, resprectively. We will describe the capabilities of the KIN, including these unique modes, and outline how they were key in our discovery that dust does not appear to be created in the outburst as in previous models, but instead was created between nova events. We also show how these first results from the KIN are consistent with Spitzer data. The KIN data show evidence of enhanced neutral atomic hydrogen emission and atomic metals including silicon located in the inner spatial regime near the white dwarf (WD) relative to the outer regime. There are also nebular emission lines and evidence of hot silicate dust in the outer spatial region, centered at approximately 17 AU from the WD, that are not found in the inner regime. The KIN and Spitzer data suggest that these emissions were excited by the nova flash in the outer spatial regime before the blast wave reached these regions. We describe the present results in terms of a new model for dust creation in recurrent novae that includes an increase in density in the plane of the orbit of the two stars created by a spiral shock wave caused by the motion of the stars through the cool wind of the red giant star. These data show the power and potential of the nulling technique which has been developed for the detection of Earth-like planets around nearby stars for the Terrestrial Planet Finder and Darwin missions.

Barry, Richard K.↗
DOE OSTI2021

Distributed and communication-efficient solutions to linear equations with special sparse structure

In this paper we report two distributed and communication-efficient algorithms based on the multi-agent system are proposed to solve a system of linear equations with the Laplacian sparse system matrix. One algorithm is based on the gradient descent method in optimization. In this algorithm, the agents only share partial information instead of all of their collective state vectors to save significant communication. The other algorithm is obtained by approximating Newton’s method for a faster convergence rate. Although it requires twice as much communication as the first one, it is still communication-efficient given the low dimension of the information shared among agents. The convergence at a linear rate is proved for both algorithms, and a comprehensive comparison of their convergence rate, communication burden, and computation costs is also performed. The proposed algorithms can be applied to various systems to solve those problems that can be modeled as a system of linear equations with a Laplacian sparse system matrix. Simulation results with the electric power system illustrate their effectiveness.

42 ENGINEERING↗
NASA NTRS2008

On Combining Thermal-Infrared and Radio-Occultation Data of Saturn's Atmosphere

Radio-occultation and thermal-infrared measurements are complementary investigations for sounding planetary atmospheres. The vertical resolution afforded by radio occultations is typically approximately 1 km or better, whereas that from infrared sounding is often comparable to a scale height. On the other hand, an instrument like CIRS can easily generate global maps of temperature and composition, whereas occultation soundings are usually distributed more sparsely. The starting point for radio-occultation inversions is determining the residual Doppler-shifted frequency, that is the shift in frequency from what it would be in the absence of the atmosphere. Hence the positions and relative velocities of the spacecraft, target atmosphere, and DSN receiving station must be known to high accuracy. It is not surprising that the inversions can be susceptible to sources of systematic errors. Stratospheric temperature profiles on Titan retrieved from Cassini radio occultations were found to be very susceptible to errors in the reconstructed spacecraft velocities (approximately equal to 1 mm/s). Here the ability to adjust the spacecraft ephemeris so that the profiles matched those retrieved from CIRS limb sounding proved to be critical in mitigating this error. A similar procedure can be used for Saturn, although the sensitivity of its retrieved profiles to this type of error seems to be smaller. One issue that has appeared in inverting the Cassini occultations by Saturn is the uncertainty in its equatorial bulge, that is, the shape in its iso-density surfaces at low latitudes. Typically one approximates that surface as a geopotential surface by assuming a barotropic atmosphere. However, the recent controversy in the equatorial winds, i.e., whether they changed between the Voyager (1981) era and later (after 1996) epochs of Cassini and some Hubble observations, has made it difficult to know the exact shape of the surface, and it leads to uncertainties in the retrieved temperature profiles of one to a few kelvins. This propagates into errors in the retrieved helium abundance, which makes use of thermal-infrared spectra and synthetic spectra computed with retrieved radio-occultation temperature profiles. The highest abundances are retrieved with the faster Voyager-era winds, but even these abundances are somewhat smaller than those retrieved from the thermal-infrared data alone (albeit with larger formal errors). The helium abundance determination is most sensitive to temperatures in the upper troposphere. Further progress may include matching the radio-occultation profiles with those from CIRS limb sounding in the upper stratosphere.

Flasar, F. M.↗