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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2016

Normalization and Implementation of Three Gravitational Acceleration Models

Unlike the uniform density spherical shell approximations of Newton, the consequence of spaceflight in the real universe is that gravitational fields are sensitive to the asphericity of their generating central bodies. The gravitational potential of an aspherical central body is typically resolved using spherical harmonic approximations. However, attempting to directly calculate the spherical harmonic approximations results in at least two singularities that must be removed to generalize the method and solve for any possible orbit, including polar orbits. Samuel Pines, Bill Lear, and Robert Gottlieb developed three unique algorithms to eliminate these singularities. This paper documents the methodical normalization of two of the three known formulations for singularity-free gravitational acceleration (namely, the Lear and Gottlieb algorithms) and formulates a general method for defining normalization parameters used to generate normalized Legendre polynomials and Associated Legendre Functions (ALFs) for any algorithm. A treatment of the conventional formulation of the gravitational potential and acceleration is also provided, in addition to a brief overview of the philosophical differences between the three known singularity-free algorithms.

Eckman, Randy A.↗
NASA NTRS1975

Evaluation of geopotential and luni-solar perturbations by a recursive algorithm

The disturbing functions due to the geopotential and Luni-solar attractions are linear and bilinear forms in spherical harmonics. Making use of recurrence relations for the solid spherical harmonics and their derivatives, recurrence formulas are obtained for high degree terms as function of lower degree for any term of those disturbing functions and their derivative with respect to any element. The equations obtained are effective when a numerical integration of the equations of motion is appropriate. In analytical theories, they provide a fast way of obtaining high degree terms starting from initial very simple functions.

Giacaglia, G. E. O.↗
NASA NTRS1993

Venus gravity and topography: 60th degree and order model

We have combined the most recent Pioneer Venus Orbiter (PVO) and Magellan (MGN) data with the earlier 1978-1982 PVO data set to obtain a new 60th degree and order spherical harmonic gravity model and a 120th degree and order spherical harmonic topography model. Free-air gravity maps are shown over regions where the most marked improvement has been obtained (Ishtar-Terra, Alpha, Bell and Artemis). Gravity versus topography relationships are presented as correlations per degree and axes orientation.

Konopliv, A. S.↗
NASA NTRS1999

Backus Effect on a Perpendicular Errors in Harmonic Models of Real vs. Synthetic Data

Measurements of geomagnetic scalar intensity on a thin spherical shell alone are not enough to separate internal from external source fields; moreover, such scalar data are not enough for accurate modeling of the vector field from internal sources because of unmodeled fields and small data errors. Spherical harmonic models of the geomagnetic potential fitted to scalar data alone therefore suffer from well-understood Backus effect and perpendicular errors. Curiously, errors in some models of simulated 'data' are very much less than those in models of real data. We analyze select Magsat vector and scalar measurements separately to illustrate Backus effect and perpendicular errors in models of real scalar data. By using a model to synthesize 'data' at the observation points, and by adding various types of 'noise', we illustrate such errors in models of synthetic 'data'. Perpendicular errors prove quite sensitive to the maximum degree in the spherical harmonic expansion of the potential field model fitted to the scalar data. Small errors in models of synthetic 'data' are found to be an artifact of matched truncation levels. For example, consider scalar synthetic 'data' computed from a degree 14 model. A degree 14 model fitted to such synthetic 'data' yields negligible error, but amplifies 4 nT (rmss) added noise into a 60 nT error (rmss); however, a degree 12 model fitted to the noisy 'data' suffers a 492 nT error (rmms through degree 12). Geomagnetic measurements remain unaware of model truncation, so the small errors indicated by some simulations cannot be realized in practice. Errors in models fitted to scalar data alone approach 1000 nT (rmss) and several thousand nT (maximum).

Voorhies, C. V.↗
DOE OSTI2023

Quantifying molecular deformation in polymer melts by a generalized Zimm plot approach

The Zimm plot has been widely used to characterize the molecular dimensions of polymers from small-angle scattering experiments, where the reciprocal intensity is analyzed as a function of the square of the magnitude of the scattering wavevector Q. Here, this work explores the benefits of analyzing the reciprocal scattering intensity from deformed polymers, extending the original Zimm plot to anisotropic materials. In the small-angle limit, a tensorial extension of the Guinier law is found for the gyration tensor and the reciprocal single-chain structure factor. In the high-Q limit, application of the spherical harmonic expansion technique to the reciprocal structure factor permits direct model-independent analysis of spatially dependent molecular deformation of polymers. Additionally, the contributions from high-order spherical harmonics become insignificant in the reciprocal-intensity representation. The proposed generalized Zimm plot approach is demonstrated computationally with the affine deformation model and the Rouse model, and experimentally with small-angle neutron scattering measurements of deformed polystyrene melts.

