Search NASA⌕ Search

SEARCH · Search NASA

Results for “generalized 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 73 records · Page 4

The non-hydrostatic figures of the terrestrial planets

Solid state creep being exponentially dependent on temperature must dominate the mechanical behavior of the mantles of terrestrial planets beneath their lithospheres. General arguments suggest that the lithospheres of the Moon and Mars are about 200 km thick; the Earth, Venus and Mercury much less. Short wavelength gravity anomalies are explained by the finite strength of the lithosphere: the lunar mascons being an example. The good correlation of the Venus and Mars gravity anomalies with topography up to spherical harmonics of degrees 10-15 is in striking contrast to the lack of correlation between the long wavelength components of the geoid and the continent-ocean distribution or even the plates. Attempts have been made to explain the former correlations by isostatic models but the depths of compensation seem implausible. Low degree harmonics of the gravity fields of the terrestrial planets as is certainly the case in the Earth must arise from the density variations driving solid state convection. In the case of Venus the less dense differentiated materials of the highlands seems to be positioned over the singular points of the convection pattern. Thus the correlated gravity field does not arise from the highlands but from the density difference in the convecting interior. In the Earth lack of correlation seems to arise from the fact that the plates have moved relative to the convection pattern the last 100 M yr.

Runcorn, S. K.↗

A Theoretical Analysis of the Effects of Fuel Motion on Airplane Dynamics

The general equations of motion for an airplane with a number of spherical fuel tanks are derived. The motion of the fuel is approximated by the motion of solid pendulums. The same type of derivation and equations are shown to apply to any type of fuel tank where the motion of the fuel may be represented in terms of undamped harmonic oscillators. Motions are calculated for a present-day high-speed airplane and a free-flying airplane model with two spherical tanks in the symmetry plane.

Schy, Albert A↗

Flyby Characterization of Lower-Degree Spherical Harmonics Around Small Bodies

Interest in studying small bodies has grown significantly in the last two decades, and there are a number of past, present, and future missions. These small body missions challenge navigators with significantly different kinds of problems than the planets and moons do. The small bodies' shape is often irregular and their gravitational field significantly weak, which make the designing of a stable orbit a complex dynamical problem. In the initial phase of spacecraft rendezvous with a small body, the determination of the gravitational parameter and lower-degree spherical harmonics are of crucial importance for safe navigation purposes. This motivates studying how well one can determine the total mass and lower-degree spherical harmonics in a relatively short time in the initial phase of the spacecraft rendezvous via flybys. A quick turnaround for the gravity data is of high value since it will facilitate the subsequent mission design of the main scientific observation campaign. We will present how one can approach the problem to determine a desirable flyby geometry for a general small body. We will work in the non-dimensional formulation since it will generalize our results across different size/mass bodies and the rotation rate for a specific combination of gravitational coefficients.

Takahashi, Yu↗

Diffusion model for lightning radiative transfer

A one-speed Boltzmann transport theory, with diffusion approximations, is applied to study the radiative transfer properties of lightning in optically thick thunderclouds. Near-infrared (lambda = 0.7774 micrometers) photons associated with a prominent oxygen emission triplet in the lightning spectrum are considered. Transient and spatially complex lightning radiation sources are placed inside a rectangular parallelepiped thundercloud geometry and the effects of multiple scattering are studied. The cloud is assumed to be composed of a homogeneous collection of identical spherical water droplets, each droplet a nearly conservative, anisotropic scatterer. Conceptually, we treat the thundercloud like a nuclear reactor, with photons replaced by neutrons, and utilize standard one-speed neutron diffusion techniques common in nuclear reactor analyses. Valid analytic results for the intensity distribution (expanded in spherical harmonics) are obtained for regions sufficiently far from sources. Model estimates of the arrival-time delay and pulse width broadening of lightning signals radiated from within the cloud are determined and the results are in good agreement with both experimental data and previous Monte Carlo estimates. Additional model studies of this kind will be used to study the general information content of cloud top lightning radiation signatures.

Koshak, William J.↗

