Search NASA⌕ Search

SEARCH · Search NASA

Results for “Spherical harmonic analysis”

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
NASA NTRS1971

Lunar gravity derived from long-period satellite motion, a proposed method

A method was devised to determine the spherical harmonic coefficients of the lunar gravity field. The method consists of a two-step data reduction and estimation process. Pseudo-Doppler data were generated simulating two different lunar orbits. The analysis included the perturbing effects of the L1 lunar gravity field, the earth, the sun, and solar radiation pressure. Orbit determinations were performed on these data and long-period orbital elements were obtained. The Kepler element rates from these solutions were used to recover L1 lunar gravity coefficients. Overall results of the experiment show that lunar gravity coefficients can be accurately determined and that the method is dynamically consistent with long-period perturbation theory.

Ferrari, A. J.↗
NASA NTRS1979

Global accuracy estimates of point and mean undulation differences obtained from gravity disturbances, gravity anomalies and potential coefficients

Through the method of truncation functions, the oceanic geoid undulation is divided into two constituents: an inner zone contribution expressed as an integral of surface gravity disturbances over a spherical cap; and an outer zone contribution derived from a finite set of potential harmonic coefficients. Global, average error estimates are formulated for undulation differences, thereby providing accuracies for a relative geoid. The error analysis focuses on the outer zone contribution for which the potential coefficient errors are modeled. The method of computing undulations based on gravity disturbance data for the inner zone is compared to the similar, conventional method which presupposes gravity anomaly data within this zone.

Jekeli, C.↗
NASA NTRS1989

High resolution analysis of satellite gradiometry

Satellite gravity gradiometry is a technique now under development which, by the middle of the next decade, may be used for the high resolution charting from space of the gravity field of the earth and, afterwards, of other planets. Some data analysis schemes are reviewed for getting detailed gravity maps from gradiometry on both a global and a local basis. It also presents estimates of the likely accuracies of such maps, in terms of normalized spherical harmonics expansions, both using gradiometry alone and in combination with data from a Global Positioning System (GPS) receiver carried on the same spacecraft. It compares these accuracies with those of current and future maps obtained from other data (conventional tracking, satellite-satellite tracking, etc.), and also with the spectra of various signals of geophysical interest.

Colombo, O. L.↗
NASA NTRS1994

The isostatic state of Mead crater

We have analyzed high-resolution Magellan Doppler tracking data over Mead crater, using both line-of-sight and spherical harmonic methods, and have found a negative gravity anomaly of about 4-5 mgal (at spacecraft altitude, 182 km). This is consistent with no isostatic compensation of the present topography; the uncertainty in the analysis allows perhaps as much as 30% compensation at shallow dpeths (approximately 25 km). This is similar to observations of large craters on Earth, which are not generally compensated, but contrasts with at least some lunar basins which are inferred to have large Moho uplifts and corresponding positive Bouguer anomalies. An uncompensated load of this size requires a lithosphere with an effective elastic lithosphere thickness greater than 30 km. In order for the crust-mantle boundary not to have participated in the deformation associated with the collapse of the transient cavity during the creation of the crater, the yield strength near the top of the mantle must have been significantly higher on Earth and Venus than on the Moon at the time of basin formation. This might be due to increased strength against frictional sliding at the higher confining pressures within the larger planets. Alternatively, the thinner crusts of Earth and Venus compared to that of the Moon may result in higher creep strength of the upper mantle at shallower depths.

Banerdt, W. B.↗
NASA NTRS2003

Measurement and Interpretation of Temporal Variations of the Earths Gravity Field Using GPS and SLR Data

