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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Study of geopotential error models used in orbit determination error analysis

The uncertainty in the geopotential model is currently one of the major error sources in the orbit determination of low-altitude Earth-orbiting spacecraft. The results of an investigation of different geopotential error models and modeling approaches currently used for operational orbit error analysis support at the Goddard Space Flight Center (GSFC) are presented, with emphasis placed on sequential orbit error analysis using a Kalman filtering algorithm. Several geopotential models, known as the Goddard Earth Models (GEMs), were developed and used at GSFC for orbit determination. The errors in the geopotential models arise from the truncation errors that result from the omission of higher order terms (omission errors) and the errors in the spherical harmonic coefficients themselves (commission errors). At GSFC, two error modeling approaches were operationally used to analyze the effects of geopotential uncertainties on the accuracy of spacecraft orbit determination - the lumped error modeling and uncorrelated error modeling. The lumped error modeling approach computes the orbit determination errors on the basis of either the calibrated standard deviations of a geopotential model's coefficients or the weighted difference between two independently derived geopotential models. The uncorrelated error modeling approach treats the errors in the individual spherical harmonic components as uncorrelated error sources and computes the aggregate effect using a combination of individual coefficient effects. This study assesses the reasonableness of the two error modeling approaches in terms of global error distribution characteristics and orbit error analysis results. Specifically, this study presents the global distribution of geopotential acceleration errors for several gravity error models and assesses the orbit determination errors resulting from these error models for three types of spacecraft - the Gamma Ray Observatory, the Ocean Topography Experiment, and the Cosmic Background Explorer.

Yee, C.↗

NASA NTRS2014

High-resolution Local Gravity Model of the South Pole of the Moon from GRAIL Extended Mission Data

We estimated a high-resolution local gravity field model over the south pole of the Moon using data from the Gravity Recovery and Interior Laboratory's extended mission. Our solution consists of adjustments with respect to a global model expressed in spherical harmonics. The adjustments are expressed as gridded gravity anomalies with a resolution of 1/6deg by 1/6deg (equivalent to that of a degree and order 1080 model in spherical harmonics), covering a cap over the south pole with a radius of 40deg. The gravity anomalies have been estimated from a short-arc analysis using only Ka-band range-rate (KBRR) data over the area of interest. We apply a neighbor-smoothing constraint to our solution. Our local model removes striping present in the global model; it reduces the misfit to the KBRR data and improves correlations with topography to higher degrees than current global models.

South Pole↗

NASA NTRS1996

Seasonal Variations of the Earth's Gravitational Field: An Analysis of Atmospheric Pressure, Ocean Tidal, and Surface Water Excitation

Monthly mean gravitational field parameters (denoted here as C(sub even)) that represent linear combinations of the primarily even degree zonal spherical harmonic coefficients of the Earth's gravitational field have been recovered using LAGEOS I data and are compared with those derived from gridded global surface pressure data of the National meteorological center (NMC) spanning 1983-1992. The effect of equilibrium ocean tides and surface water variations are also considered. Atmospheric pressure and surface water fluctuations are shown to be the dominant cause of observed annual C(sub even) variations. Closure with observations is seen at the 1sigma level when atmospheric pressure, ocean tide and surface water effects are include. Equilibrium ocean tides are shown to be the main source of excitation at the semiannual period with closure at the 1sigma level seen when both atmospheric pressure and ocean tide effects are included. The inverted barometer (IB) case is shown to give the best agreement with the observation series. The potential of the observed C(sub even) variations for monitoring mass variations in the polar regions of the Earth and the effect of the land-ocean mask in the IB calculation are discussed.

Dong, D,↗

NASA NTRS1993

An inversion of geoid and topography for mantle and crustal structure on Mars

Mars has the largest amplitude geoid anomalies and surface topography known on the terrestrial planets. A number of prior studies have analyzed Martian gravity anomalies and topography in terms of isostasy and flexure of the crust and lithosphere. Other studies have emphasized the role of mantle convection in producing gravity anomalies and topography in some regions of Mars. Geoid and topography observations for simultaneous estimates of density anomalies in the crust and mantle of Mars are inverted. In performing this study, a recent degree 50 spherical harmonic expansion of the Martian gravity field (GMM-l) and a corresponding resolution expansion of the USGS Mars topography model are used. However, our analysis is restricted to harmonic degrees up to L equals 25, which are better determined than the higher harmonics. This provides a half-wavelength horizontal resolution of 425 km.

Kiefer, Walter↗

NASA NTRS1993

Mars Gravity Field From Dual Satellite Observations

