SEARCH · Search NASA
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.
To assess the resolution of the local gravity field from kinematic measurements, a state model for motion in the gravity field of the earth is formulated. The resulting set of equations can accommodate gravity gradients, specific force, acceleration, velocity and position as input data and can take into account approximation errors as well as sensor errors.
The planetary ephemerides approximation for radar astronomy is discussed, and, in particular, the effect of this approximation on the performance of the programmable local oscillator (PLO) used in Goldstone Solar System Radar is presented. Four different approaches are considered and it is shown that the Gram polynomials outperform the commonly used technique based on Chebyshev polynomials. These methods are used to analyze the mean square, the phase error, and the frequency tracking error in the presence of the worst case Doppler shift that one may encounter within the solar system. It is shown that in the worst case the phase error is under one degree and the frequency tracking error less than one hertz when the frequency to the PLO is updated every millisecond.
To optimize the performance of any system, the sensitivity derivatives of the system's output variables with respect to its input variables must be readily available. It is also desirable that these derivatives be inexpensive to calculate as the optimization process requires many evaluations of the output variables and their derivatives. Optimization methods that have been developed for use in automated structural design programs may not be extended for use in integrated multidisciplinary design programs until adequate means of calculating accurate sensitivity derivatives of complex, internally coupled systems have been developed. Until the development of the Global Sensitivity Equations (GSE), the only method of determining the sensitivity derivatives of coupled systems has been by using finite differences. Analytical or semi-analytical derivatives do not exist as there is no analytical solution to the coupled problem. Also, difficulties arise because the finite difference method is expensive as the system has to iterate to a converged solution for each incremental input variable. The method may not be accurate, and the choice of the input variable increment may cause the difference in the output variable to be insignificant compared to computer numerical error if the choice is too small, or the process may not predict the true value of the output variable if the increment is too large. The GSE allow the system's sensitivity derivatives to be calculated as functions of the component subsystem's (local) sensitivity derivatives. These local sensitivity derivatives are calculated from specifically decoupled subsystems, whereas the GSE account for total system coupling. Since the subsystems are decoupled, it may be possible for the local derivatives to be calculated by analytical or semi-analytical methods, which generally reduce cost and improve accuracy. Several academic problems have been solved using GSE and have demonstrated encouraging results.
This dissertation addresses various issues for model classes of hyperbolic conservation laws. The basic approach developed in this work employs a new family of adaptive, hp-version, finite element methods based on a special discontinuous Galerkin formulation for hyperbolic problems. The discontinuous Galerkin formulation admits high-order local approximations on domains of quite general geometry, while providing a natural framework for finite element approximations and for theoretical developments. The use of hp-versions of the finite element method makes possible exponentially convergent schemes with very high accuracies in certain cases; the use of adaptive hp-schemes allows h-refinement in regions of low regularity and p-enrichment to deliver high accuracy, while keeping problem sizes manageable and dramatically smaller than many conventional approaches. The use of discontinuous Galerkin methods is uncommon in applications, but the methods rest on a reasonable mathematical basis for low-order cases and has local approximation features that can be exploited to produce very efficient schemes, especially in a parallel, multiprocessor environment. The place of this work is to first and primarily focus on a model class of linear hyperbolic conservation laws for which concrete mathematical results, methodologies, error estimates, convergence criteria, and parallel adaptive strategies can be developed, and to then briefly explore some extensions to more general cases. Next, we provide preliminaries to the study and a review of some aspects of the theory of hyperbolic conservation laws. We also provide a review of relevant literature on this subject and on the numerical analysis of these types of problems.
We present measurements of the column densities of interstellar D I and O I made with the Far Ultraviolet Spectroscopic Explorer (FUSE), and of H I made with the International Ultraviolet Explorer (IUE) toward the sdOB star Feige 110 [(l,b) = (74.09 deg., - 59.07 deg.); d = 179(sup +265, sub -67) pc; Z = -154(sup +57, Sub -227 pc). Our determination of the D I column density made use of curve of growth fitting and profile fitting analyses, while our O I column density determination used only curve of growth techniques. The H I column density was estimated by fitting the damping wings of the interstellar Ly(lpha) profile. We find log N(D I) = 15.47 +/- 0.06, log N(O I) = 16.73 +/- 0.10, and log N(H I) = 20.14(sup +0.13, sub -0.20) (all errors 2(sigma)). This implies D/H = (2.14 +/- 0.82) x 10(esp -5), D/O = (5.50(sup + 1.64, sub -133)) x 10(exp -2), and O/H = (3.89 +/- 1.67) x 10(exp -4). Taken with the FUSE results reported in companion papers and previous measurements of the local interstellar medium, this suggests the possibility of spatial variability in D/H for sight lines exceeding approx. 100 pc. This result may constrain models which characterize the mixing time and length scales of material in the local interstellar medium.
The Nimbus-7 and Earth Probe Total Ozone Mapping Spectrometer (TOMS) data have been reprocessed with a new retrieval algorithm, (Version 8) and an updated calibration procedure. These data have been systematically compared to total ozone data from Brewer and Dobson spectrophotometers for 73 individual ground stations. The comparisons were made as a function of latitude, solar zenith angle, reflectivity and total ozone. Results show that the accuracy of the TOMS retrieval'is much improved when aerosols are present in the atmosphere, when snow/ice and sea glint are present, and when ozone in the northern hemisphere is extremely low. TOMS overpass data are derived from the single TOMS best match measurement, almost always located within one degree of the ground station and usually made within an hour of local noon. The version 8 Earth Probe TOMS ozone values have decreased by an average of about 1% due to a much better understanding of the calibration of the instrument. The remaining differences between TOMS and ground stations suggest that there are still small errors in the TOMS retrievals. But if TOMS is used as a transfer standard to compare ground stations, the large station-to-station differences suggest the possibility of significant instrument errors at some ground stations.
Originally developed for the Moderate Resolution Imaging Spectroradiometer (MODIS) in polar, sun-synchronous low-earth orbit (LEO), the Dark Target (DT) aerosol retrieval algorithm relies on the assumption of a Surface Reflectance Parameterization (SRP) over land surfaces. Specifically for vegetated and dark-soiled surfaces, values of surface reflectance in blue and red visible-wavelength bands are assumed to be nearly linearly related to each other and to the value in a shortwave infrared (SWIR) wavelength band. This SRP also includes dependencies on scattering angle and a normalized difference vegetation index computed from two SWIR bands (NDVISWIR). As the DT retrieval algorithm is being ported to new sensors to continue and expand the aerosol data record, we assess whether the MODIS-assumed SRP can be used for these sensors. Here, we specifically assess SRP for the Advanced Baseline Imager (ABI) aboard, the Geostationary Operational Environmental Satellite (GOES)-16/East (ABIE). First, we find that using MODIS-based SRP leads to higher biases and artificial diurnal signatures in aerosol optical depth (AOD)retrievals from ABIE. The primary reason appears to be that geostationary orbit (GEO) encounters an entirely different set of observation geometry than does LEO, primarily with regards to solar angles coupled with fixed view angles. Therefore, we have developed a new SRP for GEO that draws the angular shape of the surface bidirectional reflectance. We also introduce modifications to the parametrization of both red-SWIR and blue-red spectral relationships to include additional information.The revised Red-SWIR SRP includes solar zenith angle, NDVISWIR, and land-type percentage from an ancillary database. The blue-red SRP adds dependencies on the scattering angle and NDVISWIR. The new SRPs improve the AOD retrieval of ABIE in terms of overall less bias and mitigation of the overestimation around local noon. The average bias of DT AOD compared to AERONET AOD shows a reduction from 0.08 to 0.03, while the bias of local solar noon decreases from 0.12 to 0.03.The agreement between DT and AERONET AOD is established through regression slope of 1.06 and y-intercept of 0.01with correlation coefficient is0.74. By using the new SRP, the percentage of data falling within the expected error range (±0.05 + 15%) is notably risen from 54% to 78%.
A two-step advection scheme of the Lax-Wendroff type is derived which has accuracy and phase characteristics similar to that of a third-order scheme. The scheme is exactly third-order accurate in time and space for uniform flow. The new scheme is compared with other currently used methods, and is shown to simulate well the advection of localized disturbances with steep gradients. The scheme is derived for constant flow and generalized to two-dimensional nonuniform flow.
The investigation was focused on studying the residual X-ray emission from possible counterparts of enigmatic gamma-ray transients near the Galactic plane. These sources (discovered by our group at Columbia) are highly variable gamma-ray sources with no radio-loud spectrally flat blazar in their approx 1 deg. diameter error box. BSAX Narrow Field Instruments pointed the gamma-ray error box of the prominent non-blazar Energetic Gamma Ray Experiment Telescope (EGRET) source GRO 1838-04 twice on April 17, 1997, and on April 20, 1997. Most of the EGRET error box was imaged and several faint X-ray sources were detected. BSAX Wide Field Camera data were also obtained for a variety of sources near the Galactic plane, including GRO J1838-04.
Flames in a Mach 0.3 atmospheric pressure laboratory burner rig were doped with sea salt, NaS04, and NaCl, respectively, in an effort to validate theoretical dew point predictions made by a local thermochemical equilibrium (LTCE) method of predicting condensation temperatures of sodium sulfate in flame environments. Deposits were collected on cylindrical platinum targets placed in the combustion products, and the deposition was studied as a function of collector temperature. Experimental deposition onset temperatures checked within experimental error with LTCE-predicted temperatures. A multicomponent mass transfer equation was developed to predict the rate of deposition of Na2SO4(c) via vapor transport at temperatures below the deposition onset temperature. Agreement between maximum deposition rates predicted by this chemically frozen boundary layer (CFBL) theory and those obtained in the seeded laboratory burner experiments is good.
The study has three goals. The first is to demonstrate that processes exist whose intrinsic nonlocal nature cannot be represented by local approximations. The second is to elucidate the physical nature and origins of these nonlocal processes. The third is to suggest that the methods and results described here may prove useful in constructing semiempirical models of the chromosphere by means more efficient than trial and error. Matrices are computed that describe the effect of a temperature perturbation at an arbitrary point in the loop on density, hydrogen ionized fraction, total radiative loss rate, and radiative loss rate of selected hydrogen lines and continua at all other points. It is found that the dominant nonlocal radiative transfer effects can be separated into flux divergence coefficient effects and upper level population effects. The former are most important when the perturbation takes place in a region of significant opacity. Upper level population effects arise in both optically thick and thin regions in response to nonlocal density, ionization, and interlocking effects.
A method for obtaining field enhancement factors at specific points on an aircraft utilizing a small scale model was evaluated by measuring several canonical shapes. Comparison of the form factors obtained by analytical means with measurements indicate that the experimental method has serious flaws. Errors of 200 to 300 percent were found between analytical values and measured values. As a result of the study, the analytical method is not recommended for calibration of field meters located on aircraft, and should not be relied upon in any application where the local spatial derivatives of the electric field on the model are large over the dimensions of the sensing probe.
The advantages and disadvantages of three different approaches to solving the problem of the radiometric correction of synthetic aperture radar (SAR) images of varying terrain heights are presented. The first approach involves registration of a digital elevation model (DEM) of the terrain to the image, determination of the local elevation and incidence angles, and appropriate radiometric correction. The second approach uses a DEM generated from interferometric SAR data to derive the elevation and incidence angle maps. In the third approach, a monopulse technique is employed to determine the elevation angle only. The relative errors in radiometric correction between these approaches are assessed. Calibration errors are estimated using corner reflectors deployed within some of the scenes imaged by the Jet Propulsion Laboratory airborne SAR (JPL AIRSAR).
To improve snowpack estimates in Community Land Model version 4 (CLM4), the Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover fraction (SCF) was assimilated into the Community Land Model version 4 (CLM4) via the Data Assimilation Research Testbed (DART). The interface between CLM4 and DART is a flexible, extensible approach to land surface data assimilation. This data assimilation system has a large ensemble (80-member) atmospheric forcing that facilitates ensemble-based land data assimilation. We use 40 randomly chosen forcing members to drive 40 CLM members as a compromise between computational cost and the data assimilation performance. The localization distance, a parameter in DART, was tuned to optimize the data assimilation performance at the global scale. Snow water equivalent (SWE) and snow depth are adjusted via the ensemble adjustment Kalman filter, particularly in regions with large SCF variability. The root-mean-square error of the forecast SCF against MODIS SCF is largely reduced. In DJF (December-January-February), the discrepancy between MODIS and CLM4 is broadly ameliorated in the lower-middle latitudes (2345N). Only minimal modifications are made in the higher-middle (4566N) and high latitudes, part of which is due to the agreement between model and observation when snow cover is nearly 100. In some regions it also reveals that CLM4-modeled snow cover lacks heterogeneous features compared to MODIS. In MAM (March-April-May), adjustments to snowmove poleward mainly due to the northward movement of the snowline (i.e., where largest SCF uncertainty is and SCF assimilation has the greatest impact). The effectiveness of data assimilation also varies with vegetation types, with mixed performance over forest regions and consistently good performance over grass, which can partly be explained by the linearity of the relationship between SCF and SWE in the model ensembles. The updated snow depth was compared to the Canadian Meteorological Center (CMC) data. Differences between CMC and CLM4 are generally reduced in densely monitored regions.
Reflected shortwave (SW) solar radiations at the top of atmosphere from Clouds and the Earth’s Radiant Energy System (CERES), Modern Era-Retrospective analysis for Research and Applications version 2 (MERRA-2), and ECMWF Reanalysis 5th Generation (ERA5) are examined to better under-stand their differences in spatial and temporal variations (seasonal and diurnal cycle time-scale) with respect to the observations from the Earth Polychromatic Imaging Camera (EPIC) on Deep Space Climate Observatory (DSCOVR) satellite. Comparisons between two reanalyses (MERRA-2 and ERA5) and EPIC reveal that MERRA-2 has a generally larger deviation from EPIC than ERA5, in terms of the SW radiance and diurnal variability in all seasons, which can be attributed to larger cloud biases in MERRA-2. MERRA-2 produces more ice/liquid water content than ERA5 over the tropical warm pool, leading to positive SW biases in cloud and radiance, while both reanalyses underestimate the observed SW radiance from EPIC in the stratus-topped region off the western coast of US/Mexico in the boreal summer. Himalaya/Tibet region in the boreal spring/summer and the midlatitude Southern Hemisphere in the boreal winter are the regions where MERRA-2 and ERA5 deviate largely from EPIC but their deviations have the opposite sign. Vertical structures of cloud ice/liquid water content explain reasonably well these contrasting differences between the two reanalyses. As two independent observa-ions, CERES and EPIC agree well with each other in terms of the SW radiance maps, showing 2-3% mean absolute errors over the tropical-midlatitudes. The CERES-EPIC consistency further confirms that the reanalyses still have challenges in representing the SW flux and its global distribution. In the CERES-EPIC observation differences, CERES slightly overestimates the diurnal cycle (as a function of local solar time) of the observed EPIC irradiance in the morning and underestimates in the after-noon, while the opposite is the case in the reanalyses.
Object localization and its application in tele-autonomous systems are studied. Two object localization algorithms are presented together with the methods of extracting several important types of object features. The first algorithm is based on line-segment to line-segment matching. Line range sensors are used to extract line-segment features from an object. The extracted features are matched to corresponding model features to compute the location of the object. The inputs of the second algorithm are not limited only to the line features. Featured points (point to point matching) and featured unit direction vectors (vector to vector matching) can also be used as the inputs of the algorithm, and there is no upper limit on the number of the features inputed. The algorithm will allow the use of redundant features to find a better solution. The algorithm uses dual number quaternions to represent the position and orientation of an object and uses the least squares optimization method to find an optimal solution for the object's location. The advantage of using this representation is that the method solves for the location estimation by minimizing a single cost function associated with the sum of the orientation and position errors and thus has a better performance on the estimation, both in accuracy and speed, than that of other similar algorithms. The difficulties when the operator is controlling a remote robot to perform manipulation tasks are also discussed. The main problems facing the operator are time delays on the signal transmission and the uncertainties of the remote environment. How object localization techniques can be used together with other techniques such as predictor display and time desynchronization to help to overcome these difficulties are then discussed.
An analytical technique is presented for the deriving Mars atmosphere opacity from data furnished by the IR Thermal Mapper (IRTM) instruments on the Viking landers. Each IRTM features four telescopes for viewing Mars in six spectral bands and seven spatial fields of view. The present study focuses on 20 microns measurements of the thermal inertia, 0.3-3.0 microns data for the surface albedo, and, finally, the 7, 9, and 15 microns data for the actual opacity derivations. The model takes into account the atmospheric temperature profile, scattering by dust, the surface emissivity and the variation of the surface/atmosphere temperature discontinuity. Analyses are carried out of the possible error factors and the model is used to generate opacity histories for the two Viking Lander sites. The results, when compared with other astronomical data, indicate that the method is sufficiently accurate to use opacity data to analyze the processes of storm genesis and to study local storms.
We present precise localizations of AXP 1E184-045 and SGR 1627-41 with Chandra. We obtained new infrared observations of SGR 1627-41 and reanalyzed archival observations of AXP 1E1841-045 in order to refine their positions and search for infrared counterparts. A faint source is detected inside the error circle of AXP 1E1841-045. In the case of SGR 1627-41, several sources are located within the error radius of the X-ray position, and we discuss the likelihood of one of them being the counterpart. We compare the properties of our candidates to those of other known anomalous X-ray pulsar (AXP) and soft gamma repeater (SGR) counterparts. We find that the counterpart candidates for SGR 1627-41 and SGR 1806-20 would have to be intrinsically much brighter than AXPs in order to have counterparts detectable with the observational limits currently available for these sources. To confirm the reported counterpart of SGR 1806-20, we obtained new infrared observations during the 2003 July burst activation of the source. No brightening of the suggested counterpart is detected, implying that the counterpart of SGR 1806-20 remains yet to be identified.