This investigation has had four main thrusts: 1) The analysis of seasonal variations of the Earth's gravitational field using Lageos 1 and 2 SLR data and comparisons to geophysical models. We have estimated the annual variation of the gravity field via a spherical harmonic expansion complete to degree and order 4. We have also constructed a similar model using models of the annual variation in the gravity field due to atmospheric, hydrologic, and ocean mass redistribution. These three models, when combined together, are in excellent agreement with the variations observed by satellite laser ranging. An article on these results was published in the journal Geophysical Research Letters. 2) The second thrust of our investigation has been to analyze the output of a Global Climate Model (GCM) to determine if the GRACE gravity mission can be expected to detect climate change signals. Working with Gary Russell at the Goddard Institute for Space Studies (GISS), we have determined that there are several large secular signals that GRACE might be able to detect, including secular changes in snow cover, sea ice, polar ice, ocean mass, and other variables. It is possible that some of these signals could be detected with 5 years of GRACE measurements - its hard to judge this because the interannual variability in the GCM, which could mask the climate signals, is unreliable. Certainly a follow-on GRACE mission could detect these signals when compared to the data from the initial GRACE mission.). An article on these results will be published in the journal Journal of Geophysical Research. 3) In the last year of the investigation, we developed a new technique for analyzing temporal gravity variations using "geophysical fingerprints", which was successfully demonstrated on 20 years of satellite laser ranging data [Nerem et al., 20031. 4]. We also participated in a workshop on future satellite gravity measurements, which resulted in paper on measuring ocean mass variations using GRACE [Nerem et al., 20031 and on using laser interferometry for future gravity missions [Bender et al., 20031].

Nerem, R. Steven↗
NASA NTRS2017

B-plane Evolution Under Highly Non-Keplerian Dynamics

Juno is the second of a series of New Frontiers missions and was launched in 2011. The spacecraft was set en route to Jupiter, and its Jupiter Orbit Insertion occurred on 2016-07-04 after a five-year cruise in deep space. The mission phase just prior to this large maneuver is called approach. As navigation team, we were responsible for precise orbit determination of the Juno spacecraft to ensure a successful orbit insertion. To evaluate the navigation performance we employed the B-plane mapping, where the incoming hyperbolic velocity defines the plane perpendicular to it, and the time to go to hit the plane will directly map into the uncertainties of the target point. The B-plane is a convenient tool when assessing the navigation performance as it is not affected by non-linearities due to the gravitational pull of the targeted body. However, as the Jupiter Orbit Insertion approached we started to realize that the B-plane mapping uncertainties were significantly larger than expected. We found that the B-plane mapping time had tremendous influence on the covariance inflation in the B-plane. Because of the close approach distance during the Jupiter Orbit Insertion and the strong interaction with J2 spherical harmonic coefficient, the Juno dynamics were too far from a Keplerian one assumed for the B-plane mapping. We first discuss the analytical approach of the B-plane mapping uncertainty and perform a numerical analysis to isolate the components of the Juno dynamics that cause the B-plane covariance inflation.

Bordi, John↗
NASA NTRS1993

Free-air and Bouguer gravity anomalies and the Martian crustal dichotomy

Free-air and Bouguer gravity anomalies from a 50x50 field, derived from re-analysis of Viking Orbiter and Mariner 9 tracking data and using a 50x50 expansion of the current Mars topography and the GSFC degree 50 geoid as the equipotential reference surface, with the Martian crustal dichotomy are compared. The spherical harmonic topography used has zero mean elevation, and differs from the USGS maps by about 2 km. In this field the dichotomy boundary in eastern Mars lies mostly at -1 to -2 km elevation. Bouguer gravity anomalies are shown on a map of Noachian, Hesperian, and Amazonian age terrains, simplified from current geologic maps. The map is centered at 300 deg W to show the continuity of the dichotomy boundary. Contour interval is 100 mgals. Gravity and topography were compared along approximately 40 profiles oriented parallel to the dichotomy boundary topographic gradient, to determine how the geophysical character of the boundary changes along its length and what this implies for its origin and development.

Frey, Herbert↗
NASA NTRS1980

Deconvolution and analysis of wide-angle longwave radiation data from Nimbus 6 Earth radiation budget experiment for the first year