Two strategies are examined in this paper to improve the spherical harmonic representations of the Martian gravitational field. Continous observation over a week-long data arc and a 30th degree gravitational field are assumed for analysis.

Mars↗

NASA NTRS1999

The Lunar Crustal Thickness from Analysis of the Lunar Prospector Gravity and Clementine Topography Datasets

The Lunar Prospector spacecraft has mapped the gravity field of the Moon to a level of resolution never achieved before, and a spherical harmonic representation to degree and order 100 is available. When combined with the topography dataset produced by the Clementine mission, the resulting Bouguer anomaly map is interpreted to model the thickness of the lunar crust. Such models are crucial to understanding the lunar thermal history and the formation of geological features such as mascon basins, several more of which have been newly discovered from this dataset. A two-layer planetary model was used to compute the variations of the depth to the lunar Moho. The thickness values ranged from near 0 to 120 km. There is significant agreement with previous work using the Clementine gravitational field data with differences in specific locations such as South Pole-Aitken Basin, for example.

Asmar, S.↗

NASA NTRS2004

Magnectic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynamical hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth s core and uncertainty. If this agrees with the seismologic value, then the hypotheses pass the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth s core, this yields a JGR-PI, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale change.In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at length- scale l/k, if one presumes kinetic energy is converted in either eddy- overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magneto- geostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core-mantle boundary. The resulting theoretical spectra for the core- source field and its SV are far more compatible with observation. The conversion time-scale of order 120 years is pseudo-scale-invariant. Magnetic spectra of other planets may differ; however, if a transition to non-conducting fluid hydrogen in Jupiter acts as barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte↗

NASA NTRS1982

Electrostatic disturbance forces on a 3-axis drag-free sensor

The electrostatic analysis of a multiple-capacitance 3-axis drag-free sensor is presented in this paper. The instrument consists of a proof-mass (a dense metallic ball) floating freely inside a spherical cavity enclosed by the sensor plates and the shield. Since the ball and the cavity are not necessarily concentric, the problem in three-dimensional potential theory for electrostatics is solved by the method of boundary perturbations and specifically in terms of spherical harmonics. The capacitance outputs of the instrument and the electrostatic forces acting on the system are derived as non-linear functions of the ball position, ball charge and the sensor plate potentials. The instrument sensitivity and cross-coupling effects are discussed. The analysis may also be useful for electrostatic gyros and suspensions.

Hagstrom, T.↗

NASA NTRS1995

A new method of presentation the large-scale magnetic field structure on the Sun and solar corona

The large-scale photospheric magnetic field, measured at Stanford, has been analyzed in terms of surface harmonics. Changes of the photospheric field which occur within whole solar rotation period can be resolved by this analysis. For this reason we used daily magnetograms of the line-of-sight magnetic field component observed from Earth over solar disc. We have estimated the period during which day-to-day full disc magnetograms must be collected. An original algorithm was applied to resolve time variations of spherical harmonics that reflect time evolution of large-scale magnetic field within solar rotation period. This method of magnetic field presentation can be useful enough in lack of direct magnetograph observations due to sometimes bad weather conditions. We have used the calculated surface harmonics to reconstruct the large-scale magnetic field structure on the source surface near the sun - the origin of heliospheric current sheet and solar wind streams. The obtained results have been compared with spacecraft in situ observations and geomagnetic activity. We tried to show that proposed technique can trace shon-time variations of heliospheric current sheet and short-lived solar wind streams. We have compared also our results with those obtained traditionally from potential field approximation and extrapolation using synoptic charts as initial boundary conditions.

Ponyavin, D. I.↗

NASA NTRS1994

Radiation from advanced solid rocket motor plumes

The overall objective of this study was to develop an understanding of solid rocket motor (SRM) plumes in sufficient detail to accurately explain the majority of plume radiation test data. Improved flowfield and radiation analysis codes were developed to accurately and efficiently account for all the factors which effect radiation heating from rocket plumes. These codes were verified by comparing predicted plume behavior with measured NASA/MSFC ASRM test data. Upon conducting a thorough review of the current state-of-the-art of SRM plume flowfield and radiation prediction methodology and the pertinent data base, the following analyses were developed for future design use. The NOZZRAD code was developed for preliminary base heating design and Al2O3 particle optical property data evaluation using a generalized two-flux solution to the radiative transfer equation. The IDARAD code was developed for rapid evaluation of plume radiation effects using the spherical harmonics method of differential approximation to the radiative transfer equation. The FDNS CFD code with fully coupled Euler-Lagrange particle tracking was validated by comparison to predictions made with the industry standard RAMP code for SRM nozzle flowfield analysis. The FDNS code provides the ability to analyze not only rocket nozzle flow, but also axisymmetric and three-dimensional plume flowfields with state-of-the-art CFD methodology. Procedures for conducting meaningful thermo-vision camera studies were developed.

