Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spherical harmonics”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1975

Winter bulge and diurnal variations in hydrogen inferred from AE-C composition measurements

The atomic hydrogen distribution at 250 km during December 1974 solstice was inferred, considering charge exchange equilibrium, from Atmosphere Explorer-C measurements of n(H(+)), n(O(+) and N(O). An empirical model, derived from the observations by least suqare analysis in terms of spherical harmonics, has the following characteristics: (1) n(H) increases by as much as a factor of two between the summer and winter hemispheres, (2) the n(H) diurnal variation is largest at the equator and (3) the diurnal variation is larger in the winter hemisphere than in the summer. Similar analysis of the gas temperature derived from n(N2) measurements reveals that all n(H) and Tg spherical harmonic coefficients are anticorrelated. Both the diurnal and latitudinal (annual) n(H) and Tg amplitudes are in substantial agreement with the zero flux condition, in which exospheric flow dominates the hydrogen distribution. The observed diurnal phase lag of n(H) with respect to Tg is about one hour, agreeing with theory.

Brinton, H. C.↗
NASA NTRS1985

The magnetic field at the core-mantle boundary

Models of the geomagnetic field are, in general, produced from a least-squares fit of the coefficients in a truncated spherical harmonic expansion to the available data. Downward continuation of such models to the core-mantle boundary (CMB) is an unstable process: the results are found to be critically dependent on the choice of truncation level. Modern techniques allow this fundamental difficulty to be circumvented. The method of stochastic inversion is applied to modeling the geomagnetic field. Prior information is introduced by requiring that the spectrum of spherical harmonic coefficients to fall-off in a particular manner which is consistent with the Ohmic heating in the core having a finite lower bound. This results in models with finite errors in the radial field at the CMB. Curves of zero radial field can then be determined and integrals of the radial field over patches on the CMB bounded by these null-flux curves calculated. With the assumption of negligible magnetic diffusion in the core; frozen-flux hypothesis, these integrals are time-invariant.

Bloxham, J.↗
NASA NTRS1992

Properties of the gravity fields of terrestrial planets

The properties of the gravity fields of the earth, Mars, and Venus, as expressed by spherical harmonic coefficients, are examined, using the harmonic expansions of the respective planetary topographies reported by Balmino et al. (1973), Bills and Ferrari (1978), and Bills and Kobrick (1985). The items examined include the spectral magnitudes and slopes of the gravity coefficients; the correlations between gravity and topography; and the correlations among different gravity harmonics, expressed by axiality and angularity. It was found that Venus differs from the other two planets in its great apparent depths of compensation, indicating a tectonics dominated by a stiff upper mantle. In addition, Venus has less activity deep in the mantle than do earth or Mars. Mars is marked by large gravity irregularities, as well as by their axial symmetry on a global scale. Although earth is probably the most peculiar planet, spherical harmonics do not bring out its varied characteristics. It is clearly a more active planet than Venus, with activity deep in the mantle. The lower magnitude of its higher harmonics is considered to be due to water recycled to the upper mantle.

Kaula, William M.↗
NASA NTRS1990

Determination of the ocean circulation using Geosat altimetry

A spherical harmonic model of the sea surface topography complete to degree and order 10 and a model of the earth's geopotential field complete to degree and order 50 have been obtained in a simultaneous solution using Geosat altimeter data and tracking data from 14 different satellites. The sea surface topography model compares well with oceanographic models computed using hydrographic data and ship drift data. Currently, errors in the estimated gravity field model limit the determination of the spherical harmonic coefficients of the general ocean circulation to degrees 10 and lower, corresponding to a minimum wavelength of 4000 km. Error analysis indicates that the correlation between the geoid and the sea surface topography model is less than 0.2, indicating good separation of the geoid and the sea surface topography at wavelengths of 4000 km or longer. Estimates of the scale factor for the significant wave height (H1/3), which is used to compute the electromagnetic bias correction and the bias for the Geosat altimeter, are obtained. The estimate of the H1/3 correction is 3.6 + or - 1.5 percent, and the height bias estimate is zero.

Nerem, R. S.↗
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 NTRS1993

The absorption of p-modes by sunspots - Variations with degree and order

A spherical harmonic decomposition of the p-modes into inward and outward propagating waves is employed to investigate the absorption of solar p-modes by an isolated sunspot. The absorption coefficient (averaged over frequency and azimuthal order) is found to increase with increasing horizontal wavenumber k over the range 0-0.8/Mm. For larger horizontal wavenumbers, in the range 0.8-1.5/Mm, the absorption coefficient decreases with increasing k. The absorption along each individual p-mode ridge tends to peak at an intermediate value of the spherical harmonic degree in the range 200-400. The highest absorption is found along the p(1) ridge, and the absorption decreases with increasing radial order.

Bogdan, Thomas J.↗
NASA NTRS2011

Weekly Solutions of Time-Variable Gravity from 1993 to 2010

The GRACE mission has been highly successful in determining the time-variable gravity field of the Earth, producing monthly or even more frequent solutions (cf. 10-day) solutions using both spherical harmonics and mascons. However the GRACE time series only commences in 2002 - 2003 and a gap of several years may occur in the series before a GRACE follow-on satellite is launched. Satellites tracked by SLR and DORIS have also been used to study time variations in the Earth's gravitational field. These include (most recently) the solutions of Cox and Chao (2002), Cheng et al. (2004, 2007) and Lemoine et al. (2007). In this paper we discuss the development of a new time series of low degree spherical harmonic fields based on the available SLR, DORIS and GPS data. We develop simultaneous solutions for both the geocenter and the low degree harmonics up to 5x5. The solutions integrate data from SLR geodetic satellites (e.g., Lageos1, Lageos2, Starlette, Stella, Ajisai, Larets, Westpac), altimetry satellites (TOPEX/Poseidon, Envisat, Jason-1, Jason-2), and satellites tracked solely by DORIS (e.g. SPOT2-5). We discuss some pertinent aspects of the satellite-specific modeling. We include altimeter crossovers in the weekly solutions where feasible and time permits. The resulting geocenter time series is compared with geophysical model predictions and other independently-derived solutions. Over the GRACE time period the fidelity and consistency with the GRACE solutions are presented.

Lemoine, F.↗
NASA NTRS1986

Transfer function analysis of thermospheric perturbations

Applying perturbation theory, a spectral model in terms of vectors spherical harmonics (Legendre polynomials) is used to describe the short term thermospheric perturbations originating in the auroral regions. The source may be Joule heating, particle precipitation or ExB ion drift-momentum coupling. A multiconstituent atmosphere is considered, allowing for the collisional momentum exchange between species including Ar, O2, N2, O, He and H. The coupled equations of energy, mass and momentum conservation are solved simultaneously for the major species N2 and O. Applying homogeneous boundary conditions, the integration is carred out from the Earth's surface up to 700 km. In the analysis, the spherical harmonics are treated as eigenfunctions, assuming that the Earth's rotation (and prevailing circulation) do not significantly affect perturbations with periods which are typically much less than one day. Under these simplifying assumptions, and given a particular source distribution in the vertical, a two dimensional transfer function is constructed to describe the three dimensional response of the atmosphere. In the order of increasing horizontal wave numbers (order of polynomials), this transfer function reveals five components. To compile the transfer function, the numerical computations are very time consuming (about 100 hours on a VAX for one particular vertical source distribution). However, given the transfer function, the atmospheric response in space and time (using Fourier integral representation) can be constructed with a few seconds of a central processing unit. This model is applied in a case study of wind and temperature measurements on the Dynamics Explorer B, which show features characteristic of a ringlike excitation source in the auroral oval. The data can be interpreted as gravity waves which are focused (and amplified) in the polar region and then are reflected to propagate toward lower latitudes.

Mayr, H. G.↗
DOE OSTI2021

Quantification of Deformation-Induced Concentration Fluctuations in Polymeric Liquids by Small-Angle Neutron Scattering

This study demonstrates how the large concentration fluctuations of polymeric liquids under deformation and flow can be quantitatively studied by applying the spherical harmonic expansion technique to small-angle neutron scattering. Using a series of binary polymer blends, it is shown that the emergence of the so-called butterfly patterns is caused by the change of sign in the leading anisotropic component of the small-angle spectrum, when the scattering is dominated by intermolecular correlation associated with flow-induced demixing. The increasing spatial fluctuations of concentration are evidenced by the enhancement of the isotropic component of the scattering spectrum in the zero-angle limit and peak shift of the leading anisotropic coefficients toward low Q. Additionally, the spherical harmonic expansion framework permits real-space analysis of anisotropic scattering length density correlation in a convenient form. The methodology described in this work provides a concrete venue for quantitative studies of phase behavior of fluids under deformation and flow via small-angle scattering techniques, where mainly qualitative approaches were previously employed.

Chemical structure↗
NASA NTRS1986

Statistics of the geomagnetic secular variation for the past 5Ma

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗
NASA NTRS1988

Statistics of the geomagnetic secular variation for the past 5 m.y

A new statistical model is proposed for the geomagnetic secular variation over the past 5Ma. Unlike previous models, the model makes use of statistical characteristics of the present day geomagnetic field. The spatial power spectrum of the non-dipole field is consistent with a white source near the core-mantle boundary with Gaussian distribution. After a suitable scaling, the spherical harmonic coefficients may be regarded as statistical samples from a single giant Gaussian process; this is the model of the non-dipole field. The model can be combined with an arbitrary statistical description of the dipole and probability density functions and cumulative distribution functions can be computed for declination and inclination that would be observed at any site on Earth's surface. Global paleomagnetic data spanning the past 5Ma are used to constrain the statistics of the dipole part of the field. A simple model is found to be consistent with the available data. An advantage of specifying the model in terms of the spherical harmonic coefficients is that it is a complete statistical description of the geomagnetic field, enabling us to test specific properties for a general description. Both intensity and directional data distributions may be tested to see if they satisfy the expected model distributions.

Constable, C. G.↗
NASA NTRS1993

A high resolution gravity model for Venus - GVM-1

A spherical harmonic model of the gravitational field of Venus complete to degree and order 50 has been developed using the S-band Doppler tracking data of the Pioneer Venus Orbiter (PVO) collected between 1979 and 1982. The short wavelengths of this model could only be resolved near the PVO periapse location (about 14 deg N latitude), therefore a priori constraints were applied to the model to bias poorly observed coefficients towards zero. The resulting model has a half-wavelength resolution of 400 km near the PVO periapse location, but the resolution degrades to greater than 1000 km near the poles. This gravity model correlates well with a degree 50 spherical harmonic expansion of the Venus topography derived from a combination of Magellan and PVO data. New tracking data from Magellan's gravity mission should provide some improvement to this model, although a complete model of the Venusian gravity field will depend on tracking of Magellan after the circularization of its orbit using aerobraking.

Nerem, R. S.↗
NASA NTRS2004

Unique Approaches to Analysis of Time-Variable Gravity from GRACE

We have developed an innovative analysis strategy for analysis of GRACE data. We have developed a capability to recover local/regional gravity changes using non-global functional representations (Le. surface anomalies vs. global spherical harmonics) h m the GRACE data. Our approach can take regularly or irregularly shaped regions, populate them with surface anomaly blocks of suitable area and solve for the resulting mass flux with respect to a mean field. The surface mass or gravity anomalies benefit from the application of spatial and temporal constraints to add stability to the solution. In this paper we discuss the analysis of four months of GRACE Level 1B data (accelerometry, intersatellite data, attitude information and precise orbits) from July to October 2003, recently released to the GRACE Science Team. We compare and contrast this local approach to gravity recovery, with the more conventional approach using global spherical harmonics. We review simulations of this technique which allow us to pinpoint optimum strategies for applications of this local gravity recovery approach.

Lemoine, Frank G.↗
NASA NTRS1998

Mesogranulation as A Distinct Scale of Convection in the Sun

We present evidence for the existence of mesogranulation as a scale of convection distinct from granulation and supergranulation through analysis of full-disk Doppler velocity images of the Sun collected by the Michelson Doppler Imager (MDI) aboard the NASA/ESA Solar and Heliospheric Observatory (SOHO). Our analysis procedures isolate nearly steady flows in the solar photosphere and yield power spectra of convection for spherical harmonic degrees up to I = 1000. Each spectrum exhibits an obvious supergranulation peak at I approximately 130 and a broad secondary peak at I approximately 600 with a distinct break in the spectrum between these peaks at I approximately 300. We believe that this secondary peak is a signature of mesogranulation with typical cell diameters of about 7 Mm. Our standard analysis procedure is to first remove the p-mode oscillation signal by averaging individual Dopplergrams over 17-minute intervals. Next, by fitting to standard functional forms we remove Doppler signals due to the motion of the spacecraft, the convective blueshift, solar rotation including differential rotation, and the meridional circulation in order to produce Dopplergrams dominated by convective motions. By mapping these processed images onto heliographic coordinates and projecting onto spherical harmonics, we produce a power spectrum of solar convection for each 17-minute period. We construct synthetic images and pass them through the same analysis procedure in order to determine the actual solar convection spectrum that reproduces the analyzed results. We find that a small but increasing percentage of high-degree convective power is lost in the analysis as we approach the limit of resolution of the detector but'that the broad, mesogranulation peak at I approximately 600 must be included in the convection spectrum of the synthetic images.

Bachmann, Kurt T.↗
NASA NTRS2010

GLGM-3: A Degree-ISO Lunar Gravity Model from the Historical Tracking Data of NASA Moon Orbiters

In preparation for the radio science experiment of the Lunar Reconnaissance Orbiter (LRO) mission, we analyzed the available radio tracking data of previous NASA lunar orbiters. Our goal was to use these historical observations in combination with the new low-altitude data to be obtained by LRO. We performed Precision Orbit Determination on trajectory arcs from Lunar Orbiter 1 in 1966 to Lunar Prospector in 1998, using the GEODYN II program developed at NASA Goddard Space Flight Center. We then created a set of normal equations and solved for the coefficients of a spherical harmonics expansion of the lunar gravity potential up to degree and order 150. The GLGM-3 solution obtained with a global Kaula constraint (2.5 x 10(exp -4)/sq l) shows good agreement with model LP150Q from the Jet Propulsion Laboratory, especially over the nearside. The levels of data fit with both gravity models are very similar (Doppler RMS of approx.0.2 and approx. 1-2 mm/s in the nominal and extended phases, respectiVely). Orbit overlaps and uncertainties estimated from the covariance matrix also agree well. GLGM-3 shows better correlation with lunar topography and admittance over the nearside at high degrees of expansion (l > 100), particularly near the poles. We also present three companion solutions, obtained with the same data set but using alternate inversion strategies that modify the power law constraint and expectation of the individual spherical harmonics coefficients. We give a detailed discussion of the performance of this family of gravity field solutions in terms of observation fit, orbit quality, and geophysical consistency.

Mazarico, E.↗
NASA NTRS1993

Nearly steady flows in GONG prototype data

Doppler velocity images obtained with the GONG prototype instrument were analyzed to measure the nearly steady photospheric flows. The data consists of 88 images each of velocity, intensity, and modulation obtained at 20:00 UT on 88 days from July 1992 to February 1994. Each velocity image was temporally filtered to remove the p-mode oscillations, masked to exclude active regions, and then analyzed using spherical harmonics and orthogonal functions as described by Hathaway (1992). The spectral coefficients show very consistent results for the entire time interval with some evidence of year-to-year variations. The rotation profile agrees well with previous results and exhibits a north-south asymmetry that reverses sign during the 20 month interval. The residual rotation velocities exhibit structures with amplitudes of approximately 5 m/s that may be related to torsional oscillations. The meridional circulation is directed from the equator toward the poles with a peak velocity in the photosphere of approximately 50 m/s. The higher order components are very weak but indicate a divergent flow from the mid-latitudes (opposite that found for the June 1989 data). The convective limb shift is well fit by a 3rd order polynomial. The convection spectrum has a prominent peak at spherical harmonic degrees of l approximately 150 with very little signal in the low degree modes. Analysis of this signal shows that there is no evidence for giant cell convection at the level of approximately 10 m/s for all modes up to l = 32.

Hathaway, David H.↗
NASA NTRS1980

Global gravity field to degree and order 30 from Geos 3 satellite altimetry and other data

A model of the geopotential field in spherical harmonics to degree and order 30 is obtained from Geos 3 satellite to sea surface altimetry data, terrestrial gravity measurements and satellite perturbation analysis. A general perturbation solution is employed for the calculation of the orbits of 10 satellites based on satellite laser ranging data, and 1 deg x 1 deg surface gravity data are used to compute 550 km x 550 km block anomalies by means of autocovariance analysis. Altimeter-determined sea-surface heights, which are taken as the geoid, are averaged for each 1 deg x 1 deg ocean surface area and treated by autocovariance analysis to obtain 550 x 550 km block undulations. Observation and normal equations are formed from the altimeter and surface gravity data, which together cover 1635 out of 1654 possible surface elements, and are combined with the available satellite-derived normal equations to obtain a solution for the spherical harmonics coefficients. In addition, a value of 6,378,138.23 + or - 1.3 m is obtained for the earth's semimajor axis.

Gaposchkin, E. M.↗