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 685 records · Page 38

A Numerical Method for Obtaining Monoenergetic Neutron Flux Distributions and Transmissions in Multiple-Region Slabs

This method is investigated for semi-infinite multiple-slab configurations of arbitrary width, composition, and source distribution. Isotropic scattering in the laboratory system is assumed. Isotropic scattering implies that the fraction of neutrons scattered in the i(sup th) volume element or subregion that will make their next collision in the j(sup th) volume element or subregion is the same for all collisions. These so-called "transfer probabilities" between subregions are calculated and used to obtain successive-collision densities from which the flux and transmission probabilities directly follow. For a thick slab with little or no absorption, a successive-collisions technique proves impractical because an unreasonably large number of collisions must be followed in order to obtain the flux. Here the appropriate integral equation is converted into a set of linear simultaneous algebraic equations that are solved for the average total flux in each subregion. When ordinary diffusion theory applies with satisfactory precision in a portion of the multiple-slab configuration, the problem is solved by ordinary diffusion theory, but the flux is plotted only in the region of validity. The angular distribution of neutrons entering the remaining portion is determined from the known diffusion flux and the remaining region is solved by higher order theory. Several procedures for applying the numerical method are presented and discussed. To illustrate the calculational procedure, a symmetrical slab ia vacuum is worked by the numerical, Monte Carlo, and P(sub 3) spherical harmonics methods. In addition, an unsymmetrical double-slab problem is solved by the numerical and Monte Carlo methods. The numerical approach proved faster and more accurate in these examples. Adaptation of the method to anisotropic scattering in slabs is indicated, although no example is included in this paper.

Schneider, Harold↗

An Intercomparison of the Dynamical Cores of Global Atmospheric Circulation Models for Mars

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Center and San Jose State University, Department of Meteorology. The focus of this JRI has been to evaluate the dynamical 'cores' of two global atmospheric circulation models for Mars that are in operation at the NASA Ames Research Center. The two global circulation models in use are fundamentally different: one uses spherical harmonics in its horizontal representation of field variables; the other uses finite differences on a uniform longitude-latitude grid. Several simulations have been conducted to assess how the dynamical processors of each of these circulation models perform using identical 'simple physics' parameterizations. A variety of climate statistics (e.g., time-mean flows and eddy fields) have been compared for realistic solstitial mean basic states. Results of this research have demonstrated that the two Mars circulation models with completely different spatial representations and discretizations produce rather similar circulation statistics for first-order meteorological fields, suggestive of a tendency for convergence of numerical solutions. Second and higher-order fields can, however, vary significantly between the two models.

Hollingsworth, Jeffery L.↗

An Intercomparison of the Dynamical Cores of Global Atmospheric Circulation Models for Mars

This is a Final Report for a Joint Research Interchange (JRI) between NASA Ames Research Cen- ter and San Jose State University, Department of Meteorology. The focus of this JRI has been to evaluate the dynamical "cores" of two global atmospheric circulation models for Mars that are in operation at the NASA Ames Research Center. ne two global circulation models in use are fundamentally different: one uses spherical harmonics in its horizontal representation of field variables; the other uses finite differences on a uniform longitude-latitude grid. Several simulations have been conducted to assess how the dynamical processors of each of these circulation models perform using identical "simple physics" parameterizations. A variety of climate statistics (e.g., time-mean flows and eddy fields) have been compared for realistic solstitial mean basic states. Results of this research have demonstrated that the two Mars circulation models with completely different spatial representations and discretizations produce rather similar circulation statistics for first-order meteorological fields, suggestive of a tendency for convergence of numerical solutions. Second and higher-order fields can, however, vary significantly between the two models.

Hollingsworth, Jeffery L.↗

A Search for r-Modes from 1825 to the Present

Global oscillations (r-modes) of the Sun's outer convective envelope with periods approximately 1 month and longer have been detected in several short data strings of several years duration. To test whether r-modes might persist beyond one 11 year cycle, the daily sunspot numbers from 1825 to the present were analyzed. Good evidence, but confidence level less than 3sigma, was found for most of the 14 r-modes with spherical harmonic index lambda less than or equal to 5 that can exist in the presence of solar differential rotation. The characteristic rotation rate of almost every such r-mode was detected, displaced systematically from its expected value by only 0.15%. If this probable detection is real, then most low harmonic r-modes have lifetimes exceeding one century and the rotation of the Sun's outer layers varies by less than 0.05%, except possibly at solar minimum.

Wolff, Charles L.↗

An Improved Magnetic Field Model for Jupiter's Inner Magnetosphere Using a Microsignature of Amalthea

