Search NASA⌕ Search

SEARCH · Search NASA

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 109 records · Page 6

The barotropic normal modes in certain shear flows and the traveling waves in the atmosphere

It is shown analytically and numerically that in certain shear flows the linearized nondivergent barotropic vorticity equation has a limited number of neutral normal modes. The latitudinal structures of these shear flows can be expressed as polynomials of the sine of latitude. The first few such shear flows resemble the gross features of the zonal winds in the atmosphere of the earth at different times and altitudes. The spatial structures of the neutral normal modes in these shear flows are spherical harmonics, and, as a consequence, these modes are also the exact solutions of the fully nonlinear equation because the nonlinear interaction term vanishes identically. The spatial structures of the observed 5-, 4-, 2-, and 16-day free traveling waves in the atmosphere are often identified with the spherical harmonics with indices of (m, n) = ( 1, 2), (2, 3), (3, 3), and ( 1, 4), which are known previously as the neutral normal modes of the nondivergent barotropic vorticity equation in a motionless background state. Our results could explain why these free traveling waves can survive the shearing effects of zonal flows that are far different from rest because these spherical harmonics are also normal modes in certain shear flows that resemble the observations of the atmosphere.

Chen, Ping↗

Use of non-spherical gravity harmonics for relative motion GN/C (category 2)

Detailed analysis of the Automatic Rendezvous and Capture problem indicate a need for three different regions of mathematical description for the GN&C algorithms: (1) multi-vehicle orbital mechanics to the rendezvous interface point, i.e., within 100 nm; (2) relative motion solutions (such as Clohessy-Wiltshire type) from the far-field to the near-field interface, i.e., within 1 nm and; (3) close proximity motion - the near-field motion where the relative differences in the gravitational and orbit inertial accelerations can be neglected from the equations of motion. Limit boundaries to these regions can be precisely defined by further analysis and will be functions of the tracking measurement accuracies and the computer resources available for the solution of the algorithms. This paper analyzes the relative motion in Regions 2 and 3 above and present the derivation and discussion of the general case of non-spherical gravitational perturbed relative motion. Mathematical deviations from the numerically integrated spherical gravity case and solutions from the Clohessy-Wiltshire equations are presented in the analysis. Based upon this preliminary analysis, it is recommended that further efforts be used to assess the relative position and velocity differences in Region 2 due to non-spherical gravity harmonics and that viable GN&C algorithms be developed to include these gravity perturbations (especially the effects of the first gravity harmonic, J2).

Henderson, David M.↗

A New Model of Jupiter's Magnetic Field From Juno's First Nine Orbits

A spherical harmonic model of the magnetic field of Jupiter is obtained from vector magnetic field observations acquired by the Juno spacecraft during its first nine polar orbits about the planet. Observations acquired during eight of these orbits provide the first truly global coverage of Jupiter's magnetic field with a coarse longitudinal separation of ~45 deg between perijoves. The magnetic field is represented with a degree 20 spherical harmonic model for the planetary ("internal") field, combined with a simple model of the magnetodisc for the field ("external") due to distributed magnetospheric currents. Partial solution of the underdetermined inverse problem using generalized inverse techniques yields a model ("Juno Reference Model through Perijove 9") of the planetary magnetic field with spherical harmonic coefficients well determined through degree and order 10, providing the first detailed view of a planetary dynamo beyond Earth.

Connerney, J. E. P.↗

Scalar magnetic anomaly maps of Earth derived from POGO and Magsat data

A new Polar Orbit Geophysical Observatory (POGO) scalar magnetic anomaly map at 400 km altitude is presented which consists of spherical harmonics of degree 15-60. On the basis of the common features of this map with two new Magsat anomaly maps, dawn and dusk, two scalar magnetic anomaly maps of the Earth are presented using two selection criteria with different levels of stringency. These selection criteria suppress the noncrustal components of the original maps by different amounts. The more stringent selection criteria seek to eliminate as much contamination as possible, at the expense of suppressing some anomaly signal. This map is represented by spherical harmonics of degree 15-60. The less stringent selection criteria seek to retain as much crustal signal as possible, at the expense of also retaining some contaminating fields. This map is represented by spherical harmonics of degree 15-65. The resulting two maps are highly correlated with degree correlation coefficients greater than 0.8.

Arkani-Hamed, Jafar↗

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.

Venus venus topgraphy harmonic analysis spherical ↗

Dynamic ocean-tide effects on Earth's rotation