Farmer, Richard C.↗

NASA NTRS1977

Onboard magnetic field modeling for Solar Maximum Mission /SMM/

Analysis and simulation results are presented for magnetic field models for use in attitude acquisition onboard Solar Maximum Mission (SMM). A study was made of the degree of the spherical harmonic expansion of the magnetic field required, considering mission requirements, modeling errors, and magnetometer quantization and biases. It is shown that a fifth-degree field is sufficient to provide two-degree roll angle determination accuracy with a residual magnetic bias of 10 milligauss. Also, a spherical harmonic expansion for the McIlwain L-parameter is included for the first time. This parameter will be telemetered to ground with experimental data. The fifth-degree expansion will provide the L-parameter to within two percent of accepted values. The additional onboard computational burden is the storage of 36 coefficients and an increase of about 15% in computation time. Prototype flight code was developed which is anticipated to require about 2000 bytes of core storage and 30 milliseconds of computation time per orbit point on the NSSC-1 computer.

Headrick, R. D.↗

NASA NTRS1981

Voyager 1 assessment of Jupiter's planetary magnetic field

An estimate of Jupiter's planetary magnetic field is obtained from the Voyager 1 observations of the Jovian magnetosphere. An explicit model for the magnetodisc current system is combined with a spherical harmonic model of the planetary field with both sets of parameters determined simultaneously using a nonlinear generalized inverse methodology. The resulting model fits the observations extremely well throughout the analysis interval (r 20 Jovian radii). The Jovian internal field model obtained from the Voyager 1 data is very similar to the octopole Pioneer 11 models. The best fitting magnetodisc lies in the centrifugal equator, 2/3 of the way between the rotational and magnetic equators, as appropriate for centrifugal loading of the magnetosphere by a cold plasma.

Connerney, J. E. P.↗

NASA NTRS1983

Inference of variations in the gravity field from satellite-to-satellite range rate

An analytic scheme for inferring variations of the gravity field from satellite-to-satellite range rate (low-low) is developed. As a test, it is applied to a pair of satellites in polar orbit, at altitude 160 km and spacing 100 km, with 72 data points per revolution. An assumed gravity field of tesseral spherical harmonics up to the eighth degree is completely recovered in three iterations over 64 revolutions. It is apparent that data points at regular intervals enable the use of data analysis techniques that avoid massive matrix inversions.

Kaula, W. M.↗

NASA NTRS1999

An Improved 360 Degree and Order Model of Venus Topography

We present an improved 360 degree and order spherical harmonic solution for Venus' topography. The new model uses the most recent set of Venus altimetry data with spacecraft positions derived from a recent high resolution gravity model. Geometric analysis indicates that the offset between the center of mass and center of figure of Venus is about 10 times smaller than that for the Earth, the Moon, or Mars. Statistical analyses confirm that the RMS topography follows a power law over the central part of the spectrum. Compared to the previous topography model, the new model is more highly correlated with Venus' harmonic gravity field.

Rappaport, Nicole J.↗

NASA NTRS2004

Magnetic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynamical hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth's core and uncertainty. If this agrees with the seismologic value, then the hypothesis passes the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth's core, this yields a generalized Stevenson-McLeod spectrum for the core-source field, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale range. In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at length-scale 1/k, if one presumes kinetic energy is converted in either eddy-overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small-scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magnetogeostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core-mantle boundary. The resulting theoretical spectra for the core-source field and its SV are far more compatible with observation. The conversion time-scale of order 120 years is pseudo-scale-invariant. Magnetic spectra of other planets may differ; however, if a transition to non-conducting fluid hydrogen in Jupiter acts as a barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte V.↗

NASA NTRS2004

Magnetic Probing of Core Geodynamics