Observation of a particle absorption microsignature of the Jovian satellite Amalthea during the Pioneer II close flyby of Jupiter on December 3, 1974, has been described by McKibben et al. The microsignature was also observed by the University of Iowa/Pioneer II instrument in the distribution of protons in the kinetic energy range 0.61 less than T(sub p) less than 3.41 MeV but has not been previously reported. The finer time resolution and superior data quality of the latter observations provide a fresh basis for assessing the accuracy with which various published magnetic field models of Jupiter's magnetic field models describe Jupiter's magnetic field for radial distances of the order of or less than several planetary radii. The expected time of occurrence of the minimum of Amalthea's microsignature was calculated for each model and compared to the time of the observed minimum. The discrepancy between these two times was significant for each of the published models, but it could be reduced to zero by interpolating between two of the closely related models. The spherical harmonic coefficients of this combination model are tabulated. Also included is a novel method for calculating the intensity-time profile of the microsignature. In a later paper, the improved model is utilized to analyze the energetic particle measurements in Jupiter's innermost magnetosphere by the Galileo entry probe in December 1995.

Randall, Bruce A.↗

Seasonal Variations of the Earth's Gravitational Field: An Analysis of Atmospheric Pressure, Ocean Tidal, and Surface Water Excitation

Monthly mean gravitational field parameters (denoted here as C(sub even)) that represent linear combinations of the primarily even degree zonal spherical harmonic coefficients of the Earth's gravitational field have been recovered using LAGEOS I data and are compared with those derived from gridded global surface pressure data of the National meteorological center (NMC) spanning 1983-1992. The effect of equilibrium ocean tides and surface water variations are also considered. Atmospheric pressure and surface water fluctuations are shown to be the dominant cause of observed annual C(sub even) variations. Closure with observations is seen at the 1sigma level when atmospheric pressure, ocean tide and surface water effects are include. Equilibrium ocean tides are shown to be the main source of excitation at the semiannual period with closure at the 1sigma level seen when both atmospheric pressure and ocean tide effects are included. The inverted barometer (IB) case is shown to give the best agreement with the observation series. The potential of the observed C(sub even) variations for monitoring mass variations in the polar regions of the Earth and the effect of the land-ocean mask in the IB calculation are discussed.

Dong, D,↗

A Global Ocean Tide Model From TOPEX/POSEIDON Altimetry: GOT99.2

Goddard Ocean Tide model GOT99.2 is a new solution for the amplitudes and phases of the global oceanic tides, based on over six years of sea-surface height measurements by the TOPEX/POSEIDON satellite altimeter. Comparison with deep-ocean tide-gauge measurements show that this new tidal solution is an improvement over previous global models, with accuracies for the main semidiurnal lunar constituent M2 now below 1.5 cm (deep water only). The new solution benefits from use of prior hydrodynamic models, several in shallow and inland seas as well as the global finite-element model FES94.1. This report describes some of the data processing details involved in handling the altimetry, and it provides a comprehensive set of global cotidal charts of the resulting solutions. Various derived tidal charts are also provided, including tidal loading deformation charts, tidal gravimetric charts, and tidal current velocity (or transport) charts. Finally, low-degree spherical harmonic coefficients are computed by numerical quadrature and are tabulated for the major short-period tides; these are useful for a variety of geodetic and geophysical purposes, especially in combination with similar estimates from satellite laser ranging.

Ray, Richard D.↗

Lunar and Solar Torques on the Oceanic Tides

Brosche and Seiler recently suggested that direct lunar and solar tidal torques on the oceanic tides play a significant role in the earth's short-period angular momentum balance ("short-period" here meaning daily and sub-daily). We reexamine that suggestion here, concentrating on axial torques and hence on variations in rotation rate. Only those spherical harmonic components of the ocean tide having the same degree and order as the tidal potential induce nonzero torques. Prograde components (those moving in the same direction as the tide-generating body) produce the familiar secular braking of the earth's rotation. Retrograde components, however, produce rapid variations in UTI at twice the tidal frequency. There also exist interaction torques between tidal constituents, e.g. solar torques on lunar tides. They generate UTI variations at frequencies equal to the sums and differences of the original tidal frequencies. We give estimates of the torques and angular momentum variations for each of the important regimes, secular to quarter-diurnal. For the M(sub 2) potential acting on the M(sub 2) ocean tide, we find an associated angular momentum variation of amplitude 3 x 10(exp 19) N m. This is 5 to 6 orders of magnitude smaller than the angular momentum variations associated with tidal currents. We conclude that these torques do not play a significant role in the short-period angular momentum balance.

Ray, Richard D.↗

Gravitational Constraints on the Internal Structure of Ganymede