The isostatic state of Mead crater

We have analyzed high-resolution Magellan Doppler tracking data over Mead crater, using both line-of-sight and spherical harmonic methods, and have found a negative gravity anomaly of about 4-5 mgal (at spacecraft altitude, 182 km). This is consistent with no isostatic compensation of the present topography; the uncertainty in the analysis allows perhaps as much as 30% compensation at shallow dpeths (approximately 25 km). This is similar to observations of large craters on Earth, which are not generally compensated, but contrasts with at least some lunar basins which are inferred to have large Moho uplifts and corresponding positive Bouguer anomalies. An uncompensated load of this size requires a lithosphere with an effective elastic lithosphere thickness greater than 30 km. In order for the crust-mantle boundary not to have participated in the deformation associated with the collapse of the transient cavity during the creation of the crater, the yield strength near the top of the mantle must have been significantly higher on Earth and Venus than on the Moon at the time of basin formation. This might be due to increased strength against frictional sliding at the higher confining pressures within the larger planets. Alternatively, the thinner crusts of Earth and Venus compared to that of the Moon may result in higher creep strength of the upper mantle at shallower depths.

Banerdt, W. B.↗

A numerical modeling study of the interaction between the tides and the circulation forced by high-latitude plasma convection

A spectral, time-varying thermospheric general circulation model has been used to study the nonlinear interaction at high latitudes between the tides propagating into the thermosphere from below and the circulation induced by magnetospheric forcing and in situ solar heating. The model is discrete in the vertical with 27 layers spaced by half a scale height. In the horizontal, the fields are expanded in a series of spherical harmonics using a triangular truncation at wave number 31, equivalent to a homogeneous global resolution with a minimum wavelength of 1270 km. A hypothetical uniform grid point model would require a horizontal spacing of 417 km to describe the same minimum wavelength. In the high-latitude F region the tides affect the dusk vortex of the neutral flow very little, but the dawn vortex is either suppressed or amplified dependent upon the universal time and tidal phase. In the E region neutral flow, both the dusk and dawn vortices are shifted in local time by the tides, again as a function of universal time and tidal phase. At dusk a nonlinear amplification of the sunward winds occurs for certain combinations of parameters, and at dawn the winds may be completely suppressed. Below 120 km altitude the magnetospheric forcing creates a single cyclonic vortex which is also sensitive to the high-latitude tidal structure.

Mikkelsen, I. S.↗

Method for predicting impulsive noise generated by wind turbine rotors

Large wind turbines can generate both broad band and impulsive noises. These noises can be controlled by proper choice of rotor design parameters such as rotor location with respect to the supporting tower, tower geometry and tip speed. A method was developed to calculate the impulsive noise generated when the wind turbine blade experiences air forces that are periodic functions of the rotational frequency. This phenomenon can occur when the blades operate in the wake of the support tower and the nonuniform velocity field near the ground due to wind shear. Results from this method were compared with measured sound spectra taken at locations of one to two rotor diameters from the DOE/NASA Mod-1 wind turbine. The calculated spectra generally agreed with the measured data in both the amplitude of the predominant harmonics and the roll of rate with frequency. Measured sound pressure levels far from the Mod-1 (15 rotor diameters), however, were higher than predicted. Simultaneous measurements in the near and far field indicated that the propagation effects could enhance the sound levels by more than 10 dB above that expected by spherical dispersion. These propagation effects are believed to be due to terrain and atmospheric characteristics of the Mod-1 site.

Viterna, L. A.↗

Solar oscillations - A method for deriving nonlinear effects

The frequencies of solar oscillations are so closely spaced that nonlinear interactions among modes are probable. The rate of interaction is proportional to an integral involving the eigenfunctions of the interacting modes, which usually are known only numerically. An approximation in which the eigenfunctions are strictly sinusoidal functions of a suitably defined radial variable, and the numerical details of stellar structure are banished to a coefficient in the integrand, are here explored. The physical assumptions are the same as in the asymptotic approximation of p- or g-modes. This method should allow a general investigation as to likely nonlinear interactions and which modes may participate in such interactions. The coupling of two p-modes by possible large-scale internal magnetic fields is introduced as an anisotropic pressure response to a displacement. Pairs of modes differing in frequency by less than the fraction magnetic/thermal pressure are strongly coupled, and energy appears to oscillate slowly between the two associated spherical harmonics. Potentially, upper limits may be derived for internal magnetic fields.