To better understand geomagnetic theory and observation, we can use spatial magnetic spectra for the main field and secular variation to test core dynmcal hypotheses against seismology. The hypotheses lead to theoretical spectra which are fitted to observational spectra. Each fit yields an estimate of the radius of Earth's core and uncertainty. If this agrees with the seismologic value, then the hypothes pass the test. A new way to obtain theoretical spectra extends the hydromagnetic scale analysis of Benton to scale-variant field and flow. For narrow scale flow and a dynamically weak field by the top of Earth's core, this yields a generalized Stevenson-McLeod spectrum for the core-source field, and a secular variation spectrum modulated by a cubic polynomial in spherical harmonic degree n. The former passes the tests. The latter passes many tests, but does not describe rapid dipole decline and quadrupole rebound; some tests suggest it is a bit hard, or rich in narrow scale change. In a core geodynamo, motion of the fluid conductor does work against the Lorentz force. This converts kinetic into magnetic energy which, in turn, is lost to heat via Ohmic dissipation. In the analysis at lentgh-scale l/k, if one presumes kinetic energy is converted in either eddy- overturning or magnetic free-decay time-scales, then Kolmogorov or other spectra in conflict with observational spectra can result. Instead, the rate work is done roughly balances the dissipation rate, which is consistent with small scale flow. The conversion time-scale depends on dynamical constraints. These are summarized by the magneto-geostrophic vertical vorticity balance by the top of the core, which includes anisotropic effects of rotation, the magnetic field, and the core- mantle boundary. The resulting theoretical spectra for the core-source field and its SV are far more compatible with observation. The conversion time-scale of order l20 years is pseudo-scale-invarient. Magnetic spectra of other planets may differ; however, if a transition to non-conductmg fluid hydrogen in Jupiter acts as barrier to vertical flow, as well as current, then the shape of the jovi-magnetic spectrum could be remarkably Earth-like.

Voorhies, Coerte V.↗

NASA NTRS2010

High-Speed Solution of Spacecraft Trajectory Problems Using Taylor Series Integration

It has been known for some time that Taylor series (TS) integration is among the most efficient and accurate numerical methods in solving differential equations. However, the full benefit of the method has yet to be realized in calculating spacecraft trajectories, for two main reasons. First, most applications of Taylor series to trajectory propagation have focused on relatively simple problems of orbital motion or on specific problems and have not provided general applicability. Second, applications that have been more general have required use of a preprocessor, which inevitably imposes constraints on computational efficiency. The latter approach includes the work of Berryman et al., who solved the planetary n-body problem with relativistic effects. Their work specifically noted the computational inefficiencies arising from use of a preprocessor and pointed out the potential benefit of manually coding derivative routines. In this Engineering Note, we report on a systematic effort to directly implement Taylor series integration in an operational trajectory propagation code: the Spacecraft N-Body Analysis Program (SNAP). The present Taylor series implementation is unique in that it applies to spacecraft virtually anywhere in the solar system and can be used interchangeably with another integration method. SNAP is a high-fidelity trajectory propagator that includes force models for central body gravitation with N X N harmonics, other body gravitation with N X N harmonics, solar radiation pressure, atmospheric drag (for Earth orbits), and spacecraft thrusting (including shadowing). The governing equations are solved using an eighth-order Runge-Kutta Fehlberg (RKF) single-step method with variable step size control. In the present effort, TS is implemented by way of highly integrated subroutines that can be used interchangeably with RKF. This makes it possible to turn TS on or off during various phases of a mission. Current TS force models include central body gravitation with the J2 spherical harmonic, other body gravitation, thrust, constant atmospheric drag from Earth's atmosphere, and solar radiation pressure for a sphere under constant illumination. The purpose of this Engineering Note is to demonstrate the performance of TS integration in an operational trajectory analysis code and to compare it with a standard method, eighth-order RKF. Results show that TS is 16.6 times faster on average and is more accurate in 87.5% of the cases presented.

Scott, James R.↗

NASA NTRS1979

Development and analysis of a twelfth degree and order gravity model for Mars

Satellite geodesy techniques previously applied to artificial earth satellites have been extended to obtain a high-resolution gravity field for Mars. Two-way Doppler data collected by 10 Deep Space Network (DSN) stations during Mariner 9 and Viking 1 and 2 missions have been processed to obtain a twelfth degree and order spherical harmonic model for the martian gravitational potential. The quality of this model was evaluated by examining the rms residuals within the fit and the ability of the model to predict the spacecraft state beyond the fit. Both indicators show that more data and higher degree and order harmonics will be required to further refine our knowledge of the martian gravity field. The model presented shows much promise, since it resolves local gravity features which correlate highly with the martian topography. An isostatic analysis based on this model, as well as an error analysis, shows rather complete compensation on a global (long wavelength) scale. Though further model refinements are necessary to be certain, local (short wavelength) features such as the shield volcanos in Tharsis appear to be uncompensated. These are interpreted to place some bounds on the internal structure of Mars.

Christensen, E. J.↗