One year of longwave radiation data from July 1975 through June 1976 from the Nimbus 6 satellite Earth radiation budget experiment is analyzed by representing the radiation field by a spherical harmonic expansion. The data are from the wide field of view instrument. Contour maps of the longwave radiation field and spherical harmonic coefficients to degree 12 and order 12 are presented for a 12 month data period.

Bess, T. D.↗
NASA NTRS1992

Analysis of crustal structure of Venus utilizing residual Line-of-Sight (LOS) gravity acceleration and surface topography data. A trial of global modeling of Venus gravity field using harmonic spline method

To construct Venus' gravity disturbance field (or gravity anomaly) with the spacecraft-observer line of site (LOS) acceleration perturbation data, both a global and a local approach can be used. The global approach, e.g., spherical harmonic coefficients, and the local approach, e.g., the integral operator method, based on geodetic techniques are generally not the same, so that they must be used separately for mapping long wavelength features and short wavelength features. Harmonic spline, as an interpolation and extrapolation technique, is intrinsically flexible to both global and local mapping of a potential field. Theoretically, it preserves the information of the potential field up to the bound by sampling theorem regardless of whether it is global or local mapping, and is never bothered with truncation errors. The improvement of harmonic spline methodology for global mapping is reported. New basis functions, a singular value decomposition (SVD) based modification to Parker & Shure's numerical procedure, and preliminary results are presented.

Fang, Ming↗
NASA NTRS1976

Spherical harmonic representation of the gravity field of Mars using a short-arc technique

Planetary gravity fields are primarily determined from Doppler tracking of orbiting spacecraft which are typically in high eccentricity orbits. Thus, the non-central gravitational components perturb the orbit primarily in the vicinity of periapsis. To find a computationally efficient technique and to reduce the influence of unmodeled forces on the spacecraft, an analysis has been performed of the effectiveness of processing only data in the vicinity of periapsis using multiple short arcs. Mariner 9 data are processed for the study, using 1, 2, 3 and 4 hour arcs and estimating gravity fields from 4th through 7th degree. The resulting gravity fields are inter-compared and also compared with previously published results.

Daniels, E. F.↗
NASA NTRS1995

Upper Thermosphere Winds and Temperatures in the Geomagnetic Polar Cap: Solar Cycle, Geomagnetic Activity, and Interplanetary Magnetic Field Dependencies

Ground-based Fabry-Perot interferometers located at Thule, Greenland (76.5 deg. N, 69.0 deg. W, lambda = 86 deg.) and at Sondre Stromfjord, Greenland (67.0 deg. N, 50.9 deg. W, lambda = 74 deg.) have monitored the upper thermospheric (approx. 240-km altitude) neutral wind and temperature over the northern hemisphere geomagnetic polar cap since 1983 and 1985, respectively. The thermospheric observations are obtained by determining the Doppler characteristics of the (OI) 15,867-K (630.0-nm) emission of atomic oxygen. The instruments operate on a routine, automatic, (mostly) untended basis during the winter observing seasons, with data coverage limited only by cloud cover and (occasional) instrument failures. This unique database of geomagnetic polar cap measurements now extends over the complete range of solar activity. We present an analysis of the measurements made between 1985 (near solar minimum) and 1991 (near solar maximum), as part of a long-term study of geomagnetic polar cap thermospheric climatology. The measurements from a total of 902 nights of observations are compared with the predictions of two semiempirical models: the Vector Spherical Harmonic (VSH) model of Killeen et al. (1987) and the Horizontal Wind Model (HWM) of Hedin et al. (1991). The results are also analyzed using calculations of thermospheric momentum forcing terms from the Thermosphere-ionosphere General Circulation Model TGCM) of the National Center for Atmospheric Research (NCAR). The experimental results show that upper thermospheric winds in the geomagnetic polar cap have a fundamental diurnal character, with typical wind speeds of about 200 m/s at solar minimum, rising to up to about 800 m/s at solar maximum, depending on geomagnetic activity level. These winds generally blow in the antisunward direction, but are interrupted by episodes of modified wind velocity and altered direction often associated with changes in the orientation of the Interplanetary Magnetic Field (IMF). The central polar cap (greater than approx. 80 magnetic latitude) antisunward wind speed is found to be a strong function of both solar and geomagnetic activity. The polar cap temperatures show variations in both solar and geomagnetic activity, with temperatures near 800 K for low K(sub p) and F(sub 10.7) and greater than about 2000 K for high K(sub p) and F(sub 10.7). The observed temperatures are significantly greater than those predicted by the mass spectrometer/incoherent scatter model for high activity conditions. Theoretical analysis based on the NCAR TIGCM indicates that the antisunward upper thermospheric winds, driven by upstream ion drag, basically 'coast' across the polar cap. The relatively small changes in wind velocity and direction within the polar cap are induced by a combination of forcing terms of commensurate magnitude, including the nonlinear advection term, the Coriolis term, and the pressure gradient force term. The polar cap thennospheric thermal balance is dominated by horizontal advection, and adiabatic and thermal conduction terms.

