SEARCH · Engineering Papers
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.
Topics addressed include: global gravity model development; methods for approximation of the gravity field; gravity field measuring techniques; global gravity field applications and requirements in geophysics and oceanography; and future gravity missions.
A steady WKB model of gravity wave propagation including convective adjustment is used to investigate approximations used in various gravity-wave parameterization schemes. First, it is shown that estimates of the wave breaking height assuming a single horizontal wavenumber gravity wave can lead to errors if the topography is not sinusoidal. Second, the model results show that the assumption that wave growth ceases with the onset of convection or shear instability is an oversimplification. Since convection appears first over a very limited spatial region of the wave field, the wave is initially unaffected by turbulent mixing. However, when the convection zone spreads over a large portion of the wave field the amplitude is constrained. Estimates of the heat flux by breaking gravity waves are used to develop a simple parameterization of the vertical diffusion in terms of the Reynolds stress.
A nearby analytic model of a gravity wave driven middle atmosphere circulation is formulated. The simplified model represents the one-dimensional interaction of a single gravity wave mode with a zonal-mean flow. Solutions of this model are shown to agree very well with numerical results from a more complete mean-flow/gravity-wave model (essentially the model developed by Holton, 1982). A single nondimensional parameter, the ratio of a mean flow radiative forcing parameter, and the gravity wave momentum flux incident into the middle atmosphere from below largely determine the solution of the model. For typical middle atmosphere parameter values, an approximate solution can be obtained, and this solution permits the parameter dependence of the circulation to be characterized very simply. The gravity wave driven circulation extends downward from the level (the breaking level, approximately) where the momentum flux equals the mean flow radiative forcing.
The GEOS-C spacecraft is scheduled to carry onboard a radar altimeter for the purpose of measuring the geoid undulations in oceanic areas. An independently derived geoid map will provide a valuable complement to these experiments. A detailed gravimetric geoid is presented for the Atlantic and northeast Pacific Ocean areas based upon a combination of the Goddard Space Flight Center GEM-6 earth model and surface 1 deg x 1 deg gravity data. As part of this work a number of satellite derived gravity models were evaluated to establish the model which best represented the long wave length features of the geoid in the above mentioned area. Comparisons of the detailed geoid with the astrogeodetic data provided by the National Ocean Survey and dynamically derived tracking station heights indicate that the accuracy of this combined geoid is on the order of 2 meters or better where data was dense and 5 to 7 meters where data was less dense.
The ARTEMIS mission, part of the THEMIS extended mission, is the first to fly spacecraft in the Earth-Moon Lissajous regions. In 2009, two of the five THEMIS spacecraft were redeployed from Earth-centered orbits to arrive in Earth-Moon Lissajous orbits in late 2010. Starting in August 2010, the ARTEMIS P1 spacecraft executed numerous stationkeeping maneuvers, initially maintaining a lunar L2 Lissajous orbit before transitioning into a lunar L1 orbit. The ARTEMIS P2 spacecraft entered a L1 Lissajous orbit in October 2010. In April 2011, both ARTEMIS spacecraft will suspend Lissajous stationkeeping and will be maneuvered into lunar orbits. The success of the ARTEMIS mission has allowed the science team to gather unprecedented magnetospheric measurements in the lunar Lissajous regions. In order to effectively perform lunar Lissajous stationkeeping maneuvers, the ARTEMIS operations team has provided orbit determination solutions with typical accuracies on the order of 0.1 km in position and 0.1 cm/s in velocity. The ARTEMIS team utilizes the Goddard Trajectory Determination System (GTDS), using a batch least squares method, to process range and Doppler tracking measurements from the NASA Deep Space Network (DSN), Berkeley Ground Station (BGS), Merritt Island (MILA) station, and United Space Network (USN). The team has also investigated processing of the same tracking data measurements using the Orbit Determination Tool Kit (ODTK) software, which uses an extended Kalman filter and recursive smoother to estimate the orbit. The orbit determination results from each of these methods will be presented and we will discuss the advantages and disadvantages associated with using each method in the lunar Lissajous regions. Orbit determination accuracy is dependent on both the quality and quantity of tracking measurements, fidelity of the orbit force models, and the estimation techniques used. Prior to Lissajous operations, the team determined the appropriate quantity of tracking measurements that would be needed to meet the required orbit determination accuracies. Analysts used the Orbit Determination Error Analysis System (ODEAS) to perform covariance analyses using various tracking data schedules. From this analysis, it was determined that 3.5 hours of DSN TRK-2-34 range and Doppler tracking data every other day would suffice to meet the predictive orbit knowledge accuracies in the Lissajous region. The results of this analysis are presented. Both GTDS and ODTK have high-fidelity environmental orbit force models that allow for very accurate orbit estimation in the lunar Lissajous regime. These models include solar radiation pressure, Earth and Moon gravity models, third body gravitational effects from the Sun, and to a lesser extent third body gravitational effects from Jupiter, Venus, Saturn, and Mars. Increased position and velocity uncertainties following each maneuver, due to small execution performance errors, requires that several days of post-maneuver tracking data be processed to converge on an accurate post-maneuver orbit solution. The effects of maneuvers on orbit determination accuracy will be presented, including a comparison of the batch least squares technique to the extended Kalman filter/smoother technique. We will present the maneuver calibration results derived from processing post-maneuver tracking data. A dominant error in the orbit estimation process is the uncertainty in solar radiation pressure and the resultant force on the spacecraft. An estimation of this value can include many related factors, such as the uncertainty in spacecraft reflectivity and surface area which is a function of spacecraft orientation (spin-axis attitude), uncertainty in spacecraft wet mass, and potential seasonal variability due to the changing direction of the Sun line relative to the Earth-Moon Lissajous reference frame. In addition, each spacecraft occasionally enters into Earth or Moon penumbra or umbra and these shadow crossings reduche solar radiation force for several hours. The effects of these events on orbit determination accuracy will be presented. In order to plan for upcoming stationkeeping maneuvers, the maneuver planning team must take the current orbit estimate, propagate it forward to the planned maneuver time, and determine the optimal maneuver to maintain the Lissajous orbit for one or more revolutions. The propagation is performed using a Runge-Kutta 7/8 integrator and typically the position and velocity uncertainty increases with propagation time, increasing the overall uncertainty of the orbit state at the maneuver execution time. The effect of orbit knowledge uncertainty on stationkeeping operations will be presented.
The development of appropriate modeling and adjustment procedures for the estimation of harmonic coefficients of the geopotential, from surface gravity data was studied, in order to provide an optimum way of utilizing the terrestrial gravity information in combination solutions currently developed at NASA/Goddard Space Flight Center, for use in the TOPEX/POSEIDON mission. The mathematical modeling was based on the fundamental boundary condition of the linearized Molodensky boundary value problem. Atmospheric and ellipsoidal corrections were applied to the surface anomalies. Terrestrial gravity solutions were found to be in good agreement with the satellite ones over areas which are well surveyed (gravimetrically), such as North America or Australia. However, systematic differences between the terrestrial only models and GEMT1, over extended regions in Africa, the Soviet Union, and China were found. In Africa, gravity anomaly differences on the order of 20 mgals and undulation differences on the order of 15 meters, over regions extending 2000 km in diameter, occur. Comparisons of the GEMT1 implied undulations with 32 well distributed Doppler derived undulations gave an RMS difference of 2.6 m, while corresponding comparison with undulations implied by the terrestrial solution gave RMS difference on the order of 15 m, which implies that the terrestrial data in that region are substantially in error.
TOPEX/POSEIDON sea surface height measurements are examined for quantitative consistency with known elements of the oceanic general circulation and its variability. Project-provided corrections were accepted but are at tested as part of the overall results. The ocean was treated as static over each 10-day repeat cycle and maps constructed of the absolute sea surface topography from simple averages in 2 deg x 2 deg bins. A hybrid geoid model formed from a combination of the recent Joint Gravity Model-2 and the project-provided Ohio State University geoid was used to estimate the absolute topography in each 10-day period. Results are examined in terms of the annual average, seasonal average, seasonal variations, and variations near the repeat period. Conclusion are as follows: the orbit error is now difficult to observe, having been reduced to a level at or below the level of other error sources; the geoid dominates the error budget of the estimates of the absolute topography; the estimated seasonal cycle is consistent with prior estimates; shorter-period variability is dominated on the largest scales by an oscillation near 50 days in spherical harmonics Y(sup m)(sub 1)(theta, lambda) with an amplitude near 10 cm, close to the simplest alias of the M(sub 2) tide. This spectral peak and others visible in the periodograms support the hypothesis that the largest remaining time-dependent errors lie in the tidal models. Though discrepancies attribute to the geoid are within the formal uncertainties of the good estimates, removal of them is urgent for circulation studies. Current gross accuracy of the TOPEX/POSEIDON mission is in the range of 5-10 cm, distributed overbroad band of frequencies and wavenumbers. In finite bands, accuracies approach the 1-cm level, and expected improvements arising from extended mission duration should reduce these numbers by nearly an order of magnitude.
The NASA GSFC DORIS analysis center has processed data from January 1993 to December 2014 and provided 1141 weekly solutionsin the form of normal equations for incorporation into the DORIS solution for ITRF2014. The solution time series, designated asgscwd26, were based on tracking data to eleven DORIS satellites divided generally into seven-day arcs. With respect to the ITRF2008 submission (Le Bail et al., 2010), the measurement model was updated to model the beacon frequency variations at certain DORIS sites,to apply the DORIS antenna phase law for the Starec and Alcatel antennae, and to apply the antenna offset corrections in the NASAGSFC orbit determination software rather than using the data-supplied corrections. We show that computing the antenna offset correctionsin the orbit determination software is superior to using the offset corrections that are supplied with the DORIS data, and that thisimproves the RMS of fit for SPOT-2, Envisat, SPOT-4, and SPOT-5. The updates for the force model included: (1) the development of improved nonconservative force modeling for SPOT-2, SPOT-3, SPOT-5, Envisat, and HY-2A, and (2) the application of an updatedstatic gravity model based on GRACE and GOCE data, and weekly models of the variations in the low degree gravity field deducedindependently from tracking by Satellite Laser Ranging (SLR) and DORIS. The post-ITRF2008 DORIS coordinate WRMS afterthe launch of Envisat and SPOT-5 is improved from 11.20 to 12.45 mm with ITRF2008 (Le Bail et al., 2010), to between 8.50 and9.99 mm with the gscwd26 SINEX solution. The application of the DORIS antenna phase laws shifts the DORIS scale wrt DPOD2008by +6.0 mm from 19930103 to 20020606, and by +11.4 mm from 20020613 to 20111030. The application of more detailed models of time-variable gravity reduces the slopes in the Helmert transformation parameters Tx, and Ty (w.r.t. DPOD2008) after 2005. The annual amplitude in these parameters is reduced from 3.2 mm (for Tx), 4.1 mm (for Ty), to 1.7 mm (for Tx) and 2.8 mm (for Ty).
Measurement of the gravity field of the Moon has a high science priority because of its implications for the internal structure and thermal history of the Moon, and it has a high priority for future exploration activities because of the influence of lunar gravity on spacecraft navigation and orbit maintenance. The current state of knowledge in the lunar gravity field (and the uncertainty in the knowledge) is based primarily on data accumulated from the Lunar Orbiter and Apollo programs. Data are sparse and emphasize the equatorial band (+/- 30 deg) on the near side of the Moon. There are no tracking data on the far side and only the Lunar Orbiter 5 provides a small amount of high inclination data. A host of gravity models developed from different combinations of tracking data have large discrepancies in their predictions of spacecraft motion and orbit lifetimes. There are also large disagreements in the Mercator projections of the gravity acceleration from each model, especially on the far side, where the contours tend to have no obvious relationship with the local topography. The science and engineering requirements for global gravity field mapping will be satisfied with continuous radio metric tracking of Lunar Scout 1 in a low polar orbit using the Deep Space Network and Lunar Scout 2 in a high elliptical orbit. The gravity field of the Moon will be mapped during the Scout Program using a two spacecraft concept. In the two spacecraft concept, one spacecraft is placed in a high altitude eccentric orbit while the second spacecraft is in a low altitude polar orbit. The gravity experiment requires a radio frequency that will permit two-way Doppler tracking between the spacecraft and the Deep Space Network (DSN). Both spacecraft carry NASA standard transponder systems for data transmission to Earth as well as for tracking and orbit determination. Data sufficient to produce a gravity field map could be acquired within one month with this system.
Thermal balance analysis of temperature distributions in cylindrical tube considered to model gravity gradient attitude control rod isolated in space
This paper examines the accuracy of recent Mars gravity models and the importance of perturbations due to the Mars radiation pressure and the Martian moons, Phobos and Deimos, on the trajectories of Mars orbiters. A linear orbit perturbation theory is used to characterize the patterns of gravity field near resonances for the Viking and Mariner 9 spacecraft. These resonances are shown to have considerable power and their potential for contributing to Mars gravity solutions is emphasized. It is shown that some of the same resonance orders which appear in the Viking orbits, dominate the radial orbit error spectrum for Mars Observer. Results of orbit determination simulations at the Goddard Space Flight Center show that the perturbations caused by the Martian moons and the Mars radiation pressure are larger than 0.1 mm/s, the expected precision of the Mars Observer Doppler tracking data. Tests with the Viking Doppler data indicate that best analysis of these data mandates the inclusion of the Phobos gravitational perturbation in the modeling of Viking spacecraft trajectories.
An improved model for the Earth's gravity field, TEG-1, was determined using data sets from fourteen satellites, spanning the inclination ranges from 15 to 115 deg, and global surface gravity anomaly data. The satellite measurements include laser ranging data, Doppler range-rate data, and satellite-to-ocean radar altimeter data measurements, which include the direct height measurement and the differenced measurements at ground track crossings (crossover measurements). Also determined was another gravity field model, TEG-1S, which included all the data sets in TEG-1 with the exception of direct altimeter data. The effort has included an intense scrutiny of the gravity field solution methodology. The estimated parameters included geopotential coefficients complete to degree and order 50 with selected higher order coefficients, ocean and solid Earth tide parameters, Doppler tracking station coordinates and the quasi-stationary sea surface topography. Extensive error analysis and calibration of the formal covariance matrix indicate that the gravity field model is a significant improvement over previous models and can be used for general applications in geodesy.
Pioneer Venus orbiter (PVO) gravity profiles over Aphrodite Terra are compared with profiles predicted from models of thermal isostasy, mantle convection, and Airy compensation. Topography profiles across Aphrodite are compared to model spreading-ridge profiles in order to further assess this model. Airy compensation depths and convection-layer thicknesses are greater under eastern Aphrodite than western Aphrodite. Compensation depths in the east are greater than most estimates of lithospheric thickness, suggesting that this part of the ridge is dynamically supported. In parts of western Aphrodite, the spreading-ridge model gravity provides a better fit to the data than either Airy compensation or mantle convection. Best-fit spreading rates are between 0.3 and 1.6 cm/yr. Airy compensation and mantle convection cannot be distinguished in most places using only PVO data.
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.
The current known geodetic parameters of the planets and their moons are reviewed. A 3:2 spin-orbit resonance has been calculated for Mercury, a planet for which previous mass estimates are suggested to be inaccurate by 30 pct. The Pioneer Venus Orbiter data indicated a Venus diameter of 6051 km, a geoid surface highly-correlated with the surface topography, and regional gravity anomalies. Reflective mirrors on the surface and data from the Lunar Orbiter 4 have produced values for the principal polar moment of inertia homogeneity factor, the earth-moon mass ratio, a fifth degree and order spherical harmonic gravity model, the lunar acceleration, and the lunar potential Love number. The Mars gravity field is now estimated up to 12th degree and order spherical harmonics from Mariner 9 and Viking Orbiters 1 and 2, which also obtained new topographical data. New gravity fields and mass are being calculated for the Jovian and Saturnian moons. The outer planets' radii and rotational periods have recently been revised to greater accuracy. New data will be forthcoming from the Space Telescope and the planned Venus Radar Mapper.
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.
Follow-on missions to provide continuity in the observation of the sea surface topography once the successful TOPEX/POSEIDON (T/P) oceanographic satellite mission has ended are discussed. Candidates include orbits which follow the ground tracks of T/P GEOSAT or ERS-1. The T/P precision ephemerides, estimated to be near 3 cm root-mean-square, demonstrate the radial orbit accuracy that can be achieved at 1300 km altitude. However, the radial orbit accuracy which can be achieved for a mission at the 800 km altitudes of GEOSAT and ERS-1 has not been established, and achieving an accuracy commensurate with T/P will pose a great challenge. This investigation focuses on the radial orbit accuracy that can be achieved for a mission in the GEOSAT orbit. Emphasis is given to characterizing the effects of force model errors on the estimated radial orbit accuracy, particularly those due to gravity and drag. The importance of global, continuous tracking of the satellite for reduction in these sources of orbit error is demonstrated with simulated GPS tracking data. A gravity tuning experiment is carried out to show how the effects of gravity error may be reduced. Assuming a GPS flight receiver with a full-sky tracking capability, the simulation results indicate that a 5 cm radial orbit accuracy for an altimeter satellite in GEOSAT orbit should be achievable during low-drag atmospheric conditions and after an acceptable tuning of the gravity model.
A technique is described for improving the positioning accuracy of earth survey satellites. Tracking data from Landsat-1, taken in contiguous spans equal in length to the repeat period of 18 days, have been reduced in a least-squares orbit-determination scheme to estimate two pairs of lumped resonance coefficients. Seven separate data spans yielded resonance coefficient pairs consistent to better than 10% of their values, based on the GEM-1 gravity model. Constraint equations were developed. The resonance terms were embedded in the GEM-8 gravity field to yield further improvements in orbital fit to the tracking data, and tests suggest a 20 m position capability for Landsat orbits when the resonance model is adopted. The lumped terms may be used to improve the perturbation model for all Landsat orbits, and the approach is applicable to the orbit of any satellite whose repeat period is long enough for geopotential resonance to have a significant effect on its motion.