Before the arrival of the Galileo spacecraft in the jovian system, there was little information on the interior structure of jupiter's largest moon, Ganymede. Its mean density (1,940 kg/cu m), determined by the Pioneer and Voyager spacecraft, implies a composition that is roughly 60% rock and 40% ice, which could be uniformly mixed or differentiated into a rocky core and icy mantle. Here we report measurements by the Galileo spacecraft of Ganymede's overall density and the spherical harmonics, J(sub 2) and C(sub 22) of its gravitational field. These data show clearly that Ganymede has differentiated into a core and mantle. Combined with the recent discovery of an intrinsic magnetic field, our gravity results suggest that Ganymede has a metallic core of radius 400-1,300 km surrounded by a silicate mantle, which is in turn enclosed by an ice she approx.800 km thick. Depending on whether the core is pure iron or an alloy of iron and iron sulphide, it could account for as little as 1.4% or as much as one-third of the total mass. If the ice were stripped away, Ganymede could look much like lo' in terms of its size and internal mass distribution.

Anderson, J. D.↗

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

Mantle Q at Period 12.4 Hours: New Estimates from Satellite Tracking and Altimetry

The most convincing estimates of mantle Q at periods of many hours have historically come from extrapolating seismic and free-oscillation estimates via some assumed frequency dependence, sometimes contrained by estimates from the Chandler Wobble. At the semidiurnal tidal period, direct estimates of Q have been difficult to obtain because of the dominating signals of the ocean tides, which account for more than 95!k of the tidal energy dissipation. But knowledge of the ocean tides has been rapidly improving, primarily owing to satellite altimetry, and in 1996 we reported (NATURE, 381, 595-7) an estimate of solid-earth tidal energy dissipation and mantle Q based on combining satellite altimeter measurements with tracking observations of tidally induced satellite orbit perturbations. Tidal estimates from both reveal a small systematic difference in the quadrature component of the degree-2, order-2 spherical harmonic coefficients, which we attribute to a small lag in the earth's body tide. The formalism accounts for this lag in both the altimeter and tracking solutions and also accounts for a very small contribution from the lunar atmospheric tide. Since this original report, both altimeter and tracking estimates have improved. Recent solutions for the body-tide lag at the M2 period are 0.20 +/- 0.09 degrees, implying an energy dissipation of 100 +/- 50 gigawatts and a solid-earth Q of 300. Further new solutions will be discussed, as will the prospects for significantly reducing error bars and for obtaining estimates from other tides in the diurnal band.

Ray, Richard D.↗

Mars Gravity Field Model Development from Mars Global Surveyor Tracking Data

Since Feb. 99 the MGS spacecraft has been in a near circular orbit at about 400 km altitude. The MGS has been regularly tracked by the Deep Space Network (DSN) at X-band and for a 3 week period in February was tracked almost continuously for an intensive gravity modeling activity that would form the basis of the orbital computations for the rest of the mission. The data collected during this calibration period and the earlier SPO and Hiatus periods have now been used to develop a new gravity field model for Mars that is showing considerable new detail in both the northern and southern hemispheres. Until February no data at 400 km altitude or lower had been acquired on any previous mission south of about 35S and all the previous data were of significantly lower quality. Low altitude data (-170 km) were obtained over the higher latitudes of the northern hemisphere during the SPO periods but because of the high eccentricity of the orbit nothing of similar quality was obtainable for the southern hemisphere. The new models are of spherical harmonic degree and order 70 or higher and are suggesting large anomalies are often associated with the large impact features. Gravity data have also been obtained over both the northern and southern polar ice caps. The MGS orbit quality resulting from the use of these newer models is better than any previous Mars missions and is approaching the ten's of meter level that had been hoped would be eventually realizable.

Lemoine, F. G.↗

Preliminary Correlations of Gravity and Topography from Mars Global Surveyor

The Mars Global Surveyor (MGS) spacecraft is currently in a 400-km altitude polar mapping orbit and scheduled to begin global mapping of Mars in March of 1999. Doppler tracking data collected in this Gravity Calibration Orbit prior to the nominal mapping mission combined with observations from the MGS Science Phasing Orbit in Spring - Summer 1999 and the Viking and mariner 9 orbiters has led to preliminary high resolution gravity fields. Spherical harmonic expansions have been performed to degree and order 70 and are characterized by the first high spatial resolution coverage of high latitudes. Topographic mapping by the Mars Orbiter Laser Altimeter on MGS is providing measurements of the height of the martian surface with sub-meter vertical resolution and 5-30 m absolute accuracy. Data obtained during the circular mapping phase are expected to provide the first high resolution measurements of surface heights in the southern hemisphere. The combination of gravity and topography measurements provides information on the structure of the planetary interior, i.e. the rigidity and distribution of internal density. The observations can also be used to address the mechanisms of support of surface topography. Preliminary results of correlations of gravity and topography at long planetary wavelengths will be presented and the implications for internal structure will be addressed.