Killeen, T. L.↗
NASA NTRS1990

Gravity field improvement using GPS data from Topex/Poseidon - A covariance analysis

A covariance analysis is performed using a realistic scenario for processing 10 days of GPS data, to obtain the expected improvement to the GEM-T2 gravity field. The gravity bin technique has been refined to compute the covariance matrix associated with the spherical harmonic gravity field. It is shown that the GPS data from one ten-day arc of Topex/Poseidon with no a priori can improve medium degree and order (3-26) sigmas for the parameters in the GEM-T2 gravity field by more than an order of magnitude.

Bertiger, Willy I.↗
NASA NTRS2011

Rapid Calculation of Spacecraft Trajectories Using Efficient Taylor Series Integration

A variable-order, variable-step Taylor series integration algorithm was implemented in NASA Glenn's SNAP (Spacecraft N-body Analysis Program) code. SNAP is a high-fidelity trajectory propagation program that can propagate the trajectory of a spacecraft about virtually any body in the solar system. The Taylor series algorithm's very high order accuracy and excellent stability properties lead to large reductions in computer time relative to the code's existing 8th order Runge-Kutta scheme. Head-to-head comparison on near-Earth, lunar, Mars, and Europa missions showed that Taylor series integration is 15.8 times faster than Runge- Kutta on average, and is more accurate. These speedups were obtained for calculations involving central body, other body, thrust, and drag forces. Similar speedups have been obtained for calculations that include J2 spherical harmonic for central body gravitation. The algorithm includes a step size selection method that directly calculates the step size and never requires a repeat step. High-order Taylor series integration algorithms have been shown to provide major reductions in computer time over conventional integration methods in numerous scientific applications. The objective here was to directly implement Taylor series integration in an existing trajectory analysis code and demonstrate that large reductions in computer time (order of magnitude) could be achieved while simultaneously maintaining high accuracy. This software greatly accelerates the calculation of spacecraft trajectories. At each time level, the spacecraft position, velocity, and mass are expanded in a high-order Taylor series whose coefficients are obtained through efficient differentiation arithmetic. This makes it possible to take very large time steps at minimal cost, resulting in large savings in computer time. The Taylor series algorithm is implemented primarily through three subroutines: (1) a driver routine that automatically introduces auxiliary variables and sets up initial conditions and integrates; (2) a routine that calculates system reduced derivatives using recurrence relations for quotients and products; and (3) a routine that determines the step size and sums the series. The order of accuracy used in a trajectory calculation is arbitrary and can be set by the user. The algorithm directly calculates the motion of other planetary bodies and does not require ephemeris files (except to start the calculation). The code also runs with Taylor series and Runge-Kutta used interchangeably for different phases of a mission.

Scott, James R.↗
NASA NTRS2017