Wentzel, Donat G.↗

Upper limit on solar interior rotation

Individual 5-min p-mode oscillations of spherical harmonic degree 0-2 and radial order 16-26 are revealed by the Solar Maximum Mission Spacecraft's Active Cavity Radiometer Irradiance Monitor power spectra for total solar irradiance flux variation. Rotationally split p-, g- and f-modes have been identified in the temporal power spectrum of Bos and Hill's (1983) limb-darkening data, and from these splittings, two rotational curves have been deduced which imply a solar gravitational quadrupole moment sufficiently large to preclude agreement between general relativity and planetary motion observations.

Woodard, M.↗

The Ohio State 1991 geopotential and sea surface topography harmonic coefficient models

The computation is described of a geopotential model to deg 360, a sea surface topography model to deg 10/15, and adjusted Geosat orbits for the first year of the exact repeat mission (ERM). This study started from the GEM-T2 potential coefficient model and it's error covariance matrix and Geosat orbits (for 22 ERMs) computed by Haines et al. using the GEM-T2 model. The first step followed the general procedures which use a radial orbit error theory originally developed by English. The Geosat data was processed to find corrections to the a priori geopotential model, corrections to a radial orbit error model for 76 Geosat arcs, and coefficients of a harmonic representation of the sea surface topography. The second stage of the analysis took place by doing a combination of the GEM-T2 coefficients with 30 deg gravity data derived from surface gravity data and anomalies obtained from altimeter data. The analysis has shown how a high degree spherical harmonic model can be determined combining the best aspects of two different analysis techniques. The error analysis was described that has led to the accuracy estimates for all the coefficients to deg 360. Significant work is needed to improve the modeling effort.

Rapp, Richard H.↗

Optimum data weighting and error calibration for estimation of gravitational parameters

A new technique was developed for the weighting of data from satellite tracking systems in order to obtain an optimum least squares solution and an error calibration for the solution parameters. Data sets from optical, electronic, and laser systems on 17 satellites in GEM-T1 (Goddard Earth Model, 36x36 spherical harmonic field) were employed toward application of this technique for gravity field parameters. Also, GEM-T2 (31 satellites) was recently computed as a direct application of the method and is summarized here. The method employs subset solutions of the data associated with the complete solution and uses an algorithm to adjust the data weights by requiring the differences of parameters between solutions to agree with their error estimates. With the adjusted weights the process provides for an automatic calibration of the error estimates for the solution parameters. The data weights derived are generally much smaller than corresponding weights obtained from nominal values of observation accuracy or residuals. Independent tests show significant improvement for solutions with optimal weighting as compared to the nominal weighting. The technique is general and may be applied to orbit parameters, station coordinates, or other parameters than the gravity model.

Lerch, F. J.↗

Estimation and Validation of Oceanic Mass Circulation from the GRACE Mission

Since the launch of the Gravity Recovery And Climate Experiment (GRACE) in March 2002, the Earth's surface mass variations have been monitored with unprecedented accuracy and resolution. Compared to the classical spherical harmonic solutions, global high-resolution mascon solutions allows the retrieval of mass variations with higher spatial and temporal sampling (2 degrees and 10 days). We present here the validation of the GRACE global mascon solutions by comparing mass estimates to a set of about 100 ocean bottom pressure (OSP) records, and show that the forward modelling of continental hydrology prior to the inversion of the K-band range rate data allows better estimates of ocean mass variations. We also validate our GRACE results to OSP variations modelled by different state-of-the-art ocean general circulation models, including ECCO (Estimating the Circulation and Climate of the Ocean) and operational and reanalysis from the MERCATOR project.