36 MATERIALS SCIENCE↗
NASA NTRS1992

Simulations of lunar gravity field determination for Lunar Observer

The current plan for the Lunar Observer (LO) mission is to launch in the late 1990s and insert LO into a 100 km polar circular mapping orbit. However, prior to the mapping orbit, LO will be placed in a gravity calibration orbit (GCO) at a higher altitude to determine the gravity field of the moon. This paper examines the abilities of two GCO orbits (at 200 and 500 km altitudes) to recover a high degree and order gravity truth model that includes spherical harmonics and mascons by estimating different degree and order gravity fields with spherical harmonics only. This is achieved by comparing radial accelerations from the true and estimated models at the mapping altitude and by comparing trajectory propagations. For the gravity fields estimated (up to 30th degree and order), the 500 km GCO was just as successful as the 200 km GCO in determining the gravity field.

Konopliv, Alex↗
NASA NTRS1982

Equivalent source modeling of the main field using MAGSAT data

An iterative least squares estimation algorithm with the capability for including a priori statistical information was implemented to recover multiple magnetic dipole models of the Earth's main magnetic field. The dipoles are fixed to a specified radius at or below the core-mantle boundary and centered on equal area blocks. The algorithm can solve for dipole magnitudes only (fixed orientations), or allow full freedom of orientation and solve for vector components. External field parameters and observatory anomaly biases can also be estimated simultaneously. Time dependence is modeled using first time derivatives for dipole vector components. Single-epoch and time dependent dipole models are derived using MAGSAT and observatory annul means data. Equivalent spherical harmonic representation are computed in closed form from the dipole models and compared with truncated spherical harmonic models estimated in the standard way from the same data sets. In particular, a 21 deg spatial resolution model based on 93 dipoles was computed based on observatory annual means data and a selected MAGSAT data set and was compared with candidate IGRF 1975 models and their 1980 secular variation.

Mayhew, M. A.↗
NASA NTRS1994

A 360-degree and -order model of Venus topography

This report presents the most recent spherical harmonic topography model of Venus developed at Jet Propulsion Laboratory. It was produced by a spherical harmonic analysis of the most complete set of Magellan altimetry data, augmented by Pioneer Venus and Venera data. The harmonic coefficients of the topography were computed to degree and order 360. Compared to previous topography models, this one has the highest correlation with the gravity field of Venus.

Rappaport, Nicole↗
DOE OSTI2022

Global Barotropic Tide Modeling Using Inline Self‐Attraction and Loading in MPAS‐Ocean

Abstract We examine ocean tides in the barotropic version of the Model for Prediction Across Scales (MPAS‐Ocean), the ocean component of the Department of Energy Earth system model. We focus on four factors that affect tidal accuracy: self‐attraction and loading (SAL), model resolution, details of the underlying bathymetry, and parameterized topographic wave drag. The SAL term accounts for the tidal loading of Earth's crust and the self‐gravitation of the ocean and the load‐deformed Earth. A common method for calculating SAL is to decompose mass anomalies into their spherical harmonic constituents. Here, we compare a scalar SAL approximation versus an inline SAL using a fast spherical harmonic transform package. Wave drag accounts for energy lost by breaking internal tides that are produced by barotropic tidal flow over topographic features. We compare a series of successively finer quasi‐uniform resolution meshes (62.9, 31.5, 15.7, and 7.87 km) to a variable resolution (45 to 5 km) configuration. We ran MPAS‐Ocean in a single‐layer barotropic mode forced by five tidal constituents. The 45 to 5 km variable resolution mesh obtained the best total root‐mean‐square error (5.4 cm) for the deep ocean ( 1,000 m) tide compared to TPXO8 and ran twice as fast as the quasi‐uniform 8 km mesh, which had an error of 5.8 cm. This error is comparable to those found in other forward (non‐assimilative) ocean tide models. In future work, we plan to use MPAS‐Ocean to study tidal interactions with other Earth system components, and the tidal response to climate change.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2023

Constraining the Structure under Lunar Impact Basins with Gravity

The lunar gravity field is used to estimate and constrain the depth of mass anomalies under 19 major lunar impact basins. We use radial gravitational spectra, consisting of accelerations computed either per spherical harmonic degree or cumulatively, at surface locations to obtain the distribution of the gravity signal with spherical harmonic degree and, by implication, to the likely depth below the surface. The results provide estimates for the maximum likely depths of the primary component to the mass anomalies under 19 basins. We find that the maximum depths of the primary source of mascon gravity on the lunar nearside are deeper than the depths for those on the farside when South Pole–Aitken (SPA) is excluded. All basin mass anomalies on the lunar nearside are in the mantle. The maximum depth of the primary source of the mass anomalies is 200 km beneath the surface. The upper 20 km under all basins is largely devoid of anomalies, reflecting predominantly mixing and relaxation associated with impact melt combined with ejecta fallback, as well as homogenization associated with post-basin formation impact bombardment. Except for SPA, all basin anomalies merge with the deep interior at ∼150 km or below, indicating the depth penetration of disruption of the density structure of the lunar interior associated with impact bombardment.