Drag-Free Control and Drag Force Recovery of Small Satellites

Drag-free satellites provide autonomous precision orbit determination, accurately map the static and time varying components of Earth's mass distribution, aid in our understanding of the fundamental force of gravity, and will ultimately open up a new window to our universe through the detection and observation of gravitational waves. At the heart of this technology is a gravitational reference sensor, which (a) contains and shields a free-floating proof mass from all non-gravitational forces, and (b) precisely measures the position of the test mass inside the sensor. Thus, both test mass and spacecraft follow a pure geodesic in spacetime. By tracking the position of a low Earth orbiting drag-free satellite we can directly determine the detailed shape of geodesics and through analysis, the higher order harmonics of the Earths geopotential. This paper explores two different drag-free control systems on small satellites. The first drag-free control system is a continuously compensated single thruster 3-unit CubeSat with a suspension-free spherical proof-mass. A feedback control system commands the thruster and Attitude and Determination Control System to fly the tender spacecraft with respect to the test mass. The spheres position is sensed with a LED-based differential optical shadow sensor, its electric charge controlled by photoemission using UV LEDs, and the spacecraft position is maintained with respect to the sphere using an ion electrospray propulsion system. This configuration is the most fuel-efficient drag-free system possible today. The second drag-free control system is an electro-statically suspended cubical proof-mass that is operated with a low duty cycle, limiting suspension force noise over brief, known time intervals on a small GRACE-II -like satellite. The readout is performed using a laser interferometer, which is immune to the dynamic range limitations of voltage references. This system eliminates the need for a thruster, enabling drag-free control systems for passive satellites. In both cases, the test mass position, GPS tracking data, and commanded actuation, either thrust or suspension system, can be analyzed to estimate the 3-axis drag forces acting on the satellite. The data produces the most precise maps of upper atmospheric drag forces and with additional information, detailed models that describe the dynamics of the upper atmosphere and its impact on all satellites that orbit the Earth. This paper highlights the history, applications, design, laboratory technology development and highly detailed simulation results of each control system.

Small Satellites↗
NASA NTRS1975

The complex magnetic field of Jupiter

An analysis of the characteristics of the magnetic field of the planet Jupiter is presented. The data were obtained during the flight of Pioneer 11 space probe, using a high field triaxial fluxgate magnetometer. The data are analyzed in terms of traditional Schmitt normalized spherical harmonic expansion fitted to the observations in a least squares sense. Tables of data and graphs are provided to summarize the findings.

Acuna, M. H.↗
NASA NTRS1982

An analysis of the structure of Saturn's magnetic field using charged particle absorption signatures

A new technique is derived for determining the structure of Saturn's magnetic field. This technique uses the observed positions of charged particle signatures due to the satellites and rings of Saturn to determine the parameters of an axially symmetric, spherical harmonic model of the magnetic field using the method of least squares. Absorption signatures observed along the Pioneer 11, Voyager 1, and Voyager 2 spacecraft trajectories are used to derive values for the orientation of the magnetic symmetry axis relative to Saturn's axis of rotation, the axial displacement of the center of the magnetic dipole from the center of Saturn, and the magnitude of the external field component. Comparing these results with the magnetic field model parameters deduced from analyses of magnetometer data leads to a preference for models that incorporate a northward offset of the dipole center by about 0.05 Saturn radii.

Chenette, D. L.↗
NASA NTRS1993

Theory and operation of a three-gate time-of-flight velocity analyzer

Theoretical considerations and test results are presented for a new-type velocity analyzer for incident fast neutral particles, positive ions, and negative ions. Velocity analysis is carried out by means of a pulsed, three-gate time-of-flight (TOF) technique capable of eliminating alias velocities (harmonics) to sixth order. In addition the design and operation are presented of a four-element ion lens system, with small spherical and chromatic aberrations, suitable for interfacing a large-diameter ion beam from the TOF section with a subsequent mass analyzer.

Martus, K. E.↗