This article develops 'broad-band' Liouville equations which are capable of determining the effects on the rotation of the Earth of a periodic excitation even at frequencies as high as semi-diurnal; these equations are then used to predict the rotational effects of altimetric, numerical and 32-constituent spherical harmonic ocean-tide models. The rotational model includes a frequency-dependent decoupled core, the effects of which are especially marked near retrograde diurnal frequencies; and a fully dynamic oceanic response, whose effects appear to be minor despite significant frequency dependence. The model also includes solid-earth effects which are frequency dependent as the result of both anelasticity at long periods and the fluid-core resonance at nearly diurnal periods. The effects of both tidal inertia and relative angular momentum on Earth rotation (polar motion, length of day, 'nutation' and Universal Time) are presented for 32 long- and short-period ocean tides determined as solutions to the author's spherical harmonic tide theory. The lengthening of the Chandler wobble period by the pole tide is also re-computed using the author's full theory. Additionally, using the spherical harmonic theory, tidal currents and their effects on rotation are determined for available numerical and altimetric tide height models. For all models, we find that the effects of tidal currents are at least as important as those of tide height for diurnal and semi-diurnal constituents.

Dickman, S. R.↗

Jupiter's magnetic field, magnetosphere, and interaction with the solar wind - Pioneer 11

Measurements of the magnetic field vector were obtained continuously throughout the encounter of the spacecraft with the planet Jupiter. Effects of Jupiter on the solar wind are considered along with the characteristics of the magnetopause at both low and intermediate latitudes, the three basic regions within the magnetosphere, and a spherical harmonic analysis of the Pioneer 11 measurements. The spherical harmonic representation has been used to derive contours of the magnetic field strength at the surface of Jupiter.

Smith, E. J.↗

Some new methods in geomagnetic field modeling applied to the 1960 - 1980 epoch

The utilization of satellite and surface data together permitted the incorporation of a solution for the anomaly field at each observatory. The residuals of the observatory measurements to such models is commensurate with the actual measurment accuracy. Incorporation of the anomaly estimation enabled the inclusion of stable time derivatives of the spherical harmonic coefficients up to the third derivative. A spherical harmonic model is derived with degree and order 13 in its constant and first time derivative terms, six in its second derivative terms and four in its third derivative terms.

Langel, R. A.↗

Unique determination of model coronal magnetic fields using photospheric observations

We show that the nonradial field-boundary condition (or the line-of-sight boundary condition) for the Laplacian-like equation developed by Bogdan and Low (1986) is sufficient to uniquely determine the model coronal magnetic field provided the electric currents are horizontal (or zero, the current-free case) at the solar surface as well as in the solar atmosphere between the photosphere and the source surface. The derived recursion formulae for the spherical harmonic coefficients can be used to determine the spherical harmonic coefficients in the solutions of the horizontal current models very efficiently.

Zhao, Xuepu↗

Anisotropic Developments for Homogeneous Shear Flows

The general decomposition of the spectral correlation tensor R(sub ij)(k) by Cambon et al. (J. Fluid Mech., 202, 295; J. Fluid Mech., 337, 303) into directional and polarization components is applied to the representation of R(sub ij)(k) by spherically averaged quantities. The decomposition splits the deviatoric part H(sub ij)(k) of the spherical average of R(sub ij)(k) into directional and polarization components H(sub ij)(sup e)(k) and H(sub ij)(sup z)(k). A self-consistent representation of the spectral tensor in the limit of weak anisotropy is constructed in terms of these spherically averaged quantities. The directional polarization components must be treated independently: models that attempt the same representation of the spectral tensor using the spherical average H(sub ij)(k) alone prove to be inconsistent with Navier-Stokes dynamics. In particular, a spectral tensor consistent with a prescribed Reynolds stress is not unique. The degree of anisotropy permitted by this theory is restricted by realizability requirements. Since these requirements will be less severe in a more accurate theory, a preliminary account is given of how to generalize the formalism of spherical averages to higher expansion of the spectral tensor. Directionality is described by a conventional expansion in spherical harmonics, but polarization requires an expansion in tensorial spherical harmonics generated by irreducible representations of the spatial rotation group SO(exp 3). These expansions are considered in more detail in the special case of axial symmetry.

Cambon, Claude↗

Normalization of Gravitational Acceleration Models

Unlike the uniform density spherical shell approximations of Newton, the con- sequence of spaceflight in the real universe is that gravitational fields are sensitive to the nonsphericity of their generating central bodies. The gravitational potential of a nonspherical central body is typically resolved using spherical harmonic approximations. However, attempting to directly calculate the spherical harmonic approximations results in at least two singularities which must be removed in order to generalize the method and solve for any possible orbit, including polar orbits. Three unique algorithms have been developed to eliminate these singularities by Samuel Pines [1], Bill Lear [2], and Robert Gottlieb [3]. This paper documents the methodical normalization of two1 of the three known formulations for singularity-free gravitational acceleration (namely, the Lear [2] and Gottlieb [3] algorithms) and formulates a general method for defining normalization parameters used to generate normalized Legendre Polynomials and 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.↗

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.↗

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.↗

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.↗

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.↗

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↗

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.↗