Lunar gravitational field↗
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 NTRS2010

VLBI2020: From Reality to Vision

The individual apparent motions of distant radio sources are believed to be caused by the effect of intrinsic structure variations of the active galactic nuclei (AGN). However, some cosmological models of the expanded Universe predict that systematic astrometric proper motions of distant quasars do not vanish as the radial distance from the observer to the quasar grows. These systematic effects can even increase with the distance, making it possible to measure them with high-precision astrometric techniques like VLBI. The Galactocentric acceleration of the Solar System barycenter may cause a secular aberration drift with a magnitude of 4 uas/yr. The Solar System motion relative to the cosmic microwave background produces an additional dipole effect, proportional to red shift. We analyzed geodetic VLBI data spanning from 1979 until 2009 to estimate the vector spherical harmonics in the expansion of the vector field of the proper motion of 687 radio sources. The dipole and quadrupole vector spherical harmonics were estimated with an accuracy of 1-5 as/yr. We have shown that over the next decade the geodetic VLBI may approach the level of accuracy on which the cosmological models of the Universe could be tested. Hence, it is important to organize a dedicated observational program to increase the number of measured proper motions to 3000.

Titov, Oleg↗
NASA NTRS1989

Downward continuation of the free-air gravity anomalies to the ellipsoid using the gradient solution and terrain correction: An attempt of global numerical computations

The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies have to be analytically downward continued from the earth's surface to a sphere-at least to the ellipsoid. The goal is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models. The basic method for the downward continuation is the gradient solution (the g sub 1 term). The terrain correction was also computed because of the role it can play as a correction term when calculating harmonic coefficients from surface gravity data. The fast Fourier transformation was applied to the computations.

Wang, Y. M.↗
NASA NTRS1992

Excitation of pressure modes in common envelopes

The excitation of oscillatory modes by a low-mass star orbiting inside a common envelope with a more massive star is investigated, with emphasis on adiabatic high spherical harmonic degree (l is much greater than 1) modes propagating outward in the envelope. The dominant oscillatory modes are those for which the spherical harmonic order is large, and thus the amplitudes are large close to the equatorial plane and small closer to the poles. A secondary of mass about 1 percent of the primary mass excites modes with relative surface amplitudes of a few x 10 percent. Even a Jupiter-like brown dwarf, when it is very deep in the envelope of an asymptotic giant branch star, can cause perturbations of relative surface amplitudes of about 10 percent near the equatorial plane. If mass loss is influenced by oscillation, then this mechanism can lead to higher mass loss in the equatorial plane. In case of an asymptotic giant branch primary, this might lead to the formation of an elliptical planetary nebula.

Soker, Noam↗
NASA NTRS1982

The latitude dependence of the variance of zonally averaged quantities

Geometric characteristics of the spherical earth are shown to be responsible for the increase of variance with latitude of zonally averaged meteorological statistics. An analytic model is constructed to display the effect of a spherical geometry on zonal averages, employing a sphere labeled with radial unit vectors in a real, stochastic field expanded in complex spherical harmonics. The variance of a zonally averaged field is found to be expressible in terms of the spectrum of the vector field of the spherical harmonics. A maximum variance is then located at the poles, and the ratio of the variance to the zonally averaged grid-point variance, weighted by the cosine of the latitude, yields the zonal correlation typical of the latitude. An example is provided for the 500 mb level in the Northern Hemisphere compared to 15 years of data. Variance is determined to increase north of 60 deg latitude.

North, G. R.↗
NASA NTRS1983

Planetary geodesy

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.

Sjogren, W. L.↗
NASA NTRS1980

Equivalent source modeling of the main field using Magsat data

A graphic software package was developed to plot the dipole positions for a particular model on a world map. An arrow is drawn at each dipole in the direction of the horizontal magnetization vector, with length proportional to the horizontal magnitude. A contouring package represents the radical component of the magnetization vector. Based on quiet Magsat data, several dipole models were derived. These are being verified against each other and against MGST (6/80) data set as to goodness of fit and stability of solution. The spatial power spectra computed from spherical harmonic expansions of the dipole models is being analyzed relative to crustal and core content. The power spectra for one of the dipole models and a spherical harmonic model based on Magsat data through degree and order 23 is presented.

Source record↗