Zuber, M. T.↗

Co-Seismic Mass Dislocation and Its Effect on Earth's Rotation and Gravity

Mantle processes often involve large-scale mass transport, ranging from mantle convection, tectonic motions, glacial isostatic adjustment, to tides, atmospheric and oceanic loadings, volcanism and seismicity. On very short time scale of less than an hour, co-seismic event, apart from the "shaking" that is the earthquake, leaves behind permanent (step-function-like) dislocations in the crust and mantle. This redistribution of mass changes the Earth's inertia tensor (and hence Earth's rotation in both length-of-day and polar motion), and the gravity field (in terms of spherical harmonic Stokes coefficients). The question is whether these effects are large enough to be of any significance. In this paper we report updated calculation results. The calculation uses the normal mode summation scheme, applied to 15,814 major earthquakes that occurred during 1976-1998, according to source mechanism solutions given by the Harvard Central Moment Tensor catalog. Compared to the truly large ones earlier in the century, the earthquakes we study are individually all too small to have left any discernible signature in geodetic records of Earth rotation or global gravity field. However, their collective effects continue to exhibit an extremely strong statistical tendencies. For example, earthquakes conspire to decrease J(sub 2) and J(sub 22) while shortening LOD, resulting in a rounder and more compact Earth. Strong tendency is also seen in the earthquakes trying to "nudge" the Earth rotation pole towards about 140 degree E, roughly opposite to the observed polar drift direction. The geophysical significance and implications will be further studied.

Chao, Benjamin F.↗

The Lunar Crustal Thickness from Analysis of the Lunar Prospector Gravity and Clementine Topography Datasets

The Lunar Prospector spacecraft has mapped the gravity field of the Moon to a level of resolution never achieved before, and a spherical harmonic representation to degree and order 100 is available. When combined with the topography dataset produced by the Clementine mission, the resulting Bouguer anomaly map is interpreted to model the thickness of the lunar crust. Such models are crucial to understanding the lunar thermal history and the formation of geological features such as mascon basins, several more of which have been newly discovered from this dataset. A two-layer planetary model was used to compute the variations of the depth to the lunar Moho. The thickness values ranged from near 0 to 120 km. There is significant agreement with previous work using the Clementine gravitational field data with differences in specific locations such as South Pole-Aitken Basin, for example.

Asmar, S.↗

Geoid Recovery Using Geophysical Inverse Theory Applied to Satellite to Satellite Tracking Data

This report describes a new method for determination of the geopotential, or the equivalent geoid. It is based on Satellite-to-Satellite Tracking (SST) of two co-orbiting low earth satellites separated by a few hundred kilometers. The analysis is aimed at the GRACE Mission, though it is generally applicable to any SST data. It is proposed that the SST be viewed as a mapping mission. That is, the result will be maps of the geoid or gravity, as contrasted with determination of spherical harmonics or Fourier coefficients. A method has been developed, based on Geophysical Inverse Theory (GIT), that can provide maps at a prescribed (desired) resolution and the corresponding error map from the SST data. This computation can be done area by area avoiding simultaneous recovery of all the geopotential information. The necessary elements of potential theory, celestial mechanics, and Geophysical Inverse Theory are described, a computation architecture is described, and the results of several simulations presented. Centimeter accuracy geoids with 50 to 100 km resolution can be recovered with a 30 to 60 day mission.

Gaposchkin, E. M.↗

Resonant Third-Degree Diurnal Tides in the Seas Off Western Europe

Third-degree diurnal tides are estimated from long time series of sea level measurements at three North Atlantic tide gauges. Although their amplitudes are only a few mm or less, their admittances are far larger than those of second-degree diurnal tides, just as Cartwright discovered for the M(sub 1) constituent. The tides are evidently resonantly enhanced owing to high spatial correlation between the third-degree spherical harmonic of the tidal potential and a near-diurnal oceanic normal mode that is most pronounced in the North Atlantic. By estimating the ocean tidal response across the diurnal band (5 tidal constituents plus nodal modulations), the period and Q of this mode and one nearby mode are estimated.

Ray, Richard D.↗

Orsted Initial Field Model

Magnetic measurements taken by the Orsted satellite during geomagnetic quiet conditions around January 1, 2000 have been used to derive a spherical harmonic model of the Earth's magnetic field for epoch 2000.0. The maximum degree and order of the model is 19 for internal, and 2 for external, source fields; however, coefficients above degree 14 may not be robust. Such detailed models exist for only one previous epoch, 1980. Achieved rms misfit is 2 nT for the scalar intensity and 4 nT for the vector components perpendicular to the magnetic field. This model is of higher detail than the IGRF 2000, which for scientific purposes related to the Orsted mission it supersedes.

Olsen, N.↗