Boy, J.-P.↗

Experimental study of outdoor propagation of spherically speading periodic acoustic waves of finite amplitude

The outdoor propagation of spherically spreading sound waves of finite amplitude was investigated. The main purpose of the experiments was to determine the extent to which the outdoor environment, mainly random inhomogeneity of the medium, affects finite amplitude propagation. Periodic sources with fundamental frequencies in the range 6 to 8 kHz and source levels SPLlm from 140 to 149 dB were used. The sources were an array of 7 to 10 horn drivers and a siren. The propagation path was vertical and parallel to an 85 m tower, whose elevator carried the traveling microphone. The general conclusions drawn from the experimental results were as follows. The inhomogeneities caused significant fluctuations in the instantaneous acoustic signal, but with sufficient time averaging of the measured harmonic levels, the results were comparable to results expected for propagation in a quiet medium. Propagation data for the fundamental of the siren approached within 1 dB of the weak shock saturation levels. Extra attenuation on the order of 8 dB was observed. The measurements generally confirmed the predictions of several theoretical models. The maximum propagation distance was 36 m. The narrowbeam arrays were much weaker sources. Nonlinear propagation distortion was produced, but the maximum value of extra attenuation measured was 1.5 dB. The maximum propagation distance was 76 m. The behavior of the asymetric waveforms received in one experiment qualitatively suggested that beam type diffraction effects were present. The role of diffraction of high intensity sound waves in radiation from a single horn was briefly investigated.

Theobald, M. A.↗

Topography of the Moon from the Clementine Lidar

Range measurements from the lidar instrument carried aboard the Clementine spacecraft have been used to produce an accurate global topographic model of the Moon. This paper discusses the function of the lidar; the acquisition, processing, and filtering of observations to produce a global topographic model; and the determination of parameters that define the fundamental shape of the Moon. Our topographic model: a 72nd degree and order spherical harmonic expansion of lunar radii, is designated Goddard Lunar Topography Model 2 (GLTM 2). This topographic field has an absolute vertical accuracy of approximately 100 m and a spatial resolution of 2.5 deg. The field shows that the Moon can be described as a sphere with maximum positive and negative deviations of approx. 8 km, both occurring on the farside, in the areas of the Korolev and South Pole-Aitken (S.P.-Aitken) basins. The amplitude spectrum of the topography shows more power at longer wavelengths as compared to previous models, owing to more complete sampling of the surface, particularly the farside. A comparison of elevations derived from the Clementine lidar to control point elevations from the Apollo laser altimeters indicates that measured relative topographic heights generally agree to within approx. 200 in over the maria. While the major axis of the lunar gravity field is aligned in the Earth-Moon direction, the major axis of topography is displaced from this line by approximately 10 deg to the cast and intersects the farside 24 deg north of the equator. The magnitude of impact basin topography is greater than the lunar flattening (approx. 2 km) and equatorial ellipticity (approx. 800 m), which imposes a significant challenge to interpreting the lunar figure. The floors of mare basins are shown to lie close to an equipotential surface, while the floors of unflooded large basins, except for S.P.-Aitken, lie above this equipotential. The radii of basin floors are thus consistent with a hydrostatic mechanism for the absence of significant farside maria except for S.P.-Aitken, whose depth and lack of mare require significant internal compositional and/or thermal heterogeneity. A macroscale surface roughness map shows that roughness at length scales of 10(exp 1) - 10(exp 2) km correlates with elevation and surface age.

Smith, David E.↗

A general low frequency acoustic radiation capability for NASTRAN

A new capability called NASHUA is described for calculating the radiated acoustic sound pressure field exterior to a harmonically-excited arbitrary submerged 3-D elastic structure. The surface fluid pressures and velocities are first calculated by coupling a NASTRAN finite element model of the structure with a discretized form of the Helmholtz surface integral equation for the exterior fluid. After the fluid impedance is calculated, most of the required matrix operations are performed using the general matrix manipulation package (DMAP) available in NASTRAN. Far field radiated pressures are then calculated from the surface solution using the Helmholtz exterior integral equation. Other output quantities include the maximum sound pressure levels in each of the three coordinate planes, the rms and average surface pressures and normal velocities, the total radiated power and the radiation efficiency. The overall approach is illustrated and validated using known analytic solutions for submerged spherical shells subjected to both uniform and nonuniform applied loads.

Everstine, G. C.↗

Results of dynamic testing of GP-B spherical gyroscopes

Laboratory tests of the spherical electrostatically levitated cryogenically cooled coated gyroscope being developed for the Gravity Probe B (GP-B) spacecraft (Bardas et al., 1986) are reported. Spin speed and the dc components of the trapped magnetic field are measured with three orthogonal pickup loops attached to SQUID detectors as the levitated gyro is brought up to speed by an He gas jet. Data on the spin-vector time history, mass unbalance, higher rotor-shape harmonics, and spin-vector position are presented in extensive graphs and characterized in detail, and a mathematical model of the electrostatic suspension torques is derived. Prototype gyro 86-4 is found to have mass unbalance within the range required for the GP-B mission (to detect the geodetic and motional effects predicted by general relativity theory).

Keiser, G. M.↗

Viscoelastic responses of a hard transition zone - Effects on postglacial uplifts and rotational signatures

A five-layer viscoelastic spherical model is used to calculate the transient displacements of postglacial rebound, the induced polar motions, and the temporal variations of the geopotential up to degree 8 of the zonal coefficients. Two models - one with two viscoelastic layers separated at 670 km, and the other with three layers in which a hard garnet layer lies between the upper and lower mantle - are compared. Forward modeling shows that it may be possible to discern the presence of a hard garnet layer with a viscosity of at least ten times greater than the upper mantle, on the basis of uplift data near the center of the former Laurentide ice-sheet and from polar wander and j2 data. Temporal variations of higher gravity harmonics, such as j6 and j8, can potentially place even tighter constraints on the rheological properties of the hard transition zone. A lower mantle viscosity between 2 and 4 x 10 to the 22nd Pa is generally preferred in models with a garnet layer which may be as large as 50 times more viscous than the upper mantle.

Spada, Giorgio↗

On Lateral Viscosity Contrast in the Mantle and the Rheology of Low-Frequency Geodynamics

Mantle-wide heterogeneity is largely controlled by deeply penetrating thermal convective currents. These thermal currents are likely to produce significant lateral variation in rheology, and this can profoundly influence overall material behaviour. How thermally related lateral viscosity variations impact models of glacio-isostatic and tidal deformation is largely unknown. An important step towards model improvement is to quantify, or bound, the actual viscosity variations that characterize the mantle. Simple scaling of viscosity to shear-wave velocity fluctuations yields map-views of long- wavelength viscosity variation. These give a general quantitative description and aid in estimating the depth dependence of rheological heterogeneity throughout the mantle. The upper mantle is probably characterized by two to four orders of magnitude variation (peak-to-peak). Discrepant time-scales for rebounding Holocene shorelines of Hudson Bay and southern Iceland are consistent with this characterization. Results are given in terms of a local average viscosity ratio, (Delta)eta(bar)(sub i), of volumetric concentration, phi(sub i). For the upper mantle deeper than 340 km the following reasonable limits are estimated for (delta)eta(bar) approx. equal 10(exp -2): 0.01 less than or equal to phi less than or equal to 0.15. A spectrum of ratios (Delta)eta(bar)(sub i) less than 0.1 at concentration level eta(sub i) approx. equal 10(exp -6) - 10(exp -1) in the lower mantle implies a spectrum of shorter time-scale deformational response modes for second-degree spherical harmonic deformations of the Earth. Although highly uncertain, this spectrum of spatial variation allows a purely Maxwellian viscoelastic rheology simultaneously to explain all solid tidal dispersion phenomena and long-term rebound-related mantle viscosity. Composite theory of multiphase viscoelastic media is used to demonstrate this effect.

Ivins, Erik R.↗